A systematically improved high quality genome and transcriptome of the human blood fluke Schistosoma mansoni

Genome-wide transcriptome analysis of Echinococcus multilocularis larvae and germinative cell cultures reveals genes involved in parasite stem cell function

The lethal zoonosis alveolar echinococcosis is caused by tumour-like growth of the metacestode stage of the tapeworm Echinococcus multilocularis within host organs. We previously demonstrated that metacestode proliferation is exclusively driven by somatic stem cells (germinative cells), which are the only mitotically active parasite cells that give rise to all differentiated cell types. The Echinococcus gene repertoire required for germinative cell maintenance and differentiation has not been characterised so far. We herein carried out Illumina sequencing on cDNA from Echinococcus metacestode vesicles, from metacestode tissue depleted of germinative cells, and from Echinococcus primary cell cultures. We identified a set of ~1,180 genes associated with germinative cells, which contained numerous known stem cell markers alongside genes involved in replication, cell cycle regulation, mitosis, meiosis, epigenetic modification, and nucleotide metabolism. Interestingly, we also identified 44 stem cell associated transcription factors that are likely involved in regulating germinative cell differentiation and/or pluripotency. By in situ hybridization and pulse-chase experiments, we also found a new general Echinococcus stem cell marker, EmCIP2Ah, and we provide evidence implying the presence of a slow cycling stem cell sub-population expressing the extracellular matrix factor Emkal1. RNA-Seq analyses on primary cell cultures revealed that metacestode-derived Echinococcus stem cells display an expanded differentiation capability and do not only form differentiated cell types of the metacestode, but also cells expressing genes specific for protoscoleces, adult worms, and oncospheres, including an ortholog of the schistosome praziquantel target, EmTRPMPZQ. Finally, we show that primary cell cultures contain a cell population expressing an ortholog of the tumour necrosis factor α receptor family and that mammalian TNFα accelerates the development of metacestode vesicles from germinative cells. Taken together, our analyses provide a robust and comprehensive characterization of the Echinococcus germinative cell transcriptome, demonstrate expanded differentiation capability of metacestode derived stem cells, and underscore the potential of primary germinative cell cultures to investigate developmental processes of the parasite. These data are relevant for studies into the role of Echinococcus stem cells in parasite development and will facilitate the design of anti-parasitic drugs that specifically act on the parasite germinative cell compartment.

Small Molecule Ligands of the BET-like Bromodomain, SmBRD3, Affect Schistosoma mansoni Survival, Oviposition, and Development

Schistosomiasis is a disease affecting >200 million people worldwide, but its treatment relies on a single agent, praziquantel. To investigate new avenues for schistosomiasis control, we have conducted the first systematic analysis of bromodomain-containing proteins (BCPs) in a causative species, Schistosoma mansoni. Having identified 29 putative bromodomains (BRDs) in 22 S. mansoni proteins, we selected SmBRD3, a tandem BRD-containing BCP that shows high similarity to the human bromodomain and extra terminal domain (BET) family, for further studies. Screening 697 small molecules identified the human BET BRD inhibitor I-BET726 as a ligand for SmBRD3. An X-ray crystal structure of I-BET726 bound to the second BRD of SmBRD3 [SmBRD3(2)] enabled rational design of a quinoline-based ligand (15) with an ITC Kd = 364 ± 26.3 nM for SmBRD3(2). The ethyl ester pro-drug of compound 15 (compound 22) shows substantial effects on sexually immature larval schistosomula, sexually mature adult worms, and snail-infective miracidia in ex vivo assays.

Lentiviral Transduction-based CRISPR/Cas9 Editing of Schistosoma mansoni Acetylcholinesterase

Background

Recent studies on CRISPR/Cas9-mediated gene editing in Schistosoma mansoni have shed new light on the study and control of this parasitic helminth. However, the gene editing efficiency in this parasite is modest.

Methods

To improve the efficiency of CRISPR/Cas9 genome editing in schistosomes, we used lentivirus, which has been effectively used for gene editing in mammalian cells, to deliver plasmid DNA encoding Cas9 nuclease, a sgRNA targeting acetylcholinesterase (SmAChE) and a mCherry fluorescence marker into schistosomes.

Results

MCherry fluorescence was observed in transduced eggs, schistosomula, and adult worms, indicating that the CRISPR components had been delivered into these parasite stages by lentivirus. In addition, clearly changed phenotypes were observed in SmAChE-edited parasites, including decreased SmAChE activity, reduced hatching ability of edited eggs, and altered behavior of miracidia hatched from edited eggs. Next-generation sequencing analysis demonstrated that the lentiviral transduction-based CRISPR/Cas9 gene modifications in SmAChE-edited schistosomes were homology-directed repair predominant but with much lower efficiency than that obtained using electroporation (data previously published by our laboratory) for the delivery of CRISPR components.

Conclusion

Taken together, electroporation is more efficient than lentiviral transduction in the delivery of CRISPR/Cas9 into schistosomes for programmed genome editing. The exploration of tactics for enhancing CRISPR/Cas9 gene editing provides the basis for the future improvement of programmed genome editing in S. mansoni.

The stage- and sex-specific transcriptome of the human parasite Schistosoma mansoni

