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.

Diverging patterns of introgression from Schistosoma bovis across S. haematobium African lineages.

Hybridization is a fascinating evolutionary phenomenon that raises the question of how species maintain their integrity. Inter-species hybridization occurs between certain Schistosoma species that can cause important public health and veterinary issues. In particular hybrids between Schistosoma haematobium and S. bovis associated with humans and animals respectively are frequently identified in Africa. Recent genomic evidence indicates that some S. haematobium populations show signatures of genomic introgression from S. bovis. Here, we conducted a genomic comparative study and investigated the genomic relationships between S. haematobium, S. bovis and their hybrids using 19 isolates originating from a wide geographical range over Africa, including samples initially classified as S. haematobium (n = 11), S. bovis (n = 6) and S. haematobium x S. bovis hybrids (n = 2). Based on a whole genomic sequencing approach, we developed 56,181 SNPs that allowed a clear differentiation of S. bovis isolates from a genomic cluster including all S. haematobium isolates and a natural S. haematobium-bovis hybrid. All the isolates from the S. haematobium cluster except the isolate from Madagascar harbored signatures of genomic introgression from S. bovis. Isolates from Corsica, Mali and Egypt harbored the S. bovis-like Invadolysin gene, an introgressed tract that has been previously detected in some introgressed S. haematobium populations from Niger. Together our results highlight the fact that introgression from S. bovis is widespread across S. haematobium and that the observed introgression is unidirectional.

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.

Precision mapping of snail habitat provides a powerful indicator of human schistosomiasis transmission.

Recently, the World Health Organization recognized that efforts to interrupt schistosomiasis transmission through mass drug administration have been ineffective in some regions; one of their new recommended strategies for global schistosomiasis control emphasizes targeting the freshwater snails that transmit schistosome parasites. We sought to identify robust indicators that would enable precision targeting of these snails. At the site of the world's largest recorded schistosomiasis epidemic-the Lower Senegal River Basin in Senegal-intensive sampling revealed positive relationships between intermediate host snails (abundance, density, and prevalence) and human urogenital schistosomiasis reinfection (prevalence and intensity in schoolchildren after drug administration). However, we also found that snail distributions were so patchy in space and time that obtaining useful data required effort that exceeds what is feasible in standard monitoring and control campaigns. Instead, we identified several environmental proxies that were more effective than snail variables for predicting human infection: the area covered by suitable snail habitat (i.e., floating, nonemergent vegetation), the percent cover by suitable snail habitat, and size of the water contact area. Unlike snail surveys, which require hundreds of person-hours per site to conduct, habitat coverage and site area can be quickly estimated with drone or satellite imagery. This, in turn, makes possible large-scale, high-resolution estimation of human urogenital schistosomiasis risk to support targeting of both mass drug administration and snail control efforts.

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.

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.

Nanopore Sequencing Resolves Elusive Long Tandem-Repeat Regions in Mitochondrial Genomes.

Long non-coding, tandem-repetitive regions in mitochondrial (mt) genomes of many metazoans have been notoriously difficult to characterise accurately using conventional sequencing methods. Here, we show how the use of a third-generation (long-read) sequencing and informatic approach can overcome this problem. We employed Oxford Nanopore technology to sequence genomic DNAs from a pool of adult worms of the carcinogenic parasite, Schistosoma haematobium, and used an informatic workflow to define the complete mt non-coding region(s). Using long-read data of high coverage, we defined six dominant mt genomes of 33.4 kb to 22.6 kb. Although no variation was detected in the order or lengths of the protein-coding genes, there was marked length (18.5 kb to 7.6 kb) and structural variation in the non-coding region, raising questions about the evolution and function of what might be a control region that regulates mt transcription and/or replication. The discovery here of the largest tandem-repetitive, non-coding region (18.5 kb) in a metazoan organism also raises a question about the completeness of some of the mt genomes of animals reported to date, and stimulates further explorations using a Nanopore-informatic workflow.

Multihost Transmission of Schistosoma mansoni in Senegal, 2015-2018.

