Different metazoan parasites, different transcriptomic responses, with new insights on parasitic castration by digenetic trematodes in the schistosome vector snail Biomphalaria glabrata.

BACKGROUND: Gastropods of the genus Biomphalaria (Family Planorbidae) are exploited as vectors by Schistosoma mansoni, the most common causative agent of human intestinal schistosomiasis. Using improved genomic resources, overviews of how Biomphalaria responds to S. mansoni and other metazoan parasites can provide unique insights into the reproductive, immune, and other systems of invertebrate hosts, and their responses to parasite challenges. RESULTS: Using Illumina-based RNA-Seq, we compared the responses of iM line B. glabrata at 2, 8, and 40 days post-infection (dpi) to single infections with S. mansoni, Echinostoma paraensei (both digenetic trematodes) or Daubaylia potomaca (a nematode parasite of planorbid snails). Responses were compared to unexposed time-matched control snails. We observed: (1) each parasite provoked a distinctive response with a predominance of down-regulated snail genes at all time points following exposure to either trematode, and of up-regulated genes at 8 and especially 40dpi following nematode exposure; (2) At 2 and 8dpi with either trematode, several snail genes associated with gametogenesis (particularly spermatogenesis) were down-regulated. Regarding the phenomenon of trematode-mediated parasitic castration in molluscs, we define for the first time a complement of host genes that are targeted, as early as 2dpi when trematode larvae are still small; (3) Differential gene expression of snails with trematode infection at 40dpi, when snails were shedding cercariae, was unexpectedly modest and revealed down-regulation of genes involved in the production of egg mass proteins and peptide processing; and (4) surprisingly, D. potomaca provoked up-regulation at 40dpi of many of the reproduction-related snail genes noted to be down-regulated at 2 and 8dpi following trematode infection. Happening at a time when B. glabrata began to succumb to D. potomaca, we hypothesize this response represents an unexpected form of fecundity compensation. We also document expression patterns for other Biomphalaria gene families, including fibrinogen domain-containing proteins (FReDs), C-type lectins, G-protein coupled receptors, biomphalysins, and protease and protease inhibitors. CONCLUSIONS: Our study is relevant in identifying several genes involved in reproduction that are targeted by parasites in the vector snail B. glabrata and that might be amenable to manipulation to minimize their ability to serve as vectors of schistosomes.

Metabolism of FAD, FMN and riboflavin (vitamin B2) in the human parasitic blood fluke Schistosoma mansoni.

BACKGROUND: Schistosomiasis is a parasitic disease caused by trematodes of the genus Schistosoma. The intravascular worms acquire the nutrients necessary for their survival from host blood. Since all animals are auxotrophic for riboflavin (vitamin B2), schistosomes too must import it to survive. Riboflavin is an essential component of the coenzymes flavin mononucleotide (FMN) and flavin adenine dinucleotide (FAD); these support key functions of dozens of flavoenzymes. METHODS: Here, using a combination of metabolomics, enzyme kinetics and in silico molecular analysis, we focus on the biochemistry of riboflavin and its metabolites in Schistosoma mansoni (Sm). RESULTS: We show that when schistosomes are incubated in murine plasma, levels of FAD decrease over time while levels of FMN increase. We show that live schistosomes cleave exogenous FAD to generate FMN and this ability is significantly blocked when expression of the surface nucleotide pyrophosphatase/phosphodiesterase ectoenzyme SmNPP5 is suppressed using RNAi. Recombinant SmNPP5 cleaves FAD with a Km of 178 ± 5.9 µM and Kcat/Km of 324,734 ± 36,347 M- 1.S- 1. The FAD-dependent enzyme IL-4I1 drives the oxidative deamination of phenylalanine to produce phenylpyruvate and H2O2. Since schistosomes are damaged by H2O2, we determined if SmNPP5 could impede H2O2 production by blocking IL-4I1 action in vitro. We found that this was not the case; covalently bound FAD on IL-4I1 appears inaccessible to SmNPP5. We also report that live schistosomes can cleave exogenous FMN to generate riboflavin and this ability is significantly impeded when expression of a second surface ectoenzyme (alkaline phosphatase, SmAP) is suppressed. Recombinant SmAP cleaves FMN with a Km of 3.82 ± 0.58 mM and Kcat/Km of 1393 ± 347 M- 1.S- 1. CONCLUSIONS: The sequential hydrolysis of FAD by tegumental ecto-enzymes SmNPP5 and SmAP can generate free vitamin B2 around the worms from where it can be conveniently imported by the recently described schistosome riboflavin transporter SmaRT. Finally, we identified in silico schistosome homologs of enzymes that are involved in intracellular vitamin B2 metabolism. These are riboflavin kinase (SmRFK) as well as FAD synthase (SmFADS); cDNAs encoding these two enzymes were cloned and sequenced. SmRFK is predicted to convert riboflavin to FMN while SmFADS could further act on FMN to regenerate FAD in order to facilitate robust vitamin B2-dependent metabolism in schistosomes.

