Skip to main content

Population genetic structure and phenotypic diversity of Aspidodera raillieti (Nematoda: Heterakoidea), a parasite of Didelphini marsupials in Brazil’s South and Southeast Atlantic Forest



The population genetics of parasites may be influenced by host specificity, life cycle, host geographical range, evolutionary history, and host population structure. The nematode Aspidodera raillieti infects different marsupial and rodent hosts in the Nearctic and Neotropical regions, implying a gene flow among populations. However, niche diversification of the main hosts of A. raillieti in superimposed areas may provide conditions for population genetic structuring within this parasite species. We examined the genetic structuring of A. raillieti infecting three marsupial species co-occurring along the South and Southeast Brazilian Atlantic Forest, a hotspot of biodiversity.


We employed morphometric analyses and partial mitochondrial cytochrome c oxidase I gene sequences (MT-CO1) to characterize populations via phylogenetic and phylogeographic analyses.


Among 175 A. raillieti specimens recovered from the marsupial hosts Didelphis aurita, D. albiventris, and Philander quica, we identified 99 MT-CO1 haplotypes forming four haplogroups and four clades in networks and phylogenetic trees, respectively. Clades I and II encompassed parasites of D. albiventris from the South region, clade III comprised parasites of D. aurita from the South and Southeast regions, and clade IV encompassed parasites of D. aurita and D. albiventris from the South and Southeast regions and parasites of P. quica from the South region. High genetic differentiation between clades, with a high fixation index and greater genetic variation in the analysis of molecular variance (AMOVA), indicated low gene flow between clades. Haplotypes shared among host species revealed a lack of host specificity. A significant correlation in the Mantel test suggested parasite isolation by distance, while there was no evidence of geographical structure between populations. Negative neutrality test values for clades III and IV suggested recent population expansion. Morphometric differentiation between A. raillieti specimens recovered from different host species, as well as from different localities, was more evident in males.


The genetic structure of A. raillieti populations in the South and Southeast Atlantic Forest resulted from historical events rather than from current geographical distribution or host specificity. We also demonstrate morphometric variation associated with host species and localities, suggesting phenotypic plasticity to host attributes and to spatial variables.

Graphical Abstract


The nematodes of the family Aspidoderidae Freitas, 1956, superfamily Heterakoidea Railliet & Henry, 1912, comprise Aspidodera Railliet & Henry, 1912; Paraspidodera Travassos, 1914; Lauroia Proença, 1938; and Nematomystes Sutton, Durette-Desset & Chabaud, 1980 [1,2,3]. The genus Aspidodera is morphologically characterized by cuticular expansion forming a cephalic hood with cordons arranged in six longitudinal loops at the anterior end, three lips, an oesophagus with a terminal bulb, a ventral sucker on males, a pair of spicules and a gubernaculum, and a posterior region ending in a digitiform projection [4, 5]. Aspidodera spp. are distributed along the southern Nearctic and the entire Neotropical region, parasitizing the cecum and the large intestines of mammals of the orders Didelphimorphia, Cingulata, Pilosa, and Rodentia [6].

Aspidodera raillieti Travassos, 1913 [7] is widely distributed in the Americas, occurring in the USA, Mexico, Guatemala, Costa Rica, Panama, Trinidad and Tobago, Colombia, Suriname, French Guiana, Brazil, Bolivia, Paraguay, and Argentina [2, 4, 5, 8]. Some species of Aspidodera were considered junior synonyms of A. raillieti because the morphological and morphometric differences used to propose these new species are not diagnostic characters and thus do not allow distinguishing them from A. raillieti. Among them are A. harwoodi Chandler, 1932 described from D. virginiana Kerr, 1792 in the USA; A. diazungriai Masf-Pallares and Vergara, 1971, parasitizing Didelphis pernigra J. A. Allen, 1900 in Paraguay; and A. vicentei Kohn, Fernandes and Mello, 1982, described from Nectomys squamipes Brants, 1827 in Brazil [5]. These findings suggest that the morphological and morphometric variability of A. raillieti might be associated with its geographical distribution and its hosts.

Aspidodera raillieti has been found in all Brazilian biomes in areas of Amazonia, Atlantic Forest, Caatinga, Cerrado, Pampa, and Pantanal. This nematode was found parasitizing different species of marsupials, including Didelphis aurita Wied-Neuwied, 1826; Didelphis albiventris Lund, 1840; Didelphis marsupialis Linnaeus, 1758; Philander quica (Temminck, 1824); Philander opossum (Linnaeus, 1758); Chironectes minimus (Zimmermann, 1780); Metachirus myosurus (É. Geoffroy St.-Hilaire, 1803); Caluromys lanatus (Olfers, 1818); Marmosops incanus (Lund, 1840); and a single rodent species, Nectomys squamipes (Brants, 1827) [9,10,11]. Four marsupial species of the tribe Didelphini occur in sympatry and are found infected by A. raillieti in areas of the Brazilian Atlantic Forest: D. aurita, D. albiventris, P. quica, and C. minimus [12]. The Atlantic Forest is among the most altered Brazilian biomes due to urban expansion and economic activities, retaining only approximately 12% of its original coverage [13,14,15,16].

Parasite populations may have their genetic structure influenced by several factors, such as the degree of host specificity, effective population size, geographical distance between populations, host dispersal ability, evolutionary history, host population structures, and life cycle complexity [16,17,18,19]. The low degree of host specificity exhibited by A. raillieti presupposes significant gene flow between populations of different host species, despite this nematode having a monoxenous life cycle [10]. Therefore, the A. raillieti population structure is expected to depend on host movements [17, 20,21,22], geographical distance, and/or historical processes [18].

A study on the ecological modelling of D. aurita and D. albiventris, based on climatic niches throughout their geographical distributions, showed that these marsupials might explore different niches in areas where they co-occur, mostly in areas of forest-grassland mosaics [23]. In addition, D. aurita may act as a biotic barrier for D. albiventris, where this habitat mosaic is not available for their coexistence [23]. These different host ecological characteristics may promote favourable conditions for the emergence of genetic structuring patterns in parasite populations, which is consistent with the use of different host species.

However, the distribution of P. quica overlaps with that of D. aurita and partially overlaps with that of D. albiventris [24, 25]. Philander quica and D. aurita have niche overlap and compete with each other [26, 27]; both species occur in humid forested areas [27], and predation upon P. quica by D. aurita occurs [27, 28], which may favour parasite transmission from one host species to another. These aspects may promote gene flow among the parasite populations of different hosts.

Nevertheless, another study demonstrated that populations of D. albiventris exhibit patterns of isolation by distance when comparing disjunction distributions between South and Southeast Brazil [29]. The same pattern was observed for D. aurita at a regional scale in Southeast Brazil [30, 31]. Likewise, this pattern may occur in parasite populations.

In this context, we aimed to examine the population genetic structure of the nematode A. raillieti, a parasite of the marsupials D. aurita, D. albiventris, and P. quica in South and Southeast Atlantic Forest localities. We hypothesized that A. raillieti populations are genetically structured as a function of their host species and/or geographical distances.


Host sampling

This study was conducted in eight localities in the Brazilian Atlantic Forest, from the state of Rio Grande do Sul in Brazil’s South Region to the state of Espírito Santo in the Southeast Region (Fig. 1, Tables 1 and 2). These localities include different natural forest formations, comprising dense ombrophilous forest, mixed ombrophilous forest, and semi-deciduous seasonal forest (Table 1) [14].

Fig. 1
figure 1

Sampling localities of A. raillieti recovered from marsupial hosts in southern and southeastern Brazil. The geometric shapes represent different host species, and shape colours indicate their respective municipalities. The background map displays elevation values in hypsometric tints