The flatworm Schistosoma mansoni is an important but neglected pathogen that causes the disease schistosomiasis in millions of people worldwide. The parasite has a complex life cycle, undergoing sexual reproduction in a mammalian host and asexual replication in a snail host. Understanding the molecular mechanisms that the parasite uses to transition between hosts and develop into dimorphic reproductively competent adults may reveal new strategies for control. We present the first comprehensive transcriptomic analysis of S. mansoni, from eggs to sexually naïve worms. Focusing on eight life stages spanning free-living water-borne and parasitic stages from both intermediate and definitive hosts, we have generated deep RNA-seq data for five replicates per group for a total of 75 data sets. The data were produced using a single approach to increase the accuracy of stage-to-stage comparisons and made accessible via a user-friendly tool to visualise and explore gene expression ( https://lifecycle.schisto.xyz/ ). These data are valuable for understanding the biology and sex-specific development of schistosomes and the interpretation of complementary genomic and functional genetics studies.

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.

A report on the complete mitochondrial genome of the trematode Azygia robusta Odhner, 1911, its new definitive host from the Russian Far East, and unexpected phylogeny of Azygiidae within Digenea, as inferred from mitogenome sequences

New data on the complete mitochondrial genome of Azygia robusta (Azygiidae) were obtained by the next-generation sequencing (NGS) approach. The mitochondrial DNA (mtDNA) of A. robusta had a length of 13 857 bp and included 12 protein-coding genes, two ribosomal genes, 22 transfer RNA genes, and two non-coding regions. The nucleotide sequences of the complete mitochondrial genomes of two A. robusta specimens differed from each other by 0.12 ± 0.03%. Six of 12 protein-coding genes demonstrated intraspecific variation. The difference between the nucleotide sequences of the complete mitochondrial genomes of A. robusta and Azygia hwangtsiyui was 26.95 ± 0.35%; the interspecific variation of protein-coding genes between A. robusta and A. hwangtsiyui ranged from 20.5 ± 0.9% (cox1) to 30.7 ± 1.2% (nad5). The observed gene arrangement in the mtDNA sequence of A. robusta was identical to that of A. hwangtsiyui. Codon usage and amino acid frequencies were highly similar between A. robusta and A. hwangtsiyui. The results of phylogenetic analyses based on mtDNA protein-coding regions showed that A. robusta is closely related to A. hwangtsiyui (belonging to the same suborder, Azygiida) that formed a distinct early-diverging branch relative to all other Digenea. A preliminary morphological analysis of paratypes of the two azygiid specimens studied showed visible morphological differences between them. The specimen extracted from Sakhalin taimen (Parahucho perryi) was most similar to A. robusta. Thus, we here provide the first record of a new definitive host, P. perryi, for A. robusta and also molecular characteristics of the trematode specimens.

Targeted insertion and reporter transgene activity at a gene safe harbor of the human blood fluke, Schistosoma mansoni

The identification and characterization of genomic safe harbor sites (GSHs) can facilitate consistent transgene activity with minimal disruption to the host cell genome. We combined computational genome annotation and chromatin structure analysis to predict the location of four GSHs in the human blood fluke, Schistosoma mansoni, a major infectious pathogen of the tropics. A transgene was introduced via CRISPR-Cas-assisted homology-directed repair into one of the GSHs in the egg of the parasite. Gene editing efficiencies of 24% and transgene-encoded fluorescence of 75% of gene-edited schistosome eggs were observed. The approach advances functional genomics for schistosomes by providing a tractable path for generating transgenics using homology-directed, repair-catalyzed transgene insertion. We also suggest that this work will serve as a roadmap for the development of similar approaches in helminths more broadly.

What about the Cytoskeletal and Related Proteins of Tapeworms in the Host's Immune Response? An Integrative Overview

Recent advances have increased our understanding of the molecular machinery in the cytoskeleton of mammalian cells, in contrast to the case of tapeworm parasites, where cytoskeleton remains poorly characterized. The pertinence of a better knowledge of the tapeworm cytoskeleton is linked to the medical importance of these parasitic diseases in humans and animal stock. Moreover, its study could offer new possibilities for the development of more effective anti-parasitic drugs, as well as better strategies for their surveillance, prevention, and control. In the present review, we compile the results of recent experiments on the cytoskeleton of these parasites and analyze how these novel findings might trigger the development of new drugs or the redesign of those currently used in addition to supporting their use as biomarkers in cutting-edge diagnostic tests.

Molecular characterization of Smtdc-1 and Smddc-1 discloses roles as male-competence factors for the sexual maturation of Schistosoma mansoni females

Introduction

Schistosomes are the only mammalian flatworms that have evolved separate sexes. A key question of schistosome research is the male-dependent sexual maturation of the female since a constant pairing contact with a male is required for the onset of gonad development in the female. Although this phenomenon is long known, only recently a first peptide-based pheromone of males was identified that contributes to the control of female sexual development. Beyond this, our understanding of the molecular principles inducing the substantial developmental changes in a paired female is still rudimentary.

Objectives

Previous transcriptomic studies have consistently pointed to neuronal genes being differentially expressed and upregulated in paired males. These genes included Smp_135230 and Smp_171580, both annotated as aromatic-L-amino-acid decarboxylases (DOPA decarboxylases). Here, we characterized both genes and investigated their roles in male-female interaction of S. mansoni.

Methodologies/findings

Sequence analyses indicated that Smp_135230 represents an L-tyrosine decarboxylase (Smtdc-1), whereas Smp_171580 represents a DOPA decarboxylase (Smddc-1). By qRT-PCR, we confirmed the male-specific and pairing-dependent expression of both genes with a significant bias toward paired males. RNA-interference experiments showed a strong influence of each gene on gonad differentiation in paired females, which was enhanced by double knockdown. Accordingly, egg production was significantly reduced. By confocal laser scanning microscopy, a failure of oocyte maturation was found in paired knockdown females. Whole-mount in situ hybridization patterns exhibited the tissue-specific occurrence of both genes in particular cells at the ventral surface of the male, the gynecophoral canal, which represents the physical interface of both genders. These cells probably belong to the predicted neuronal cluster 2 of S. mansoni.

Conclusion

Our results suggest that Smtdc-1 and Smddc-2 are male-competence factors that are expressed in neuronal cells at the contact zone between the genders as a response of pairing to subsequently control processes of female sexual maturation.

Characterization of Schistosome Sox Genes and Identification of a Flatworm Class of Sox Regulators

Schistosome helminths infect over 200 million people across 78 countries and are responsible for nearly 300,000 deaths annually. However, our understanding of basic genetic pathways crucial for schistosome development is limited. The sex determining region Y-box 2 (Sox2) protein is a Sox B type transcriptional activator that is expressed prior to blastulation in mammals and is necessary for embryogenesis. Sox expression is associated with pluripotency and stem cells, neuronal differentiation, gut development, and cancer. Schistosomes express a Sox-like gene expressed in the schistosomula after infecting a mammalian host when schistosomes have about 900 cells. Here, we characterized and named this Sox-like gene SmSOXS1. SmSoxS1 protein is a developmentally regulated activator that localizes to the anterior and posterior ends of the schistosomula and binds to Sox-specific DNA elements. In addition to SmSoxS1, we have also identified an additional six Sox genes in schistosomes, two Sox B, one SoxC, and three Sox genes that may establish a flatworm-specific class of Sox genes with planarians. These data identify novel Sox genes in schistosomes to expand the potential functional roles for Sox2 and may provide interesting insights into early multicellular development of flatworms.

Using ChEMBL to Complement Schistosome Drug Discovery

Schistosomiasis is one of the most important neglected tropical diseases. Until an effective vaccine is registered for use, the cornerstone of schistosomiasis control remains chemotherapy with praziquantel. The sustainability of this strategy is at substantial risk due to the possibility of praziquantel insensitive/resistant schistosomes developing. Considerable time and effort could be saved in the schistosome drug discovery pipeline if available functional genomics, bioinformatics, cheminformatics and phenotypic resources are systematically leveraged. Our approach, described here, outlines how schistosome-specific resources/methodologies, coupled to the open-access drug discovery database ChEMBL, can be cooperatively used to accelerate early-stage, schistosome drug discovery efforts. Our process identified seven compounds (fimepinostat, trichostatin A, NVP-BEP800, luminespib, epoxomicin, CGP60474 and staurosporine) with ex vivo anti-schistosomula potencies in the sub-micromolar range. Three of those compounds (epoxomicin, CGP60474 and staurosporine) also demonstrated potent and fast-acting ex vivo effects on adult schistosomes and completely inhibited egg production. ChEMBL toxicity data were also leveraged to provide further support for progressing CGP60474 (as well as luminespib and TAE684) as a novel anti-schistosomal compound. As very few compounds are currently at the advanced stages of the anti-schistosomal pipeline, our approaches highlight a strategy by which new chemical matter can be identified and quickly progressed through preclinical development.

Transforming growth factor-β signalling regulates protoscolex formation in the Echinococcus multilocularis metacestode

The lethal zoonosis alveolar echinococcosis (AE) is caused by tumor-like, infiltrative growth of the metacestode larval stage of the tapeworm Echinococcus multilocularis. We previously showed that the metacestode is composed of posteriorized tissue and that the production of the subsequent larval stage, the protoscolex, depends on re-establishment of anterior identities within the metacestode germinative layer. It is, however, unclear so far how protoscolex differentiation in Echinococcus is regulated. We herein characterized the full complement of E. multilocularis TGFβ/BMP receptors, which is composed of one type II and three type I receptor serine/threonine kinases. Functional analyzes showed that all Echinococcus TGFβ/BMP receptors are enzymatically active and respond to host derived TGFβ/BMP ligands for activating downstream Smad transcription factors. In situ hybridization experiments demonstrated that the Echinococcus TGFβ/BMP receptors are mainly expressed by nerve and muscle cells within the germinative layer and in developing brood capsules. Interestingly, the production of brood capsules, which later give rise to protoscoleces, was strongly suppressed in the presence of inhibitors directed against TGFβ/BMP receptors, whereas protoscolex differentiation was accelerated in response to host BMP2 and TGFβ. Apart from being responsive to host TGFβ/BMP ligands, protoscolex production also correlated with the expression of a parasite-derived TGFβ-like ligand, EmACT, which is expressed in early brood capsules and which is strongly expressed in anterior domains during protoscolex development. Taken together, these data indicate an important role of TGFβ/BMP signalling in Echinococcus anterior pole formation and protoscolex development. Since TGFβ is accumulating around metacestode lesions at later stages of the infection, the host immune response could thus serve as a signal by which the parasite senses the time point at which protoscoleces must be produced. Overall, our data shed new light on molecular mechanisms of host-parasite interaction during AE and are relevant for the development of novel treatment strategies.

Proteases and their inhibitors involved in Schistosoma mansoni egg-host interaction revealed by comparative transcriptomics with Fasciola hepatica eggs

Schistosoma mansoni eggs are the main causative agents of the pathological manifestations of schistosomiasis. The eggs are laid in the host bloodstream, then they migrate through the intestinal wall into the lumen. However, a significant proportion of the eggs become lodged in the liver, where they cause inflammation and fibrosis. In this study, we focus on a specific group of proteins expressed by the egg, namely proteases and their inhibitors. These molecules are often involved in schistosome-host interactions, but are still unexplored in the egg stage. Using RNA-seq and comparative transcriptomics of immature and mature S. mansoni eggs, we mapped the portfolio of proteases and their inhibitors, and determined their gene expression levels. In addition, we compared these data with gene expression of proteases and their inhibitors in Fasciola hepatica eggs. Fasciola hepatica eggs served as a useful comparative model, as they do not migrate through tissues and inflict pathology. We detected transcription of 135 and 117 proteases in S. mansoni and F. hepatica eggs, respectively, with 87 identified as orthologous between the two species. In contrast, we observed only four orthologous inhibitors out of 21 and 16 identified in S. mansoni and F. hepatica eggs, respectively. Among others, we measured high and developmentally regulated levels of expression of metalloproteases in S. mansoni eggs, specifically aminopeptidase N1, endothelin-converting enzyme 1, and several leishmanolysin-like peptidases. We identified highly transcribed protease inhibitors serpin and alpha-2-macroglobulin that are unique to S. mansoni eggs, and antistasin-like inhibitor in F. hepatica eggs. This study provides new insights into the portfolio of proteases and inhibitors expressed by S. mansoni with potential roles in egg tissue migration, stimulation of angiogenesis, and interaction with host blood and immunity.

A cercarial invadolysin interferes with the host immune response and facilitates infection establishment of Schistosoma mansoni

Schistosoma mansoni employs immune evasion and immunosuppression to overcome immune responses mounted by its snail and human hosts. Myriad immunomodulating factors underlie this process, some of which are proteases. Here, we demonstrate that one protease, an invadolysin we have termed SmCI-1, is released from the acetabular glands of S. mansoni cercaria and is involved in creating an immunological milieu favorable for survival of the parasite. The presence of SmCI-1 in the cercarial stage of S. mansoni is released during transformation into the schistosomula. SmCI-1 functions as a metalloprotease with the capacity to cleave collagen type IV, gelatin and fibrinogen. Additionally, complement component C3b is cleaved by this protease, resulting in inhibition of the classical and alternative complement pathways. Using SmCI-1 knockdown cercariae, we demonstrate that SmCI-1 protects schistosomula from complement-mediated lysis in human plasma. We also assess the effect of SmCI-1 on cytokine release from human peripheral blood mononuclear cells, providing compelling evidence that SmCI-1 promotes an anti-inflammatory microenvironment by enhancing production of IL-10 and suppressing the production of inflammatory cytokines like IL-1B and IL-12p70 and those involved in eosinophil recruitment and activation, like Eotaxin-1 and IL-5. Finally, we utilize the SmCI-1 knockdown cercaria in a mouse model of infection, revealing a role for SmCI-1 in S. mansoni survival.

Trans-splicing in the cestode Hymenolepis microstoma is constitutive across the life cycle and depends on gene structure and composition

Spliced leader (SL) trans-splicing is a key process during mRNA maturation of many eukaryotes, in which a short sequence (SL) is transferred from a precursor SL-RNA into the 5' region of an immature mRNA. This mechanism is present in flatworms, in which it is known to participate in the resolution of polycistronic transcripts. However, most trans-spliced transcripts are not part of operons, and it is not clear if this process may participate in additional regulatory mechanisms in this group. In this work, we present a comprehensive analysis of SL trans-splicing in the model cestode Hymenolepis microstoma. We identified four different SL-RNAs which are indiscriminately trans-spliced to 622 gene models. SL trans-splicing is enriched in constitutively expressed genes and does not appear to be regulated throughout the life cycle. Operons represented at least 20% of all detected trans-spliced gene models, showed conservation to those of the cestode Echinococcus multilocularis, and included complex loci such as an alternative operon (processed as either a single gene through cis-splicing or as two genes of a polycistron). Most insertion sites were identified in the 5' untranslated region (UTR) of monocistronic genes. These genes frequently contained introns in the 5' UTR, in which trans-splicing used the same acceptor sites as cis-splicing. These results suggest that, unlike other eukaryotes, trans-splicing is associated with internal intronic promoters in the 5' UTR, resulting in transcripts with strong splicing acceptor sites without competing cis-donor sites, pointing towards a simple mechanism driving the evolution of novel SL insertion sites.

A major locus confers triclabendazole resistance in Fasciola hepatica and shows dominant inheritance

Fasciola hepatica infection is responsible for substantial economic losses in livestock worldwide and poses a threat to human health in endemic areas. The mainstay of control in livestock and the only drug licenced for use in humans is triclabendazole (TCBZ). TCBZ resistance has been reported on every continent and threatens effective control of fasciolosis in many parts of the world. To date, understanding the genetic mechanisms underlying TCBZ resistance has been limited to studies of candidate genes, based on assumptions of their role in drug action. Taking an alternative approach, we combined a genetic cross with whole-genome sequencing to localise a ~3.2Mbp locus within the 1.2Gbp F. hepatica genome that confers TCBZ resistance. We validated this locus independently using bulk segregant analysis of F. hepatica populations and showed that it is the target of drug selection in the field. We genotyped individual parasites and tracked segregation and reassortment of SNPs to show that TCBZ resistance exhibits Mendelian inheritance and is conferred by a dominant allele. We defined gene content within this locus to pinpoint genes involved in membrane transport, (e.g. ATP-binding cassette family B, ABCB1), transmembrane signalling and signal transduction (e.g. GTP-Ras-adenylyl cyclase and EGF-like protein), DNA/RNA binding and transcriptional regulation (e.g. SANT/Myb-like DNA-binding domain protein) and drug storage and sequestration (e.g. fatty acid binding protein, FABP) as prime candidates for conferring TCBZ resistance. This study constitutes the first experimental cross and genome-wide approach for any heritable trait in F. hepatica and is key to understanding the evolution of drug resistance in Fasciola spp. to inform deployment of efficacious anthelmintic treatments in the field.

CRISPR interference for sequence-specific regulation of fibroblast growth factor receptor A in Schistosoma mansoni

Employing the flatworm parasite Schistosoma mansoni as a model, we report the first application of CRISPR interference (CRISPRi) in parasitic helminths for loss-of-function studies targeting the SmfgfrA gene which encodes the stem cell marker, fibroblast growth factor receptor A (FGFRA). SmFGFRA is essential for maintaining schistosome stem cells and critical in the schistosome-host interplay. The SmfgfrA gene was targeted in S. mansoni adult worms, eggs and schistosomula using a catalytically dead Cas9 (dCas9) fused to a transcriptional repressor KRAB. We showed that SmfgfrA repression resulted in considerable phenotypic differences in the modulated parasites compared with controls, including reduced levels of SmfgfrA transcription and decreased protein expression of SmFGFRA, a decline in EdU (thymidine analog 5-ethynyl-2'-deoxyuridine, which specifically stains schistosome stem cells) signal, and an increase in cell apoptosis. Notably, reduced SmfgfrA transcription was evident in miracidia hatched from SmfgfrA-repressed eggs, and resulted in a significant change in miracidial behavior, indicative of a durable repression effect caused by CRISPRi. Intravenous injection of mice with SmfgfrA-repressed eggs resulted in granulomas that were markedly reduced in size and a decline in the level of serum IgE, emphasizing the importance of SmFGFRA in regulating the host immune response induced during schistosome infection. Our findings show the feasibility of applying CRISPRi for effective, targeted transcriptional repression in schistosomes, and provide the basis for employing CRISPRi to selectively perturb gene expression in parasitic helminths on a genome-wide scale.

Single-sex schistosomiasis: a mini review

Schistosomiasis is a neglected tropical disease caused by dioecious blood flukes of the genus Schistosoma and second to malaria as a parasitic disease with significant socio-economic impacts. Mating is essential for maturation of male and female schistosomes and for females to lay of eggs, which are responsible for the pathogenesis and propagation of the life cycle beyond the mammalian host. Single-sex schistosomes, which do not produce viable eggs without mating, have been overlooked given the symptomatic paucity of the single-sex schistosomiasis and limited diagnostic toolkit. Besides, single-sex schistosomes are less sensitive to praziquantel. Therefore, these issues should be considered to achieve the elimination of this infection disease. The aim of this review is to summarize current progress in research of single-sex schistosomes and host-parasite interactions.

Schistosoma mansoni phosphoglycerate mutase: a glycolytic ectoenzyme with thrombolytic potential

Schistosomiasis is a debilitating parasitic disease caused by intravascular flatworms called schistosomes (blood flukes) that affects >200 million people worldwide. Proteomic analysis has revealed the surprising presence of classical glycolytic enzymes – typically cytosolic proteins – located on the extracellular surface of the parasite tegument (skin). Immunolocalization experiments show that phosphoglycerate mutase (PGM) is widely expressed in parasite tissues and is highly expressed in the tegument. We demonstrate that live Schistosoma mansoni parasites express enzymatically active PGM on their tegumental surface. Suppression of PGM using RNA interference (RNAi) diminishes S. mansoni PGM (SmPGM) gene expression, protein levels, and surface enzyme activity. Sequence comparisons place SmPGM in the cofactor (2,3-bisphosphoglycerate)-dependent PGM (dPGM) family. We have produced recombinant SmPGM (rSmPGM) in an enzymatically active form in Escherichia coli. The Michaelis-Menten constant (K_m) of rSmPGM for its glycolytic substrate (3-phosphoglycerate) is 0.85 mM ± 0.02. rSmPGM activity is inhibited by the dPGM-specific inhibitor vanadate. Here, we show that rSmPGM not only binds to plasminogen but also promotes its conversion to an active form (plasmin) in vitro. This supports the hypothesis that host-interactive tegumental proteins (such as SmPGM), by enhancing plasmin formation, may help degrade blood clots around the worms in the vascular microenvironment and thus promote parasite survival in vivo.

Profiling the serum proteome during Schistosoma mansoni infection in the BALB/c mice: A focus on the altered lipid metabolism as a key modulator of host-parasite interactions

Schistosomiasis represents a condition in which every aspect of the disease, starting from skin invasion of the cercariae to egg laying by adult worms, incites a tissue response from the vertebrate host. This response, whether acute or chronic, leads to the appearance of reporter molecules of tissue injury in bodily fluids that could be surveyed as markers for disease diagnosis, status and prognosis. In this scenario, the serum proteome associated with a schistosome infection remains poorly explored; particularly by the use of high-throughput mass spectrometric instrumentation. In this study, we aimed to comparatively examine the serum proteome of control versus infected BALB/c mice, spanning the interval between the onset of egg laying and the peak of the acute phase of infection. Compositional analysis of the sera, using one dimensional reversed-phase fractionation of tryptic peptides coupled to mass spectrometry, allowed identification of 453 constituents. Among these, over 30% (143 molecules) were differentially present comparing sera from infected and non-infected mice, as revealed by quantitative label-free shotgun approach. The majority of proteins exhibiting altered levels was categorised as belonging to immune response (acute phase-related proteins) followed by those linked to lipid transport and metabolism. Inspection of the lipid profile from control and infected individuals demonstrated more pronounced and significant alterations in triglycerides, VLDL and HDL fractions (p<0,001), attesting for a disturbance in circulating lipid molecules, and suggesting a key role in host-parasite interactions. Our findings provide a global view of the serum proteome in the context of experimental schistosomiasis during the acute phase of infection. It contributes by listing key molecules that could be monitored to inform on the associated inflammatory disease status. We hope it will shed light into uncovered aspects of the Schistosoma mansoni parasitism in the vertebrate host, particularly those related to modulation of the lipid metabolism mediating immune responses.

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.

Chromosome-level genome assembly defines female-biased genes associated with sex determination and differentiation in the human blood fluke Schistosoma japonicum

Schistosomiasis is a neglected tropical disease of humans caused by blood flukes of the genus Schistosoma, the only dioecious parasitic flatworm. Although aspects of sex determination, differentiation and reproduction have been studied in some Schistosoma species, almost nothing is known for Schistosoma japonicum, the causative agent of schistosomiasis japonica. This mainly reflects the lack of high-quality genomic and transcriptomic resources for this species. As current genomes for S. japonicum are highly fragmented, we assembled and report a chromosome-level reference genome (seven autosomes, the Z-chromosome and partial W-chromosome), achieving a substantially enhanced gene annotation. Utilizing this genome, we discovered that the sex chromosomes of S. japonicum and its congener S. mansoni independently suppressed recombination during evolution, forming five and two evolutionary strata, respectively. By exploring the W-chromosome and sex-specific transcriptomes, we identified 35 W-linked genes and 257 female-preferentially transcribed genes (FTGs) from our chromosomal assembly and uncovered a signature for sex determination and differentiation in S. japonicum. These FTGs clustering within autosomes or the Z-chromosome exhibit a highly dynamic transcription profile during the pairing of female and male schistosomula, thereby representing a critical phase for the maturation of the female worms and suggesting distinct layers of regulatory control of gene transcription at this development stage. Collectively, these data provide a valuable resource for further functional genomic characterization of S. japonicum, shed light on the evolution of sex chromosomes in this highly virulent human blood fluke, and provide a pathway to identify novel targets for development of intervention tools against schistosomiasis.

Docking-Based Virtual Screening Enables Prioritizing Protein Kinase Inhibitors With In Vitro Phenotypic Activity Against Schistosoma mansoni

Schistosomiasis is a parasitic neglected disease with praziquantel (PZQ) utilized as the main drug for treatment, despite its low effectiveness against early stages of the worm. To aid in the search for new drugs to tackle schistosomiasis, computer-aided drug design has been proved a helpful tool to enhance the search and initial identification of schistosomicidal compounds, allowing fast and cost-efficient progress in drug discovery. The combination of high-throughput in silico data followed by in vitro phenotypic screening assays allows the assessment of a vast library of compounds with the potential to inhibit a single or even several biological targets in a more time- and cost-saving manner. Here, we describe the molecular docking for in silico screening of predicted homology models of five protein kinases (JNK, p38, ERK1, ERK2, and FES) of Schistosoma mansoni against approximately 85,000 molecules from the Managed Chemical Compounds Collection (MCCC) of the University of Nottingham (UK). We selected 169 molecules predicted to bind to SmERK1, SmERK2, SmFES, SmJNK, and/or Smp38 for in vitro screening assays using schistosomula and adult worms. In total, 89 (52.6%) molecules were considered active in at least one of the assays. This approach shows a much higher efficiency when compared to using only traditional high-throughput in vitro screening assays, where initial positive hits are retrieved from testing thousands of molecules. Additionally, when we focused on compound promiscuity over selectivity, we were able to efficiently detect active compounds that are predicted to target all kinases at the same time. This approach reinforces the concept of polypharmacology aiming for “one drug-multiple targets”. Moreover, at least 17 active compounds presented satisfactory drug-like properties score when compared to PZQ, which allows for optimization before further in vivo screening assays. In conclusion, our data support the use of computer-aided drug design methodologies in conjunction with high-throughput screening approach.

Schistosoma mansoni Fibroblast Growth Factor Receptor A Orchestrates Multiple Functions in Schistosome Biology and in the Host-Parasite Interplay

Stem cells play significant roles in driving the complex life cycle of Schistosoma mansoni. Fibroblast growth factor (FGF) receptor A (SmFGFRA) is essential for maintaining the integrity of schistosome stem cells. Using immunolocalization, we demonstrated that SmFGFRA was distributed abundantly in germinal/stem cells of different S. mansoni life stages including eggs, miracidia, cercariae, schistosomula and adult worms. Indeed, SmFGFRA was also localized amply in embryonic cells and in the perinuclear region of immature eggs; von Lichtenberg's layer and the neural mass of mature eggs; the ciliated surface and neural mass of miracidia; the tegument cytosol of cercariae, schistosomula and adult worms; and was present in abundance in the testis and vitellaria of adult worms of S. mansoni. The distribution pattern of SmFGFRA illustrates the importance of this molecule in maintaining stem cells, development of the nervous and reproductive system of schistosomes, and in the host-parasite interplay. We showed SmFGFRA can bind human FGFs, activating the mitogen activated protein kinase (MAPK) pathway of adult worms in vitro. Inhibition of FGF signaling by the specific tyrosine kinase inhibitor BIBF 1120 significantly reduced egg hatching ability and affected the behavior of miracidia hatched from the treated eggs, emphasizing the importance of FGF signaling in driving the life cycle of S. mansoni. Our findings provide increased understanding of the complex schistosome life cycle and host-parasite interactions, indicating components of the FGF signaling pathway may represent promising targets for developing new interventions against schistosomiasis.

Evolution of sexual systems, sex chromosomes and sex-linked gene transcription in flatworms and roundworms

Many species with separate male and female individuals (termed ‘gonochorism’ in animals) have sex-linked genome regions. Here, we investigate evolutionary changes when genome regions become completely sex-linked, by analyses of multiple species of flatworms (Platyhelminthes; among which schistosomes recently evolved gonochorism from ancestral hermaphroditism), and roundworms (Nematoda) which have undergone independent translocations of different autosomes. Although neither the evolution of gonochorism nor translocations fusing ancestrally autosomal regions to sex chromosomes causes inevitable loss of recombination, we document that formerly recombining regions show genomic signatures of recombination suppression in both taxa, and become strongly genetically degenerated, with a loss of most genes. Comparisons with hermaphroditic flatworm transcriptomes show masculinisation and some defeminisation in schistosome gonad gene expression. We also find evidence that evolution of sex-linkage in nematodes is accompanied by transcriptional changes and dosage compensation. Our analyses also identify sex-linked genes that could assist future research aimed at controlling some of these important parasites.

Transitions between hermaphroditic and separate sexes are relatively understudied in animals compared to pants. Here, Wang et al. reconstruct the evolution of separate sexes in the flatworms and complex changes of sex chromosomes in the roundworms.

Molecular fusion events in carcinogenic organisms: a bioinformatics study for the detection of fused proteins between viruses, bacteria and eukaryotes

Molecular fusion events have a prominent role in the initial steps of carcinogenesis. In this study, a bioinformatics analysis was performed between four organisms that are known to induce cancer development in humans: two viruses, Human Herpesvirus 4, and Human T-cell leukaemia virus, one bacterium, Helicobacter Pylori, and one trematode, Schistosoma mansoni. The annotated proteomes from these organisms were analysed using the SAFE software to identify protein fusion events, which may provide insight into protein function similarities and possible merging events during the course of evolution. Based on the results, five fused proteins with very similar functions were detected, whereas proteins with different functions that might act in the same molecular complex or biochemical pathway were not found. Thus, this study analysed the above four well-known cancer-related organisms with de novo bioinformatics programs and provided useful information on protein fusion events, hopefully leading to deeper understanding of carcinogenenesis.

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.

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.

RNA-Guided AsCas12a- and SpCas9-Catalyzed Knockout and Homology Directed Repair of the Omega-1 Locus of the Human Blood Fluke, Schistosoma mansoni

The efficiency of the RNA-guided AsCas12a nuclease of Acidaminococcus sp. was compared with SpCas9 from Streptococcus pyogenes, for functional genomics in Schistosoma mansoni. We deployed optimized conditions for the ratio of guide RNAs to the nuclease, donor templates, and electroporation parameters, to target a key schistosome enzyme termed omega-1. Programmed cleavages catalyzed by Cas12a and Cas9 resulted in staggered- and blunt-ended strand breaks, respectively. AsCas12a was more efficient than SpCas9 for gene knockout, as determined by TIDE analysis. CRISPResso2 analysis confirmed that most mutations were deletions. Knockout efficiency of both nucleases markedly increased in the presence of single-stranded oligodeoxynucleotide (ssODN) template. With AsCas12a, ssODNs representative of both the non-CRISPR target (NT) and target (T) strands were tested, resulting in KO efficiencies of 15.67, 28.71, and 21.43% in the SpCas9 plus ssODN, AsCas12a plus NT-ssODN, and AsCas12a plus T-ssODN groups, respectively. Trans-cleavage against the ssODNs by activated AsCas12a was not apparent in vitro. SpCas9 catalyzed more precise transgene insertion, with knock-in efficiencies of 17.07% for the KI_Cas9 group, 14.58% for KI_Cas12a-NT-ssODN, and 12.37% for KI_Cas12a-T-ssODN. Although AsCas12a induced fewer mutations per genome than SpCas9, the phenotypic impact on transcription and expression of omega-1 was similar for both nucleases.

Schistosoma mansoni α-N-acetylgalactosaminidase (SmNAGAL) regulates coordinated parasite movement and egg production

α-galactosidase (α-GAL) and α-N-acetylgalactosaminidase (α-NAGAL) are two glycosyl hydrolases responsible for maintaining cellular homeostasis by regulating glycan substrates on proteins and lipids. Mutations in the human genes encoding either enzyme lead to neurological and neuromuscular impairments seen in both Fabry- and Schindler/Kanzaki- diseases. Here, we investigate whether the parasitic blood fluke Schistosoma mansoni, responsible for the neglected tropical disease schistosomiasis, also contains functionally important α-GAL and α-NAGAL proteins. As infection, parasite maturation and host interactions are all governed by carefully-regulated glycosylation processes, inhibiting S. mansoni's α-GAL and α-NAGAL activities could lead to the development of novel chemotherapeutics. Sequence and phylogenetic analyses of putative α-GAL/α-NAGAL protein types showed Smp_089290 to be the only S. mansoni protein to contain the functional amino acid residues necessary for α-GAL/α-NAGAL substrate cleavage. Both α-GAL and α-NAGAL enzymatic activities were higher in females compared to males (p<0.05; α-NAGAL > α-GAL), which was consistent with smp_089290's female biased expression. Spatial localisation of smp_089290 revealed accumulation in parenchymal cells, neuronal cells, and the vitellaria and mature vitellocytes of the adult schistosome. siRNA-mediated knockdown (>90%) of smp_089290 in adult worms significantly inhibited α-NAGAL activity when compared to control worms (siLuc treated males, p<0.01; siLuc treated females, p<0.05). No significant reductions in α-GAL activities were observed in the same extracts. Despite this, decreases in α-NAGAL activities correlated with a significant inhibition in adult worm motility as well as in egg production. Programmed CRISPR/Cas9 editing of smp_089290 in adult worms confirmed the egg reduction phenotype. Based on these results, Smp_089290 was determined to act predominantly as an α-NAGAL (hereafter termed SmNAGAL) in schistosome parasites where it participates in coordinating movement and oviposition processes. Further characterisation of SmNAGAL and other functionally important glycosyl hydrolases may lead to the development of a novel anthelmintic class of compounds.

Comparative characterization of microRNAs of Schistosoma japonicum from SCID mice and BALB/c mice: Clues to the regulation of parasite growth and development

Schistosomiasis, caused by a parasite with a wide range of mammalian hosts, remains one of the most prevailing parasitic diseases in the world. While numerous studies have reported that the growth and reproduction of schistosomes in immunodeficient mice was significantly retarded, the underlying molecular mechanisms have yet to be revealed. In this study, we comparatively analyzed the microRNA expression of Schistosoma japonicum derived from SCID and BALB/c mice on the 35^th day post-infection by high-throughput RNA sequencing as prominent morphological abnormalities had been observed in schistosomes from SCID mice when compared with those from BALB/c mice. The results revealed that more than 72% and 61% of clean reads in the small RNA libraries of female and male schistosomes, respectively, could be mapped to the selected miRs in the miRBase or the sequences of species-specific genomes. Further analysis identified 122 miRNAs using TPM >0.01 as the threshold value, including 75 known and 47 novel miRNAs, 96 of which were commonly expressed across all the four tested schistosome libraries. Comparative analysis of the libraries of schistosomes from SCID and BALB/c mice identified 15 differentially expressed miRNAs (5 up-regulated and 10 down-regulated) among females and 16 among males (9 up-regulated and 7 down-regulated). Integrated analysis of the two sets of differentially expressed miRNAs of female and male worms identified 2 miRNAs (sja-miR-3488 and sja-miR-novel_29) that overlapped between female and male datasets. Prediction of miRNA targets and Gene Ontology (GO) term enrichment analysis of the predicted target genes revealed that these genes were involved in some important biological processes, such as nucleic acid metabolic process, macromolecule modification, and cellular aromatic compound metabolic process. The predicted target genes were further matched to the differentially expressed genes in male and female schistosomes from the above two hosts, obtaining 7 genes that may be responsible for regulating the growth, development and sex maturation of schistosomes. Taken together, this study provides the first identification of differentially expressed miRNAs in schistosomes from SCID and BALB/c mice. These miRNAs and their predicted target mRNAs are probably involved in the regulation of development, growth, and maturation of schistosomes. Therefore, this study expands our understanding of schistosome development regulation and host-parasite relationship, and also provides a valuable set of potential anti-schistosomal targets for prevention and control of schistosomiasis.

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

[Figure: see text].

Codon usage in the flatworm Schistosoma mansoni is shaped by the mutational bias towards A+T and translational selection, which increases GC-ending codons in highly expressed genes

Schistosoma mansoni is a trematode flatworm that parasitizes humans and produces a disease called bilharzia. At the genomic level, it is characterized by a low genomic GC content and an "isochore-like" structure, where GC-richest regions, mainly placed at the extremes of the chromosomes, are interspersed with low GC-regions. Furthermore, the GC-richest regions are at the same time the gene-richest, and where the most heavily expressed genes are placed. Taking these features into account, we decided to reanalyze the codon usage of this flatworm. Our results show that a) when all genes are considered together, the strong mutational bias towards A + T leads to a predominance of A/T-ending codons, b) a multivariate analysis discriminates between highly and lowly expressed genes, c) the sequences expressed at highest levels display a significant increase in G/C-ending codons, d) when comparing the molecular distances with a closely related species the synonymous distance in highly expressed genes is significantly lower than in lowly expressed sequences. Therefore, we conclude that despite previous results, which were performed with a small sample of genes, codon usage in S. mansoni is the result of two forces that operate in opposite directions: while mutational bias leads to a predominance of A/T codons, translational selection, working at the level of speed, increment G/C ending triplets.

