Stem cell heterogeneity drives the parasitic life cycle of Schistosoma mansoni

Anterior-posterior polarity signals differentially regulate regeneration-competence of the tapeworm Hymenolepis diminuta

Competence to regenerate lost tissues varies widely across species. The rat tapeworm, Hymenolepis diminuta, undergoes continual cycles of shedding and regenerating thousands of reproductive segments to propagate the species. Despite its prowess, H. diminuta can only regenerate posteriorly from a singular tissue: the neck or germinative region (GR). What cells and signaling pathways restrict regeneration competence to the GR? In this study, we show that the head regulates regeneration-competence by promoting maintenance of the GR and inhibiting proglottid formation in a distance-dependent manner. Anterior-posterior (A-P) patterning within the GR provide local signals that mediate these head-dependent responses. βcat1 is necessary for stem cell maintenance, proliferation and proglottidization. On the other hand, sfrp is necessary for maintaining the GR at its proper length. Our study demonstrates that the head organizes a balance of pro- and anti-regeneration signals that must be integrated together and therefore control competence to regenerate.

RNAi screening of uncharacterized genes identifies promising druggable targets in Schistosoma japonicum

Schistosomiasis affects more than 250 million people worldwide and is one of the neglected tropical diseases. Currently, the treatment of schistosomiasis relies on a single drug-praziquantel-which has led to increasing pressure from drug resistance. Therefore, there is an urgent need to find new treatments. The development of genome sequencing has provided valuable information for understanding the biology of schistosomes. In the genome of Schistosoma japonicum, approximately 11% of the protein-coding sequences are uncharacterized genes (UGs) annotated as "hypothetical protein" or "protein of unknown function." These poorly understood genes have been unjustifiably neglected, although some may be essential for the survival of the parasites and serve as potential drug targets. In this study, we systematically mined the highly expressed UGs in both genders of this parasite throughout key developmental stages in their mammalian host, using our previously published S. japonicum genome and RNA-seq data. By employing in vitro RNA interference (RNAi), we screened 126 UGs that lack homologs in Homo sapiens and identified 8 that are essential for the parasite vitality. We further investigated two UGs, Sjc_0002003 and Sjc_0009272, which resulted in the most severe phenotypes. Fluorescence in situ hybridization demonstrated that both genes were expressed throughout the body without sex bias. Silencing either Sjc_0002003 or Sjc_0009272 reduced the cell proliferation in the body. Furthermore, in vivo RNAi indicated both genes are required for the growth and survival of the parasites in the mammalian host. For Sjc_0002003, we further characterize the underlying molecular cause of the observed phenotype. Through RNA-seq analysis and functional studies, we revealed that silencing Sjc_0002003 reduces the expression of a series of intestinal genes, including Sjc_0007312 (hypothetical protein), Sjc_0008276 (vha-17), Sjc_0002942 (PLA2G15), and Sjc_0003646 (SJCHGC09134 protein), leading to gut dilation. Our work highlights the importance of UGs in schistosomes as promising targets for drug development in the treatment of the schistosomiasis.

Application of microphysiological systems to unravel the mechanisms of schistosomiasis egg extravasation

Despite decades of control efforts, the prevalence of schistosomiasis remains high in many endemic regions, posing significant challenges to global health. One of the key factors contributing to the persistence of the disease is the complex life cycle of the Schistosoma parasite, the causative agent, which involves multiple stages of development and intricate interactions with its mammalian hosts and snails. Among the various stages of the parasite lifecycle, the deposition of eggs and their migration through host tissues is significant, as they initiate the onset of the disease pathology by inducing inflammatory reactions and tissue damage. However, our understanding of the mechanisms underlying Schistosoma egg extravasation remains limited, hindering efforts to develop effective interventions. Microphysiological systems, particularly organ-on-a-chip systems, offer a promising approach to study this phenomenon in a controlled experimental setting because they allow the replication of physiological microenvironments in vitro. This review provides an overview of schistosomiasis, introduces the concept of organ-on-a-chip technology, and discusses its potential applications in the field of schistosomiasis research.

The retinoic acid family-like nuclear receptor SmRAR identified by single-cell transcriptomics of ovarian cells controls oocyte differentiation in Schistosoma mansoni

Studies on transcription regulation in platyhelminth development are scarce, especially for parasitic flatworms. Here, we employed single-cell transcriptomics to identify genes involved in reproductive development in the trematode model Schistosoma mansoni. This parasite causes schistosomiasis, a major neglected infectious disease affecting >240 million people worldwide. The pathology of schistosomiasis is closely associated with schistosome eggs deposited in host organs including the liver. Unlike other trematodes, schistosomes exhibit distinct sexes, with egg production reliant on the pairing-dependent maturation of female reproductive organs. Despite this significance, the molecular mechanisms underlying ovary development and oocyte differentiation remain largely unexplored. Utilizing an organ isolation approach for S. mansoni, we extracted ovaries of paired females followed by single-cell RNA sequencing (RNA-seq) with disassociated oocytes. A total of 1967 oocytes expressing 7872 genes passed quality control (QC) filtering. Unsupervised clustering revealed four distinct cell clusters: somatic, germ cells and progeny, intermediate and late germ cells. Among distinct marker genes for each cluster, we identified a hitherto uncharacterized transcription factor of the retinoic acid receptor family, SmRAR. Functional analyses of SmRAR and associated genes like Smmeiob (meiosis-specific, oligonucleotide/oligosaccharide binding motif (OB) domain-containing) demonstrated their pairing-dependent and ovary-preferential expression and their decisive roles in oocyte differentiation of S. mansoni.

Regeneration in the absence of canonical neoblasts in an early branching flatworm

The remarkable regenerative abilities of flatworms are closely linked to neoblasts - adult pluripotent stem cells that are the only division-competent cell type outside of the reproductive system. Although the presence of neoblast-like cells and whole-body regeneration in other animals has led to the idea that these features may represent the ancestral metazoan state, the evolutionary origin of both remains unclear. Here we show that the catenulid Stenostomum brevipharyngium, a member of the earliest-branching flatworm lineage, lacks conventional neoblasts despite being capable of whole-body regeneration and asexual reproduction. Using a combination of single-nuclei transcriptomics, in situ gene expression analysis, and functional experiments, we find that cell divisions are not restricted to a single cell type and are associated with multiple fully differentiated somatic tissues. Furthermore, the cohort of germline multipotency genes, which are considered canonical neoblast markers, are not expressed in dividing cells, but in the germline instead, and we experimentally show that they are neither necessary for proliferation nor regeneration. Overall, our results challenge the notion that canonical neoblasts are necessary for flatworm regeneration and open up the possibility that neoblast-like cells may have evolved convergently in different animals, independent of their regenerative capacity.

