Whole genome amplification and exome sequencing of archived schistosome miracidia

Genomic data reveal a north-south split and introgression history of blood fluke (Schistosoma haematobium) populations from across Africa

The human parasitic fluke, Schistosoma haematobium hybridizes with the livestock parasite S. bovis in the laboratory, but the extent of hybridization in nature is unclear. We analyzed 34.6 million single nucleotide variants in 162 samples from 18 African countries, revealing a sharp genetic discontinuity between northern and southern S. haematobium. We found no evidence for recent hybridization. Instead the data reveal admixture events that occurred 257-879 generations ago in northern S. haematobium populations. Fifteen introgressed S. bovis genes are approaching fixation in northern S. haematobium with four genes potentially driving adaptation. We identified 19 regions that were resistant to introgression; these were enriched on the sex chromosomes. These results (i) demonstrate strong barriers to gene flow between these species, (ii) indicate that hybridization may be less common than currently envisaged, but (iii) reveal profound genomic consequences of interspecific hybridization between schistosomes of medical and veterinary importance.

Hybrid-Capture Target Enrichment in Human Pathogens: Identification, Evolution, Biosurveillance, and Genomic Epidemiology

High-throughput sequencing (HTS) has revolutionised the field of pathogen genomics, enabling the direct recovery of pathogen genomes from clinical and environmental samples. However, pathogen nucleic acids are often overwhelmed by those of the host, requiring deep metagenomic sequencing to recover sufficient sequences for downstream analyses (e.g., identification and genome characterisation). To circumvent this, hybrid-capture target enrichment (HC) is able to enrich pathogen nucleic acids across multiple scales of divergences and taxa, depending on the panel used. In this review, we outline the applications of HC in human pathogens-bacteria, fungi, parasites and viruses-including identification, genomic epidemiology, antimicrobial resistance genotyping, and evolution. Importantly, we explored the applicability of HC to clinical metagenomics, which ultimately requires more work before it is a reliable and accurate tool for clinical diagnosis. Relatedly, the utility of HC was exemplified by COVID-19, which was used as a case study to illustrate the maturity of HC for recovering pathogen sequences. As we unravel the origins of COVID-19, zoonoses remain more relevant than ever. Therefore, the role of HC in biosurveillance studies is also highlighted in this review, which is critical in preparing us for the next pandemic. We also found that while HC is a popular tool to study viruses, it remains underutilised in parasites and fungi and, to a lesser extent, bacteria. Finally, weevaluated the future of HC with respect to bait design in the eukaryotic groups and the prospect of combining HC with long-read HTS.

Population genomics of helminth parasites

Next generation sequencing technologies have facilitated a shift from a few targeted loci in population genetic studies to whole genome approaches. Here, we review the types of questions and inferences regarding the population biology and evolution of parasitic helminths being addressed within the field of population genomics. Topics include parabiome, hybridization, population structure, loci under selection and linkage mapping. We highlight various advances, and note the current trends in the field, particularly a focus on human-related parasites despite the inherent biodiversity of helminth species. We conclude by advocating for a broader application of population genomics to reflect the taxonomic and life history breadth displayed by helminth parasites. As such, our basic knowledge about helminth population biology and evolution would be enhanced while the diversity of helminths in itself would facilitate population genomic comparative studies to address broader ecological and evolutionary concepts.

A review of the genetic determinants of praziquantel resistance in Schistosoma mansoni: Is praziquantel and intestinal schistosomiasis a perfect match?

Schistosomiasis is a neglected tropical disease (NTD) caused by parasitic trematodes belonging to the Schistosoma genus. The mainstay of schistosomiasis control is the delivery of a single dose of praziquantel (PZQ) through mass drug administration (MDA) programs. These programs have been successful in reducing the prevalence and intensity of infections. Due to the success of MDA programs, the disease has recently been targeted for elimination as a public health problem in some endemic settings. The new World Health Organization (WHO) treatment guidelines aim to provide equitable access to PZQ for individuals above two years old in targeted areas. The scale up of MDA programs may heighten the drug selection pressures on Schistosoma parasites, which could lead to the emergence of PZQ resistant schistosomes. The reliance on a single drug to treat a disease of this magnitude is worrying should drug resistance develop. Therefore, there is a need to detect and track resistant schistosomes to counteract the threat of drug resistance to the WHO 2030 NTD roadmap targets. Until recently, drug resistance studies have been hindered by the lack of molecular markers associated with PZQ resistance. This review discusses recent significant advances in understanding the molecular basis of PZQ action in S. mansoni and proposes additional genetic determinants associated with PZQ resistance. PZQ resistance will also be analyzed in the context of alternative factors that may decrease efficacy within endemic field settings, and the most recent treatment guidelines recommended by the WHO.