Characterizing genetic variation on the Z chromosome in Schistosoma japonicum reveals host-parasite co-evolution.

BACKGROUND: Schistosomiasis is a neglected tropical disease that afflicts millions of people worldwide; it is caused by Schistosoma, the only dioecious flukes with ZW systems. Schistosoma japonicum is endemic to Asia; the Z chromosome of S. japonicum comprises one-quarter of the entire genome. Detection of positive selection using resequencing data to understand adaptive evolution has been applied to a variety of pathogens, including S. japonicum. However, the contribution of the Z chromosome to evolution and adaptation is often neglected. METHODS: We obtained 1,077,526 high-quality SNPs on the Z chromosome in 72 S. japonicum using re-sequencing data publicly. To examine the faster Z effect, we compared the sequence divergence of S. japonicum with two closely related species, Schistosoma haematobium and S. mansoni. Genetic diversity was compared between the Z chromosome and autosomes in S. japonicum by calculating the nucleotide diversity (π) and Dxy values. Population structure was also assessed based on PCA and structure analysis. Besides, we employed multiple methods including Tajima's D, FST, iHS, XP-EHH, and CMS to detect positive selection signals on the Z chromosome. Further RNAi knockdown experiments were performed to investigate the potential biological functions of the candidate genes. RESULTS: Our study found that the Z chromosome of S. japonicum showed faster evolution and more pronounced genetic divergence than autosomes, although the effect may be smaller than the variation among genes. Compared with autosomes, the Z chromosome in S. japonicum had a more pronounced genetic divergence of sub-populations. Notably, we identified a set of candidate genes associated with host-parasite co-evolution. In particular, LCAT exhibited significant selection signals within the Taiwan population. Further RNA interference experiments suggested that LCAT is necessary for S. japonicum survival and propagation in the definitive host. In addition, we identified several genes related to the specificity of the intermediate host in the C-M population, including Rab6 and VCP, which are involved in adaptive immune evasion to the host. CONCLUSIONS: Our study provides valuable insights into the adaptive evolution of the Z chromosome in S. japonicum and further advances our understanding of the co-evolution of this medically important parasite and its hosts.

Contribution of parasite and host genotype to immunopathology of schistosome infections.

BACKGROUND: The role of pathogen genotype in determining disease severity and immunopathology has been studied intensively in microbial pathogens including bacteria, fungi, protozoa and viruses but is poorly understood in parasitic helminths. The medically important blood fluke Schistosoma mansoni is an excellent model system to study the impact of helminth genetic variation on immunopathology. Our laboratory has demonstrated that laboratory schistosome populations differ in sporocyst growth and cercarial production in the intermediate snail host and worm establishment and fecundity in the vertebrate host. Here, we (i) investigate the hypothesis that schistosome genotype plays a significant role in immunopathology and related parasite life history traits in the vertebrate mouse host and (ii) quantify the relative impact of parasite and host genetics on infection outcomes. METHODS: We infected BALB/c and C57BL/6 mice with four different laboratory schistosome populations from Africa and the Americas. We quantified disease progression in the vertebrate host by measuring body weight and complete blood count (CBC) with differential over a 12-week infection period. On sacrifice, we assessed parasitological (egg and worm counts, fecundity), immunopathological (organ measurements and histopathology) and immunological (CBC with differential and cytokine profiles) characteristics to determine the impact of parasite and host genetics. RESULTS: We found significant variation between parasite populations in worm numbers, fecundity, liver and intestine egg counts, liver and spleen weight, and fibrotic area but not in granuloma size. Variation in organ weight was explained by egg burden and intrinsic parasite factors independent of egg burden. We found significant variation between infected mouse lines in cytokine levels (IFN-γ, TNF-α), eosinophils, lymphocytes and monocyte counts. CONCLUSIONS: This study showed that both parasite and host genotype impact the outcome of infection. While host genotype explains most of the variation in immunological traits, parasite genotype explains most of the variation in parasitological traits, and both host and parasite genotypes impact immunopathology outcomes.