Schistosome esophageal gland factor MEG-8.2 drives host cell lysis and interacts with host immune proteins

Schistosomes are blood flukes that ingest large amounts of host blood during their intra-mammalian stage. The ingested blood contains leukocytes that can be harmful, yet the parasites survive inside the host for decades, reflecting superb immune evasion mechanisms that remain poorly understood. Our previous work discovered that FoxA, a forkhead transcription factor, drives the production of the esophageal gland, an anterior digestive organ essential for degrading the ingested leukocytes and for in vivo survival. However, a comprehensive molecular makeup of the esophageal gland remains unclear. Importantly, which of the esophageal gland factors are responsible for degrading the ingested leukocytes, their mechanism of action, and how such a function relates to parasite survival and immune evasion remains unknown. Here, we identify additional esophageal gland genes by taking a comparative transcriptomics approach to identify transcripts altered in foxA knockdown adult schistosomes. A targeted RNAi screen coupled with biochemistry reveals that specific domains of the micro-exon gene MEG-8.2, can drive host cell lysis in a concentration-dependent manner. Using pull-down assays coupled with mass spectrometry, we discover that MEG-8.2 interacts with several host membrane and extracellular proteins that play important roles in activating innate and/or adaptive immunity. Together, our findings suggest a dual role for MEG-8.2 in effectively lysing the ingested cells in the esophageal lumen and interacting with specific host proteins to neutralize or suppress the host immunity. These findings lay an important foundation for exploiting esophageal gland factors to treat schistosomiasis.

Flexible use of conserved motif vocabularies constrains genome access in cell type evolution

Cell types evolve into a hierarchy with related types grouped into families. How cell type diversification is constrained by the stable separation between families over vast evolutionary times remains unknown. Here, integrating single-nucleus multiomic sequencing and deep learning, we show that hundreds of sequence features (motifs) divide into distinct sets associated with accessible genomes of specific cell type families. This division is conserved across highly divergent, early-branching animals including flatworms and cnidarians. While specific interactions between motifs delineate cell type relationships within families, surprisingly, these interactions are not conserved between species. Consistently, while deep learning models trained on one species can predict accessibility of other species' sequences, their predictions frequently rely on distinct, but synonymous, motif combinations. We propose that long-term stability of cell type families is maintained through genome access specified by conserved motif sets, or 'vocabularies', whereas cell types diversify through flexible use of motifs within each set.

A single-cell atlas of the miracidium larva of Schistosoma mansoni reveals cell types, developmental pathways, and tissue architecture

Schistosoma mansoni is a parasitic flatworm that causes the major neglected tropical disease schistosomiasis. The miracidium is the first larval stage of the life cycle. It swims and infects a freshwater snail, transforms into a mother sporocyst, where its stem cells generate daughter sporocysts that give rise to human-infective cercariae larvae. To understand the miracidium at cellular and molecular levels, we created a whole-body atlas of its ~365 cells. Single-cell RNA sequencing identified 19 transcriptionally distinct cell clusters. In situ hybridisation of tissue-specific genes revealed that 93% of the cells in the larva are somatic (57% neural, 19% muscle, 13% epidermal or tegument, 2% parenchyma, and 2% protonephridia) and 7% are stem. Whereas neurons represent the most diverse somatic cell types, trajectory analysis of the two main stem cell populations indicates that one of them is the origin of the tegument lineage and the other likely contains pluripotent cells. Furthermore, unlike the somatic cells, each of these stem populations shows sex-biased transcriptional signatures suggesting a cell-type-specific gene dosage compensation for sex chromosome-linked loci. The miracidium represents a simple developmental stage with which to gain a fundamental understanding of the molecular biology and spatial architecture of schistosome cells.

Regeneration in the absence of canonical neoblasts in an early branching flatworm

The remarkable regenerative abilities of flatworms are closely linked to neoblasts - adult pluripotent stem cells that are the only division-competent cell type outside of the reproductive system. Although the presence of neoblast-like cells and whole-body regeneration in other animals has led to the idea that these features may represent the ancestral metazoan state, the evolutionary origin of both remains unclear. Here we show that the catenulid Stenostomum brevipharyngium, a member of the earliest-branching flatworm lineage, lacks conventional neoblasts despite being capable of whole-body regeneration and asexual reproduction. Using a combination of single-nuclei transcriptomics, in situ gene expression analysis, and functional experiments, we find that cell divisions are not restricted to a single cell type and are associated with multiple fully differentiated somatic tissues. Furthermore, the cohort of germline multipotency genes, which are considered canonical neoblast markers, are not expressed in dividing cells, but in the germline instead, and we experimentally show that they are neither necessary for proliferation nor regeneration. Overall, our results challenge the notion that canonical neoblasts are necessary for flatworm regeneration and open up the possibility that neoblast-like cells may have evolved convergently in different animals, independent of their regenerative capacity.

Neoblast-like stem cells of Fasciola hepatica

The common liver fluke (Fasciola hepatica) causes the disease fasciolosis, which results in considerable losses within the global agri-food industry. There is a shortfall in the drugs that are effective against both the adult and juvenile life stages within the mammalian host, such that new drug targets are needed. Over the last decade the stem cells of parasitic flatworms have emerged as reservoirs of putative novel targets due to their role in development and homeostasis, including at host-parasite interfaces. Here, we investigate and characterise the proliferating cells that underpin development in F. hepatica. We provide evidence that these cells are capable of self-renewal, differentiation, and are sensitive to ionising radiation- all attributes of neoblasts in other flatworms. Changes in cell proliferation were also noted during the early stages of in vitro juvenile growth/development (around four to seven days post excystment), which coincided with a marked reduction in the nuclear area of proliferating cells. Furthermore, we generated transcriptomes from worms following irradiation-based ablation of neoblasts, identifying 124 significantly downregulated transcripts, including known stem cell markers such as fgfrA and plk1. Sixty-eight of these had homologues associated with neoblast-like cells in Schistosoma mansoni. Finally, RNA interference mediated knockdown of histone h2b (a marker of proliferating cells), ablated neoblast-like cells and impaired worm development in vitro. In summary, this work demonstrates that the proliferating cells of F. hepatica are equivalent to neoblasts of other flatworm species and demonstrate that they may serve as attractive targets for novel anthelmintics.

