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

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

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

Schistosoma haematobium and Schistosoma bovis first generation hybrids undergo gene expressions changes consistent with species compatibility and heterosis

When two species hybridize, the two parental genomes are brought together and some alleles might interact for the first time. To date, the extent of the transcriptomic changes in first hybrid generations, along with their functional outcome constitute an important knowledge gap, especially in parasite species. Here we explored the molecular and functional outcomes of hybridization in first-generation hybrids between the blood fluke parasites Schistosoma haematobium and S. bovis. Through a transcriptomic approach, we measured gene expression in both parental species and hybrids. We described and quantified expression profiles encountered in hybrids along with the main biological processes impacted. Up to 7,100 genes fell into a particular hybrid expression profile (intermediate between the parental expression levels, over-expressed, under-expressed, or expressed like one of the parental lines). Most of these genes were different depending on the direction of the parental cross (S. bovis mother and S. haematobium father or the reverse) and depending on the sex. For a given sex and cross direction, the vast majority of genes were hence unassigned to a hybrid expression profile: either they were differentially expressed genes but not typical of any hybrid expression profiles or they were not differentially expressed neither between hybrids and parental lines nor between parental lines. The most prevalent profile of gene expression in hybrids was the intermediate one (24% of investigated genes). These results suggest that transcriptomic compatibility between S. haematobium and S. bovis remains quite high. We also found support for an over-dominance model (over- and under-expressed genes in hybrids compared to parental lines) potentially associated with heterosis. In females in particular, processes such as reproductive processes, metabolism and cell interactions as well as signaling pathways were indeed affected. Our study hence provides new insight on the biology of Schistosoma hybrids with evidences supporting compatibility and heterosis.

Schistosomiasis diagnosis: Challenges and opportunities for elimination

OVERVIEW

The roadmap adopted by the World Health Organization (WHO) for eliminating neglected tropical diseases aims to eliminate schistosomiasis, as a public health concern, by 2030. While progress has been made towards reducing schistosomiasis morbidity control in several sub-Saharan African countries, there is still more that needs to be done. Proper surveillance using accurate diagnostics with acceptable sensitivity and specificity is essential for evaluating the success of all efforts against schistosomiasis. Microscopy, despite its low sensitivity, remains the gold standard approach for diagnosing the disease. Although many efforts have been made to develop new diagnostics based on circulating parasite proteins, genetic markers, schistosome egg morphology, and their paramagnetic properties, none has been robust enough to replace microscopy. This review highlights common diagnostic approaches for detecting schistosomiasis in field and clinical settings, major challenges, and provides new and novel opportunities and diagnosis pathways that will be critical in supporting elimination of schistosomiasis.

METHODS

We searched for relevant and reliable published literature from PubMed, Scopus, google scholar, and Web of science. The search strategies were primarily determined by subtopic, and hence the following words were used (schistosom*, diagnosis, Kato-Katz, antibody test, circulating antigen, POC-CCA, UCP-LF-CAA, molecular diagnostics, nucleic acid amplification test, microfluidics, lab-on a disk, lab-on chip, recombinase polymerase amplification (RPA), LAMP, portable sequencer, nanobody test, identical multi-repeat sequences, diagnostic TPPs, REASSURED, extraction free), and Boolean operators AND and/OR were used to refine the searching capacity. Due to the global public health nature of schistosomiasis, we also searched for reliable documents, reports, and research papers published by international health organizations, World Health Organization (WHO), and Center for Disease control and Elimination.

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.

A single locus determines praziquantel response in Schistosoma mansoni

We previously performed a genome-wide association study (GWAS) to identify the genetic basis of praziquantel (PZQ) response in schistosomes, identifying two quantitative trait loci situated on chromosomes 2 and 3. We reanalyzed this GWAS using the latest (version 10) genome assembly showing that a single locus on chromosome 3, rather than two independent loci, determines drug response. These results reveal that PZQ response is monogenic and demonstrates the importance of high-quality genomic information.

Trematode Genomics and Proteomics

Trematode infections stand out as one of the frequently overlooked tropical diseases, despite their wide global prevalence and remarkable capacity to parasitize diverse host species and tissues. Furthermore, these parasites hold significant socio-economic, medical, veterinary and agricultural implications. Over the past decades, substantial strides have been taken to bridge the information gap concerning various "omic" tools, such as proteomics and genomics, in this field. In this edition of the book, we highlight recent progress in genomics and proteomics concerning trematodes with a particular focus on the advances made in the past 5 years. Additionally, we present insights into cutting-edge technologies employed in studying trematode biology and shed light on the available resources for exploring the molecular facets of this particular group of parasitic helminths.

