Whole-genome sequencing of Schistosoma mansoni reveals extensive diversity with limited selection despite mass drug administration

Characterizing genetic variation on the Z chromosome in Schistosoma japonicum reveals host-parasite co-evolution

BACKGROUND

Schistosomiasis is a neglected tropical disease that afflicts millions of people worldwide; it is caused by Schistosoma, the only dioecious flukes with ZW systems. Schistosoma japonicum is endemic to Asia; the Z chromosome of S. japonicum comprises one-quarter of the entire genome. Detection of positive selection using resequencing data to understand adaptive evolution has been applied to a variety of pathogens, including S. japonicum. However, the contribution of the Z chromosome to evolution and adaptation is often neglected.

METHODS

We obtained 1,077,526 high-quality SNPs on the Z chromosome in 72 S. japonicum using re-sequencing data publicly. To examine the faster Z effect, we compared the sequence divergence of S. japonicum with two closely related species, Schistosoma haematobium and S. mansoni. Genetic diversity was compared between the Z chromosome and autosomes in S. japonicum by calculating the nucleotide diversity (π) and Dxy values. Population structure was also assessed based on PCA and structure analysis. Besides, we employed multiple methods including Tajima's D, FST, iHS, XP-EHH, and CMS to detect positive selection signals on the Z chromosome. Further RNAi knockdown experiments were performed to investigate the potential biological functions of the candidate genes.

RESULTS

Our study found that the Z chromosome of S. japonicum showed faster evolution and more pronounced genetic divergence than autosomes, although the effect may be smaller than the variation among genes. Compared with autosomes, the Z chromosome in S. japonicum had a more pronounced genetic divergence of sub-populations. Notably, we identified a set of candidate genes associated with host-parasite co-evolution. In particular, LCAT exhibited significant selection signals within the Taiwan population. Further RNA interference experiments suggested that LCAT is necessary for S. japonicum survival and propagation in the definitive host. In addition, we identified several genes related to the specificity of the intermediate host in the C-M population, including Rab6 and VCP, which are involved in adaptive immune evasion to the host.

CONCLUSIONS

Our study provides valuable insights into the adaptive evolution of the Z chromosome in S. japonicum and further advances our understanding of the co-evolution of this medically important parasite and its hosts.

Development of non-sedating antischistosomal benzodiazepines

The neglected tropical disease schistosomiasis infects over 200 million people worldwide and is treated with just one broad spectrum antiparasitic drug (praziquantel). Alternative drugs are needed in the event of emerging praziquantel resistance or treatment failure. One promising lead that has shown efficacy in animal models and a human clinical trial is the benzodiazepine meclonazepam, discovered by Roche in the 1970's. Meclonazepam was not brought to market because of dose-limiting sedative side effects. However, the human target of meclonazepam that causes sedation (GABAARs) are not orthologous to the parasite targets that cause worm death. Therefore, we were interested in whether the structure of meclonazepam could be modified to produce antiparasitic benzodiazepines that do not cause host sedation. We synthesized 18 meclonazepam derivatives with modifications at different positions on the benzodiazepine ring system and tested them for in vitro antiparasitic activity. This identified five compounds that progressed to in vivo screening in a murine model, two of which cured parasite infections with comparable potency to meclonazepam. When these two compounds were administered to mice that were run on the rotarod test, both were less sedating than meclonazepam. These findings demonstrate the proof of concept that meclonazepam analogs can be designed with an improved therapeutic index, and point to the C3 position of the benzodiazepine ring system as a logical site for further structure-activity exploration to further optimize this chemical series.

Therapeutic Potential of Natural Products in the Treatment of Schistosomiasis