Winners vs. losers: Schistosoma mansoni intestinal and liver eggs exhibit striking differences in gene expression and immunogenicity.

The eggs of the blood fluke Schistosoma mansoni are the main cause of the clinical manifestations of chronic schistosomiasis. After laying, the egg "winners" attach to the endothelium of the mesenteric vein and, after a period of development, induce the growth of a small granuloma, which facilitates their passage to the intestinal lumen. Egg "losers" carried by the bloodstream to non-specific tissues also undergo full development and induce large granuloma formation, but their life ends there. Although these trapped eggs represent a dead end in the parasite life cycle, the vast majority of studies attempting to describe the biology of the S. mansoni eggs have studied these liver-trapped "losers" instead of migrating intestinal "winners". This raises the fundamental question of how these eggs differ. With robust comparative transcriptomic analysis performed on S. mansoni eggs isolated 7 weeks post infection, we show that gene expression is critically dependent on tissue localization, both in the early and late stages of development. While mitochondrial genes and venom allergen-like proteins are significantly upregulated in mature intestinal eggs, well-described egg immunomodulators IPSE/alpha-1 and omega-1, together with micro-exon genes, are predominantly expressed in liver eggs. In addition, several proteases and protease inhibitors previously implicated in egg-host interactions display clear tissue-specific gene expression patterns. These major differences in gene expression could be then reflected in the observed different ability of liver and intestinal soluble egg antigens to elicit host immune responses and in the shorter viability of miracidia hatched from liver eggs. Our comparative analysis provides a new perspective on the biology of parasite's eggs in the context of their development and tissue localization. These findings could contribute to a broader and more accurate understanding of parasite eggs interactions with the host, which have historically been often restricted to liver eggs and sometimes inaccurately generalized.

Dual Strategy to Design New Agents Targeting Schistosoma mansoni: Advancing Phenotypic and SmCB1 Inhibitors for Improved Efficacy.

In this study, we have identified and optimized two lead structures from an in-house screening, with promising results against the parasitic flatworm Schistosoma mansoni and its target protease S. mansoni cathepsin B1 (SmCB1). Our correlation analysis highlighted the significance of physicochemical properties for the compounds' in vitro activities, resulting in a dual approach to optimize the lead structures, regarding both phenotypic effects in S. mansoni newly transformed schistosomula (NTS), adult worms, and SmCB1 inhibition. The optimized compounds from both approaches ("phenotypic" vs "SmCB1" approach) demonstrated improved efficacy against S. mansoni NTS and adult worms, with 2h from the "SmCB1" approach emerging as the most potent compound. 2h displayed nanomolar inhibition of SmCB1 (Ki = 0.050 µM) while maintaining selectivity toward human off-target cathepsins. Additionally, the greatly improved efficacy of compound 2h toward S. mansoni adults (86% dead worms at 10 µM, 68% at 1 µM, 35% at 0.1 µM) demonstrates its potential as a new therapeutic agent for schistosomiasis, underlined by its improved permeability.

Combining a transcriptomic approach and a targeted metabolomics approach for deciphering the molecular bases of compatibility phenotype in the snail Biomphalaria glabrata toward Schistosoma mansoni.