Genome-wide analysis of Schistosoma mansoni reveals limited population structure and possible praziquantel drug selection pressure within Ugandan hot-spot communities

Populations within schistosomiasis control areas, especially those in Africa, are recommended to receive regular mass drug administration (MDA) with praziquantel (PZQ) as the main strategy for controlling the disease. The impact of PZQ treatment on schistosome genetics remains poorly understood, and is limited by a lack of high-resolution genetic data on the population structure of parasites within these control areas. We generated whole-genome sequence data from 174 individual miracidia collected from both children and adults from fishing communities on islands in Lake Victoria in Uganda that had received either annual or quarterly MDA with PZQ over four years, including samples collected immediately before and four weeks after treatment. Genome variation within and between samples was characterised and we investigated genomic signatures of natural selection acting on these populations that could be due to PZQ treatment. The parasite population on these islands was more diverse than found in nearby villages on the lake shore. We saw little or no genetic differentiation between villages, or between the groups of villages with different treatment intensity, but slightly higher genetic diversity within the pre-treatment compared to post-treatment parasite populations. We identified classes of genes significantly enriched within regions of the genome with evidence of recent positive selection among post-treatment and intensively treated parasite populations. The differential selection observed in post-treatment and pre-treatment parasite populations could be linked to any reduced susceptibility of parasites to praziquantel treatment.

Genomic analysis of a parasite invasion: Colonization of the Americas by the blood fluke Schistosoma mansoni

Schistosoma mansoni, a snail-borne, blood fluke that infects humans, was introduced into the Americas from Africa during the Trans-Atlantic slave trade. As this parasite shows strong specificity to the snail intermediate host, we expected that adaptation to South American Biomphalaria spp. snails would result in population bottlenecks and strong signatures of selection. We scored 475,081 single nucleotide variants in 143 S. mansoni from the Americas (Brazil, Guadeloupe and Puerto Rico) and Africa (Cameroon, Niger, Senegal, Tanzania, and Uganda), and used these data to ask: (i) Was there a population bottleneck during colonization? (ii) Can we identify signatures of selection associated with colonization? (iii) What were the source populations for colonizing parasites? We found a 2.4- to 2.9-fold reduction in diversity and much slower decay in linkage disequilibrium (LD) in parasites from East to West Africa. However, we observed similar nuclear diversity and LD in West Africa and Brazil, suggesting no strong bottlenecks and limited barriers to colonization. We identified five genome regions showing selection in the Americas, compared with three in West Africa and none in East Africa, which we speculate may reflect adaptation during colonization. Finally, we infer that unsampled populations from central African regions between Benin and Angola, with contributions from Niger, are probably the major source(s) for Brazilian S. mansoni. The absence of a bottleneck suggests that this is a rare case of a serendipitous invasion, where S. mansoni parasites were pre-adapted to the Americas and able to establish with relative ease.

Integrating genomic and epidemiologic data to accelerate progress toward schistosomiasis elimination