In West Africa, Schistosoma spp. are capable of infecting multiple definitive hosts, a lifecycle feature that may complicate schistosomiasis control. We characterized the evolutionary relationships among multiple Schistosoma mansoni isolates collected from snails (intermediate hosts), humans (definitive hosts), and rodents (definitive hosts) in Senegal. On a local scale, diagnosis of S. mansoni infection ranged 3.8%-44.8% in school-aged children, 1.7%-52.6% in Mastomys huberti mice, and 1.8%-7.1% in Biomphalaria pfeifferi snails. Our phylogenetic framework confirmed the presence of multiple S. mansoni lineages that could infect both humans and rodents; divergence times of these lineages varied (0.13-0.02 million years ago). We propose that extensive movement of persons across West Africa might have contributed to the establishment of these various multihost S. mansoni clades. High S. mansoni prevalence in rodents at transmission sites frequented by humans further highlights the implications that alternative hosts could have on future public health interventions.

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.

Genetic diversity of Biomphalaria pfeifferi, the intermediate host of Schistosoma mansoni in Shamva district, Zimbabwe: role on intestinal schistosomiasis transmission.

The fresh water snail Biomphalaria pfeifferi is the intermediate host for Schistosoma mansoni, which causes human intestinal schistosomiasis in Zimbabwe. Despite the medical importance of this intermediate host, there are no current data on its molecular characterization in Zimbabwe. In 2016, human water contact sites were identified in four communities in Madziwa area, Shamva district, Zimbabwe. The survey sites were recorded and mapped using a global positioning system. A 655 bp region of the mitochondrial cytochrome oxidase subunit I gene was amplified in 70 B. pfeifferi snails. The sequence data were analysed to determine the relationships between the individual snails, their inter, intra population diversity and structure. Overall, four unique cox1 haplotypes, with a haplotype diversity of 0.608, were identified in the snails. One haplotype spanned across most of the sites. There was no clear geographical clustering of haplotypes. The mean diversity among the haplotypes was very low (0.009), while the net divergence among the collection sites ranged from 0.000 to 0.026. The diversity within and between the sites was 0.017 and 0.012 respectively. This data advances our knowledge of the understanding of the population structure of B. pfeifferi in Madziwa area, Zimbabwe, with the high occurrence of one haplotype indicating the possibility of a recent bottleneck followed by population expansion. The population genetic structure of B. pfeifferi snails described here has provided an opportunity to investigate the contribution of snail genetics to variation in disease burden; and development of control strategies that exploit genetic differences in susceptibility to parasites.

An update on non-invasive urine diagnostics for human-infecting parasitic helminths: what more could be done and how?

Reliable diagnosis of human helminth infection(s) is essential for ongoing disease surveillance and disease elimination. Current WHO-recommended diagnostic assays are unreliable in low-endemic near-elimination settings and typically involve the invasive, onerous and potentially hazardous sampling of bodily fluids such as stool and blood, as well as tissue via biopsy. In contrast, diagnosis by use of non-invasive urine sampling is generally painless, more convenient and low risk. It negates the need for specialist staff, can usually be obtained immediately upon request and is better accepted by patients. In some instances, urine-based diagnostic assays have also been shown to provide a more reliable diagnosis of infection when compared to traditional methods that require alternative and more invasive bodily samples, particularly in low-endemicity settings. Given these relative benefits, we identify and review current research literature to evaluate whether non-invasive urine sampling is currently exploited to its full potential in the development of diagnostic tools for human helminthiases. Though further development, assessment and validation are needed before their routine use in control programmes, low-cost, rapid and reliable assays capable of detecting transrenal helminth-derived antigens and cell-free DNA show excellent promise for future use at the point-of-care in high-, medium- and even low-endemicity elimination settings.

Analytical and Clinical Assessment of a Portable, Isothermal Recombinase Polymerase Amplification (RPA) Assay for the Molecular Diagnosis of Urogenital Schistosomiasis.

Accurate diagnosis of urogenital schistosomiasis is crucial for disease surveillance and control. Routine diagnostic methods, however, lack sensitivity when assessing patients with low levels of infection still able to maintain pathogen transmission. Therefore, there is a need for highly sensitive diagnostic tools that can be used at the point-of-care in endemic areas. Recombinase polymerase amplification (RPA) is a rapid and sensitive diagnostic tool that has been used to diagnose several pathogens at the point-of-care. Here, the analytical performance of a previously developed RPA assay (RT-ShDra1-RPA) targeting the Schistosoma haematobium Dra1 genomic region was assessed using commercially synthesised S. haematobium Dra1 copies and laboratory-prepared samples spiked with S. haematobium eggs. Clinical performance was also assessed by comparing diagnostic outcomes with that of a reference diagnostic standard, urine-egg microscopy. The RT-ShDra1-RPA was able to detect 1 × 101 copies of commercially synthesised Dra1 DNA as well as one S. haematobium egg within laboratory-spiked ddH2O samples. When compared with urine-egg microscopy, the overall sensitivity and specificity of the RT-ShDra1-RPA assay was 93.7% (±88.7-96.9) and 100% (±69.1-100), respectively. Positive and negative predictive values were 100% (±97.5-100) and 50% (±27.2-72.8), respectively. The RT-ShDra1-RPA therefore shows promise as a rapid and highly sensitive diagnostic tool able to diagnose urogenital schistosomiasis at the point-of-care.

