Skip to main content

Spatial distance between sites of sampling associated with genetic variation among Neospora caninum in aborted bovine foetuses from northern Italy



Neospora caninum, a coccidian protozoan, represents an important cause of bovine abortion. Available N. caninum strains show considerable variation in vitro and in vivo, including different virulence in cattle. To which extent sexual recombination, which is possible in the intestines of domestic dogs and closely related carnivores as definitive hosts, contributes to this variation is not clear yet.


Aborted bovine foetuses were collected between 2015 and early 2019 from Italian Holstein Friesian dairy herds suffering from reproductive problems. A total of 198 samples were collected from 165 intensive farms located in Lombardy, northern Italy. N. caninum samples were subjected to multilocus-microsatellite genotyping using ten previously established microsatellite markers. In addition to our own data, those from a recent study providing data on five markers from other northern Italian regions were included and analysed.


Of the 55 samples finally subjected to genotyping, 35 were typed at all or 9 out of 10 loci and their individual multilocus-microsatellite genotype (MLMG) determined. Linear regression revealed a statistically significant association between the spatial distance of the sampling sites with the genetic distance of N. caninum MLMGs (P < 0.001). Including data from this and a previous North Italian study into eBURST analysis revealed that several of N. caninum MLMGs from northern Italy separate into four groups; most of the samples from Lombardy clustered in one of these groups. Principle component analysis revealed similar clusters and confirmed MLMG groups identified by eBURST. Variations observed between MLMGs were not equally distributed over all loci, but predominantly observed in MS7, MS6A, or MS10.


Our findings confirm the concept of local N. caninum subpopulations. The geographic distance of sampling was associated with the genetic distance as determined by microsatellite typing. Results suggest that multi-parental recombination in N. caninum is a rare event, but does not exclude uniparental mating. More comprehensive studies on microsatellites in N. caninum and related species like Toxoplasma gondii should be undertaken, not only to improve genotyping capabilities, but also to understand possible functions of these regions in the genomes of these parasites.

Graphical Abstract


Neospora caninum is a protozoan coccidian parasite closely related to Toxoplasma gondii and Besnoitia besnoiti. It has a worldwide distribution and causes foetal losses or stillbirth in livestock, especially cattle [1,2]. Domestic dogs and other phylogenetically closely related wild carnivores like coyotes, wolves, or dingoes are the only known definitive hosts of N. caninum [1,3,4,5], i.e. hosts, in which sexual recombination can take place. The main mode of transmission in cattle, which are important intermediate hosts of the parasite, seems to be endogenous vertical transmission [1,2], i.e. transmission from a persistently or chronically infected dam to her offspring [6]. Exogenous vertical transmission to a foetus or an unborn calf occurs in cattle that became infected during pregnancy after the ingestion of oocysts shed by a definitive host [6,7]. Transmission via oocysts by definitive hosts, most likely farm dogs, seems to occur frequently in cattle populations. These infections are the likely cause of epidemic abortions in cattle herds [8]. Such abortion storms were reported to be associated with both (i) low avidity, i.e. a recently established antibody response against N. caninum [8,9], and (ii) an identical microsatellite pattern in N. caninum detected in several foetuses of a single herd [10].

A recent study on the genomes of 50 N. caninum isolates collected worldwide from a wide range of hosts using 19 linked and unlinked genetic markers showed that there is only a single genotype of N. caninum worldwide [11]. Moreover, whole genome sequencing of seven isolates from two continents revealed less than 104 bi-allelic single nucleotide polymorphisms (SNPs), which is very little compared to the situation in the genome of the closely related apicomplexan parasite T. gondii with > 106 SNPs between the compared strains [11]. Since > 50% of the SNP clustered in six haploblocks which had been already partially observed in a former study [12], it was concluded that uni-parental reproduction together with a non-sexual expansion formed the actual worldwide N. caninum population rather than non-sexual expansion alone [11]. The extremely limited genetic diversity was explained by genetic bottlenecks during the domestication of cattle in the Near East probably about 10,000 years ago [11]. Breeding and moving specific cattle breeds including the Holstein Friesian breed examined here during the following centuries may also have contributed [13]. The domestication of dogs may have formed another bottleneck, which shaped the N. caninum population that exists today [11].

Nevertheless, N. caninum is less uniform than the previously mentioned study [11] suggests. In fact, different N. caninum isolates show large differences in vitro [14,15,16] and in vivo [17,18,19], also in cattle [20,21,22,23]. The pioneering work of Regidor-Cerrillo et al. [24] contributed a number of microsatellite markers that allow fingerprinting N. caninum isolates or DNAs and undertake population studies.

We aimed at genotyping N. caninum in aborted bovine foetuses from Lombardy, one of the most important dairy cattle production areas in Italy [25]. In this region, where the Italian Holstein Friesian breed prevails, N. caninum was suspected as an important cause of abortion (MTM, unpublished data). Moreover, a high N. caninum seroprevalence was revealed in cattle from northern Italy [26,27]. We therefore determined the proportion of N. caninum PCR-positive aborted foetuses in this area and characterised the available isolates by multi-locus microsatellite genotyping.

A previous worldwide study [28], a South-American study [29], a local study from Spain [30], and a recent study conducted in northern Italy [31] suggested the existence of N. caninum sub-populations specific for particular countries or regions. The North-Italian study conducted predominantly in two areas (Piedmont and Veneto-Trento) close to our study area suggested at least three sub-populations. One aim was to find out how N. caninum isolates in aborted foetuses from Lombardy fit into the pattern of these sub-populations. To understand the reasons for genetic differences among N. caninum of one region, we used detailed geographic information on the sites, where N. caninum-positive foetal material had been sampled to assess, whether genetic differences were associated with the spatial distances of farms. A third aim was to look at farms with repeated submission of N. caninum-positive foetuses and to investigate, if the microsatellite data would support either uni-parental or multi-parental sexual reproduction.

Material and methods

Study area

All farms, from which aborted foetuses had been sampled, were located in the “Bassa Padana” in the Po valley. This area comprises the territory of the provinces of Lodi, Cremona, Mantova, Pavia, and the south of the provinces of Milano and Brescia (Fig. 1). This is one of the largest dairy production areas in Italy and stands out for the high density of cattle farms, mainly based on Italian Holstein Friesian under the intensive production system. In particular, the Italian National Zootechnical Registry counted 5446 dairy farms hosting 1,071,164 animals in the Lombardy region, corresponding to 20.7% and 40.8% of all farms and animals in Italy, with an annual production of 5,215,408 tons of milk (; last access 06.08.2020). Most farms in the area host between 100 and 500 animals, but several farms hold more than 500 cattle (National Zootechnical Database,; last access 06.08.2020). In these dairy cattle farms, mainly animals born in the farm are used for replacement.