It is estimated that 250 million people worldwide are affected by schistosomiasis. Disease transmission is related to the poor sanitation and hygiene habits that affect residents of impoverished regions in tropical and subtropical countries. The main species responsible for causing disease in humans are Schistosoma Mansoni, S. japonicum, and S. haematobium, each with different geographic distributions. Praziquantel is the drug predominantly used to treat this disease, which offers low effectiveness against immature and juvenile parasite forms. In addition, reports of drug resistance prompt the development of novel therapeutic approaches. Natural products represent an important source of new compounds, especially those obtained from plant sources. This review compiles data from several in vitro and in vivo studies evaluating various compounds and essential oils derived from plants with cercaricidal and molluscicidal activities against both juvenile and adult forms of the parasite. Finally, this review provides an important discussion on recent advances in molecular and computational tools deemed fundamental for more rapid and effective screening of new compounds, allowing for the optimization of time and resources.

Praziquantel inhibits Caenorhabditis elegans development and species-wide differences might be cct-8-dependent

Anthelmintic drugs are used to treat parasitic roundworm and flatworm infections in humans and other animals. Caenorhabditis elegans is an established model to investigate anthelmintics used to treat roundworms. In this study, we use C. elegans to examine the mode of action and the mechanisms of resistance against the flatworm anthelmintic drug praziquantel (PZQ), used to treat trematode and cestode infections. We found that PZQ inhibited development and that this developmental delay varies by genetic background. Interestingly, both enantiomers of PZQ are equally effective against C. elegans, but the right-handed PZQ (R-PZQ) is most effective against schistosome infections. We conducted a genome-wide association mapping with 74 wild C. elegans strains to identify a region on chromosome IV that is correlated with differential PZQ susceptibility. Five candidate genes in this region: cct-8, znf-782, Y104H12D.4, Y104H12D.2, and cox-18, might underlie this variation. The gene cct-8, a subunit of the protein folding complex TRiC, has variation that causes a putative protein coding change (G226V), which is correlated with reduced developmental delay. Gene expression analysis suggests that this variant correlates with slightly increased expression of both cct-8 and hsp-70. Acute exposure to PZQ caused increased expression of hsp-70, indicating that altered TRiC function might be involved in PZQ responses. To test if this variant affects development upon exposure to PZQ, we used CRISPR-Cas9 genome editing to introduce the V226 allele into the N2 genetic background (G226) and the G226 allele into the JU775 genetic background (V226). These experiments revealed that this variant was not sufficient to explain the effects of PZQ on development. Nevertheless, this study shows that C. elegans can be used to study PZQ mode of action and resistance mechanisms. Additionally, we show that the TRiC complex requires further evaluation for PZQ responses in C. elegans.

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.

Diagnosis of helminths depends on worm fecundity and the distribution of parasites within hosts

Helminth transmission and morbidity are dependent on the number of mature parasites within a host; however, observing adult worms is impossible for many natural infections. An outstanding challenge is therefore relating routine diagnostics, such as faecal egg counts, to the underlying worm burden. This relationship is complicated by density-dependent fecundity (egg output per worm reduces due to crowding at high burdens) and the skewed distribution of parasites (majority of helminths aggregated in a small fraction of hosts). We address these questions for the carcinogenic liver fluke Opisthorchis viverrini, which infects approximately 10 million people across Southeast Asia, by analysing five epidemiological surveys (n = 641) where adult flukes were recovered. Using a mechanistic model, we show that parasite fecundity varies between populations, with surveys from Thailand and Laos demonstrating distinct patterns of egg output and density-dependence. As the probability of observing faecal eggs increases with the number of mature parasites within a host, we quantify diagnostic sensitivity as a function of the worm burden and find that greater than 50% of cases are misdiagnosed as false negative in communities close to elimination. Finally, we demonstrate that the relationship between observed prevalence from routine diagnostics and true prevalence is nonlinear and strongly influenced by parasite aggregation.

An adaptive phase II/III safety and efficacy randomized controlled trial of single day or three-day fixed-dose albendazole-ivermectin co-formulation versus albendazole for the treatment of Trichuris trichiura and other STH infections. ALIVE trial protocol