Table 1 Host species used in the present study, the total number of host individuals analysed, and specimens (N) per host species, with their respective locality abbreviations, municipalities and states, geographical coordinates, heights in metres above sea level (m.a.s.l.), and locality descriptions
Table 2 Distances between localities where hosts infected by A. raillieti were collected (in kilometres)

Marsupials were captured using Tomahawk Live Trap (Hazelhurst, WI) model 201 traps (16″ × 5″ × 5″) baited with a mixture of peanut butter, banana, oats, and bacon.

Helminth recovery

The digestive tract of marsupials was screened for parasites, and both the large intestine and the cecum were examined for the presence of specimens of A. raillieti. Organs were placed separately in Petri dishes, washed twice in physiological saline solution (NaCl 0.85%) to remove tissue debris, and stored in 70% ethanol solution. For examination, nematodes were clarified in 25% glycerin alcohol. Measurements and drawings were produced with the aid of a camera lucida attached to a Nikon Eclipse E200 MV R microscope. Specimens were randomly selected for measurements. The number of specimens used for morphometric analyses was based on the parasite burden (Table 3).

Table 3 Female, male, and total numbers of marsupial hosts and A. raillieti specimens measured from each locality and region

Discriminant analysis of principal components

We performed a discriminant analysis of principal components (DAPC) [32] to compare helminth morphometric differences considering localities and host species (except for Cariacica, Espírito Santo [CAR-ES] due to the low sample size of hosts and helminths recovered).

DAPC is a robust method for describing variations between defined groups, selecting principal components (PCs) that explain the greatest variation between groups while minimizing the variation within each group. We used the cross-validation optimization procedure to identify the ideal number of PCs to be retained by DAPC and selected the components associated with the lowest root mean squared error [33]. Finally, we determined the percentage of A. raillieti specimens correctly classified within their original group. Thus, we used DAPC to evaluate the results from the genetic analyses. DAPC was performed using the package ‘adegenet’ [34] in the R software environment, version 4.0.2 [35].

Genomic DNA isolation, amplification, and sequencing

