- Research
- Open access
- Published:
A cyst-forming coccidian with large geographical range infecting forest and commensal rodents: Sarcocystis muricoelognathis sp. nov.
Parasites & Vectors volume 17, Article number: 135 (2024)
Abstract
Background
The geographic distribution and host-parasite interaction networks of Sarcocystis spp. in small mammals in eastern Asia remain incompletely known.
Methods
Experimental infections, morphological and molecular characterizations were used for discrimination of a new Sarcocystis species isolated from colubrid snakes and small mammals collected in Thailand, Borneo and China.
Results
We identified a new species, Sarcocystis muricoelognathis sp. nov., that features a relatively wide geographic distribution and infects both commensal and forest-inhabiting intermediate hosts. Sarcocystis sporocysts collected from rat snakes (Coelognathus radiatus, C. flavolineatus) in Thailand induced development of sarcocysts in experimental SD rats showing a type 10a cyst wall ultrastructure that was identical with those found in Rattus norvegicus from China and the forest rat Maxomys whiteheadi in Borneo. Its cystozoites had equal sizes in all intermediate hosts and locations, while sporocysts and cystozoites were distinct from other Sarcocystis species. Partial 28S rRNA sequences of S. muricoelognathis from M. whiteheadi were largely identical to those from R. norvegicus in China but distinct from newly sequenced Sarcocystis zuoi. The phylogeny of the nuclear 18S rRNA gene placed S. muricoelognathis within the so-called S. zuoi complex, including Sarcocystis attenuati, S. kani, S. scandentiborneensis and S. zuoi, while the latter clustered with the new species. However, the phylogeny of the ITS1-region confirmed the distinction between S. muricoelognathis and S. zuoi. Moreover, all three gene trees suggested that an isolate previously addressed as S. zuoi from Thailand (KU341120) is conspecific with S. muricoelognathis. Partial mitochondrial cox1 sequences of S. muricoelognathis were almost identical with those from other members of the group suggesting a shared, recent ancestry. Additionally, we isolated two partial 28S rRNA Sarcocystis sequences from Low’s squirrel Sundasciurus lowii that clustered with those of S. scandentiborneensis from treeshews.
Conclusions
Our results provide strong evidence of broad geographic distributions of rodent-associated Sarcocystis and host shifts between commensal and forest small mammal species, even if the known host associations remain likely only snapshots of the true associations.
Graphical Abstract
Background
Sarcocystis spp. of rodents represent a highly diversified group of tissue cyst-forming coccidia, for which the patterns of species diversification, geographical distribution and host associations remain poorly known. Over the last 5 decades, various species of Sarcocystis have been observed in the striated musculature of commensal Rattus spp. and related rodent species, including the widespread Sarcocystis cymruensis [1,2,3,4,5], S. kani [6], S. murinotechis [7], S. ratti [5, 8], S. singaporensis, S. villivillosi, S. zamani [2, 9], S. sulawesiensis [2, 10], S. zuoi [11] and a Sarcocystis sp. from masked owls [12]. Sarcocystis rodentifelis [13, 14] is considered synonymous with S. cymruensis [3]. In this group, Sarcocystis spp. with snake-rodent life cycle are among the most prevalent Apicomplexan parasites found in rodents in the Indo-Australian Archipelago [2, 10, 15], also showing prominent presence in environmental samples [16].
Sarcocystis murinotechis has been a somewhat enigmatic parasite, because since its early description from infected elapid snakes (Notechis ater; now, Notechis scutatus [Peters, 1861], see [17]) and rodents in Australia, i.e. Rattus rattus, R. norvegicus, R. lutreolus, Pseudomys higginsi and Mastacomys fuscus [7, 18], no molecular information has become available. Furthermore, the ultrastructure of its sarcocyst shares some similarities with Sarcocystis spp. from Eastern and Southeast Asia, such as S. zuoi from China, which develops in colubrid snakes as definitive host [11], and a species described as S. murinotechis-like that was observed in rats in Thailand and neighboring countries [2, 15]. Thus, more morphological as well as molecular data were needed to discriminate between these species. Watthanakaiwan et al. [19] amplified nuclear rRNA gene sequences from a Sarcocystis species observed in Rattus tiomanicus in Thailand, which appeared to be molecularly close to S. zuoi and was subsequently cited as this species although the authors did not provide morphological evidence for its cyst wall structure. In a recent study, we found that one of the 18S rRNA gene sequences (KU341120) published by Watthanakaiwan et al. [19] was genetically distinct from S. zuoi and hypothesized that it could belong to another species [20].
Here, we synthesized comparative morphological and molecular evidence from new and previously collected specimens from Thailand, Borneo and China and insights from transmission experiments using sporocysts from naturally infected colubrid (rat) snakes in Thailand to clarify the status of the S. murinotechis-like species in Thailand and its potential distribution in other countries of the region. We describe a new Sarcocystis species and discuss its phylogenetic position and host relationships within the so-called S. zuoi complex, which is an expanding group of molecularly cryptic Sarcocystis spp. that infects colubrid snakes and small mammals in Asia [6].
Methods
Sample collection and examination
Sarcocystis oocysts/sporocysts were recovered from the intestines of two road-killed, adult colubrid snakes, one copperhead rat snake, Coelognathus radiatus, collected in the Bangkok metropolitan area, and one yellow striped snake, Coelognathus flavolineatus, which had been collected in Chumphon Province in southern Thailand. Intestinal scrapings containing sporocysts were further purified on a discontinuous sucrose density gradient as described previously [6].
The rodents from Borneo examined for this study included two Maxomys whiteheadi, three Niviventer cremoriventer and four Sundasciurus lowii that were live-trapped between 2012 and 2013 using locally made drop-door, wire-mesh traps in different habitat types from mature forest to urban habitats in the Kota Kinabalu and Ranau Districts in Sabah, Malaysia. Species identification was based on phenotypic characteristics according to Aplin et al. [21] and Musser and Carleton [22]. Captured animals were transferred into small cotton bags and killed by cervical dislocation. Samples of striated musculature from diaphragm and latero-ventral abdominal wall (Musculus obliquus externus) or M. quadriceps femoris were fixed and stored in 70% ethanol. Histological preparation of sarcocysts in muscle samples and subsequent processing for ultrastructural examination have been described [20].
In China, a total of 189 Norway rats were captured between 2017 and 2021 by live trapping on farmland in Anning and Nanjian Prefectures, both located in southwestern China. After the animals had been killed with CO2, their corpses were examined in a laboratory for presence of sarcocysts. This included light microscopical inspection of fresh preparations of musculature from tongue, esophagus, diaphragm, various skeletal muscles and the heart. Single sarcocysts were separated from surrounding tissue using dissection needles and subsequently processed for transmission electron microscopy (TEM) as described [23].
Light microscopic size measurements of parasite stages and statistical analysis
The sizes of sporocysts and cystozoites of the new Sarcocystis isolates under study were measured and compared with other species by using one-way ANOVA as previously described [6].
Husbandry and experimental infection of laboratory rats
All experiments involving the transmission of Sarcocystis sporocysts to a rodent host were performed using outbred Sprague-Dawley (SD) rats, which were obtained from the National Laboratory Animal Centre of Mahidol University in Thailand. SD rats were kept under coccidia-free conditions in isolated and climate-controlled rooms of the Agricultural Zoology Research Group at the Department of Agriculture in Bangkok, Thailand. Two individuals were housed together per cage (Makrolon polycarbonate: 60 × 38 × 20 cm) under controlled temperature and humidity conditions, exposed to a 12:12 h light/dark cycle and ad libitum diet and water supply.
Sporocysts isolated from one Coelognathus radiatus and one C. flavolineatus were used to inoculate three male, 3-month-old SD rats each per isolate; additionally, two rats each were inoculated with PBS and kept in the same room serving as negative controls. The animals were inoculated with about 1000 purified sporocysts each using a stomach tube. Afterwards, the animals were monitored daily for potential signs of disease. All six parasite-inoculated rats and the four negative control rats were killed after 5 months with CO2 and necropsied. We examined fresh muscle tissue from esophagus, tongue, diaphragm, muscles of the abdominal wall, femoral muscles and masseter for presence of sarcocysts under a light microscope. The size of freshly released cystozoites was determined after rupturing individual sarcocysts in phosphate-buffered saline at room temperature. Samples of muscle tissue were fixed and processed for electron microscopy as described previously [24].
Gene sequencing and sequence comparisons using BLAST
Only samples from Borneo and China were available for molecular studies, as samples taken during the laboratory infection experiments in Thailand had been destroyed during storage. For samples from Borneo, DNA extraction was performed as described in detail previously [20]; 25 µg samples of striated musculature were used that had been checked for presence of sarcocysts by LM and TEM. Nuclear 28S rRNA gene sequences were amplified following the protocol of Prakas et al. [25] using the primers KL-P1F and KL-P1R (Table 1). Raw sequences were manually checked and refined using MEGA-6 [26]. For samples from China, four sarcocysts derived from different individual Norway rats were used for molecular analysis. DNA was extracted from each sarcocyst using a commercial kit (TIanamp Genomic DNA, Tiangen Biotech®, Beijing, China) according to the manufacturer’s instructions. Four genes, nuclear 18S rRNA, 28S rRNA, ITS1-5.8S-ITS-2 and mitochondrial cox1, were selected for the characterization of this species. The primers used are listed in Table 1. PCR amplifications were performed as previously described [23]. PCR products were gel purified using an E.Z.N.A.® Gel Extraction Kit (Omega Bio-Tek, Inc., Georgia, USA) and ligated to pCE2 TA/Blunt-Zero vector using a 5-min TA/Blunt-Zero Cloning Kit (Vazyme Biotech Co., Ltd., Nanjing, China) according to the product manuals. The ligated vectors were transformed to Trelief® 5α Chemically Competent Cell (Tsingke Biotechnology Co., Ltd., Beijing, China). The positive bacterial colonies were selected for bi-directional sequencing on an ABI PRISM TM 3730 XL DNA Analyzer.
Initial GenBank screening and pairwise similarity comparisons of the amplified sequences were performed using the web-based Basic Local Alignment Search Tool (BLASTn) of the National Center for Biotechnology Information of the National Institutes of Health, USA.
Phylogenetic analysis
Alignment of rRNA sequences and identification of gene boundaries
Before alignment of the 18S and 28S rRNA sequences and the Internal Transcribed Spacer 1 (ITS1), we identified their gene boundaries based on the predicted secondary structure of the rRNA of Toxoplasma gondii as published by Gagnon et al. [27] and truncated the ends if necessary. For the 18S rRNA gene, we additionally employed secondary structure modeling using SSU-ALIGN software (version 0.1.1, February 2016) as described previously [6]. Full-length 28S rRNA sequences were truncated from the 3’-end to a maximum of 2500 nucleotides to accommodate the analysis limits of the alignment software.
The new ITS1 sequences described here were carefully inspected and excised from the amplified ITS1-5.8S-ITS2 gene complex. For sequence comparisons, 25 Apicomplexan ITS1 sequences were selected that were labeled as being complete by GenBank. Especially, we checked whether they started after the 18 s rRNA [28] and were consistent with the start motifs that followed the terminal motif (3ʹ-end) of the nuclear 18S rRNA gene of T. gondii: 5ʹ-TCATTCA-3ʹ. After this motif, ITS1 sequences typically started with 5ʹ-CACA-3ʹ (e.g. Sarcocystis cf. zuoi, KU341120), 5ʹ-AACG-3ʹ (e.g. S. neurona, AF081944), 5ʹ-CACC-3ʹ (e.g. S. hirsuta, KT901231) and 5ʹ-CACG-3ʹ (e.g. T. gondii, AY259045), or similar. To identify the terminal nucleotides of an ITS1 sequence, we searched for the potential presence of the start motif of the 5.8S rRNA of T. gondii (5ʹ-AAATTTT’-3ʹ) according to Gagnon et al. [27], which was slightly modified in Sarcocystis spp. to 5ʹ-ACAATTTT-3ʹ and included in a stretch of about 230 nucleotides (nt) of the highly conserved 5.8S rRNA gene. Accordingly, we truncated all newly sequenced ITS1 sequences before this start motif, which resulted in ITS1 sequences that were 629–675 nt long. We did not alter the ends of the other, publicly available ITS1 sequences of tissue-cyst forming coccidia, which were analyzed as they were, except in the case of the four ITS1 sequences that had been designated as ‘Sarcocystis zuoi’ from Thailand (KU341118/19/20/21) and deposited as gene complex, requiring excision for further alignment.
The 18S, 28S rRNA and ITS1 sequences of the Sarcocystis isolates under investigation were aligned with 38, 25 and 25 published Apicomplexan sequences available at GenBank [29], respectively, using the multiple sequence alignment algorithm of the ‘R-Coffee’ web server, which considers the predicted secondary structure of RNA [30]. The alignment procedure was repeated several times, while we finally selected three independent alignments with the smallest proportion of ambiguous positions for each gene. Conserved areas were inspected visually for potential misalignments.
Phylogenetic reconstruction of gene trees
Reconstruction of phylogenetic trees in the case of the 18S and 28S rRNA genes was performed by the maximum likelihood (ML) and minimum evolution (ME) methods as implemented in Mega-X [31]. The ME criterion was used because it provided a distance-based approach to tree-building in contrast to the character-based ML [32]. As implemented in Mega-X, ME also allowed for distance estimation under heterogeneous substitution patterns among lineages [33], a condition that applies to coccidian parasites [34]. Our initial analyses of the 18S had shown that ML could not resolve (into dichotomies) the clade containing the new Sarcocystis species and S. zuoi. Initial tree(s) for the heuristic search under ML were obtained automatically by applying neighbor-joining and BioNJ algorithms. The ME tree was searched using the close-neighbor-interchange (CNI) algorithm at a search level of 2; the neighbor-joining algorithm was used to generate the initial tree. For analysis of the ITS1 locus, initial tree(s) were obtained by applying neighbor-joining and BioNJ algorithms to a matrix of pairwise distances estimated using the maximum composite likelihood (MCL) approach. For all analyses, 1000 bootstrap replicates were performed.
The appropriate models for nucleotide substitution rates were determined by ML analysis using MEGA-X, selecting for each gene the model with the lowest BIC (Bayesian information criterion) and AICc (Akaike information criterion, corrected) values. For analysis of a longer stretch of the nuclear 28S rRNA gene (comparison of 1300–1400 homologous characters), the Tamura-Nei substitution model was used with a discrete gamma distribution rate variation (four categories; + G, parameter = 0.6874) while the rate variation model allowed for some sites to be invariable (+ I; 32.39% of sites). For shorter 28S rRNA sequences (360–400 characters), the Kimura two-parameter model was employed with a discrete gamma distribution (four categories; + G, parameter = 0.2299). For the nuclear 18S rRNA gene under ML analysis, we selected the Tamura three-parameter substitution model with a discrete gamma distribution (four categories; + G, parameter = 0.1537) and some sites being evolutionarily invariable ([+ I], 38.77% sites). For ME phylogenetic analysis of the 18S rRNA, the rate variation among sites was modeled with a gamma distribution (shape parameter = 0.16); the differences in the composition bias among sequences were considered in evolutionary comparisons. The ITS1 sequences were analyzed using the Hasegawa-Kishino-Yano substitution model under a discrete gamma distribution to model evolutionary rate differences among sites (four categories; + G, parameter = 4.2949).
For comparison of other Sarcocystis spp. with the novel Sarcocystis isolates, we included all sequences in tree building that showed a high degree of sequence similarity during initial GenBank screens with BLASTn. Second, we selected sequences of known Sarcocystis spp. parasitizing snakes as definitive hosts and representative taxa of major lineages of the Sarcocystinae and Toxoplasmatinae as described [6, 35]. Eimeriid coccidia served as outgroup. Accession numbers of all sequences are given in the phylogenetic trees of the respective genes.
Results
Sarcocyst development in experimentally infected rats
The sporocysts of Sarcocystis isolated from the two ratsnake species in Thailand, Coelognathus radiatus and C. flavolineatus, were indistinguishable in size, measuring on average 10.07 (± 0.39; s.d.) × 7.18 (± 0.25) µm (n = 20) and 10.28 (± 0.80) × 7.51 (± 0.73) µm (n = 21), respectively (Fig. 1a). After inoculation of six SD rats with both sporocyst isolates, we observed sarcocysts in the striated musculature of all six rats that were examined 5 months after infection. In fresh tissue preparations, sarcocysts of both isolates were numerous and morphologically identical, measuring 400–700 × 70–90 µm (length × diameter; n = 25 from four rats). We observed sarcocysts in the striated musculature of abdomen, thorax, limbs, diaphragm and esophagus; their form was cigar-shaped, while longer cysts could show slightly pointed ends. Live sarcocysts, which had been freed from host tissue, showed a rim of palisade-like, villar protrusions that were between 8 and 10 µm long and about 2–4 µm wide at their base; however, protrusions at the tip of a cyst could be longer, measuring between 10 and 16 µm with a base that was on average 3.6 µm wide (range: 2.5–4.1 µm, n = 15). The interior of native sarcocysts appeared to be septated (Fig. 1b). Banana-shaped cystozoites, freshly released from ruptured, native sarcocysts, were on average 5.9 (± 0.6; s.d.) µm long and 2.2 (± 0.4) µm in diameter (n = 21; both isolates).
Ultrastructurally, the mean diameter of the sarcocyst wall (= primary cyst wall [PCW] plus ground substance [GS]) was 10.16 (± 2.34) µm (n = 20), while the villar protrusions measured on average 8.65 (± 2.05; n = 20; two cysts, one each from each isolate) µm in length and 2.50 (± 0.61) µm in diameter. The electron-dense, knob-like structures of the PCW were about 150 nm wide at the base of the villar protrusions, whereas their width was about 75 nm in the more apical parts, where knobs could fuse to form a continuous layer of PCW (Fig. 1c). The layer of GS underneath the villar protrusions was thin, on average about 270 nm wide (Fig. 1c, e). Notably, cystozoites of the isolates from Thailand only possessed two rhoptries, which was evident from longitudinal as well as cross sections (Fig. 1d, e). In cross sections through the anterior third of cystozoites, we counted on average 60 micronemes (± 14; n = 6). These micronemes were loosely scattered throughout the anterior half of the cell; they appeared short and stumpy and were about 200 nm long. Dense granules were present but rarely observed. Metrocytes occurred regularly, more frequently towards the ends of a sarcocyst; apparently, they divided by endodyogeny since they contained only two zoites in all observed cases (Fig. 1e). As shown in Fig. 1f, large cyst wall protrusions could attain a length well over 10 µm.
The four negative control rats examined 5 months after infection were free of any Sarcocystis infection. All rats remained healthy throughout the experiment.
Morphology of sarcocysts from wild rodents in Borneo and China
Sarcocyst-infected muscle tissue from Borneo was isolated from two wild-caught Maxomys whiteheadi (Whitehead’s spiny rat), which showed numerous sarcocysts with a striated wall in histological sections of musculature of the abdominal wall. Because cysts were mainly in the cross-sectional plane, we can only provide their mean diameter, which was 85.6 (± 17.1) µm (n = 5). At higher magnification (Fig. 2a), the cyst wall consisted of broad, palisade-like protrusions which were on average 5.1 (± 2.5) µm long (n = 10). At the tip of one cyst, protrusions were longer, measuring 6.8–9.1 µm. However, because most of the protrusions in histological sections appeared to be bent or compressed, their actual size was probably larger, as measurements from live sarcocysts from other sampling locations suggested. The protrusions of the sarcocyst wall were indistinguishable from the samples from Thailand and China in terms of ultrastructure (Fig. 2b). Although preservation of ultrastructural details of the Bornean sample was relatively low, the typical, knobbed and electron-dense lining of the PCW of the protrusions was clearly visible including presence of electron-light granular material filling the core. On average, protrusions were 7005 (± 1235) nm long (n = 6); their diameter in the apical part ranged between 1940 and 3690 nm, at the base between 1625 and 2610 nm (n = 10); thus, some of the protrusions appeared enlarged towards the apical part. Two fully stretched cystozoites measured 4300–4940 × 1060–1259 nm, while all observed cystozoites possessed only two rhoptries like the isolates from Thailand (Fig. 2c). Cystozoites in H&E-stained histological sections measured 4.6 (± 0.4) × 1.2 (± 0.2) µm (n = 10).
Muscle cysts found in 14 out of 189 (7.4%) Norway rats from southwestern China were identical to those observed for the Sarcocystis isolates from Thailand and Borneo (Fig. 2e). There were numerous sarcocysts in the skeletal muscles, but none in the heart. Under the light microscope, live sarcocysts measured 770–1580 × 87–123 μm in size (n = 20). The cysts had palisade-like villar protrusions measuring 8.5–11 μm in length and 1.7–2.4 μm in width (n = 30) (Fig. 2d); they were filled with numerous banana-shape cystozoites measuring 6.5 (± 0.6; s.d.) × 2.1 (± 0.15) μm in size (n = 25). At the ultrastructural level, the primary cyst wall contained numerous vertical, palisade-like protrusions, measuring 6.5–10.8 × 1.6–3.1 μm in size (n = 30), which contained electron-light granular matter in their core (Fig. 2e). The PCW consisted of a knobbed, electron-dense layer and showed minute undulations over its entire surface. A layer of ground substances with 200–300 nm thickness was located immediately underneath the primary sarcocyst wall.
Morphological comparison of the novel Sarcocystis sp. with closely related taxa
The mean (± s.d.) sporocyst length and diameter of both isolates of Sarcocystis sp. from Thailand given as 10.07 (± 0.39) × 7.18 (± 0.25) µm (n = 20) (from Coelognathus radiatus) and 10.28 (± 0.80) × 7.51 (± 0.73) µm (n = 21) (from C. flavolineatus) compared to other studies as follows: Sarcocystis clethrionomyelaphis isolated in Germany (11.5 [± 0.6] × 8.7 [± 0.6] μm, n = 12; [36]), S. clethrionomyelaphis from China (11.0 [± 0.8] × 8.0 [± 0.5] μm, n = 25; [37]), S. zuoi from China (10.8 [± 0.7] × 8.0 [± 0.7] μm, n = 25; [11]), S. attenuati isolated in China (9.9 [± 0.4] × 6.6 [± 0.2] μm, n = 30; [23]), S. kani from Thailand (11.8 [± 0.5] × 7.5 [± 0.2] μm, n = 14; [6]) and S. murinotechis from Australia (11.0 [± 0.5] × 7.25 [± 0.7] μm, n = 10; [7]) (Fig. 3). Sporocyst lengths of Sarcocystis sp. from C. radiatus were significantly shorter than in S. clethrionomyelaphis, S. kani, S. murinotechis and S. zuoi (Holm-Sidak post hoc P < 0.001 following one-way ANOVA test with F(7, 149) = 22.5, P < 0.001) but not in the case of S. attenuati (P = 0.321) and Sarcocystis sp. isolated from C. flavolineatus (P = 0.489). The mean sporocyst diameter of Sarcocystis sp. from C. radiatus was significantly larger than that measured for S. attenuati and smaller than for S. clethrionomyelaphis, S. kani and S. zuoi (post hoc P < 0.001 following ANOVA with F(7, 149) = 30.4, P < 0.001) whereas no apparent difference was found in the measured diameters of S. murinotechis (P = 0.738) and Sarcocystis sp. isolated from C. flavolineatus (P = 0.136); the power of all performed tests with alpha = 0.05 was 1.
We also examined the size of live cystozoites of the new Sarcocystis sp. from Thailand and China (mean length × diameter in µm and s.d.: Thailand, 5.9 [± 0.6] × 2.2 [± 0.4], n = 21; China, 6.5 [± 0.6] × 2.1 [± 0.15], n = 25) in comparison with S. clethrionomyelaphis from China (9.5 [± 1.0] × 2.7 [± 0.1], n = 40; same authors as above), S. kani (7.7 [± 0.4] × 2.5 [± 0.2], n = 7), S. zuoi (8.5 [± 0.8] × 2.8 [± 0.1], n = 35), S. attenuati (9.2 [± 0.7] × 2.5 [± 0.3], n = 40) and S. eothenomysi (6.8 [± 0.4] × 3.1 [± 0.3], n = 23; [38]); the latter species represented lineage S2. In comparison, the isolate of Sarcocystis sp. from Thailand was indistinguishable in size from the isolate of Sarcocystis sp. from China but was significantly different from all other examined species: The cystozoite length of Sarcocystis sp. from Thailand was significantly shorter than those of all studied species (post hoc Holm-Sidak test P < 0.05 with ANOVA F(6, 170) = 80.565, P < 0.001) but indistinguishable in length from those in R. norvegicus in China (P = 0.074). Likewise, cystozoite diameters of Sarcocystis sp. from Thailand were significantly smaller than those from all other species (post hoc P < 0.05 with ANOVA F(6, 170) = 50.6, P ≤ 0.001) but indistinguishable in size from those in R. norvegicus from China (P = 0.187); power of all performed tests = 1 (alpha = 0.05).
Because there were no data on the size of cystozoites of S. murinotechis available, we re-examined the length of its cystozoites taken from an image of a sarcocyst in muscle tissue published by Munday et al. [18] that was associated with S. murinotechis [7]. In this sarcocyst, cystozoites of S. murinotechis were 7.2 (± 0.8) µm long (n = 12). In our own histological sections of sarcocysts of Sarcocystis sp. from rats inoculated with sporocysts from Coelognathus radiatus and C. flavolineatus, cystozoites measured 4.2 (± 0.5) µm (n = 14) and 4.5 (± 0.9) µm (n = 6) in length, respectively. When we compared these samples of Sarcocystis sp. (plus cystozoite length in tissue sections from Borneo, see previous section) with data on S. kani [6], S. scandentiborneensis [20] and the measurements of S. murinotechis, cystozoites of Sarcocystis sp. from Thailand (ex C. radiatus) and Borneo (ex M. whiteheadi) were significantly shorter than S. kani, S. murinotechis and S. scandentiborneensis (Holm-Sidak’s P < 0.001 with ANOVA F(5, 136) = 35.3, P < 0.001) but not different in intraspecific comparison: Sarcocystis sp. ex C. radiatus versus Sarcocystis sp. ex M. whiteheadi from Borneo; P = 0.198; Sarcocystis sp. ex C. radiatus versus Sarcocystis sp. ex C. flavolineatus; P = 0.360; power of all performed tests = 1 (alpha = 0.05).
Taken together, these comparisons provided substantial evidence that the new species showed significant morphological disparity (including its sporocysts) to S. murinotechis, S. zuoi and other Sarcocystis spp., while the cystozoites of Sarcocystis sp. from Thailand, Borneo and China were similar in size.
Molecular evidence for the new Sarcocystis sp. from intermediate hosts in Borneo and China and other novel Sarcocystis sequences from forest rodents
In China, targeted genetic markers were successfully sequenced from four individuals of R. norvegicus, including the nuclear 18S and 28S rRNA genes, the ITS1-5.8S-ITS2 rRNA gene complex and mitochondrial cox1 (Additional file 1: Table S1). From the samples in Borneo, we were able to amplify three partial 28S rRNA sequences from cyst-infected musculature of two specimens of M. whiteheadi, of which the sequence E369.13-28S (GenBank OR979643) stood out because of its reasonable length (610 nt) useful for phylogenetic analysis. The other two sequences from a second specimen of M. whiteheadi were too short for tree building (isolates E389.13A-28S, 142 bp, and E389.13B-28S, 137 bp) but 100% identical with isolate E369.13-28S (Additional file 1: Table S1), potentially indicating that these sequences also belonged to the new Sarcocystis sp. Moreover, sarcocyst morphology (LM and TEM) in the second forest rat was identical to that in the first specimen (Fig. 2b, c). We also amplified a sequence from another forest-dwelling rat, Niviventer cremoriventer (OR979644, isolate E388.13A-28S, 413 nt), which was 99.03% identical to the sequences of Sarcocystis sp. from Borneo and China but showed the highest identity score (99.76%) with sequence E357.13-KLP1F (OR979647; Additional file 1: Table S1), a partial 28S rRNA sequence that was amplified from the same sarcocyst DNA (isolated from a treeshrew) which was deposited as holotype of Sarcocystis scandentiborneensis [20]. Although S. scandentiborneensis was so far known to develop in treeshrews only, we amplified two more putative 28S rRNA gene sequences of this species from muscle samples of rodents: here, two specimens of Low’s squirrel Sundasciurus lowii, while these sequences (E346.13-28S, E367.13-28S) showed the highest identity score (99.82%) with E357.13-KLP1F of S. scandentiborneensis in pairwise comparisons (Additional file 1: Table S1). The light microscopic appearance and ultrastructure of the corresponding sarcocysts from the squirrels also corresponded to features of S. scandentiborneensis, corroborating the link to that species. Presumably intraspecific sequence comparisons at around 600-bp alignment length of the Sarcocystis sp. 28S rRNA E369.13-28S from M. whiteheadi with isolate Rn-cl2.3-28S from Norway rats in China (including the other haplotypes) and KU341120 from R. tiomanicus in Thailand all resulted in 99% identity. Although sequence S. scandentiborneensis E357.13-KLP1F also showed 99% identity with E369.13-28S of the new Sarcocystis sp. from Whitehead’s spiny rat, the alignment length under pairwise BLAST analysis was around 60 bp shorter, apparently too short to include the bp differences that would have separated the two species by the similarity score. This was indicated by the result that another (putative) sequence of S. scandentiborneensis that was similar in length (614 bp, isolate E367-13-28S, OR979645) only showed a similarity score of 98.37% in pairwise alignment with sequence E369.13-28S; both sequences included the expansion segment D2 at the 3’-end. Moreover, character-based phylogenetic reconstruction of the 28S rRNA tree under the ML criterion using the same sequences could distinguish well between the two species (see below). Importantly, BLAST comparisons of Sarcocystis sp. from Borneo and China with a newly sequenced, full-length 28S rRNA of S. zuoi (OR979648, isolate RN-cl2-3-28S) revealed a sequence identity of only 97% indicating that Sarcocystis sp. was different from S. zuoi at the 28S rRNA locus (Additional file 1: Table S1).
Since we could not examine our own DNA samples of Sarcocystis sp. from Thailand, we regarded isolate KU341120 from Thailand as the most likely candidate to be conspecific with the former species. Apart from the high identity scores between the 28S rRNA sequence KU341120 and the partial, homologous sequences of Sarcocystis sp. from Borneo and China (in the latter comparison over a stretch of about 1500 bp), alignment of the 18S rRNA gene sequence KU341120 with the new, full-length 18S rRNA from China (Sarcocystis sp. haplotypes Zhils1-18S and Zhils2-18S) resulted in identity scores > 99% over a length of 1700 alignment positions (Additional file 1: Table S1). This included the highly variable domains V2 and V4 according to the predicted secondary structure of the T. gondii reference sequence [27]. Taken together, these results strongly suggested that the 28S and 18S rRNA gene sequences deposited under KU341120 represented DNA of Sarcocystis sp. from Thailand. In contrast to the results on the 28S rRNA gene mentioned above, Sarcocystis sp. 18S rRNA haplotypes Zhils1-18S and Zhils2-18S showed the highest identity scores with the type sequences of S. zuoi (JQ029112, JQ029113), which made the two taxa indistinguishable at the 18S rRNA locus (Additional file 1: Table S1). However, one must consider here that the S. zuoi sequences were truncated [20] so that a potential divergence between the taxa may have been concealed.
We observed the greatest sequence variability in the ITS1-5.8S-ITS2 gene complex, whereby the most relevant (and available) sequences for comparison with Sarcocystis sp. from China (isolates Cl312-ITS/5.8S, Cl41-ITS/5.8S, Cl43-ITS/5.8S, Cl45-ITS/5.8S) were KU341120 and those from S. attenuati. Sequence identity of Sarcocystis sp. with KU341120 ranged between 92 and 93% (1268–1276 alignment positions), whereas this value was 89% for comparisons with various haplotypes of S. attenuati over 712–736 alignment positions (Additional file 1: Table S1). Although we had amplified various sequences of the ITS complex from S. zuoi in China, the amplicons only included ITS1 and 5.8S rRNA and are therefore not listed in Additional file 1: Table S1. However, we used the excised ITS1 of S. zuoi for phylogenetic analysis (below).
Four new sequences of the cox1 gene of Sarcocystis sp. (only one sequence deposited: PP033596, isolate 11526FR1-Cox1) were not only completely identical among each other but also highly similar to S. kani, S. scandentiborneensis, S. attenuati and S. zuoi (Additional file 1: Table S1): The new Sarcocystis sp. differed from S. kani in only three out of 927 homologous positions, in five out of 1333 bp compared to S. attenuati (MZ889673), in three out of 976 bp in the case of S. scandentiborneensis and in eight out of 978 bp compared to the newly sequenced cox1 gene of S. zuoi (PP033597, isolate SF1-Rm1011-Cox1). Its divergence compared to a recently discovered Sarcocystis from Rattus argentiventer in Sumatra (Sarcocystis sp. 2, ON989200) was higher with 20 out of 986 bp changes. Since all these taxa are considered part of the molecularly cryptic S. zuoi complex of species [6], the new results here show that the latter two species exhibited the greatest divergence at the cox1 locus within this group, although overall variability was relatively low. Furthermore, it confirmed for cox1 (as for the 28S rRNA) that Sarcocystis sp. isolated from Norway rats in China was different from the original S. zuoi isolated from the same intermediate host in China. Sequence divergence of Sarcocystis sp. compared to S. singaporensis and S. zamani from pythonid hosts was higher (Additional file 1: Table S1). The latter two species and the S. zuoi—complex are included in the so-called S1-lineage of snakehost Sarcocystis spp. in the 18S rRNA phylogenetic tree [35], whereby the cox1 gene tree topology is similar to the 18S rRNA except that members of S. zuoi complex form a polytomic clade [6]. The same polytomy also resulted after phylogenetic analysis of cox1 of the new Sarcocystis species, which is why we do not show the cox1 tree here. The new cox1 sequences contained the barcode area, which in the case of Sarcocystis sp. started with a complete, in the case of S. zuoi incomplete helix 1, as exemplified for S. scandentiborneensis [20].
Phylogenies of the 28S, 18S rRNA genes and ITS1 region
Phylogenetic reconstruction of the 28S rRNA gene tree was performed considering longer sequences (comparing about 1400 homologous positions) as well as shorter sequences to accommodate the samples from Borneo (400–500 sites). Maximum likelihood (ML) analysis with 29 longer sequences, including Sarcocystis sp. from Norway rats in China, sequence KU341120 isolated from R. tiomanicus in Thailand and S. zuoi from R. norvegicus in China revealed that the former two taxa clustered together while S. zuoi was a sister species included in a separate clade that also contained the sequences KU341118/19/21 from R. tiomanicus in Thailand (Fig. 4a). The Sarcocystis species from pythonid snakes (S. singaporensis, S. zamani) branched off basally from the monophyletic clade shared with the Sarcocystis spp. of interest. Sarcocystis pantherophisi, the only representative of the S2 lineage for this marker, branched off separately from the S1 cluster with high branch support. Subsequently, we used S. pantherophisi as outgroup in a smaller tree of 16 28S rRNA sequences that included the shorter amplicons of the samples from Borneo. As shown in Fig. 4b, all sequences of Sarcocystis sp. from China and KU341120 from Thailand formed a monophyletic clade that also included isolate E369-13-28S from M. whiteheadi in Borneo. Hence, all three geographic locations of the novel species were combined in one monophyletic clade with relatively high branch support, which was congruent with the high identity scores observed in the BLAST analyses (Additional file 1: Table S1). A sister clade contained the two sequences retrieved from the squirrel Sundasciurus lowii that clustered with S. scandentiborneensis (isolate E357.13-28S). Similar to the topology of the 28S rRNA tree based on longer sequences, S. zuoi clustered with sequences KU341118/19/21, apart from the new Sarcocystis species.
A different tree topology resulted after reconstruction of the 18S rRNA tree (albeit after inclusion of additional members of the S. zuoi complex, such as S. attenuati and S. kani), where Sarcocystis sp. and sequence KU341120 clustered with S. zuoi in a fashion that could not unequivocally distinguish between these taxa: Character-based ML analysis produced a polytomic clade with very short branch lengths, while the distance-based ME algorithm resulted in a dichotomic topology of this clade which suggested that isolate KU341120 and Sarcocystis sp. assumed an ancestral position relative to S. zuoi (Fig. 5): There was no speciation time/branching (sensu [39]) of Sarcocystis sp. after its split from S. zuoi. Similar to the 28S rRNA tree, the clade containing the members of the S. zuoi complex showed good bootstrap support and species of the S2 snakehost lineage (e.g. S. pantherophisi, Sarcocystis nesbitti, S. atheridis) were part of a separate cluster of sequences.
Since the 18S rRNA tree could not unequivocally distinguish between S. muricoelognathis and S. zuoi, and the cox1 locus was largely uninformative regarding the Sarcocystis species of interest (Additional file 1: Table S1), we also reconstructed the phylogeny of the ITS1 locus (after excision from the ITS1-5.8S-ITS2 complex deposited under OR977565-68). We hoped to find additional evidence for the divergence between Sarcocystis sp. and S. zuoi, which apparently use the same intermediate hosts in the wild in China. Although alignment of the ITS1 sequences was demanding because of their high variability (i.e. relatively few conserved positions shared between species), secondary structure-guided alignment and subsequent ML analysis resulted in a tree topology that was well in agreement with the 28S rRNA tree: The ITS1 sequences of Sarcocystis sp. and Sarcocystis isolate KU341120 clustered together, while the newly sequenced ITS1 of S. zuoi (OR990322-25) formed a separate clade with high branch support (Fig. 6). Similarly, S. singaporensis and S. zamani branched off basally from the shared S1 clade, and S. pantherophisi was located on a separate branch.
Taxonomic summary of Sarcocystis muricoelognathis sp. nov.
(Sarcocystis sp.: Figs. 1, 2, 3, 4, 5, 6, Additional file 1: Table S1).
Diagnosis
Oocyst containing two sporocysts with four sporozoites each; sporocysts measuring on average 10.1–10.3 × 7.2–7.5 µm; sarcocysts microscopic with striated cyst wall in fresh preparations (similar to S. singaporensis but cyst wall protrusions not resting on stalks), cigar-shaped with slightly pointed ends in older cysts; in skeletal musculature except heart; native sarcocysts in rats infected experimentally for 5 months measured 400–700 × 70–90 µm, in wild-caught rats 770–1580 × 87–123 μm; cyst wall with palisade-like, villar protrusions, 8 to 11 µm long and 2–4 µm wide at their base; protrusions at the tips of the sarcocyst often longer, 10–16 µm; banana-shaped, live cystozoites on average 5.9–6.5 µm long and 2.1–2.2 µm in diameter; ultrastructurally, villar protrusions palisade-like, attached upright to the cyst (if not bent by preparation), sometimes slightly broadened towards the apical ends, delimited by a knob-like lining and filled with fine, electron-light granular material; cystozoites with only two rhoptries throughout distribution range.
Type natural definitive hosts
Coelognathus radiatus (holotype), Coelognathus flavolineatus (paratype).
Type experimental intermediate host
Rattus norvegicus (holotype)
Type natural intermediate hosts
Maxomys whiteheadi (paratype); Rattus norvegicus (paratype)
Type localities
Holotype: central Bangkok, Thailand; paratypes: Chumphon, southern Thailand; Poring Hot Springs, Northern Borneo, Sabah, Malaysia; Anning and Nanjian Prefectures, China.
Other locations and natural hosts not included in type series
Sarcocysts with the same ultrastructure recorded from Rattus argentiventer and Bandicota indica in Thailand [2] and reported from Rattus rattus in Singapore [15].
Etymology
Species name in reference to the intermediate and definitive hosts.
Type specimens deposited
Holotype deposited at the Department of Parasitology, University of Hohenheim, Emil-Wolff-Str. 34, 70599 Stuttgart-Hohenheim, Germany, accession number DP/R41C961764: Resin-embedded sarcocyst from Rattus norvegicus infected with sporocysts from Coelognathus radiatus; paratype from Borneo deposited at the Sabah Parks Museum, Kinabalu Park, Sabah, Malaysia, accession number SP/EAC/00102: sarcocyst in musculature from Maxomys whiteheadi; paratypes from China deposited at the Zoological Specimen Museum of Yunnan University, Kunming, accession number Prot202302: formalin-fixed tissue containing sarcocysts as well as photomicrographs.
Sequences deposited at GenBank
18S rRNA: OR976520, OR976521; 28S rRNA: isolates from China OR979639, OR979640, OR979641, OR979642, and Borneo OR979643; ITS1-5.8S-ITS2: OR977565, OR977566, OR977567, OR977568; cox1: PP033596.
ZooBank Registration
This article has been submitted to ZooBank under the Life Science Identifier (LSID) urn:lsid:zoobank.org:pub:E86A1FEF-6189-4E28-A627-11ABFAA635AC. The LSID for the new species name Sarcocystis muricoelognathis is urn:lsid:zoobank.org:act:000AF647-D4AC-4609-914D-61C1449A2DAA.
Discussion
Taxonomic remarks
Since the ultrastructure of the sarcocyst wall is in many cases a useful criterion for distinguishing between species of Sarcocystis [40], we show here that S. muricoelognathis is clearly distinct from two similar species, S. zuoi and S. murinotechis. In the case of S. zuoi, its villar protrusions of the cyst wall are bent by a certain angle, which is also visible in live sarcocysts when the surrounding host tissue has been removed, and exhibit a reticulate base [11]. In contrast, S. muricoelognathis shows upright protrusions with a largely non-reticulate base. Compared to S. murinotechis from Australia [7], the two former species have markedly longer protrusions (roughly twice as long). Although only two electron micrographs of the sarcocysts of S. murinotechis exist in the original species description, the accompanying text characterizes them as possessing “short, broad protrusions” that were about 6 µm in length [7]. Relatively short protrusions are also visible in a histological section of a sarcocyst of S. cf. murinotechis from Rattus lutreolus [18] supporting the observations on S. murinotechis published by these authors 2 years later. All in all, we propose that the sarcocyst ultrastructure of S. muricoelognathis could be categorized under group 10a of the common cyst wall classification system for Sarcocystis [40]. Additionally, the sizes of the transmission stages, sporocysts and cystozoites, were significantly different among the three species discussed above. Notably, sporocysts of S. murinotechis exhibited a shape index different from S. muricoelognathis and S. zuoi: Size differences of sporocysts can be regarded a suitable indicator for the potential presence of different Sarcocystis species [41]. Regarding intraspecific variability, cystozoites of S. muricoelognathis were remarkably constant in size independent of geographic location and intermediate host and consistently showed two rhoptries at the ultrastructural level. These characteristics underscore the uniqueness of S. muricoelognathis compared to other species with striated cyst walls, which we have discussed previously [20]. At the light microscopic level, sarcocysts of S. muricoelognathis may be easily confused with S. singaporensis, because both species are microscopic, show a striated cyst wall with relatively long protrusions, exhibit pointed sarcocyst ends and share similar geographic distributions and rodent hosts. However, the cyst wall protrusions of the latter species rest on typical stalks, which appear as a thin, basal layer of the striated wall at lower magnification in live and histological preparations [2]. Molecularly, S. muricoelognathis can be unequivocally identified using its 28S rRNA and ITS1 sequences, whereby the partial 28S rRNA gene (≥ 600 bp) is already sufficient for differential diagnosis if the variable expansion segments D1 and D2 are included (as discussed below).
In a taxonomic context, we think we provided substantial molecular evidence for the proposition that the rRNA sequences from Thailand deposited under accession number KU341120 [19] do not belong to S. zuoi but represent S. muricoelognathis. Although a high degree of sequence identity between partial sequences alone may be insufficient for species discrimination, we must add here that the 28S rRNA gene sequence of S. muricoelognathis from Borneo included the variable expansion segments D1 and D2 of the rRNA molecule according to the predicted secondary structure of the reference sequence of T. gondii (L25635, [27]). Especially, domain D2 has been shown to be highly variable [42]. Moreover, the overlapping parts of the alignments of partial sequence KU341120 with the (almost) full-length 28S rRNA sequences of S. muricoelognathis from China included expansion segments D1–D6, which, in view of an identity score of 99%, strongly suggested that these sequences belong to the same species. Furthermore, both taxa clustered together in the ITS1 phylogenetic tree, which is increasingly used for discriminating between closely related Sarcocystis spp. (e.g. [43]) and shows much greater evolutionary rates than other rRNA markers [28]. Given that the sequences under KU341120 have been included in major rRNA databases under the species name ‘Sarcocystis zuoi’ and are frequently cited in the context of Sarcocystis of rodents (e.g. [44, 45]), it will probably take some time until any annotation in line with our proposal could be expected. Moreover, based on the observations here and a previous study [6], it can be safely assumed that the sequences KU341118/19/21 also do not represent S. zuoi but belong to a different species, because they clustered in a clade apart from the original S. zuoi sequences.
Close phylogenetic relationships between S. muricoelognathis and its sister species
Sarcocystis muricoelognathis is a new member of what we call the S. zuoi complex, a group of Sarcocystis species cycling between colubrid snakes and small mammals in eastern Asia that shows remarkably close genetic relationships but significant morphological disparity regarding sarcocyst ultrastructure as well as the sizes of its transmission stages [6]. Close genetic relatedness was especially apparent among S. muricoelognathis, S. attenuati [23], S. kani [6] and S. scandentiborneensis [20], from which the former differed in ≤ 5-bp changes within the partial mitochondrial cox1 gene. The cox1 gene sequences included the barcode fragment, which is used globally for discrimination between species of the tree of life [46]. Although the newly sequenced cox1 gene of S. zuoi, the indicator species of this complex, showed a few additional bp changes, its identity score also reached 99% compared to the species above. This suggests incomplete lineage sorting (ILS) at this locus and a recent common ancestor of these taxa. Mitochondrial DNA shows a relatively rapid rate of mutation, which makes it suitable as a marker for more recent evolutionary history of natural populations [47, 48]. In contrast, the 28S rRNA and ITS1 trees could distinguish well between S. muricoelognathis and its sister species, which was indicated by a mostly dichotomic tree resolution and absence of ILS, which can be recognized by polyphyly or paraphyly in fast evolving species [49, 50]. Although the 18S rRNA gene was also useful for discriminating between various taxa of the S. zuoi complex, it could not unequivocally resolve the divergence between S. zuoi and S. muricoelognathis in this study for the reason outlined earlier. In the future, a full-length 18S rRNA of S. zuoi should be sequenced to aid in the further study of this group of Sarcocystis spp. The fast-evolving ITS1 marker (e.g. compare scale bars of Figs. 5 and 6) produced results congruent with the 28S and 18S rRNA gene trees. We can confirm that the alignment of ITS sequences was challenging because of their high variability and often incorrect sequence boundaries in public databases [51]. We are confident that careful determination of starting and end points of ITS1 and using a secondary structure-guided sequence alignment (because ITS folds into stems, bulges and loops; [52]) helped find the correct tree.
The close genetic relatedness among the species of the S. zuoi complex as apparent in the cox1 phylogeny [6] and revealed here for S. muricoelognathis as well was also reflected in the relatively short branch lengths of this clade in the 18S rRNA gene tree. This contrasted with longer branch lengths of the species of the S2 snakehost lineage (sensu [35]), such as S. pantherophisi, S. atheridis or S. nesbitti. Differences in branch lengths can indicate variable evolutionary rates among the coccidia [34] but could also be a result of different evolutionary histories [53]. In view of this, we hypothesized that the Sarcocystis species of the S. zuoi complex may undergo radiation that followed the documented recent radiations of their intermediate and definitive hosts in Asia [6].
Host shifting of Sarcocystis between commensal and forest rodents in the biodiversity hotspot Borneo
Remarkably, S. muricoelognathis exhibits a relatively wide distribution range in Asia (China, Thailand, northern Borneo and probably also Singapore) and apparently utilizes forest and commensal rodents as intermediate hosts. A wide distribution range is in line with other Sarcocystis species from colubrid snakes, notably, with S. clethrionomyelaphis, which was originally isolated in Germany from Zamenis longissimus [36] and later in China from Elaphe taeniura [37]. Also, the geographical range of S. zuoi was previously already assumed to include China, Thailand, Malaysia and Japan [11, 19, 54, 55], respectively. However, we show here that the four 18S rRNA gene sequences from Thailand previously designated as ‘S. zuoi’ are not conspecific with S. zuoi as discussed above. Furthermore, other sequences of the same gene with the same species designation collected in Malaysia and Japan may not be S. zuoi, because they were relatively short and not accompanied by a morphological description of a sarcocyst [54, 55]. Hence, based on current evidence S. zuoi records are from China only; however, that does not exclude the possibility that this species has a wider distribution than currently anticipated.
The island of Borneo with its dipterocarp tropical forests is a global biodiversity hotspot for mammalian species [56]. In a previous study, we showed that S. scandentiborneensis is a parasite infecting treeshrews in Borneo [20]. Here, we found first molecular and morphological indications that the intermediate host range of S. scandentiborneensis includes the rat Niviventer cremoriventer and the squirrel Sundasciurus lowii. Although these findings require further confirmation, especially in the case of Niviventer, they hint at the possibility that the intermediate host range of this species in its forest habitat might be larger than previously anticipated. Our observation that S. muricoelognathis occurs in commensal and forest rodents points in the same direction; we observed the new Sarcocystis species in the rat Maxomys whiteheadi, which is frequently found in primary and logged forest in Borneo but exhibits limited habitat overlap with commensal Rattus spp. [57]. While S. muricoelognathis is a prevalent encounter in commensal rats in China (this study), Thailand [2] and on Singaporean islands [15], the newly reported Sarcocystis species with snake-rodent life cycle in forest habitats question the origin of associations with rodents and the possible history of host shifting events. In this context, the study of Paperna et al. [15] is interesting insofar as they showed that infection rates with S. singaporensis of the commensal rodent Rattus rattus were higher in forest habitats on islands that were less disturbed by human activity, i.e. that were more likely to harbor the pythonid definitive host at densities that could maintain the parasite’s life cycle. A follow-up study in Singapore about 12 years later confirmed the trend that parasite transmission was negatively affected by expanding urbanization [58]. Furthermore, a recent phylogenetic study on the evolution of dietary preferences of colubrid snakes found that Coelognathus flavolineatus and C. radiatus are prototypes of a dietary generalist that occupied tropical forests in Asia [59]. The evolutionary history of rat snakes traces back to the forests of Southeast Asia, where their ancestors originated around 37 million years ago [60, 61]. Moreover, although Sarcocystis singaporensis is one of the most prevalent and widely distributed coccidian species in commensal rats of the Indo-Australian Archipelago [35], its intermediate host spectrum on Java and Sulawesi in Indonesia appears to be greatly extended to a variety of forest rodents [10].
This growing body of insight suggests that transmission cycles between snakes and rodents are closed in both forest and anthropogenically modified habitats, questioning whether commensal rodents or those confined to natural forest habitats were the original intermediate hosts. In this context, it was hypothesized that infection of phylogenetically older intermediate hosts like shrews and treeshrews by phylogenetically younger species of the S. zuoi complex (i.e. S. attenuati and S. scandentiborneensis, respectively) might be the result of host shifting or adaptive processes rather than co-speciation [6]. Our finding that S. scandentiborneensis also likely occurs in sciurid rodents highlights the possibility that such a host shift could have occurred from the Rodentia to the Scandentia. For eimeriid coccidia infecting rodents, the occurrence of repeated host switches among closely related genetic lineages has been established [62]. Furthermore, the paradigm of Sarcocystis spp. being highly host-specific for their intermediate hosts might require a more differentiated analysis, since deviations from this generalization have been observed for Sarcocystis spp. infecting wild ruminants and farm animals [63, 64], members of the lineage clustering with Sarcocystis neurona (e.g. S. canis, [43]) and, early on, Sarcocystis spp. with snake-rodent life cycle [65].
Conclusions
Sarcocystis muricoelognathis sp. nov. cycles between colubrid snakes as definitive and forest and commensal rodents as intermediate hosts, further increasing the number of species of the so-called S. zuoi complex, members of which appear to have expanded intermediate host ranges. Although we anticipate that these host associations are likely only snapshots of the true associations, our results provide evidence of broad geographic distributions of rodent-associated Sarcocystis and host shifts between commensal and forest small mammals. Future work is required to better understand the possible eco-evolutionary pathways of host shifting of Sarcocystis species, which also deserves attention from a contemporary conservation and disease control perspective.
Availability of data and materials
All relevant data of this study are contained within the manuscript, while type materials of the new Sarcocystis species have been deposited in the institutions as indicated in the taxonomic summary section. All new DNA sequences are publicly available at GenBank under the following accession numbers: 18S rRNA: Sarcocystis sp. isolates Zhils1-18S, Zhils2-18S (OR976520, OR976521); 28S rRNA: Sarcocystis sp. isolates Rn-cl1.3-28S, Rn-cl2.3-28S, Rn-cl3.1-28S, Rn-cl3.3-28S (OR979639, OR979640, OR979641, OR979642); Sarcocystis sp. isolate E369-13-28S (OR979643); Sarcocystis cf. scandentiborneensis isolates E388-13A-28S, E367-13-28S, E346-13-28S (OR979644, OR979645, OR979646); Sarcocystis scandentiborneensis isolate E357-13-28S (OR979647); Sarcocystis zuoi isolate Rn-cl2.3-28S (OR979648); ITS1-5.8S rRNA-ITS2: Sarcocystis sp. isolates Cl312-ITS-5.8S, Cl41-ITS-5.8S, Cl43-ITS-5.8S, Cl45-ITS-5.8S (OR977565, OR977566, OR977567, OR977568); ITS1: Sarcocystis zuoi isolates cl1.2, cl2.1, cl2.2, cl3.3 (OR990322, OR990323, OR990324, OR990325); cox1: Sarcocystis sp. isolate 11526FR1-Cox1 (PP033596); Sarcocystis zuoi isolate SF1-Rm1011-Cox1 (PP033597).
Abbreviations
- ANOVA:
-
Analysis of variance
- cox1 :
-
Mitochondrial Cytochrome oxidase subunit 1 gene
- ES:
-
Eukaryotic expansion segments of the rRNA
- ITS:
-
Internal transcribed spacer region
- ME:
-
Minimum evolution criterion
- ML:
-
Maximum likelihood criterion
- PCR:
-
Polymerase chain reaction
- s.d.:
-
Standard deviation
- s.e.:
-
Standard error of the mean
- SD rats:
-
Sprague-Dawley rats
- 18S, 5.8S and, 28 S rRNA :
-
Nuclear 18S, 5.8S and 28S ribosomal RNA genes
References
Ashford RW. Sarcocystis cymruensis n. sp., a parasite of rats Rattus norvegicus and cats Felis catus. Ann Trop Med Parasitol. 1978;72:37–43.
Jäkel T, Khoprasert Y, Sorger I, Kliemt D, Seehabutr V, Suasa-ard K, et al. Sarcosporidiasis in rodents from Thailand. J Wildl Dis. 1997;33:860–7.
Hu JJ, Liao JY, Meng Y, Guo YM, Chen XW, Zuo YX. Identification of Sarcocystis cymruensis in wild Rattus flavipectus and Rattus norvegicus from Peoples Republic of China and its transmission to rats and cats. J Parasitol. 2011;97:421–4.
Antunes Murata FH, Cerqueira-Cezar CK, Thompson PC, Tiwari K, Mowery JD, Verma SK, et al. Sarcocystis cymruensis: discovery in Western Hemisphere in the Brown rat (Rattus norvegicus) from Grenada, West Indies: redescription, molecular characterization, and transmission to IFN-gamma gene knockout mice via sporocysts from experimentally infected domestic cat (Felis catus). Parasitol Res. 2018;117:1195–204.
Zeng H, Guo Y, Ma C, Deng S, Hu J, Zhang Y. Redescription and molecular characterization of sarcocysts of Sarcocystis cymruensis from Norway rats (Rattus norvegicus) and Sarcocystis ratti from black rats (R. rattus) in China. Parasitol Res. 2020;119:3785–91.
Jäkel T, Raisch L, Richter S, Wirth M, Birenbaum D, Ginting S, et al. Morphological and molecular phylogenetic characterization of Sarcocystis kani sp. nov. and other novel, closely related Sarcocystis spp. infecting small mammals and colubrid snakes in Asia. Int J Parasitol Parasites Wildl. 2023;22:184–98.
Munday BL, Mason RW. Sarcocystis and related organisms in Australian wildlife: III. Sarcocystis murinotechis sp.n. life cycle in rats (Rattus, Pseudomys and Mastocomys spp.) and tiger snakes (Notechis ater). J Wildl Dis. 1980;16:83–8.
Prakas P, Kirillova V, Gavarane I, Gravele E, Butkauskas D, Rudaityte-Lukosiene E, et al. Morphological and molecular description of Sarcocystis ratti n. sp. from the black rat (Rattus rattus) in Latvia. Parasitol Res. 2019;118:2689–94.
Beaver PC, Maleckar JR. Sarcocystis singaporensis Zaman and Colley, (1975) 1976, Sarcocystis villivillosi sp. n., and Sarcocystis zamani sp. n.: development, morphology, and persistence in the laboratory rat, Rattus norvegicus. J Parasitol. 1981;67:241–56.
O’Donoghue PJ, Watts CH, Dixon BR. Ultrastructure of Sarcocystis spp. (Protozoa: Apicomplexa) in rodents from North Sulawesi and West Java, Indonesia. J Wildl Dis. 1987;23:225–32.
Hu JJ, Meng Y, Guo YM, Liao JY, Song JL. Completion of the life cycle of Sarcocystis zuoi, a parasite from the Norway rat, Rattus norvegicus. J Parasitol. 2012;98:550–3.
Munday BL. An isosporan parasite of masked owls producing sarcocysts in rats. J Wildl Dis. 1983;19:146–7.
Grikieniene J, Arnastauskiene T, Kutkiene L. On some disregarded ways of sarcosporidians circulation and remarks about systematics of the genus Sarcocystis Lankester, 1882 with the description of the new species from rodents. Ekologija. 1993;1:16–24.
Slapeta JR, Modry D, Votypka J, Jirku M, Koudela B, Lukes J. Multiple origin of the dihomoxenous life cycle in sarcosporidia. Int J Parasitol. 2001;31:413–7.
Paperna I, Peh KSH, Martelli P, Koh LP, Sodhi NS. Factors affecting Sarcocystis infection of rats on small tropical islands. Ecol Res. 2004;19:475–83.
Lee FCH. Finding Sarcocystis spp. on the Tioman Island: 28S rRNA gene next-generation sequencing reveals nine new Sarcocystis species. J Water Health. 2019;17:416–27.
Uetz P, Freed P, Aguilar R, Reyes F, Hošek J. The Reptile Database. http://www.reptile-database.org. 2022.
Munday BL, Mason RW, Hartley WJ, Presidente PJ, Obendorf D. Sarcocystis and related organisms in Australian wildlife: I. Survey findings in mammals. J Wildl Dis. 1978;14:417–33.
Watthanakaiwan V, Sukmak M, Hamarit K, Kaolim N, Wajjwalku W, Muangkram Y. Molecular characterization of the ribosomal DNA unit of Sarcocystis singaporensis, Sarcocystis zamani and Sarcocystis zuoi from rodents in Thailand. J Vet Med Sci. 2017;79:1412–8.
Ortega Pérez P, Wibbelt G, Brinkmann A, Galindo Puentes JA, Tuh FYY, Lakim MB, et al. Description of Sarcocystis scandentiborneensis sp. nov. from treeshrews (Tupaia minor, T. tana) in northern Borneo with annotations on the utility of COI and 18S rDNA sequences for species delineation. Int J Parasitol Parasites Wildl. 2020;12:220–31.
Aplin KP, Brown PR, Jacob J, Krebs CJ, Singleton GR. Field methods for rodent studies in Asia and the Indo-Pacific. Canberra: ACIAR; 2003.
Musser G, Carleton M. Superfamily muroidea. In: Wilson DE, Reeder DM, editors. Mammal species of the world: a taxonomic and geographic reference. Baltimore: Johns Hopkins University; 2005. p. 894–1531.
Hu J, Sun J, Guo Y, Zeng H, Zhang Y, Tao J. Infection of the Asian gray shrew Crocidura attenuata (Insectivora: Soricidae) with Sarcocystis attenuati n. sp. (Apicomplexa: Sarcocystidae) in China. Parasit Vectors. 2022;15:13.
Jäkel T. Cyclic transmission of Sarcocystis gerbilliechis n. sp. by the Arabian saw-scaled viper, Echis coloratus, to rodents of the subfamily Gerbillinae. J Parasitol. 1995;81:626–31.
Prakas P, Oksanen A, Butkauskas D, Sruoga A, Kutkiene L, Svazas S, et al. Identification and intraspecific genetic diversity of Sarcocystis rileyi from ducks, Anas spp., Lithuania and Finland. J Parasitol. 2014;100:657–61.
Tamura K, Stecher G, Peterson D, Filipski A, Kumar S. MEGA6: molecular evolutionary genetics analysis version 6.0. Mol Biol Evol. 2013;30:2725–9.
Gagnon S, Bourbeau D, Levesque RC. Secondary structures and features of the 18S, 5.8S and 26S ribosomal RNAs from the Apicomplexan parasite Toxoplasma gondii. Gene. 1996;173:129–35.
Bass D, Richards TA, Matthai L, Marsh V, Cavalier-Smith T. DNA evidence for global dispersal and probable endemicity of protozoa. BMC Evol Biol. 2007;7:162.
Clark K, Karsch-Mizrachi I, Lipman DJ, Ostell J, Sayers EW. GenBank. Nucleic Acids Res. 2017;44:D67-72.
Tommaso PD, Moretti S, Xenarios I, Orobitg M, Montanyola A, Chang JM, et al. T-Coffee: a web server for the multiple sequence alignment of protein and RNA sequences using structural information and homology extension. Nucleic Acids Res. 2011;39:W13-7.
Kumar S, Stecher G, Li M, Knyaz C, Tamura K. MEGA X: molecular evolutionary genetics analysis across computing platforms. Mol Biol Evol. 2018;35:1547–9.
Rzhetsky A, Nei M. A simple method for estimating and testing minimum evolution trees. Mol Biol Evol. 1992;9:945–67.
Tamura K, Kumar S. Evolutionary distance estimation under heterogeneous substitution pattern among lineages. Mol Biol Evol. 2002;19:1727–36.
Morrison DA, Bornstein S, Thebo P, Wernery U, Kinne J, Mattsson JG. The current status of the small subunit rRNA phylogeny of the coccidia (Sporozoa). Int J Parasitol. 2004;34:501–14.
Wassermann M, Raisch L, Lyons JA, Natusch DJD, Richter S, Wirth M, et al. Examination of Sarcocystis spp. of giant snakes from Australia and Southeast Asia confirms presence of a known pathogen—Sarcocystis nesbitti. PLoS ONE. 2017;12:e0187984.
Matuschka F-R. Sarcocystis clethrionomyelaphis n. sp. from snakes of the genus Elaphe and different voles of the family Arvicolidae. J Parasitol. 1986;72:226–31.
Hu JJ, Liu TT, Liu Q, Esch GW, Chen JQ. Sarcocystis clethrionomyelaphis Matuschka, 1986 (Apicomplexa: Sarcocystidae) infecting the large oriental vole Eothenomys miletus (Thomas) (Cricetidae: Microtinae) and its phylogenetic relationships with other species of Sarcocystis Lankester, 1882. Syst Parasitol. 2015;91:273–9.
Hu JJ, Liu Q, Yang YF, Esch GW, Guo YM, Zou FC. Sarcocystis eothenomysi n. sp. (Apicomplexa: Sarcocystidae) from the large oriental vole Eothenomys miletus (Thomas) (Cricetidae: Microtinae) from Anning, China. Syst Parasitol. 2014;89:73–81.
Hedges SB, Marin J, Suleski M, Paymer M, Kumar S. Tree of life reveals clock-like speciation and diversification. Mol Biol Evol. 2015;32:835–45.
Dubey JP, Calero-Bernal R, Rosenthal BM, Speer CA, Fayer R. Sarcocystosis of animals and humans. Boca Raton: CRC Press; 2016.
McAllister CT, Upton SJ, Trauth SE, Dixon JR. Coccidian parasites (Apicomplexa) from snakes in the southcentral and southwestern United States: new host and geographic records. J Parasitol. 1995;81:63–8.
Mugridge NB, Morrison DA, Jäkel T, Heckeroth AR, Tenter AM, Johnson AM. Effects of sequence alignment and structural domains of ribosomal DNA on phylogeny reconstruction for the protozoan family Sarcocystidae. Mol Biol Evol. 2000;17:1842–53.
Leger JS, Chen Y, Sakamaki K, Mena A, Raverty SA, Rotstein D, et al. Fatal hepatic sarcocystosis in three captive and one free-ranging pinniped. Int J Parasitol Parasites Wildl. 2023;22:249–54.
Canova V, Helman E, Robles MDR, Abba AM, More G. First report of Sarcocystis spp. (Apicomplexa, Sarcocystidae) in Lagostomus maximus (Desmarest, 1917) (Rodentia, Chinchillidae) in Argentina. Int J Parasitol Parasites Wildl. 2023;20:180–6.
Prakas P, Stirke V, Sneideris D, Rakauskaite P, Butkauskas D, Balciauskas L. Protozoan parasites of Sarcocystis spp. in rodents from commercial Orchards. Animals. 2023. https://doi.org/10.3390/ani13132087.
Pentinsaari M, Salmela H, Mutanen M, Roslin T. Molecular evolution of a widely-adopted taxonomic marker (COI) across the animal tree of life. Sci Rep. 2016;6:35275.
Brown WM, George M Jr, Wilson AC. Rapid evolution of animal mitochondrial DNA. Proc Natl Acad Sci USA. 1979;76:1967–71.
Norman JA, Blackmore CJ, Rourke M, Christidis L. Effects of mitochondrial DNA rate variation on reconstruction of Pleistocene demographic history in a social avian species, Pomatostomus superciliosus. PLoS ONE. 2014;9:e106267.
Belfiore NM, Liu L, Moritz C. Multilocus phylogenetics of a rapid radiation in the genus Thomomys (Rodentia: Geomyidae). Syst Biol. 2008;57:294–310.
Maddison WP, Knowles LL. Inferring phylogeny despite incomplete lineage sorting. Syst Biol. 2006;55:21–30.
Bengtsson-Palme J, Ryberg M, Hartmann M, Branco S, Wang Z, Godhe A, et al. Improved software detection and extraction of ITS1 and ITS2 from ribosomal ITS sequences of fungi and other eukaryotes for analysis of environmental sequencing data. Methods Ecol Evol. 2013;4:914–9.
Rampersad SN. ITS1, 5.8S and ITS2 secondary structure modelling for intra-specific differentiation among species of the Colletotrichum gloeosporioides sensu lato species complex. SpringerPlus. 2014;3:684.
Molnar K, Ostoros G, Dunams-Morel D, Rosenthal BM. Eimeria that infect fish are diverse and are related to, but distinct from, those that infect terrestrial vertebrates. Infect Genet Evol. 2012;12:1810–5.
Lau YL, Chang PY, Subramaniam V, Ng YH, Mahmud R, Ahmad AF, et al. Genetic assemblage of Sarcocystis spp. In Malaysian snakes. Parasit Vectors. 2013;6:257.
Abe N, Matsubara K, Tamukai K, Miwa Y, Takami K. Molecular evidence of Sarcocystis species in captive snakes in Japan. Parasitol Res. 2015;114:3175–9.
Struebig MJ, Wilting A, Metcalfe K, Kramer-Schadt S. Targeted conservation to safeguard a biodiversity hotspot from climate and land-cover change. Curr Biol. 2015;25:372–8.
Wells K, Kalko EK, Lakim MB, Pfeiffer M. Effects of rain forest logging on species richness and assemblage composition of small mammals in Southeast Asia. J Biogeogr. 2007;34:1087–99.
Devan-Song A, Luz S, Mathew A, Low MR, Bickford DP. Pythons, parasites, and pests: anthropogenic impacts on Sarcocystis (Sarcocystidae) transmission in a multi-host system. Biotropica. 2017;49:706–15.
Barends JM, Maritz B. Dietary specialization and habitat shifts in a clade of Afro-Asian colubrid snakes (Colubridae: Colubrinae). Ichthyol Herpetol. 2022;110:278–91.
Burbrink FT, Lawson R. How and when did old world ratsnakes disperse into the new world? Mol Phylogenet Evol. 2007;43:173–89.
Rage J-C. Fossil record. In: Seigel Richard A, Collins Joseph T, Novak Susan S, editors. Snakes: ecology and evolutionary biology. New York: McGraw-Hill Publishing; 1987. p. 51–76.
Macova A, Hoblikova A, Hypsa V, Stanko M, Martinu J, Kvicerova J. Mysteries of host switching: diversification and host specificity in rodent-coccidia associations. Mol Phylogenet Evol. 2018;127:179–89.
Dahlgren SS, Gjerde B. Molecular characterization of five Sarcocystis species in red deer (Cervus elaphus), including Sarcocystis hjorti n. sp., reveals that these species are not intermediate host specific. Parasitology. 2010;137:815–40.
Marandykina-Prakiene A, Butkauskas D, Gudiskis N, Juozaityte-Ngugu E, Bagdonaite DL, Kirjusina M, et al. Sarcocystis species richness in sheep and goats from Lithuania. Vet Sci. 2013. https://doi.org/10.3390/vetsci10080520.
Matuschka FR. Reptiles as intermediate and/or final hosts of Sarcosporidia. Parasitol Res. 1987;73:22–32.
Barta JR, Martin DS, Liberator PA, Dashkevicz M, Anderson JW, Feighner SD, et al. Phylogenetic relationships among eight Eimeria species infecting domestic fowl inferred using complete small subunit ribosomal DNA sequences. J Parasitol. 1997;83:262–71.
Fischer S, Odening K. Characterization of bovine Sarcocystis species by analysis of their 18S ribosomal DNA sequences. J Parasitol. 1998;84:50–4.
Fenger CK, Granstrom DE, Langemeier JL, Stamper S, Donahue JM, Patterson JS, et al. Identification of opossums (Didelphis virginiana) as the putative definitive host of Sarcocystis neurona. J Parasitol. 1995;81:916–9.
Mugridge NB, Morrison DA, Johnson AM, Luton K, Dubey JP, Votypka J, et al. Phylogenetic relationships of the genus Frenkelia: a review of its history and new knowledge gained from comparison of large subunit ribosomal ribonucleic acid gene sequences. Int J Parasitol. 1999;29:957–72.
Acknowledgements
We thank the Sabah Biodiversity Council (Majlis Biodiversity Sabah) for their support in obtaining research permits. Regarding the work in Thailand, we are grateful to the team of the Agricultural Zoology Research Group at the Department of Agriculture in Bangkok for their permission to use their laboratory for microscopic examinations and support in animal experimentation. In addition, the excellent technical assistance of Doris Krummow, Dagmar Viertel and Susanne Auls is much appreciated.
Funding
The study was supported by partial funding from the German Federal Ministry of Economic Cooperation and Development (BMZ) (Project No. 2002.2156.4; TJ for German International Cooperation), the Natural Sciences Foundation of China (Grant 32260119), the Science and Technology Planning Project in Key Areas of Yunnan Province, China (Grant 202303AP140004) and Graduate Innovation Project of Yunnan University (Grant ZC-23234473).
Author information
Authors and Affiliations
Contributions
JH and TJ conceived the study, analyzed the data and drafted the manuscript. TQ and PP conducted the molecular experiments and morphological analyses, while TJ contributed to the phylogenetic analyses. TQ, GW and TJ performed the ultrastructural studies; YK and SG performed the animal experimentation. ML, KW, JH and TJ did the field sampling. All authors reviewed and approved the final manuscript.
Corresponding authors
Ethics declarations
Ethics approval and consent to participate
For sampling in Borneo, biological resource access and export permits were issued by the Sabah Biodiversity Centre (JKM-MBS.1000–2/2(35), JKM-MBS.1000–2/2(63)). Collecting specimens in China was approved by the Animal Ethics Committee of Yunnan University (permission number ynucae20180002). In Thailand, animal experimentation was permitted under the bilateral collaboration agreement between the Department of Agriculture and German International Cooperation (GIZ) (Project number 2002.2156.4). Animal experimentation was performed according to the ‘Ethical Principles and Guidelines for the Use of Animals’ as published by the National Research Council of Thailand (NRC 1999), the ‘Guide for the Care and Use of Laboratory Animals’ of the National Academy of Sciences USA (1996), and guidelines by the American Veterinary Medical Association (https://www.avma.org).
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.
Pairwise sequence comparisons using the novel sequences of Sarcocystis sp., S. scandentiborneensis and S. zuoi and their association with sarcocyst morphology where possible (remarks).
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
Qin, T., Ortega-Perez, P., Wibbelt, G. et al. A cyst-forming coccidian with large geographical range infecting forest and commensal rodents: Sarcocystis muricoelognathis sp. nov.. Parasites Vectors 17, 135 (2024). https://doi.org/10.1186/s13071-024-06230-8
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s13071-024-06230-8