Identification of long noncoding RNAs in Schistosoma mansoni and Schistosoma japonicum

Dynamic profiles of lncRNAs reveal a functional natural antisense RNA that regulates the development of Schistosoma japonicum

Schistosomes are flatworm parasites that undergo a complex life cycle involving two hosts. The regulation of the parasite's developmental processes relies on both coding RNAs and non-coding RNAs. However, the roles of non-coding RNAs, including long non-coding RNAs (lncRNAs) in schistosomes remain largely unexplored. Here we conduct advanced RNA sequencing on male and female S. japonicum during their pairing and reproductive development, resulting in the identification of nearly 8,000 lncRNAs. This extensive dataset enables us to construct a comprehensive co-expression network of lncRNAs and mRNAs, shedding light on their interactions during the crucial reproductive stages within the mammalian host. Importantly, we have also revealed a specific lncRNA, LNC3385, which appears to play a critical role in the survival and reproduction of the parasite. These findings not only enhance our understanding of the dynamic nature of lncRNAs during the reproductive phase of schistosomes but also highlight LNC3385 as a potential therapeutic target for combating schistosomiasis.

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.

Discovery of long non-coding RNAs in the liver fluke, Fasciola hepatica

Long non-coding (lnc)RNAs are a class of eukaryotic RNA that do not code for protein and are linked with transcriptional regulation, amongst a myriad of other functions. Using a custom in silico pipeline we have identified 6,436 putative lncRNA transcripts in the liver fluke parasite, Fasciola hepatica, none of which are conserved with those previously described from Schistosoma mansoni. F. hepatica lncRNAs were distinct from F. hepatica mRNAs in transcript length, coding probability, exon/intron composition, expression patterns, and genome distribution. RNA-Seq and digital droplet PCR measurements demonstrated developmentally regulated expression of lncRNAs between intra-mammalian life stages; a similar proportion of lncRNAs (14.2%) and mRNAs (12.8%) were differentially expressed (p<0.001), supporting a functional role for lncRNAs in F. hepatica life stages. While most lncRNAs (81%) were intergenic, we identified some that overlapped protein coding loci in antisense (13%) or intronic (6%) configurations. We found no unequivocal evidence for correlated developmental expression within positionally correlated lncRNA:mRNA pairs, but global co-expression analysis identified five lncRNA that were inversely co-regulated with 89 mRNAs, including a large number of functionally essential proteases. The presence of micro (mi)RNA binding sites in 3135 lncRNAs indicates the potential for miRNA-based post-transcriptional regulation of lncRNA, and/or their function as competing endogenous (ce)RNAs. The same annotation pipeline identified 24,141 putative lncRNAs in F. gigantica. This first description of lncRNAs in F. hepatica provides an avenue to future functional and comparative genomics studies that will provide a new perspective on a poorly understood aspect of parasite biology.

Roles of microRNAs and Long Non-Coding RNAs Encoded by Parasitic Helminths in Human Carcinogenesis

Infectious agents such as viruses, bacteria, and parasites can lead to cancer development. Infection with the helminthic parasite Schistosoma haematobium can cause cancer of the urinary bladder in humans, and infection with the parasites Clonorchis sinensis and Opisthorchis viverrini can promote cholangiocarcinoma. These three pathogens have been categorized as “group 1: carcinogenic to humans” by the International Agency for Research on Cancer (IARC). Additionally, the parasite Schistosoma japonicum has been associated with liver and colorectal cancer and classified as “group 2B: possibly carcinogenic to humans”. These parasites express regulatory non-coding RNAs as microRNAs (miRNAs) and long non-coding RNAs (lncRNAs), which modulate genic expression in different biological processes. In this review, we discuss the potential roles of miRNAS and lncRNAs encoded by helminthic parasites that are classified by the IARC as carcinogenic and possibly carcinogenic to humans. The miRNAs of these parasites may be involved in carcinogenesis by modulating the biological functions of the pathogen and the host and by altering microenvironments prone to tumor growth. miRNAs were identified in different host fluids. Additionally, some miRNAs showed direct antitumoral effects. Together, these miRNAs show potential for use in future therapeutic and diagnostic applications. LncRNAs have been less studied in these parasites, and their biological effects in the parasite–host interaction are largely unknown.

Genome-wide identification of circular RNAs in adult Schistosoma japonicum

Circular RNAs (circRNAs) are a class of novel, widespread, covalently closed RNAs that have played an essential role in animal gene regulation. To systematically explore circRNAs in the blood fluke Schistosoma japonicum, we performed RNA sequencing and bioinformatics analysis, and found that hundreds of circRNAs showed gender-associated expression. Among these identified circRNAs, more than 77.54% and 74.73% were putatively derived from the exon region of the genome and some circRNAs showed gender-associated expressions. The functional prediction of circRNAs (circ_003826 and circ_004690) showed potential binding sites and possibly acted as the sponge to regulate microRNAs (miRNAs) sja-miR-1, sja-miR-133 and sja-miR-3504. Altogether, these findings demonstrated that S. japonicum also contains circRNAs, which may have potential regulatory roles during schistosome development.

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.

Detection of Schistosoma mansoni long non-coding RNAs in the infected C57BL/6 mouse liver

Long non-coding RNAs (lncRNAs) perform several types of regulatory functions and have been recently explored in the genus Schistosoma. Although sequencing and bioinformatics approaches have demonstrated the presence of hundreds of lncRNAs and microRNAs (miRNAs) in this genus, information regarding their abundance, characteristics, and potential functions linked to Schistosoma mansoni biology and parasite-host interaction is limited. Our objectives in the present study were to verify whether 15 previously identified S. mansoni lncRNAs are detectable in the host liver. In addition, we assess whether these lncRNAs are present in the S. mansoni infective form and the stages inside the definitive host. The detection of these 15 S. mansoni lncRNAs and a long terminal repeat (LTR) retrotransposon Saci 4 was performed in the eggs, cercariae, and 3.5-h schistosomula. All lncRNAs were found to be expressed in these stages; some of the lncRNAs were found in the livers of the infected C57BL/6 mice. In conclusion, S. mansoni lncRNAs were detected in host livers and quantified. Furthermore, many of the lncRNAs analyzed showed differential expression in the larval stages, indicating that they play a stage-specific regulatory role.

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.

Post-genomic progress in helminth parasitology

Helminth parasitology is an important discipline, which poses often unique technical challenges. One challenge is that helminth parasites, particularly those in humans, are often difficult to obtain alive and in sufficient quantities for study; another is the challenge of studying these organisms in vitro – no helminth parasite life cycle has been fully recapitulated outside of a host. Arguably, the key issue retarding progress in helminth parasitology has been a lack of experimental tools and resources, certainly relative to the riches that have driven many parasitologists to adopt free-living model organisms as surrogate systems. In response to these needs, the past 10–12 years have seen the beginnings of helminth parasitology's journey into the ‘omics’ era, with the release of abundant sequencing resources, and the functional genomics tools with which to test biological hypotheses. To reflect this progress, the 2019 Autumn Symposium of the British Society for Parasitology was held in Queen's University Belfast on the topic of ‘post-genomic progress in helminth parasitology’. This issue presents examples of the current state of play in the field, while this editorial summarizes how genomic datasets and functional genomic tools have stimulated impressive recent progress in our understanding of parasite biology.

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.

Step-by-Step Bioinformatics Analysis of Schistosoma mansoni Long Non-coding RNA Sequences

In the last few years, long non-coding RNAs (lncRNAs) have been widely studied in humans, and their relevance for physiological and pathological conditions has been demonstrated. In parasites, there are only a few works, such as in Plasmodium falciparum, where it was shown that an lncRNA regulates the expression of a gene associated with immune system evasion, also indicating the relevance of understanding the role of this class of RNAs in parasites. In Schistosoma mansoni, in the last 2 years, there were four published articles related to the annotation of lncRNAs in different life cycle stages using RNA-Seq libraries. In order to make this process of lncRNA identification and annotation more accessible to biologists with no bioinformatics training, considering the growing number of S. mansoni RNA-Seq libraries publicly available from different sources, such as ovary tissues from bi-sex and single-sex infections, and the potential of lncRNAs as therapeutic targets, we provide this step-by-step protocol of lncRNA identification and quantification. This guide includes the download of RNA-Seq libraries from a public database and reads processing and mapping against the genome, transcript reconstruction, novel lncRNA identification, transcripts expression level determination, and the identification of differentially expressed lncRNAs.

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.

High-Throughput RNA Sequencing Reveals NDUFC2-AS lncRNA Promotes Adipogenic Differentiation in Chinese Buffalo (Bubalus bubalis L)

The buffalo (Bubalus bubalis L.) is prevalent in China and the increasing demand for meat production has changed its role from being a beast of burden to a meat source. The low fat deposition level has become one of the main barriers for its use in meat production. It is urgent to reveal factors involved in fat deposition in buffalo. This study performed RNA sequencing to investigate both long noncoding RNAs (lncRNAs) and mRNAs of adipose tissues in young and adult buffalos. A total of 124 lncRNAs and 2008 mRNAs showed differential expression patterns between young and adult samples. Coexpression analysis and functional enrichment revealed 585 mRNA–lncRNA pairs with potential function in fat deposition. After validation by qRT-PCR, we focused on a lncRNA transcribed from the ubiquinone oxidoreductase subunit C2 (NDUFC2) antisense (AS) strand which showed high correlation with thyroid hormone responsive protein (THRSP). NDUFC2-AS lncRNA is highly expressed in adipose tissue and maturation adipocytes and mainly exists in the nucleus. Functional assays demonstrated that NDUFC2-AS lncRNA promotes adipogenic differentiation by upregulating the expression levels of THRSP and CCAAT enhancer binding protein alpha (C/EBPα) in buffalo. These results indicate that NDUFC2-AS lncRNA promotes fat deposition in buffalo.

Long Non-Coding RNAs in the Regulation of Gene Expression: Physiology and Disease

The identification of RNAs that are not translated into proteins was an important breakthrough, defining the diversity of molecules involved in eukaryotic regulation of gene expression. These non-coding RNAs can be divided into two main classes according to their length: short non-coding RNAs, such as microRNAs (miRNAs), and long non-coding RNAs (lncRNAs). The lncRNAs in association with other molecules can coordinate several physiological processes and their dysfunction may impact in several pathologies, including cancer and infectious diseases. They can control the flux of genetic information, such as chromosome structure modulation, transcription, splicing, messenger RNA (mRNA) stability, mRNA availability, and post-translational modifications. Long non-coding RNAs present interaction domains for DNA, mRNAs, miRNAs, and proteins, depending on both sequence and secondary structure. The advent of new generation sequencing has provided evidences of putative lncRNAs existence; however, the analysis of transcriptomes for their functional characterization remains a challenge. Here, we review some important aspects of lncRNA biology, focusing on their role as regulatory elements in gene expression modulation during physiological and disease processes, with implications in host and pathogens physiology, and their role in immune response modulation.