Computational analysis of spliced leader trans-splicing in the regenerative flatworm Macrostomum lignano reveals its prevalence in conserved and stem cell related genes

In eukaryotes, trans-splicing is a process of nuclear pre-mRNA maturation where two different RNA molecules are joined together by the spliceosomal machinery utilizing mechanisms similar to cis-splicing. In diverse taxa of lower eukaryotes, spliced leader (SL) trans-splicing is the most frequent type of trans-splicing, when the same sequence derived from short small nuclear RNA molecules, called SL RNAs, is attached to the 5’ ends of different non-processed pre-mRNAs. One of the functions of SL trans-splicing is processing polycistronic pre-mRNA molecules transcribed from operons, when several genes are transcribed as one pre-mRNA molecule. However, only a fraction of trans-spliced genes reside in operons, suggesting that SL trans-splicing must also have some other, less understood functions. Regenerative flatworms are informative model organisms which hold the keys to understand the mechanism of stem cell regulation and specialization during regeneration and homeostasis. Their ability to regenerate is fueled by the division and differentiation of the adult somatic stem cell population called neoblasts. Macrostomum lignano is a flatworm model organism where substantial technological advances have been achieved in recent years, including the development of transgenesis. Although a large fraction of genes in M. lignano were estimated to be SL trans-spliced, SL trans-splicing was not studied in detail in M. lignano before. Here, we performed the first comprehensive study of SL trans-splicing in M. lignano. By reanalyzing the existing genome and transcriptome data of M. lignano, we estimate that 30 % of its genes are SL trans-spliced, 15 % are organized in operons, and almost 40 % are both SL trans-spliced and in operons. We annotated and characterized the sequence of SL RNA and characterized conserved cis- and SL transsplicing motifs. Finally, we found that a majority of SL trans-spliced genes are evolutionarily conserved and significantly over-represented in neoblast-specific genes. Our findings suggest an important role of SL trans-splicing in the regulation and maintenance of neoblasts in M. lignano.