Fig. 1
figure 1

Genetic distance between N. caninum isolates from bovine foetuses collected in Lombardy, northern Italy, is associated with the geographical distance of the sampling sites. a Graphical representation of a linear regression model to associate Bruvo’s genetic distance of N. caninum in bovine foetuses and spatial distance of sampling sites (km). b Dendrogram representing the results of K-means clustering using Bruvo’s distance with non-parametric bootstrapping; the analysis identified at least four groups (different colours of sample designation, consisting of a string of sample ID, village name, and farm number), separated by bootstrap values > 50. c Map showing farm locations, using the same colours as for the groups in the dendrogram

The sites have an average altitude of about 69 m above sea level. The climate is continental, typical of the Po valley, with hot muggy summers with a few thunderstorms and cold and foggy winters with some snow. There is a large annual thermal excursion with a mean maximum temperature of 25–28 °C and a mean minimum temperature of − 1 to – 2 °C. Rainfall is distributed over the year with peaks in spring and autumn. The average total annual rainfall ranges between 700 and 1200 mm.

Sample collection

Aborted bovine foetuses were collected between 2015 and early 2019 from Italian Holstein Friesian dairy herds suffering from reproductive problems; previous unpublished diagnostic results (i.e. by serology or PCR) had suggested endemic N. caninum associated abortion problems in these herds. Overall, 198 samples were collected from 165 intensive farms located in the Po valley (Provinces of Cremona, Lodi, Mantova, and Milano). Abortion events had been notified by the farm veterinarian, foetuses were maintained refrigerated and collected directly from the farms or submitted to the laboratory. In the laboratory, a pool of organs was prepared from each aborted foetus, which included brain, lung, and liver. Pools were used to streamline the diagnostic procedure and to keep costs at a minimum. Target organs, i.e. brain, lung, and liver, were isolated from aborted foetuses, and about 10 g of each organ was pooled into a 50-ml tube and subsequently homogenized mechanically after the addition of 10 ml of 1 × phosphate buffered saline using a laboratory homogenizer (Omni GLH 850, OMNI International, Kennesaw, GA, USA) with Soft Tissue Omni Tips Plastic Probes (Omni GLH 850, OMNI International, Kennesaw, GA, USA) to avoid cross-contamination. Homogenized pools were stored at – 20 °C until DNA extraction. The geographical localization (municipality and province as well as the geographic coordinates) of each farm with abortion episodes was recorded along with the data for each aborted foetus.

DNA extraction

Genomic DNA was extracted from tissue homogenates of aborted foetuses using a commercial kit (NucleoSpin® Tissue, Macherey-Nagel, Düren, Germany), according to the manufacturer’s instructions. To optimize DNA extraction from aborted foetuses, a 200 µl aliquot of the homogenate was pre-lysed with buffer T1 and proteinase K as recommended for the NucleoSpin Tissue Kit, by scaling up eight times the volumes used for the initial digestion. After digestion (56 °C, overnight), 230 µl of the final suspension (1840 µl) was taken and the protocol of the NucleoSpin Tissue kit followed as recommended. A negative extraction control was included in each batch of ~ 24 samples, which was subsequently tested by real-time PCR in the same way as the samples from aborted foetuses were analysed.

Screening N. caninum real-time PCR

Foetal DNA samples were initially subjected to a N. caninum-specific real-time PCR with the primers and probes targeting the Nc5 gene as previously described [32,33]. Primers were used at a final concentration of 600 nM and the fluorogenic probe (NeoProbe) at 100 nM. Ten μl of 2 × iQ Supermix (Bio-Rad Laboratories GmbH, Feldkirchen, Germany) was added to each reaction. For PCR analysis, 5 μl of DNA from each sample was used in a 20 μl reaction on a 96-well reaction plate. The thermal cycling protocol consisted of an initial polymerase activation and DNA denaturation at 95 °C for 5 min, followed by 46 cycles of amplification including denaturation at 95 °C for 10 s and annealing/extension at 58 °C for 30 s. NC-1 [34] was included as a positive control, while water PCR Reagent (Sigma-Aldrich) was used as a negative control; negative extraction but no endogenous controls were included.

Microsatellite typing

Samples confirmed as positive for N. caninum were subjected to multilocus-microsatellite typing using ten microsatellite markers (MS1B, 2, 3, 5, 6A, 6B, 7, 10, 12, and 21). Nested-PCR (n-PCR) techniques were used for the amplification of microsatellite-containing regions from N. caninum DNA from aborted foetuses. Microsatellites MS2 and MS10 were amplified and sequenced using previously described primers and protocols [35]. The microsatellites MS1B, 3, 5, 6A, 6B, 7, 12, and 21 were analysed by nested-PCR and fragment lengths determined by capillary electrophoresis [10]. Forward primers of the internal PCRs were labelled with either 6-FAM, HEX, or NED dyes to allow the simultaneous length determination of three different microsatellites per capillary electrophoresis. In addition, previously sequenced microsatellite-containing fragments from the in vitro-grown N. caninum NC-1 strain [36] were included in the amplifications, the sequencing, and sizing analyses as a positive control.

For all n-PCRs, primers were used at a final concentration of 0.5 µM and dNTPs at 250 µM each (Amershan Biosciences, Piscataway, NJ, USA). DyNAzyme II DNA polymerase (Finnzymes, Espoo, Finland) was added at 1 U/25 µl with the provided buffer. The reaction mix was supplemented with bovine serum albumin at a concentration of 20 µg/ml; 1.5 µl of genomic DNA or 1 µl of amplification product was used as template for the external and internal PCR amplification steps, respectively. Water PCR Reagent (Sigma-Aldrich) was used as a negative control. DNA from cell culture-derived N. caninum tachyzoites (NC-1) was used as a positive control. The PCR reactions were performed in a thermal cycler (Eppendorf Mastercycler, Personal Thermal Cycler). The external PCR was performed with an initial denaturation step of 94 °C for 5 min, followed by 35 cycles of denaturation (1 min at 94 °C), annealing (1 min at 60 °C for all MS, except for MS2: 50 °C, and MS10: 54 °C) and extension (1 min at 72 °C), and a final extension step at 72 °C for 10 min. The internal PCR consisted of an initial denaturation of 94 °C for 5 min, followed by 35 cycles of denaturation: 1 min at 94 °C, annealing: 1 min at 60 °C (MS2, 10) or 30 s at 52 °C (MS1B, 3, 5, 6A, 21), 54 °C (MS6B, 12) or 60 °C (MS7) and extension: 1 min at 72 °C (MS2, 10), 30 s at 70 °C (MS1B, 3, 5, 6A, 6B, 12, 21) or 65 °C (MS7), and a final extension step at 72 ° C (MS2,10), 70 °C (MS1B, 3, 5, 6A, 6B, 12, 21) or 65 °C (MS7) for 10 min.

