Identify the critical protein-coding genes and long noncoding RNAs in cardiac myxoma

Cardiac myxoma (CM) is the most common benign cardiac tumor which is mostly sporadic. Increasing evidence show that protein-coding genes (PCGs) and long noncoding RNAs (lncRNAs) play important roles in the pathology processes of multiple cancers. However, the functional roles and regulatory mechanisms of RNAs interaction in CM are still unclear. In this study, we investigated three pairs of surgically excised CM by high throughput sequencing and screened a set of PCGs and lncRNAs which were differentially expressed and could serve as expression markers in CM. By constructing protein-protein interactions (PPI) and lncRNA-mRNA coexpressing network, we screened out a CM-related hub lncRNA-mRNA modules, which were enriched in different pathways such as MAPK and TGF-beta whose imbalance were validated by q-PCR. In addition, we identified a specific dysregulated competing endogenous RNA (ceRNA) network in CM by integrating lncRNA-miRNA-mRNA interactions. These results will help us to understand the interaction mechanisms of RNAs in CM and provide novel PCGs and lncRNAs as potential therapeutic targets for CM.     

Integrated analysis of co‐expression and ceRNA network identifies five lncRNAs as prognostic markers for breast cancer

Long non‐coding RNAs (lncRNAs), which competitively bind miRNAs to regulate target mRNA expression in the competing endogenous RNAs (ceRNAs) network, have attracted increasing attention in breast cancer research. We aim to find more effective therapeutic targets and prognostic markers for breast cancer. LncRNA, mRNA and miRNA expression profiles of breast cancer were downloaded from TCGA database. We screened the top 5000 lncRNAs, top 5000 mRNAs and all miRNAs to perform weighted gene co‐expression network analysis. The correlation between modules and clinical information of breast cancer was identified by Pearson's correlation coefficient. Based on the most relevant modules, we constructed a ceRNA network of breast cancer. Additionally, the standard Kaplan‐Meier univariate curve analysis was adopted to identify the prognosis of lncRNAs. Ultimately, a total of 23 and 5 modules were generated in the lncRNAs/mRNAs and miRNAs co‐expression network, respectively. According to the Green module of lncRNAs/mRNAs and Blue module of miRNAs, our constructed ceRNA network consisted of 52 lncRNAs, 17miRNAs and 79 mRNAs. Through survival analysis, 5 lncRNAs (AL117190.1, COL4A2‐AS1, LINC00184, MEG3 and MIR22HG) were identified as crucial prognostic factors for patients with breast cancer. Taken together, we have identified five novel lncRNAs related to prognosis of breast cancer. Our study has contributed to the deeper understanding of the molecular mechanism of breast cancer and provided novel insights into the use of breast cancer drugs and prognosis.     

Long noncoding RNA MEG3 play an important role in osteosarcoma development through sponging microRNAs

More and more evidence indicate long noncoding RNAs (lncRNAs) as competing endogenous RNAs (ceRNAs) to indirectly regulate messenger RNAs (mRNAs) by acting as microRNA (miRNA) sponges, which represents a novel layer of gene regulation that plays a critical role in the development of cancers. However, functional roles and regulatory mechanisms of lncRNA-mediated ceRNAs network in osteosarcoma are still largely unknown. Here, we comprehensively compared the expression profiles of mRNAs, lncRNAs, and miRNAs between osteosarcoma and normal samples from the Gene Expression Omnibus (GEO) to elaborate related latent mechanisms. Two lncRNAs, ie, LINC01560 and MEG3, were identified to be aberrantly expressed. Importantly, MEG3 was considered as a promising diagnostic biomarker and therapeutic target for patients with osteosarcoma according to the Kaplan-Meier analysis of another independent osteosarcoma data set from the Cancer Genome Atlas (P = 0.05). Eventually, we successfully established a dysregulated lncRNA-related ceRNA network, including one osteosarcoma-specific lncRNA, three miRNAs and four mRNAs. In conclusion, this study should be beneficial for improving our understanding of the lncRNA-mediated ceRNA regulatory mechanisms in the pathogenesis of osteosarcoma and providing it with novel candidate diagnostic and therapeutic biomarkers.     

Identification of functional lncRNAs in atrial fibrillation by integrative analysis of the lncRNA-mRNA network based on competing endogenous RNAs hypothesis

