- Research
- Open Access
- Published:
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 volume 15, Article number: 203 (2022)
Abstract
Background
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.
Methods
We employed morphometric analyses and partial mitochondrial cytochrome c oxidase I gene sequences (MT-CO1) to characterize populations via phylogenetic and phylogeographic analyses.
Results
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.
Conclusion
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

Background
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.
Methods
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].
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).
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 (https://plataformas.fiocruz.br/). 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.
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 3.5.2.2 [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.
Results
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).
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).
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).
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).
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).
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).
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).
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).
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).
Discussion
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].
Conclusion
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.
Abbreviations
- m.a.s.l.:
-
Meters above sea level
- CAR-ES:
-
Cariacica, Espírito Santo
- PET-RJ:
-
Petrópolis, Rio de Janeiro
- RIO-RJ:
-
Rio de Janeiro, Rio de Janeiro
- PTY-RJ:
-
Paraty, Rio de Janeiro
- SGS-MG:
-
São Gonçalo do Sapucaí, Minas Gerais
- CUR-PR:
-
Curitiba, Paraná
- SAI-SC:
-
Santo Amaro da Imperatriz, Santa Catarina
- POA-RS:
-
Porto Alegre, Rio Grande do Sul
- NaCl 0.85%:
-
Physiological saline solution
- DAPC:
-
Discriminant analysis of principal components
- PCs:
-
Principal components
- MT-CO1:
-
Mitochondrially encoded cytochrome c oxidase I
- PCR:
-
Polymerase chain reaction
- PhyPA:
-
Phylogenetic method
- ML:
-
Maximum likelihood
- IB:
-
Bayesian inference
- SMS:
-
Smart Model Selection
- AIC:
-
Akaike information criterion
- aLRT:
-
Approximate likelihood-ratio test
- BP:
-
Bootstrap percentage
- GTR:
-
General time-reversible nucleotide substitution model
- I:
-
proportion of invariable sites
- G:
-
gamma-distributed heterogeneity of substitution rates among sites
- MCMC:
-
Markov chain Monte Carlo
- BPP:
-
Bayesian posterior probabilities
- H:
-
Numbers of haplotypes
- S:
-
Polymorphic sites
- Hd:
-
Haplotype diversity
- π:
-
Nucleotide diversity
- AMOVA:
-
Analysis of molecular variance
- F st :
-
Fixation index
- IBD:
-
Isolation by distance
- AICc:
-
Corrected Akaike information criterion
- ESS:
-
Effective Sample Size
References
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.
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.
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.
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.
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.
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.
Travassos L. Sobre as espécies brasileiras da subfamília Heterakinae Railliet & Henry, 1912. Mem Inst Oswaldo Cruz. 1913;5:271–318.
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.
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.
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.
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.
Gardner AL, editor. Mammals of South America, volume 1: marsupials, xenarthrans, shrews, and bats, vol. 2. Chicago: University of Chicago Press; 2008.
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.
Marques MC, Grelle CE, editors. The Atlantic Forest: History, biodiversity, threats and opportunities of the mega-diverse forest. Berlin: Springer Nature; 2021.
Mittermeier RA, Grill PR, Hoffman M, Pilgrim J. Hotspots revisited: Earth’s biologically richest and most endangered terrestrial ecoregions Cemex. Mexico City. 2004.
Myers N, Mittermeier RA, Mittermeier CG, Da Fonseca GA, Kent J. Biodiversity hotspots for conservation priorities. Nature. 2000;403:853–8.
Blouin MS, Yowell CA, Courtney CH, Dame JB. Host movement and the genetic structure of populations of parasitic nematodes. Genetics. 1995;141:1007–14.
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.
Huyse T, Poulin R, Theron A. Speciation in parasites: a population genetics approach. Trends Parasitol. 2005;21:469–75.
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.
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.
Criscione CD, Poulin R, Blouin MS. Molecular ecology of parasites: elucidating ecological and microevolutionary processes. Mol Ecol. 2005;14:2247–57.
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.
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.
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.
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.
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.
Wilson DE. Opossum predation: Didelphis on Philander. J Mammal. 1970;51:386–7.
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.
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.
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.
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.
Jombart T, Collins C. A tutorial for discriminant analysis of principal components (DAPC) using adegenet, Rvignette. 2015;1–37.
Jombart T. adegenet: a R package for the multivariate analysis of genetic markers. Bioinform. 2008;24:1403–5.
R Core Team (2020). R: A language and environment for statistical computing. R Foundation Statistical Computing. Vienna, Austria. 2020.
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.
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.
Abascal F, Zardoya R, Telford MJ. TranslatorX: multiple alignment of nucleotide sequences guided by amino acid translations. Nucleic Acids Res. 2010;38:W7–13.
Edgar RC. MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 2004;32:1792–7.
Maddison WP, Maddison DR. Mesquite: a modular system for evolutionary analysis, 2021.
Xia X. DAMBE6: new tools for microbial genomics, phylogenetics, and molecular evolution. J Hered. 2017;108:431–7.
Xia X, Lemey P. Assessing substitution saturation with DAMBE. Phylogenet Handb. 2009;2:615–30.
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.
Xia X, Xie Z. DAMBE: software package for data analysis in molecular biology and evolution. J Hered. 2001;92:371–3.
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.
Lefort V, Longueville JE, Gascuel O. SMS: smart model selection in PhyML. Mol Biol Evol. 2017;34:2422–4.
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.
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.
Rambaut A, Suchard MA, Xie D, Drummond AJ. Tracer v1. 6. Computer program and documentation distributed by the author. 2014.
Leigh JW, Bryant D. PopART: Full-feature software for haplotype network construction. Methods Ecol Evol. 2015;6:1110–6.
Bandelt H, Forster P, Röhl A. Median-joining networks for inferring intraspecific phylogenies. Mol Biol Evol. 1999;16:37–48.
Librado P, Rozas J. DnaSP v5: a software for comprehensive analysis of DNA polymorphism data. Bioinform. 2009;25:1451–2.
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.
Wright S. Evolution and the genetics of populations, vol. 4. Variability within and among natural populations. 1978.
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.
Tajima F. Statistical method for testing the neutral mutation hypothesis by DNA polymorphism. Genetics. 1989;123:585–95.
Fu Y. Statistical tests of neutrality of mutations against population growth, hitchhiking and background selection. Genetics. 1997;147:915–25.
Mantel N. The detection of disease clustering and a generalized regression approach. Cancer Res. 1967;27:209–20.
Paradis E, Claude J, Strimmer K. APE: analyses of phylogenetics and evolution in R language. Bioinform. 2004;20:289–90.
Swofford DL. PAUP*. Phylogenetic analysis using parsimony (*and other methods). Version 4. Sunderland: Sinauer Associates; 2003.
Nychka D, Furrer R, Sain S. (2017). Fields: Tools for Spatial Data (R Package).
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.
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.
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.
Freeland JR, Kirk H, Petersen S. Molecular ecology. 2nd ed. Hoboken: Wiley-Blackwell; 2011.
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.
Hartl DL, Clark AG. Princípios de Genética de Populações. 4th ed. Porto Alegre: Artmed; 2010.
Rull V. Speciation timing and neotropical biodiversity: the Tertiary-Quaternary debate in the light of molecular phylogenetic evidence. Mol Ecol. 2008;17:2722–9.
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.
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.
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.
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.
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.
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.
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).
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.
Cerqueira R. The distribution of Didelphis in South America (Polyprotodontia, Didelphidae). J Biogeogr. 1985;135–145.
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.
Hewitt G. The genetic legacy of the Quaternary ice ages. Nature. 2000;405:907–13.
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.
Avise JC. Phylogeography: the history and formation of species. Harvard: Harvard University Press; 2000.
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.
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.
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.
Acknowledgements
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.
Funding
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) #16.03.37.062, 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
Contributions
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
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 http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
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). https://doi.org/10.1186/s13071-022-05288-6
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s13071-022-05288-6
Keywords
- Aspidoderidae
- Didelphis
- Philander
- AMOVA
- Fixation index
- Mantel test
- Morphology
- MT-CO1
- Phylogeny
- Phylogeography