The amplification products were visualized after electrophoresis in 2% agarose gels stained with ethidium bromide using a 100 bp DNA ladder (Invitrogen GmbH, Karlsruhe, Germany) as a size standard.

MS2 and MS10 were Sanger-sequenced because these microsatellite markers consist of three sub-loci (or motifs) each. This situation could only be resolved by sequencing and not by capillary electrophoresis. For sequencing MS2 and MS10, bands of the expected size were excised from the agarose gels and purified with a commercial kit (NucleoSpin® Gel and PCR Clean-up, Macherey-Nagel, Düren, Germany), following the manufacturer’s instructions. Purified amplification products were then cloned into a commercially available vector (pGEM®-T Easy Vector System I, Promega, Mannheim, Germany) and used to transform chemically competent Escherichia coli (OneShot TOP10, Thermo Fisher Scientific, Langenselbold, Germany). The transformed E. coli were cultivated, and the plasmid DNA was subsequently collected using a commercial kit (QIAprep Spin Miniprep Kit, Qiagen, Hilden, Germany), according to the manufacturer’s instructions. Sequencing was finally performed using the BigDye Terminator v1.1 Cycle Seq. Kit (Thermo Fisher Scientific, Langenselbold, Germany) and passage with NucleoSEQ Columns (Macherey-Nagel, Düren, Germany) for cleaning up nucleic acids, in an ABI 3130 capillary sequencer (Thermo Fisher Scientific, Langenselbold, Germany).

For all other MS, N. caninum-positive samples were analysed on an ABI 3130 capillary sequencer (Thermo Fisher Scientific, Langenselbold, Germany). A ROX dye-labelled standard (GeneScan™ 500 ROX™ dye Size Standard, Thermo Fisher Scientific, Langenselbold, Germany) was included in each analysis as a size reference. The results were analysed using the Geneious 1.11.5 software (Biomatters, Inc., 2365 Northside Dr., Suite 560, San Diego, CA 92108, USA).

Statistical analysis