The global community has adopted ambitious goals to eliminate schistosomiasis as a public health problem, and new tools are needed to achieve them. Mass drug administration programs, for example, have reduced the burden of schistosomiasis, but the identification of hotspots of persistent and reemergent transmission threaten progress toward elimination and underscore the need to couple treatment with interventions that reduce transmission. Recent advances in DNA sequencing technologies make whole-genome sequencing a valuable and increasingly feasible option for population-based studies of complex parasites such as schistosomes. Here, we focus on leveraging genomic data to tailor interventions to distinct social and ecological circumstances. We consider two priority questions that can be addressed by integrating epidemiological, ecological, and genomic information: (1) how often do non-human host species contribute to human schistosome infection? and (2) what is the importance of locally acquired versus imported infections in driving transmission at different stages of elimination? These questions address processes that can undermine control programs, especially those that rely heavily on treatment with praziquantel. Until recently, these questions were difficult to answer with sufficient precision to inform public health decision-making. We review the literature related to these questions and discuss how whole-genome approaches can identify the geographic and taxonomic sources of infection, and how such information can inform context-specific efforts that advance schistosomiasis control efforts and minimize the risk of reemergence.

Genetic analysis of praziquantel response in schistosome parasites implicates a transient receptor potential channel

[Figure: see text].

Genetic architecture of transmission stage production and virulence in schistosome parasites

Both theory and experimental data from pathogens suggest that the production of transmission stages should be strongly associated with virulence, but the genetic bases of parasite transmission/virulence traits are poorly understood. The blood fluke Schistosoma mansoni shows extensive variation in numbers of cercariae larvae shed and in their virulence to infected snail hosts, consistent with expected trade-offs between parasite transmission and virulence. We crossed schistosomes from two populations that differ 8-fold in cercarial shedding and in their virulence to Biomphalaria glabrata snail hosts, and determined four-week cercarial shedding profiles in F0 parents, F1 parents and 376 F2 progeny from two independent crosses in inbred snails. Sequencing and linkage analysis revealed that cercarial production is polygenic and controlled by five QTLs (i.e. Quantitative Trait Loci). These QTLs act additively, explaining 28.56% of the phenotypic variation. These results demonstrate that the genetic architecture of key traits relevant to schistosome ecology can be dissected using classical linkage mapping approaches.

Patterns of relatedness and genetic diversity inferred from whole genome sequencing of archival blood fluke miracidia (Schistosoma japonicum)

Genomic approaches hold great promise for resolving unanswered questions about transmission patterns and responses to control efforts for schistosomiasis and other neglected tropical diseases. However, the cost of generating genomic data and the challenges associated with obtaining sufficient DNA from individual schistosome larvae (miracidia) from mammalian hosts have limited the application of genomic data for studying schistosomes and other complex macroparasites. Here, we demonstrate the feasibility of utilizing whole genome amplification and sequencing (WGS) to analyze individual archival miracidia. As an example, we sequenced whole genomes of 22 miracidia from 11 human hosts representing two villages in rural Sichuan, China, and used these data to evaluate patterns of relatedness and genetic diversity. We also down-sampled our dataset to test how lower coverage sequencing could increase the cost effectiveness of WGS while maintaining power to accurately infer relatedness. Collectively, our results illustrate that population-level WGS datasets are attainable for individual miracidia and represent a powerful tool for ultimately providing insight into overall genetic diversity, parasite relatedness, and transmission patterns for better design and evaluation of disease control efforts.

Parasite Population Genetic Contributions to the Schistosomiasis Consortium for Operational Research and Evaluation within Sub-Saharan Africa

Analyses of the population genetic structure of schistosomes under the "Schistosomiasis Consortium for Operational Research and Evaluation" (SCORE) contrasting treatment pressure scenarios in Tanzania, Niger, and Zanzibar were performed to provide supplementary critical information with which to evaluate the impact of these large-scale control activities and guide how activities could be adjusted. We predicted that population genetic analyses would reveal information on a range of important parameters including, but not exclusive to, recruitment and transmission of genotypes, occurrence of hybridization events, differences in reproductive mode, and degrees of inbreeding, and hence, the evolutionary potential, and responses of parasite populations under contrasting treatment pressures. Key findings revealed that naturally high levels of gene flow and mixing of the parasite populations between neighboring sites were likely to dilute any effects imposed by the SCORE treatment arms. Furthermore, significant inherent differences in parasite fecundity were observed, independent of current treatment arm, but potentially of major impact in terms of maintaining high levels of ongoing transmission in persistent "biological hotspot" sites. Within Niger, naturally occurring Schistosoma haematobium/Schistosoma bovis viable hybrids were found to be abundant, often occurring in significantly higher proportions than that of single-species S. haematobium infections. By examining parasite population genetic structures across hosts, treatment regimens, and the spatial landscape, our results to date illustrate key transmission processes over and above that which could be achieved through standard parasitological monitoring of prevalence and intensity alone, as well as adding to our understanding of Schistosoma spp. life history strategies in general.

Snail-Related Contributions from the Schistosomiasis Consortium for Operational Research and Evaluation Program Including Xenomonitoring, Focal Mollusciciding, Biological Control, and Modeling

The Schistosomiasis Consortium for Operational Research and Evaluation (SCORE) was created in 2008 to answer questions of importance to program managers working to reduce the burden of schistosomiasis in Africa. In the past, intermediate host snail monitoring and control was an important part of integrated schistosomiasis control. However, in Africa, efforts to control snails have declined dramatically over the last 30 years. A resurgence of interest in the control of snails has been prompted by the realization, backed by a World Health Assembly resolution (WHA65.21), that mass drug administration alone may be insufficient to achieve schistosomiasis elimination. SCORE has supported work on snail identification and mapping and investigated how xenomonitoring techniques can aid in the identification of infected snails and thereby identify potential transmission areas. Focal mollusciciding with niclosamide was undertaken in Zanzibar and Côte d'Ivoire as a part of elimination studies. Two studies involving biological control of snails were conducted: one explored the association of freshwater riverine prawns and snail hosts in Côte d'Ivoire and the other assessed the current distribution of Procambarus clarkii, the invasive Louisiana red swamp crayfish, in Kenya and its association with snail hosts and schistosomiasis transmission. SCORE also supported modeling studies on the importance of snail control in achieving elimination and a meta-analysis of the impact of molluscicide-based snail control programs on human schistosomiasis prevalence and incidence. SCORE's snail control studies contributed to increased investment in building capacity, and specimens collected during SCORE research deposited in the Schistosomiasis Collections at the Natural History Museum (SCAN) will provide a valuable resource for the years to come.

Oxamniquine resistance alleles are widespread in Old World Schistosoma mansoni and predate drug deployment

Do mutations required for adaptation occur de novo, or are they segregating within populations as standing genetic variation? This question is key to understanding adaptive change in nature, and has important practical consequences for the evolution of drug resistance. We provide evidence that alleles conferring resistance to oxamniquine (OXA), an antischistosomal drug, are widespread in natural parasite populations under minimal drug pressure and predate OXA deployment. OXA has been used since the 1970s to treat Schistosoma mansoni infections in the New World where S. mansoni established during the slave trade. Recessive loss-of-function mutations within a parasite sulfotransferase (SmSULT-OR) underlie resistance, and several verified resistance mutations, including a deletion (p.E142del), have been identified in the New World. Here we investigate sequence variation in SmSULT-OR in S. mansoni from the Old World, where OXA has seen minimal usage. We sequenced exomes of 204 S. mansoni parasites from West Africa, East Africa and the Middle East, and scored variants in SmSULT-OR and flanking regions. We identified 39 non-synonymous SNPs, 4 deletions, 1 duplication and 1 premature stop codon in the SmSULT-OR coding sequence, including one confirmed resistance deletion (p.E142del). We expressed recombinant proteins and used an in vitro OXA activation assay to functionally validate the OXA-resistance phenotype for four predicted OXA-resistance mutations. Three aspects of the data are of particular interest: (i) segregating OXA-resistance alleles are widespread in Old World populations (4.29–14.91% frequency), despite minimal OXA usage, (ii) two OXA-resistance mutations (p.W120R, p.N171IfsX28) are particularly common (>5%) in East African and Middle-Eastern populations, (iii) the p.E142del allele has identical flanking SNPs in both West Africa and Puerto Rico, suggesting that parasites bearing this allele colonized the New World during the slave trade and therefore predate OXA deployment. We conclude that standing variation for OXA resistance is widespread in S. mansoni.

It has been argued that drug resistance is unlikely to spread rapidly in helminth parasites infecting humans. This is based, at least in part, on the premise that resistance mutations are rare or absent within populations prior to treatment, and take a long time to reach appreciable frequencies because helminth parasite generation time is long. This argument is critically dependent on the starting frequency of resistance alleles–if high levels of “standing variation” for resistance are present prior to deployment of treatment, resistance may spread rapidly. We examined frequencies of oxamniquine resistance alleles present in Schistosoma mansoni from Africa and the Middle East where oxamniquine has seen minimal use. We found that oxamniquine resistance alleles are widespread in the Old World, ranging from 4.29% in the Middle East to 14.91% in East African parasite populations. Furthermore, we show that resistance alleles from West African and the Caribbean schistosomes share a common origin, suggesting that these alleles travelled to the New World with S. mansoni during the transatlantic slave trade. Together, these results demonstrate extensive standing variation for oxamniquine resistance. Our results have important implications for both drug treatment policies and drug development efforts, and demonstrate the power of molecular surveillance approaches for guiding helminth control.

Ancient Hybridization and Adaptive Introgression of an Invadolysin Gene in Schistosome Parasites

Introgression among parasite species has the potential to transfer traits of biomedical importance across species boundaries. The parasitic blood fluke Schistosoma haematobium causes urogenital schistosomiasis in humans across sub-Saharan Africa. Hybridization with other schistosome species is assumed to occur commonly, because genetic crosses between S. haematobium and livestock schistosomes, including S. bovis, can be staged in the laboratory, and sequencing of mtDNA and rDNA amplified from microscopic miracidia larvae frequently reveals markers from different species. However, the frequency, direction, age, and genomic consequences of hybridization are unknown. We hatched miracidia from eggs and sequenced the exomes from 96 individual S. haematobium miracidia from infected patients from Niger and the Zanzibar archipelago. These data revealed no evidence for contemporary hybridization between S. bovis and S. haematobium in our samples. However, all Nigerien S. haematobium genomes sampled show hybrid ancestry, with 3.3-8.2% of their nuclear genomes derived from S. bovis, providing evidence of an ancient introgression event that occurred at least 108-613 generations ago. Some S. bovis-derived alleles have spread to high frequency or reached fixation and show strong signatures of directional selection; the strongest signal spans a single gene in the invadolysin gene family (Chr. 4). Our results suggest that S. bovis/S. haematobium hybridization occurs rarely but demonstrate profound consequences of ancient introgression from a livestock parasite into the genome of S. haematobium, the most prevalent schistosome species infecting humans.

Evaluation of DNA Extraction Methods on Individual Helminth Egg and Larval Stages for Whole-Genome Sequencing

Whole-genome sequencing is being rapidly applied to the study of helminth genomes, including de novo genome assembly, population genetics, and diagnostic applications. Although late-stage juvenile and adult parasites typically produce sufficient DNA for molecular analyses, these parasitic stages are almost always inaccessible in the live host; immature life stages found in the environment for which samples can be collected non-invasively offer a potential alternative; however, these samples typically yield very low quantities of DNA, can be environmentally resistant, and are susceptible to contamination, often from bacterial or host DNA. Here, we have tested five low-input DNA extraction protocols together with a low-input sequencing library protocol to assess the feasibility of whole-genome sequencing of individual immature helminth samples. These approaches do not use whole-genome amplification, a common but costly approach to increase the yield of low-input samples. We first tested individual parasites from two species spotted onto FTA cards-egg and L1 stages of Haemonchus contortus and miracidia of Schistosoma mansoni-before further testing on an additional five species-Ancylostoma caninum, Ascaridia dissimilis, Dirofilaria immitis, Strongyloides stercoralis, and Trichuris muris-with an optimal protocol. A sixth species-Dracunculus medinensis-was included for comparison. Whole-genome sequencing followed by analyses to determine the proportion of on- and off-target mapping revealed successful sample preparations for six of the eight species tested with variation both between species and between different life stages from some species described. These results demonstrate the feasibility of whole-genome sequencing of individual parasites, and highlight a new avenue toward generating sensitive, specific, and information-rich data for the diagnosis and surveillance of helminths.

