We derive a simple tensor algebraic expression of the modified Eshelby tensor for a spherical inclusion embedded in an arbitrarily anisotropic matrix in terms of three tensor quantities (the fourth-order identity tensor, the elastic stiffness tensor, and the Eshelby tensor) and two scalar quantities (the inclusion radius and interfacial spring constant), when the interfacial damage is modelled as a linear-spring layer of vanishing thickness. We validate the expression for a triclinic crystal involving 21 independent elastic constants against finite element analysis.