The Schistosoma mansoni nuclear receptor FTZ-F1 maintains esophageal gland function via transcriptional regulation of meg-8.3

Schistosomes infect over 200 million of the world's poorest people, but unfortunately treatment relies on a single drug. Nuclear hormone receptors are ligand-activated transcription factors that regulate diverse processes in metazoans, yet few have been functionally characterized in schistosomes. During a systematic analysis of nuclear receptor function, we found that an FTZ-F1-like receptor was essential for parasite survival. Using a combination of transcriptional profiling and chromatin immunoprecipitation (ChIP), we discovered that the micro-exon gene meg-8.3 is a transcriptional target of SmFTZ-F1. We found that both Smftz-f1 and meg-8.3 are required for esophageal gland maintenance as well as integrity of the worm's head. Together, these studies define a new role for micro-exon gene function in the parasite and suggest that factors associated with the esophageal gland could represent viable therapeutic targets.

Genetic architecture of transmission stage production and virulence in schistosome parasites

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

Temperature-induced reorganisation of Schistocephalus solidus (Cestoda) proteome during the transition to the warm-blooded host

The protein composition of the cestode Schistocephalus solidus was measured in an experiment simulating the trophic transmission of the parasite from a cold-blooded to a warm-blooded host. The first hour of host colonisation was studied in a model experiment, in which sticklebacks Gasterosteus aculeatus infected with S. solidus were heated at 40°C for 1 h. As a result, a decrease in the content of one tegument protein was detected in the plerocercoids of S. solidus. Sexual maturation of the parasites was initiated in an experiment where S. solidus larvae were taken from fish and cultured in vitro at 40°C for 48 h. Temperature-independent changes in the parasite proteome were investigated by incubating plerocercoids at 22°C for 48 h in culture medium. Analysis of the proteome allowed us to distinguish the temperature-induced genes of S. solidus, as well as to specify the molecular markers of the plerocercoid and adult worms. The main conclusion of the study is that the key enzymes of long-term metabolic changes (glycogen consumption, protein production, etc.) in parasites during colonisation of a warm-blooded host are induced by temperature.

