Genome-wide identification and functional annotation of Plasmodium falciparum long noncoding RNAs from RNA-seq data

Novel insights into the role of long non-coding RNA in the human malaria parasite, Plasmodium falciparum

The complex life cycle of Plasmodium falciparum requires coordinated gene expression regulation to allow host cell invasion, transmission, and immune evasion. Increasing evidence now suggests a major role for epigenetic mechanisms in gene expression in the parasite. In eukaryotes, many lncRNAs have been identified to be pivotal regulators of genome structure and gene expression. To investigate the regulatory roles of lncRNAs in P. falciparum we explore the intergenic lncRNA distribution in nuclear and cytoplasmic subcellular locations. Using nascent RNA expression profiles, we identify a total of 1768 lncRNAs, of which 718 (~41%) are novels in P. falciparum. The subcellular localization and stage-specific expression of several putative lncRNAs are validated using RNA-FISH. Additionally, the genome-wide occupancy of several candidate nuclear lncRNAs is explored using ChIRP. The results reveal that lncRNA occupancy sites are focal and sequence-specific with a particular enrichment for several parasite-specific gene families, including those involved in pathogenesis and sexual differentiation. Genomic and phenotypic analysis of one specific lncRNA demonstrate its importance in sexual differentiation and reproduction. Our findings bring a new level of insight into the role of lncRNAs in pathogenicity, gene regulation and sexual differentiation, opening new avenues for targeted therapeutic strategies against the deadly malaria parasite.

The role of long noncoding RNAs in malaria parasites

The human malaria parasites, including Plasmodium falciparum, persist as a major cause of global morbidity and mortality. The recent stalling of progress toward malaria elimination substantiates a need for novel interventions. Controlled gene expression is central to the parasite's numerous life cycle transformations and adaptation. With few specific transcription factors (TFs) identified, crucial roles for chromatin states and epigenetics in parasite transcription have become evident. Although many chromatin-modifying enzymes are known, less is known about which factors mediate their impacts on transcriptional variation. Like those of higher eukaryotes, long noncoding RNAs (lncRNAs) have recently been shown to have integral roles in parasite gene regulation. This review aims to summarize recent developments and key findings on the role of lncRNAs in P. falciparum.

Cryptosporidium parvum hijacks a host's long noncoding RNA U90926 to evade intestinal epithelial cell-autonomous antiparasitic defense

Cryptosporidium is a zoonotic apicomplexan parasite that infects the gastrointestinal epithelium and other mucosal surfaces in humans. It is an important opportunistic pathogen in AIDS patients and a leading cause of infectious diarrhea and diarrheal-related death in children worldwide. The intestinal epithelial cells provide the first line of defense against Cryptosporidium infection and play a central role in activating and regulating the host's antiparasitic response. Increasing evidence suggests that long noncoding RNAs (lncRNAs) participate in host-pathogen interactions and play a regulatory role in the pathogenesis of diseases but the underlying molecular mechanisms are not fully understood. We previously identified a panel of host lncRNAs that are upregulated in murine intestinal epithelial cells following Cryptosporidium infection, including U90926. We demonstrate here that U90926 is acting in a pro-parasitic manner in regulating intestinal epithelial cell-autonomous antiparasitic defense. Inhibition of U90926 resulted in a decreased infection burden of the parasite while overexpression of U90926 showed an increase in infection burden in cultured murine intestinal epithelial cells. Induction of U90926 suppressed transcription of epithelial defense genes involved in controlling Cryptosporidium infection through epigenetic mechanisms. Specifically, transcription of Aebp1, which encodes the Aebp1 protein, a potent modulator of inflammation and NF-κB signaling, was suppressed by U90926. Gain- or loss-of-function of Aebp1 in the host's epithelial cells caused reciprocal alterations in the infection burden of the parasite. Interestingly, Cryptosporidium carries the Cryptosporidium virus 1 (CSpV1), a double-stranded (ds) RNA virus coding two dsRNA fragments, CSpV1-dsRdRp and CSpV1-dsCA. Both CSpV1-dsRdRp and CSpV1-dsCA can be delivered into infected cells as previously reported. We found that cells transfected with in vitro transcribed CSpV1-dsCA or CSpV1-dsRdRp displayed an increased level of U90926, suggesting that CSpV1 is involved in the upregulation of U90926 during Cryptosporidium infection. Our study highlights a new strategy by Cryptosporidium to hijack a host lncRNA to suppress epithelial cell-autonomous antiparasitic defense and allow for a robust infection.

A manually curated annotation characterises genomic features of P. falciparum lncRNAs

BACKGROUND

Important regulation occurs at the level of transcription in Plasmodium falciparum and growing evidence suggests that these apicomplexan parasites have complex regulatory networks. Recent studies implicate long noncoding RNAs (lncRNAs) as transcriptional regulators in P. falciparum. However, due to limited research and the lack of necessary experimental tools, our understanding of their role in the malaria-causing parasite remains largely unelucidated. In this work, we address one of these limitations, the lack of an updated and improved lncRNA annotation in P. falciparum.

RESULTS

We generated long-read RNA sequencing data and integrated information extracted and curated from multiple sources to manually annotate lncRNAs. We identified 1119 novel lncRNAs and validated and refined 1250 existing annotations. Utilising the collated datasets, we generated evidence-based ranking scores for each annotation and characterised the distinct genomic contexts and features of P. falciparum lncRNAs. Certain features indicated subsets with potential biological significance such as 25 lncRNAs containing multiple introns, 335 lncRNAs lacking mutations in piggyBac mutagenic studies and lncRNAs associated with specific biologic processes including two new types of lncRNAs found proximal to var genes.

CONCLUSIONS

