CeRNASeek: an R package for identification and analysis of ceRNA regulation

Competing endogenous RNAs in head and neck squamous cell carcinoma: a review

Our understanding of RNA biology has evolved with recent advances in research from it being a non-functional product to molecules of the genome with specific regulatory functions. Competitive endogenous RNA (ceRNA), which has gained prominence over time as an essential part of post-transcriptional regulatory mechanism, is one such example. The ceRNA biology hypothesis states that coding RNA and non-coding RNA co-regulate each other using microRNA (miRNA) response elements. The ceRNA components include long non-coding RNAs, pseudogene and circular RNAs that exert their effect by interacting with miRNA and regulate the expression level of its target genes. Emerging evidence has revealed that the dysregulation of the ceRNA network is attributed to the pathogenesis of various cancers, including the head and neck squamous cell carcinoma (HNSCC). This is the most prevalent cancer developed from the mucosal epithelium in the lip, oral cavity, larynx and pharynx. Although many efforts have been made to comprehend the cause and subsequent treatment of HNSCC, the morbidity and mortality rate remains high. Hence, there is an urgent need to understand the holistic progression of HNSCC, mediated by ceRNA, that can have immense relevance in identifying novel biomarkers with a defined therapeutic intervention. In this review, we have made an effort to highlight the ceRNA biology hypothesis with a focus on its involvement in the progression of HNSCC. For the identification of such ceRNAs, we have additionally highlighted a number of databases and tools.

Construction of ceRNA network and key gene screening in cervical squamous intraepithelial lesions

BACKGROUND

This study aimed to construct an endogenous competition network for cervical squamous intraepithelial lesions using differential gene screening.

METHODS

GSE149763 was used to screen differentially expressed long non-coding RNAs (lncRNAs) and mRNAs to predict correlated microRNAs (miRNAs). The correlated miRNAs and GSE105409 were used to screen differentially expressed miRNAs for differential co-expression analysis, and the co-expressed differentially expressed miRNAs were used to predict correlated mRNAs. Differentially expressed mRNAs, miRNAs, and lncRNAs were visualized, and differential gene screening, enrichment, and pathway analysis were performed.

RESULTS

The ceRNA network of cervical squamous intraepithelial was successfully established and a potential differentially expressed network was identified. The key genes were VEGFA and FOS, and the key pathway was the MAPK signaling pathway.

CONCLUSIONS

The differential expression and potential effects of the lncRNA BACH1-IT1/miR-140-5p/VEGFA axis, key genes, VEGFA and FOS, and MAPK signaling in CIN were clarified, and the occurrence and potential effects of CIN were further clarified. The underlying molecular mechanism provides a certain degree of reference for subsequent treatments and experimental research.

miRspongeR 2.0: an enhanced R package for exploring miRNA sponge regulation

Summary

MicroRNA (miRNA) sponges influence the capability of miRNA-mediated gene silencing by competing for shared miRNA response elements and play significant roles in many physiological and pathological processes. It has been proved that computational or dry-lab approaches are useful to guide wet-lab experiments for uncovering miRNA sponge regulation. However, all of the existing tools only allow the analysis of miRNA sponge regulation regarding a group of samples, rather than the miRNA sponge regulation unique to individual samples. Furthermore, most existing tools do not allow parallel computing for the fast identification of miRNA sponge regulation. Here, we present an enhanced version of our R/Bioconductor package, miRspongeR 2.0. Compared with the original version introduced in 2019, this package extends the resolution of miRNA sponge regulation from the multi-sample level to the single-sample level. Moreover, it supports the identification of miRNA sponge networks using parallel computing, and the construction of sample-sample correlation networks. It also provides more computational methods to infer miRNA sponge regulation and expands the ground truth for validation. With these new features, we anticipate that miRspongeR 2.0 will further accelerate the research on miRNA sponges with higher resolution and more utilities.

Availability and implementation

http://bioconductor.org/packages/miRspongeR/.

Supplementary information

Supplementary data are available at Bioinformatics Advances online.

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.

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.