Schistosoma mansoni Larval Extracellular Vesicle protein 1 (SmLEV1) is an immunogenic antigen found in EVs released from pre-acetabular glands of invading cercariae

Extracellular Vesicles (EVs) are an integral component of cellular/organismal communication and have been found in the excreted/secreted (ES) products of both protozoan and metazoan parasites. Within the blood fluke schistosomes, EVs have been isolated from egg, schistosomula, and adult lifecycle stages. However, the role(s) that EVs have in shaping aspects of parasite biology and/or manipulating host interactions is poorly defined. Herein, we characterise the most abundant EV-enriched protein in Schistosoma mansoni tissue-migrating schistosomula (Schistosoma mansoni Larval Extracellular Vesicle protein 1 (SmLEV1)). Comparative sequence analysis demonstrates that lev1 orthologs are found in all published Schistosoma genomes, yet homologs are not found outside of the Schistosomatidae. Lifecycle expression analyses collectively reveal that smlev1 transcription peaks in cercariae, is male biased in adults, and is processed by alternative splicing in intra-mammalian lifecycle stages. Immunohistochemistry of cercariae using a polyclonal anti-recombinant SmLEV1 antiserum localises this protein to the pre-acetabular gland, with some disperse localisation to the surface of the parasite. S. mansoni-infected Ugandan fishermen exhibit a strong IgG1 response against SmLEV1 (dropping significantly after praziquantel treatment), with 11% of the cohort exhibiting an IgE response and minimal levels of detectable antigen-specific IgG4. Furthermore, mice vaccinated with rSmLEV1 show a slightly reduced parasite burden upon challenge infection and significantly reduced granuloma volumes, compared with control animals. Collectively, these results describe SmLEV1 as a Schistosomatidae-specific, EV-enriched immunogen. Further investigations are now necessary to uncover the full extent of SmLEV1's role in shaping schistosome EV function and definitive host relationships.