The insights and the annotation presented in this study will serve as valuable tools for researchers seeking to understand the role of lncRNAs in parasite biology through both bioinformatics and experimental approaches.

LINC01526 Promotes Proliferation and Metastasis of Gastric Cancer by Interacting with TARBP2 to Induce GNG7 mRNA Decay

Simple Summary

Many long noncoding RNAs play an important role in gastric cancer progression. In this study, we focused on LINC01526. Through expression and functional analyses, we obtained a preliminary understanding of the pro-cancer role of LINC01526 in gastric cancer. Furthermore, RNA pull-down and RNA immunoprecipitation chip assays demonstrated that LINC01526 interacts with TARBP2, an RNA-binding protein controlling mRNA stability. Moreover, TARBP2 could bind and destabilize GNG7 transcripts. Finally, the rescue assay disclosed that LINC01526 promoted gastric cancer progression by interacting with TARBP2, leading to the degradation of GNG7 mRNA.

Abstract

Gastric cancer is the most common malignancy of the human digestive system. Long noncoding RNAs (lncRNAs) influence the occurrence and development of gastric cancer in multiple ways. However, the function and mechanism of LINC01526 in gastric cancer remain unknown. Herein, we investigated the function of LINC01526 with respect to the malignant progression of gastric cancer. We found that LINC01526 was upregulated in gastric cancer cells and tissues. The function experiments in vitro and the Xenograft mouse model in vivo proved that LINC01526 could promote gastric cancer cell proliferation and migration. Furthermore, LINC01526 interacted with TAR (HIV-1) RNA-binding protein 2 (TARBP2) and decreased the mRNA stability of G protein gamma 7 (GNG7) through TARBP2. Finally, the rescue assay showed that downregulating GNG7 partially rescued the cell proliferation inhibited by LINC01526 or TARBP2 silencing. In summary, LINC01526 promoted gastric cancer progression by interacting with TARBP2, which subsequently degraded GNG7 mRNA. This study not only explores the role of LINC01526 in gastric cancer, but also provides a laboratory basis for its use as a new biomarker for diagnosis and therapeutic targets.

Emerging biology of noncoding RNAs in malaria parasites

In eukaryotic organisms, noncoding RNAs (ncRNAs) have been implicated as important regulators of multifaceted biological processes, including transcriptional, posttranscriptional, and epigenetic regulation of gene expression. In recent years, it is becoming clear that protozoan parasites encode diverse ncRNA transcripts; however, little is known about their cellular functions. Recent advances in high-throughput “omic” studies identified many novel long ncRNAs (lncRNAs) in apicomplexan parasites, some of which undergo splicing, polyadenylation, and encode small proteins. To date, only a few of them are characterized, leaving a big gap in our understanding regarding their origin, mode of action, and functions in parasite biology. In this review, we focus on lncRNAs of the human malaria parasite Plasmodium falciparum and highlight their cellular functions and possible mechanisms of action.

Non-coding RNAs in malaria infection

Malaria is one of the most severe infectious diseases affecting humans and it is caused by protozoan pathogens of the species Plasmodium (spp.). The malaria parasite Plasmodium is characterized by a complex, multistage life cycle that requires tight gene regulation which allows for host invasion and defense against host immune responses. Unfortunately, the mechanisms regulating gene expression during Plasmodium infection remain largely elusive, though several lines of evidence implicate a major involvement of non-coding RNAs (ncRNAs). The ncRNAs have been found to play a key role in regulating transcriptional and post-transcriptional events in a broad range of organisms including Plasmodium. In Plasmodium ncRNAs have been shown to regulate key events in the multistage life cycle and virulence ability. Here we review recent progress involving ncRNAs (microRNAs, long non-coding RNAs, and circular RNAs) and their role as regulators of gene expression during Plasmodium infection in human hosts with focus on the possibility of using these molecules as biomarkers for monitoring disease status. We also discuss the surprising function of ncRNAs in mediating the complex interplay between parasite and human host and future perspectives of the field. This article is categorized under: RNA in Disease and Development > RNA in Disease.

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.

A Bioinformatics Approach to Identifying Potential Biomarkers for Cryptosporidium parvum: A Coccidian Parasite Associated with Fetal Diarrhea

Cryptosporidium parvum (C. parvum) is a protozoan parasite known for cryptosporidiosis in pre-weaned calves. Animals and patients with immunosuppression are at risk of developing the disease, which can cause potentially fatal diarrhoea. The present study aimed to construct a network biology framework based on the differentially expressed genes (DEGs) of C. parvum infected subjects. In this way, the gene expression profiling analysis of C. parvum infected individuals can give us a snapshot of actively expressed genes and transcripts under infection conditions. In the present study, we have analyzed microarray data sets and compared the gene expression profiles of the patients with the different data sets of the healthy control. Using a network medicine approach to identify the most influential genes in the gene interaction network, we uncovered essential genes and pathways related to C. parvum infection. We identified 164 differentially expressed genes (109 up- and 54 down-regulated DEGs) and allocated them to pathway and gene set enrichment analysis. The results underpin the identification of seven significant hub genes with high centrality values: ISG15, MX1, IFI44L, STAT1, IFIT1, OAS1, IFIT3, RSAD2, IFITM1, and IFI44. These genes are associated with diverse biological processes not limited to host interaction, type 1 interferon production, or response to IL-gamma. Furthermore, four genes (IFI44, IFIT3, IFITM1, and MX1) were also discovered to be involved in innate immunity, inflammation, apoptosis, phosphorylation, cell proliferation, and cell signaling. In conclusion, these results reinforce the development and implementation of tools based on gene profiles to identify and treat Cryptosporidium parvum-related diseases at an early stage.

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.

Full-Length Transcriptome Analysis of Plasmodium falciparum by Single-Molecule Long-Read Sequencing

Malaria, an infectious disease caused by Plasmodium parasites, still accounts for amounts of deaths annually in last decades. Despite the significance of Plasmodium falciparum as a model organism of malaria parasites, our understanding of gene expression of this parasite remains largely elusive since lots of progress on its genome and transcriptome are based on assembly with short sequencing reads. Herein, we report the new version of transcriptome dataset containing all full-length transcripts over the whole asexual blood stages by adopting a full-length sequencing approach with optimized experimental conditions of cDNA library preparation. We have identified a total of 393 alternative splicing (AS) events, 3,623 long non-coding RNAs (lncRNAs), 1,555 alternative polyadenylation (APA) events, 57 transcription factors (TF), 1,721 fusion transcripts in P. falciparum. Furthermore, the shotgun proteome was performed to validate the full-length transcriptome of P. falciparum. More importantly, integration of full-length transcriptomic and proteomic data identified 160 novel small proteins in lncRNA regions. Collectively, this full-length transcriptome dataset with high quality and accuracy and the shotgun proteome analyses shed light on the complex gene expression in malaria parasites and provide a valuable resource for related functional and mechanistic researches on P. falciparum genes.

Analysis of Long Non-Coding RNA in Cryptosporidium parvum Reveals Significant Stage-Specific Antisense Transcription

Cryptosporidium is a protist parasite that has been identified as the second leading cause of moderate to severe diarrhea in children younger than two and a significant cause of mortality worldwide. Cryptosporidium has a complex, obligate, intracellular but extra cytoplasmic lifecycle in a single host. How genes are regulated in this parasite remains largely unknown. Long non-coding RNAs (lncRNAs) play critical regulatory roles, including gene expression across a broad range of organisms. Cryptosporidium lncRNAs have been reported to enter the host cell nucleus and affect the host response. However, no systematic study of lncRNAs in Cryptosporidium has been conducted to identify additional lncRNAs. In this study, we analyzed a C. parvum in vitro strand-specific RNA-seq developmental time series covering both asexual and sexual stages to identify lncRNAs associated with parasite development. In total, we identified 396 novel lncRNAs, mostly antisense, with 86% being differentially expressed. Surprisingly, nearly 10% of annotated mRNAs have an antisense transcript. lncRNAs occur most often at the 3' end of their corresponding sense mRNA. Putative lncRNA regulatory regions were identified and many appear to encode bidirectional promoters. A positive correlation between lncRNA and upstream mRNA expression was observed. Evolutionary conservation and expression of lncRNA candidates was observed between C. parvum, C. hominis and C. baileyi. Ten C. parvum protein-encoding genes with antisense transcripts have P. falciparum orthologs that also have antisense transcripts. Three C. parvum lncRNAs with exceptional properties (e.g., intron splicing) were experimentally validated using RT-PCR and RT-qPCR. This initial characterization of the C. parvum non-coding transcriptome facilitates further investigations into the roles of lncRNAs in parasite development and host-pathogen interactions.

Long non-coding RNAs towards precision medicine in gastric cancer: early diagnosis, treatment, and drug resistance

Gastric cancer is a deadly disease and remains the third leading cause of cancer-related death worldwide. The 5-year overall survival rate of patients with early-stage localized gastric cancer is more than 60%, whereas that of patients with distant metastasis is less than 5%. Surgical resection is the best option for early-stage gastric cancer, while chemotherapy is mainly used in the middle and advanced stages of this disease, despite the frequently reported treatment failure due to chemotherapy resistance. Therefore, there is an unmet medical need for identifying new biomarkers for the early diagnosis and proper management of patients, to achieve the best response to treatment. Long non-coding RNAs (lncRNAs) in body fluids have attracted widespread attention as biomarkers for early screening, diagnosis, treatment, prognosis, and responses to drugs due to the high specificity and sensitivity. In the present review, we focus on the clinical potential of lncRNAs as biomarkers in liquid biopsies in the diagnosis and prognosis of gastric cancer. We also comprehensively discuss the roles of lncRNAs and their molecular mechanisms in gastric cancer chemoresistance as well as their potential as therapeutic targets for gastric cancer precision medicine.

The cryptic unstable transcripts are associated with developmentally regulated gene expression in blood-stage Plasmodium falciparum

The tight gene expression regulation controls the development and pathogenesis of human malaria parasite Plasmodium falciparum throughout the complex life cycle. Recent studies have revealed the pervasive nascent transcripts in the genome of P. falciparum, suggesting the existence of a hidden transcriptome involved in the dynamic gene expression. However, the landscape and related biological functions of nascent non-coding RNAs (ns-ncRNAs) are still poorly explored. Here we profiled the transcription dynamics of nascent RNAs by rRNA-depleted and stranded RNA sequencing over the course of 48-h intraerythrocytic developmental cycle (IDC). We identified the genome-wide sources of a total of 2252 ns-ncRNAs, mostly originating from intergenic and untranslated regions of annotated genes. By integrating the nascent RNA abundances with ATAC-seq and ChIP-seq analysis, we uncovered the euchromatic microenvironment surrounding the ns-ncRNA loci, and revealed a positive correlation between ns-ncRNAs and corresponding mRNA abundances. Finally, by gene knock-down strategy, we showed that the cooperation of RNA exosome catalytic subunit PfDis3 and PfMtr4 cofactor played a major role in ns-ncRNAs degradation. Collectively, this study contributes to understanding of the potential roles of short-lived nascent ncRNAs in regulating gene expression in malaria parasites.

Plasma mEV levels in Ghanain malaria patients with low parasitaemia are higher than those of healthy controls, raising the potential for parasite markers in mEVs as diagnostic targets