A mounting body of evidence has suggested that long noncoding RNAs (lncRNAs) play critical roles in human diseases by acting as competing endogenous RNAs (ceRNAs). However, the functions and ceRNA mechanisms of lncRNAs in atrial fibrillation (AF) remain to date unclear. In this study, we constructed an AF-related lncRNA-mRNA network (AFLMN) based on ceRNA theory, by integrating probe reannotation pipeline and microRNA (miRNA)-target regulatory interactions. Two lncRNAs with central topological properties in the AFLMN were first obtained. By using bidirectional hierarchical clustering, we identified two modules containing four lncRNAs, which were significantly enriched in many known pathways of AF. To elucidate the ceRNA interactions in certain disease or normal condition, the dysregulated lncRNA-mRNA crosstalks in AF were further analyzed, and six hub lncRNAs were obtained from the network. Furthermore, random walk analysis of the AFLMN suggested that lncRNA RP11-296O14.3 may function importantly in the pathological process of AF. All these eight lncRNAs that were identified from previous steps (RP11-363E7.4, GAS5, RP11-410L14.2, HAGLR, RP11-421L21.3, RP11-111K18.2, HOTAIRM1, and RP11-296O14.3) exhibited a strong diagnostic power for AF. The results of our study provide new insights into the functional roles and regulatory mechanisms of lncRNAs in AF, and facilitate the discovery of novel diagnostic biomarkers or therapeutic targets.     

Analysis of competing endogenous RNA network identifies a poorly differentiated cancer-specific RNA signature for hepatocellular carcinoma

Plenty of evidence has suggested that long noncoding RNAs (lncRNAs) play a vital role in competing endogenous RNA (ceRNA) networks. Poorly differentiated hepatocellular carcinoma (PDHCC) is a malignant phenotype. This paper aimed to explore the effect and the underlying regulatory mechanism of lncRNAs on PDHCC as a kind of ceRNA. Additionally, prognosis prediction was assessed. A total of 943 messenger RNAs (mRNAs), 86 miRNAs, and 468 lncRNAs that were differentially expressed between 137 PDHCCs and 235 well-differentiated HCCs were identified. Thereafter, a ceRNA network related to the dysregulated lncRNAs was established according to bioinformatic analysis and included 29 lncRNAs, 9 miRNAs, and 96 mRNAs. RNA-related overall survival (OS) curves were determined using the Kaplan-Meier method. The lncRNA ARHGEF7-AS2 was markedly correlated with OS in HCC (P = .041). Moreover, Cox regression analysis revealed that patients with low ARHGEF7-AS2 expression were associated with notably shorter survival time (P = .038). In addition, the area under the curve values of the lncRNA signature for 1-, 3-, and 5-year survival were 0.806, 0.741, and 0.701, respectively. Furthermore, a lncRNA nomogram was established, and the C-index of the internal validation was 0.717. In vitro experiments were performed to demonstrate that silencing ARHGEF7-AS2 expression significantly promoted HCC cell proliferation and migration. Taken together, our findings shed more light on the ceRNA network related to lncRNAs in PDHCC, and ARHGEF7-AS2 may be used as an independent biomarker to predict the prognosis of HCC.     

LNRLMI: Linear neighbour representation for predicting lncRNA‐miRNA interactions

LncRNA and miRNA are key molecules in mechanism of competing endogenous RNAs(ceRNA), and their interactions have been discovered with important roles in gene regulation. As supplementary to the identification of lncRNA‐miRNA interactions from CLIP‐seq experiments, in silico prediction can select the most potential candidates for experimental validation. Although developing computational tool for predicting lncRNA‐miRNA interaction is of great importance for deciphering the ceRNA mechanism, little effort has been made towards this direction. In this paper, we propose an approach based on linear neighbour representation to predict lncRNA‐miRNA interactions (LNRLMI). Specifically, we first constructed a bipartite network by combining the known interaction network and similarities based on expression profiles of lncRNAs and miRNAs. Based on such a data integration, linear neighbour representation method was introduced to construct a prediction model. To evaluate the prediction performance of the proposed model, k‐fold cross validations were implemented. As a result, LNRLMI yielded the average AUCs of 0.8475 ± 0.0032, 0.8960 ± 0.0015 and 0.9069 ± 0.0014 on 2‐fold, 5‐fold and 10‐fold cross validation, respectively. A series of comparison experiments with other methods were also conducted, and the results showed that our method was feasible and effective to predict lncRNA‐miRNA interactions via a combination of different types of useful side information. It is anticipated that LNRLMI could be a useful tool for predicting non‐coding RNA regulation network that lncRNA and miRNA are involved in.     