Crystallographic approach to fragment-based hit discovery against Schistosoma mansoni purine nucleoside phosphorylase

Several Schistosoma species cause Schistosomiasis, an endemic disease in 78 countries that is ranked second amongst the parasitic diseases in terms of its socioeconomic impact and human health importance. The drug recommended for treatment by the WHO is praziquantel (PZQ), but there are concerns associated with PZQ, such as the lack of information about its exact mechanism of action, its high price, its effectiveness - which is limited to the parasite's adult form - and reports of resistance. The parasites lack the de novo purine pathway, rendering them dependent on the purine salvage pathway or host purine bases for nucleotide synthesis. Thus, the Schistosoma purine salvage pathway is an attractive target for the development of necessary and selective new drugs. In this study, the purine nucleotide phosphorylase II (PNP2), a new isoform of PNP1, was submitted to a high-throughput fragment-based hit discovery using a crystallographic screening strategy. PNP2 was crystallized and crystals were soaked with 827 fragments, a subset of the Maybridge 1000 library. X-ray diffraction data was collected and structures were solved. Out of 827-screened fragments we have obtained a total of 19 fragments that show binding to PNP2. Fourteen of these fragments bind to the active site of PNP2, while five were observed in three other sites. Here we present the first fragment screening against PNP2.

The function of twister ribozyme variants in non-LTR retrotransposition in Schistosoma mansoni