Intestinal cell diversity and treatment responses in a parasitic nematode at single cell resolution

BACKGROUND

Parasitic nematodes, significant pathogens for humans, animals, and plants, depend on diverse organ systems for intra-host survival. Understanding the cellular diversity and molecular variations underlying these functions holds promise for developing novel therapeutics, with specific emphasis on the neuromuscular system's functional diversity. The nematode intestine, crucial for anthelmintic therapies, exhibits diverse cellular phenotypes, and unraveling this diversity at the single-cell level is essential for advancing knowledge in anthelmintic research across various organ systems.

RESULTS

Here, using novel single-cell transcriptomics datasets, we delineate cellular diversity within the intestine of adult female Ascaris suum, a parasitic nematode species that infects animals and people. Gene transcripts expressed in individual nuclei of untreated intestinal cells resolved three phenotypic clusters, while lower stringency resolved additional subclusters and more potential diversity. Clusters 1 and 3 phenotypes displayed variable congruence with scRNA phenotypes of C. elegans intestinal cells, whereas the A. suum cluster 2 phenotype was markedly unique. Distinct functional pathway enrichment characterized each A. suum intestinal cell cluster. Cluster 2 was distinctly enriched for Clade III-associated genes, suggesting it evolved within clade III nematodes. Clusters also demonstrated differential transcriptional responsiveness to nematode intestinal toxic treatments, with Cluster 2 displaying the least responses to short-term intra-pseudocoelomic nematode intestinal toxin treatments.

CONCLUSIONS

This investigation presents advances in knowledge related to biological differences among major cell populations of adult A. suum intestinal cells. For the first time, diverse nematode intestinal cell populations were characterized, and associated biological markers of these cells were identified to support tracking of constituent cells under experimental conditions. These advances will promote better understanding of this and other parasitic nematodes of global importance, and will help to guide future anthelmintic treatments.

Single-cell transcriptomics of the human parasite Schistosoma mansoni first intra-molluscan stage reveals tentative tegumental and stem-cell regulators

Schistosomiasis is a major Neglected Tropical Disease, caused by the infection with blood flukes in the genus Schistosoma. To complete the life cycle, the parasite undergoes asexual and sexual reproduction within an intermediate snail host and a definitive mammalian host, respectively. The intra-molluscan phase provides a critical amplification step that ensures a successful transmission. However, the cellular and molecular mechanisms underlying the development of the intra-molluscan stages remain poorly understood. Here, single cell suspensions from S. mansoni mother sporocysts were produced and sequenced using the droplet-based 10X Genomics Chromium platform. Six cell clusters comprising two tegument, muscle, neuron, parenchyma and stem/germinal cell clusters were identified and validated by in situ hybridisation. Gene Ontology term analysis predicted key biological processes for each of the clusters, including three stem/germinal sub-clusters. Furthermore, putative transcription factors predicted for stem/germinal and tegument clusters may play key roles during parasite development and interaction with the intermediate host.

Purification and transcriptomic characterization of proliferative cells of Mesocestoides corti selectively affected by irradiation

Flatworms depend on stem cells for continued tissue growth and renewal during their life cycles, making these cells valuable drug targets. While neoblasts are extensively characterized in the free-living planarian Schmidtea mediterranea, and similar stem cells have been characterized in the trematode Schistosoma mansoni, their identification and characterization in cestodes is just emerging. Since stem cells are generally affected by irradiation, in this work we used this experimental approach to study the stem cells of the model cestode Mesocestoides corti. We found that gamma irradiation produces a dose-dependent decrease in proliferative cells, requiring higher doses than in other flatworms to completely abolish proliferation. The treatment results in the downregulation of candidate marker genes. Transcriptomic studies reveal that several genes downregulated after irradiation are conserved with other flatworms, and are related to cell cycle, DNA replication and repair functions. Furthermore, proliferative cells were isolated by cell sorting and also characterized transcriptomically. We found that the set of genes characteristic of proliferative cells agrees well with those downregulated during irradiation, and have a significant overlap with those expressed in planarian neoblasts or S. mansoni stem cells. Our study highlights that conserved mechanisms of stem cell biology may be functional in flatworms, suggesting that these could be relevant targets to evaluate in the control of parasitic species.

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.

Rhodopsin orphan GPCR20 interacts with neuropeptides and directs growth, sexual differentiation, and egg production in female Schistosoma mansoni

IMPORTANCE

Schistosomes cause schistosomiasis, one of the neglected tropical diseases as defined by the WHO. For decades, the treatment of schistosomiasis relies on a single drug, praziquantel. Due to its wide use, there is justified fear of resistance against this drug, and a vaccine is not available. Besides its biological relevance in signal transduction processes, the class of G protein-coupled receptors (GPCRs) is also well suited for drug design. Against this background, we characterized one GPCR of Schistosoma mansoni, SmGPCR20, at the molecular and functional level. We identified two potential neuropeptides (NPPs) as ligands, SmNPP26 and SmNPP40, and unraveled their roles, in combination with SmGPCR20, in neuronal processes controlling egg production, oogenesis, and growth of S. mansoni females. Since eggs are closely associated with the pathogenesis of schistosomiasis, our results contribute to the understanding of processes leading to egg production in schistosomes, which is under the control of pairing in this exceptional parasite.

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.

Schistosomiasis

Schistosomiasis is a major cause of morbidity in the world and almost 800 million people worldwide are at risk for schistosomiasis; it is second only to malaria as a major infectious disease. Globally, it is estimated that the disease affects more than 250 million people in 78 countries of the world and is responsible for some 280,000-500,000 deaths each year. The three major schistosomes infecting humans are Schistosoma mansoni, S. japonicum, and S. haematobium. This chapter covers a wide range of aspects of schistosomiasis, including basic biology of the parasites, epidemiology, immunopathology, treatment, control, vaccines, and genomics/proteomics. In this chapter, the reader will understand the significant toll this disease takes in terms of mortality and morbidity. A description of the various life stages of schistosomes is presented, which will be informative for both those unfamiliar with the disease and experienced scientists. Clinical and public health aspects are addressed that cover acute and chronic disease, diagnosis, current treatment regimens and alternative drugs, and schistosomiasis control programs. A brief overview of genomics and proteomics is included that details recent advances in the field that will help scientists investigate the molecular biology of schistosomes. The reader will take away an appreciation for general aspects of schistosomiasis and the current research advances.