Biomphalaria glabrata is a freshwater snail and the obligatory intermediate host of Schistosoma mansoni parasite, the etiologic agent of intestinal Schistosomiasis, in South America and Caribbean. Interestingly in such host-parasite interactions, compatibility varies between populations, strains or individuals. This observed compatibility polymorphism is based on a complex molecular-matching-phenotype, the molecular bases of which have been investigated in numerous studies, notably by comparing between different strains or geographical isolates or clonal selected snail lines. Herein we propose to decipher the constitutive molecular support of this interaction in selected non-clonal resistant and susceptible snail strain originating from the same natural population from Brazil and thus having the same genetic background. Thanks to a global RNAseq transcriptomic approach on whole snail, we identified a total of 328 differentially expressed genes between resistant and susceptible phenotypes among which 129 were up-regulated and 199 down-regulated. Metabolomic studies were used to corroborate the RNAseq results. The activation of immune genes and specific metabolic pathways in resistant snails might provide them with the capacity to better respond to parasite infection.

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.

Multiplexed CRISPR-Cas9 protocol for large transgene integration into the Schistosoma mansoni genome.

Precise, on-target CRISPR-Cas9 genome editing has been shown in Schistosoma mansoni, involving both non-homology end joining and homology-directed repair pathways. Here, we present a multiplexed CRISPR-Cas9 protocol for large transgene integration into the S. mansoni genome. We describe steps for deploying multiplexed ribonucleoprotein complexes (RNPs) and donor DNA preparation. We then detail procedures for introducing RNPs into schistosome eggs by square-wave electroporation in the presence of a 5' phosphorothioate-modified double-stranded donor transgene. For complete details on the use and execution of this protocol, please refer to Ittiprasert et al. (2023).1.

Laboratory and field validation of the recombinase polymerase amplification assay targeting the Schistosoma mansoni mitochondrial minisatellite region (SmMIT-RPA) for snail xenomonitoring for schistosomiasis.

Improvements in diagnostics for schistosomiasis in both humans and snail hosts are priorities to be able to reach the World Health Organization (WHO) goal of eliminating the disease as a public health problem by 2030. In this context, molecular isothermal amplification tests, such as Recombinase Polymerase Amplification (RPA), are promising for use in endemic areas at the point-of-need for their accuracy, robustness, simplicity, and time-effectiveness. The developed recombinase polymerase amplification assay targeting the Schistosoma mansoni mitochondrial minisatellite region (SmMIT-RPA) was used to detect S. mansoni DNA from both laboratory and field Biomphalaria snails. Laboratory snails were experimentally infected and used at one, seven, and 28 days post-exposure (dpe) to 10 S. mansoni miracidia to provide samples in the early pre-patent infection stage. Field samples of Biomphalaria spp. were collected from the Mucuri Valley and Jequitinhonha Valley regions in the state of Minas Gerais, Brazil, which are endemic for S. mansoni. The sensitivity and specificity of the SmMIT-RPA assay were analysed and compared with existing loop-mediated isothermal amplification (LAMP), PCR-based methods, parasitological examination of the snails, and nucleotide sequencing. The SmMIT-RPA assay was able to detect S. mansoni DNA in the experimentally infected Biomphalaria glabrata as early as one dpe to 10 miracidia. It also detected S. mansoni infections (55.5% prevalence) in the field samples with the highest accuracy (100% sensitivity and specificity) compared with the other molecular tests used as the reference. Results from this study indicate that the SmMIT-RPA assay is a good alternative test to be used for snail xenomonitoring of S. mansoni due to its high sensitivity, accuracy, and the possibility of detecting early pre-patent infection. Its simplicity and portability also make it a suitable methodology in low-resource settings.

Subtelomeric plasticity contributes to gene family expansion in the human parasitic flatworm Schistosoma mansoni.

BACKGROUND: The genomic region that lies between the telomere and chromosome body, termed the subtelomere, is heterochromatic, repeat-rich, and frequently undergoes rearrangement. Within this region, large-scale structural changes enable gene diversification, and, as such, large multicopy gene families are often found at the subtelomere. In some parasites, genes associated with proliferation, invasion, and survival are often found in these regions, where they benefit from the subtelomere's highly plastic, rapidly changing nature. The increasing availability of complete (or near complete) parasite genomes provides an opportunity to investigate these typically poorly defined and overlooked genomic regions and potentially reveal relevant gene families necessary for the parasite's lifestyle. RESULTS: Using the latest chromosome-scale genome assembly and hallmark repeat richness observed at chromosome termini, we have identified and characterised the subtelomeres of Schistosoma mansoni, a metazoan parasitic flatworm that infects over 250 million people worldwide. Approximately 12% of the S. mansoni genome is classified as subtelomeric, and, in line with other organisms, we find these regions to be gene-poor but rich in transposable elements. We find that S. mansoni subtelomeres have undergone extensive interchromosomal recombination and that these sites disproportionately contribute to the 2.3% of the genome derived from segmental duplications. This recombination has led to the expansion of subtelomeric gene clusters containing 103 genes, including the immunomodulatory annexins and other gene families with unknown roles. The largest of these is a 49-copy plexin domain-containing protein cluster, exclusively expressed in the tegument-the tissue located at the host-parasite physical interface-of intramolluscan life stages. CONCLUSIONS: We propose that subtelomeric regions act as a genomic playground for trial-and-error of gene duplication and subsequent divergence. Owing to the importance of subtelomeric genes in other parasites, gene families implicated in this subtelomeric expansion within S. mansoni warrant further characterisation for a potential role in parasitism.

