henhaobio
Sign in

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.

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

  1. 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…

  2. 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

  3. 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…

  4. 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…

  5. 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.

Similarity Map
SelectedSame fieldCross-field
THISAjakaye '23Kincaid-Smi '19Huguenin '22Platt '25Platt '19Luviano '26
80% semantic matchMicrobiology
O. G. Ajakaye, E. Enabulele, J. Balogun+2 more·bioRxiv·Jul 24, 2023