Construction of lncRNA-associated ceRNA networks to identify prognostic lncRNA biomarkers for glioblastoma

Long noncoding RNAs (lncRNAs) serve as competitive endogenous RNAs (ceRNAs) that play significant regulatory roles in the pathogenesis of tumors. However, the role of lncRNAs, especially the lncRNA-related ceRNA regulatory network, in glioblastoma (GBM) has not been fully elucidated. The goal of the current study was to construct lncRNA-microRNA-mRNA-related ceRNA networks for further investigation of their mechanism of action in GBM. We downloaded data from The Cancer Genome Atlas (TCGA) and the Gene Expression Omnibus (GEO) databases and identified differential lncRNAs, microRNAs (miRNAs), and messenger RNAs (mRNAs) associated with GBM. A ceRNA network was constructed and analyzed to examine the relationship between lncRNAs and patients’ overall survival. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGGs) were used to analyze the related mRNAs to indirectly explain the mechanism of action of lncRNAs. The potential effective drugs for the treatment of GBM were identified using the connectivity map (CMap). After integrated analysis, we obtained a total of 210 differentially expressed lncRNAs, 90 differentially expressed miRNAs, and 2508 differentially expressed mRNAs (DEmRNAs) from the TCGA and GEO databases. Using these differential genes, we constructed a lncRNA-associated ceRNA network. Six lncRNAs in the ceRNA network were associated with the overall survival of patients with GBM. Through KEGG analysis, it was found that the DEmRNAs involved in the network are related to cancer-associated pathways, for instance, mitogen-activated protein kinase and Ras signaling pathways. CMap analysis revealed four small-molecule compounds that could be used as drugs for the treatment of GBM. In this study, a multi-database joint analysis was used to construct a lncRNA-related ceRNA network to help identify the regulatory functions of lncRNAs in the pathogenesis of GBM.     

Comprehensive analysis of a ceRNA network reveals potential prognostic cytoplasmic lncRNAs involved in HCC progression

The aberrant expression of long noncoding RNAs (lncRNAs) has drawn increasing attention in the field of hepatocellular carcinoma (HCC) biology. In the present study, we obtained the expression profiles of lncRNAs, microRNAs (miRNAs), and messenger RNAs (mRNAs) in 371 HCC tissues and 50 normal tissues from The Cancer Genome Atlas (TCGA) and identified hepatocarcinogenesis-specific differentially expressed genes (DEGs, log fold change ≥ 2, FDR < 0.01), including 753 lncRNAs, 97 miRNAs, and 1,535 mRNAs. Because the specific functions of lncRNAs are closely related to their intracellular localizations and because the cytoplasm is the main location for competitive endogenous RNA (ceRNA) action, we analyzed not only the interactions among these DEGs but also the distributions of lncRNAs (cytoplasmic, nuclear or both). Then, an HCC-associated deregulated ceRNA network consisting of 37 lncRNAs, 10 miRNAs, and 26 mRNAs was constructed after excluding those lncRNAs located only in the nucleus. Survival analysis of this network demonstrated that 15 lncRNAs, 3 miRNAs, and 16 mRNAs were significantly correlated with the overall survival of HCC patients (p < 0.01). Through multivariate Cox regression and lasso analysis, a risk score system based on 13 lncRNAs was constructed, which showed good discrimination and predictive ability for HCC patient survival time. This ceRNA network-construction approach, based on lncRNA distribution, not only narrowed the scope of target lncRNAs but also provided specific candidate molecular biomarkers for evaluating the prognosis of HCC, which will help expand our understanding of the ceRNA mechanisms involved in the early development of HCC.     

Analysing the relationship between lncRNA and protein‐coding gene and the role of lncRNA as ceRNA in pulmonary fibrosis