Background: Soil-transmitted helminths (STH) are targeted for control through mass drug-administration campaigns to prevent morbidity affecting at-risk groups in endemic regions. Although broadly successful, the use of albendazole and mebendazole achieved variable progress, with deficiencies against Trichuris trichiura and a predictable low efficacy against Strongyloides stercoralis. Novel drug combinations offer a potential solution, providing they can be delivered safely and maintain efficacy against all STH species. Here we present the protocol of a clinical trial to evaluate a fixed-dose combination (FDC) tablet containing albendazole and ivermectin that will be compared against albendazole against STH . Methods: An adaptive phase II/III randomized controlled trial will be undertaken in STH endemic sites in Ethiopia, Kenya and Mozambique to evaluate an oral FDC of 400 mg albendazole and either 9- or 18 mg ivermectin. FDC will be administered as a single dose or single doses over three-consecutive days and assessed against a single dose of 400 mg albendazole. In the phase II trial, 126 T. trichiura-infected children weighting 15 to 45 kg will be treated in a dose-escalation manner to determine safety objectives. In the phase III trial, 1097 participants aged 5 to 18 years old infected with T. trichiura, hookworm and S. stercoralis will be recruited to determine safety and efficacy. The trial will be open-label with blinded outcome assessors. Cure rate measured 21-days after-treatment in duplicate Kato-Katz is the primary efficacy outcome. Secondary objectives include efficacy evaluation by quantitative polymerase chain reaction (PCR) as an outcome measurement, description of pharmacokinetic parameters, palatability and acceptability evaluations, and monitoring of anthelmintic resistance. Conclusions: This trial with registrational goals seeks to evaluate an innovative fixed-dose combination of albendazole and ivermectin co-formulated tablets, with the goal of providing an anthelmintic regimen with improved efficacy and spectrum of coverage against STH. ClinicalTrials.gov registration: NCT05124691 (18/11/2021).

Praziquantel: An update on the mechanism of its action against schistosomiasis and new therapeutic perspectives

Praziquantel (PZQ) is the drug of choice for the treatment of all forms of schistosomiasis, although its mechanisms of action are not completely understood. PZQ acts largely on adult worms. This narrative literature review describes what is known about the mechanisms of action of PZQ against schistosomes from in vitro and in vivo studies and highlights the molecular targets in parasites and immune responses induced in definitive hosts by this drug. Moreover, new therapeutic uses of PZQ are discussed. Studies have demonstrated that in addition to impacting voltage-operated Ca2 + channels, PZQ may interact with other schistosome molecules, such as myosin regulatory light chain, glutathione S-transferase, and transient receptor potential channels. Following PZQ administration, increased T regulatory type 1 (Tr1) cell differentiation and decreased inflammation were observed, indicating that PZQ promotes immunoregulatory pathways. Although PZQ is widely used in mass drug administration schemes, the existence of resistant parasites has not been proven; however, it is a concern that should be constantly investigated in human populations. In addition, we discuss studies that evaluate health applications of PZQ (other than helminth infection), such as its effect in cancer therapy and its adjuvant action in vaccines against viruses.

Genomic landscape of drug response reveals mediators of anthelmintic resistance

Like other pathogens, parasitic helminths can rapidly evolve resistance to drug treatment. Understanding the genetic basis of anthelmintic drug resistance in parasitic nematodes is key to tracking its spread and improving the efficacy and sustainability of parasite control. Here, we use an in vivo genetic cross between drug-susceptible and multi-drug-resistant strains of Haemonchus contortus in a natural host-parasite system to simultaneously map resistance loci for the three major classes of anthelmintics. This approach identifies new alleles for resistance to benzimidazoles and levamisole and implicates the transcription factor cky-1 in ivermectin resistance. This gene is within a locus under selection in ivermectin-resistant populations worldwide; expression analyses and functional validation using knockdown experiments support that cky-1 is associated with ivermectin survival. Our work demonstrates the feasibility of high-resolution forward genetics in a parasitic nematode and identifies variants for the development of molecular diagnostics to combat drug resistance in the field.

Potential drivers for schistosomiasis persistence: Population genetic analyses from a cluster-randomized urogenital schistosomiasis elimination trial across the Zanzibar islands