Functional genomics in Spiralia

Our understanding of the mechanisms that modulate gene expression in animals is strongly biased by studying a handful of model species that mainly belong to three groups: Insecta, Nematoda and Vertebrata. However, over half of the animal phyla belong to Spiralia, a morphologically and ecologically diverse animal clade with many species of economic and biomedical importance. Therefore, investigating genome regulation in this group is central to uncovering ancestral and derived features in genome functioning in animals, which can also be of significant societal impact. Here, we focus on five aspects of gene expression regulation to review our current knowledge of functional genomics in Spiralia. Although some fields, such as single-cell transcriptomics, are becoming more common, the study of chromatin accessibility, DNA methylation, histone post-translational modifications and genome architecture are still in their infancy. Recent efforts to generate chromosome-scale reference genome assemblies for greater species diversity and optimise state-of-the-art approaches for emerging spiralian research systems will address the existing knowledge gaps in functional genomics in this animal group.

Fractionation of a sex-inducing substance from flatworms using open-column chromatography and reverse-phase high-performance liquid chromatography

A substance that sexualizes planarians, an ancestral group of parasitic flatworms, is widely present in planarians and parasitic flatworms. Here, we present a protocol for extracting and purifying the active fraction with sex-inducing activity. We describe steps for homogenization of flatworms, sample concentration, open-column chromatography, and reverse-phase high-performance liquid chromatography. We then detail a feeding bioassay to confirm sex-inducing activity. The obtained active fraction may positively affect parasitic flatworm sexual maturation and can be tested by adding it into the culture media. For complete details on the use and execution of this protocol, please refer to Sekii et al. (2023).1.

Schistosoma mansoni excretory-secretory products induce protein kinase signalling, hyperkinesia, and stem cell proliferation in the opposite sex

Adult male and female schistosomes in copula dwell within human blood vessels and lay eggs that cause the major Neglected Tropical Disease human schistosomiasis. How males and females communicate to each other is poorly understood; however, male-female physical interaction is known to be important. Here, we investigate whether excretory-secretory products (ESPs), released into the external milieu by mature Schistosoma mansoni, might induce responses in the opposite sex. We demonstrate that ESPs adhere to the surface of opposite sex worms inducing the activation of extracellular signal-regulated kinase (ERK) and p38 mitogen-activated protein kinase (p38 MAPK) pathways, particularly in the parasite tegument. Furthermore, we show that mature worms stimulated signalling in juvenile worms. Strikingly, we demonstrate that ESPs from the opposite sex promote stem cell proliferation, in an ERK- and p38 MAPK-dependent manner, in the tegument and within the testes of males, and the ovaries and vitellaria of females. Hyperkinesia also occurs following opposite sex ESP exposure. Our findings support the hypothesis that male and female schistosomes may communicate over distance to modulate key processes underlying worm development and disease progression, opening unique avenues for schistosomiasis control.

New technologies to study helminth development and host-parasite interactions

How parasites develop and survive, and how they stimulate or modulate host immune responses are important in understanding disease pathology and for the design of new control strategies. Microarray analysis and bulk RNA sequencing have provided a wealth of data on gene expression as parasites develop through different life-cycle stages and on host cell responses to infection. These techniques have enabled gene expression in the whole organism or host tissue to be detailed, but do not take account of the heterogeneity between cells of different types or developmental stages, nor the spatial organisation of these cells. Single-cell RNA-seq (scRNA-seq) adds a new dimension to studying parasite biology and host immunity by enabling gene profiling at the individual cell level. Here we review the application of scRNA-seq to establish gene expression cell atlases for multicellular helminths and to explore the expansion and molecular profile of individual host cell types involved in parasite immunity and tissue repair. Studying host-parasite interactions in vivo is challenging and we conclude this review by briefly discussing the applications of organoids (stem-cell derived mini-tissues) to examine host-parasite interactions at the local level, and as a potential system to study parasite development in vitro. Organoid technology and its applications have developed rapidly, and the elegant studies performed to date support the use of organoids as an alternative in vitro system for research on helminth parasites.

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.

Planarians to schistosomes: an overview of flatworm cell-types and regulators

Schistosomiasis remains a major neglected tropical disease that afflicts over 200 million people globally. Schistosomes, the aetiological agent of schistosomiasis, are parasitic flatworms that propagate between molluscan and mammalian hosts. Inside the mammalian host, schistosomes rapidly grow over 100-fold in size and develop into a sexually mature male or female that thrives in the bloodstream for several decades. Recent work has identified schistosome stem cells as the source that drives parasite transmission, reproduction and longevity. Moreover, studies have begun to uncover molecular programmes deployed by stem cells that are essential for tissue development and maintenance, parasite survival and immune evasion. Such programmes are reminiscent of neoblast-driven development and regeneration of planarians, the free-living flatworm relative of schistosomes. Over the last few decades, research in planarians has employed modern functional genomic tools that significantly enhanced our understanding of stem cell-driven animal development and regeneration. In this review, we take a broad stroke overview of major flatworm organ systems at the cellular and molecular levels. We summarize recent advances on genetic regulators that play critical roles in differentiation and maintenance of flatworm cell types. Finally, we provide perspectives on how investigation of basic parasite biology is critical to discovering new approaches to battle schistosomiasis.

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.

Sex-inducing effects toward planarians widely present among parasitic flatworms

Various parasitic flatworms infect vertebrates for sexual reproduction, often causing devastating diseases in their hosts. Consequently, flatworms are of great socioeconomic and biomedical importance. Although the cessation of parasitic flatworm sexual reproduction is a major target of anti-parasitic drug design, little is known regarding bioactive compounds controlling flatworm sexual maturation. Using the planarian Dugesia ryukyuensis, we observed that sex-inducing substances found in planarians are also widespread in parasitic flatworms, such as monogeneans and flukes (but not in tapeworms). Reverse-phase HPLC analysis revealed the sex-inducing substance(s) eluting around the tryptophan retention time in the fluke Calicophoron calicophorum, consistent with previous studies on the planarian Bipalium nobile, suggesting that the substance(s) is likely conserved among flatworms. Moreover, six of the 18 ovary-inducing substances identified via transcriptome and metabolome analyses are involved in purine metabolism. Our findings provide a basis for understanding and modifying the life cycles of various parasitic flatworms.