This study sought to measure medium-sized extracellular vesicles (mEVs) in plasma, when patients have low Plasmodium falciparum early in infection. We aimed to define the relationship between plasma mEVs and: (i) parasitaemia, (ii) period from onset of malaria symptoms until seeking medical care (patient delay, PD), (iii) age and (iv) gender. In this cross-sectional study, n = 434 patients were analysed and Nanosight Tracking Analysis (NTA) used to quantify mEVs (vesicles of 150–500 nm diameter, isolated at 15,000 × g, β-tubulin-positive and staining for annexin V, but weak or negative for CD81). Overall plasma mEV levels (1.69 × 10¹⁰ mEVs mL⁻¹) were 2.3-fold higher than for uninfected controls (0.51 × 10¹⁰ mEVs mL⁻¹). Divided into four age groups, we found a bimodal distribution with 2.5- and 2.1-fold higher mEVs in infected children (<11 years old [yo]) (median:2.11 × 10¹⁰ mEVs mL⁻¹) and the elderly (>45 yo) (median:1.92 × 10¹⁰ mEVs mL⁻¹), respectively, compared to uninfected controls; parasite density varied similarly with age groups. There was a positive association between mEVs and parasite density (r = 0.587, p < 0.0001) and mEVs were strongly associated with PD (r = 0.919, p < 0.0001), but gender had no effect on plasma mEV levels (p = 0.667). Parasite density was also exponentially related to patient delay. Gender (p = 0.667) had no effect on plasma mEV levels. During periods of low parasitaemia (PD = 72h), mEVs were 0.93-fold greater than in uninfected controls. As 75% (49/65) of patients had low parasitaemia levels (20–500 parasites µL⁻¹), close to the detection limits of microscopy of Giemsa-stained thick blood films (5–150 parasites µL⁻¹), mEV quantification by NTA could potentially have early diagnostic value, and raises the potential of Pf markers in mEVs as early diagnostic targets.

Involvement of Cryptosporidium parvum Cdg7_FLc_1000 RNA in the Attenuation of Intestinal Epithelial Cell Migration via Trans-Suppression of Host Cell SMPD3

Intestinal infection by Cryptosporidium parvum causes inhibition of epithelial turnover, but underlying mechanisms are unclear. Previous studies demonstrate that a panel of parasite RNA transcripts of low protein-coding potential are delivered into infected epithelial cells. Using in vitro and in vivo models of intestinal cryptosporidiosis, we report here that host delivery of parasite Cdg7_FLc_1000 RNA results in inhibition of epithelial cell migration through suppression of the gene encoding sphingomyelinase 3 (SMPD3). Delivery of Cdg7_FLc_1000 into infected cells promotes the histone methyltransferase G9a-mediated H3K9 methylation in the SMPD3 locus. The DNA-binding transcriptional repressor, PR domain zinc finger protein 1, is required for the assembly of Cdg7_FLc_1000 into the G9a complex and associated with the enrichment of H3K9 methylation at the gene locus. Pathologically, nuclear transfer of Cryptosporidium parvum Cdg7_FLc_1000 RNA is involved in the attenuation of intestinal epithelial cell migration via trans-suppression of host cell SMPD3.

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.

Identification and comparison of long non-coding RNAs in the silk gland between domestic and wild silkworms

Under long-term artificial selection, the domestic silkworm (Bombyx mori) has increased its silk yield tremendously in comparison with its wild progenitor, Bombyx mandarina. However, the molecular mechanism of silk yield increase is still unknown. Comparative analysis of long non-coding RNAs (lncRNAs) may provide some insights into understanding this phenotypic variation. In this study, using RNA sequencing technology data of silk gland in domestic and wild silkworms, we identified 599 lncRNAs in the silk gland of the silkworm. Compared with protein-coding genes, the silk gland lncRNA genes tend to have fewer exon numbers, shorter transcript length and lower GC-content. Moreover, we found that three lncRNA genes are significantly and differentially expressed between domestic and wild silkworms. The potential targets of two differentially expressed lncRNAs (DELs) (dw4sg_0040 and dw4sg_0483) and the expression-correlated genes with the two DELs are mainly enriched in the related processes of silk protein translation. This implies that these DELs may affect the phenotypic variation in silk yield between the domestic and wild silkworms through the post-transcriptional regulation of silk protein.

Identification of long noncoding RNAs in Schistosoma mansoni and Schistosoma japonicum

Schistosomiasis is a major parasitic disease caused by 3 principal species of schistosome. Studies of schistosome transcriptomes have focused on protein-coding transcripts and although miRNAs are attracting increased attention, few reports have concerned the long noncoding RNAs (lncRNAs). These have been shown to play key roles in the regulation of gene expression through interactions with mRNAs, proteins and miRNAs. In this study, we first identified lncRNAs from RNA-seq data in Schistosoma mansoni and Schistosoma japonicum: 3247 and 3033 potential lncRNAs were found in these two species respectively. ChIP-seq analysis to determine H3K4me3 profiles along the gene regions corresponding to lncRNAs showed that in 12% of cases this mark was enriched in regions proximal to the transcription start sites, supporting their validity as actively transcribed genes. Besides, the sequence conservation of lncRNAs between schistosome species was much lower than that of mRNAs, but higher than that of the randomly selected genomic sequences, which is consistent with that in mammals. Our results demonstrate that lncRNAs form a significant part of the schistosome transcriptome and suggest that they play an important role in the biology of the parasite.

Identification of human Phosphatidyl Inositol 5-Phosphate 4-Kinase as an RNA binding protein that is imported into Plasmodium falciparum

Plasmodium falciparum is a causative agent for malaria and has a complex life cycle in human and mosquito hosts. Translation repression of specific set of mRNA has been reported in gametocyte stages of this parasite. A conserved element present in the 3'UTR of some of these transcripts was identified. Biochemical studies have identified components of the RNA storage and/or translation inhibitor complex but it is not yet clear how the complex is specifically recruited on the RNA targeted for translation regulation. We used the 3'UTR region of translationally regulated transcripts to identify Phosphatidyl-inositol 5-phosphate 4-kinase (PIP4K2A) as the protein that associates with these RNAs. We further show that recombinant PIP4K2A has the RNA binding activity and can associate specifically with Plasmodium 3'UTR RNAs. Immunostainings show that hPIP4K2A is imported into the Plasmodium parasite from RBC. These results identify a novel RNA binding role for PIP4K2A that may play a role in Plasmodium propagation.

Trans-suppression of host CDH3 and LOXL4 genes during Cryptosporidium parvum infection involves nuclear delivery of parasite Cdg7_FLc_1000 RNA

Intestinal infection by Cryptosporidium parvum causes significant alterations in the gene expression profile in host epithelial cells. Previous studies demonstrate that a panel of parasite RNA transcripts of low protein-coding potential are delivered into infected host cells and may modulate host gene transcription. Using in vitro models of human intestinal cryptosporidiosis, we report here that trans-suppression of the cadherin 3 (CDH3) and lysyl oxidase like 4 (LOXL4) genes in human intestinal epithelial cells following C. parvum infection involves host delivery of the Cdg7_FLc_1000 RNA, a C. parvum RNA that has been previously demonstrated to be delivered into the nuclei of infected host cells. Downregulation of CDH3 and LOXL4 genes was detected in host epithelial cells following C. parvum infection or in cells expressing the parasite Cdg7_FLc_1000 RNA. Knockdown of Cdg7_FLc_1000 attenuated the trans-suppression of CDH3 and LOXL4 genes in host cells induced by infection. Interestingly, Cdg7_FLc_1000 was detected to be recruited to the promoter regions of both CDH3 and LOXL4 gene loci in host cells following C. parvum infection. Host delivery of Cdg7_FLc_1000 promoted the PH domain zinc finger protein 1 (PRDM1)-mediated H3K9 methylation associated with trans-suppression in the CDH3 gene locus, but not the LOXL4 gene. Therefore, our data suggest that host delivery of Cdg7_FLc_1000 causes CDH3 trans-suppression in human intestinal epithelial cells following C. parvum infection through PRDM1-mediated H3K9 methylation in the CDH3 gene locus, whereas Cdg7_FLc_1000 induces trans-suppression of the host LOXL4 gene through H3K9/H3K27 methylation-independent mechanisms.

Trans-suppression of defense DEFB1 gene in intestinal epithelial cells following Cryptosporidium parvum infection is associated with host delivery of parasite Cdg7_FLc_1000 RNA

To counteract host immunity, Cryptosporidium parvum has evolved multiple strategies to suppress host antimicrobial defense. One such strategy is to reduce the production of the antimicrobial peptide beta-defensin 1 (DEFB1) by host epithelial cells but the underlying mechanisms remain unclear. Recent studies demonstrate that a panel of parasite RNA transcripts of low protein-coding potential are delivered into infected host cells and may modulate host gene transcription. Using in vitro models of intestinal cryptosporidiosis, in this study, we analyzed the expression profile of host beta-defensin genes in host cells following infection. We found that C. parvum infection caused a significant downregulation of the DEFB1 gene. Interestingly, downregulation of DEFB1 gene was associated with host delivery of Cdg7_FLc_1000 RNA transcript, a C. parvum RNA that has previously demonstrated to be delivered into the nuclei of infected host cells. Knockdown of Cdg7_FLc_1000 in host cells could attenuate the trans-suppression of host DEFB1 gene and decreased the parasite burden. Therefore, our data suggest that trans-suppression of DEFB1 gene in intestinal epithelial cells following C. parvum infection involves host delivery of parasite Cdg7_FLc_1000 RNA, a process that may be relevant to the epithelial defense evasion by C. parvum at the early stage of infection.

Nuclear delivery of parasite Cdg2_FLc_0220 RNA transcript to epithelial cells during Cryptosporidium parvum infection modulates host gene transcription

Intestinal infection by the zoonotic protozoan, Cryptosporidium parvum, causes significant alterations in the gene expression profile in host epithelial cells. The molecular mechanisms of how C. parvum may modulate host cell gene transcription and the pathological significance of such alterations are largely unclear. Previous studies demonstrate that a panel of parasite RNA transcripts are delivered into infected host cells and may modulate host gene transcription. Using in vitro models of intestinal cryptosporidiosis, in this study, we analyzed the impact of host delivery of C. parvum Cdg2_FLc_0220 RNA transcript on host gene expression profile. We found that alterations in host gene expression profile following C. parvum infection were partially associated with the nuclear delivery of Cdg2_FLc_0220. Specifically, we identified a total of 46 overlapping upregulated genes and 8 overlapping downregulated genes in infected cells and cells transfected with Full-Cdg2_FLc_0220. Trans-suppression of the DAZ interacting zinc finger protein 1 like (DZIP1L) gene, the top overlapping downregulated gene in host cells following C. parvum infection and cells transfected with Full-Cdg2_FLc_0220, was mediated by G9a, independent of PRDM1. Cdg2_FLc_0220-mediated trans-suppression of the DZIP1L gene was independent of H3K9 and H3K27 methylation. Data from this study provide additional evidence that delivery of C. parvum Cdg2_FLc_0220 RNA transcript in infected epithelial cells modulates the transcription of host genes, contributing to the alterations in the gene expression profile in host epithelial cells during C. parvum infection.

Characterization of P. vivax blood stage transcriptomes from field isolates reveals similarities among infections and complex gene isoforms

Our understanding of the structure and regulation of Plasmodium vivax genes is limited by our inability to grow the parasites in long-term in vitro cultures. Most P. vivax studies must therefore rely on patient samples, which typically display a low proportion of parasites and asynchronous parasites. Here, we present stranded RNA-seq data generated directly from a small volume of blood from three Cambodian vivax malaria patients collected before treatment. Our analyses show surprising similarities of the parasite gene expression patterns across infections, despite extensive variations in parasite stage proportion. These similarities contrast with the unique gene expression patterns observed in sporozoites isolated from salivary glands of infected Colombian mosquitoes. Our analyses also indicate that more than 10% of P. vivax genes encode multiple, often undescribed, protein-coding sequences, potentially increasing the diversity of proteins synthesized by blood stage parasites. These data also greatly improve the annotations of P. vivax gene untranslated regions, providing an important resource for future studies of specific genes.

Delivery of parasite Cdg7_Flc_0990 RNA transcript into intestinal epithelial cells during Cryptosporidium parvum infection suppresses host cell gene transcription through epigenetic mechanisms

Cryptosporidial infection causes dysregulated transcription of host genes key to intestinal epithelial homeostasis, but the underlying mechanisms remain obscure. Previous studies demonstrate that several Cryptosporidium parvum (C. parvum) RNA transcripts are selectively delivered into epithelial cells during host cell invasion and may modulate gene transcription in infected cells. We report here that C. parvum infection suppresses the transcription of LRP5, SLC7A8, and IL33 genes in infected intestinal epithelium. Trans-suppression of these genes in infected host cells is associated with promoter enrichment of suppressive epigenetic markers (i.e., H3K9me3). Cdg7_FLc_0990, a C. parvum RNA that has previously demonstrated to be delivered into the nuclei of infected epithelial cells, is recruited to the promoter regions of LRP5, SLC7A8, and IL33 genes. Cdg7_FLc_0990 appears to be recruited to their promoter regions together with G9a, a histone methyltransferase for H3K9 methylation. The PR domain zinc finger protein 1, a G9a-interacting protein, is required for the assembly of Cdg7_FLc_0990 to the G9a complex and gene-specific enrichment of H3K9 methylation. Our data demonstrate that cryptosporidial infection induces epigenetic histone methylations in infected cells through nuclear transfer of parasite Cdg7_Flc_0990 RNA transcript, resulting in transcriptional suppression of the LRP5, SLC7A8, and IL33 genes.

Delivery of Parasite RNA Transcripts Into Infected Epithelial Cells During Cryptosporidium Infection and Its Potential Impact on Host Gene Transcription

Cryptosporidium parvum is an important opportunistic parasite pathogen for immunocompromised individuals and a common cause of diarrhea in young children. Previous studies have identified a panel of RNA transcripts of very low protein-coding potential in C. parvum. Using an in vitro model of human intestinal cryptosporidiosis, we report here that some of these C. parvum RNA transcripts were selectively delivered into the nuclei of host epithelial cells during C. parvum infection. Nuclear delivery of several such parasitic RNAs, including Cdg7_FLc_0990, involved heat-shock protein 70-mediated nuclear importing mechanism. Overexpression of Cdg7_FLc_0990 in intestinal epithelial cells resulted in significant changes in expression levels of specific genes, with significant overlapping with alterations in gene expression profile detected in host cells after C. parvum infection. Our data demonstrate that C. parvum transcripts of low protein-coding potential are selectively delivered into epithelial cells during infection and may modulate gene transcription in infected host cells.

Epigenetic landscapes underlining global patterns of gene expression in the human malaria parasite, Plasmodium falciparum

The dynamic chromatin landscape displaying combinatorial complexity of the epigenome impacts gene expression that underlies many events of differentiation and cell cycle progression. In the past few years, epigenetic mechanisms have emerged as important processes involved in the tight gene regulation in malaria parasites, Plasmodium spp. Focusing predominantly on Plasmodium falciparum, the species associated with the most severe form of the disease, many advances have been made in our understanding of the interaction between transcriptional regulation and epigenetic mechanisms as the pivotal processes in regulating life cycle progression, host parasite interactions and parasite adaptation to the host environment. This review focuses on the epigenome and its effect on transcriptional regulation in P. falciparum, highlighting its unique, evolutionary diverse features.

BmncRNAdb: a comprehensive database of non-coding RNAs in the silkworm, Bombyx mori

BACKGROUND

Long non-coding RNAs (lncRNAs) may play critical roles in a wide range of developmental processes of higher organisms. Recently, lncRNAs have been widely identified across eukaryotes and many databases of lncRNAs have been developed for human, mouse, fruit fly, etc. However, there is rare information about them in the only completely domesticated insect, silkworm (Bombyx mori).

DESCRIPTION

In this study, we systematically scanned lncRNAs using the available silkworm RNA-seq data and public unigenes. Finally, we identified and collected 6281 lncRNAs in the silkworm. Besides, we also collected 1986 microRNAs (miRNAs) from previous studies. Then, we organized them into a comprehensive and web-based database, BmncRNAdb. This database offers a user-friendly interface for data browse and online analysis as well as the three online tools for users to predict the target genes of lncRNA or miRNA.

CONCLUSIONS

We have systematically identified and collected the silkworm lncRNAs and constructed a comprehensive database of the silkworm lncRNAs and miRNAs. This work gives a glimpse into lncRNAs of the silkworm and lays foundations for the ncRNAs study of the silkworm and other insects in the future. The BmncRNAdb is freely available at http://gene.cqu.edu.cn/BmncRNAdb/index.php .

HP1BP3, a Chromatin Retention Factor for Co-transcriptional MicroRNA Processing

Recent studies suggest that the microprocessor (Drosha-DGCR8) complex can be recruited to chromatin to catalyze co-transcriptional processing of primary microRNAs (pri-miRNAs) in mammalian cells. However, the molecular mechanism of co-transcriptional miRNA processing is poorly understood. Here we find that HP1BP3, a histone H1-like chromatin protein, specifically associates with the microprocessor and promotes global miRNA biogenesis in human cells. Chromatin immunoprecipitation (ChIP) studies reveal genome-wide co-localization of HP1BP3 and Drosha and HP1BP3-dependent Drosha binding to actively transcribed miRNA loci. Moreover, HP1BP3 specifically binds endogenous pri-miRNAs and facilitates the Drosha/pri-miRNA association in vivo. Knockdown of HP1BP3 compromises pri-miRNA processing by causing premature release of pri-miRNAs from the chromatin. Taken together, these studies suggest that HP1BP3 promotes co-transcriptional miRNA processing via chromatin retention of nascent pri-miRNA transcripts. This work significantly expands the functional repertoire of the H1 family of proteins and suggests the existence of chromatin retention factors for widespread co-transcriptional miRNA processing.

Expression of the Plasmodium berghei actin II gene is controlled by elements in a long genomic region

Plasmodium parasites have two actin isoforms. Actin I is ubiquitously expressed, while the second actin isoform is expressed in the sexual stages and ookinetes. Reverse genetic analysis revealed two phenotypes in parasites lacking the protein: a block in male gametogenesis (exflagellation) and a second phenotype in oocyst development, dependent upon the expression of the gene in female gametocytes. Here, we report that the genetic complementation of two independent mutants lacking actin II does not fully restore wild-type function. Constructs were integrated in the c-rrna locus, previously used for expression of transgenes, in order to determine the dependence of expression on actin II flanking genomic regions. Partial restoration of male gametogenesis was achieved when the transgene contained, in addition to the coding region, 1.2 kb upstream of the actin II open reading frame. Another transgene, which comprised 2.7 kb of actin II 5' flanking regions and the cognate 3' downstream sequence, fully restored exflagellation. However, in both complemented strains, oocyst development was severely impaired compared to the WT. These data suggest that male gametocyte expression of actin II is dependent upon extensive flanking regions, while female expression requires even longer genomic sequences for correct expression of the gene.

Genome-Wide Identification and Characterization of Long Non-Coding RNAs from Mulberry (Morus notabilis) RNA-seq Data

Numerous sources of evidence suggest that most of the eukaryotic genome is transcribed into protein-coding mRNAs and also into a large number of non-coding RNAs (ncRNAs). Long ncRNAs (lncRNAs), a group consisting of ncRNAs longer than 200 nucleotides, have been found to play critical roles in transcriptional, post-transcriptional, and epigenetic gene regulation across all kingdoms of life. However, lncRNAs and their regulatory roles remain poorly characterized in plants, especially in woody plants. In this paper, we used a computational approach to identify novel lncRNAs from a published RNA-seq data set and analyzed their sequences and expression patterns. In total, 1133 novel lncRNAs were identified in mulberry, and 106 of these lncRNAs displayed a predominant tissue-specific expression in the five major tissues investigated. Additionally, functional predictions revealed that tissue-specific lncRNAs adjacent to protein-coding genes might play important regulatory roles in the development of floral organ and root in mulberry. The pipeline used in this study would be useful for the identification of lncRNAs obtained from other deep sequencing data. Furthermore, the predicted lncRNAs would be beneficial towards an understanding of the variations in gene expression in plants.

Analysis of nucleosome positioning landscapes enables gene discovery in the human malaria parasite Plasmodium falciparum

Background

Plasmodium falciparum, the deadliest malaria-causing parasite, has an extremely AT-rich (80.7 %) genome. Because of high AT-content, sequence-based annotation of genes and functional elements remains challenging. In order to better understand the regulatory network controlling gene expression in the parasite, a more complete genome annotation as well as analysis tools adapted for AT-rich genomes are needed. Recent studies on genome-wide nucleosome positioning in eukaryotes have shown that nucleosome landscapes exhibit regular characteristic patterns at the 5’- and 3’-end of protein and non-protein coding genes. In addition, nucleosome depleted regions can be found near transcription start sites. These unique nucleosome landscape patterns may be exploited for the identification of novel genes. In this paper, we propose a computational approach to discover novel putative genes based exclusively on nucleosome positioning data in the AT-rich genome of P. falciparum.

Results

Using binary classifiers trained on nucleosome landscapes at the gene boundaries from two independent nucleosome positioning data sets, we were able to detect a total of 231 regions containing putative genes in the genome of Plasmodium falciparum, of which 67 highly confident genes were found in both data sets. Eighty-eight of these 231 newly predicted genes exhibited transcription signal in RNA-Seq data, indicative of active transcription. In addition, 20 out of 21 selected gene candidates were further validated by RT-PCR, and 28 out of the 231 genes showed significant matches using BLASTN against an expressed sequence tag (EST) database. Furthermore, 108 (47 %) out of the 231 putative novel genes overlapped with previously identified but unannotated long non-coding RNAs. Collectively, these results provide experimental validation for 163 predicted genes (70.6 %). Finally, 73 out of 231 genes were found to be potentially translated based on their signal in polysome-associated RNA-Seq representing transcripts that are actively being translated.

Conclusion