A Whole Transcriptome Analysis in Peripheral Blood Suggests That Energy Metabolism and Inflammation Are Involved in Major Depressive Disorder

INTRODUCTION

Previous studies on transcriptional profiles suggested dysregulation of multiple RNA species in major depressive disorder (MDD). However, the interaction between different types of RNA was neglected. Therefore, integration of different RNA species in transcriptome analysis would be helpful for interpreting the functional readout of the transcriptome in MDD.

METHODS

A whole transcriptome sequencing were performed on the peripheral blood of 15 patients with MDD and 15 matched healthy controls (HCs). The differential expression of miRNAs, lncRNAs, circRNAs, and mRNAs was examined between MDD and HCs using empirical analysis of digital gene expression data in R (edgeR). Weighted correlation network analysis (WGCNA) was used to identify RNA co-expression modules associated with MDD. A ceRNA network was constructed for interpretation of interactions between different RNA species. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses were conducted to explore potential biological mechanisms associated with MDD.

RESULTS

Multiple RNAs and co-expression modules were identified to be significantly dysregulated in MDD compared to HCs. Based on the differential RNAs, a ceRNA network that were dysregulated in MDD were constructed. The pathway networks that related to oxidative phosphorylation and the chemokine signaling were found to be associated with MDD.

CONCLUSION

Our results suggested that the processes of energy metabolism and inflammation may be involved in the pathophysiology of MDD.

Identification of CDK2-Related Immune Forecast Model and ceRNA in Lung Adenocarcinoma, a Pan-Cancer Analysis

Background

Tumor microenvironment (TME) plays important roles in different cancers. Our study aimed to identify molecules with significant prognostic values and construct a relevant Nomogram, immune model, competing endogenous RNA (ceRNA) in lung adenocarcinoma (LUAD).

Methods

“GEO2R,” “limma” R packages were used to identify all differentially expressed mRNAs from Gene Expression Omnibus (GEO) and The Cancer Genome Atlas (TCGA) databases. Genes with P-value <0.01, LogFC>2 or <-2 were included for further analyses. The function analysis of 250 overlapping mRNAs was shown by DAVID and Metascape software. By UALCAN, Oncomine and R packages, we explored the expression levels, survival analyses of CDK2 in 33 cancers. “Survival,” “survminer,” “rms” R packages were used to construct a Nomogram model of age, gender, stage, T, M, N. Univariate and multivariate Cox regression were used to establish prognosis-related immune forecast model in LUAD. CeRNA network was constructed by various online databases. The Genomics of Drug Sensitivity in Cancer (GDSC) database was used to explore correlations between CDK2 expression and IC50 of anti-tumor drugs.

Results

A total of 250 differentially expressed genes (DEGs) were identified to participate in many cancer-related pathways, such as activation of immune response, cell adhesion, migration, P13K-AKT signaling pathway. The target molecule CDK2 had prognostic value for the survival of patients in LUAD (P = 5.8e-15). Through Oncomine, TIMER, UALCAN, PrognoScan databases, the expression level of CDK2 in LUAD was higher than normal tissues. Pan-cancer analysis revealed that the expression, stage and survival of CDK2 in 33 cancers, which were statistically significant. Through TISIDB database, we selected 13 immunodepressants, 21 immunostimulants associated with CDK2 and explored 48 genes related to these 34 immunomodulators in cBioProtal database (P < 0.05). Gene Set Enrichment Analysis (GSEA) and Metascape indicated that 49 mRNAs were involved in PUJANA ATM PCC NETWORK (ES = 0.557, P = 0, FDR = 0), SIGNAL TRANSDUCTION (ES = –0.459, P = 0, FDR = 0), immune system process, cell proliferation. Forest map and Nomogram model showed the prognosis of patients with LUAD (Log-Rank = 1.399e-08, Concordance Index = 0.7). Cox regression showed that four mRNAs (SIT1, SNAI3, ASB2, and CDK2) were used to construct the forecast model to predict the prognosis of patients (P < 0.05). LUAD patients were divided into two different risk groups (low and high) had a statistical significance (P = 6.223e-04). By “survival ROC” R package, the total risk score of this prognostic model was AUC = 0.729 (SIT1 = 0.484, SNAI3 = 0.485, ASB2 = 0.267, CDK2 = 0.579). CytoHubba selected ceRNA mechanism medicated by potential biomarkers, 6 lncRNAs-7miRNAs-CDK2. The expression of CDK2 was associated with IC50 of 89 antitumor drugs, and we showed the top 20 drugs with P < 0.05.

