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.

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.

Achieving Elimination as a Public Health Problem for Schistosoma mansoni and S. haematobium: When Is Community-Wide Treatment Required?

The World Health Organization (WHO) has set elimination as a public health problem (EPHP) as a goal for schistosomiasis. As the WHO treatment guidelines for schistosomiasis are currently under revision, we investigate whether school-based or community-wide treatment strategies are required for achieving the EPHP goal. In low- to moderate-transmission settings with good school enrolment, we find that school-based treatment is sufficient for achieving EPHP. However, community-wide treatment is projected to be necessary in certain high-transmission settings as well as settings with low school enrolment. Hence, the optimal treatment strategy depends on setting-specific factors such as the species present, prevalence prior to treatment, and the age profile of infection.

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.

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.

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.

Interrupting seasonal transmission of Schistosoma haematobium and control of soil-transmitted helminthiasis in northern and central Côte d'Ivoire: a SCORE study protocol.

BACKGROUND: To achieve a world free of schistosomiasis, the objective is to scale up control and elimination efforts in all endemic countries. Where interruption of transmission is considered feasible, countries are encouraged to implement a comprehensive intervention package, including preventive chemotherapy, information, education and communication (IEC), water, sanitation and hygiene (WASH), and snail control. In northern and central Côte d'Ivoire, transmission of Schistosoma haematobium is seasonal and elimination might be achieved. In a cluster-randomised trial, we will assess different treatment schemes to interrupt S. haematobium transmission and control soil-transmitted helminthiasis over a 3-year period. We will compare the impact of (i) arm A: annual mass drug administration (MDA) with praziquantel and albendazole before the peak schistosomiasis transmission season; (ii) arm B: annual MDA after the peak schistosomiasis transmission season; (iii) arm C: two yearly treatments before and after peak schistosomiasis transmission; and (iv) arm D: annual MDA before peak schistosomiasis transmission, coupled with chemical snail control using niclosamide. METHODS/DESIGN: The prevalence and intensity of S. haematobium and soil-transmitted helminth infections will be assessed using urine filtration and Kato-Katz thick smears, respectively, in six administrative regions in northern and central parts of Côte d'Ivoire. Once a year, urine and stool samples will be collected and examined from 50 children aged 5-8 years, 100 children aged 9-12 years and 50 adults aged 20-55 years in each of 60 selected villages. Changes in S. haematobium and soil-transmitted helminth prevalence and intensity will be assessed between years and stratified by intervention arm. In the 15 villages randomly assigned to intervention arm D, intermediate host snails will be collected three times per year, before niclosamide is applied to the selected freshwater bodies. The snail abundance and infection rates over time will allow drawing inference on the force of transmission. DISCUSSION: This cluster-randomised intervention trial will elucidate whether in an area with seasonal transmission, the four different treatment schemes can interrupt S. haematobium transmission and control soil-transmitted helminthiasis. Lessons learned will help to guide schistosomiasis control and elimination programmes elsewhere in Africa. TRIAL REGISTRATION: ISRCTN ISRCTN10926858 . Registered 21 December 2016. Retrospectively registered.

Sensory Protein Kinase Signaling in Schistosoma mansoni Cercariae: Host Location and Invasion.

Schistosoma mansoni cercariae display specific behavioral responses to abiotic/biotic stimuli enabling them to locate and infect the definitive human host. Here we report the effect of such stimulants on signaling pathways of cercariae in relation to host finding and invasion. Cercariae exposed to various light/temperature regimens displayed modulated protein kinase C (PKC), extracellular signal-regulated kinase (ERK) and p38 mitogen-activated protein kinase (p38 MAPK) activities, with distinct responses at 37 °C and intense light/dark, when compared to 24 °C under normal light. Kinase activities were localized to regions including the oral sensory papillae, acetabular ducts, tegument, acetabular glands, and nervous system. Furthermore, linoleic acid modulated PKC and ERK activities concurrent with the temporal release of acetabular gland components. Attenuation of PKC, ERK, and p38 MAPK activities significantly reduced gland component release, particularly in response to linoleic acid, demonstrating the importance of these signaling pathways to host penetration mechanisms.

Diagnostics for schistosomiasis in Africa and Arabia: a review of present options in control and future needs for elimination.

Within the World Health Organization 2012-2020 roadmap for control and elimination of schistosomiasis, the scale-up of mass drug administration with praziquantel is set to change the epidemiological landscape across Africa and Arabia. Central in measuring progress is renewed emphasis upon diagnostics which operate at individual, community and environmental levels by assessing reductions in disease, infections and parasite transmission. However, a fundamental tension is revealed between levels for present diagnostic tools, and methods applied in control settings are not necessarily adequate for application in elimination scenarios. Indeed navigating the transition from control to elimination needs careful consideration and planning. In the present context of control, we review current options for diagnosis of schistosomiasis at different levels, highlighting several strengths and weaknesses therein. Future challenges in elimination are raised and we propose that more cost-effective diagnostics and clinical staging algorithms are needed. Using the Kingdom of Saudi Arabia as a contemporary example, embedding new diagnostic methods within the primary care health system is discussed with reference to both urogenital and intestinal schistosomiasis.

Accelerating snail vector genomics.

The three most important genera of snails for the transmission of schistosomes are Bulinus, Biomphalaria and Oncomelania. Each of these genera, found in two distantly related families, includes species that act as the intermediate host for one of the three most widespread schistosome species infecting humans, Schistosoma haematobium, S. mansoni and S. japonicum, respectively. An important step in the fight against schistosomiasis in Asia has been taken with the publication of the article "Chromosome-level genome assembly of Oncomelania hupensis: the intermediate snail host of Schistosoma japonicum", which means that genomes for all three major genera, including species across three continents, are now available in the public domain. This includes the first genomes of African snail vectors, namely Biomphalaria sudanica, Bi. pfeifferi and Bulinus truncatus, as well as high-quality chromosome level assemblies for South American Bi. glabrata. Most importantly, the wealth of new genomic and transcriptomic data is helping to establish the specific molecular mechanisms that underly compatibility between snails and their schistosomes, which although diverse and complex, may help to identify potential targets dictating host parasite interactions that can be utilised in future transmission control strategies. This new work on Oncomelania hupensis and indeed studies on other snail vectors, which provide deep insights into the genome, will stimulate research that may well lead to new and much needed control interventions.