Differences in N. caninum-specific real-time PCR results (Ct values) in relation to the number of microsatellite markers typed were assessed using the command “wilcox.test” from the package “stats” in R, version 3.5.3 (R Foundation for Statistical Computing, Vienna, Austria; Microsatellite marker-specific differences in the number of repeats between N. caninum isolates from different regions were assessed using the command “pairwise.wilcox.test” from the package “stats” in R, version 3.5.3 “BH adjusted” option because of multiple testing [37].

For calculating the relative genetic distance between microsatellite genotypes, we used a method that specifically permits analysing organisms with different polyploidy levels and takes stepwise mutational processes into account [38]. Pairwise Bruvo’s genetic distance was calculated using microsatellite data of all DNA samples from aborted foetuses in Lombardy, for which a complete MLMG was obtained (n = 20). Since MS10 and MS2 were combined microsatellite markers, which consist of three sub-loci (or motifs) each, data on these sub-loci were included separated into the data set. Calculations were done by R, version 3.5.3., employing “Poppr”, version 2, a R package designed for the analysis of populations with mixed modes of reproduction [39], by using the command “bruvo.dist”.

To visualize the relatedness of individual MLMGs, a dendrogram was established using the “Poppr” package representing the results of K-means clustering using Bruvo’s distance with non-parametric bootstrapping, a process that consists of randomly sampling loci with replacement, recalculating the dendrogram, and grossing up the bootstrap support (in percent of success).

For testing of a standardized index of association (\(I_{{\text{A}}}^{{\text{S}}}\)) by multilocus linkage disequilibrium (LD) among different genotypes including nine loci the LIAN v 3.7 web interface program ( was applied as described [28,40].

Testing population differentiation was done by likelihood ratio G‐statistic [41], using the function “test.between” from the R-package “HierFstat” [42].

Geographic distance based on the geographic coordinates of the farms were calculated (straight-line distance, accounting for curvature of the earth) using the command “spDists” of the package “sp”, version 1.4-2, in R. To determine the relatedness of Bruvo’s distance, geographic distance and the time between samplings linear regression was performed using the “lm” command in R, version 3.5.3 from the package “stats”.

The eBURST software [43] was used to generate networks based on MLMGs using the double-locus option (DVL—at least 7 shared loci of 9). Samples with missing data were excluded from the analysis. As the eBURST analysis was performed on samples from this study and those from a previous study, the analysis was restricted to those microsatellite loci that were used in both studies. Since MS10 is a combined microsatellite marker, which consists of three sub-loci (or motifs), these sub-loci were analysed separately. The same data (MLMGs data collected in the present study on foetuses from Lombardy and data from a previous study) were also analysed by a Principal Coordinates Analysis (PCoA) using the command “dudi.pco” employing the R package “Poppr”, version 2.

Figures were assembled using R, version 3.5.3 or 4.0.0 (packages “ggplot2”, “reshape” and “scales”).


Overall N. caninum DNA findings in aborted foetuses

Out of 198 aborted foetuses 55 were positive for N. caninum by real-time PCR, yielding a prevalence of 27.8% (presence of the parasite in pooled tissue homogenates of brain, lung, and liver); 43 farms recorded at least one positive foetus (26.1%). Overall, 55 samples from 43 herds of foetuses collected were subjected to microsatellite typing.

Relationship between real-time PCR results and success in genotyping

The majority of the 55 DNA samples subjected to multilocus microsatellite typing were typed at all 10 (n = 20) or at 9 (n = 15) loci. Nine DNA samples could be typed at eight (n = 2), seven (n = 2) or six (n = 6) loci. The remaining ten DNA samples could only be typed at five (n = 3), four (n = 3), three (n = 2), two (n = 1), or one (n = 1) of the loci.

The numbers of genotyped loci were related to the amount of parasitic DNA in the sample as reflected by the Ct values determined by real-time PCR (Table 1).

Table 1 N. caninum DNA content as reflected by the Ct value in samples of aborted foetuses and success in genotyping

Microsatellite typing on individual farms at different time points revealed differences in the MLMGs

In seven farms, more than one foetus had been sampled and detected N. caninum DNA typed in the study period from 2015 to 2018. Most of the differences in typing comprised only one or two repeat units in individual microsatellites (Table 2). However, in five farms, more prominent differences per locus (i.e. differences of 3 or even more repeat units per locus) were observed. In herd 6, i.e. the farm with the largest number of foetuses analysed, 48.1% of the comparisons of individual loci revealed differences (Tables 2, 3). These more prominent differences affected the loci MS7 and MS21 (chromosome 7a), MS10 (MS10.2, chromosome 8), and MS6A (chromosome 10) (Table 3).

Table 2 Comparison of the MLMG patterns in dairy cattle farms in northern Italy, where more than one Neospora caninum-positive foetus had been sampled and genotyped during the study period
Table 3 MLMG patterns on dairy cattle farms from northern Italy where more than one Neospora caninum-positive foetus had been sampled during the study period, which showed differences (bold numbers) in length for ≥ 3 repeats at microsatellite loci on various chromosomes (Chr2, 7a, 8, 9, 10, 12)

Relationship between the spatial distance of sampling sites and the genetic distance of N. caninum isolates

For 20 of the N. caninum-positive samples of bovine foetuses from Lombardy, which had been typed at the complete set of microsatellite loci (n = 10, this study), pairwise Bruvo’s genetic distances were calculated. Bruvo’s genetic distance and the geographical distance between sampling sites of the individual N. caninum-positive foetuses (i.e. based on the geographical coordinates of 17 farms) were tested for correlation. Linear regression revealed that Bruvo’s genetic distance correlated statistically significantly with the geographical distance between the sampling sites (Fig. 1a; Table 4). The model (Model 1, Table 4) had an adjusted R2 of 19.5%. Including the number of days between the sampling of foetuses into the model did not improve the model significantly (Model 2, Table 4).

Table 4 Linear regression models to characterize the association between genetic distance (Bruvos’s genetic distance) of N. caninum isolates and the spatial distance of the sampling sites or the time between dates of sampling

Bruvo’s genetic distance with non-parametric bootstrapping allowed the identification of at least four groups, separated by bootstrap values > 50 (Fig. 1b). The farms, from which the largest group of samples originated (Fig. 1b, red), were located with two exceptions central-north in the Cremona district (Fig. 1c, red). The farms, from which the remaining samples had been derived, were located in the far north (Fig. 1c, green), in the south of the Cremona district (Fig. 1c, blue), or in the neighbouring district of Mantova (Fig. 1c, yellow).

An eBURST and PCoA analysis revealed different groups or clusters of N. caninum MLMGs in northern Italy

In addition to the microsatellite data established in this study (Lombardy, n = 25), further data on 50 N. caninum isolates from northern Italy covering mainly other regions than Lombardy (Piedmont, n = 17; Veneto-Trento, n = 6, including also one goat isolate), Lombardy (n = 1), and a further North-Italian bovine N. caninum isolate (n = 1) were analysed. These additional data were available from a recent study [31]. Since the latter study had employed a set of microsatellite markers that overlapped only partially with the one we used, our analysis was restricted to MS6A, MS6B, MS10, MS12, and MS21, for which data were available from the present and the previous study. Linkage disequilibrium (LD) was assessed for the entire population excluding the goat sample; results (\(I_{{\text{A}}}^{{\text{S}}}\) = 0.0411, VD = 1.8249, L = 1.5508, = 0.0411) indicated LD because VD > L [40].

Network analysis using eBURST with double locus variation (DLV) among six loci and three MS10 motifs or sub-loci (MS10.1-3) revealed that many N. caninum MLMGs from Lombardy (n = 11/25; this study) clustered separately (eBURST G4) from those obtained from Piedmont (n = 10/17, eBURST G1, G2, G3) or Veneto-Trento (n = 3/6, eBURST G1, G2) (Fig. 2a). The grouping of isolates from Piedmont and Veneto-Trento had already been assessed by others [31]. This previous grouping of particular N. caninum isolates into the eBURST groups G1, G2, and G3 matched perfectly (n = 10/10) with our grouping (Additional file 1: Table S1). Genotyping of an Italian isolate, for which the province of origin was not known (violet dot), revealed results close to most isolates from Lombardy (group G4). The full representation of the MLMST network using the MST option in eBURST, including the complete data set (n = 49 N. caninum samples), shows the G4 samples in the centre with most of the Veneto-Trento samples on the right side (G1) and the Piedmont samples both left (G2, G3) and right (G1) in the network (Fig. 2b).

Fig. 2
figure 2

eBURST analysis of 49 northern Italian Neospora caninum samples. Specimens from Lombardy (red dots) tend to cluster separately from those from Piedmont (green dots) or Veneto-Trento (blue dots). a Analysis using the Double Locus Variant (eBURST DLV) option shows that N. caninum multilocus genotypes can be separated into four groups (G1–4). Genotyping of an Italian isolate for which the province of origin was not known (violet dot) revealed results close to most isolates from Lombardy (group G4). b The full MST option in eBURST, including the complete data set (n = 50 N. caninum samples) shows the location of the individual groups within a network of all samples. The analysis was restricted to microsatellite markers available from this and a previous study [31]. Only samples that could be typed for all microsatellite markers were included. Moreover, MS10, which combines variation in three separate motifs (sub-loci), was analysed per each motif individually. Groups G1–G3 resemble the grouping reported in [31], while G4 represents a new group including 11 of 25 bovine N. caninum samples from Lombardy that were added to the analysis by this study

Investigating N. caninum-positive samples (n = 48, excluding a sample for which the province of origin was unknown and the goat sample) by PCoA revealed a clear Axis 1 separation of the samples from Piedmont (Fig. 3) and Lombardy (Fig. 3). While n = 17/17 samples from Piedmont were located on the left side, all from Lombardy (n = 17/25) were on the right side (Fig. 3). The majority (n = 3/5) of samples from Veneto-Trento (Fig. 3) were left or right in the PCoA graph.

Fig. 3
figure 3

Principle coordinate analysis (PCoA) separates north Italian N. caninum based on their MLMGs into clusters. Data from the present and a previous study conducted in cattle from northern Italy [31] were used for this analysis. The two dominant eigenvalues of a PCoA based on the allele-sharing coefficient separated all N. caninum-positive samples isolated from Piedmont (PMN) and most from Veneto-Trento (VEN), so that they clustered on the left side of the PCoA graph, while most of the samples from Lombardy (LOM) were located on the right side of the PCoA graph. Groups G1–4 (represented by circles with numbers 1–4), as determined by eBURST, were located clearly in separate positions. While G1–3 was found on the left side, G4 (exclusively Lombardy samples) was localised on the right side of the PCoA graph

When the eBURST grouping (DLV option) was compared to the PCoA result, eBURST G1 and eBURST G3 were located separately in the upper left quarter of the PCoA graph, while eBURST G2 was located in the lower left part (Fig. 3). Neospora caninum samples that had been separated by eBURST into G4 (exclusively Lombardy samples) clustered in both the upper and lower corners on the right side of the PCoA graph (Fig. 3).

Differences between the N. caninum populations in northern Italy were assessed using FST analysis values (F) and Nei’s unbiased genetic distance (D) (Table 5). Both results for F and D suggest that there is a statistically significant genetic difference between N. caninum from Piedmont and Lombardy as well as between Lombardy and Veneto-Trento (Table 5).

Table 5 Genetic differences between northern Italian N. caninum

The graphical representation of the heterogeneity in numbers of repeat units, stratified for various loci, suggested differences between bovine N. caninum sampled in the North Italian regions of Lombardy, Piedmont, and Veneto-Trento (Fig. 4). The range between the lowest and highest number of repeats was largest for MS10.2, followed by MS6A, MS5, and MS7. However, the pairwise Wilcoxon test (BH adjusted P value) revealed statistically significant differences between regions only for MS5, MS6B, MS7, MS10.1, and MS10.2 (Fig. 4).

Fig. 4
figure 4

Differences in the heterogeneity of the numbers of repeat units characteristic for Neospora caninum microsatellite loci examined in the present and a previous study conducted in cattle from northern Italy [31], i.e. in the regions Lombardy (LOM), Piedmont (PMN), and Veneto-Trento (VEN). The range between the lowest and highest number of repeats was largest for MS10.2 (17 repeats), followed by MS6A (13), MS5 (9) and MS7 (9). For the remaining repeats, the range was four (MS12, MS21, MS10.1, MS10.3) or three repeats (MS6B). Pairwise Wilcoxon test (BH adjusted P value) revealed statistically significant differences in MS5 (LOM vs VEN), MS6B (LOM vs PMN, LOM vs VEN), MS7 (LOM vs PMN), MS10.1 (LOM vs VEN, PMN vs VEN) and MS10.2 (LOM vs PMN, LOM vs VEN)


This study focused on microsatellite multilocus genotyping of N. caninum in bovine foetuses sampled in Lombardy, northern Italy, between 2015 and early 2019 from dairy herds suffering from reproductive problems. The farms selected for this study had a background of N. caninum infections in the past as the results of serological testing of a limited number of cattle per herd had indicated; 27.8% of aborted foetuses (55/198) sampled in 43/165 farms tested positive for N. caninum, which shows that this parasite is present in many herds of this region. Because pooled samples of brain, lung, and liver had been used to streamline the analytic process, the diagnostic sensitivity might have been reduced and the analysis of individual and further organs (including also cotyledons) could have increased the proportion of positive findings. Since it was not the aim of this study to assess the importance of N. caninum as a cause of bovine abortion, we did not collect information on differential diagnoses, i.e. on potential other infectious or non-infectious causes of bovine abortion. It is therefore possible that not all N. caninum-positive abortions observed herein were caused by neosporosis. In addition to PCR, further diagnostic steps would have been necessary to confirm neosporosis, including histological and serological examinations as well as the careful exclusion of other causes of abortion [44].

A similar study had also focused on MLMGs in N. caninum in cattle and a goat from northern Italy [31]. Both beef (Piedmontese breed) and dairy (Italian Holstein Friesian breed) as well as crossbreeds or local cattle breeds were included, whereas our work was focused on Italian Holstein Friesian cattle under intensive farming conditions. Regarding the geographical area, this previous study had concentrated on bovine foetuses from the regions of Piedmont and Veneto-Trento (n = 38) and included only a few samples (n = 2) from Lombardy. In addition, the previous study lacked detailed data on the geographic origin of samples. The regions of Piedmont and Veneto-Trento neighbour our study area; Piedmont is located west and Veneto-Trento east from Lombardy. On the basis of microsatellite typing results, it has been hypothesized that local N. caninum subpopulations exist in Europe [30,32].

Entirely unexpected was our finding that geographic distance in such a small area as Lombardy (< 30,000 km²) may at least partially explain the MLMG-based genetic distance of N. caninum isolates. This suggests that MLMGs remain relatively constant over time in particular herds and only slowly disperse in an area, probably due to animal trade. As it can be assumed that trade between herds was mostly local, this may have resulted in the observed correlation of the geographic distance with the MLMG-based genetic distance between N. caninum isolates. Within herds, N. caninum seems to be transmitted endogenously mainly, i.e. from latently infected dams to their offspring. However, even an exogenous transmission via oocysts shedding definitive hosts (e.g. farm dogs or wolves) could explain the conserved microsatellite patterns we observed, if the definitive hosts were infected only with a single strain, i.e. the local one, which results in a so-called uniparental mating of N. caninum [11]. Both endogenous and exogenous transmission without or only limited contact to neighbouring herds or infectious material might have resulted in local N. caninum sub-populations.

All cattle analysed in our study in Lombardy belonged to the Italian Holstein Friesian breed, which was founded in the late twentieth century by the importation of cattle from The Netherlands and North America [45,46]. This breed originates from the north of The Netherlands and Germany and became the predominant dairy breed world-wide; present Holstein Friesian cattle represent crosses of the Dutch Friesian and North-American Holstein lines [45,46]. It has been hypothesized that individual N. caninum strains may have been imported together with the Holstein Friesian cattle breed into Italy [31]. Although MLMGs of N. caninum from northern Italy showed a relationship to those from Germany [31], it is almost impossible to clarify, whether the N. caninum genotypes now observed in Italian Holstein Friesian cattle originated from Italy (i.e. from local cattle breeds) or were introduced by the importation of cattle from abroad. Attempts to find genetic differences between N. caninum isolates from Holstein Friesian and from local breeds (i.e. the Piedmontese breed), which could have supported the hypothesis of importation, failed [31]. Interestingly, it has been observed that Argentinian N. caninum isolates are still related to those from Spain, even long after the importation of the first Iberian cattle into South America in the fifteenth century [28].

In a few herds, we were able to sample more than a single foetus, and in four of these herds, several months or even a period of years separated the sampling dates. In contrast to findings in bovine herds with epidemic abortions (i.e. in herds where a point source exposure to N. caninum oocysts shed by a dog was assumed [10]), particular microsatellite loci differed strongly (i.e. by 3 or even more repeats). It may be discussed whether these larger differences in single loci (MS7, MS6A) or sub-loci (MS10.2) might be an indication for sexual recombination in a definitive host (probably by a farm dog) in the past. However, typing details (Table 2) revealed that related loci on the same chromosome (like MS1B, MS6B) or even sub-loci (e.g. MS10.1 or MS10.2) were not affected. Thus, these differences are best explained by the loss or gain of repeat units, which is typical for microsatellites and caused by point or polymerase template-slippage mutations [47,48]. Nevertheless, sexual recombination events in these herds cannot be completely ruled out, if this recombination occurred with two or more separate strains (i.e. multi-parental) with very similar microsatellite patterns, e.g. the microsatellite pattern of the local N. caninum population.

Our findings in individual herds with more than one N. caninum MLMG, but also the overall analysis of the N. caninum MLMGs in northern Italy, suggest that particular microsatellite loci may be more prone to variation than others. In our study, especially the sub-locus MS10.2, followed by the loci MS6A, MS5 and MS7, showed the strongest variation in the number of repeats. It is well known that several factors contribute to these differences among microsatellite loci, including repeat number, the sequence of the repeat motif, the length of the repeat unit, and flanking sequences (reviewed in [49,50,51]). Thus, high variability in particular microsatellite loci may indicate that these loci are of limited importance for the parasite. On the other hand, the observation of changes in particular microsatellite loci raises questions regarding the possible effects of such variations on the affected parasite.

Microsatellites are highly abundant in the noncoding DNA of all eukaryotic genomes [52] and changes in these loci might be “neutral”, but microsatellites may also be located in coding regions and variation can be associated with an altered phenotype (reviewed [50]). A recent study on different Plasmodium spp. affecting humans suggested that a significant proportion (one fifth to one third) of microsatellite-related sequences are related to coding sequences. Based on gene ontology, the respective coding sequences can be involved in molecular functions like binding or in biologic processes such as metabolism or reproduction [53]. Although we are far from understanding the biological relevance of microsatellites for N. caninum, it would be intriguing if future studies could comprehensively address the microsatellites of N. caninum to gain more knowledge on their function. It has so far not been possible to link particular microsatellite patterns to particular traits, e.g. virulence for foetuses [54]. Nevertheless, virulence differences have been observed among N. caninum strains also in its main intermediate host, i.e. cattle [15,17,18,20,21,22,23], and it would be fascinating to know to which extent differences in microsatellite loci contribute to virulence. Recently, highly virulent isolates of N. caninum were shown to express a subset of particular secreted proteins in more abundance [55]. The reasons for differential expression between strains could at least in part be due to differences in microsatellite loci as shown for other eukaryotic cells and organisms [52].

In T. gondii, an apicomplexan parasite closely related to N. caninum, the situation is similar. Currently, relatively conserved microsatellite loci are used to differentiate dominant clonal lineages of T. gondii, while also a number of other microsatellite loci are known, which were called finger-printing microsatellite loci. They can be used to differentiate strains of a single type on a local level [56].


Our findings confirm the concept of local N. caninum sub-populations. For the first time, we could show a correlation between the genetic distances of N. caninum isolates based on MLMGs and the geographic distances of the places, where the isolates had been obtained. Our results confirm that sexual recombination in N. caninum is a rare event. Possible reasons for this might be that endogenous vertical transmission is dominating and that the chance for a definitive host to feed on intermediate hosts with a mixed infection of viable and non-related N. caninum strains is extremely low. More comprehensive studies on microsatellites in N. caninum and related species like T. gondii should be undertaken, not only to improve genotyping capabilities, but also to understand the possible function of these regions in the genomes of these important parasites.

Availability of data and materials

Data supporting the conclusions of this article are included within the article and its additional files. The raw datasets used and analyzed during the present study are available as Additional file 1: Table S1 or from the corresponding author upon reasonable request.



Multilocus-microsatellite genotype




Principle coordinate analysis


Double locus variant








  1. Dubey JP, Hemphill A, Calero-Bernal R, Schares G. Neosporosis in animals. Boca Rotan: CRC Press; 2017.

    Book  Google Scholar 

  2. Dubey JP, Schares G, Ortega-Mora LM. Epidemiology and control of neosporosis and Neospora caninum. Clin Microbiol Rev. 2007;20:323–67.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  3. Dubey JP, Jenkins MC, Rajendran C, Miska K, Ferreira LR, Martins J, et al. Gray wolf (Canis lupus) is a natural definitive host for Neospora caninum. Vet Parasitol. 2011;181:382–7.

    Article  PubMed  CAS  Google Scholar 

  4. Gondim LF, McAllister MM, Pitt WC, Zemlicka DE. Coyotes (Canis latrans) are definitive hosts of Neospora caninum. Int J Parasitol. 2004;34:159–61.

    Article  PubMed  Google Scholar 

  5. King JS, Slapeta J, Jenkins DJ, Al-Qassab SE, Ellis JT, Windsor PA. Australian dingoes are definitive hosts of Neospora caninum. Int J Parasitol. 2010;40:945–50.

    Article  PubMed  Google Scholar 

  6. Trees AJ, Williams DJ. Endogenous and exogenous transplacental infection in Neospora caninum and Toxoplasma gondii. Trends Parasitol. 2005;21:558–61.

    Article  PubMed  Google Scholar 

  7. McCann CM, McAllister MM, Gondim LF, Smith RF, Cripps PJ, Kipar A, et al. Neospora caninum in cattle: experimental infection with oocysts can result in exogenous transplacental infection, but not endogenous transplacental infection in the subsequent pregnancy. Int J Parasitol. 2007;37:1631–9.

    Article  PubMed  CAS  Google Scholar 

  8. McAllister MM, Wallace RL, Björkman C, Gao L, Firkins LD. A probable source of Neospora caninum infection in an abortion outbreak in dairy cows. Bovine Pract. 2005;39:69–74.

    Google Scholar 

  9. Schares G, Bärwald A, Staubach C, Söndgen P, Rauser M, Schröder R, et al. p38-avidity-ELISA: examination of herds experiencing epidemic or endemic Neospora caninum-associated bovine abortion. Vet Parasitol. 2002;106:293–305.

    Article  PubMed  CAS  Google Scholar 

  10. Basso W, Schares S, Minke L, Bärwald A, Maksimov A, Peters M, et al. Microsatellite typing and avidity analysis suggest a common source of infection in herds with epidemic Neospora caninum-associated bovine abortion. Vet Parasitol. 2010;173:24–31.

    Article  PubMed  CAS  Google Scholar 

  11. Khan A, Fujita AW, Randle N, Regidor-Cerrillo J, Shaik JS, Shen K, et al. Global selective sweep of a highly inbred genome of the cattle parasite Neospora caninum. Proc Natl Acad Sci USA. 2019;116:22764–73.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  12. Calarco L, Barratt J, Ellis J. Genome wide Identification of mutational hotspots in the apicomplexan parasite Neospora caninum and the implications for virulence. Genome Biol Evol. 2018;10:2417–31.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  13. Felius M, Beerling M-L, Buchanan DS, Theunissen B, Koolmees PA, Lenstra JA. On the history of cattle genetic resources. Diversity. 2014;6:705–50.

    Article  Google Scholar 

  14. Regidor-Cerrillo J, Gomez-Bautista M, Sodupe I, Aduriz G, Alvarez-Garcia G, Del PI, et al. In vitro invasion efficiency and intracellular proliferation rate comprise virulence-related phenotypic traits of Neospora caninum. Vet Res. 2011;42:41.

    Article  PubMed  PubMed Central  Google Scholar 

  15. Dellarupe A, Regidor-Cerrillo J, Jimenez-Ruiz E, Schares G, Unzaga JM, Venturini MC, et al. Comparison of host cell invasion and proliferation among Neospora caninum isolates obtained from oocysts and from clinical cases of naturally infected dogs. Exp Parasitol. 2014;145:22–8.

    Article  PubMed  CAS  Google Scholar 

  16. Jimenez-Pelayo L, Garcia-Sanchez M, Regidor-Cerrillo J, Horcajo P, Collantes-Fernandez E, Gomez-Bautista M, et al. Differential susceptibility of bovine caruncular and trophoblast cell lines to infection with high and low virulence isolates of Neospora caninum. Parasites Vect. 2017;10:463.

    Article  CAS  Google Scholar 

  17. Rojo-Montejo S, Collantes-Fernandez E, Regidor-Cerrillo J, Alvarez-Garcia G, Marugan-Hernandez V, Pedraza-Diaz S, et al. Isolation and characterization of a bovine isolate of Neospora caninum with low virulence. Vet Parasitol. 2009;159:7–16.

    Article  PubMed  Google Scholar 

  18. Regidor-Cerrillo J, Gomez-Bautista M, Del PI, Jimenez-Ruiz E, Aduriz G, Ortega-Mora LM. Influence of Neospora caninum intra-specific variability in the outcome of infection in a pregnant BALB/c mouse model. Vet Res. 2010;41:52.

    Article  PubMed  PubMed Central  Google Scholar 

  19. Dellarupe A, Regidor-Cerrillo J, Jimenez-Ruiz E, Schares G, Unzaga JM, Venturini MC, et al. Clinical outcome and vertical transmission variability among canine Neospora caninum isolates in a pregnant mouse model of infection. Parasitology. 2014;141:356–66.

    Article  PubMed  CAS  Google Scholar 

  20. Rojo-Montejo S, Collantes-Fernandez E, Blanco-Murcia J, Rodriguez-Bertos A, Risco-Castillo V, Ortega-Mora LM. Experimental infection with a low virulence isolate of Neospora caninum at 70 days gestation in cattle did not result in foetopathy. Vet Res. 2009;40:5.

    Article  CAS  Google Scholar 

  21. Chryssafidis AL, Canton G, Chianini F, Innes EA, Madureira EH, Gennari SM. Pathogenicity of Nc-Bahia and Nc-1 strains of Neospora caninum in experimentally infected cows and buffaloes in early pregnancy. Parasitol Res. 2014;113:1521–8.

    Article  PubMed  Google Scholar 

  22. Regidor-Cerrillo J, Arranz-Solis D, Benavides J, Gomez-Bautista M, Castro-Hermida JA, Mezo M, et al. Neospora caninum infection during early pregnancy in cattle: how the isolate influences infection dynamics, clinical outcome and peripheral and local immune responses. Vet Res. 2014;45:10.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  23. Jimenez-Pelayo L, Garcia-Sanchez M, Vazquez P, Regidor-Cerrillo J, Horcajo P, Collantes-Fernandez E, et al. Early Neospora caninum infection dynamics in cattle after inoculation at mid-gestation with high (Nc-Spain7)- or low (Nc-Spain1H)-virulence isolates. Vet Res. 2019;50:1.

    Article  CAS  Google Scholar 

  24. Regidor-Cerrillo J, Pedraza-Diaz S, Gomez-Bautista M, Ortega-Mora LM. Multilocus microsatellite analysis reveals extensive genetic diversity in Neospora caninum. J Parasitol. 2006;92:517–24.

    Article  PubMed  CAS  Google Scholar 

  25. Zucali M, Tamburini A, Sandrucci A, Bava L. Global warming and mitigation potential of milk and meat production in Lombardy (Italy). J Clean Prod. 2017;153:474–82.

    Article  CAS  Google Scholar 

  26. Otranto D, Llazari A, Testini G, Traversa D, di Regalbono AF, Badan M, et al. Seroprevalence and associated risk factors of neosporosis in beef and dairy cattle in Italy. Vet Parasitol. 2003;118:7–18.

    Article  PubMed  Google Scholar 

  27. Sala G, Gazzonis A, Boccardo A, Coppoletta E, Galasso C, Manfredi MT, et al. Using beef-breed semen in seropositive dams for the control of bovine neosporosis. Prev Vet Med. 2018;161:127–33.

    Article  PubMed  Google Scholar 

  28. Regidor-Cerrillo J, Diez-Fuertes F, Garcia-Culebras A, Moore DP, Gonzalez-Warleta M, Cuevas C, et al. Genetic diversity and geographic population structure of bovine Neospora caninum determined by microsatellite genotyping analysis. PLoS ONE. 2013;8:e72678.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  29. Cabrera A, Fresia P, Berna L, Silveira C, Macias-Rioseco M, Arevalo AP, et al. Isolation and molecular characterization of four novel Neospora caninum strains. Parasitol Res. 2019;118:3535–42.

    Article  PubMed  Google Scholar 

  30. Pedraza-Diaz S, Marugan-Hernandez V, Collantes-Fernandez E, Regidor-Cerrillo J, Rojo-Montejo S, Gomez-Bautista M, et al. Microsatellite markers for the molecular characterization of Neospora caninum: application to clinical samples. Vet Parasitol. 2009;166:38–46.

    Article  PubMed  CAS  Google Scholar 

  31. Regidor-Cerrillo J, Horcajo P, Ceglie L, Schiavon E, Ortega-Mora LM, Natale A. Genetic characterization of Neospora caninum from Northern Italian cattle reveals high diversity in European N. caninum populations. Parasitol Res. 2020;119:1353–62.

    Article  PubMed  Google Scholar 

  32. Constantin EM, Schares G, Grossmann E, Sauter K, Romig T, Hartmann S. Untersuchungen zur Rolle des Rotfuchses (Vulpes vulpes) als möglicher Endwirt von Neospora caninum. Berl Munch Tierarztl Wochenschr. 2011;124:148–53.

    PubMed  Google Scholar 

  33. Legnani S, Pantchev N, Forlani A, Zini E, Schares G, Balzer J, et al. Emergence of cutaneous neosporosis in a dog receiving immunosuppressive therapy: molecular identification and management. Vet Dermatol. 2016;27:49-e14.

    Article  PubMed  Google Scholar 

  34. Dubey JP, Hattel AL, Lindsay DS, Topper MJ. Neonatal Neospora caninum infection in dogs: isolation of the causative agent and experimental transmission. J Am Vet Med Assoc. 1988;193:1259–63.

    PubMed  CAS  Google Scholar 

  35. Basso W, Schares S, Barwald A, Herrmann DC, Conraths FJ, Pantchev N, et al. Molecular comparison of Neospora caninum oocyst isolates from naturally infected dogs with cell culture-derived tachyzoites of the same isolates using nested polymerase chain reaction to amplify microsatellite markers. Vet Parasitol. 2009;160:43–50.

    Article  PubMed  CAS  Google Scholar 

  36. Dubey JP, Dorough KR, Jenkins MC, Liddell S, Speer CA, Kwok OCH, et al. Canine neosporosis: clinical signs, diagnosis, treatment and isolation of Neospora caninum in mice and cell culture. Int J Parasitol. 1998;28:1293–304.

    Article  PubMed  CAS  Google Scholar 

  37. Benjamini Y, Hochberg Y. Controlling the false discovery rate: a practical and powerful approach to multiple testing. J Roy Stat Soc Ser B (Methodol). 1995;57:289–300.

    Google Scholar 

  38. Bruvo R, Michiels NK, D’Souza TG, Schulenburg H. A simple method for the calculation of microsatellite genotype distances irrespective of ploidy level. Mol Ecol. 2004;13:2101–6.

    Article  PubMed  CAS  Google Scholar 

  39. Kamvar ZN, Tabima JF, Grunwald NJ. Poppr: an R package for genetic analysis of populations with clonal, partially clonal, and/or sexual reproduction. PeerJ. 2014;2:e281.

    Article  PubMed  PubMed Central  Google Scholar 

  40. Haubold B, Hudson RR. LIAN 3.0: detecting linkage disequilibrium in multilocus data. Bioinformatics. 2000;16:847–9.

    Article  PubMed  CAS  Google Scholar 

  41. Goudet J. HierFstat, a package for R to compute and test hierarchical F-statistics. Mol Ecol Notes. 2005;5:184–6.

    Article  Google Scholar 

  42. de Meeûs T, Goudet J. A step-by-step tutorial to use HierFstat to analyse populations hierarchically structured at multiple levels. Infect Genet Evol. 2007;7:731–5.

    Article  PubMed  CAS  Google Scholar 

  43. Feil EJ, Li BC, Aanensen DM, Hanage WP, Spratt BG. eBURST: inferring patterns of evolutionary descent among clusters of related bacterial genotypes from multilocus sequence typing data. J Bacteriol. 2004;186:1518–30.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  44. Dubey JP, Schares G. Diagnosis of bovine neosporosis. Vet Parasitol. 2006;140:1–34.

    Article  PubMed  CAS  Google Scholar 

  45. Mancini G, Gargani M, Chillemi G, Nicolazzi EL, Marsan PA, Valentini A, et al. Signatures of selection in five Italian cattle breeds detected by a 54K SNP panel. Mol Biol Rep. 2014;41:957–65.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  46. Buchanan DS. Breeds of dairy cattle (major Bos taurus breeds). In: Reference module in food science. The Netherlands: Elsevier; 2016.

  47. Schlötterer C, Tautz D. Slippage synthesis of simple sequence DNA. Nucleic Acids Res. 1992;20:211–5.

    Article  PubMed  PubMed Central  Google Scholar 

  48. Viguera E, Canceill D, Ehrlich SD. Replication slippage involves DNA polymerase pausing and dissociation. EMBO J. 2001;20:2587–95.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  49. Ellegren H. Microsatellite mutations in the germline: implications for evolutionary inference. Trends Genet. 2000;16:551–8.

    Article  PubMed  CAS  Google Scholar 

  50. Schlötterer C. Evolutionary dynamics of microsatellite DNA. Chromosoma. 2000;109:365–71.

    Article  PubMed  Google Scholar 

  51. Estoup A, Jarne P, Cornuet JM. Homoplasy and mutation model at microsatellite loci and their consequences for population genetics analysis. Mol Ecol. 2002;11:1591–604.

    Article  PubMed  CAS  Google Scholar 

  52. Bagshaw ATM. Functional mechanisms of microsatellite DNA in eukaryotic genomes. Genome Biol Evol. 2017;9:2428–43.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  53. Mathema VB, Nakeesathit S, White NJ, Dondorp AM, Imwong M. Genome-wide microsatellite characteristics of five human Plasmodium species, focusing on Plasmodium malariae and P. ovale curtisi. Parasite. 2020;27:34.

    Article  PubMed  PubMed Central  Google Scholar 

  54. Dubey JP, Schares G. Neosporosis in animals—the last five years. Vet Parasitol. 2011;180:90–108.

    Article  PubMed  CAS  Google Scholar 

  55. Roman LRS, Horcajo P, Regidor-Cerrillo J, Fernandez-Escobar M, Collantes-Fernandez E, Gutierrez-Blazquez D, et al. Comparative tachyzoite proteome analyses among six Neospora caninum isolates with different virulence. Int J Parasitol. 2020;50:377–88.

    Article  CAS  Google Scholar 

  56. Ajzenberg D, Collinet F, Aubert D, Villena I, Darde ML, French ToxoBs network g, et al. The rural-urban effect on spatial genetic structure of type II Toxoplasma gondii strains involved in human congenital toxoplasmosis, France 2002–2009. Infect Genet Evol. 2015;36:511–6.

    Article  PubMed  Google Scholar 

Download references


The authors are grateful to the veterinarians (particularly C. Basevi, C. Barisani, S. Bonfanti, M. Colombo, E. Fumagalli, G. Gelati, R. Landriscina, P. Bossi, C. Santambrogio, R. Zanchetta) who helped in the collection of bovine aborted foetuses. We thank Patrick Wysocki for providing the map of farm locations.


Open Access funding enabled and organized by Projekt DEAL. LV was funded by the German Academic Exchange Service (DAAD) Short-Term Research Grant (no. 57440917).

Author information

Authors and Affiliations



LV, MTM and GS designed the study. LV, SAZ and MTM collected the samples. LV, FH, MT and GS collected the data. LV, ALG, MM, CL, MT, PM and GS performed experiments and analysed the samples. GS and PM statistically analysed the data. LV, PM and GS interpreted the data. GS, LV, PM, FJC and MTM made major contributions to the writing of the manuscript. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Gereon Schares.

Ethics declarations

Ethics approval and consent to participate

All procedures were approved by the Institutional Animal Care and Use Committee of Università degli Studi di Milano (“Organismo Preposto al Benessere degli Animali,” Prot. no. OPBA_34_2017).

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1

: Table S1. Microsatellite data on Neospora caninum from norther Italy, including also data on intermediate host species, region, province, sampling site, sampling date, real-time PCR result, eBURST designation, and reference.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Villa, L., Maksimov, P., Luttermann, C. et al. Spatial distance between sites of sampling associated with genetic variation among Neospora caninum in aborted bovine foetuses from northern Italy. Parasites Vectors 14, 47 (2021).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: