Long non-coding RNAs are essential for Schistosoma mansoni pairing-dependent adult worm homeostasis and fertility.

The trematode parasite Schistosoma mansoni causes schistosomiasis, which affects over 200 million people worldwide. Schistosomes are dioecious, with egg laying depending on the females' obligatory pairing with males. Long non-coding RNAs (lncRNAs) are transcripts longer than 200 nucleotides with low or no protein-coding potential that have been involved in other species with reproduction, stem cell maintenance, and drug resistance. In S. mansoni, we recently showed that the knockdown of one lncRNA affects the pairing status of these parasites. Here, we re-analyzed public RNA-Seq data from paired and unpaired adult male and female worms and their gonads, obtained from mixed-sex or single-sex cercariae infections, and found thousands of differentially expressed pairing-dependent lncRNAs among the 23 biological samples that were compared. The expression levels of selected lncRNAs were validated by RT-qPCR using an in vitro unpairing model. In addition, the in vitro silencing of three selected lncRNAs showed that knockdown of these pairing-dependent lncRNAs reduced cell proliferation in adult worms and their gonads, and are essential for female vitellaria maintenance, reproduction, and/or egg development. Remarkably, in vivo silencing of each of the three selected lncRNAs significantly reduced worm burden in infected mice by 26 to 35%. Whole mount in situ hybridization experiments showed that these pairing-dependent lncRNAs are expressed in reproductive tissues. These results show that lncRNAs are key components intervening in S. mansoni adult worm homeostasis, which affects pairing status and survival in the mammalian host, thus presenting great potential as new therapeutic target candidates.

The Aedes aegypti peritrophic matrix controls arbovirus vector competence through HPx1, a heme-induced peroxidase.

Aedes aegypti mosquitoes are the main vectors of arboviruses. The peritrophic matrix (PM) is an extracellular layer that surrounds the blood bolus. It acts as an immune barrier that prevents direct contact of bacteria with midgut epithelial cells during blood digestion. Here, we describe a heme-dependent peroxidase, hereafter referred to as heme peroxidase 1 (HPx1). HPx1 promotes PM assembly and antioxidant ability, modulating vector competence. Mechanistically, the heme presence in a blood meal induces HPx1 transcriptional activation mediated by the E75 transcription factor. HPx1 knockdown increases midgut reactive oxygen species (ROS) production by the DUOX NADPH oxidase. Elevated ROS levels reduce microbiota growth while enhancing epithelial mitosis, a response to tissue damage. However, simultaneous HPx1 and DUOX silencing was not able to rescue bacterial population growth, as explained by increased expression of antimicrobial peptides (AMPs), which occurred only after double knockdown. This result revealed hierarchical activation of ROS and AMPs to control microbiota. HPx1 knockdown produced a 100-fold decrease in Zika and dengue 2 midgut infection, demonstrating the essential role of the mosquito PM in the modulation of arbovirus vector competence. Our data show that the PM connects blood digestion to midgut immunological sensing of the microbiota and viral infections.

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.

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.

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.

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.

Long non-coding RNA levels can be modulated by 5-azacytidine in Schistosoma mansoni.

Schistosoma mansoni is a flatworm that causes schistosomiasis, a neglected tropical disease that affects more than 200 million people worldwide. There is only one drug indicated for treatment, praziquantel, which may lead to parasite resistance emergence. The ribonucleoside analogue 5-azacytidine (5-AzaC) is an epigenetic drug that inhibits S. mansoni oviposition and ovarian development through interference with parasite transcription, translation and stem cell activities. Therefore, studying the downstream pathways affected by 5-AzaC in S. mansoni may contribute to the discovery of new drug targets. Long non-coding RNAs (lncRNAs) are transcripts longer than 200 nucleotides with low or no protein coding potential that have been involved in reproduction, stem cell maintenance and drug resistance. We have recently published a catalog of lncRNAs expressed in S. mansoni life-cycle stages, tissues and single cells. However, it remains largely unknown if lncRNAs are responsive to epigenetic drugs in parasites. Here, we show by RNA-Seq re-analyses that hundreds of lncRNAs are differentially expressed after in vitro 5-AzaC treatment of S. mansoni females, including intergenic, antisense and sense lncRNAs. Many of these lncRNAs belong to co-expression network modules related to male metabolism and are also differentially expressed in unpaired compared with paired females and ovaries. Half of these lncRNAs possess histone marks at their genomic loci, indicating regulation by histone modification. Among a selected set of 8 lncRNAs, half of them were validated by RT-qPCR as differentially expressed in females, and some of them also in males. Interestingly, these lncRNAs are also expressed in other life-cycle stages. This study demonstrates that many lncRNAs potentially involved with S. mansoni reproductive biology are modulated by 5-AzaC and sheds light on the relevance of exploring lncRNAs in response to drug treatments in parasites.