The World Health Organization's revised NTD Roadmap and the newly launched Guidelines target elimination of schistosomiasis as a public health problem in all endemic areas by 2030. Key to meeting this goal is elucidating how selective pressures imposed by interventions shape parasite populations. Our aim was to identify any differential impact of a unique cluster-randomized tri-armed elimination intervention (biannual mass drug administration (MDA) applied alone or in association with either mollusciciding (snail control) or behavioural change interventions) across two Zanzibarian islands (Pemba and Unguja) on the population genetic composition of Schistosoma haematobium over space and time. Fifteen microsatellite loci were used to analyse individual miracidia collected from infected individuals across islands and intervention arms at the start (2012 baseline: 1,522 miracidia from 176 children; 303 from 43 adults; age-range 6-75, mean 12.7 years) and at year 5 (2016: 1,486 miracidia from 146 children; 214 from 25 adults; age-range 9-46, mean 12.4 years). Measures of genetic diversity included allelic richness (Ar), Expected (He) and Observed heterozygosity (Ho), inbreeding coefficient (FST), parentage analysis, estimated worm burden, worm fecundity, and genetic sub-structuring. There was little evidence of differential selective pressures on population genetic diversity, inbreeding or estimated worm burdens by treatment arm, with only the MDA+snail control arm within Unguja showing trends towards reduced diversity and altered inbreeding over time. The greatest differences overall, both in terms of parasite fecundity and genetic sub-structuring, were observed between the islands, consistent with Pemba's persistently higher mean infection intensities compared to neighbouring Unguja, and within islands in terms of infection hotspots (across three definitions). These findings highlight the important contribution of population genetic analyses to elucidate extensive genetic diversity and biological drivers, including potential gene-environmental factors, that may override short term selective pressures imposed by differential disease control strategies. Trial Registration: ClinicalTrials.gov ISRCTN48837681.

Prospects for genomic surveillance for selection in schistosome parasites

Schistosomiasis is a neglected tropical disease caused by multiple parasitic Schistosoma species, and which impacts over 200 million people globally, mainly in low- and middle-income countries. Genomic surveillance to detect evidence for natural selection in schistosome populations represents an emerging and promising approach to identify and interpret schistosome responses to ongoing control efforts or other environmental factors. Here we review how genomic variation is used to detect selection, how these approaches have been applied to schistosomes, and how future studies to detect selection may be improved. We discuss the theory of genomic analyses to detect selection, identify experimental designs for such analyses, and review studies that have applied these approaches to schistosomes. We then consider the biological characteristics of schistosomes that are expected to respond to selection, particularly those that may be impacted by control programs. Examples include drug resistance, host specificity, and life history traits, and we review our current understanding of specific genes that underlie them in schistosomes. We also discuss how inherent features of schistosome reproduction and demography pose substantial challenges for effective identification of these traits and their genomic bases. We conclude by discussing how genomic surveillance for selection should be designed to improve understanding of schistosome biology, and how the parasite changes in response to selection.

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.

Improving translational power in antischistosomal drug discovery

Schistosomiasis is a poverty-associated tropical disease caused by blood dwelling trematodes that threaten approximately 10% of the world population. Praziquantel, the sole drug currently available for treatment, is insufficient to eliminate the disease and the clinical drug development pipeline is empty. Here, we review the characteristics of the patent Schistosoma mansoni mouse model used for in vivo antischistosomal drug discovery, highlighting differences in the experimental set-up across research groups and their potential influence on experimental results. We explore the pharmacokinetic/pharmacodynamic relationship of selected drug candidates, showcasing opportunities to improve the drug profile to accelerate the transition from the early drug discovery phase to new clinical candidates.

Scratching the Itch: Updated Perspectives on the Schistosomes Responsible for Swimmer's Itch around the World

