This work investigates the transient effects of a single subsurface inclusion over the pressure, film thickness, and von Mises stress in a line elastohydrodynamic lubrication contact. Results are obtained with a fully-coupled finite element model for either a stiff or a soft inclusion moving at the speed of the surface. Two cases analyzed consider the inclusion moving either at the same speed as the mean velocity of the lubricant or moving slower. Two additional cases investigate reducing either the size of the inclusion or its stiffness differential with respect to the matrix. It is shown that the well-known two-wave elastohydrodynamic lubrication mechanism induced by surface features is also applicable to the inclusions. Also, that the effects of the inclusion become weaker both when its size is reduced and when its stiffness approaches that of the matrix. A direct comparison with predictions by the semi-analytical model of Morales-Espejel et al. ( Proc IMechE, Part J: J Engineering Tribology 2017; 231) shows reasonable qualitative agreement. Quantitatively some differences are observed which, after accounting for the semi-analytical model's simplicity, physical agreement, and computational efficiency, may then be considered as reasonable for engineering applications.