PLoS PathogensEvolutionary BiologyPeer-reviewed journal article
Misleading inference of schistosome epidemiology from ribosomal internal transcribed spacer (ITS) and mitochondrial DNA
E.E. Enabuele, Roy N. PlattiD, Ehizogie Egbeobauwaye AdeyemiiD, M.S.O. AisieniD, Oluwaremilekun Grace AjakayeiD, Mahmud Umar Ali, Ebube Charles AmaechiiD, Tolulope Ebenezer AtalabiiD, Timothy AutaiD, Oluwaseun Bunmi AwosoluiD, Adamu Galadima Dagona, Omoyemwen Edo-Taiwo, Ejikeugwu ChikaiD, C. IgbeneghuiD, Victor Stephen NjomiD, Marian Onwude AgbuguiiD, Mary-Kate N. Orji, Funso Olu Peter Oyinloye, Esther Oyemade, Habibat J. Ozemoka, Christopher R. Pam, Uchenna Iyioku UgahiD, Jenna M. HulkeiD, Grace A. Arya, Tim AndersoniD
Abstract
The nuclear, internal transcribed spacer (ITS) and mitochondrial cox1 markers are widely used to differentiate Schistosoma haematobium from its livestock counterparts, S. bovis and S. curassoni. Schistosoma isolated from humans with ITS and cox1 alleles from livestock parasites are typically inferred to be zoonotic infections and those with heterozygous ITS alleles (suggesting mixed species ancestry) are classified as recent hybrids. These classifications assume that the ITS and cox1 markers accurately reflect genome-wide ancestry. Here, we evaluated the reliability of this classification scheme by genotyping ITS and cox1 from 132 parasites isolated from human urine, and from 37 adult schistosomes collected from cattle at 14 Nigerian locations. We also genome sequenced each sample to empirically determine livestock schistosome ancestry. ITS/cox1 genotyping suggested extensive recent hybridization and zoonotic infection. Among parasites from humans, 10.1% carried both S. curassoni and S. haematobium ITS, consistent with F1 or early generation hybrids, 21% had livestock schistosome markers at both cox1 and ITS suggesting zoonotic infection, while 13.7% carried S. bovis cox1 alongside mixed S. curassoni and S. haematobium ITS, suggesting more complex ancestry. Genome sequencing revealed a very different picture. All parasites from humans formed a tight cluster regardless of ITS or cox1 genotype, while all worms from cattle were well differentiated. We found no schistosomes containing 50% livestock parasite ancestry consistent with F1s. Instead, we observed regionally varying levels of S. bovis introgression, with modest levels in southern Nigeria (mean = 4.9%) and low levels in northern Nigeria (mean = 0.06%). These results demonstrate that: (i) two-locus genotyping is uninformative for detecting zoonotic infection or recent hybridization between S. haematobium and livestock schistosomes and (ii) previous data generated using this approach requires reinterpretation. These findings reveal the limitations of widely-used approaches for documenting zoonotic infection and hybridization between S. haematobium and livestock schistosome species. Author SummaryMolecular markers provide useful tools for investigating epidemiology of closely related pathogens, but care is required when interpreting the data generated. This study critically examines the interpretation of two widely used molecular markers - the maternally-inherited mitochondrial cox1 locus and ribosomal internal transcribed spacer (ITS) - that are widely used for molecular epidemiology studies of schistosome parasites of humans and livestock. Data from these two markers have been used to infer recent hybridization between livestock and human parasites, and to identify livestock-to-human transmission. To examine the accuracy of epidemiological inferences from the cox1/ITS genotypes we collected parasites from both cattle and people in Nigeria. We genotyped each parasite using cox1 and ITS, as well as using whole-genome sequencing. Cox1/ITS genotyping suggested high rates of hybridization and livestock-to-human transmission. By contrast, whole-genome data show that parasites sampled from humans and cattle are distinct and correspond to independent species. The discrepant results demonstrate that genotyping just two markers provide insufficient resolution and can lead to erroneous conclusions about parasite epidemiology and poor management decisions. We suggest that future work should sample multiple markers from each parasite, using amplicon sequencing or other cost-effective approaches, to improve our understanding of schistosome epidemiology.
- organism: human
- organism: mouse
- assay: whole-genome sequencing
- assay: DNA methylation
- research focus: benchmarking
- research focus: phylogenetics / evolution
- research focus: chromatin & epigenetic regulation
- approach: benchmark study
- approach: statistical model
In-depth summary
Claude Sonnet 4.6 · xhigh effort · 8 Nov 2026 · from the full text of version 2 (17 Aug 2026) · CC BY 4.0
Whole-genome SNVs from Nigerian schistosomes show cox1/ITS genotyping wrongly flags human-derived parasites as recent hybrids or zoonotic infections; they are S. haematobium with low livestock introgression.
Two-locus cox1/ITS genotyping is widely used to infer frequent hybridization and livestock-to-human transmission, which feeds into One Health control proposals, yet its classification schemes were never validated with genomic data. This work tests them against whole-genome data from the same parasites.
Key results
ShownTwo-marker cox1/ITS classification labelled most human-derived miracidia as non-pure S. haematobium, including early-generation, late-generation, complex and zoonotic-like profiles.
90/132 (68.2%) discordant or heterozygous; 42/132 (31.8%) pure S. haematobium. Of the 90: 14 (10.6%) early-generation, 25 (18.9%) late-generation, 19 (14.4%) complex, 28 (21.2%) homozygous S. curassoni ITS with S. bovis mitochondria. All 37 cattle worms were S. bovis. · Categorization of Schistosoma samples using…
ShownGenome-wide PCA and a coalescent phylogeny place all human-derived miracidia in one cluster, separate from livestock schistosomes, whatever their ITS or cox1 class.
PCA of 278,855 unlinked SNVs: PC1 62.2%, PC2 22.6%, one tight human cluster. SVDquartets on 1,061,546 SNVs gave monophyletic species clades, 12.74% of quartets incompatible, 100% support at species nodes. · Principal component analysis (PCA); Nuclear Phylogenetic analysis
ShownSupervised Admixture found only low livestock ancestry in human miracidia, including those that cox1/ITS called F1-like or fully livestock-derived. No recent hybrid or zoonotic infection was detected.
Livestock ancestry in humans 0-18.5%. 14 putative early hybrids: mean 0.04% (0.002%-0.7%). 28 with S. curassoni ITS plus S. bovis cox1: mean 7.17% (0.002%-18.4%). Lab F1 control: 49% S. bovis, 51% S. haematobium. · Admixture; ITS and mtDNA do not identify F1 hybrids or…
ShownSampling location predicted nuclear S. haematobium ancestry almost entirely. ITS genotype added little and mitotype added nothing, so the two markers do not report genome-wide introgression.
Geography-only OLS: R² = 0.927, p = 4.6 × 10⁻⁶⁶. Adding ITS and mitotype: ΔR² = 0.004, model R² = 0.931. ITS partial R² = 0.051 (p = 0.048); mitotype partial R² <0.001 (p = 0.999). · Do ITS genotype and mitotype predict…
ShownLivestock nuclear introgression was largely confined to southern Nigerian sites, with Osun an exception where none was detected. Introgressed S. bovis mitochondria were also far more common in the south. Levels were similar among parasites within a location.
ANOVA across locations F8,123 = 194.97, p = 4.6 x 10^-66. Ebonyi 13.2%-18.4% (mean 16.0%); none detected in Borno, Kebbi, Kano, Katsina, Osun. S. bovis mitotype: 68 of 73 southern vs 5 of 59 northern parasites. · Biogeography of introgression in miracidia from…
Methods, data and reproducibility
Single miracidia and adult worms were genotyped by Sanger ITS and assembled-mitogenome cox1, then compared with 1,877,149 whole-genome SNVs analysed by PCA, supervised Admixture, SVDquartets and OLS regression.
- Miracidia hatched from eggs in urine of infected people, whole-genome sequenced and ITS Sanger genotyped · Schistosoma haematobium (human-derived, Nigeria) · n = 132 miracidia from 53 patients, 9 states · PRJNA635500; PRJNA561522
- Adult schistosomes from cattle at abattoirs, whole-genome sequenced and ITS Sanger genotyped · Schistosoma bovis (cattle, Nigeria) · n = 37 worms from 12 cattle, 4 abattoirs · PRJNA635500; PRJNA561522
- Previously published schistosome genomes: reference species, natural S. bovis x S. curassoni hybrids, lab F1 S. haematobium x S. bovis hybrid, S. margrebowiei outgroup · Schistosoma spp. · n = 34 sequences
- Reference panel for supervised Admixture (three samples per species, outside Nigeria) · S. haematobium, S. bovis, S. curassoni · n = nine samples · SRR11907383, SRR11907395, SRR11907458, ERR3012902, ERR5908628, ERR5919562, SRR11907426, SRR11907431, SRR11907524
- ITS1-5.8S-ITS2 Sanger genotypes at five diagnostic sites; reference chromatogram profiles for peak calling · Schistosoma spp. · n = 169 samples · GenBank PP546795, PP546579, PP546674 (reference profiles)
- Filtered autosomal SNV call set mapped to the S. haematobium reference genome · Schistosoma spp. · n = 1,877,149 variants in 203 individuals · GCF_000699445.3 (reference genome)
- Methods & tools
- Trimmomatic v0.39, BBMap v38.18, GATK v4.5.0.0 HaplotypeCaller, VCFtools v0.1.16, PLINK v1.90b6.2, Admixture v1.3.0 (supervised), SVDquartets v1.49 in PAUP*, GetOrganelle v1.7.7.1, RAxML-NG v1.2.2, Tracy (Indigo algorithm), GenMap v1.3.0, Ordinary least squares regression
- Software & environment
- Conda v24.11.3 environments; environment YAML files and Jupyter notebooks archived in the GitHub repository and on Zenodo.
- Assay conditions & compute
- Standard desktop with at least 4 processors and 16 GB RAM, or a single node with 96 cores and 1 TB RAM.
- Reproducibility gaps
- Per-sample metadata and accessions are given in Supplemental Table 2 rather than in the main text; wall-clock runtimes are not reported.
Limitations (4)
- Geographic scope is Nigeria only; applicability to other West African countries or other parts of Africa is undetermined. · stated by the authors
- Whole-genome sequencing covered 132 human-derived miracidia; deeper sampling might still reveal true hybridization or zoonotic transmission events. · stated by the authors
- The cause of the ITS/mtDNA versus genome discrepancy (lineage sorting versus ancestral gene flow) is discussed but not tested with the data. · noted in this summary
- The supervised Admixture panel has only three samples per species from outside Nigeria, so small livestock-ancestry estimates (for example below 1%) may be sensitive to panel choice. · noted in this summary
Read this if you use cox1/ITS genotyping to infer schistosome hybridization or zoonotic transmission and want genome-scale evidence on how reliable those inferences are.
- Citations
- 1
- FWCI
- 6.14
- Field percentile
- top 3%
Reproducibility ledger
Nobody has reported trying to reproduce this paper yet. Did you replicate an experiment or run the code? Your report saves the next lab a week.
Sign in to file a reproducibility report.
Related papers
- Extant interspecific hybridization among trematodes within the Schistosoma haematobium species complex in Nigeria
O. G. Ajakaye, E. Enabulele, … M. E. Grigg·bioRxiv·Jul 24, 2023
- Morphological and Genomic Characterisation of the Hybrid Schistosome Infecting Humans In Europe Reveals a Complex Admixture Between Schistosoma haematobium and Schistosoma bovis Parasites.
Julien Kincaid-Smith, Alan Tracey, … Ève Toulza·bioRxiv·Jan 3, 2019
- MALDI-TOF : A new tool for the identification of Schistosoma cercariae and detection of hybrids
A. Huguenin, J. Kincaid-Smith, … H. Ferte·bioRxiv·Jun 14, 2022
- Genomic data reveal a north-south split and introgression history of blood fluke (Schistosoma haematobium) populations from across Africa
Roy N. Platt, Egie Elisha Enabulele, … Timothy J. C. Anderson·Nature CommunicationsPeer-reviewed·Apr 13, 2025
- Ancient hybridization and adaptive introgression of an invadolysin gene in schistosome parasites
R. N. Platt, M. McDew-White, … T. J. C. Anderson·bioRxiv·Apr 22, 2019
- Transcriptomic plasticity in hybrid schistosomes can contribute to their zoonotic potential
Nelia Luviano, Eglantine Mathieu‐Bégné, … Eve Toulza·Philosophical Transactions of the Royal Society B Biological SciencesPeer-reviewed·Jan 8, 2026