Within the framework of the Landau-Ginzburg-Devonshire approach we explore the impact of elastic anisotropy, electrostriction, flexoelectric couplings, and mismatch strain on the domain structure morphology in ferroelectric core-shell nanoparticles of spherical shape. We perform finite element modelling (FEM) for multiaxial ferroelectric nanoparticle cores covered with an elastically-isotropic soft or elastically-anisotropic rigid paraelectric shell, with and without mismatch strains. In the case of a core covered with a soft shell, the FEM results show that at room temperature a single polarization vortex with a dipolar kernel can be stable if the electrostriction coupling is weak. With increasing anisotropic electrostriction coupling, the vortex disappears and is replaced by complex flux-closure structures, which are formed in the equatorial plane and transform into an elongated vortex with a central 180-degree domain wall near the core poles. This complex domain morphology develops in the core due to the anisotropic electrostriction, and the flexoelectric coupling leads to an additional curvature and twist of the polarization isosurfaces. In contrast to this, FEM performed for a core covered with a rigid shell shows that, at room temperature, the anisotropic elastic properties of the shell can stabilize vortex-like structures with three flux-closure domains, which gradually cross in the equatorial plane of the core and transform into 120-degree domains near the core poles. The flexoelectric coupling leads to a noticeable curling of the flux-closure domain walls. A mismatch strain compensates the curling of the flux-closure domains in the core confined by the elastically-anisotropic rigid shell. Our analysis reveals different types of topological defects, namely Bloch point structures (BPS) and Ising lines, that form in a ferroelectric core covered with a soft or rigid shell.