Although most studies of digenetic trematodes of the family Schistosomatidae dwell on representatives causing human schistosomiasis, the majority of the 130 identified species of schistosomes infect birds or non-human mammals. The cercariae of many of these species can cause swimmer's itch when they penetrate human skin. Recent years have witnessed a dramatic increase in our understanding of schistosome diversity, now encompassing 17 genera with eight more lineages awaiting description. Collectively, schistosomes exploit 16 families of caenogastropod or heterobranch gastropod intermediate hosts. Basal lineages today are found in marine gastropods and birds, but subsequent diversification has largely taken place in freshwater, with some reversions to marine habitats. It seems increasingly likely that schistosomes have on two separate occasions colonized mammals. Swimmer's itch is a complex zoonotic disease manifested through several different routes of transmission involving a diversity of different host species. Swimmer's itch also exemplifies the value of adopting the One Health perspective in understanding disease transmission and abundance because the schistosomes involved have complex life cycles that interface with numerous species and abiotic components of their aquatic environments. Given the progress made in revealing their diversity and biology, and the wealth of questions posed by itch-causing schistosomes, they provide excellent models for implementation of long-term interdisciplinary studies focused on issues pertinent to disease ecology, the One Health paradigm, and the impacts of climate change, biological invasions and other environmental perturbations.

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.

Worms and bugs of the gut: the search for diagnostic signatures using barcoding, and metagenomics-metabolomics

Gastrointestinal (GI) helminth infections cause significant morbidity in both humans and animals worldwide. Specific and sensitive diagnosis is central to the surveillance of such infections and to determine the effectiveness of treatment strategies used to control them. In this article, we: (i) assess the strengths and limitations of existing methods applied to the diagnosis of GI helminth infections of humans and livestock; (ii) examine high-throughput sequencing approaches, such as targeted molecular barcoding and shotgun sequencing, as tools to define the taxonomic composition of helminth infections; and (iii) discuss the current understanding of the interactions between helminths and microbiota in the host gut. Stool-based diagnostics are likely to serve as an important tool well into the future; improved diagnostics of helminths and their environment in the gut may assist the identification of biomarkers with the potential to define the health/disease status of individuals and populations, and to identify existing or emerging anthelmintic resistance.

A genetic TRP down the channel to praziquantel resistance

The anthelmintic praziquantel (PZQ) is an essential tool in controlling schistosomiasis, so reports of reduced PZQ efficacy are of great public health concern. Le Clec'h et al. recently identified a gene responsible for PZQ resistance in experimentally selected resistant Schistosoma mansoni. The importance of this locus in natural infections remains to be established.

An Overview of Transcriptional Responses of Schistosome-Susceptible (M line) or -Resistant (BS-90) Biomphalaria glabrata Exposed or Not to Schistosoma mansoni Infection

Background

We seek to provide a comprehensive overview of transcriptomics responses of immune-related features of the gastropod Biomphalaria glabrata (Bg) following exposure to Schistosoma mansoni (Sm), a trematode causing human schistosomiasis. Responses of schistosome-susceptible (M line, or SUS) and -resistant (BS-90, or RES) Bg strains were characterized following exposure to Sm for 0.5, 2, 8 or 40 days post-exposure (dpe).

Methods

RNA-Seq and differential expression analysis were undertaken on 56 snails from 14 groups. We considered 7 response categories: 1) constitutive resistance factors; 2) constitutive susceptibility factors; 3) generalized stress responses; 4) induced resistance factors; 5) resistance factors suppressed in SUS snails; 6) suppressed/manipulated factors in SUS snails; and 7) tolerance responses in SUS snails. We also undertook a gene co-expression network analysis. Results from prior studies identifying schistosome resistance/susceptibility factors were examined relative to our findings.

Results

A total of 792 million paired-end reads representing 91.2% of the estimated 31,985 genes in the Bg genome were detected and results for the 7 categories compiled and highlighted. For both RES and SUS snails, a single most supported network of genes with highly correlated expression was found.

Conclusions