Long non‐coding RNAs (lncRNAs) are involved in various pathophysiologic processes and human diseases. However, their dynamics and corresponding functions in pulmonary fibrosis remain poorly understood. In this study, portions of lncRNAs adjacent or homologous to protein‐coding genes were determined by searching the UCSC genome bioinformatics database. This was found to be potentially useful for exploring lncRNA functions in disease progression. Previous studies showed that competing endogenous RNA (ceRNA) hypothesis is another method to predict lncRNA function. However, little is known about the function of ceRNA in pulmonary fibrosis. In this study, we selected two differentially expressed lncRNAs MRAK088388 and MRAK081523 to explore their regulatory mechanisms. MRAK088388 and MRAK081523 were analysed as long‐intergenic non‐coding RNAs (lincRNAs), and identified as orthologues of mouse lncRNAs AK088388 and AK081523, respectively. qRT‐PCR and in situ hybridization (ISH) showed that they were significantly up‐regulated, and located in the cytoplasm of interstitial lung cells. We also showed that MRAK088388 and N4bp2 had the same miRNA response elements (MREs) for miR‐200, miR‐429, miR‐29, and miR‐30, whereas MRAK081523 and Plxna4 had the same MREs for miR‐218, miR‐141, miR‐98, and let‐7. Moreover, the expression levels of N4bp2 and Plxna4 significantly increased in fibrotic rats, and were highly correlated with those of MRAK088388 and MRAK081523, respectively. Among their shared miRNAs, miR‐29b‐3p and let‐7i‐5p decreased in the model group, and were negatively correlated with the expression of MRAK088388 and MRAK081523, respectively. MRAK088388 and MRAK081523 could regulate N4bp2 and Plxna4 expression by sponging miR‐29b‐3p and let‐7i‐5p, respectively, and possessed regulatory functions as ceRNAs. Thus, our study may provide insights into the functional interactions of lncRNA, miRNA and mRNA, and lead to new theories for the pathogenesis and treatment of pulmonary fibrosis.     

Cancer Specific Long Noncoding RNAs Show Differential Expression Patterns and Competing Endogenous RNA Potential in Hepatocellular Carcinoma

Long noncoding RNAs (lncRNAs) regulate gene expression by acting with microRNAs (miRNAs). However, the roles of cancer specific lncRNA and its related competitive endogenous RNAs (ceRNA) network in hepatocellular cell carcinoma (HCC) are not fully understood. The lncRNA profiles in 372 HCC patients, including 372 tumor and 48 adjacent non-tumor liver tissues, from The Cancer Genome Atlas (TCGA) and NCBI GEO omnibus (GSE65485) were analyzed. Cancer specific lncRNAs (or HCC related lncRNAs) were identified and correlated with clinical features. Based on bioinformatics generated from miRcode, starBase, and miRTarBase, we constructed an lncRNA-miRNA-mRNA network (ceRNA network) in HCC. We found 177 cancer specific lncRNAs in HCC (fold change ≥ 1.5, P < 0.01), 41 of them were also discriminatively expressed with gender, race, tumor grade, AJCC tumor stage, and AJCC TNM staging system. Six lncRNAs (CECR7, LINC00346, MAPKAPK5-AS1, LOC338651, FLJ90757, and LOC283663) were found to be significantly associated with overall survival (OS, log-rank P < 0.05). Collectively, our results showed the lncRNA expression patterns and a complex ceRNA network in HCC, and identified a complex cancer specific ceRNA network, which includes 14 lncRNAs and 17 miRNAs in HCC.     

LncRNA‐PAGBC acts as a microRNA sponge and promotes gallbladder tumorigenesis

Long noncoding RNAs (lncRNAs) play roles in the development and progression of many cancers; however, the contributions of lncRNAs to human gallbladder cancer (GBC) remain largely unknown. In this study, we identify a group of differentially expressed lncRNAs in human GBC tissues, including prognosis‐associated gallbladder cancer lncRNA (lncRNA‐PAGBC), which we find to be an independent prognostic marker in GBC. Functional analysis indicates that lncRNA‐PAGBC promotes tumour growth and metastasis of GBC cells. More importantly, as a competitive endogenous RNA (ceRNA), lncRNA‐PAGBC competitively binds to the tumour suppressive microRNAs miR‐133b and miR‐511. This competitive role of lncRNA‐PAGBC is required for its ability to promote tumour growth and metastasis and to activate the AKT/mTOR pathway. Moreover, lncRNA‐PAGBC interacts with polyadenylate binding protein cytoplasmic 1 (PABPC1) and is stabilized by this interaction. This work provides novel insight on the molecular pathogenesis of GBC.     

Reconstruction and analysis of the aberrant lncRNA‐miRNA‐mRNA network based on competitive endogenous RNA in CESC

