HumanViCe: host ceRNA network in virus infected cells in human

Construction and Analysis of ceRNA Networks Reveal the Key Genes Associated with Bovine Herpesvirus Type 1 Infection

Background

Virus infection can cause the changes of lncRNA expression levels to regulate the interaction between virus and host, but the relationship between BHV-1 infection and lncRNA has not been reported.

Methods

In this study, in order to reveal the molecular mechanism of RNA in BoHV-1 infection, the Madin-Darby bovine kidney (MDBK) cells were infected with BoHV-1, transcriptome sequencing were performed by next-generation sequencing at 18 h or 24 h or 33 h of viral infection and then based on the competitive endogenous RNA (ceRNA) theory, lncRNA-miRNA-mRNA networks were constructed using these high-throughput sequencing data. The network analysis, Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis were performed for functional annotation and exploration of ncRNA ceRNAs in BoHV-1 infection.

Results

The results showed that 48 lncRNAs, 123 mRNAs and 20 miRNAs as differentially expressed genes, and the mitogen activated protein kinase (MAPK) pathway and calcium signaling pathway were significantly enriched in the ceRNA network. Some differentially expressed lncRNA genes were randomly selected for verification by RT-qPCR, and the results showed that their expression trend was consistent with the results of transcriptome sequencing data.

Conclusion

This study revealed that BoHV-1 infection can affect the expression of RNAs in MDBK cells and the regulation of ceRNA network to carry out corresponding biological functions in the host, but further experimental studies are still necessary to prove the hub genes function in ceRNA network and the molecular mechanism in BoHV-1 infection.

Circular RNAs in hepatocellular carcinoma: biogenesis, function, and pathology

Hepatocellular carcinoma (HCC) is one of the most common causes of cancer-related death worldwide. Both genetic and environmental factors through a multitude of underlying molecular mechanisms participate in the pathogenesis of HCC. Recently, numerous studies have shown that circular RNAs (circRNAs), an emerging class of non-coding RNAs characterized by the presence of covalent bonds linking 3' and 5' ends, play an important role in the initiation and progression of cancers, including HCC. In this review, we outline the current status of the field of circRNAs, with an emphasis on the functions and mechanisms of circRNAs in HCC and its microenvironment. We also summarize and discuss recent advances of circRNAs as biomarkers and therapeutic targets. These efforts are anticipated to throw new insights into future perspectives about circRNAs in basic, translational and clinical research, eventually advancing the diagnosis, prevention and treatment of HCC.

Integrative Multiomics and Regulatory Network Analyses Uncovers the Role of OAS3, TRAFD1, miR-222-3p, and miR-125b-5p in Hepatitis E Virus Infection

The hepatitis E virus (HEV) is a long-ignored virus that has spread globally with time. It ranked 6th among the top risk-ranking viruses with high zoonotic spillover potential; thus, considering its viral threats is a pressing priority. The molecular pathophysiology of HEV infection or the underlying cause is limited. Therefore, we incorporated an unbiased, systematic methodology to get insights into the biological heterogeneity associated with the HEV. Our study fetched 93 and 2016 differentially expressed genes (DEGs) from chronic HEV (CHEV) infection in kidney-transplant patients, followed by hub module selection from a weighted gene co-expression network (WGCN). Most of the hub genes identified in this study were associated with interferon (IFN) signaling pathways. Amongst the genes induced by IFNs, the 2'-5'-oligoadenylate synthase 3 (OAS3) protein was upregulated. Protein-protein interaction (PPI) modular, functional enrichment, and feed-forward loop (FFL) analyses led to the identification of two key miRNAs, i.e., miR-222-3p and miR-125b-5p, which showed a strong association with the OAS3 gene and TRAF-type zinc finger domain containing 1 (TRAFD1) transcription factor (TF) based on essential centrality measures. Further experimental studies are required to substantiate the significance of these FFL-associated genes and miRNAs with their respective functions in CHEV. To our knowledge, it is the first time that miR-222-3p has been described as a reference miRNA for use in CHEV sample analyses. In conclusion, our study has enlightened a few budding targets of HEV, which might help us understand the cellular and molecular pathways dysregulated in HEV through various factors. Thus, providing a novel insight into its pathophysiology and progression dynamics.

CircRNA/miRNA/mRNA axis participates in the progression of partial bladder outlet obstruction

BACKGROUND

More and more evidence showed that circRNA/miRNA/mRNA axis played a vital role in the pathogenesis of some diseases. However, the role of circRNA/miRNA/mRNA axis in partial bladder outlet obstruction (pBOO) remains unknown. Our study aimed to explore the complex regulatory mechanism of circRNA/miRNA/mRNA axis in pBOO.

METHODS

The pBOO rat model was established, and the bladder tissues were collected for mRNA sequencing. The differentially expressed mRNAs were analyzed by high-throughput sequencing, and the GO and KEGG analysis of the differentially expressed mRNAs were performed. Competing endogenous RNAs (ceRNAs) analysis identified the potential regulation function of circRNA/miRNA/mRNA axis in pBOO. qRT-PCR detected the expression of circRNA/miRNA/mRNA. miRanda software was performed to predict the relationship between circRNA and miRNA, miRNA and mRNA.

RESULTS

Compared with the sham group, a total of 571 mRNAs were differentially expressed in the pBOO group, of which 286 were up-regulated and 285 were down-regulated. GO analysis showed that the mRNAs were mainly involved in cellular process, single-organism process, and cell, etc. KEGG analysis showed that the enriched signaling pathways were metabolic pathways, cell adhesion molecules (CAMs), and HTLV-I infection, etc. Based on the previous transcriptome data and differentially expressed circRNAs, we drew the ceRNA network regulation diagram. qRT-PCR results confirmed that chr3:113195876|113197193/rno-miR-30c-1-3p/Gata4, chr1:126188351|126195625/rno-miR-153-5p/Diaph3, and chr9:81258380|81275269/rno-miR-135b-5p/Pigr axis may have ceRNA function. miRanda confirmed there have the binding sites of circRNA/miRNA/mRNA axis.