We isolated genomic DNA from mid-section fragments of each adult specimen of A. raillieti using the QIAamp DNA Mini Kit (QIAGEN, Hilden, Germany) according to the manufacturer’s instructions. Before DNA isolation, each specimen was morphologically characterized and subsequently washed in distilled water. Amplifications of the mitochondrial cytochrome c oxidase subunit I gene (MT-CO1) by polymerase chain reaction (PCR) were performed individually for each A. raillieti specimen using the primer cocktail described by Prosser et al. [36]. Each reaction contained 12.5 μl of PCR Master Mix (Promega Corporation) (50 units/ml Taq DNA polymerase, 400 μM dATP, 400 μM dGTP, 400 μM dCTP, 400 μM dTTP, 3 mM MgCl2), 0.5 μl of each primer cocktail (10 μM of a three-forward primer mix and 10 μM of a three-reverse primer mix), 1 μl of genomic DNA, and ultrapure water in a total volume of 25 μl. Thermal cycling conditions were 94 °C for 1 min; five cycles at 94 °C for 40 s, 45 °C for 40 s, and 72 °C for 1 min; 35 cycles at 94 °C for 40 s, 51 °C for 40 s, and 72 °C for 1 min; and a final extension at 72 °C for 5 min [36]. Successfully amplified amplicons were purified using the Illustra GFX PCR DNA and Gel Band Purification Kit (GE Healthcare, Little Chalfont, Bucks, UK) following the manufacturer’s protocol and then cycle sequenced using the Big Dye Terminator v3.1 Cycle Sequencing Kit (Applied Biosystems, Carlsbad, CA, USA), individually for each cocktail primer for better accuracy. Cycle-sequenced product precipitation, formamide resuspension, and sequencing were conducted at the Capillary Sequencing (SANGER) Platform, P01-001-RPT/FIOCRUZ ( Sequencing was performed using the 96-capillary 3730xl DNA Analyzer (Applied Biosystems).

Molecular phylogenetic and phylogeographic analyses

For each specimen, DNA sequencing reads were assembled into contigs and edited for ambiguities using the Geneious version 9.1.8 bioinformatics software platform [37], resulting in a consensus sequence.

All sequences obtained in this study were deposited in GenBank (accession numbers OL676808–OL676982) (Additional file 1: Table S1). Two datasets were used in this study. The first dataset, used for phylogenetic inferences, included the sequences we generated and those available in GenBank (Additional file 1: Table S1), including sequences from other species of the family Aspidoderidae (Table 4). As an out-group, we added a sequence of Lauroia trinidadensis (Aspidoderidae). Since Nematomystes and A. raillieti formed a monophyletic group in previous studies [2, 8], we included Nematomystes spp. to the ingroup in our phylogenetic analyses. The second dataset, used for phylogeographic analyses, included only our 175 sequences of A. raillieti.

Table 4 Species, GenBank accession number, geographical locality, host, and references of Aspidoderidae GenBank sequences used in this study

At each dataset, we aligned the MT-CO1 sequences using the TranslatorX web server [38], employing amino acid translations to align protein-coding nucleotide sequences, using the MUSCLE algorithm [39]. The resulting alignments were manually trimmed of poorly aligned regions using the Mesquite software package, version 3.61 [40]. The presence of pseudogenes was checked using a phylogenetic method (PhyPA) based on pairwise alignment [41]. Substitution saturation in the matrices was assessed via the Xia test [42, 43]. Both tests were conducted using the DAMBE software package, version 6.4.79 [44].

Phylogenetic analyses under maximum likelihood (ML) as optimality criteria were generated using the PhyML, version 3.0 online web server [45]. Evolutionary model selection was implemented with SMS (Smart Model Selection) [46] in PhyML using the Akaike information criterion (AIC). Node support was assessed by the approximate likelihood-ratio test (aLRT) and by non-parametric bootstrap percentages (ML-BP) after 1000 replicates.

Bayesian inference (BI) analyses were performed using MrBayes, version 3.2.6 [47], executed on XSEDE through the CIPRES Science Gateway [48]. Independent GTR + I + G (general time-reversible nucleotide substitution model, with a proportion of invariable sites and gamma distribution of rates among sites) models were used for each codon position, with unlinking of base frequencies and parameters. Markov chain Monte Carlo (MCMC) sampling was performed for 10,000,000 generations with four simultaneous chains in two runs. Node support was assessed by Bayesian posterior probabilities (BPP), calculated from trees sampled every 100 generations, after removing the first 25% ‘burn-in’ generations. Sampling adequacy was assessed using the program Tracer, version 1.6 [49], to calculate the effective sample sizes (ESSs) of the parameters. Values above 200 effectively independent samples were considered robust.

Haplotype networks were inferred using the program PopART, version 1.7 [50], under the median-joining method [51]. We used DNAsp, version 5.10.1 [52], to organize A. raillieti sequences into groups according to (1) the clades recovered in the ML and BI phylogenetic trees, (2) the host species, and (3) the geographical localities. Additionally, using DNAsp, the genetic diversity of each group was calculated by the numbers of haplotypes (H), polymorphic sites (S), haplotype diversity (Hd), and nucleotide diversity (π).

Population genetic analyses

Analysis of molecular variance (AMOVA) [53] and fixation index (Fst) [54] were calculated using the Arlequin software package, version [55]. We used AMOVA to analyse genetic variability between and within previously defined groups and Fst to measure levels of genetic differentiation between groups. We also used the program Arlequin to assess deviation from neutrality using Tajima’s D [56] and Fu’s Fs [57] tests for each of the previously determined groups.

The Mantel test [58] was used to verify the correlation between genetic distances and geographical distances of A. raillieti specimens, as well as the correlation between genetic distances and elevation differences from all geographical localities studied. The genetic distance matrix was calculated in the package ‘ape’ [59] using the evolutionary model selected in the automated model selection feature of the program PAUP*, version 4.0a167 [60], under AICc. The geographical distance matrix was built from the geographical coordinates of the studied localities using the package ‘fields’ [61]. The elevation difference matrix was generated from the elevation of each locality using the package ‘vegan’ [62]. All procedures of the Mantel test were computed within the R software environment, version 4.0.2 [35]. All statistical analyses were performed at the 5% significance level.


Discriminant analysis of principal components

The number of PCs retained in the DAPC varied according to the investigated group. For the comparison between localities, 11 and 14 PCs were retained for female and male A. raillieti, respectively. For the comparison between host species, eight and 14 PCs were retained for female and male A. raillieti, respectively. The proportion of male specimens correctly classified to their original group was approximately 83% for locality and 75% for host species. The proportion of female specimens correctly classified was 75% for localities and 67% for host species.

We observed greater morphometric variation in parasites among localities than among host species (Fig. 2). Furthermore, this result was more evident among male A. raillieti, revealing three clusters, each cluster formed by specimens with high morphometric proximity. One of these clusters was formed by the locality of Porto Alegre, Rio Grande do Sul (POA-RS); the second was formed by the localities of São Gonçalo do Sapucaí, Minas Gerais (SGS-MG), Santo Amaro da Imperatriz, Santa Catarina (SAI-SC), Curitiba, Paraná (CUR-PR), and Petrópolis, Rio de Janeiro (PET-RJ); and the third was formed by Rio de Janeiro, Rio de Janeiro (RIO-RJ) and Paraty, Rio de Janeiro (PTY-RJ) (Fig. 2a).

Fig. 2
figure 2

a, c Population clusters for males and females. b, d Aspidodera raillieti specimens based on DAPC along with discriminant analysis (DA) eigenvalues, showing morphometric variations between a, b localities and c, d host species

Among the variables that best discriminated among locality groups for males were length and width of the sucker, length and width of the hood, length of the cordon, and length of the caudal spine (Fig. 3a). For females, the variables were cephalic hood length, cephalic hood width, cephalic cordon length, bulb length, and bulb width (Fig. 4a). The variables that best discriminated parasites between host species for males were length and width of the sucker, width of the hood, and length of the caudal spine (Fig. 3b). For females, these variables included the distance of the nerve ring to the anterior end, length and width of the bulb, and length of the oesophagus (Fig. 4b).

Fig. 3
figure 3

Morphometric variable contributions to A. raillieti male specimens. a Variables that best discriminated locality groups. b Variables that best discriminated host species groups

Fig. 4
figure 4

Morphometric variable contributions to A. raillieti female specimens. a Variables that best discriminated locality groups. b Variables that best discriminated host species groups

Phylogenetic analyses and genetic diversity

We successfully sequenced adult worms recovered from 17 D. aurita, nine D. albiventris, and four P. quica hosts collected in eight different localities. We obtained 489 bp MT-CO1 gene consensus sequences from 175 A. raillieti specimens. The first matrix, used for phylogenetic inferences, resulted in 183 taxa and 480 sites. The second matrix, used for phylogeographic analyses, included only our 175 sequences of A. raillieti with 489 sites. No long branch was observed in the PhyPA method, indicating the absence of pseudogenes. Xia’s test provided no evidence for substitution saturation in any data matrix.

The ML best-fit model chosen by SMS in PhyML under AIC was GTR + G, with four substitution rate categories and gamma shape parameter α = 0.095, resulting in a tree with lnL = −3030.323236 score. The BI sampling, after 25% ‘burn-in’, resulted in a mean estimated marginal likelihood of –3129.3479 (standard deviation = 18.2293; median = −3128.849). The ESSs were robust for all parameters.

Tree topologies from ML and BI analyses were similar, with some variations in nodes and support values. Aspidodera raillieti specimens formed a monophyletic group in both phylogenies (BP-ML = 0.42, aLRT = 1.00, BPP = 0.95). Within A. raillieti, four main monophyletic groups were recovered and identified as clades I, II, III, and IV (Fig. 5a, b).

Fig. 5
figure 5

a, b Phylogenetic trees of partial MT-CO1 sequences of A. raillieti from this study and of aspidoderid species from GenBank. The Lauroia trinidadensis sequence was added as an out-group. Clades of A. raillieti are collapsed: clade I pink, clade II green, clade III blue, clade IV yellow. a Bayesian inference topology. Node values are BPP. b Maximum likelihood topology. Node values are ML-BP and aLRT > 0.50 support. c Median-joining network of partial MT-CO1 sequence haplotypes of A. raillieti from this study. Circle sizes are proportional to haplotype frequencies, and colours represent the localities where each haplotype occurs. Lines circling groups of haplotypes delimit the four clades recovered in our phylogenetic analyses. Circles identified by numbers in the haplotype network represent haplotypes shared between localities. Each hatch mark along the lines connecting haplotypes in the median-joining networks represents one mutation

Clade I (BP-ML = 0.74, aLRT = 1.00, BPP = 1.00) comprised A. raillieti parasites of D. albiventris from POA-RS (Additional file 2: Fig. S1). Clade II (BP-ML = 0.93, aLRT = 1.00, BPP = 1.00) comprised parasites of D. albiventris from POA-RS (Additional file 3: Fig. S2). Clade III (BP-ML < 0.5, aLRT = 1.00, BPP = 1.00) comprised parasites of D. aurita from CAR-ES, RIO-RJ, PTY-RJ, and CUR-PR (Additional file 4: Fig. S3). Finally, clade IV (BP-ML = 0.66, aLRT = 1.00, BPP = 1.00) comprised parasites of D. aurita from RIO-RJ, PTY-RJ, PET-RJ, and CUR-PR; parasites of D. albiventris from SGS-MG, CUR-PR, and POA-RS; and parasites of P. quica from SAI-SC (Additional file 5: Fig. S4, Fig. 6).

Fig. 6
figure 6

Geographical distribution of the clades identified in this study and the number of haplotypes for each clade by studied localities

In the BI topology (Fig. 5a), clade I was recovered as sister to clades II, III, and IV (BPP = 0.95), whereas clades II and III formed a monophyletic group (BPP = 0.99), sister to clade IV (BPP = 0.66). In the ML topology (Fig. 5b), clades II and III formed a monophyletic group (aLRT = 1.00, BP-ML = 0.80), sister to a monophyletic group formed by clades I and IV (aLRT = 1.00, BP-ML > 0.5).

Haplotype networks

Among our 175 partial MT-CO1 gene sequences (489 bp) of A. raillieti, we identified 99 haplotypes with 114 polymorphic sites. These 99 haplotypes were grouped into four haplogroups, corresponding to the clades recovered in the phylogenies, separated by genetic distances of 20 to 30 mutational steps (Fig. 5c). The molecular diversity indices of groups were separated by the clades recovered in phylogenies, host species, and localities. All groups studied (clades, host species, and localities) had prominent levels of haplotype diversity but low levels of nucleotide diversity (Table 5).

Table 5 Genetic diversity indices for A. raillieti sequence groups divided by clades, host species, and localities

The localities RIO-RJ and PTY-RJ included clades III and IV haplotypes. The locality POA-RS included clades I, II, and IV haplotypes. The locality CUR-PR had only one clade III haplotype, while all others had clade IV haplotypes. Haplotypes from the localities CAR-ES (clade III), PET-RJ (clade IV), SGS-MG (clade IV) and SAI-SC (clade IV) clustered in only one haplogroup per locality (Fig. 5c). Some haplotypes were shared between localities. Haplotypes 9, 13, 21, and 79 were shared between the RIO-RJ and PTY-RJ localities. Haplotype 91 was shared between PET-RJ and SGS-MG. Haplotype 16 was shared between PTY-RJ and CUR-PR. Haplotype 64 was shared between CUR-PR and SAI-SC. Finally, haplotype 53 was shared between POA-RS and SAI-SC (Fig. 5c).

We also observed haplotype sharing between host species. Haplotype 91 was shared between D. aurita (PET-RJ) and D. albiventris (SGS-MG). Haplotype 64 was shared between D. aurita (CUR-PR) and P. quica (SAI-SC). Haplotype 53 was shared between D. albiventris (POA-RS) and P. quica (SAI-SC) (Fig. 7).

Fig. 7
figure 7

Median-joining network for partial MT-CO1 sequence haplotypes of 175 A. raillieti sequences from the present study. Circle sizes are proportional to haplotype frequencies, and the colours represent the hosts (D. aurita, D. albiventris, and P. quica) in which each haplotype occurs. Each hatch mark along the lines connecting haplotypes represents one mutation. Only haplotypes shared between hosts have their identifications represented in the haplotype network

Population structure

The AMOVA result for clade groups revealed higher variation among clades, which represented 80.64% of the total variation, while the variation within clades represented 19.36%. The genetic variation for the host species groups revealed less variation among host species, which represented 28.50% of the variation, while the variation within host species represented 71.50%. The genetic variation for locality groups revealed higher variation among localities, which represented 57.52% of the variation, while the variation within localities represented 42.48% (Table 6).

Table 6 Analysis of molecular variance (AMOVA) of A. raillieti among and within clades recovered in phylogenetic inferences, host species, and geographical localities

Considering the fixation index, all Fst values revealed significant genetic differences between clades (P < 0.05), host species (P < 0.01) (Table 7), and locality groups (P < 0.05). No difference was observed between RIO-RJ and CAR-ES or between PTY-RJ and CAR-ES (Table 8). The Mantel tests executed using the Kimura two-parameter nucleotide substitution model, calculated using PAUP, indicated a significant and positive correlation between genetic distance and geographical distance (N = 175, r = 0.27, P = 0.001) (Fig. 8a) and a significant and positive correlation between genetic distance and elevation difference (N = 175, r = 0.32, P = 0.001) (Fig. 8b).

Table 7 Pairwise Fst values of A. raillieti sequences between clades and host species groups
Table 8 Pairwise Fst values of A. raillieti sequences between geographical locality groups
Fig. 8
figure 8

a Mantel test of correlation between genetic divergences and geographical distances; b Mantel test of correlation between genetic divergences and elevation differences

Historical demography

In Tajima’s D and Fu’s Fs neutrality tests, calculated for the clades, D values (−2.02) were significant only for clade III (P = 0.005), while Fs values were significant for both clade III (−15.09) and clade IV (−24.39) (P = 0.000). For the host species groups, the D and Fs values were significant only for the group of A. raillieti parasites of P. quica (−1.545, P = 0.05 and −8.333, P = 0.0004, respectively). For the locality groups, D values were significant only for the localities PTY-RJ (−1.819, P = 0.024), PET-RJ (−1.654, P = 0.027) and CUR-PR (−1.746, P = 0.015), whereas Fs values were significant for the localities CUR-PR (−6.942, P = 0.007), SGS-MG (−6.931, P = 0.000), and SAI-SC (−8.333, P = 0.003).


Morphometric comparison between A. raillieti specimens from different localities and hosts.

The DAPC for A. raillieti specimen groups associated with different hosts (D. albiventris, D. aurita, and P. quica) indicated morphometric differentiation for both females and males, however, with partial overlapping of some specimens from different host species.

Among the variables that best discriminated A. raillieti specimens associated with host species and localities, our results identified the cephalic hood width for male and female specimens, the sucker length and width for females, and the bulb length and width for females. This morphometric variability may be a consequence of adaptations to environmental conditions, as observed in the trematode Echinostoma paraensei Lie & Basch, 1967 [63].

When analysing DAPC morphometric differentiation between A. raillieti specimen groups associated with different localities, female nematodes had greater overlap between localities, while males were differentiated into three groups. The first group was composed of male specimens from the localities RIO-RJ and PTY-RJ, the second was composed of specimens from POA-RS, and the third was composed of specimens from PET-RJ, SGS-MG, CUR-PR, and SAI-SC.

Our findings indicated the influence of host species, as well as locality, on the morphometry of A. raillieti, thus suggesting phenotypic plasticity regarding host attributes and spatial variables [64].

Population structure of A. raillieti

Population structure can be influenced by different evolutionary forces. Among them, gene flow is considered of fundamental importance, as it allows the exchange of genetic information between populations, homogenizing the variation among them [65]. Population genetic studies of parasitic helminths have shown that gene exchange between populations is strongly influenced by the movement of their vertebrate [17, 22] and invertebrate hosts [66]. Aspidodera raillieti is a parasite with a monoxenous life cycle [10], depending on its mammalian hosts for dispersion. Thus, its population structure is expected to be dependent on the movement and encounter of mammalian hosts.

Our phylogenetic trees and networks recovered four clades of A. raillieti. These results provide evidence that these lineages experienced past events that contributed to the genetic divergence observed between clades, since the genetic structuring observed in these helminths was not related to geographical distance, as evidenced by the presence of divergent clades in the same locality.

The genetic structure observed using AMOVA indicated a greater genetic variation among clades (interpopulation) than within clades (intrapopulation). The highly significant Fst values, indicating high genetic differentiation between clades [67], concurred with the AMOVA results, also indicating limited gene flow between the four clades. However, analysing other genetic markers from independent loci would be necessary to deliver a clearer picture of the evolutionary history of A. raillieti.

No geographical structuring was observed, since the AMOVA results had similar percentages of genetic variation among and within localities (57.52% and 42.48%, respectively). This was congruent with Fst, with high values both between geographically distant localities and between closer localities. Nevertheless, the Mantel test showed a significant and positive correlation between genetic and geographical distances, indicating isolation by distance (IBD). The disagreement between the population structure analyses and the Mantel test may be a consequence of the co-occurrence of haplotypes from different clades in the same locality, possibly due to historical processes. The phylogeographic patterns of species can be affected by different factors, such as dispersal or vicariant events, which can promote differentiation between populations [65]. In addition, it has been postulated that the current distributions of several lineages of mammals, birds, and amphibians in South America originated from several mechanisms, such as Quaternary climatic oscillations and Tertiary orogenic events [68, 69]. However, to identify which event led to the divergence between the lineages recovered in our analyses, a well-calibrated molecular clock would be necessary to estimate the divergence times between them compared to known events [65].

The high genetic diversity observed in A. raillieti has also been identified in populations of the nematode Heligmosomoides polygyrus (Dujardin, 1845), a parasite of the forest rodent Apodemus sylvaticus (Linnaeus, 1758). To understand the phylogeographic pattern of H. polygyrus from different localities in Europe, Nieberding et al. [70] studied the genetic structure of the populations of this monoxenous nematode. The authors found a high number of haplotypes for the MT-CYB gene, totalling 126 haplotypes from 136 sequenced specimens. Five main groups were observed, both in phylogenetic reconstructions and haplotype networks, which showed a high degree of genetic divergence, being separated by a genetic distance of 18 to 35 mutational steps, as observed between A. raillieti clades (20 to 30 mutational steps), with some haplotypes co-occurring in some localities.

Similarly, the infective free-living stages of H. polygyrus and A. raillieti have no specialized structures for dispersal; thus, gene flow between populations depends mainly on host movements and social behaviour. However, H. polygurus is a host specialist and has a wide distribution congruent with its host, partially reflecting the phylogeographic history of its host [70, 71]. In contrast, A. raillieti is a host-generalist parasite that is able to infect several marsupials and one rodent species, which makes its phylogeographic patterns potentially more complex, requiring studies with broader geographical ranges and other host species to better understand its genetic population structure.

Analysing A. raillieti by host species, the phylogenetic and phylogeographic results showed no evidence of population specificity for host species, unlike the morphometric analyses, since nematode haplotypes from different marsupial species were in the same clade. Additionally, some haplotypes were shared between different host species. Corroborating this pattern, no structuring was observed in A. raillieti associated with host species, as AMOVA indicated low genetic variation between the specimens recovered from each host species (interhost). The significant Fst values between host species indicated moderate genetic differentiation. Moreover, this genetic differentiation may result partially from geographical distances among localities and from divergent clades.

Some ecological characteristics of the host species studied promote unfavourable conditions for the emergence of parasite population genetic structuring patterns. Didelphis aurita and D. albiventris are omnivorous frugivores, while P. quica is an omnivorous insectivore, and all have overlapping diets [72]. They also have the same locomotor habits, as both are scansorial [73, 74]. However, D. albiventris is a habitat generalist, while the other two species mostly occur in forested regions [75, 76]. In addition, both D. aurita and D. albiventris occur in abundance in degraded areas [77], unlike P. quica [29]. Although Cáceres et al. [23] have shown that in areas of sympatry, D. aurita and D. albiventris explore different niches, this barrier may not have been sufficient to prevent gene flow between populations of A. raillieti in these hosts. We also observed a significant correlation between genetic divergence and elevation differences. The geographical distribution of D. albiventris is larger than that of D. aurita, as the former is a more habitat generalist and has greater climatic tolerance than the latter, which is also reflected in the elevation.

Moreover, we expected to find greater genetic differentiation between parasites recovered from P. quica and D. albiventris than between parasites recovered from P. quica and D. aurita, as the distribution of P. quica overlaps that of D. aurita but not completely that of D. albiventris [12, 29]. As D. aurita and P. quica have niche overlap, compete [27, 28], and may have intraguild predation [30], all these characteristics may favour parasite gene flow between these host populations. However, less genetic differentiation was found between A. raillieti specimens recovered from P. quica and D. albiventris than between A. raillieti specimens recovered from P. quica and D. aurita. This may be because A. raillieti samples of P. quica were collected from a single locality, closer to D. albiventris than to D. aurita localities.

As observed in A. raillieti populations, the nematode Trichostrongylus axei (Cobbold, 1879), a host-generalist parasite that infects multiple sympatric wild ungulates, showed no evidence of genetic structure associated with host species [20]. The authors proposed that T. axei populations would be structured due to the degree of spatial niche partitioning between hosts.

López-Caballero et al. [78] performed a study on the genetic divergence of populations of the acanthocephalan Oligacanthorhynchus microcephalus (Rudolphi, 1819) parasitizing three definitive hosts of the tribes Didelphini, Didelphis marsupialis, D. virginiana, and Philander opossum from different localities in Mexico. Phylogenetic analyses demonstrated a similar pattern to that found for A. raillieti, in which the specimens of O. microcephalus were grouped into three main clades, which were not correlated either with definitive host species or with geographical distributions. The lack of population structuring was attributed to several aspects, including host natural histories, dispersal abilities, sympatries, overlapping diets, and the fact that the arthropod intermediate hosts of O. microcephalus are distributed throughout the entire geographical range of this parasite.

Demographic history of Aspidodera raillieti

Climate changes, such as Pleistocene glaciations, promoted the retraction of tropical forests, forming refuges and the subsequent expansion of these forests due to climate amelioration. Populations from refuge areas that undergo postglacial demographic expansion have consequent genetic signatures [79, 80]. These demographic fluctuations can be detected by some analyses, such as the neutrality tests used in this study [56, 57]. Significant negative D or Fs values in neutrality tests suggest a population undergoing purifying selection or expansion, characterized by an excess of rare alleles [67].

In our neutrality tests, parasitic specimens from the PTY-RJ, PET-RJ, CUR-PR, SGS-MG, and SAI-SC localities had significant negative values, some of which were negative for Tajima’s D or Fu’s Fs. This expansion signature was congruent with the genetic diversity indices, showing high haplotype diversity and low nucleotide diversity for the MT-CO1 gene for all groups studied (clades, localities, and hosts). These results also showed that although there were many haplotypes, they differed from each other by only a few nucleotide substitutions. This pattern is consistent with a rapidly expanding population from a small effective population size [81].

Since A. raillieti has a Neotropical and partially Nearctic distribution, future studies should include specimens from other biomes, encompassing its entire distribution range, to better understand the evolutionary history of this parasite. It would also be necessary to include nuclear genetic markers from independent loci to verify whether the phylogeographic pattern observed for the MT-CO1 gene is corroborated. Additionally, the inclusion of a time scale to estimate divergence times between clades would make it possible to verify congruence between cladogenesis and palaeogeographical and climatic events [71, 82].


Based on our results, we concluded that the genetic structure of A. raillieti populations in the South and Southeast Atlantic Forest was likely associated with historical events, such as past climate changes, and not with the host species D. aurita, D. albiventris, and P. quica or with the current geographical distribution of this parasitic nematode. We also observed greater morphometric variation than molecular structuring associated with host species and localities, suggesting phenotypic plasticity related to host functional traits, as well as to spatial variables.

Availability of data and materials

All the new A. raillieti MT-CO1 sequences obtained for this study were deposited in NCBI GenBank under accession numbers OL676808-OL676982.



Meters above sea level


Cariacica, Espírito Santo


Petrópolis, Rio de Janeiro


Rio de Janeiro, Rio de Janeiro


Paraty, Rio de Janeiro


São Gonçalo do Sapucaí, Minas Gerais


Curitiba, Paraná


Santo Amaro da Imperatriz, Santa Catarina


Porto Alegre, Rio Grande do Sul

NaCl 0.85%:

Physiological saline solution


Discriminant analysis of principal components


Principal components


Mitochondrially encoded cytochrome c oxidase I


Polymerase chain reaction


Phylogenetic method


Maximum likelihood


Bayesian inference


Smart Model Selection


Akaike information criterion


Approximate likelihood-ratio test


Bootstrap percentage


General time-reversible nucleotide substitution model


proportion of invariable sites


gamma-distributed heterogeneity of substitution rates among sites


Markov chain Monte Carlo


Bayesian posterior probabilities


Numbers of haplotypes


Polymorphic sites


Haplotype diversity


Nucleotide diversity


Analysis of molecular variance

F st :

Fixation index


Isolation by distance


Corrected Akaike information criterion


Effective Sample Size


  1. Chagas-Moutinho VA, Sant’Anna V, Oliveira-Menezes A, De Souza W. New Aspidoderidae species parasite of Didelphis aurita (Mammalia: Didelphidae): A light and scanning electron microscopy approach. Acta Trop. 2014;130:162–6.

    Article  CAS  PubMed  Google Scholar 

  2. Jiménez FA, Carreno RA, Gardner SL. Aspidodera kinsellai n. sp. (Nematoda: Heterakoidea) from nine-banded armadillos in Middle America with notes on phylogeny and host-parasite biogeography. J Parasitol. 2013;99:1056–61.

    Article  PubMed  Google Scholar 

  3. Jiménez-Ruiz FA, Gardner SL, Noronha D, Pinto RM. The systematic position of Lauroiinae Skrjabin and Schikhobalova, 1951 (Nemata: Heterakoidea: Aspidoderidae), as revealed by the analysis of traits used in its diagnosis. Cladistics. 2008;24:459–76.

    Article  PubMed  Google Scholar 

  4. Chagas-Moutinho VA, Oliveira-Menezes A, Cárdenas MQ, Lanfredi RM. Further description of Aspidodera raillieti (Nematoda: Aspidoderidae) from Didelphis marsupialis (Mammalia: Didelphidae) by light and scanning electron microscopy. Parasitol Res. 2007;101:1331–6.

    Article  CAS  PubMed  Google Scholar 

  5. Santos CP, Lent H, Gomes DC. The genus Aspidodera Railliet and Henry.1912 (Nematoda: Heterakoidea): revision, new synonyms and key for species. Rev Bras Biol. 1990;50:1017–31.

    Google Scholar 

  6. Jiménez-Ruiz FA, Gardner SL, Varela-Stokes AS. Aspidoderidae from North America, with the description of new species of Aspidodera (Nematoda: Heterakoidea). J Parasitol. 2006;92:847–54.

    Article  PubMed  Google Scholar 

  7. Travassos L. Sobre as espécies brasileiras da subfamília Heterakinae Railliet & Henry, 1912. Mem Inst Oswaldo Cruz. 1913;5:271–318.

    Article  Google Scholar 

  8. Jiménez FA, Gardner SL, Navone G, Ortí G. Four events of host switching in Aspidoderidae (Nematoda) involve convergent lineages of mammals. J Parasitol. 2012;98:1166–75.

    Article  PubMed  Google Scholar 

  9. Cirino BS, Costa SFD, Maldonado A, Gentile R. First study on the helminth community structure of the neotropical marsupial Metachirus myosuros (Didelphimorphia, Didelphidae). Rev Bras Parasitol Vet. 2020;29:e005420.

    Article  CAS  PubMed  Google Scholar 

  10. Jiménez FA, Catzeflis F, Gardner SL. Structure of parasite component communities of didelphis marsupials: insights from a comparative study. J Parasitol. 2011;97:779–87.

    Article  PubMed  Google Scholar 

  11. Vicente JJ, Gomes DC, Araujo-Filho NA. Alguns helmintos de marsupiais e roedores da Ilha Grande, Estado do Rio de Janeiro. Atas da Sociedade de Biologia do Rio de Janeiro, Rio de Janeiro. 1982;26:3–4.

    Google Scholar 

  12. Gardner AL, editor. Mammals of South America, volume 1: marsupials, xenarthrans, shrews, and bats, vol. 2. Chicago: University of Chicago Press; 2008.

    Google Scholar 

  13. Colombo AF, Joly CA. Brazilian Atlantic Forest lato sensu: the most ancient Brazilian forest, and a biodiversity hotspot, is highly threatened by climate change. Braz J Biol. 2010;70:697–708.

    Article  CAS  PubMed  Google Scholar 

  14. Marques MC, Grelle CE, editors. The Atlantic Forest: History, biodiversity, threats and opportunities of the mega-diverse forest. Berlin: Springer Nature; 2021.

    Google Scholar 

  15. Mittermeier RA, Grill PR, Hoffman M, Pilgrim J. Hotspots revisited: Earth’s biologically richest and most endangered terrestrial ecoregions Cemex. Mexico City. 2004.

  16. Myers N, Mittermeier RA, Mittermeier CG, Da Fonseca GA, Kent J. Biodiversity hotspots for conservation priorities. Nature. 2000;403:853–8.

    Article  CAS  PubMed  Google Scholar 

  17. Blouin MS, Yowell CA, Courtney CH, Dame JB. Host movement and the genetic structure of populations of parasitic nematodes. Genetics. 1995;141:1007–14.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  18. Gorton MJ, Kasl EL, Detwiler JT, Criscione CD. Testing local-scale panmixia provides insights into the cryptic ecology, evolution, and epidemiology of metazoan animal parasites. Parasitology. 2012;139:981–97.

    Article  PubMed  Google Scholar 

  19. Huyse T, Poulin R, Theron A. Speciation in parasites: a population genetics approach. Trends Parasitol. 2005;21:469–75.

    Article  PubMed  Google Scholar 

  20. Archie EA, Ezenwa VO. Population genetic structure and history of a generalist parasite infecting multiple sympatric host species. Int J Parasitol. 2011;41:89–98.

    Article  PubMed  Google Scholar 

  21. Barrett LG, Thrall PH, Burdon JJ, Linde CC. Life history determines genetic structure and evolutionary potential of host-parasite interactions. Trends Ecol Evol. 2008;23:678–85.

    Article  PubMed  PubMed Central  Google Scholar 

  22. Criscione CD, Poulin R, Blouin MS. Molecular ecology of parasites: elucidating ecological and microevolutionary processes. Mol Ecol. 2005;14:2247–57.

    Article  CAS  PubMed  Google Scholar 

  23. Cáceres NC, de Moraes WM, Melo GL, Meloro C, Sponchiado J, Carvalho RDS, et al. Which factors determine spatial segregation in the South American opossums (Didelphis aurita and D. albiventris)? An ecological niche modelling and geometric morphometrics approach. PLoS ONE. 2016;11:e0157723.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  24. Santori RT, Astúa D, Finotti R, Lessa LG, Cerqueira R. Inter and intraspecific differences in food resources selection and use in captive Philander frenatus and Didelphis aurita (Didelphimorphia; Didelphidae). Oecol Aust. 2015;19:102.

    Article  Google Scholar 

  25. Ceotto P, Finotti R, Santori R, Cerqueira R. Diet variation of the marsupials Didelphis aurita and Philander frenatus (Didelphimorphia, Didelphidae) in a rural area of Rio de Janeiro state, Brazil. Mastozool Neotrop. 2009;16:49–58.

    Google Scholar 

  26. Voss RS, Díaz-Nieto JF, Jansa SA. A revision of Philander (Marsupialia: Didelphidae), part 1: P. quica, P. canus, and a new species from Amazonia. Am Mus Novit. 2018;3891:1–70.

    Article  Google Scholar 

  27. Moura MC, Vieira MV, Cerqueira R. Occasional intraguild predation structuring small mammal assemblages: the marsupial Didelphis aurita in the Atlantic Forest of Brazil. Austral Ecol. 2009;34:481–9.

    Article  Google Scholar 

  28. Wilson DE. Opossum predation: Didelphis on Philander. J Mammal. 1970;51:386–7.

    Article  Google Scholar 

  29. Sousa LC, Gontijo CM, Botelho HA, Fonseca CG. Mitochondrial genetic variability of Didelphis albiventris (Didelphimorphia, Didelphidae) in Brazilian localities. Genet Mol Biol. 2012;35:522–9.

    Article  PubMed  PubMed Central  Google Scholar 

  30. Gentile R, Kajin M. Long-term empirical studies of small mammals: the contribution of Dr. Rui Cerqueira to population biology. Oecol Aust. 2015;19:1–15.

    Article  Google Scholar 

  31. Gentile R. Biologia de populações de Didelphis aurita (Didelphimorphia, Didelphidae): Dinâmica e variabilidade genética de populações na Floresta Atlântica do Estado do Rio de Janeiro. PhD thesis, Universidade Federal do Rio de Janeiro, Rio de Janeiro, Brasil, 2000.

  32. Jombart T, Devillard S, Balloux F. Discriminant analysis of principal components: a new method for the analysis of genetically structured populations. BMC Genet. 2010;11:94.

    Article  PubMed  PubMed Central  Google Scholar 

  33. Jombart T, Collins C. A tutorial for discriminant analysis of principal components (DAPC) using adegenet, Rvignette. 2015;1–37.

  34. Jombart T. adegenet: a R package for the multivariate analysis of genetic markers. Bioinform. 2008;24:1403–5.

    Article  CAS  Google Scholar 

  35. R Core Team (2020). R: A language and environment for statistical computing. R Foundation Statistical Computing. Vienna, Austria. 2020.

  36. Prosser SW, Velarde-Aguilar MG, León-Règagnon V, Hebert PD. Advancing nematode barcoding: a primer cocktail for the cytochrome C oxidase subunit I gene from vertebrate parasite nematodes. Mol Ecol Resour. 2013;13:1108–15.

    CAS  PubMed  Google Scholar 

  37. Kearse M, Moir R, Wilson A, Stones-Havas S, Cheung M, Sturrock S, et al. Geneious Basic: an integrated and extendable desktop software platform for the organization and analysis of sequence data. Bioinform. 2012;28:1647–9.

    Article  Google Scholar 

  38. Abascal F, Zardoya R, Telford MJ. TranslatorX: multiple alignment of nucleotide sequences guided by amino acid translations. Nucleic Acids Res. 2010;38:W7–13.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  39. Edgar RC. MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 2004;32:1792–7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  40. Maddison WP, Maddison DR. Mesquite: a modular system for evolutionary analysis, 2021.

  41. Xia X. DAMBE6: new tools for microbial genomics, phylogenetics, and molecular evolution. J Hered. 2017;108:431–7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  42. Xia X, Lemey P. Assessing substitution saturation with DAMBE. Phylogenet Handb. 2009;2:615–30.

    Article  Google Scholar 

  43. Xia X, Xie Z, Salemi M, Chen L, Wang Y. An index of substitution saturation and its application. Mol Phylogenet Evol. 2003;26:1–7.

    Article  CAS  PubMed  Google Scholar 

  44. Xia X, Xie Z. DAMBE: software package for data analysis in molecular biology and evolution. J Hered. 2001;92:371–3.

    Article  CAS  PubMed  Google Scholar 

  45. Guindon S, Dufayard JF, Lefort V, Anisimova M, Hordijk W, Gascuel O. New algorithms and methods to estimate maximum-likelihood phylogenies: assessing the performance of PhyML 3.0. Syst Biol. 2010;59:307–21.

    Article  CAS  PubMed  Google Scholar 

  46. Lefort V, Longueville JE, Gascuel O. SMS: smart model selection in PhyML. Mol Biol Evol. 2017;34:2422–4.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  47. Ronquist F, Teslenko M, Van Der Mark P, Ayres DL, Darling A, Höhna S, et al. MrBayes 3.2: efficient Bayesian phylogenetic inference and model choice across a large model space. Syst Biol. 2012;61:539–42.

    Article  PubMed  PubMed Central  Google Scholar 

  48. Miller MA, Pfeiffer W, Schwartz T. (2010) "Creating the CIPRES Science Gateway for inference of large phylogenetic trees" in Proceedings of the Gateway Computing Environments Workshop (GCE), 14 Nov. 2010, New Orleans, LA pp 1–8.

  49. Rambaut A, Suchard MA, Xie D, Drummond AJ. Tracer v1. 6. Computer program and documentation distributed by the author. 2014.

  50. Leigh JW, Bryant D. PopART: Full-feature software for haplotype network construction. Methods Ecol Evol. 2015;6:1110–6.

    Article  Google Scholar 

  51. Bandelt H, Forster P, Röhl A. Median-joining networks for inferring intraspecific phylogenies. Mol Biol Evol. 1999;16:37–48.

    Article  CAS  PubMed  Google Scholar 

  52. Librado P, Rozas J. DnaSP v5: a software for comprehensive analysis of DNA polymorphism data. Bioinform. 2009;25:1451–2.

    Article  CAS  Google Scholar 

  53. Excoffier L, Smouse PE, Quattro JM. Analysis of molecular variance inferred from metric distances among DNA haplotypes: application to human mitochondrial DNA restriction data. Genetics. 1992;131:479–91.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  54. Wright S. Evolution and the genetics of populations, vol. 4. Variability within and among natural populations. 1978.

  55. Excoffier L, Laval G, Schneider S. Arlequin (version 3.0): an integrated software package for population genetics data analysis. Evolutionary bioinformatics online. 2005;1:117693430500100003.

  56. Tajima F. Statistical method for testing the neutral mutation hypothesis by DNA polymorphism. Genetics. 1989;123:585–95.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  57. Fu Y. Statistical tests of neutrality of mutations against population growth, hitchhiking and background selection. Genetics. 1997;147:915–25.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  58. Mantel N. The detection of disease clustering and a generalized regression approach. Cancer Res. 1967;27:209–20.

    CAS  PubMed  Google Scholar 

  59. Paradis E, Claude J, Strimmer K. APE: analyses of phylogenetics and evolution in R language. Bioinform. 2004;20:289–90.

    Article  CAS  Google Scholar 

  60. Swofford DL. PAUP*. Phylogenetic analysis using parsimony (*and other methods). Version 4. Sunderland: Sinauer Associates; 2003.

    Google Scholar 

  61. Nychka D, Furrer R, Sain S. (2017). Fields: Tools for Spatial Data (R Package).

  62. Oksanen J, Guillaume Blanchet F, Friendly M, Kindt R, Legendre P, McGlinn D, et al. vegan: Community Ecology Package. R package version. 2020;2.5–7.

  63. Júnior AM, Zeitone BK, Amado LA, Amado IF, Machado-Silva JR, Lanfredi RM. Biological variation between two Brazilian geographical isolates of Echinostoma paraensei. J Helminthol. 2005;79:345–51.

    Article  Google Scholar 

  64. Cardoso TS, Andreazzi CS, Junior AM, Gentile R. Functional traits shape small mammal-helminth network: patterns and processes in species interactions. Parasitology. 2021;148:947–55.

    Article  PubMed  Google Scholar 

  65. Freeland JR, Kirk H, Petersen S. Molecular ecology. 2nd ed. Hoboken: Wiley-Blackwell; 2011.

    Book  Google Scholar 

  66. Peng J, He ZP, Zhang S, Lun ZR, Wu ZD, Fan CK. Phylogeography of Angiostrongylus cantonensis (Nematoda: Angiostrongylidae) in southern China and some surrounding areas. PLoS Negl Trop Dis. 2017;11:e0005776.

    Article  PubMed  PubMed Central  Google Scholar 

  67. Hartl DL, Clark AG. Princípios de Genética de Populações. 4th ed. Porto Alegre: Artmed; 2010.

    Google Scholar 

  68. Rull V. Speciation timing and neotropical biodiversity: the Tertiary-Quaternary debate in the light of molecular phylogenetic evidence. Mol Ecol. 2008;17:2722–9.

    Article  PubMed  Google Scholar 

  69. Turchetto-Zolet AC, Pinheiro F, Salgueiro F, Palma-Silva C. Phylogeographical patterns shed light on evolutionary process in South A merica. Mol Ecol. 2013;22:1193–213.

    Article  CAS  PubMed  Google Scholar 

  70. Nieberding C, Libois R, Douady CJ, Morand S, Michaux JR. Phylogeography of a nematode (Heligmosomoides polygyrus) in the western Palearctic region: persistence of northern cryptic populations during ice ages? Mol Ecol. 2005;14:765–79.

    Article  CAS  PubMed  Google Scholar 

  71. Nieberding C, Morand S, Libois R, Michaux JR. A parasite reveals cryptic phylogeographic history of its host. Proc R Soc Lond B Biol Sci. Series B: Biological Sciences. 2004;271:2559–2568.

  72. Santori RT, Astúa de Moraes D, Grelle CE, Cerqueira R. Natural diet at a restinga forest and laboratory food preferences of the opossum Philander frenata in Brazil. Stud Neotrop Fauna Environ. 1997;32:12–6.

    Article  Google Scholar 

  73. Cáceres NC, Araujo Monteiro-Filho EL. The common opossum, Didelphis aurita, as a seed disperser of several plants in southern Brazil. Ciênc Cult. 2000;52:41–4.

    Google Scholar 

  74. Cunha AA, Vieira MV. Support diameter, incline, and vertical movements of four didelphid marsupials in the Atlantic forest of Brazil. J Zool. 2002;258:419–26.

    Article  Google Scholar 

  75. Cáceres NC, Machado AF. Spatial, dietary and temporal niche dimensions in ecological segregation of two sympatric, congeneric marsupial species. Open Ecol J. 2013;6(1).

  76. Cáceres NC. Food habits and seed dispersal by the white-eared opossum, Didelphis albiventris, in southern Brazil. Stud Neotrop Fauna Environ. 2002;37:97–104.

    Article  Google Scholar 

  77. Cerqueira R. The distribution of Didelphis in South America (Polyprotodontia, Didelphidae). J Biogeogr. 1985;135–145.

  78. López-Caballero J, Mata-López R, García-Varela M, Pérez-Ponce de León G. Genetic divergence of Oligacanthorhynchus microcephalus (Acanthocephala: Archiacanthocephala: Oligacanthorhynchidae), parasite of three species of opossum (Mammalia: Didelphidae) across Central and Southeastern Mexico. Comp Parasitol. 2015;82:175–86.

    Article  Google Scholar 

  79. Hewitt G. The genetic legacy of the Quaternary ice ages. Nature. 2000;405:907–13.

    Article  CAS  PubMed  Google Scholar 

  80. Lessa EP, Cook JA, Patton JL. Genetic footprints of demographic expansion in North America, but not Amazonia, during the Late Quaternary. PNAS. 2013;100:10331–4.

    Article  CAS  Google Scholar 

  81. Avise JC. Phylogeography: the history and formation of species. Harvard: Harvard University Press; 2000.

    Book  Google Scholar 

  82. Li L, Lü L, Nadler SA, Gibson DI, Zhang LP, Chen HX, et al. Molecular phylogeny and dating reveal a terrestrial origin in the early carboniferous for ascaridoid nematodes. Syst Biol. 2018;67:888–900.

    Article  PubMed  Google Scholar 

  83. Sikes RS. Guidelines of the American Society of Mammalogists for the use of wild mammals in research and education. J Mammal. 2016;2016:663–88.

    Article  Google Scholar 

  84. CFMV, 2013. Guia Brasileiro de Boas Práticas em Eutanásia em Animais - Conceitos e Procedimentos Recomendados. CFMV, Conselho Federal de Medicina Veterinária, Brasília, DF.

Download references


We would like to thank the staff and students of the Laboratory of Biology and Parasitology of Reservoir Wild Mammals (LABPMR/IOC/FIOCRUZ) who helped in the field and lab work. We are grateful to Dr Paulo D’Andrea for providing his ICMBio licence, Dr Rui Cerqueira for his FATMA and INEA authorizations, and Dr Gilberto Gazêta for his IEMA authorization. We are indebted to Fiocruz Atlantic Forest Campus (CFMA), State Institute of the Environment (INEA), Serra dos Órgãos National Park (PARNASO), Pedra Branca State Park (PEPB), Serra do Tabuleiro State Park (PEST), and Duas Bocas Biological Reserve (RBDB) for providing us permission to collect and local facilities. We especially thank the DNA Sequencing Platform of the Oswaldo Cruz Institute (RPT/FIOCRUZ) for assistance with sample sequencing, Ricardo Baptista Schmidt for imaging support, and Fernando de Oliveira dos Santos for helping with the maps.


This project was financially supported by the National Council for Scientific and Technological Development (CNPq) PPBioMA #457524/2012-0, CNPq #306352/2014-1, by the Foundation Carlos Chagas Filho Research Support of the State of Rio de Janeiro (FAPERJ) #E-26/111.296/2014, by the Oswaldo Cruz Institute (IOC/FIOCRUZ) #, and by the Foundation Carlos Chagas Filho Research Support of the State of Rio de Janeiro—FAPERJ (APQ1 E-26/010.001597/2019). RG received research grants from CNPq (304355/2018-6). TSC received grants from the Program for the Innovation Promotion (Inova Fiocruz), and KV received grants provided by the CNPq.

Author information

Authors and Affiliations



KV, AMJ, and RVV co-designed the study. SFC, RG, and RVV collected the biological material. KV, AMJ, and SFC identified and measured parasites. KV performed all laboratory procedures of molecular biology. TSC conducted DAPC and Mantel tests. KV and RVV conducted all the genetic analyses. KV drafted the manuscript, and AMJ, RG, and RVV revised it. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Roberto do Val Vilela.

Ethics declarations

Ethics approval and consent to participate

Euthanasia of marsupials followed the guidelines of the Animal Care and Use Committee of the American Society of Mammalogists [83] and the Brazilian Guide to Good Practices for Euthanasia in Animals by the Ethics, Bioethics, and Animal Welfare Committee of Brazil’s Federal Council of Veterinary Medicine [84]. Permits for marsupial capture and handling were issued by the Chico Mendes Institute for Biodiversity Conservation (ICMBio/SISBIO, licence numbers 13373-1, 45839-1, and 63023-4); the Environment Foundation of Santa Catarina (FATMA, Environmental Authorization AuA 043/2014/GERUC/DPEC); the State Environmental Institute of Rio de Janeiro (INEA, licence number 020/2011); and the State Institute for the Environment and Water Resources of Espirito Santo (IEMA, Research Authorization GRN n° 026A-2017). Procedures were approved by the ethics committees on animal use of the Oswaldo Cruz Foundation and the Oswaldo Cruz Institute (CEUA-FIOCRUZ, licence number LW-39/14, and CEUA-IOC, licence number L-036/2018).

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.

Municipality and state of origin, hosts, and GenBank accession number of MT-CO1 gene sequences of 175 A. raillieti specimens from this study.

Additional file 2: Figure S1.

Bayesian phylogenetic trees of partial MT-CO1 sequences of A. raillieti from this study and of aspidoderid species from GenBank. The sequence of Lauroia trinidadensis was added as an out-group. Clade I of A. raillieti is not collapsed.

Additional file 3: Figure S2.

Bayesian phylogenetic trees of partial MT-CO1 sequences of A. raillieti from this study and of aspidoderid species from GenBank. The sequence of Lauroia trinidadensis was added as an out-group. Clade II of A. raillieti is not collapsed.

Additional file 4: Figure S3.

Bayesian phylogenetic trees of partial MT-CO1 sequences of A. raillieti from this study and of aspidoderid species from GenBank. The sequence of Lauroia trinidadensis was added as an out-group. Clade III of A. raillieti is not collapsed.

Additional file 5: Figure S4.

Bayesian phylogenetic trees of partial MT-CO1 sequences of A. raillieti from this study and of aspidoderid species from GenBank. The sequence of Lauroia trinidadensis was added as an out-group. Clade IV of A. raillieti is not collapsed.

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

Verify currency and authenticity via CrossMark

Cite this article

Varella, K., Vilela, R.d., Gentile, R. et al. Population genetic structure and phenotypic diversity of Aspidodera raillieti (Nematoda: Heterakoidea), a parasite of Didelphini marsupials in Brazil’s South and Southeast Atlantic Forest. Parasites Vectors 15, 203 (2022).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI:


  • Aspidoderidae
  • Didelphis
  • Philander
  • Fixation index
  • Mantel test
  • Morphology
  • MT-CO1
  • Phylogeny
  • Phylogeography