The genome and transcriptome of the snail Biomphalaria sudanica s.l.: immune gene diversification and highly polymorphic genomic regions in an important African vector of Schistosoma mansoni.

BACKGROUND: Control and elimination of schistosomiasis is an arduous task, with current strategies proving inadequate to break transmission. Exploration of genetic approaches to interrupt Schistosoma mansoni transmission, the causative agent for human intestinal schistosomiasis in sub-Saharan Africa and South America, has led to genomic research of the snail vector hosts of the genus Biomphalaria. Few complete genomic resources exist, with African Biomphalaria species being particularly underrepresented despite this being where the majority of S. mansoni infections occur. Here we generate and annotate the first genome assembly of Biomphalaria sudanica sensu lato, a species responsible for S. mansoni transmission in lake and marsh habitats of the African Rift Valley. Supported by whole-genome diversity data among five inbred lines, we describe orthologs of immune-relevant gene regions in the South American vector B. glabrata and present a bioinformatic pipeline to identify candidate novel pathogen recognition receptors (PRRs). RESULTS: De novo genome and transcriptome assembly of inbred B. sudanica originating from the shoreline of Lake Victoria (Kisumu, Kenya) resulted in a haploid genome size of ~ 944.2 Mb (6,728 fragments, N50 = 1.067 Mb), comprising 23,598 genes (BUSCO = 93.6% complete). The B. sudanica genome contains orthologues to all described immune genes/regions tied to protection against S. mansoni in B. glabrata, including the polymorphic transmembrane clusters (PTC1 and PTC2), RADres, and other loci. The B. sudanica PTC2 candidate immune genomic region contained many PRR-like genes across a much wider genomic region than has been shown in B. glabrata, as well as a large inversion between species. High levels of intra-species nucleotide diversity were seen in PTC2, as well as in regions linked to PTC1 and RADres orthologues. Immune related and putative PRR gene families were significantly over-represented in the sub-set of B. sudanica genes determined as hyperdiverse, including high extracellular diversity in transmembrane genes, which could be under pathogen-mediated balancing selection. However, no overall expansion in immunity related genes was seen in African compared to South American lineages. CONCLUSIONS: The B. sudanica genome and analyses presented here will facilitate future research in vector immune defense mechanisms against pathogens. This genomic/transcriptomic resource provides necessary data for the future development of molecular snail vector control/surveillance tools, facilitating schistosome transmission interruption mechanisms in Africa.

A single locus determines praziquantel response in Schistosoma mansoni.

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

Expression profile analysis of the transient receptor potential (TRPM) channel, a possible target of praziquantel in Schistosoma japonicum.