CONCLUSIONS

CircRNA/miRNA/mRNA axis was involved in the progression of pBOO. Our research on the circRNA/miRNA/mRNA axis revealed new pathogenesis and treatment strategies for pBOO.

ceRNAshiny: An Interactive R/Shiny App for Identification and Analysis of ceRNA Regulation

The competing endogenous RNA (ceRNA) network is a newly discovered post-transcriptional regulation that controls both physiological and pathological progresses. Increasing research studies have been pivoted on this theory to explore the function of novel non-coding RNAs, pseudogenes, circular RNAs, and messenger RNAs. Although there are several R packages or computational tools to analyze ceRNA networks, an urgent need for easy-to-use computational tools still remains to identify ceRNA regulation. Besides, the conventional tools were mainly devoted to investigating ceRNAs in malignancies instead of those in neurodegenerative diseases. To fill this gap, we developed ceRNAshiny, an interactive R/Shiny application, which integrates widely used computational methods and databases to provide and visualize the construction and analysis of the ceRNA network, including differential gene analysis and functional annotation. In addition, demo data in ceRNAshiny could provide ceRNA network analyses about neurodegenerative diseases such as Parkinson's disease. Overall, ceRNAshiny is a user-friendly application that benefits all researchers, especially those who lack an established bioinformatic pipeline and are interested in studying ceRNA networks.

Competitive Endogenous RNA Network Activates Host Immune Response in SARS-CoV-2-, panH1N1 (A/California/07/2009)-, and H7N9 (A/Shanghai/1/2013)-Infected Cells

The global outbreak of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is still ongoing, as is research on the molecular mechanisms underlying cellular infection by coronaviruses, with the hope of developing therapeutic agents against this pandemic. Other important respiratory viruses such as 2009 pandemic H1N1 and H7N9 avian influenza virus (AIV), influenza A viruses, are also responsible for a possible outbreak due to their respiratory susceptibility. However, the interaction of these viruses with host cells and the regulation of post-transcriptional genes remains unclear. In this study, we detected and analyzed the comparative transcriptome profiling of SARS-CoV-2, panH1N1 (A/California/07/2009), and H7N9 (A/Shanghai/1/2013) infected cells. The results showed that the commonly upregulated genes among the three groups were mainly involved in autophagy, pertussis, and tuberculosis, which indicated that autophagy plays an important role in viral pathogenicity. There are three groups of commonly downregulated genes involved in metabolic pathways. Notably, unlike panH1N1 and H7N9, SARS-CoV-2 infection can inhibit the m-TOR pathway and activate the p53 signaling pathway, which may be responsible for unique autophagy induction and cell apoptosis. Particularly, upregulated expression of IRF1 was found in SARS-CoV-2, panH1N1, and H7N9 infection. Further analysis showed SARS-CoV-2, panH1N1, and H7N9 infection-induced upregulation of lncRNA-34087.27 could serve as a competitive endogenous RNA to stabilize IRF1 mRNA by competitively binding with miR-302b-3p. This study provides new insights into the molecular mechanisms of influenza A virus and SARS-CoV-2 infection.