The twister ribozyme is widely distributed over numerous organisms and is especially abundant in Schistosoma mansoni, but has no confirmed biological function. Of the 17 non-LTR retrotransposons known in S. mansoni, none have thus far been associated with ribozymes. Here we report the identification of novel twister variant (T-variant) ribozymes and their function in S. mansoni non-LTR retrotransposition. We show that T-variant ribozymes are located at the 5′ end of Perere-3 non-LTR retrotransposons in the S. mansoni genome. T-variant ribozymes were demonstrated to be catalytically active in vitro. In reporter constructs, T-variants were shown to cleave in vivo, and cleavage of T-variants was sufficient for the translation of downstream reporter genes. Our analysis shows that the T-variants and Perere-3 are transcribed together. Target site duplications (TSDs); markers of target-primed reverse transcription (TPRT) and footmarks of retrotransposition, are located adjacent to the T-variant cleavage site and suggest that T-variant cleavage has taken place inS. mansoni. Sequence heterogeneity in the TSDs indicates that Perere-3 retrotransposition is not site-specific. The TSD sequences contribute to the 5′ end of the terminal ribozyme helix (P1 stem). Based on these results we conclude that T-variants have a functional role in Perere-3 retrotransposition.

Satellite-Like W-Elements: Repetitive, Transcribed, and Putative Mobile Genetic Factors with Potential Roles for Biology and Evolution of Schistosoma mansoni