Whole-genome sequence of the bovine blood fluke Schistosoma bovis supports interspecific hybridization with S. haematobium

Mesenteric infection by the parasitic blood fluke Schistosoma bovis is a common veterinary problem in Africa and the Middle East and occasionally in the Mediterranean Region. The species also has the ability to form interspecific hybrids with the human parasite S. haematobium with natural hybridisation observed in West Africa, presenting possible zoonotic transmission. Additionally, this exchange of alleles between species may dramatically influence disease dynamics and parasite evolution. We have generated a 374 Mb assembly of the S. bovis genome using Illumina and PacBio-based technologies. Despite infecting different hosts and organs, the genome sequences of S. bovis and S. haematobium appeared strikingly similar with 97% sequence identity. The two species share 98% of protein-coding genes, with an average sequence identity of 97.3% at the amino acid level. Genome comparison identified large continuous parts of the genome (up to several 100 kb) showing almost 100% sequence identity between S. bovis and S. haematobium. It is unlikely that this is a result of genome conservation and provides further evidence of natural interspecific hybridization between S. bovis and S. haematobium. Our results suggest that foreign DNA obtained by interspecific hybridization was maintained in the population through multiple meiosis cycles and that hybrids were sexually reproductive, producing viable offspring. The S. bovis genome assembly forms a highly valuable resource for studying schistosome evolution and exploring genetic regions that are associated with species-specific phenotypic traits.

In this article we detail the assembly and functional annotation of the Schistosoma bovis genome. S. bovis is a parasitic flatworm that primarily infects bovines, with important economic consequences in affected countries. However, it is also a close relative of the human carcinogenic parasite Schistosoma haematobium which is a serious health issue in many endemic countries in Sub-Saharan Africa. The close relationship and overlapping geographical distribution of S. bovis and S. haematobium allows these to hybridise in the wild increasing their genetic diversity and presenting the risk of zoonotic transmission, i.e. the transmission from animals to humans. The hybridization between human and ruminant schistosomes is of particular interest as interspecific hybridization may have dramatic impacts on transmission rates, disease dynamics, control interventions and parasite evolution. By whole-genome sequencing and comparative genomics we present evidence that fertile hybrids are indeed present in the wild, presenting the potential risk of transmission from animal reservoirs to humans.

Advancing the multi-disciplinarity of parasitology within the British Society for Parasitology: studies of host-parasite evolution in an ever-changing world

The study of parasites typically crosses into other research disciplines and spans across diverse scales, from molecular- to populational-levels, notwithstanding promoting an understanding of parasites set within evolutionary time. Today, the 2030 Sustainable Development Goals (SDGs) help frame much of contemporary parasitological research, since parasites can be found in all ecosystems, blighting human, animal and plant health. In recognition of the multi-disciplinary nature of parasitological research, the 2017 Autumn Symposium of the British Society for Parasitology was held in London to provide a forum for novel exchange across medical, veterinary and wildlife fields of study. Whilst the meeting was devoted to the topic of parasitism, it sought to foster mutualism, the antithesis perhaps of parasitism, by forging new academic connections and social networks to exchange novel ideas. The meeting also celebrated the longstanding career of Professor David Rollinson, FLS in the award of the International Federation for Tropical Medicine Medal for his efforts spanning 40 years of parasitological research. Indeed, David has done so much to explore and promote the fascinating biology of parasitism, as exemplified by the 15 manuscripts contained within this Special Issue.

Genetic Crosses and Linkage Mapping in Schistosome Parasites

Linkage mapping - utilizing experimental genetic crosses to examine cosegregation of phenotypic traits with genetic markers - is now 100 years old. Schistosome parasites are exquisitely well suited to linkage mapping approaches because genetic crosses can be conducted in the laboratory, thousands of progeny are produced, and elegant experimental work over the last 75 years has revealed heritable genetic variation in multiple biomedically important traits such as drug resistance, host specificity, and virulence. Application of this approach is timely because the improved genome assembly for Schistosoma mansoni and developing molecular toolkit for schistosomes increase our ability to link phenotype with genotype. We describe current progress and potential future directions of linkage mapping in schistosomes.