Development of a Molecular Snail Xenomonitoring Assay to Detect Schistosoma haematobium and Schistosoma bovis Infections in their Bulinus Snail Hosts.

Schistosomiasis, a neglected tropical disease of medical and veterinary importance, transmitted through specific freshwater snail intermediate hosts, is targeted for elimination in several endemic regions in sub-Saharan Africa. Multi-disciplinary methods are required for both human and environmental diagnostics to certify schistosomiasis elimination when eventually reached. Molecular xenomonitoring protocols, a DNA-based detection method for screening disease vectors, have been developed and trialed for parasites transmitted by hematophagous insects, such as filarial worms and trypanosomes, yet few have been extensively trialed or proven reliable for the intermediate host snails transmitting schistosomes. Here, previously published universal and Schistosoma-specific internal transcribed spacer (ITS) rDNA primers were adapted into a triplex PCR primer assay that allowed for simple, robust, and rapid detection of Schistosoma haematobium and Schistosoma bovis in Bulinus snails. We showed this two-step protocol could sensitively detect DNA of a single larval schistosome from experimentally infected snails and demonstrate its functionality for detecting S. haematobium infections in wild-caught snails from Zanzibar. Such surveillance tools are a necessity for succeeding in and certifying the 2030 control and elimination goals set by the World Health Organization.

Schistosome Interactions within the Schistosoma haematobium Group, Malawi.

Molecular analysis of atypical schistosome eggs retrieved from children in Malawi revealed genetic interactions occurring between human (Schistosoma haematobium) and livestock (S. mattheei and S. bovis) schistosome species. Detection of hybrid schistosomes adds a notable new perspective to the epidemiology and control of urogenital schistosomiasis in central Africa.

Establishing the Production of Male Schistosoma mansoni Cercariae for a Controlled Human Infection Model.

To accelerate the development of novel vaccines for schistosomiasis, we set out to develop a human model for Schistosoma mansoni infection in healthy volunteers. During natural infections, female schistosomes produce eggs that give rise to morbidity. Therefore, we produced single-sex, male Schistosoma mansoni cercariae for human infection without egg production and associated pathology. Cercariae were produced in their intermediate snail hosts in accordance with the principles of good manufacturing practice (GMP). The application of GMP principles to an unconventional production process is a showcase for the controlled production of complex live challenge material in the European Union or under Food and Drug Administration guidance.

Development of a lateral flow recombinase polymerase assay for the diagnosis of Schistosoma mansoni infections.

Infection with Schistosoma mansoni causes intestinal schistosomiasis, a major health problem across Africa. The accurate diagnosis of intestinal schistosomiasis is vital to inform surveillance/control programs. Diagnosis mainly relies on microscopic detection of eggs in faecal samples but many factors affect sensitivity. Molecular diagnostics are sensitive and specific but application is limited as necessary infrastructure, financial resources and skilled personnel are often lacking in endemic settings. Recombinase Polymerase Amplification (RPA) is an isothermal DNA amplification/detection technology that is practical in nearly any setting. Here we developed a RPA lateral flow (LF) assay targeting the 28S rDNA region of S. mansoni. The 28S LF-RPA assay's lower limit of detection was 10pg DNA with the lower test parameters permitting sufficient amplification being 6 min and 25°C. Optimal assay parameters were 40-45°C and 10 min with an analytical sensitivity of 102 copies of DNA. Additionally the PCRD3 lateral flow detection cassettes proved more robust and sensitive compared to the Milenia HybriDetect strips. This 28S LF-RPA assay produces quick reproducible results that are easy to interpret, require little infrastructure and is a promising PON test for the field molecular diagnosis of intestinal schistosomiasis.

Urogenital schistosomiasis and hybridization between Schistosoma haematobium and Schistosoma bovis in adults living in Richard-Toll, Senegal.

