Atlas of Schistosoma mansoni long non-coding RNAs and their expression correlation to protein-coding genes

Expression of Transposable Elements throughout the Fasciola hepatica Trematode Life Cycle

BACKGROUND

Transposable elements (TEs) are major components of eukaryotic genomes. The extensive body of evidence suggests that although they were once considered "genomic parasites", transposons and their transcripts perform specific functions, such as regulation of early embryo development. Understanding the role of TEs in such parasites as trematodes is becoming critically important. Fasciola hepatica, a parasite affecting humans and livestock, undergoes a complex life cycle in diverse environments and hosts, and knowledge about its life cycle regulation is scarce so far.

METHODS

We summarized the data regarding the repetitive elements in F. hepatica and conducted bulk RNA-seq analysis across its life cycle stages. TE expression profiles were analyzed, focusing on differential expression and potential homology with previously described long non-coding RNAs (lncRNAs).

RESULTS

Differential expression analysis revealed stage-specific TE transcription patterns, notably peaking during egg and metacercariae stages. Some TEs showed homology with known lncRNAs and contained putative transcription factor binding sites. Interestingly, TE transcription levels were highest in eggs and metacercariae compared to adults, suggesting regulatory roles in trematode life cycle transitions.

CONCLUSIONS

These findings suggest that TEs may play roles in regulating trematode life cycle transitions. Moreover, TE homology with lncRNAs underscores their significance in gene regulation.

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.

Functional Intricacy and Symmetry of Long Non-Coding RNAs in Parasitic Infections

RNAs are a class of molecules and the majority in eukaryotes are arbitrarily termed non- coding transcripts which are broadly classified as short and long non-coding RNAs. Recently, knowledge of the identification and functions of long non-coding RNAs have continued to accumulate and they are being recognized as important molecules that regulate parasite-host interface, parasite differentiation, host responses, and disease progression. Herein, we present and integrate the functions of host and parasite long non-coding RNAs during infections within the context of epigenetic re-programming and molecular crosstalk in the course of host-parasite interactions. Also, the modular range of parasite and host long non-coding RNAs in coordinated parasite developmental changes and host immune dynamic landscapes are discussed. We equally canvass the prospects of long non-coding RNAs in disease diagnosis and prognosis. Hindsight and suggestions are offered with the aim that it will bolster our understanding for future works on host and parasite long non-coding RNAs.

In Silico Analysis of Common Long Noncoding RNAs in Schistosoma mansoni and Schistosoma haematobium

Background

Schistosomiasis caused by Schistosoma parasites is one of the most common parasitic infections worldwide. Genetic regulation of the genus Schistosoma, which has different developmental stages throughout its life, is quite complex. In these parasites, thousands of long noncoding RNAs (lncRNAs) estimated to be functional were identified. Identifying the transcripts expressed in common and detecting their functions for better understanding of the role of these lncRNAs require a comparative study.

Methods

Assembled RNA-seq datasets belonging to S. mansoni and S. haematobium were obtained from the National Center for Biotechnology. A basic local alignment search tool (BLASTN) analysis was conducted against previously constructed lncRNA library to identify the common lncRNAs between two species. LncRNAs target genes and their gene ontology annotation was performed.

Results

In S. mansoni and S. haematobium, 5132 and 3589 lncRNA transcripts were detected, respectively. These two species had 694 lncRNAs in common. A significant number of lncRNAs was determined to be transcribed from sex chromosomes. The frequently expressed lncRNAs appear to be involved in metabolic and biological regulation processes.

Conclusions

These two species share similar lncRNAs; thus, this finding is a clue that they might have similar functions. In sexual development, they especially might play important roles. Our results will provide important clues to further studies about interactions between human hosts and parasites and the infection mechanisms of Schistosoma parasites.

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.

RNA sequencing analysis of altered expression of long noncoding RNAs associated with Schistosoma japonicum infection in the murine liver and spleen

BACKGROUND