A large portion of animal and plant genomes consists of noncoding DNA. This part includes tandemly repeated sequences and gained attention because it offers exciting insights into genome biology. We investigated satellite-DNA elements of the platyhelminth Schistosoma mansoni, a parasite with remarkable biological features. Schistosoma mansoni lives in the vasculature of humans causing schistosomiasis, a disease of worldwide importance. Schistosomes are the only trematodes that have evolved separate sexes, and the sexual maturation of the female depends on constant pairing with the male. The schistosome karyotype comprises eight chromosome pairs, males are homogametic (ZZ) and females are heterogametic (ZW). Part of the repetitive DNA of S. mansoni are W-elements (WEs), originally discovered as female-specific satellite DNAs in the heterochromatic block of the W-chromosome. Based on new genome and transcriptome data, we performed a reanalysis of the W-element families (WEFs). Besides a new classification of 19 WEFs, we provide first evidence for stage-, sex-, pairing-, gonad-, and strain-specific/preferential transcription of WEs as well as their mobile nature, deduced from autosomal copies of full-length and partial WEs. Structural analyses suggested roles as sources of noncoding RNA-like hammerhead ribozymes, for which we obtained functional evidence. Finally, the variable WEF occurrence in different schistosome species revealed remarkable divergence. From these results, we propose that WEs potentially exert enduring influence on the biology of S. mansoni. Their variable occurrence in different strains, isolates, and species suggests that schistosome WEs may represent genetic factors taking effect on variability and evolution of the family Schistosomatidae.

Homology modelling, molecular dynamics simulation and docking evaluation of β-tubulin of Schistosoma mansoni

Schistosomiasis is one of the neglected diseases causing considerable morbidity and mortality throughout the world. Microtubules with its main component, tubulin play a vital role in helminthes including schistosomes. Benzimidazoles represent potential drug candidates by binding β-tubulin. The study aimed to generate a homology model for the β-tubulin of S. mansoni using the crystal structure of O visaries (Sheep) β-tubulin (PDB ID: 3N2G D) as a template, then different β-tubulin models were generated and two previously reported benzimidazole derivatives (NBTP-F and NBTP-OH) were docked to the generated models, the binding results indicated that both S. mansoni, S. haematobium were susceptible to the two NBTP derivatives. Additionally, three mutated versions of S. mansoni β-tubulin wild-type were generated and the mutation (F185Y) seems to slightly enhance the ligand binding. Dynamics simulation experiments showed S. haematobium β-tubulin is highly susceptible to the tested compounds; similar to S. mansoni, moreover, mutated models of S. mansoni β-tubulin altered its NBTPs susceptibility. Moreover, additional seven new benzimidazole derivatives were synthesized and tested by molecular docking on the generated model binding site of S. mansoni β-tubulin and were found to have good interaction inside the pocket.

Assessment of reference genes at six different developmental stages of Schistosoma mansoni for quantitative RT-PCR

Reverse-transcription quantitative real-time polymerase chain reaction (RT-qPCR) is the most used, fast, and reproducible method to confirm large-scale gene expression data. The use of stable reference genes for the normalization of RT-qPCR assays is recognized worldwide. No systematic study for selecting appropriate reference genes for usage in RT-qPCR experiments comparing gene expression levels at different Schistosoma mansoni life-cycle stages has been performed. Most studies rely on genes commonly used in other organisms, such as actin, tubulin, and GAPDH. Therefore, the present study focused on identifying reference genes suitable for RT-qPCR assays across six S. mansoni developmental stages. The expression levels of 25 novel candidates that we selected based on the analysis of public RNA-Seq datasets, along with eight commonly used reference genes, were systematically tested by RT-qPCR across six developmental stages of S. mansoni (eggs, miracidia, cercariae, schistosomula, adult males and adult females). The stability of genes was evaluated with geNorm, NormFinder and RefFinder algorithms. The least stable candidate reference genes tested were actin, tubulin and GAPDH. The two most stable reference genes suitable for RT-qPCR normalization were Smp_101310 (Histone H4 transcription factor) and Smp_196510 (Ubiquitin recognition factor in ER-associated degradation protein 1). Performance of these two genes as normalizers was successfully evaluated with females maintained unpaired or paired to males in culture for 8 days, or with worm pairs exposed for 16 days to double-stranded RNAs to silence a protein-coding gene. This study provides reliable reference genes for RT-qPCR analysis using samples from six different S. mansoni life-cycle stages.

The Search for a Schistosomiasis Vaccine: Australia's Contribution

Schistosomiasis, a neglected tropical disease caused by parasitic flatworms of the genus Schistosoma, results in considerable human morbidity in sub-Saharan Africa, in particular, but also parts of the Middle East, South America, and Southeast Asia. The anti-schistosome drug praziquantel is efficacious and safe against the adult parasites of all Schistosoma species infecting humans; however, it does not prevent reinfection and the development of drug resistance is a constant concern. The need to develop an effective vaccine is of great importance if the health of many in the developing world is to be improved. Indeed, vaccination, in combination with other public health measures, can provide an invaluable tool to achieve lasting control, leading to schistosomiasis elimination. Australia has played a leading role in schistosomiasis vaccine research over many years and this review presents an overview of some of the significant contributions made by Australian scientists in this important area.

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.

A Comparative Proteomic Analysis of Praziquantel-Susceptible and Praziquantel-Resistant Schistosoma mansoni Reveals Distinct Response Between Male and Female Animals

Schistosomiasis is a chronic neglected tropical disease saddling millions of people in the world, mainly children living in poor rural areas. Praziquantel (PZQ) is currently the only drug used for the treatment and control of this disease. However, the extensive use of this drug has brought concern about the emergence of PZQ-resistance/tolerance by Schistosoma mansoni. Studies of Schistosoma spp. genome, transcriptome, and proteome are crucial to better understand this situation. In this in vitro study, we compare the proteomes of a S. mansoni variant strain stably resistant to PZQ and isogenic to its fully susceptible parental counterpart, identifying proteins from male and female adult parasites of PZQ-resistant and PZQ-susceptible strains, exposed and not exposed to PZQ. A total of 60 Schistosoma spp. proteins were identified, some of which present or absent in either strain, which may putatively be involved in the PZQ-resistance phenomenon. These proteins were present in adult parasites not exposed to PZQ, but some of them disappeared when these adult parasites were exposed to the drug. Understanding the development of PZQ-resistance in S. mansoni is crucial to prolong the efficacy of the current drug and develop markers for monitoring the potential emergence of drug resistance.

Multiomic analysis of Schistosoma mansoni reveals unique expression profiles in cercarial heads and tails

Schistosomes require both molluscan and mammalian hosts for development. The larval cercaria exits the snail host and swims to identify and invade the mammalian host. The cercaria has two macrostructures, the head and the tail. The head invades the host, where it matures into an adult worm. The tail is lost after host invasion. Translation in the cercaria differs in each macrostructure, with higher levels of translation in the cercarial tail and little to no translational activity in the cercarial head. We compared the transcriptome and proteome of the cercarial head and tail and observed stark differences between the two macrostructures. We identified unique and differentially expressed transcripts and proteins, including ribosomal components expressed in higher levels in tails than in heads, which may explain the differences in translation levels between heads and tails. We also characterized the weak correlation between transcription and translation in infectious cercarial heads and tails.