Introduction
The family Subuluridae Travassos, Reference Travassos1914 currently comprises six subfamilies: Allodapinae Inglis, Reference Inglis1958; Labiobulurinae Quentin, Reference Quentin1969; Parasubulurinae Berghe & Vuylsteke, Reference Berghe and Vuylsteke1938; Leipoanematinae Chabaud, Reference Chabaud1957; Subulurinae Travassos, Reference Travassos1914; and Echidnonematinae Smales et al. 2020. This family is characterized by the structure of the buccal end in which the lobes of the oesophagus extend anteriorly to form a pharyngeal portion located between the true oral capsule and the oesophagus (Smales et al., Reference Smales, Elliot and Chisholm2021). Recent work has adopted scanning electron microscopy (SEM) as a complementary tool for the identification and redescription of species belonging to the superfamily Subuluroidea, contributing to a refined description of the species (Smales, Reference Smales2009; Baruš et al., Reference Baruš, Mašová, Koubková and Sitko2013; Du et al., Reference Du, Xu, Li and Li2014; Guo et al., Reference Guo, Zhang, Li and Li2019; Smales et al., Reference Smales, Elliot and Chisholm2021).
Species of the genus Subulura Molin, 1860 (Ascaridida: Subuluroidea) occur in birds, lizards and mammals, with more than 60 valid Subulura species reported worldwide (Vicente et al., Reference Vicente, Sluys, Fontes and Kiefer2000; Baruš et al., Reference Baruš, Mašová, Koubková and Sitko2013). To date, four Subulura species have been reported parasitizing marsupials from South America: Subulura interrogans Lent & Freitas, Reference Lent and Freitas1935, S. Subulura Foster, Reference Foster1939, Subulura trinitatis Wolfgang, Reference Wolfgang1951 and S. amazonica Pereira & Machado Filho, Reference Pereira and Machado Filho1968, with only one Australian species, Subulura peramelis Baylis, Reference Baylis1930, reported.
We recovered parasites of the white-bellied woolly mouse opossum, Marmosa (Micoreus) constantiae Thomas, Reference Thomas1904 in the Brazilian Amazon rainforest, municipality of Sinop, state of Mato Grosso (MT), Brazil. Examination of these parasites revealed a new species of subulurid. We describe the morphology of the species using light microscopy and SEM. In addition, we sequenced part of the mitochondrially encoded cytochrome c oxidase subunit I (MT-CO1) gene of this new species and performed a phylogenetic analysis.
Materials and methods
Marsupial collection
Fifty-three white-bellied woolly mouse opossums were captured in a fragmented landscape of the Amazon rainforest, as part of a larger study described elsewhere (De Mendonça et al., Reference De Mendonça, Colle, Freitas, Martins, Horta, Oliveira and Pacheco2020).
The phylogenetic study also included parasites of the long-furred woolly mouse opossum, Marmosa (Micoreus) demerarae (Thomas, Reference Thomas1905) from the Amazonia biome, municipality of Porto Acre, state of Acre (AC), Brazil and specimens of Primasubulura jacchi recovered from the white-tufted marmoset, Callithrix jacchus (Linnaeus, Reference Linnaeus1758), from the Atlantic Forest biome, municipality of Rio de Janeiro (RJ), state of Rio de Janeiro. All specimens were donated by Laboratório de Biologia e Parasitologia de Mamíferos Silvestres Reservatórios.
Recovery of helminths and morphological identification
All helminths were recovered and processed according to Hoffman, Reference Hoffman1987, with 66 specimens in total (33 males and 33 females). In the laboratory, the helminth specimens were morphologically characterized. First, the specimens were clarified in phenol 90%. Drawings were then produced with the aid of a camera lucida attached to a Zeiss Scope Z1 light microscope (Zeiss, Göttingen, Germany). The structures observed were measured from digital images captured by a Zeiss Axio Cam HRC (Zeiss, Germany), using Carl Zeiss AxioVision Rel. 4.7 accessory software. All measurements were in millimetres (table 1). Identification followed the nematode keys of Vicente et al. (Reference Vicente, Rodrigues, Gomes and Pinto1997) and Anderson et al. (Reference Anderson, Chabaud and Willmott2009). The prevalence, abundance and mean intensity of parasitism for each helminth species recovered were calculated according to Bush et al. (Reference Bush, Lafferty, Lotz and Shostak1997). The type specimens of S. interrogans (CHIOC: 30,334 and 31,232 a–f) and S. amazonica (CHIOC: 30,337 and 30,336) used for comparisons were borrowed from the Helminthological Collection of the Oswaldo Cruz Institute (CHIOC) of Rio de Janeiro, Brazil. Specimens of P. jacchi (Diesing, Reference Diesing1861) and Primasubuura distans (Rudolphi 1809) registered as. Subulura jacchi (CHIOC; 1108) and Subulura distans (CHIOC: 324), respectively, were analysed for comparison.
Table 1. Morphometric data on species of the genus Subulura parasites of marsupials in the Americas.
a this host may otherwise be Caluromys derbianus, based this species distribution.
For SEM, six specimens (four males and two females) were processed according to a protocol modified by Souza et al. (Reference Souza, Lopes-Torres, Garcia, Gomes, Rodrigues-Silva and Maldonado2017). The helminths were dehydrated in a 70–100% ethanol gradient. First, the samples were dehydrated in 70% ethanol for 48 h and then 80%, 90% and absolute ethanol for 20 min, at each step. Finally, the samples were critical point dried in carbon dioxide, mounted on metal stubs and coated with gold (20 nm). The specimens were then examined using a microscope model JEOL JSM-6390 (JEOL, Tokyo, Japan) at the Rudolf Barth Electron Microscopy Platform Oswaldo Cruz Institute, Fiocruz PDTIS/FIOCRUZ.
Molecular and phylogenetic analyses
Genomic DNA of one specimen was extracted using a Qiagen QIAamp DNA Extraction Kit (QIAGEN, Venlo, The Netherlands) following the manufacturer's instructions.
Partial MT-CO1 gene fragments were amplified by polymerase chain reaction (PCR) using primer cocktails described by Prosser et al. (Reference Prosser, Velarde-Aguilar, León-Règagnon and Hebert2013). PCR reactions were performed at a volume of 25 μl for each sample, with reagents in the following volumes: 12.5 μl of PCR Master Mix (Promega Corporation, Madison, WI, USA), 8.5 μl of ultrapure water, 0.5 μl of each forward and reverse primer cocktails (10 μM) and 3 μl of genomic DNA from the sample. Amplifications were performed on a Veriti Thermal Cycler (Life Technologies, CA, USA), following the parameters described by Prosser et al. (Reference Prosser, Velarde-Aguilar, León-Règagnon and Hebert2013). Amplicons (PCR products) were visualized using GelRed (Biotium), on an ultraviolet transilluminator after 1.5% agarose gel electrophoresis.
The amplified products were purified using an Illustra GFX PCR DNA and Gel Band Purification Kit (GE Healthcare, Little Chalfont, Buckinghamshire, UK) following the manufacturer's protocol. Individual primer sequences of the purified samples were obtained using a BigDye Terminator v3.1 Sequencing Kit (Applied Biosystems, USA). This resulted in six bidirectional sequencing reads for improved accuracy. Sequencing was performed using a 96-capillary 3730xl DNA Analyzer (Applied Biosystems) at the Genomic Platform, Fiocruz Technological Platforms Network Oswaldo Cruz Institute (RPT-Fiocruz) of Rio de Janeiro, Brazil.
Sequencing chromatograms were assembled into contigs and edited using the software package Geneious 9.1.8 (Kearse et al., Reference Kearse, Moir and Wilson2012), resulting in a consensus sequence. In addition to the consensus sequence of Subulura eliseae sp. n., we obtained MT-CO1 gene sequences from another specimen of S. eliseae sp. n. parasitizing M. demerarae and from a specimen of P. jacchi, parasite of C. jacchus. The DNA sequences obtained were compared using the BLASTn algorithm with sequences available in the National Center for Biotechnology Information database (http://www.ncbi.nlm.nih.gov).
Our partial MT-CO1 gene dataset included the sequence of a representative of the superfamily Subuluroidea, Aulonocephalus pennula Chandler, Reference Chandler1935. Sequences belonging to species of the superfamilies Ascaridoidea and Heterakoidea were included as outgroups (table 2).
Table 2. Species of Nematoda used in the present molecular analyses.
a State of Mato Grosso, Brazil.
b State of Acre, Brazil.
Sequences of each dataset were aligned using the ClustalW multiple sequence alignment program (Thompson et al., Reference Thompson, Higgins and Gibson1994). The resulting alignments of poorly aligned regions were manually trimmed using Mesquite, version 3.61 (Maddison & Maddison, Reference Maddison and Maddison2018). The substitution saturation level of the matrix was evaluated according to the methods of Xia et al. (Reference Xia, Xie, Salemi, Chen and wang2003) and Xia & Lemey (Reference Xia, Lemey, Lemey, Salemi and Vandamme2009) using DAMBE, version 6.4.79 (Xia & Xie, Reference Xia and Xie2001). Maximum likelihood (ML) inference was performed using the online PhyML 3.0 (Guindon et al., Reference Guindon, Dufayard, Lefort, Anisimova, Hordijk and Gascuel2010). Nucleotide substitution model selection was executed using Smart Model Selection (SMS) (Lefort et al., Reference Lefort, Longueville and Gascuel2017) in PhyML under the Akaike information criterion. Node support in ML trees was determined using the approximate branch likelihood ratio test (aLRT) (Anisimova & Gascuel, Reference Anisimova and Gascuel2006) and by nonparametric bootstrap percentages (ML-BP), with 1000 pseudoreplications, both implemented in PhyML 3.0.
For Bayesian inference (BI), MrBayes, version 3.2.6 (Ronquist et al., Reference Ronquist, Teslenko and Van DerMark2012) on XSEDE was used, using the CIPRES Science Gateway platform (Miller et al., Reference Miller, Pfeiffer and Schwartz2010). Taking into account the different evolutionary processes at each codon position of the MT-CO1 gene, the substitution models were calculated separately for each position using the automated model selection in PAUP*, version 4.0a164 (Swofford, Reference Swofford2003), under the Bayesian information criterion, and implemented with unlinking base frequencies, models and parameters. Markov chain Monte Carlo samplings were performed over 10,000,000 generations, with four simultaneous chains in two runs. Node support in the BI trees was given by Bayesian posterior probability (BPP) calculated from trees sampled every 100 generations after the removal of a burn-in fraction of 25%. To assess the adequateness of our sampling, we used Tracer, version 1.7.1 (Rambaut et al., Reference Rambaut, Suchard, Xie and Drummond2014) to calculate the effective sample size of the parameters. Values above 200 effectively independent samples were considered robust.
Results
Zoobank registration: The Life Science identifier (LSID) for S. eliseae sp. n. urn: lsid: zoobank.org:act:726EBE11-2ABA-47F0-B8C2-E262390C747D
S. eliseae sp. n.
Type-host: Marmosa (Micoreus) constantiae Thomas, Reference Thomas1904
Type-locality: Municipality of Sinop, state of Mato Grosso (MT), Brazil.
Site-of-infection: large intestine.
Prevalence: 13.20% (seven of 53 hosts collected)
Mean intensity: 9.42 (66 helminths collected from seven infected hosts)
Mean abundance: 1.24 (66 helminths collected from 53 hosts collected)
Specimen deposit: Helminthological Collection of the Oswaldo Cruz Institute (CHIOC) of Rio de Janeiro. The accession number of the holotype CHIOC number 39317a (male), accession number of the CHIOC number 39317b allotype (female) and accession numbers of paratype CHIOC number 39317c (one male and one female).
Other material studied: Specimens of S. eliseae sp. n. parasite of M. demerarae, from AC were deposited under the CHIOC number 39318 (one male and one female). Specimens of P. jacchi from the state of Rio de Janeiro were deposited under CHIOC number 39316 (one male and one female). All samples were deposited in a liquid medium.
Etymology: The new species is named in honour of Claudia Elise Xavier de Andrade, mother of the first author.
Description
Medium-sized, whitish nematodes. Cuticle finely striated. Maximum width at mid-body region. Cervical alae developed, starting at the base of cephalic plate and ending just after to posterior excretory pore (figs 1a, e and 2a). Deirids not observed. Cephalic end rounded, without developed lips. Oral aperture simple, circular, surrounded by four ovoid, submedial papillae, and two lateral amphids (figs 1b and 3b). Buccal cavity cylindrical, sclerotized walls and three small pharyngeal complex structures at the base (fig. 1c, d). Muscular oesophagus, terminating in conspicuous valved bulb (fig. 1a). Nerve ring at the first third of the body. Excretory pore slightly posterior to nerve ring (fig. 3a). Tail of both sexes conical, with digitiform tip (fig. 1f, g).
Fig. 1. Subulura elisae sp. n. from Marmosa constantiae in Brazil. (a) female anterior part, lateral view; (b, c) cephalic extremity apical view; (d) cephalic extremity lateral view; (e) male anterior part, ventral view; (f) female posterior end; (g) male posterior end, lateral view; (h) male posterior end, ventral view; (i) tip of spicule, lateral view; (j) gubernaculum, ventral view; (k) region of vulva, lateral view. Scale bars: (a, f, g, h) = 50 μm; (b, c, d, i, j, k) = 10 μm; (e) = 100 μm.
Fig. 2. Photographs of Subulura elisae sp. n. (a) vulva, lateral view; (b) gubernaculum, lateral view; (c) egg; (d) egg located in the ovijector. Scale bars: (a, b, d) = 50 μm; (c) = 10 μm.
Fig. 3. Scanning electron photomicrographs of a specimen of Subulura elisae sp. n. (a) female anterior part, excretory pore (arrowed); (b) apical view, papillae (p) and amphids (am); (c) female posterior part, anus; (d) region of vulva, lateral view.
Male (n = 4): Posterior body ventrally curved, 8.56–9.39 (9.02) long, 354–398 (377) wide. Cervical alae 990–1.046 (1.018) mm long starting at the base of cephalic plate and ending just after to posterior excretory pore. Buccal capsule (fig. 1d), 42–47 (44.2) long, 34–39 (36.5) wide. Oesophagus 1.032–1.117 (1.055) mm long. Bulb rounded (fig. 1a), 193–219 (213.5) long, 182–205 (194) wide. Nerve ring and excretory pore 276–333 (308.5), and 399–510 (455.5), from anterior extremity, respectively. Curved caudal end (fig. 1h). Distance cloaca 133–241 (200.5) to the tail end. Precloacal sucker (figs. 1h and 4a, c), elliptical, 141–179 (159.7) long, 571–678 (622) from posterior extremity. Caudal alae absent. Sessile caudal papillae, ten pairs; three pairs pre-cloacal (fig. 4c, d), two pairs ad cloacal (fig. 4d), five pairs post-cloacal (fig. 4e, online Supplementary fig. S1). Last postcloacal papillae from the tail end 65–77 (72). Spicules similar, slender, slightly curvilinear, alate, distal tips pointed, 1.376–1.475 (1.431) mm long (figs. 1i and 4b). Gubernaculum triangular (figs. 1j and 2b), 190–231 (215.5).
Fig. 4. Scanning electron photomicrographs of a specimen of Subulura elisae sp. n. (a) male posterior part, papillae (ten pairs); (b) spicule (s); (c) precloacal pseudo-sucker; ventral view, papillae (1st pair); (d) precloacal papillae (2nd and 3rd pairs); adcloacal papillae (4th and 5th pairs), ventral view; (e) postcloacal papillae (6th, 7th, 8th, 9th and 10th pairs) lateral view.
Females (n = 5): Body curved ventrally, 12.345–13.942 (13.195) mm long, 435–548 (491.4) wide. Cervical alae 1.141–1.413 (1.250) mm long. Buccal capsule 43–49 (45.2) long, 40–45 (43) wide. Oesophagus, 1.007–1.245 (1.160) mm long. Bulb, 201–231 (218.8) long, 230–261 (240.8) wide. Nerve ring and excretory pore 275–365 (330.8), and 509–569 (542.7), from anterior extremity, respectively. Uterus filled eggs, with circumvolutions extended anteriorly from the vulva to the bulb, 2.323–3.973 (2.910) mm and posteriorly reaching the vulva and/or close to the anus, 5.652–6.063 (5.872) mm. Vulva small 4.331–5.496 (5.203) mm from the anterior extremity (figs. 1k, 2a and 3d). Tail and tip tail long with 910–1.141 (1.031) mm and 70–135 (110.8) (fig. 3c), respectively. Eggs, 71–77 (74) long, 41–50 (45) wide (fig. 2c, d).
Molecular analyses
The MT-CO1 gene sequence of S. eliseae sp. n. from Sinop, MT, Brazil is 601 base pairs (bp) in length. The sequence of the S. eliseae sp. n. from Porto Acre, AC, Brazil is 620 bp in length. To increase the representativeness of subuluroid taxa in our dataset, we amplified and sequenced the partial MT-CO1 gene from the subulurid P. jacchi, obtaining a sequence 612 bp in length. The MT-CO1 sequences of S. eliseae sp. n. (MT and AC) and P. jacchi were deposited in the GenBank database (http://www.ncbi.nlm.nih.gov) (accession numbers: OM432014, OM432016 and OM432015).
The MT-CO1 gene sequences from the present study, aligned with those retrieved from GenBank, resulted in a matrix of eight taxa and 579 characters. From these, 408 characters were constant, 76 were variable, and 95 were parsimony informative. In the ML analyses, the PhyML-SMS selected GTR + G as the best-fit nucleotide substitution model, with optimized ML frequencies, four rate categories and an estimated gamma-shape parameter of α=0.216. The best log-likelihood ML-tree score was -1893.108849. In the BI analyses, PAUP* selected the substitution models TrN + I for the first position, F81 for the second position, and HKY + G for the third position. The BI mean estimated marginal likelihood was -1704.5373 and the median was -1704.272. The effective samples sizes were well above 200 for all parameters.
The ML and BI phylogenies had similar topologies, with little variation in the nodes or support values in the matrix. In all the topologies, the sequences of Subuluroidea species formed a monophyletic group (aLRT= 1.00, ML-BP = 1.00, BPP = 1.00). All analyses confirmed that S. eliseae sp. n. from Acre state and S. eliseae sp. n. from Mato Grosso state formed a monophyletic group (aLRT = 1.00, ML-BP = 0.99, BPP = 1.00), with P. jacchi a sister to the monophyletic S. eliseae sp. n., with moderate to high support (aLRT = 1.00, ML-BP = 0.60, BPP = 0.82). Aulonocephalus pennula grouped with the clade formed by P. jacchi and S. eliseae sp. n., with moderate support (aLRT = -, ML-BP = 0.97, BPP = 1.00) indicating monophyly of the Subuluridae family (fig. 5).
Fig. 5. Maximum likelihood tree (ML) based on an analysis of mitochondrial cytochrome c oxidase I sequences inferring a phylogenetic relationship between Subulura elisae sp. n. and other sequences from GenBank. Node values are aLRT, ML-BP, and BPP supports.
Discussion
Currently, the identification of Subulura species is based mainly on morphological and morphometric characters, that is, body size, lengths of the oesophagus and spicules, number and arrangement of caudal papillae, morphology of the gubernaculum and position of the vulva (Guo et al., Reference Guo, Zhang, Li and Li2019). We allocated the new species to the subfamily Subulurinae due to its complex pharyngeal structure with three pharyngeal portions, mouth opening simple and/or hexagonal and males with an elongated precloacal pseudo-sucker without a defined border.
The new species can be distinguished morphologically from the other four Subulura parasite species of marsupials by the number of caudal papillae, structure and dimensions (morphometric) and size of the spicule. Subulura eliseae sp. n. is smaller than S. lanigeri, but both species are longer than S. interrogans and S. trinitatis. It does not have caudal alae, a structure described in S. interrogans and S. trinitatis. The number of caudal papillae (ten pairs) found in S. eliseae sp. n. also differs from that found in other Subulura spp. described in South American marsupials.
In their description of S. amazonica and characterization the S. interrogans female allotype Pereira & Machado Filho (Reference Pereira and Machado Filho1968) mentioned three pairs of tooth-like pharyngeal projections. However, they seem not to have accurately documented an important Subulura characteristic. Quentin (Reference Quentin1969) studied ontogenesis of the cephalic structure of the family Subuluroidea, describing three small pharyngeal portions of a complex structure at the base of the buccal capsule of the family Subuluridae species. Presuming that the pharyngeal projections described by Pereira & Machado Filho (Reference Pereira and Machado Filho1968) are the pharyngeal portions of the complex structure described by Quentin (Reference Quentin1969), we consider that there are three pharyngeal formations in Subulura species, as suggested by Quentin (Reference Quentin1969), instead of three pairs, as suggested by Pereira & Machado Filho (Reference Pereira and Machado Filho1968) (table 2).
In conclusion we describe a new Subulura species, S. eliseae sp. n., parasitizing mouse opossums Marmosa spp. As molecular data for the family Subuluridae species are insufficient, few inferences could be drawn from our phylogenies. As expected, our two sequences of S. eliseae sp. n. from Acre and Mato Grosso states formed a strongly supported clade, sister to the sequence of P. jacchi, thus supporting the subfamily Subulurinae, sister to A. pennula, which belongs to the subfamily Allodapinae, finally supporting the family Subuluridae.
Subulura eliseae sp. n. is the fifth species of the genus to be described in marsupials from the Americas, a genus that now comprises 61 species. This is only the second Subulura and the third Subuluridae species to be molecularly characterized. Adding to this, our sequence of S. jacchi is the fourth Subuluridae species thus far to be molecularly characterized. At present, morphological knowledge about many Subulura species is scattered or poorly reported. Thus, an integrative taxonomic approach comprising more species of the genus Subulura and family Subuluridae, including morphological, ecological and molecular characters, is essential to improve our understanding of their evolutionary and ecological associations with their hosts.
Supplementary material
To view supplementary material for this article, please visit https://doi.org/10.1017/S0022149X22000244
Acknowledgements
We thank the Universidade Federal Mato Grosso for the donation of the helminths used in this study. We are grateful to Dr Marcelo Knoff for providing literature. Thanks also to Dr Sócrates Fraga da Costa Neto for collecting the specimen of Primasubulura.
Financial support
Funding for host collection was provided by the National Council for Scientific and Technological Development – CNPq (#447557/2014-9; #310352/2016) and the Foundation for Research Support of Mato Grosso State – FAPEMAT (#477017/2011). Authors BEAS, LCF, and RFBM received the funding of a scholarship from the Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES) – Brazil – Finance code 001.
Conflict of interest
None.
Ethical approval
The authors followed all applicable institutional and national laws and guidelines during this research.
Open access