The WHO considers schistosomiasis, which is controlled by the mass administration of the drug praziquantel (PZQ), to be a neglected tropical disease. Despite its clinical use for over four decades, PZQ remains the only choice of chemotherapy against this disease. Regarding the previous studies that demonstrated that PZQ activates the transient receptor potential (TRP) channel in Schistosoma mansoni (Sm.TRPMPZQ), the expression profile of the ortholog of this channel gene (Smp_246790.5) in S. japonicum (EWB00_008853) (Sj.TRPMPZQ) was analyzed. The relative expression of this gene in various stages of the parasite lifecycle was analyzed by quantitative real-time reverse transcription-PCR (qRT-PCR), and the expression of Sj.TRPMPZQ was observed by immunohistochemical staining using anti-serum against the recombinant Sj.TRPMPZQ protein. qRT-PCR revealed the significantly lower mRNA expression in the snail stage in comparison to other stages (p < 0.01). The relative quantity of the Sj.TRPMPZQ expression for paired females, unpaired males, and eggs was 60%, 56%, and 68%, respectively, in comparison to paired males that showed the highest expression (p < 0.05). Interestingly, immunostaining demonstrated that Sj.TRPMPZQ is expressed in the parenchyma which contains muscle cells, neuronal cells and tegument cells in adult worms. This may support the two major effects of PZQ-worm paralysis and tegument disruption-induced by channel activation. Moreover, the channel was expressed in both the eggshell and the miracidia inside, but could not be observed in sporocyst. These results suggest that the expression of Sj.TRPMPQZ corresponds to the known sensitivity of S. japonicum to PZQ.

Fluorescent non transgenic schistosoma to decipher host-parasite phenotype compatibility.

Schistosomiasis is considered as a significant public health problem, imposing a deeper understanding of the intricate interplay between parasites and their hosts. Unfortunately, current invasive methodologies employed to study the compatibility and the parasite development impose limitations on exploring diverse strains under various environmental conditions, thereby impeding progress in the field. In this study, we demonstrate the usefulness for the trematode parasite Schistosma mansoni, leveranging a fluorescence-imaging-based approach that employs fluorescein 5-chloromethylfluorescein diacetate (CMFDA) and 5-chloromethylfluorescein diacetate (CMAC) as organism tracker for intramolluscan studies involving the host snail Biomphalaria glabrata. These probes represent key tools for qualitatively assessing snail infections with unmatched accuracy and precision. By monitoring the fluorescence of parasites within the snail vector, our method exposes an unprecedented glimpse into the host-parasite compatibility landscape. The simplicity and sensitivity of our approach render it an ideal choice for evolutionary studies, as it sheds light on the intricate mechanisms governing host-parasite interactions. Fluorescent probe-based methods play a pivotal role in characterizing factors influencing parasite development and phenotype of compatibility, paving the way for innovative, effective, and sustainable solutions to enhance our understanding host-parasite immunobiological interaction and compatibility.

Transcriptome analysis of peripheral blood of Schistosoma mansoni infected children from the Albert Nile region in Uganda reveals genes implicated in fibrosis pathology.

Over 290 million people are infected by schistosomes worldwide. Schistosomiasis control efforts focus on mass drug treatment with praziquantel (PZQ), a drug that kills the adult worm of all Schistosoma species. Nonetheless, re-infections have continued to be detected in endemic areas with individuals living in the same area presenting with varying infection intensities. Our objective was to characterize the transcriptome profiles in peripheral blood of children between 10-15 years with varying intensities of Schistosoma mansoni infection living along the Albert Nile in Uganda. RNA extracted from peripheral blood collected from 44 S. mansoni infected (34 high and 10 low by circulating anodic antigen [CAA] level) and 20 uninfected children was sequenced using Illumina NovaSeq S4 and the reads aligned to the GRCh38 human genome. Differential gene expression analysis was done using DESeq2. Principal component analysis revealed clustering of gene expression by gender when S. mansoni infected children were compared with uninfected children. In addition, we identified 14 DEGs between S. mansoni infected and uninfected individuals, 56 DEGs between children with high infection intensity and uninfected individuals, 33 DEGs between those with high infection intensity and low infection intensity and no DEGs between those with low infection and uninfected individuals. We also observed upregulation and downregulation of some DEGs that are associated with fibrosis and its regulation. These data suggest expression of fibrosis associated genes as well as genes that regulate fibrosis in S. mansoni infection. The relatively few significant DEGS observed in children with schistosomiasis suggests that chronic S. mansoni infection is a stealth infection that does not stimulate a strong immune response.

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

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

Variants of IL6, IL10, FCN2, RNASE3, IL12B and IL17B loci are associated with Schistosoma mansoni worm burden in the Albert Nile region of Uganda.