Our results clearly indicate that nucleosome positioning data contains sufficient information for novel gene discovery. As distinct nucleosome landscapes around genes are found in many other eukaryotic organisms, this methodology could be used to characterize the transcriptome of any organism, especially when coupled with other DNA-based gene finding and experimental methods (e.g., RNA-Seq).

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-2214-9) contains supplementary material, which is available to authorized users.

Genome-wide identification and functional analysis of lincRNAs acting as miRNA targets or decoys in maize

Background

Long intergenic noncoding RNAs (lincRNAs) are endogenous non-coding RNAs (ncRNAs) that are transcribed from ‘intergenic’ regions of the genome and may play critical roles in regulating gene expression through multiple RNA-mediated mechanisms. MicroRNAs (miRNAs) are single-stranded small ncRNAs of approximately 21–24 nucleotide (nt) that are involved in transcriptional and post-transcriptional gene regulation. While miRNAs functioning as mRNA repressors have been studied in detail, the influence of miRNAs on lincRNAs has seldom been investigated in plants.

Methods

LincRNAs as miRNA targets or decoys were predicted via GSTAr.pl script with a set of rules, and lincRNAs as miRNA targets were validated by degradome data. Conservation analysis of lincRNAs as miRNA targets or decoys were conducted using BLASTN and MAFFT. The function of lincRNAs as miRNA targets were predicted via a lincRNA-mRNA co-expression network, and the function of lincRNAs as miRNA decoys were predicted according to the competing endogenous RNA (ceRNA) hypothesis.

Results

In this work, we developed a computational method and systematically predicted 466 lincRNAs as 165 miRNA targets and 86 lincRNAs as 58 miRNA decoys in maize (Zea mays L.). Furthermore, 34 lincRNAs predicted as 33 miRNA targets were validated based on degradome data. We found that lincRNAs acting as miRNA targets or decoys are a common phenomenon, which indicates that the regulated networks of miRNAs also involve lincRNAs. To elucidate the function of lincRNAs, we reconstructed a miRNA-regulated network involving 78 miRNAs, 117 lincRNAs and 8834 mRNAs. Based on the lincRNA-mRNA co-expression network and the competing endogenous RNA hypothesis, we predicted that 34 lincRNAs that function as miRNA targets and 86 lincRNAs that function as miRNA decoys participate in cellular and metabolic processes, and play role in catalytic activity and molecular binding functions.

Conclusions

This work provides a comprehensive view of miRNA-regulated networks and indicates that lincRNAs can participate in a layer of regulatory interactions as miRNA targets or decoys in plants, which will enable in-depth functional analysis of lincRNAs.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-2024-0) contains supplementary material, which is available to authorized users.

Genome-wide characterization of long intergenic non-coding RNAs (lincRNAs) provides new insight into viral diseases in honey bees Apis cerana and Apis mellifera

BACKGROUND

Long non-coding RNAs (lncRNAs) are a class of RNAs that do not encode proteins. Recently, lncRNAs have gained special attention for their roles in various biological process and diseases.

RESULTS

In an attempt to identify long intergenic non-coding RNAs (lincRNAs) and their possible involvement in honey bee development and diseases, we analyzed RNA-seq datasets generated from Asian honey bee (Apis cerana) and western honey bee (Apis mellifera). We identified 2470 lincRNAs with an average length of 1011 bp from A. cerana and 1514 lincRNAs with an average length of 790 bp in A. mellifera. Comparative analysis revealed that 5 % of the total lincRNAs derived from both species are unique in each species. Our comparative digital gene expression analysis revealed a high degree of tissue-specific expression among the seven major tissues of honey bee, different from mRNA expression patterns. A total of 863 (57 %) and 464 (18 %) lincRNAs showed tissue-dependent expression in A. mellifera and A. cerana, respectively, most preferentially in ovary and fat body tissues. Importantly, we identified 11 lincRNAs that are specifically regulated upon viral infection in honey bees, and 10 of them appear to play roles during infection with various viruses.

CONCLUSIONS

This study provides the first comprehensive set of lincRNAs for honey bees and opens the door to discover lincRNAs associated with biological and hormone signaling pathways as well as various diseases of honey bee.

Dynamic and Widespread lncRNA Expression in a Sponge and the Origin of Animal Complexity

Long noncoding RNAs (lncRNAs) are important developmental regulators in bilaterian animals. A correlation has been claimed between the lncRNA repertoire expansion and morphological complexity in vertebrate evolution. However, this claim has not been tested by examining morphologically simple animals. Here, we undertake a systematic investigation of lncRNAs in the demosponge Amphimedon queenslandica, a morphologically simple, early-branching metazoan. We combine RNA-Seq data across multiple developmental stages of Amphimedon with a filtering pipeline to conservatively predict 2,935 lncRNAs. These include intronic overlapping lncRNAs, exonic antisense overlapping lncRNAs, long intergenic nonprotein coding RNAs, and precursors for small RNAs. Sponge lncRNAs are remarkably similar to their bilaterian counterparts in being relatively short with few exons and having low primary sequence conservation relative to protein-coding genes. As in bilaterians, a majority of sponge lncRNAs exhibit typical hallmarks of regulatory molecules, including high temporal specificity and dynamic developmental expression. Specific lncRNA expression profiles correlate tightly with conserved protein-coding genes likely involved in a range of developmental and physiological processes, such as the Wnt signaling pathway. Although the majority of Amphimedon lncRNAs appears to be taxonomically restricted with no identifiable orthologs, we find a few cases of conservation between demosponges in lncRNAs that are antisense to coding sequences. Based on the high similarity in the structure, organization, and dynamic expression of sponge lncRNAs to their bilaterian counterparts, we propose that these noncoding RNAs are an ancient feature of the metazoan genome. These results are consistent with lncRNAs regulating the development of animals, regardless of their level of morphological complexity.