Proteomic profiling of serum extracellular vesicles identifies diagnostic markers for echinococcosis

Echinococcosis is a parasitic disease caused by the metacestodes of Echinococcus spp. The disease has a long latent period and is largely underdiagnosed, partially because of the lack of effective early diagnostic approaches. Using liquid chromatography-mass spectrometry, we profiled the serum-derived extracellular vesicles (EVs) of E. multilocularis-infected mice and identified three parasite-origin proteins, thioredoxin peroxidase 1 (TPx-1), transitional endoplasmic reticulum ATPase (TER ATPase), and 14-3-3, being continuously released by the parasites into the sera during the infection via EVs. Using ELISA, both TPx-1 and TER ATPase were shown to have a good performance in diagnosis of experimental murine echinococcosis as early as 10 days post infection and of human echinococcosis compared with that of control. Moreover, TER ATPase and TPx-1 were further demonstrated to be suitable for evaluation of the prognosis of patients with treatment. The present study discovers the potential of TER ATPase and TPx-1 as promising diagnostic candidates for echinococcosis.

Single-cell RNA-seq analyses show that long non-coding RNAs are conspicuously expressed in Schistosoma mansoni gamete and tegument progenitor cell populations

Schistosoma mansoni is a flatworm that causes schistosomiasis, a neglected tropical disease that affects over 200 million people worldwide. New therapeutic targets are needed with only one drug available for treatment and no vaccine. Long non-coding RNAs (lncRNAs) are transcripts longer than 200 nucleotides with low or no protein-coding potential. In other organisms, they have been shown as involved with reproduction, stem cell maintenance and drug resistance, and they tend to exhibit tissue-specific expression patterns. S. mansoni expresses thousands of lncRNA genes; however, the cell type expression patterns of lncRNAs in the parasite remain uncharacterized. Here, we have re-analyzed publicly available single-cell RNA-sequencing (scRNA-seq) data obtained from adult S. mansoni to identify the lncRNAs signature of adult schistosome cell types. A total of 8023 lncRNAs (79% of all lncRNAs) were detected. Analyses of the lncRNAs expression profiles in the cells using statistically stringent criteria were performed to identify 74 lncRNA gene markers of cell clusters. Male gamete and tegument progenitor lineages clusters contained most of the cluster-specific lncRNA markers. We also identified lncRNA markers of specific neural clusters. Whole-mount in situ hybridization (WISH) and double fluorescence in situ hybridization were used to validate the cluster-specific expression of 13 out of 16 selected lncRNA genes (81%) in the male and female adult parasite tissues; for one of these 16 gene loci, probes for two different lncRNA isoforms were used, which showed differential isoform expression in testis and ovary. An atlas of the expression profiles across the cell clusters of all lncRNAs detected in our analysis is available as a public website resource (http://verjolab.usp.br:8081). The results presented here give strong support to a tissue-specific expression and to a regulated expression program of lncRNAs in S. mansoni. This will be the basis for further exploration of lncRNA genes as potential therapeutic targets.

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.

A chromosome-level genome of the human blood fluke Schistosoma japonicum identifies the genomic basis of host-switching

The evolution and adaptation of S. japonicum, a zoonotic parasite that causes human schistosomiasis, remain unclear because of the lack of whole-genome data. We construct a chromosome-level S. japonicum genome and analyze it together with 72 samples representing six populations of the entire endemic region. We observe a Taiwan zoophilic lineage splitting from zoonotic populations ∼45,000 years ago, consistent with the divergent history of their intermediate hosts. Interestingly, we detect a severe population bottleneck in S. japonicum, largely coinciding with human history in Asia during the last glacial maximum. We identify several genomic regions underlying natural selection, including GATAD2A and Lmln, both showing remarkable differentiation among different areas. RNAi knockdown suggests association of GATAD2A with parasite development and infection in definitive hosts, while Lmln relates to the specificity of the intermediate hosts. Our study provides insights into the evolution of S. japonicum and serves as a resource for further studies.

Comparative transcriptome profiles of Schistosoma japonicum larval stages: Implications for parasite biology and host invasion

Schistosoma japonicum is prevalent in Asia with a wide mammalian host range, which leads to highly harmful zoonotic parasitic diseases. Most previous transcriptomic studies have been performed on this parasite, but mainly focus on stages inside the mammalian host. Moreover, few larval transcriptomic data are available in public databases. Here we mapped the detailed transcriptome profiles of four S. japonicum larval stages including eggs, miracidia, sporocysts and cercariae, providing a comprehensive development picture outside of the mammalian host. By analyzing the stage-specific/enriched genes, we identified functional genes associated with the biological characteristic at each stage: e.g. we observed enrichment of genes necessary for DNA replication only in sporocysts, while those involved in proteolysis were upregulated in sporocysts and/or cercariae. This data indicated that miracidia might use leishmanolysin and neprilysin to penetrate the snail, while elastase (SjCE2b) and leishmanolysin might contribute to the cercariae invasion. The expression profile of stem cell markers revealed potential germinal cell conversion during larval development. Additionally, our analysis indicated that tandem duplications had driven the expansion of the papain family in S. japonicum. Notably, all the duplicated cathepsin B-like proteases were highly expressed in cercariae. Utilizing our 3rd version of S. japonicum genome, we further characterized the alternative splicing profiles throughout these four stages. Taken together, the present study provides compressive gene expression profiles of S. japonicum larval stages and identifies a set of genes that might be involved in intermediate and definitive host invasion.

Daily rhythms in gene expression of the human parasite Schistosoma mansoni

Background

The consequences of the earth’s daily rotation have led to 24-h biological rhythms in most organisms. Even some parasites are known to have daily rhythms, which, when in synchrony with host rhythms, can optimise their fitness. Understanding these rhythms may enable the development of control strategies that take advantage of rhythmic vulnerabilities. Recent work on protozoan parasites has revealed 24-h rhythms in gene expression, drug sensitivity and the presence of an intrinsic circadian clock; however, similar studies on metazoan parasites are lacking. To address this, we investigated if a metazoan parasite has daily molecular oscillations, whether they reveal how these longer-lived organisms can survive host daily cycles over a lifespan of many years and if animal circadian clock genes are present and rhythmic. We addressed these questions using the human blood fluke Schistosoma mansoni that lives in the vasculature for decades and causes the tropical disease schistosomiasis.

Results

Using round-the-clock transcriptomics of male and female adult worms collected from experimentally infected mice, we discovered that ~ 2% of its genes followed a daily pattern of expression. Rhythmic processes included a stress response during the host’s active phase and a ‘peak in metabolic activity’ during the host’s resting phase. Transcriptional profiles in the female reproductive system were mirrored by daily patterns in egg laying (eggs are the main drivers of the host pathology). Genes cycling with the highest amplitudes include predicted drug targets and a vaccine candidate. These 24-h rhythms may be driven by host rhythms and/or generated by a circadian clock; however, orthologs of core clock genes are missing and secondary clock genes show no 24-h rhythmicity.

Conclusions

There are daily rhythms in the transcriptomes of adult S. mansoni, but they appear less pronounced than in other organisms. The rhythms reveal temporally compartmentalised internal processes and host interactions relevant to within-host survival and between-host transmission. Our findings suggest that if these daily rhythms are generated by an intrinsic circadian clock then the oscillatory mechanism must be distinct from that in other animals. We have shown which transcripts oscillate at this temporal scale and this will benefit the development and delivery of treatments against schistosomiasis.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-021-01189-9.

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.

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.

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.

SchistoCyte Atlas: A Single-Cell Transcriptome Resource for Adult Schistosomes

Schistosomes cause untold disease and disability in the developing world. Here, we introduce SchistoCyte Atlas, a web-based platform for exploring gene expression at single-cell resolution in adult Schistosoma mansoni. Similar resources accessible to non-specialists across the globe will expedite our ability to understand the biology of these devastating parasites.

Mapping single-cell atlases throughout Metazoa unravels cell type evolution

Comparing single-cell transcriptomic atlases from diverse organisms can elucidate the origins of cellular diversity and assist the annotation of new cell atlases. Yet, comparison between distant relatives is hindered by complex gene histories and diversifications in expression programs. Previously, we introduced the self-assembling manifold (SAM) algorithm to robustly reconstruct manifolds from single-cell data (Tarashansky et al., 2019). Here, we build on SAM to map cell atlas manifolds across species. This new method, SAMap, identifies homologous cell types with shared expression programs across distant species within phyla, even in complex examples where homologous tissues emerge from distinct germ layers. SAMap also finds many genes with more similar expression to their paralogs than their orthologs, suggesting paralog substitution may be more common in evolution than previously appreciated. Lastly, comparing species across animal phyla, spanning sponge to mouse, reveals ancient contractile and stem cell families, which may have arisen early in animal evolution.

Single-cell deconstruction of stem-cell-driven schistosome development

Schistosomes cause one of the most devastating neglected tropical diseases, schistosomiasis. Their transmission is accomplished through a complex life cycle with two obligate hosts and requires multiple radically different body plans specialized for infecting and reproducing in each host. Recent single-cell transcriptomic studies on several schistosome body plans provide a comprehensive map of their cell types, which include stem cells and their differentiated progeny along an intricate developmental hierarchy. This progress not only extends our understanding of the basic biology of the schistosome life cycle but can also inform new therapeutic and preventive strategies against the disease, as blocking the development of specific cell types through genetic manipulations has shown promise in inhibiting parasite survival, growth, and reproduction.

Mitochondrial state determines functionally divergent stem cell population in planaria

Mitochondrial state changes were shown to be critical for stem cell function. However, variation in the mitochondrial content in stem cells and the implication, if any, on differentiation is poorly understood. Here, using cellular and molecular studies, we show that the planarian pluripotent stem cells (PSCs) have low mitochondrial mass compared with their progenitors. Transplantation experiments provided functional validation that neoblasts with low mitochondrial mass are the true PSCs. Further, the mitochondrial mass correlated with OxPhos and inhibiting the transition to OxPhos dependent metabolism in cultured cells resulted in higher PSCs. In summary, we show that low mitochondrial mass is a hallmark of PSCs in planaria and provide a mechanism to isolate live, functionally active, PSCs from different cell cycle stages (G0/G1 and S, G2/M). Our study demonstrates that the change in mitochondrial metabolism, a feature of PSCs is conserved in planaria and highlights its role in organismal regeneration.

•

Mitochondrial state differs between stem (X1) and differentiated (Xins) cells

•

X1 cells with low MTG are enriched for pluripotent cells compared with high MTG cells

•

MTG-based sorting yields functional neoblasts from G1, S/G2/M phase of cell cycle

•

Inhibition of mitochondrial activity affects neoblast differentiation in vitro

Innovations and Advances in Schistosome Stem Cell Research

Schistosomes infect about 250 million people globally causing the devastating and persistent disease of schistosomiasis. These blood flukes have a complicated life cycle involving alternating infection of freshwater snail intermediate and definitive mammalian hosts. To survive and flourish in these diverse environments, schistosomes transition through a number of distinct life-cycle stages as a result of which they change their body plan in order to quickly adapt to each new environment. Current research suggests that stem cells, present in adults and larvae, are key in aiding schistosomes to facilitate these changes. Given the recent advances in our understanding of schistosome stem cell biology, we review the key roles that two major classes of cells play in the different life cycle stages during intramolluscan and intramammalian development; these include the germinal cells of sporocysts involved in asexual reproduction in molluscan hosts and the neoblasts of adult worms involved in sexual reproduction in human and other mammalian hosts. These studies shed considerable new light in revealing the stem cell heterogeneity driving the propagation of the schistosome life cycle. We also consider the possibility and value of establishing stem cell lines in schistosomes to advance schistosomiasis research. The availability of such self-renewable resources will provide new platforms to study stem cell behavior and regulation, and to address fundamental aspects of schistosome biology, reproductive development and survival. In turn, such studies will create new avenues to unravel individual gene function and to optimize genome-editing processes in blood flukes, which may lead to the design of novel intervention strategies for schistosomiasis.

Schistosoma mansoni alter transcription of immunomodulatory gene products following in vivo praziquantel exposure

Control of the neglected tropical disease schistosomiasis relies almost entirely on praziquantel (PZQ) monotherapy. How PZQ clears parasite infections remains poorly understood. Many studies have examined the effects of PZQ on worms cultured in vitro, observing outcomes such as muscle contraction. However, conditions worms are exposed to in vivo may vary considerably from in vitro experiments given the short half-life of PZQ and the importance of host immune system engagement for drug efficacy in animal models. Here, we investigated the effects of in vivo PZQ exposure on Schistosoma mansoni. Measurement of pro-apoptotic caspase activation revealed that worm death occurs only after parasites shift from the mesenteric vasculature to the liver, peaking 24 hours after drug treatment. This indicates that PZQ is not directly schistocidal, since PZQ’s half-life is ~2 hours in humans and ~30 minutes in mice, and focuses attention on parasite interactions with the host immune system following the shift of worms to the liver. RNA-Seq of worms harvested from mouse livers following sub-lethal PZQ treatment revealed drug-evoked changes in the expression of putative immunomodulatory and anticoagulant gene products. Several of these gene products localized to the schistosome esophagus and may be secreted into the host circulation. These include several Kunitz-type protease inhibitors, which are also found in the secretomes of other blood feeding animals. These transcriptional changes may reflect mechanisms of parasite immune-evasion in response to chemotherapy, given the role of complement-mediated attack and the host innate/humoral immune response in parasite elimination. One of these isoforms, SmKI-1, has been shown to exhibit immunomodulatory and anti-coagulant properties. These data provide insight into the effect of in vivo PZQ exposure on S. mansoni, and the transcriptional response of parasites to the stress of chemotherapy.

The disease schistosomiasis is caused by parasitic worms that live within the circulatory system. While this disease infects over 200 million people worldwide, treatment relies almost entirely on one drug, praziquantel, whose mechanism is poorly understood. In this study, we analyzed the effects of praziquantel treatment on the gene expression of parasites harvested from mice treated with praziquantel chemotherapy. Despite the rapid action of the drug on worms in vitro, we found that key outcomes in vivo (measurement of cell death and changes in gene expression) occurred relatively late (12+ hours after drug administration). We found that worms increased the expression of immunomodulatory gene products in response to praziquantel, including a Kunitz-type protease inhibitor that localized to the worm esophagus and may be secreted to the external host environment. These are an intriguing class of proteins, because they display anti-coagulant and immunomodulatory properties. Up-regulation of these gene products may reflect a parasite mechanism of immune-evasion in response to chemotherapy. This research provides insight into the mechanism of praziquantel by observing the effect of this drug on worms within the context of the host immune system.

Single-cell analysis of Schistosoma mansoni identifies a conserved genetic program controlling germline stem cell fate

Schistosomes are parasitic flatworms causing one of the most prevalent infectious diseases from which millions of people are currently suffering. These parasites have high fecundity and their eggs are both the transmissible agents and the cause of the infection-associated pathology. Given its biomedical significance, the schistosome germline has been a research focus for more than a century. Nonetheless, molecular mechanisms that regulate its development are only now being understood. In particular, it is unknown what balances the fate of germline stem cells (GSCs) in producing daughter stem cells through mitotic divisions versus gametes through meiosis. Here, we perform single-cell RNA sequencing on juvenile schistosomes and capture GSCs during de novo gonadal development. We identify a genetic program that controls the proliferation and differentiation of GSCs. This program centers around onecut, a homeobox transcription factor, and boule, an mRNA binding protein. Their expressions are mutually dependent in the schistosome male germline, and knocking down either of them causes over-proliferation of GSCs and blocks germ cell differentiation. We further show that this germline-specific regulatory program is conserved in the planarian, schistosome's free-living evolutionary cousin, but the function of onecut has changed during evolution to support GSC maintenance.

Single-cell atlas of the first intra-mammalian developmental stage of the human parasite Schistosoma mansoni

Over 250 million people suffer from schistosomiasis, a tropical disease caused by parasitic flatworms known as schistosomes. Humans become infected by free-swimming, water-borne larvae, which penetrate the skin. The earliest intra-mammalian stage, called the schistosomulum, undergoes a series of developmental transitions. These changes are critical for the parasite to adapt to its new environment as it navigates through host tissues to reach its niche, where it will grow to reproductive maturity. Unravelling the mechanisms that drive intra-mammalian development requires knowledge of the spatial organisation and transcriptional dynamics of different cell types that comprise the schistomulum body. To fill these important knowledge gaps, we perform single-cell RNA sequencing on two-day old schistosomula of Schistosoma mansoni. We identify likely gene expression profiles for muscle, nervous system, tegument, oesophageal gland, parenchymal/primordial gut cells, and stem cells. In addition, we validate cell markers for all these clusters by in situ hybridisation in schistosomula and adult parasites. Taken together, this study provides a comprehensive cell-type atlas for the early intra-mammalian stage of this devastating metazoan parasite.

A single-cell RNA-seq atlas of Schistosoma mansoni identifies a key regulator of blood feeding

Schistosomiasis is a neglected tropical disease that infects 240 million people. With no vaccines and only one drug available, new therapeutic targets are needed. The causative agents, schistosomes, are intravascular flatworm parasites that feed on blood and lay eggs, resulting in pathology. The function of the parasite's various tissues in successful parasitism are poorly understood, hindering identification of therapeutic targets. Using single-cell RNA sequencing (RNA-seq), we characterize 43,642 cells from the adult schistosome and identify 68 distinct cell populations, including specialized stem cells that maintain the parasite's blood-digesting gut. These stem cells express the gene hnf4, which is required for gut maintenance, blood feeding, and pathology in vivo. Together, these data provide molecular insights into the organ systems of this important pathogen and identify potential therapeutic targets.

The esophageal gland mediates host immune evasion by the human parasite Schistosoma mansoni

Schistosomes are parasitic flatworms that cause schistosomiasis, a neglected tropical disease affecting over 200 million people. Schistosomes develop multiple body plans while navigating their complex life cycle, which involves two different hosts: a mammalian definitive host and a molluscan intermediate host. Their survival and propagation depend upon proliferation and differentiation of stem cells necessary for parasite homeostasis and reproduction. Infective larvae released from snails carry a handful of stem cells that serve as the likely source of new tissues as the parasite adapts to life inside the mammalian host; however, the role of these stem cells during this critical life cycle stage remains unclear. Here, we characterize stem cell fates during early intramammalian development. Surprisingly, we find that the esophageal gland, an accessory organ of the digestive tract, develops before the rest of the digestive system is formed and blood feeding is initiated, suggesting a role in processes beyond nutrient uptake. To explore such a role, we examine schistosomes that lack the esophageal gland due to knockdown of a forkhead-box transcription factor, Sm-foxA, which blocks development and maintenance of the esophageal gland, without affecting the development of other somatic tissues. Intriguingly, schistosomes lacking the esophageal gland die after transplantation into naive mice, but survive in immunodeficient mice lacking B cells. We show that parasites lacking the esophageal gland are unable to lyse ingested immune cells within the esophagus before passing them into the gut. These results unveil an immune-evasion mechanism mediated by the esophageal gland, which is essential for schistosome survival and pathogenesis.

Schistosome TRP channels: An appraisal

Ion channels underlie electrical excitability in cells and are essential for a variety of functions, most notably neuromuscular and sensory activity. They are also validated targets for a preponderance of approved anthelmintic compounds. Transient receptor potential (TRP) channels constitute an ion channel superfamily whose members play important roles in sensory signaling, regulation of ion homeostasis, organellar trafficking, and other key cellular and organismal activities. Unlike most other ion channels, TRP channels are often polymodal, gated by a variety of mechanisms. Furthermore, TRP channels fall into several classes or subtypes based on sequence and structure. Until recently, there had been very little investigation of the properties and functions of TRP channels from parasitic helminths, including schistosomes, but that situation has changed in the past few years. Indeed, it is now clear that at least some schistosome TRP channels exhibit unusual pharmacological properties, and, intriguingly, both a mammalian and a schistosome TRP channel are activated by praziquantel, the current antischistosomal drug of choice. With the latest release of the Schistosoma mansoni genome database, several changes in predicted TRP channel sequences appeared, some of which were significant. This review updates and reassesses the TRP channel repertoire in S. mansoni, examines recent findings regarding these potential therapeutic targets, and provides guideposts for some of the physiological functions that may be mediated by these channels in schistosomes.

Dynamic Expression of Long Non-Coding RNAs Throughout Parasite Sexual and Neural Maturation in Schistosoma japonicum

Schistosoma japonicum is a flatworm that causes schistosomiasis, a neglected tropical disease. S. japonicum RNA-Seq analyses has been previously reported in the literature on females and males obtained during sexual maturation from 14 to 28 days post-infection in mouse, resulting in the identification of protein-coding genes and pathways, whose expression levels were related to sexual development. However, this work did not include an analysis of long non-coding RNAs (lncRNAs). Here, we applied a pipeline to identify and annotate lncRNAs in 66 S. japonicum RNA-Seq publicly available libraries, from different life-cycle stages. We also performed co-expression analyses to find stage-specific lncRNAs possibly related to sexual maturation. We identified 12,291 S. japonicum expressed lncRNAs. Sequence similarity search and synteny conservation indicated that some 14% of S. japonicum intergenic lncRNAs have synteny conservation with S. mansoni intergenic lncRNAs. Co-expression analyses showed that lncRNAs and protein-coding genes in S. japonicum males and females have a dynamic co-expression throughout sexual maturation, showing differential expression between the sexes; the protein-coding genes were related to the nervous system development, lipid and drug metabolism, and overall parasite survival. Co-expression pattern suggests that lncRNAs possibly regulate these processes or are regulated by the same activation program as that of protein-coding genes.

LncRNAs in molluscan and mammalian stages of parasitic schistosomes are developmentally-regulated and coordinately expressed with protein-coding genes

Despite the low level expression of some long noncoding RNAs (lncRNAs), the differential expression of specific lncRNAs plays important roles during the development of many organisms. Schistosomes, parasitic flatworms that are responsible for schistosomiasis, infects over 200 million people resulting in chronic disease and hundreds of thousands of deaths. Schistosomes have a complex life cycle that transitions between molluscan and mammalian hosts. In a molluscan snail host, the sporocyst stage develops over 5 weeks undergoing asexual reproduction to give rise to free-swimming and infectious cercariae that penetrate human skin and eventually mature into egg producing worms in mammals. The tight integration of the sporocyst to the snail host hepatopancreas hinders the -omics study in the molluscan stage, so the sporocyst transcriptome has only been examined for lncRNAs in immature in vitro samples. Here we analyzed the in vivo mature sporocyst transcriptome to identify 4,930 total lncRNAs between the molluscan and mammalian stages of the parasite. We further demonstrate that the lncRNAs are differentially expressed in a development-dependent manner. In addition, we constructed a co-expression correlation network between lncRNAs and protein-coding (PC) genes that was used to identify clusters of lncRNA transcripts with potential functional relevance. We also describe lncRNA-lncRNA and lncRNA-kinome correlations that identify lncRNAs with prospective roles in gene regulation. Finally, our results show clear differential expression patterns of lncRNAs in host-dependent development stages of S. mansoni and ascribe potential functional roles in development based on predicted intracellular interaction.

A single-parasite transcriptional atlas of Toxoplasma Gondii reveals novel control of antigen expression

Toxoplasma gondii, a protozoan parasite, undergoes a complex and poorly understood developmental process that is critical for establishing a chronic infection in its intermediate hosts. Here, we applied single-cell RNA-sequencing (scRNA-seq) on >5,400 Toxoplasma in both tachyzoite and bradyzoite stages using three widely studied strains to construct a comprehensive atlas of cell-cycle and asexual development, revealing hidden states and transcriptional factors associated with each developmental stage. Analysis of SAG1-related sequence (SRS) antigenic repertoire reveals a highly heterogeneous, sporadic expression pattern unexplained by measurement noise, cell cycle, or asexual development. Furthermore, we identified AP2IX-1 as a transcription factor that controls the switching from the ubiquitous SAG1 to rare surface antigens not previously observed in tachyzoites. In addition, comparative analysis between Toxoplasma and Plasmodium scRNA-seq results reveals concerted expression of gene sets, despite fundamental differences in cell division. Lastly, we built an interactive data-browser for visualization of our atlas resource.