1) Several constitutive differences in gene expression between SUS and RES snails were noted, the majority over-represented in RES; 2) There was little indication of a generalized stress response shared by SUS and RES snails at 0.5 or 2 dpe; 3) RES snails mounted a strong, multi-faceted response by 0.5 dpe that carried over to 2 dpe; 4) The most notable SUS responses were at 40 dpe, in snails shedding cercariae, when numerous features were either strongly down-regulated indicative of physiological distress or parasite manipulation, or up-regulated, suggestive of tolerance or survival-promoting effects; 5) Of 55 genes previously identified in genome wide mapping studies, 29 (52.7%) were responsive to Sm, as were many familiar resistance-associated genes (41.0%) identified by other means; 6) Both network analysis and remarkably specific patterns of expression of lectins and G protein-coupled receptors in categories 4, 6 and 7 were indicative of orchestrated responses of different suites of genes in SUS or RES snails following exposure to Sm.

Chromosome-level genome of Schistosoma haematobium underpins genome-wide explorations of molecular variation

Urogenital schistosomiasis is caused by the blood fluke Schistosoma haematobium and is one of the most neglected tropical diseases worldwide, afflicting > 100 million people. It is characterised by granulomata, fibrosis and calcification in urogenital tissues, and can lead to increased susceptibility to HIV/AIDS and squamous cell carcinoma of the bladder. To complement available treatment programs and break the transmission of disease, sound knowledge and understanding of the biology and ecology of S. haematobium is required. Hybridisation/introgression events and molecular variation among members of the S. haematobium-group might effect important biological and/or disease traits as well as the morbidity of disease and the effectiveness of control programs including mass drug administration. Here we report the first chromosome-contiguous genome for a well-defined laboratory line of this blood fluke. An exploration of this genome using transcriptomic data for all key developmental stages allowed us to refine gene models (including non-coding elements) and annotations, discover ‘new’ genes and transcription profiles for these stages, likely linked to development and/or pathogenesis. Molecular variation within S. haematobium among some geographical locations in Africa revealed unique genomic ‘signatures’ that matched species other than S. haematobium, indicating the occurrence of introgression events. The present reference genome (designated Shae.V3) and the findings from this study solidly underpin future functional genomic and molecular investigations of S. haematobium and accelerate systematic, large-scale population genomics investigations, with a focus on improved and sustained control of urogenital schistosomiasis.

More than 100 million people are infected with the carcinogenic blood fluke Schistosoma haematobium, the aetiological agent of urogenital schistosomiasis—a neglected tropical disease (NTD). In spite of its major significance, little is known about this fluke, its interactions with the human and snail intermediate hosts and the pathogenesis of the urogenital form of schistosomiasis at the molecular and biochemical levels. To enable research in these areas, we report the first chromosome-level genome and markedly enhanced gene models for S. haematobium. Comparative genomic analyses also reveal evidence of past introgression events between or among closely related schistosome species. This present reference genome for S. haematobium and the findings from this study should underpin future functional genomic and molecular investigations of S. haematobium and accelerate systematic, large-scale population genomics investigations, with a focus on improved control of urogenital schistosomiasis.

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

Whole-genome sequencing of Schistosoma mansoni reveals extensive diversity with limited selection despite mass drug administration

Control and elimination of the parasitic disease schistosomiasis relies on mass administration of praziquantel. Whilst these programmes reduce infection prevalence and intensity, their impact on parasite transmission and evolution is poorly understood. Here we examine the genomic impact of repeated mass drug administration on Schistosoma mansoni populations with documented reduced praziquantel efficacy. We sequenced whole-genomes of 198 S. mansoni larvae from 34 Ugandan children from regions with contrasting praziquantel exposure. Parasites infecting children from Lake Victoria, a transmission hotspot, form a diverse panmictic population. A single round of treatment did not reduce this diversity with no apparent population contraction caused by long-term praziquantel use. We find evidence of positive selection acting on members of gene families previously implicated in praziquantel action, but detect no high frequency functionally impactful variants. As efforts to eliminate schistosomiasis intensify, our study provides a foundation for genomic surveillance of this major human parasite.