Schistosomiasis is a chronic, debilitating infectious disease caused by members of the genus Schistosoma. Previous findings have suggested a relationship between infection with Schistosoma spp. and alterations in the liver and spleen of infected animals. Recent reports have shown the regulatory role of noncoding RNAs, such as long noncoding RNAs (lncRNAs), in different biological processes. However, little is known about the role of lncRNAs in the mouse liver and spleen during Schistosoma japonicum infection.

METHODS

In this study, we identified and investigated lncRNAs using standard RNA sequencing (RNA-Seq). The biological functions of the altered expression of lncRNAs and their target genes were predicted using bioinformatics. Ten dysregulated lncRNAs were selected randomly and validated in reverse transcription-quantitative real-time polymerase chain reaction (RT-qPCR) experiments.

RESULTS

Our study identified 29,845 and 33,788 lncRNAs from the liver and spleen, respectively, of which 212 were novel lncRNAs. We observed that 759 and 789 of the lncRNAs were differentially expressed in the respective organs. The RT-qPCR results correlated well with the sequencing data. In the liver, 657 differentially expressed lncRNAs were predicted to target 2548 protein-coding genes, whereas in the spleen 660 differentially expressed lncRNAs were predicted to target 2673 protein-coding genes. Moreover, functional annotation showed that the target genes of the differentially expressed lncRNAs were associated with cellular processes, metabolic processes, and binding, and were significantly enriched in metabolic pathways, the cell cycle, ubiquitin-mediated proteolysis, and pathways in cancer.

CONCLUSIONS

Our study showed that numerous lncRNAs were differentially expressed in S. japonicum-infected liver and spleen compared to control liver and spleen; this suggested that lncRNAs may be involved in pathogenesis in the liver and spleen during S. japonicum infection.

Molecular signatures of the rediae, cercariae and adult stages in the complex life cycles of parasitic flatworms (Digenea: Psilostomatidae)

Background

Parasitic flatworms (Trematoda: Digenea) represent one of the most remarkable examples of drastic morphological diversity among the stages within a life cycle. Which genes are responsible for extreme differences in anatomy, physiology, behavior, and ecology among the stages? Here we report a comparative transcriptomic analysis of parthenogenetic and amphimictic generations in two evolutionary informative species of Digenea belonging to the family Psilostomatidae.

Methods

In this study the transcriptomes of rediae, cercariae and adult worm stages of Psilotrema simillimum and Sphaeridiotrema pseudoglobulus, were sequenced and analyzed. High-quality transcriptomes were generated, and the reference sets of protein-coding genes were used for differential expression analysis in order to identify stage-specific genes. Comparative analysis of gene sets, their expression dynamics and Gene Ontology enrichment analysis were performed for three life stages within each species and between the two species.

Results

Reference transcriptomes for P. simillimum and S. pseudoglobulus include 21,433 and 46,424 sequences, respectively. Among 14,051 orthologous groups (OGs), 1354 are common and specific for two analyzed psilostomatid species, whereas 13 and 43 OGs were unique for P. simillimum and S. pseudoglobulus, respectively. In contrast to P. simillimum, where more than 60% of analyzed genes were active in the redia, cercaria and adult worm stages, in S. pseudoglobulus less than 40% of genes had such a ubiquitous expression pattern. In general, 7805 (36.41%) and 30,622 (65.96%) of genes were preferentially expressed in one of the analyzed stages of P. simillimum and S. pseudoglobulus, respectively. In both species 12 clusters of co-expressed genes were identified, and more than a half of the genes belonging to the reference sets were included into these clusters. Functional specialization of the life cycle stages was clearly supported by Gene Ontology enrichment analysis.

Conclusions

During the life cycles of the two species studied, most of the genes change their expression levels considerably, consequently the molecular signature of a stage is not only a unique set of expressed genes, but also the specific levels of their expression. Our results indicate unexpectedly high level of plasticity in gene regulation between closely related species. Transcriptomes of P. simillimum and S. pseudoglobulus provide high quality reference resource for future evolutionary studies and comparative analyses.