Chromosome-scale genome of the human blood fluke Schistosoma mekongi and its implications for public health

BACKGROUND

Schistosoma mekongi is a human blood fluke causing schistosomiasis that threatens approximately 1.5 million humans in the world. Nonetheless, the limited available S. mekongi genomic resources have hindered understanding of its biology and parasite-host interactions for disease management and pathogen control. The aim of our study was to integrate multiple technologies to construct a high-quality chromosome-level assembly of the S. mekongi genome.

METHODS

The reference genome for S. mekongi was generated through integrating Illumina, PacBio sequencing, 10 × Genomics linked-read sequencing, and high-throughput chromosome conformation capture (Hi-C) methods. In this study, we conducted de novo assembly, alignment, and gene prediction to assemble and annotate the genome. Comparative genomics allowed us to compare genomes across different species, shedding light on conserved regions and evolutionary relationships. Additionally, our transcriptomic analysis focused on genes associated with parasite-snail interactions in S. mekongi infection. We employed gene ontology (GO) enrichment analysis for functional annotation of these genes.

RESULTS

In the present study, the S. mekongi genome was both assembled into 8 pseudochromosomes with a length of 404 Mb, with contig N50 and scaffold N50 lengths of 1168 kb and 46,759 kb, respectively. We detected that 43% of the genome consists of repeat sequences and predicted 9103 protein-coding genes. We also focused on proteases, particularly leishmanolysin-like metalloproteases (M8), which are crucial in the invasion of hosts by 12 flatworm species. Through phylogenetic analysis, it was discovered that the M8 gene exhibits lineage-specific amplification among the genus Schistosoma. Lineage-specific expansion of M8 was observed in blood flukes. Additionally, the results of the RNA-seq revealed that a mass of genes related to metabolic and biosynthetic processes were up-regulated, which might be beneficial for cercaria production.

CONCLUSIONS

This study delivers a high-quality, chromosome-scale reference genome of S. mekongi, enhancing our understanding of the divergence and evolution of Schistosoma. The molecular research conducted here also plays a pivotal role in drug discovery and vaccine development. Furthermore, our work greatly advances the understanding of host-parasite interactions, providing crucial insights for schistosomiasis intervention strategies.

Emerging questions on the mechanisms and dynamics of 3D genome evolution in spiralians

Information on how 3D genome topology emerged in animal evolution, how stable it is during development, its role in the evolution of phenotypic novelties and how exactly it affects gene expression is highly debated. So far, data to address these questions are lacking with the exception of a few key model species. Several gene regulatory mechanisms have been proposed, including scenarios where genome topology has little to no impact on gene expression, and vice versa. The ancient and diverse clade of spiralians may provide a crucial testing ground for such mechanisms. Sprialians have followed distinct evolutionary trajectories, with some clades experiencing genome expansions and/or large-scale genome rearrangements, and others undergoing genome contraction, substantially impacting their size and organisation. These changes have been associated with many phenotypic innovations in this clade. In this review, we describe how emerging genome topology data, along with functional tools, allow for testing these scenarios and discuss their predicted outcomes.

A single locus determines praziquantel response in Schistosoma mansoni

We previously performed a genome-wide association study (GWAS) to identify the genetic basis of praziquantel (PZQ) response in schistosomes, identifying two quantitative trait loci (QTL) situated on chromosome 2 and chromosome 3. We reanalyzed this GWAS using the latest (v10) genome assembly showing that a single locus on chromosome 3, rather than two independent loci, determines drug response. These results reveal that praziquantel response is monogenic and demonstrates the importance of high-quality genomic information.

Protective human IgE responses are promoted by comparable life-cycle dependent Tegument Allergen-Like expression in Schistosoma haematobium and Schistosoma mansoni infection

Schistosoma haematobium is the most prevalent of the human-infecting schistosome species, causing significant morbidity in endemically exposed populations. Despite this, it has been relatively understudied compared to its fellow species, S. mansoni. Here we provide the first comprehensive characterization of the S. haematobium Tegument Allergen-Like protein family, a key protein family directly linked to protective immunity in S. mansoni infection. Comparable with observations for S. mansoni, parasite phylogenetic analysis and relative gene expression combined with host serological analysis support a cross-reactive relationship between S. haematobium TAL proteins, exposed to the host immune system as adult worms die, and closely related proteins, exposed during penetration by the infecting cercarial and early schistosomulae stages. Specifically, our results strengthen the evidence for host immunity driven by cross-reactivity between family members TAL3 and TAL5, establishing it for the first time for S. haematobium infection. Furthermore, we build upon this relationship to include the involvement of an additional member of the TAL protein family, TAL11 for both schistosome species. Finally, we show a close association between experience of infection and intensity of transmission and the development of protective IgE responses to these antigens, thus improving our knowledge of the mechanisms by which protective host immune responses develop. This knowledge will be critical in understanding how control efforts such as mass drug administration campaigns influence the development of host immunity and subsequent patterns of infection and disease within endemic populations.

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.