In vitro activity of aryl-thiazole derivatives against Schistosoma mansoni schistosomula and adult worms.

Schistosomiasis is caused by a trematode of the genus Schistosoma and affects over 200 million people worldwide. The only drug recommended by the World Health Organization for treatment and control of schistosomiasis is praziquantel. Development of new drugs is therefore of great importance. Thiazoles are regarded as privileged structures with a broad spectrum of activities and are potential sources of new drug prototypes, since they can act through interactions with DNA and inhibition of DNA synthesis. In this context, we report the synthesis of a series of thiazole derivatives and their in vitro schistosomicidal activity by testing eight molecules (NJ03-08; NJ11-12) containing thiazole structures. Parameters such as motility and mortality, egg laying, pairing and parasite viability by ATP quantification, which were influenced by these compounds, were evaluated during the assays. Scanning electron microscopy (SEM) was utilized for evaluation of morphological changes in the tegument. Schistosomula and adult worms were treated in vitro with different concentrations (6.25 to 50 µM) of the thiazoles for up to 5 and 3 days, respectively. After in vitro treatment for five days with 6.25 µM NJ05 or NJ07 separately, we observed a decrease of 30% in schistosomula viability, whilst treatment with NJ05+NJ07 lead to a reduction of 75% in viability measured by ATP quantitation and propidium iodide labeling. Adult worms' treatment with 50 µM NJ05, NJ07 or NJ05 + NJ07 showed decreased motility to 30-50% compared with controls. Compound NJ05 was more effective than NJ07, and adult worm viability after three days was reduced to 25% in parasites treated with 50 µM NJ05, compared with a viability reduction to 40% with 50 µM NJ07. SEM analysis showed severe alterations in adult worms with formation of bulges and blisters throughout the dorsal region of parasites treated with NJ05 or NJ07. Oviposition was extremely affected by treatment with the NJ series compounds; at concentrations of 25 µM and 50 µM, oviposition reached almost zero with NJ05, NJ07 or NJ05 + NJ07 already at day one. Tested genes involved in egg biosynthesis were all confirmed by qPCR as downregulated in females treated with 25 µM NJ05 for 2 days, with a significant reduction in expression of p14, Tyrosinase 2, p48 and fs800. NJ05, NJ07 or NJ05+NJ07 treatment of HEK293 (human embryonic cell line) and HES (human epithelial cell line) showed EC50 in the range of 18.42 to 145.20 µM. Overall, our results demonstrate that those molecules are suitable targets for further development into new drugs for schistosomiasis treatment, although progress is needed to lessen the cytotoxic effects on human cells. According to the present study, thiazole derivatives have schistosomicidal activities and may be part of a possible new arsenal of compounds against schistosomiasis.

Weighted Gene Co-Expression Analyses Point to Long Non-Coding RNA Hub Genes at Different Schistosoma mansoni Life-Cycle Stages.

Long non-coding RNAs (lncRNAs) (>200 nt) are expressed at levels lower than those of the protein-coding mRNAs, and in all eukaryotic model species where they have been characterized, they are transcribed from thousands of different genomic loci. In humans, some four dozen lncRNAs have been studied in detail, and they have been shown to play important roles in transcriptional regulation, acting in conjunction with transcription factors and epigenetic marks to modulate the tissue-type specific programs of transcriptional gene activation and repression. In Schistosoma mansoni, around 10,000 lncRNAs have been identified in previous works. However, the limited number of RNA-sequencing (RNA-seq) libraries that had been previously assessed, together with the use of old and incomplete versions of the S. mansoni genome and protein-coding transcriptome annotations, have hampered the identification of all lncRNAs expressed in the parasite. Here we have used 633 publicly available S. mansoni RNA-seq libraries from whole worms at different stages (n = 121), from isolated tissues (n = 24), from cell-populations (n = 81), and from single-cells (n = 407). We have assembled a set of 16,583 lncRNA transcripts originated from 10,024 genes, of which 11,022 are novel S. mansoni lncRNA transcripts, whereas the remaining 5,561 transcripts comprise 120 lncRNAs that are identical to and 5,441 lncRNAs that have gene overlap with S. mansoni lncRNAs already reported in previous works. Most importantly, our more stringent assembly and filtering pipeline has identified and removed a set of 4,293 lncRNA transcripts from previous publications that were in fact derived from partially processed mRNAs with intron retention. We have used weighted gene co-expression network analyses and identified 15 different gene co-expression modules. Each parasite life-cycle stage has at least one highly correlated gene co-expression module, and each module is comprised of hundreds to thousands lncRNAs and mRNAs having correlated co-expression patterns at different stages. Inspection of the top most highly connected genes within the modules' networks has shown that different lncRNAs are hub genes at different life-cycle stages, being among the most promising candidate lncRNAs to be further explored for functional characterization.