Since the construction of the Diama Dam (1985), the epidemiology of schistosomiasis along the Senegal River Basin (SRB) has been extremely dynamic with outbreaks of both intestinal and urogenital schistosomiasis. In the early 2000s, technicians reported cases of suspected urogenital schistosomiasis in adults from the local hospital in Richard-Toll, Lower SRB. The genetic analysis of schistosome miracidia isolated from 11 patients in 2012 from two neighbourhoods (Campement and Gaya) of Richard-Toll confirmed infection with Schistosoma haematobium but also S. haematobium/S. bovis hybrids. Thirty-seven per cent of the miracidia were S. bovis/S. haematobium hybrids and 63% were pure S. haematobium. The data are discussed in relation to the ongoing dynamic epidemiology of the schistosomes in Senegal and the need to treat non-target individuals.

Paragonimus and paragonimiasis in West and Central Africa: unresolved questions.

Paragonimiasis, human lung fluke disease, is a foodborne anthropozoonosis caused by the trematodes assigned to Paragonimus and is regarded by the World Health Organization as a Neglected Tropical Disease (NTD). The life cycle of this medically important parasite centres on a complex freshwater biological community that includes two intermediate hosts: a mollusc and a decapod, usually a brachyuran. Although there is a perception that the biology, symptoms, diagnosis and treatment of Paragonimus is well understood, in reality, this is not the case, especially in Africa. Much remains unknown concerning the life-cycle of the parasite, its transmission, the current epidemiology of the disease, diagnosis and the effectiveness of treatment. Furthermore, cases of paragonimiasis may be misdiagnosed as resistant tuberculosis (TB) because of the similar pulmonary symptoms and no remission after anti TB therapy. The endemic foci of human paragonimiasis in Africa have been reported mainly in the forest zones of Upper Guinea (Liberia, Guinea and Ivory Coast) and Lower Guinea (Nigeria, Cameroon, Equatorial Guinea and Gabon). Despite the perceived medical importance of paragonimiasis, relatively little attention has been paid to this NTD since its discovery in Africa in the 1960s. This review focuses on the current understanding of the life cycle and transmission of Paragonimus in Africa, discusses its diagnosis and public health importance and highlights many outstanding gaps in the knowledge that still exist for this NTD.

Occurrence of Schistosoma bovis on Pemba Island, Zanzibar: implications for urogenital schistosomiasis transmission monitoring.

The causative agent of urogenital schistosomiasis, Schistosoma haematobium, was thought to be the only schistosome species transmitted through Bulinus snails on Unguja and Pemba Island (Zanzibar, United Republic of Tanzania). For insights into the environmental risk of S. haematobium transmission on Pemba Island, malacological surveys collecting Bulinus globosus and B. nasutus, two closely related potential intermediate hosts of S. haematobium were conducted across the island in November 2016. Of 1317 B. globosus/B. nasutus collected, seven B. globosus, identified through sequencing a DNA region of the mitochondrial cytochrome oxidase subunit 1 (cox1), were observed with patent infections assumed to be S. haematobium. However, when the collected cercariae were identified through sequencing a region of the cox1 and the nuclear internal transcribed spacer (ITS1 + 2), schistosomes from five of these B. globosus collected from a single locality were in fact S. bovis. The identified presence of S. bovis raises concerns for animal health on Pemba, and complicates future transmission monitoring of S. haematobium. These results show the pertinence for not only sensitive, but also species-specific markers to be used when identifying cercariae during transmission monitoring, and also provide the first molecular confirmation for B. globosus transmitting S. bovis in East Africa.

Whole genome amplification and exome sequencing of archived schistosome miracidia.

Adult schistosomes live in the blood vessels and cannot easily be sampled from humans, so archived miracidia larvae hatched from eggs expelled in feces or urine are commonly used for population genetic studies. Large collections of archived miracidia on FTA cards are now available through the Schistosomiasis Collection at the Natural History Museum (SCAN). Here we describe protocols for whole genome amplification of Schistosoma mansoni and Schistosome haematobium miracidia from these cards, as well as real time PCR quantification of amplified schistosome DNA. We used microgram quantities of DNA obtained for exome capture and sequencing of single miracidia, generating dense polymorphism data across the exome. These methods will facilitate the transition from population genetics, using limited numbers of markers to population genomics using genome-wide marker information, maximising the value of collections such as SCAN.