A growing body of studies has demonstrated that long non‐coding RNA (lncRNA) are regarded as the primary section of the ceRNA network. This is thought to be the case owing to its regulation of protein‐coding gene expression by functioning as miRNA sponges. However, functional roles and regulatory mechanisms of lncRNA‐mediated ceRNA in cervical squamous cell carcinoma (CESC), as well as their use for potential prediction of CESC prognosis, remains unknown. The aberrant expression profiles of mRNA, lncRNA, and miRNA of 306 cervical squamous cancer tissues and three adjacent cervical tissues were obtained from the TCGA database. A lncRNA‐mRNA‐miRNA ceRNA network in CESC was constructed. Meanwhile, Gene Ontology (GO) and KEGG pathway analysis were performed using Cytoscape plug‐in BinGo and DAVID database. We identified a total of 493 lncRNA, 70 miRNA, and 1921 mRNA as differentially expressed profiles. An aberrant lncRNA‐mRNA‐miRNA ceRNA network was constructed in CESC, it was composed of 50 DElncRNA, 18 DEmiRNA, and 81 DEmRNA. According to the overall survival analysis, 3 out of 50 lncRNA, 10 out of 81 mRNA, and 1 out of 18 miRNA functioned as prognostic biomarkers for patients with CESC (P value < 0.05). We extracted the sub‐network in the ceRNA network and found that two novel lncRNA were recognized as key genes. These included lncRNA MEG3 and lncRNA ADAMTS9‐AS2. The present study provides a new insight into a better understanding of the lncRNA‐related ceRNA network in CESC, and the novel recognized ceRNA network will help us to improve our understanding of lncRNA‐mediated ceRNA regulatory mechanisms in the pathogenesis of CESC.     

To construct a ceRNA regulatory network as prognostic biomarkers for bladder cancer

Emerging evidence demonstrates that competing endogenous RNA (ceRNA) hypothesis has played a role in molecular biological mechanisms of cancer occurrence and development. But the effect of ceRNA network in bladder cancer (BC), especially lncRNA‐miRNA‐mRNA regulatory network of BC, was not completely expounded. By means of The Cancer Genome Atlas (TCGA) database, we compared the expression of RNA sequencing (RNA‐Seq) data between 19 normal bladder tissue and 414 primary bladder tumours. Then, weighted gene co‐expression network analysis (WGCNA) was conducted to analyse the correlation between two sets of genes with traits. Interactions between miRNAs, lncRNAs and target mRNAs were predicted by MiRcode, miRDB, starBase, miRTarBase and TargetScan. Next, by univariate Cox regression and LASSO regression analysis, the 86 mRNAs obtained by prediction were used to construct a prognostic model which contained 4 mRNAs (ACTC1 + FAM129A + OSBPL10 + EPHA2). Then, by the 4 mRNAs in the prognostic model, a ceRNA regulatory network with 48 lncRNAs, 14 miRNAs and 4 mRNAs was constructed. To sum up, the ceRNA network can further explore gene regulation and predict the prognosis of BC patients.     

Expression and functional analysis of lncRNAs in the hippocampus of immature rats with status epilepticus

Long non‐coding RNAs (lncRNAs) have been implicated in the regulation of gene expression at various levels. However, to date, the expression profile of lncRNAs in status epilepticus (SE) was unclear. In our study, the expression profile of lncRNAs was investigated by high‐throughput sequencing based on a lithium/pilocarpine‐induced SE model in immature rats. Furthermore, weighted correlation network analysis (WGCNA), gene ontology (GO) analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis were performed to construct co‐expression networks and establish functions of the identified hub lncRNAs in SE. The functional role of a hub lncRNA (NONRATT010788.2) in SE was investigated in an in vitro model. Our results indicated that 7082 lncRNAs (3522 up‐regulated and 3560 down‐regulated), which are involved in cell proliferation, inflammatory responses, angiogenesis and autophagy, were dysregulated in the hippocampus of immature rats with SE. Additionally, WGCNA identified 667 up‐regulated hub lncRNAs in turquoise module that were involved in apoptosis, inflammatory responses and angiogenesis via regulation of HIF‐1, p53 and chemokine signalling pathways and via inflammatory mediator regulation of TRP channels. Knockdown of an identified hub lncRNA (NONRATT010788.2) inhibited neuronal apoptosis in vitro. Taken together, our study is the first to demonstrate the expression profile and potential function of lncRNAs in the hippocampus of immature rats with SE. The defined hub lncRNAs may participate in the pathogenesis of SE via lncRNA‐miRNA‐mRNA network.