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

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.

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.

Abundant genetic variation is retained in many laboratory schistosome populations

Schistosomes are obligately sexual blood flukes that can be maintained in the laboratory using freshwater snails as intermediate and rodents as definitive hosts. The genetic composition of laboratory schistosome populations is poorly understood: whether genetic variation has been purged due to serial inbreeding or retained is unclear. We sequenced 19 - 24 parasites from each of five laboratory Schistosoma mansoni populations and compared their genomes with published exome data from four S. mansoni field populations. We found abundant genomic variation (0.897 - 1.22 million variants) within laboratory populations: these retained on average 49% (π = 3.27e-04 - 8.94e-04) of the nucleotide diversity observed in the four field parasite populations (π = 1.08e-03 - 2.2e-03). However, the pattern of variation was very different in laboratory and field populations. Tajima's D was positive in all laboratory populations except SmBRE, indicative of recent population bottlenecks, but negative in all field populations. Current effective population size estimates of laboratory populations were lower (2 - 258) compared to field populations (3,174 - infinity). The distance between markers at which linkage disequilibrium (LD) decayed to 0.5 was longer in laboratory populations (59 bp - 180 kb) compared to field populations (9 bp - 9.5 kb). SmBRE was the least variable; this parasite also shows low fitness across the lifecycle, consistent with inbreeding depression. The abundant genetic variation present in most laboratory schistosome populations has several important implications: (i) measurement of parasite phenotypes, such as drug resistance, using laboratory parasite populations will determine average values and underestimate trait variation; (ii) genome-wide association studies (GWAS) can be conducted in laboratory schistosome populations by measuring phenotypes and genotypes of individual worms; (iii) genetic drift may lead to divergence in schistosome populations maintained in different laboratories. We conclude that the abundant genetic variation retained within many laboratory schistosome populations can provide valuable, untapped opportunities for schistosome research.

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.

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.

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.

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.

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.

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.

Long non-coding RNAs as possible therapeutic targets in protozoa, and in Schistosoma and other helminths

Long non-coding RNAs (lncRNAs) emerged in the past 20 years due to massive amounts of scientific data regarding transcriptomic analyses. They have been implicated in a plethora of cellular processes in higher eukaryotes. However, little is known about lncRNA possible involvement in parasitic diseases, with most studies only detecting their presence in parasites of human medical importance. Here, we review the progress on lncRNA studies and their functions in protozoans and helminths. In addition, we show an example of knockdown of one lncRNA in Schistosoma mansoni, SmLINC156349, which led to in vitro parasite adhesion, motility, and pairing impairment, with a 20% decrease in parasite viability and 33% reduction in female oviposition. Other observed phenotypes were a decrease in the proliferation rate of both male and female worms and their gonads, and reduced female lipid and vitelline droplets that are markers for well-developed vitellaria. Impairment of female worms' vitellaria in SmLINC156349-silenced worms led to egg development deficiency. All those results demonstrate the great potential of the tools and methods to characterize lncRNAs as potential new therapeutic targets. Further, we discuss the challenges and limitations of current methods for studying lncRNAs in parasites and possible solutions to overcome them, and we highlight the future directions of this exciting field.

Germinal development in embryonic rediae of the hemiuroid digenean Bunocotyle progenetica: an ultrastructural study

Embryonic development of reproductive organs in rediae of the digenean Bunocotyle progenetica was studied using transmission electron microscopy. The germinal primordium becomes morphologically distinct in early embryos as a weakly separated cell mass with a forming cavity. It consists of undifferentiated, differentiating, and supporting cells. As embryos develop, the supporting cells form a wall around the enlarging cavity. Other cells of the germinal primordium are incorporated into the wall as solitary cells or as small cell aggregations. Those situated posteriorly give rise to an incipient germinal mass functioning during postembryonic development. Undifferentiated and differentiating cells in the middle and the anterior part of the primordium ensure a considerable growth of the cavity wall, which incorporates solitary germinal cells. In advanced embryonic rediae, these cells mature, cleave, and give rise to germinal balls, which enter the forming brood cavity. In the most mature embryonic rediae, all these early cercarial embryos reside in a brood cavity, which is lined by that time with a syncytium continuous with the supporting tissue of the incipient germinal mass. Based on our results and the literature data, we suggest that the morphogenesis of the reproductive apparatus of the daughter parthenitae in hemiuroid digeneans may be characterized by (1) emergence of an incipient brood cavity within the germinal primordium, (2) formation of the cavity lining from the cells of the germinal primordium, (3) fragmentation and uneven distribution of the germinal material of the germinal primordium around the cavity, and (4) an anticipatory development of some of this germinal material.

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.