Conclusion

In conclusion, our study identified CDK2 related immune forecast model, Nomogram model, forest map, ceRNA network, IC50 of anti-tumor drugs, to predict the prognosis and guide targeted therapy for LUAD patients.

CeNet Omnibus: an R/Shiny application to the construction and analysis of competing endogenous RNA network

Background

The competing endogenous RNA (ceRNA) regulation is a newly discovered post-transcriptional regulation mechanism and plays significant roles in physiological and pathological progress. CeRNA networks provide global views to help understand the regulation of ceRNAs. CeRNA networks have been widely used to detect survival biomarkers, select candidate regulators of disease genes, and predict long noncoding RNA functions. However, there is no software platform to provide overall functions from the construction to analysis of ceRNA networks.

Results

To fill this gap, we introduce CeNet Omnibus, an R/Shiny application, which provides a unified framework for the construction and analysis of ceRNA network. CeNet Omnibus enables users to select multiple measurements, such as Pearson correlation coefficient (PCC), mutual information (MI), and liquid association (LA), to identify ceRNA pairs and construct ceRNA networks. Furthermore, CeNet Omnibus provides a one-stop solution to analyze the topological properties of ceRNA networks, detect modules, and perform gene enrichment analysis and survival analysis. CeNet Omnibus intends to cover comprehensiveness, high efficiency, high expandability, and user customizability, and it also offers a web-based user-friendly interface to users to obtain the output intuitionally.

Conclusion

CeNet Omnibus is a comprehensive platform for the construction and analysis of ceRNA networks. It is highly customizable and outputs the results in intuitive and interactive. We expect that CeNet Omnibus will assist researchers to understand the property of ceRNA networks and associated biological phenomena. CeNet Omnibus is an R/Shiny application based on the Shiny framework developed by RStudio. The R package and detailed tutorial are available on our GitHub page with the URL https://github.com/GaoLabXDU/CeNetOmnibus.

LncRNA AC010789.1 Promotes Colorectal Cancer Progression by Targeting MicroRNA-432-3p/ZEB1 Axis and the Wnt/β-Catenin Signaling Pathway

Accumulating literatures have indicated that long non-coding RNAs (lncRNAs) are crucial molecules in tumor progression in various human cancers, including colorectal cancer (CRC). However, the clinical significance and regulatory mechanism of a vast majority of lncRNAs in CRC remain to be determined. The current study aimed to explore the function and molecular mechanism of lncRNA AC010789.1 in CRC progression. AC010789.1 found to be overexpressed in CRC tissues and cells. High expression of AC010789.1 was associated with lymph node metastasis and poor prognosis. Moreover, AC010789.1 silencing inhibited proliferation, migration, invasion and epithelial-mesenchymal transition (EMT) in vitro as well as tumorigenesis and metastasis in vivo. Mechanistically, we demonstrated that repression of AC010789.1 promoted miR-432-3p expression, and miR-432-3p directly binds to ZEB1. We then proved the anti-tumor role of miR-432-3p in CRC, showing that the inhibitory effect of AC010789.1 knockdown on CRC cells was achieved by the upregulation of miR-432-3p but downregulation of ZEB1. We also established that silencing AC010789.1 suppressed the Wnt/β-catenin signaling pathway. However, this inhibitory effect was partially counteracted by inhibition of miR-432-3p. In summary, these results reveal that silencing AC010789.1 suppresses CRC progression via miR-432-3p-mediated ZEB1 downregulation and suppression of the Wnt/β-catenin signaling pathway, highlighting a potentially promising strategy for CRC treatment.