Competing endogenous RNA network mediated by circ_3205 in SARS-CoV-2 infected cells

Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) is a new member of the Betacoronaviridae family, responsible for the recent pandemic outbreak of COVID-19. To start exploring the molecular events that follow host cell infection, we queried VirusCircBase and identified a circular RNA (circRNA) predicted to be synthesized by SARS-CoV-2, circ_3205, which we used to probe: (i) a training cohort comprised of two pools of cells from three nasopharyngeal swabs of SARS-CoV-2 infected (positive) or uninfected (negative, UCs) individuals; (ii) a validation cohort made up of 12 positive and 3 negative samples. The expression of circRNAs, miRNAs and miRNA targets was assayed through real-time PCR. CircRNA-miRNA interactions were predicted by TarpMiR, Analysis of Common Targets for circular RNAs (ACT), and STarMir tools. Enrichment of the biological processes and the list of predicted miRNA targets were retrieved from DIANA miRPath v3.0. Our results showed that the predicted SARS-CoV-2 circ_3205 was expressed only in positive samples and its amount positively correlated with that of SARS-CoV-2 Spike (S) mRNA and the viral load (r values = 0.80952 and 0.84867, Spearman's correlation test, respectively). Human (hsa) miR-298 was predicted to interact with circ_3205 by all three predictive tools. KCNMB4 and PRKCE were predicted as hsa-miR-298 targets. Interestingly, the function of both is correlated with blood coagulation and immune response. KCNMB4 and PRKCE mRNAs were upregulated in positive samples as compared to UCs (6 and 8.1-fold, p values = 0.049 and 0.02, Student's t test, respectively) and their expression positively correlated with that of circ_3205 (r values = 0.6 and 0.25, Spearman's correlation test, respectively). We propose that our results convincingly suggest that circ_3205 is a circRNA synthesized by SARS-CoV-2 upon host cell infection and that it may behave as a competitive endogenous RNA (ceRNA), sponging hsa-miR-298 and contributing to the upregulation of KCNMB4 and PRKCE mRNAs.

ViRBase v3.0: a virus and host ncRNA-associated interaction repository with increased coverage and annotation

As a means to aid in the investigation of viral infection mechanisms and identification of more effective antivirus targets, the availability of a source which continually collects and updates information on the virus and host ncRNA-associated interaction resources is essential. Here, we update the ViRBase database to version 3.0 (http://www.virbase.org/ or http://www.rna-society.org/virbase/). This update represents a major revision: (i) the total number of interaction entries is now greater than 820,000, an approximately 70-fold increment, involving 116 virus and 36 host organisms, (ii) it supplements and provides more details on RNA annotations (including RNA editing, RNA localization and RNA modification), ncRNA SNP and ncRNA-drug related information and (iii) it provides two additional tools for predicting binding sites (IntaRNA and PRIdictor), a visual plug-in to display interactions and a website which is optimized for more practical and user-friendly operation. Overall, ViRBase v3.0 provides a more comprehensive resource for virus and host ncRNA-associated interactions enabling researchers a more effective means for investigation of viral infections.

Illuminating lncRNA Function Through Target Prediction

Most of the transcribed human genome codes for noncoding RNAs (ncRNAs), and long noncoding RNAs (lncRNAs) make for the lion's share of the human ncRNA space. Despite growing interest in lncRNAs, because there are so many of them, and because of their tissue specialization and, often, lower abundance, their catalog remains incomplete and there are multiple ongoing efforts to improve it. Consequently, the number of human lncRNA genes may be lower than 10,000 or higher than 200,000. A key open challenge for lncRNA research, now that so many lncRNA species have been identified, is the characterization of lncRNA function and the interpretation of the roles of genetic and epigenetic alterations at their loci. After all, the most important human genes to catalog and study are those that contribute to important cellular functions-that affect development or cell differentiation and whose dysregulation may play a role in the genesis and progression of human diseases. Multiple efforts have used screens based on RNA-mediated interference (RNAi), antisense oligonucleotide (ASO), and CRISPR screens to identify the consequences of lncRNA dysregulation and predict lncRNA function in select contexts, but these approaches have unresolved scalability and accuracy challenges. Instead-as was the case for better-studied ncRNAs in the past-researchers often focus on characterizing lncRNA interactions and investigating their effects on genes and pathways with known functions. Here, we focus most of our review on computational methods to identify lncRNA interactions and to predict the effects of their alterations and dysregulation on human disease pathways.

Construction and Investigation of Competing Endogenous RNA Networks and Candidate Genes Involved in SARS-CoV-2 Infection

Introduction

The current COVID-19 pandemic caused by a novel coronavirus SARS-CoV-2 is a quickly developing global health crisis, yet the mechanisms of pathogenesis in COVID-19 are not fully understood.

Methods

The RNA sequencing data of SARS-CoV-2-infected cells was obtained from the Gene Expression Omnibus (GEO). The differentially expressed mRNAs (DEmRNAs), long non-coding RNAs (DElncRNAs), and microRNAs (DEmiRNAs) were identified by edgeR, and the SARS-CoV-2-associated competing endogenous RNA (ceRNA) network was constructed based on the prediction of bioinformatic databases. The Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses were conducted with the SARS-CoV-2-related DEmRNAs, and the protein–protein interaction network was also built basing on STRING database. The ROC analysis was performed for assessing the diagnostic efficiency of hub genes.

Results

The results indicated that SARS-CoV-2-related DEmRNAs were associated with the interferon signaling pathway and other antiviral processes, such as IFNL3, IFNL1 and CH25H. Our analysis suggested that lncRNA NEAT1 might regulate the host immune response through two miRNAs, hsa-miR-374-5p and hsa-miR-155-5p, which control the expression of SOCS1, IL6, IL1B, CSF1R, CD274, TLR6, and TNF. Additionally, IFI6, HRASLS2, IGFBP4 and PTN may be potential targets based on an analysis comparing the transcriptional responses of SARS-CoV-2 infection with that of other respiratory viruses.

Discussion

The unique ceRNA network identified potential non-coding RNAs and their possible targets as well as a new perspective to understand the molecular mechanisms of the host immune response to SARS-CoV-2. This study may also aid in the development of innovative diagnostic and therapeutic strategies.

Roles of ncRNAs as ceRNAs in Gastric Cancer

Although ignored in the past, with the recent deepening of research, significant progress has been made in the field of non-coding RNAs (ncRNAs). Accumulating evidence has revealed that microRNA (miRNA) response elements regulate RNA. Long ncRNAs, circular RNAs, pseudogenes, miRNAs, and messenger RNAs (mRNAs) form a competitive endogenous RNA (ceRNA) network that plays an essential role in cancer and cardiovascular, neurodegenerative, and autoimmune diseases. Gastric cancer (GC) is one of the most common cancers, with a high degree of malignancy. Considerable progress has been made in understanding the molecular mechanism and treatment of GC, but GC’s mortality rate is still high. Studies have shown a complex ceRNA crosstalk mechanism in GC. lncRNAs, circRNAs, and pseudogenes can interact with miRNAs to affect mRNA transcription. The study of the involvement of ceRNA in GC could improve our understanding of GC and lead to the identification of potential effective therapeutic targets. The research strategy for ceRNA is mainly to screen the different miRNAs, lncRNAs, circRNAs, pseudogenes, and mRNAs in each sample through microarray or sequencing technology, predict the ceRNA regulatory network, and, finally, conduct functional research on ceRNA. In this review, we briefly discuss the proposal and development of the ceRNA hypothesis and the biological function and principle of ceRNAs in GC, and briefly introduce the role of ncRNAs in the GC’s ceRNA network.

The bioinformatics toolbox for circRNA discovery and analysis

Circular RNAs (circRNAs) are a unique class of RNA molecule identified more than 40 years ago which are produced by a covalent linkage via back-splicing of linear RNA. Recent advances in sequencing technologies and bioinformatics tools have led directly to an ever-expanding field of types and biological functions of circRNAs. In parallel with technological developments, practical applications of circRNAs have arisen including their utilization as biomarkers of human disease. Currently, circRNA-associated bioinformatics tools can support projects including circRNA annotation, circRNA identification and network analysis of competing endogenous RNA (ceRNA). In this review, we collected about 100 circRNA-associated bioinformatics tools and summarized their current attributes and capabilities. We also performed network analysis and text mining on circRNA tool publications in order to reveal trends in their ongoing development.

Circular RNA in cancer development and immune regulation

Circular RNAs (circRNAs) are a class of single‐stranded RNAs with closed loop structures formed by covalent bonds of head and tail. Exploration of circRNAs is continually increasing; however, their functional relevance largely remains to be elucidated. In general, they are stable, abundant, conserved and expressed in tissue‐specific manner. These distinct properties and their diverse cellular actions indicate that circRNAs modulate transcription and translation, and may even function as translation templates. Growing evidence reveals that circRNAs contribute to various physiological and pathological processes, including the initiation and progression of cancer. In this review, we present the current knowledge about circRNAs in cancer development, as well as their potential for use as biomarkers and even therapeutic targets. CircRNA’s role in immune regulation and antitumour immunotherapy is also discussed. In addition, possible challenges in antitumour therapy are raised, and current progress and future perspectives are provided.

lncRNA CERS6-AS1 as ceRNA promote cell proliferation of breast cancer by sponging miR-125a-5p to upregulate BAP1 expression

Long noncoding RNAs (lncRNAs) can act as oncogene and tumor suppressor genes in many types of cancers including breast cancer (BC). Our previous study has indicated microRNA (miR)-125a-5p was downregulated and function as a tumor suppressor in BC. However, its upstream regulation mechanism is still unclear. In this study, we used bioinformatics algorithms, RNA pulldown assay, and dual-luciferase reports assay to predict and confirm lncRNA CERS6-AS1 interacted with miR-125a-5p. Then we found CERS6-AS1 was upregulated in BC tissues. Experimental results of tumor growth in nude mice show that CERS6-AS1 promotes tumor growth. Furthermore, CERS6-AS1 regulated BC susceptibility gene 1-associated protein 1 (BAP1) expression via sponging miR-125a-5p via Western blot analysis and quantitative polymerase chain reaction arrays. Finally, we showed that miR-125a-5p had opposing effects to those of CERS6-AS1 on BC cells, demonstrating that CERS6-AS1 may promote cell proliferation and inhibit cell apoptosis via sponging miR-125a-5p. Our results indicated CERS6-AS1 promote BC cell proliferation and inhibit cell apoptosis via sponging miR-125a-5p to upregulate BAP1 expression.

Unravelling host-pathogen interactions: ceRNA network in SARS-CoV-2 infection (COVID-19)

•

Bioinformatics analysis of mouse mRNA expression dataset for presumptive SARS-CoV-2 targets.

•

Induction of ISGs-Isg15, Oasl1, Usp18 and Ddx58 with no marked changes in the expression of IFNs.

•

No induction of ACE2 and TMPRSS2, raising implications for host factor limitations.

•

Identification of ceRNA network including miR-124-3p, Ddx58, lncRNA (Gm26917) and circRNAs (Ppp1r10, C330019G07RiK).

•

Virus regulates the expression of lnc and circRNAs, acting as sponges for miR-124-3p targeting Ddx58.

COVID-19 is a lurking calamitous disease caused by an unusual virus, SARS-CoV-2, causing massive deaths worldwide. Nonetheless, explicit therapeutic drugs or clinically approved vaccines are not available for COVID-19. Thus, a comprehensive research is crucially needed to decode the pathogenic tools, plausible drug targets, committed to the development of efficient therapy. Host-pathogen interactions via host cellular components is an emerging field of research in this respect. miRNAs have been established as vital players in host-virus interactions. Moreover, viruses have the capability to manoeuvre the host miRNA networks according to their own obligations. Besides protein coding mRNAs, noncoding RNAs might also be targeted in infected cells and viruses can exploit the host miRNA network via ceRNA effect. We have predicted a ceRNA network involving one miRNA (miR-124-3p), one mRNA (Ddx58), one lncRNA (Gm26917) and two circRNAs (Ppp1r10, C330019G07RiK) in SARS-CoV infected cells. We have identified 4 DEGs-Isg15, Ddx58, Oasl1, Usp18 by analyzing a mRNA GEO dataset. There is no notable induction of IFNs and IFN-induced ACE2, significant receptor responsible for S-protein binding mediated viral entry. Pathway enrichment and GO analysis conceded the enrichment of pathways associated with interferon signalling and antiviral-mechanism by IFN-stimulated genes. Further, we have identified 3 noncoding RNAs, playing as potential ceRNAs to the genes associated with immune mechanisms. This integrative analysis has identified noncoding RNAs and their plausible targets, which could effectively enhance the understanding of molecular mechanisms associated with viral infection. However, validation of these targets is further corroborated to determine their therapeutic efficacy.

RNAInter in 2020: RNA interactome repository with increased coverage and annotation

Research on RNA-associated interactions has exploded in recent years, and increasing numbers of studies are not limited to RNA-RNA and RNA-protein interactions but also include RNA-DNA/compound interactions. To facilitate the development of the interactome and promote understanding of the biological functions and molecular mechanisms of RNA, we updated RAID v2.0 to RNAInter (RNA Interactome Database), a repository for RNA-associated interactions that is freely accessible at http://www.rna-society.org/rnainter/ or http://www.rna-society.org/raid/. Compared to RAID v2.0, new features in RNAInter include (i) 8-fold more interaction data and 94 additional species; (ii) more definite annotations organized, including RNA editing/localization/modification/structure and homology interaction; (iii) advanced functions including fuzzy/batch search, interaction network and RNA dynamic expression and (iv) four embedded RNA interactome tools: RIscoper, IntaRNA, PRIdictor and DeepBind. Consequently, RNAInter contains >41 million RNA-associated interaction entries, involving more than 450 thousand unique molecules, including RNA, protein, DNA and compound. Overall, RNAInter provides a comprehensive RNA interactome resource for researchers and paves the way to investigate the regulatory landscape of cellular RNAs.

Long non-coding RNA and MicroRNA profiling provides comprehensive insight into non-coding RNA involved host immune responses in ALV-J-infected chicken primary macrophage

Avian leukosis virus subgroup J (ALV-J) infection can cause tumors and immunosuppression in infected chickens. Macrophages play a crucial role in host defense against invading pathogens. In the present study, whole transcriptome analysis was performed to analyze the host factors including genes, microRNA (miRNA), long non-coding RNA (lncRNA) and their regulatory network in chicken primary monocyte-derived macrophages (MDMs). In total, 128 differentially expressed (DE) lncRNAs and 15 DE miRNAs were identified in MDMs at 3 h post infection (hpi), and 30 DE lncRNAs and 8 DE miRNAs were identified in MDMs at 36 hpi during ALV-J infection. We further constructed the DE lncRNAs-mRNAs, miRNA-mRNAs and lncRNAs-miRNA-mRNAs interaction networks. The results suggested that DE lncRNAs and miRNAs are involved in the regulation of CCND3 and SOCS5 in Jak-STAT signaling pathway via ceRNA network in ALV-J-infected MDMs at 3 hpi. In addition, lncRNAs including XLOC_672329, ALDBGALG0000001429, XLOC_016500 and ALDBGALG0000000253 cis-regulating CH25H, CISH, IL-1β and CD80 respectively in MDMs at 3 hpi participated in host antiviral responses. Our findings give a comprehensive view of the connection between non-coding RNA and ALV-J in chicken primary macrophages, and provide an excellent resource for further studies of epigenetic effects on ALV-J disease resistance breeding as well as immune system and genomic researches.

Integration analysis for novel lncRNA markers predicting tumor recurrence in human colon adenocarcinoma

BACKGROUND

Numerous evidence has suggested that long non-coding RNA (lncRNA) acts an important role in tumor biology. This study focuses on the identification of novel prognostic lncRNA biomarkers predicting tumor recurrence in human colon adenocarcinoma.

METHODS

We obtained the research data from The Cancer Genome Atlas (TCGA) database. The interaction among different expressed lncRNA, miRNA and mRNA markers between colon adenocarcinoma patients with and without tumor recurrence were verified with miRcode, starBase and miRTarBase databases. We established the lncRNA-miRNA-mRNA competing endogenous RNA (ceRNA) network based on the verified association between the selected markers. We performed the functional enrichment analysis to obtain better understanding of the selected lncRNAs. Then we use multivariate logistic regression to identify the prognostic lncRNA markers with covariates. We also generated a nomogram predicting tumor recurrence risk based on the identified lncRNA biomarkers and clinical covariates.

RESULTS

We included 12,727 lncRNA, 1881 miRNA and 47,761 mRNA profiling and clinical features for 113 colon adenocarcinoma patients obtained from the TCGA database. After filtration, we used 37 specific lncRNAs, 60 miRNAs and 148 mRNAs in the ceRNA network analysis. We identified five lncRNAs as prognostic lncRNA markers predicting tumor recurrence in colon adenocarcinoma, in which four of them were identified for the first time. Finally, we generated a nomogram illustrating the association between the identified lncRNAs and the tumor recurrence risk in colon adenocarcinoma.

CONCLUSIONS

The four newly identified lncRNA biomarkers might be potential prognostic biomarkers predicting tumor recurrence in colon adenocarcinoma. We recommend that further clinical and fundamental researches be conducted on the identified lncRNA markers.

Cigarette smoke extract alters genome-wide profiles of circular RNAs and mRNAs in primary human small airway epithelial cells

As a novel kind of non-coding RNA, circular RNAs (circRNAs) were involved in various biological processes. However, the role of circRNAs in the developmental process of chronic obstructive pulmonary disease (COPD) is still unclear. In the present study, by using a cell model of COPD in primary human small airway epithelial cells (HSAECs) treated with or without cigarette smoke extract (CSE), we uncovered 4,379 previously unknown circRNAs in human cells and 903 smoke-specific circRNAs, with the help of RNA-sequencing and bioinformatic analysis. Moreover, 3,872 up- and 4,425 down-regulated mRNAs were also identified under CSE stimulation. Furthermore, a putative circRNA-microRNA-mRNA network was constructed for in-depth mechanism exploration, which indicated that differentially expressed circRNAs could influence expression of some key genes that participate in response to pentose phosphate pathway, ATP-binding cassette (ABC) transporters, glycosaminoglycan biosynthesis pathway and cancer-related pathways. Our research indicated that cigarette smoke had an influence on the biogenesis of circRNAs and mRNAs. CircRNAs might be involved in the response to CSE in COPD through the circRNA-mediated ceRNA networks.

Circular RNAs as a novel layer of regulatory mechanism in multiple sclerosis

Multiple sclerosis (MS) is believed to be an autoimmune disease of the central nervous system (CNS) in which autoreactive immune cells recognizing myelin antigens lead to demyelination and axonal injury. Mechanisms inducing and controlling the pathogenesis of MS have not been fully elucidated. Recent studies suggest an important role of epigenetic processes during the development of MS. One of the most significant discoveries in the field of epigenetic contribution to immune response has been the recognition of a group of microRNAs (miRNAs). These single-stranded non-coding RNA molecules regulate the expression of genes encoding proteins and have already been shown to be involved in pathogenesis of MS. Some miRNAs enhance generation of pro-inflammatory immune cells by promoting Th1 and Th17 pathways and others contribute to regulatory and tissue repair processes. The miRNA-dependent controlling process of autoimmune reactions is highly complex because of miRNA redundancy and multitarget nature of most of these molecules. Recently it was discovered that circular RNAs (circRNA) representing a new class of RNA possess a unique ability to control miRNAs by blocking their activity. CircRNAs are called natural miRNA "sponges" as the single circRNA molecule is able to neutralize several miRNAs and thus might determine the availability of miRNAs for their posttranscription regulation. Thus, circRNAs emerged as critical factors in epigenetic regulation of many human diseases including MS. In addition, in contrary to other RNA species they are very stable in the blood and other biological fluids and thus might be considered as a candidate for a biomarker of MS.

MicroRNA and Nonsense Transcripts as Putative Viral Evasion Mechanisms

Viral proteins encode numerous antiviral activities to modify the host immunity. In this article, we hypothesize that viral genomes and gene transcripts interfere with host gene expression using passive mechanisms to deregulate host microRNA (miRNA) activity. We postulate that various RNA viruses mimic or block binding between a host miRNA and its target transcript, a phenomenon mediated by the miRNA seed site at the 5′ end of miRNA. Virus-encoded miRNA seed sponges (vSSs) can potentially bind to host miRNA seed sites and prevent interaction with their native targets thereby relieving native miRNA suppression. In contrast, virus-encoded miRNA seed mimics (vSMs) may mediate considerable downregulation of host miRNA activity. We analyzed genomes from diverse RNA viruses for vSS and vSM signatures and found an abundance of these motifs indicating that this may be a mechanism of deceiving host immunity. Employing respiratory syncytial virus and measles virus as models, we reveal that regions surrounding vSS or vSM motifs have features characteristics of pre-miRNA templates and show that RSV viral transcripts are processed into small RNAs that may behave as vSS or vSM effectors. These data suggest that complex molecular interactions likely occur at the host-virus interface. Identifying the mechanisms in the network of interactions between the host and viral transcripts can help uncover ways to improve vaccine efficacy, therapeutics, and potentially mitigate the adverse events that may be associated with some vaccines.

Comprehensive analysis of circRNAs expression profiles in different periods of MDBK cells infected with bovine viral diarrhea virus

CircRNAs play an important regulatory role in the regulation of disease. However, we have a limited understanding of the role of circRNAs in the host's complex protective and pathological mechanisms of BVDV infection. Transcriptome analysis of circRNAs in host cells after BVDV infection may allow us to understand the biological functions of circRNAs in the regulation of BVDV infection. Here, we identified a total of 19,118 circRNAs from the MBDK cells (at 12 h, 24 h, and 48 h post-infection) infected with BVDV by using RNA-seq technology. We confirmed several circRNAs using RT-PCR and DNA sequencing, and qRT-PCR analysis was performed to identify several circRNAs expression and circRNAs resistance to RNase R digestion. GO (Gene Ontology) and KEGG (Kyoto Encyclopedia of Genes and Genomes) enrichment analysis showed that the host genes of differentially expressed circRNAs were involved in the regulation of cell proliferation, apoptosis, cycle and viral infection related signaling pathways. These results indicate that circRNA in host cells plays a broad regulatory role after BVDV infection and provides a valuable resource for studying circRNA biology in host cells after BVDV infection.

Long noncoding RNA: a crosslink in biological regulatory network

Long noncoding RNAs (lncRNAs) had been defined as a novel class of functional RNAs longer than 200 nucleotides around a decade ago. It is widely acknowledged that lncRNAs play a significant role in regulation of gene expression, but the biological and molecular mechanisms are diverse and complex, and remain to be determined. Especially, the regulatory network of lncRNAs associated with other biological molecules is still a controversial matter, thus becoming a new frontier of the studies on transcriptome. Recent advance in high-throughput sequencing technologies and bioinformatics approaches may be an accelerator to lift the mysterious veil. In this review, we will outline well-known associations between lncRNAs and other biological molecules, demonstrate the diverse bioinformatics approaches applied in prediction and analysis of lncRNA interaction and perform a case study for lncRNA linc00460 to concretely decipher the lncRNA regulatory network.

Interplay Between Long Noncoding RNAs and MicroRNAs in Cancer

In the last decade noncoding RNAs (ncRNAs) have been extensively studied in several biological processes and human diseases including cancer. microRNAs (miRNAs) are the best-known class of ncRNAs. miRNAs are small ncRNAs of around 20-22 nucleotides (nt) and are crucial posttranscriptional regulators of protein coding genes. Recently, new classes of ncRNAs, longer than miRNAs have been discovered. Those include intergenic noncoding RNAs (lincRNAs) and circular RNAs (circRNAs). These novel types of ncRNAs opened a very exciting field in biology, leading researchers to discover new relationships between miRNAs and long noncoding RNAs (lncRNAs), which act together to control protein coding gene expression. One of these new discoveries led to the formulation of the "competing endogenous RNA (ceRNA) hypothesis." This hypothesis suggests that an lncRNA acts as a sponge for miRNAs reducing their expression and causing the upregulation of miRNA targets. In this chapter we first discuss some recent discoveries in this field showing the mutual regulation of miRNAs, lncRNAs, and protein-coding genes in cancer. We then discuss the general approaches for the study of ceRNAs and present in more detail a recent computational approach to explore the ability of lncRNAs to act as ceRNAs in human breast cancer that has been shown to be, among the others, the most precise and promising.

Circular RNAs and Their Emerging Roles in Immune Regulation

Circular ribonucleic acid (RNA) molecules (circRNAs) are covalently closed loop RNA molecules with no 5' end caps or 3' poly (A) tails, which are generated by back-splicing. Originally, circRNAs were considered to be byproducts of aberrant splicing. However, in recent years, development of high-throughput sequencing has led to gradual recognition of functional circRNAs, and increasing numbers of studies have elucidated their roles in cancer, neurologic diseases, and cardiovascular disorders. Nevertheless, studies of the functions of circRNAs in the immune system are relatively scarce. In this review, we detail relevant research on the biogenesis and classification of circRNAs, describe their functional mechanisms and approaches to their investigation, and summarize recent studies of circRNA function in the immune system.

Potential regulatory role of circular RNA in idiopathic pulmonary fibrosis

Idiopathic pulmonary fibrosis (IPF) is a chronic, progressive type of interstitial pneumonia with unknown causes, poor prognosis and no effective therapy available. Circular RNAs (circRNAs), which serve as potential therapeutic targets and diagnostic biomarkers for certain diseases, represent a recent hotspot in the field of RNA research. In the present study, a total of 67 significantly dysregulated circRNAs were identified in the plasma of IPF patients by using a circRNA microarray. Among these circRNAs, 38 were upregulated, whereas 29 were downregulated. Further validation of the results by polymerase chain reaction analysis indicated that Homo sapiens (hsa)_circRNA_100906, hsa_circRNA_102100 and hsa_circRNA_102348 were significantly upregulated, whereas hsa_circRNA_101225, hsa_circRNA_104780 and hsa_circRNA_101242 were downregulated in plasma samples of IPF patients compared with those in samples from healthy controls. The majority of differentially expressed circRNAs were generated from exonic regions. The host genes of the differentially expressed circRNAs were involved in the regulation of the cell cycle, adherens junctions and RNA transport. The competing endogenous RNA (ceRNA) network of the circRNAs/micro(mi)RNAs/mRNAs indicated that circRNA-protected mRNA participated in transforming growth factor-β1, hypoxia-inducible factor-1, Wnt, Janus kinase, Rho-associated protein kinase, vascular endothelial growth factor, mitogen-activated protein kinase, Hedgehog and nuclear factor κB signalling pathways or functioned as biomarkers for pulmonary fibrosis. Furthermore, luciferase reporter assays confirmed that hsa_circRNA_100906 and hsa_circRNA_102348 directly interact with miR-324-5p and miR-630, respectively, which were downregulated in IPF patients. The present study provided a novel avenue for exploring the underlying molecular mechanisms of IPF disease.

Heat stress alters genome-wide profiles of circular RNAs in Arabidopsis

1599 novel circRNAs and 1583 heat stress-specific circRNAs were identified in Arabidopsis. Heat stress enhanced accumulation of circRNAs remarkably. Heat stress altered the sizes of circRNAs, numbers of circularized exons and alterative circularization events. A putative circRNA-mediated ceRNA networks under heat stress was established. Heat stress retards plant growth and destabilizes crop yield. The noncoding RNAs were demonstrated to be involved in plant response to heat stress. As a newly-characterized class of noncoding RNAs, circular RNAs (circRNAs) play important roles in transcriptional and post-transcriptional regulation. A few recent investigations indicated that plant circRNAs were differentially expressed under abiotic stress. However, little is known about how heat stress mediates biogenesis of circRNAs in plants. Here, we uncovered 1599 previously-unknown circRNAs and 1583 heat-specific circRNAs, by RNA-sequencing and bioinformatic analysis. Our results indicated that much more circRNAs were expressed under heat stress than in control condition. Besides, heat stress also increased the length of circRNAs, the quantity of circularized exons, and alternative circularization events. Moreover, we observed a positive correlation between expression patterns of some circRNAs and their parental genes. The prediction of ceRNA (competing endogenous RNA) networks indicated that differentially-expressed circRNAs could influence expression of many important genes, that participate in response to heat stress, hydrogen peroxide, and phytohormone signaling pathways, by interacting with the corresponding microRNAs. Together, our observations indicated that heat stress had great impacts on the biogenesis of circRNAs. Heat-induced circRNAs might participate in plant response to heat stress through the circRNA-mediated ceRNA networks.

Competitive endogenous RNA in colorectal cancer: A systematic review

Colorectal cancer is one of the most common malignant tumours. Competitive endogenous RNA (ceRNA) networks have been hypothesized, in which various RNAs regulate each other's expression using microRNA response elements (MREs). Recent evidence has highlighted the crucial regulatory roles of ceRNA networks in colorectal cancer. In this review, we summarize the present research methods as well as the currently known ceRNA competitors and targets in colorectal cancer. In addition, we discuss the significance of ceRNA and shortcomings of current studies of colorectal cancer.

Genome-wide analysis of circular RNAs in prenatal and postnatal muscle of sheep

Circular RNAs (circRNAs), a type of non-coding RNA with circular structure, were generated by back splicing and widely expressed in animals and plants. Recent studies have shown that circRNAs extensively participate in cell proliferation, cell differentiation, cell autophagy and other biological processes. However, the role and expression of circRNAs in the development and growth of muscle have not been studied in sheep. In our study, we first used RNA-seq to study the circRNAs in prenatal and postnatal longissimus dorsi muscle of sheep. A total of 6113 circRNAs were detected from the RNA-seq data. Several circRNAs were identified using reverse transcription PCR, DNA sequencing and RNase R digestion experiments. The expression levels of several circRNAs in prenatal and postnatal muscle were confirmed by Real-Time RT-PCR. The gene ontology (GO) and KEGG enrichment analysis of the host gene of the circRNAs showed that these circRNAs were mainly involved in the growth and development of muscle related signaling pathways. These circRNAs might sponge microRNAs (miRNAs) in predicted circRNA-miRNA-mRNA networks. The circRNAs expression profiles in muscle provided an important reference for the study of circRNAs in sheep.

Circular RNAs: A novel type of non-coding RNA and their potential implications in antiviral immunity

Circular RNAs (circRNAs), a novel type of non-coding RNAs (ncRNAs), are ubiquitously expressed in eukaryotic cells during post-transcriptional processes. Unlike linear RNAs, circRNAs form covalent-closed continuous loops without 5' to 3' polarities and poly (A) tails. With advances in high-throughput sequencing technology, numerous circRNAs have been identified in plants, animals and humans. Notably, circRNAs display cell-type, tissue-type and developmental-stage specific expression patterns in eukaryotic transcriptome, which reveals their significant regulatory functions in gene expression. More importantly, circRNAs serve as microRNA (miRNA) sponges and crucial regulators of gene expression. Additionally, circRNAs modulate pre-mRNA alternative splicing and possess protein-coding capacity. CircRNAs exhibit altered expression under pathological conditions and are strongly associated with the development of various human diseases. Interestingly, circRNAs can also induce antiviral immune responses. A recent study found that the delivery of circRNAs generated in vitro activates RIG-I-mediated innate immune responses and provides protection against viral infection. The antiviral dsRNA-binding proteins, NF90/NF110, act as key regulators in circRNA biogenesis. NF90/NF110 are also functional in inhibiting viral replication through binding to viral mRNAs. In this review, we provide a comprehensive overview on the classification, biogenesis and functions of circRNAs. We also discuss the critical role of circRNAs in eliciting antiviral immunity, providing evidence for the potential implications of circRNAs in antiviral therapies.

microRNA inhibitors: Natural and artificial sequestration of microRNA

MicroRNA (miRNAs) is post-transcriptional regulator of mRNA. However, the prevalence and activity of miRNA are regulated by other regulators. miRNA inhibitors are natural or artificial RNA transcripts that sequestrate miRNAs and decrease or even eliminate miRNA activity. Competing endogenous RNAs (ceRNAs) are natural and intracellular miRNA inhibitors that compete to bind to shared miRNA recognition elements (MREs) to decrease microRNA availability and relieve the repression of target RNAs. In recent years, studies have revealed that ceRNA crosstalk is involved in many pathophysiological processes and adds a new dimension to miRNA regulation. Artificial miRNA inhibitors are RNA transcripts that are synthesized via chemical and genetic methods. Artificial miRNA inhibitors can be used in miRNA loss-of-function research and gene therapies for certain diseases. In this review, we summarize the recent advances in the two different types of miRNA inhibitors.

Machine learning-based identification of endogenous cellular microRNA sponges against viral microRNAs

A "miRNA sponge" is an artificial oligonucleotide-based miRNA inhibitor containing multiple binding sites for a specific miRNA. Each miRNA sponge can bind and sequester several miRNA copies, thereby decreasing the cellular levels of the target miRNA. In addition to developing artificial miRNA sponges, scientists have sought endogenous RNA transcripts and found that long non-coding RNAs, competing endogenous RNAs, pseudogenes, circular RNAs, and coding RNAs could act as miRNA sponges under precise conditions. Here we present a computational approach for the prediction of endogenous human miRNA sponge candidates targeting viral miRNAs derived from pathogenic human viruses. Viral miRNA binding sites were predicted using a newly-developed machine learning-based method, and candidate interactions between miRNAs and sponge RNAs were experimentally validated using luciferase reporter assay, western blot analysis, and flow cytometry. We found that BX649188.1 functions as a potential natural miRNA sponge against kshv-miR-K12-7-3p.

The long non-coding RNA, SNHG6-003, functions as a competing endogenous RNA to promote the progression of hepatocellular carcinoma

The expression of long non-coding RNAs (lncRNAs) is dysregulated in hepatocellular carcinoma (HCC). However, the functions and contributions of lncRNAs remain largely unknown. Here, we identified a critical role of SNHG6-003 in HCC. We found that five SNHG6 transcripts were differentially expressed in HCC tissues while only the SNHG6-003 had an oncogenic function. Ectopic expression of SNHG6-003 in HCC cells promoted cell proliferation and induced drug resistance, whereas SNHG6-003 knockdown promoted apoptosis. Moreover, SNHG6-003 functioned as a competitive endogenous RNA (ceRNA), effectively becoming sponge for miR-26a/b and thereby modulating the expression of transforming growth factor-β-activated kinase 1 (TAK1). Importantly, expression analysis revealed that both SNHG6-003 and TAK1 were upregulated in human cancers, exhibiting a co-expression pattern. In HCC patients, high expression of SNHG6-003 closely correlated with tumor progression and shorter survival. Thus, targeting the ceRNA network involving SNHG6-003 may be used as a treatment strategy against HCC.

Endogenous microRNA sponges: evidence and controversy

The competitive endogenous RNA (ceRNA) hypothesis proposes that transcripts with shared microRNA (miRNA) binding sites compete for post-transcriptional control. This hypothesis has gained substantial attention as a unifying function for long non-coding RNAs, pseudogene transcripts and circular RNAs, as well as an alternative function for messenger RNAs. Empirical evidence supporting the hypothesis is accumulating but not without attracting scepticism. Recent studies that model transcriptome-wide binding-site abundance suggest that physiological changes in expression of most individual transcripts will not compromise miRNA activity. In this Review, we critically evaluate the evidence for and against the ceRNA hypothesis to assess the impact of endogenous miRNA-sponge interactions.

Dynamic changes in global microRNAome and transcriptome reveal complex miRNA-mRNA regulated host response to Japanese Encephalitis Virus in microglial cells

Microglia cells in the brain play essential role during Japanese Encephalitis Virus (JEV) infection and may lead to change in microRNA (miRNA) and mRNA profile. These changes may together control disease outcome. Using Affymetrix microarray platform, we profiled cellular miRNA and mRNA expression at multiple time points during viral infection in human microglial (CHME3) cells. In silico analysis of microarray data revealed a phased pattern of miRNAs expression, associated with JEV replication and provided unique signatures of infection. Target prediction and pathway enrichment analysis identified anti correlation between differentially expressed miRNA and the gene expression at multiple time point which ultimately affected diverse signaling pathways including Notch signaling pathways in microglia. Activation of Notch pathway during JEV infection was demonstrated in vitro and in vivo. The expression of a subset of miRNAs that target multiple genes in Notch signaling pathways were suppressed and their overexpression could affect JEV induced immune response. Further analysis provided evidence for the possible presence of cellular competing endogenous RNA (ceRNA) associated with innate immune response. Collectively, our data provide a uniquely comprehensive view of the changes in the host miRNAs induced by JEV during cellular infection and identify Notch pathway in modulating microglia mediated inflammation.

Identification of lncRNA-associated competing triplets reveals global patterns and prognostic markers for cancer

Recent studies have suggested that long non-coding RNAs (lncRNAs) can interact with microRNAs (miRNAs) and indirectly regulate miRNA targets though competing interactions. However, the molecular mechanisms underlying these interactions are still largely unknown. In this study, these lncRNA–miRNA–gene interactions were defined as lncRNA-associated competing triplets (LncACTs), and an integrated pipeline was developed to identify lncACTs that are active in cancer. Competing lncRNAs had sponge features distinct from non-competing lncRNAs. In the lncACT cross-talk network, disease-associated lncRNAs, miRNAs and coding-genes showed specific topological patterns indicative of their competence and control of communication within the network. The construction of global competing activity profiles revealed that lncACTs had high activity specific to cancers. Analyses of clustered lncACTs revealed that they were enriched in various cancer-related biological processes. Based on the global cross-talk network and cluster analyses, nine cancer-specific sub-networks were constructed. H19- and BRCA1/2-associated lncACTs were able to discriminate between two groups of patients with different clinical outcomes. Disease-associated lncACTs also showed variable competing patterns across normal and cancer patient samples. In summary, this study uncovered and systematically characterized global properties of human lncACTs that may have prognostic value for predicting clinical outcome in cancer patients.