Schistosomes contain divergent ligand-gated ion channels with an atypical Cys-loop motif

Ligand-gated ion channels (LGICs) are important regulators of neuromuscular function, making them attractive antiparasitic drug targets. While roundworm LGICs are targeted by several anthelmintic classes, flatworm LGICs are less studied. Chromosome-level genome assemblies have recently been released for Schistosoma flatworm species that cause the disease schistosomiasis. These have allowed us to comprehensively predict schistosome LGICs, adding to prior annotations. Analysis of LGIC sequences revealed a clade of receptors lacking cysteines at the eponymous Cys-loop region of the channel. Since these atypical channels are divergent from mammalian LGICs, they may be promising targets to treat diseases caused by parasitic flatworms.

Annotation of G-Protein Coupled Receptors in the Genomes of Parasitic Blood Flukes

Infection with Schistosoma parasitic flatworms ( Schistosoma haematobium, Schistosoma mansoni and Schistosoma japonicum ) causes the neglected tropical disease schistosomiasis. There is a need to identify new chemotherapies to treat these parasites, and G-protein coupled receptors (GPCRs) are a logical druggable targets to explore given they control key aspects of schistosome biology such as neuromuscular function and reproduction. Updated chromosome level genome assemblies for each of the three major species have recently been released. However, studies on these GPCRs require accurate, updated genome annotations. Here, we have re-annotated the GPCRs present in each of the three major schistosome species.

Improving helminth genome resources in the post-genomic era

Rapid advancement in high-throughput sequencing and analytical approaches has seen a steady increase in the generation of genomic resources for helminth parasites. Now, helminth genomes and their annotations are a cornerstone of numerous efforts to compare genetic and transcriptomic variation, from single cells to populations of globally distributed parasites, to genome modifications to understand gene function. Our understanding of helminths is increasingly reliant on these genomic resources, which are primarily static once published and vary widely in quality and completeness between species. This article seeks to highlight the cause and effect of this variation and argues for the continued improvement of these genomic resources - even after their publication - which is necessary to provide a more accurate and complete understanding of the biology of these important pathogens.

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.

Genomic evidence of contemporary hybridization between Schistosoma species

Hybridization between different species of parasites is increasingly being recognised as a major public and veterinary health concern at the interface of infectious diseases biology, evolution, epidemiology and ultimately control. Recent research has revealed that viable hybrids and introgressed lineages between Schistosoma spp. are prevalent across Africa and beyond, including those with zoonotic potential. However, it remains unclear whether these hybrid lineages represent recent hybridization events, suggesting hybridization is ongoing, and/or whether they represent introgressed lineages derived from ancient hybridization events. In human schistosomiasis, investigation is hampered by the inaccessibility of adult-stage worms due to their intravascular location, an issue which can be circumvented by post-mortem of livestock at abattoirs for Schistosoma spp. of known zoonotic potential. To characterise the composition of naturally-occurring schistosome hybrids, we performed whole-genome sequencing of 21 natural livestock infective schistosome isolates. To facilitate this, we also assembled a de novo chromosomal-scale draft assembly of Schistosoma curassoni. Genomic analyses identified isolates of S. bovis, S. curassoni and hybrids between the two species, all of which were early generation hybrids with multiple generations found within the same host. These results show that hybridization is an ongoing process within natural populations with the potential to further challenge elimination efforts against schistosomiasis.

Schistosomiasis is a chronic and debilitating major neglected tropical disease affecting both humans and livestock. Increasingly, zoonotic spillover of livestock infections, facilitated by hybridization between different Schistosoma species, is increasingly being recognised as a risk to human health. Multiple surveys conducted within endemic regions have found a high prevalence of these hybrid lineages. However, it is often unclear whether these lineages are derived from recent hybridization events, suggesting hybridization is ongoing and may be linked to anthropogenic environmental change, or simply indicators of introgression from ancient hybridization events. To understand the origin and evolution of these hybrid lineages, we produced a chromosomal-scale assembly of Schistosoma curassoni and performed whole-genome sequencing of 21 natural livestock-infective S. curassoni, S. bovis and hybridized schistosome isolates, including multi-stage sampling from the same hosts. Our analyses exclusively identified early generation hybrid lineages, including multiple unrelated generations within the same hosts, suggesting that these hybrids are viable and derived from multiple independent hybridization events.

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.