BACKGROUND: Individuals genetically susceptible to high schistosomiasis worm burden may contribute disproportionately to transmission and could be prioritized for control. Identifying genes involved may guide development of therapy. METHODOLOGY/PRINCIPAL FINDINGS: A cohort of 606 children aged 10-15 years were recruited in the Albert Nile region of Uganda and assessed for Schistosoma mansoni worm burden using the Up-Converting Particle Lateral Flow (UCP-LF) test detecting circulating anodic antigen (CAA), point-of-care Circulating Cathodic Antigen (POC-CCA) and Kato-Katz tests. Whole genome genotyping was conducted on 326 children comprising the top and bottom 25% of worm burden. Linear models were fitted to identify variants associated with worm burden in preselected candidate genes. Expression quantitative trait locus (eQTL) analysis was conducted for candidate genes with UCP-LF worm burden included as a covariate. Single Nucleotide Polymorphism loci associated with UCP-LF CAA included IL6 rs2066992 (OR = 0.43, p = 0.0006) and rs7793163 (OR = 2.0, p = 0.0007); IL21 SNP kgp513476 (OR 1.79, p = 0.0025) and IL17B SNP kgp708159 (OR = 0.35, p = 0.0028). A haplotype in the IL10 locus was associated with lower worm burden (OR = 0.53, p = 0.015) and overlapped SNPs rs1800896, rs1800871 and rs1800872. Significant haplotypes (p<0.05, overlapping significant SNP) associated with worm burden were observed in IL6 and the Th17 pathway IL12B and IL17B genes. There were significant eQTL in the IL6, IL5, IL21, IL25 and IFNG regions. CONCLUSIONS: Variants associated with S. mansoni worm burden were in IL6, FCN2, RNASE3, IL10, IL12B and IL17B gene loci. However only eQTL associations remained significant after Bonferroni correction. In summary, immune balance, pathogen recognition and Th17 pathways may play a role in modulating Schistosoma worm burden. Individuals carrying risk variants may be targeted first in allocation of control efforts to reduce the burden of schistosomiasis in the community.

Identification of Genomic Regions Implicated in Susceptibility to Schistosoma mansoni Infection in a Murine Backcross Genetic Model.

Only a small number of infected people are highly susceptible to schistosomiasis, showing high levels of infection or severe liver fibrosis. The susceptibility to schistosome infection is influenced by genetic background. To assess the genetic basis of susceptibility and identify the chromosomal regions involved, a backcross strategy was employed to generate high variation in schistosomiasis susceptibility. This strategy involved crossing the resistant C57BL/6J mouse strain with the susceptible CBA/2J strain. The resulting F1 females (C57BL/6J × CBA/2J) were then backcrossed with CBA/2J males to generate the backcross (BX) cohort. The BX mice exhibited a range of phenotypes, with disease severity varying from mild to severe disease, lacking a fully resistant group. We observed four levels of infection intensity using cluster and principal component analyses and K-means based on parasitological, pathological, and immunological trait measurements. The mice were genotyped with 961 informative SNPs, leading to the identification of 19 new quantitative trait loci (QTL) associated with parasite burden, liver lesions, white blood cell populations, and antibody responses. Two QTLs located on chromosomes 15 and 18 were linked to the number of granulomas, liver lesions, and IgM levels. The corresponding syntenic human regions are located in chromosomes 8 and 18. None of the significant QTLs had been reported previously.

Revisiting Schistosoma mansoni Micro-Exon Gene (MEG) Protein Family: A Tour into Conserved Motifs and Annotation.

Genome sequencing of the human parasite Schistosoma mansoni revealed an interesting gene superfamily, called micro-exon gene (meg), that encodes secreted MEG proteins. The genes are composed of short exons (3-81 base pairs) regularly interspersed with long introns (up to 5 kbp). This article recollects 35 S. mansoni specific meg genes that are distributed over 7 autosomes and one pair of sex chromosomes and that code for at least 87 verified MEG proteins. We used various bioinformatics tools to produce an optimal alignment and propose a phylogenetic analysis. This work highlighted intriguing conserved patterns/motifs in the sequences of the highly variable MEG proteins. Based on the analyses, we were able to classify the verified MEG proteins into two subfamilies and to hypothesize their duplication and colonization of all the chromosomes. Together with motif identification, we also proposed to revisit MEGs' common names and annotation in order to avoid duplication, to help the reproducibility of research results and to avoid possible misunderstandings.