子宫内膜癌中的lncRNA-miRNA-mRNA ceRNA网络的整合分析

子宫内膜癌(uterine corpus endometrial carcinoma, UCEC)是危害女性健康的癌症之一,但其发生发展机制尚不完全清楚。基因的异常表达在细胞癌变过程中发挥着重要的作用。本研究利用生物信息学方法对UCEC中异常表达的基因进行网络调控分析,为UCEC的机制研究及预后治疗提供理论依据。首先,利用"limma"包筛选得到差异表达的RNA分子;接着,利用R软件中的"GDCRNATools"包构建lncRNA-miRNA-mRNAceRNA网络;最后,利用DAVID(TheDatabasefor Annotation,Visualization and Integrated Discovery)、R软件中的"survival"包、MCODE(Molecular COmplex Detection)插件、StarBase数据库等对ceRNA网络中的RNA分子进行综合分析。结果显示,在UCEC中,我们鉴定了1 319个mRNA、68个lncRNA及100个miRNA为差异表达的RNA分子。同时,利用这些UCEC中的差异表达RNA分子,成功构建了ceRNA网络,结果表明,HCG11、LINC00958、LINC00667、MAGI2-AS3和AC093010.3可能具有ceRNA的功能。通过进一步对该网络成员的分析,结果发现,该网络中的mRNA分子聚集于内皮细胞增殖、细胞骨架和锌离子结合负调控等生物学过程中。生存分析显示,hsa-miR-449a、LINC00958、PKIA和DPYSL2等19个差异表达的RNA分子与UCEC患者预后显著相关。最后,利用MCODE插件对ceRNA网络进行筛选,共获得一个子网络,其成员分别是LINC00667、hsa-miR-449a、hsa-miR-34a-5p和RECK。该研究结果提示,在UCEC的发生发展中存在lncRNA、miRNA和mRNA的差异表达,并且这些RNA分子之间存在ceRNA网络模式的调控,其中一些关键分子与患者的预后相关,这将为进一步研究与理解UCEC的发病机制及预后提供理论依据。     

Long noncoding RNA LINC00958 regulates cell sensitivity to radiotherapy through RRM2 by binding to microRNA-5095 in cervical cancer

Long noncoding RNAs (lncRNAs) have been implicated in the regulation of resistance to radiotherapy in cervical cancer, which is a type of gynecological disease with high mortality in women around the world. Hence, our purpose is to delineate the involvement of LINC00958 in regulating cell sensitivity to radiotherapy in cervical cancer. LINC00958 expression in cervical cancer was assayed, followed by verification of the relationship among LINC00958, microRNA-5095 (miR-5095) and ribonucleotide reductase subunit M2 (RRM2). Hela cells were transduced with up-/downregulation of miR-5095 or RRM2, or LINC00958 silencing, respectively, and then treated with or without a 6 Gy dose of X-ray irradiation. Then the cell proliferation, apoptosis, survival fraction rate, as well as sensitivity to radiotherapy, were assessed. Finally, xenograft tumor in nude mice was established by transplanting Hela cells transfected with sh-LINC00958 and irradiated with 6 Gy of X-ray. High expression of LINC00958 was revealed in The Cancer Genome Atlas and Gene Expression Profiling Interactive Analysis, as well as in radiation-resistant patients, which was associated with lower sensitivity to radiotherapy in cervical cancer. Moreover, cervical cancer patients with higher LINC00958 expression exhibited a shorter overall survival according to Kaplan–Meier analysis. In addition, LINC00958 could regulate the expression of RRM2 by competing for miR-5095. A combination of radiotherapy with LINC00958 silencing, RRM2 downregulation or miR-5095 overexpression was found to inhibit cervical cancer cell proliferation and tumor growth, while promoting cell apoptosis both in vitro and in vivo. Collectively, our results suggest that LINC00958 could regulate RRM2 by competing to miR-5095, which regulates cell sensitivity to radiotherapy in cervical cancer.     

Long noncoding RNA HOXA-AS2 regulates the expression of SCN3A by sponging miR-106a in breast cancer

Breast cancer is the most commonly diagnosed cancer that affects women worldwide. This study aimed to investigate the competing endogenous RNAs (ceRNAs) mechanism in breast cancer. Microarray data were downloaded from the University of California Santa Cruz (UCSC) Xena database. The limma package was used to screen the differentially expressed messenger RNAs (DEMs) and differentially expressed long noncoding RNAs (DELs). Subsequently, functional analysis was performed using DAVID tool. After constructing the protein-protein interaction (PPI) network, we identified the major gene modules using the Cytoscape software. Univariate survival analysis in the survival package was performed. Finally, the ceRNA regulatory network was constructed to identify the critical genes. A total of 1380 DEMs and 345 DELs were identified in breast cancer samples compared with normal samples. Functional enrichment analysis showed that DEMs were mainly involved in cell division, and cell cycle. We screened four major gene modules and identified the hub nodes in these functional modules. Several DEMs (including FABP7, C4BPA, and LAMB3) and three long noncoding RNAs (lncRNAs) (LINC00092, SLC26A4.AS1, and COLCA1) exhibited significant correlation with patients' survival outcomes. In the ceRNA network, the lncRNA HOXA-AS2 regulated the expression level of SCN3A by interacting with hsa-miR-106a-5p. Thus, our study investigated the ceRNA mechanism in breast cancer. The results showed that lncRNA HOXA-AS2 might modulate the expression of SCN3A by sponging miR-106a in breast cancer.     

ce‐Subpathway: Identification of ceRNA‐mediated subpathways via joint power of ceRNAs and pathway topologies

Competing endogenous RNAs (ceRNAs) represent a novel mechanism of gene regulation that may mediate key subpathway regions and contribute to the altered activities of pathways. However, the classical methods used to identify pathways fail to specifically consider ceRNAs within the pathways and key regions impacted by them. We proposed a powerful strategy named ce‐Subpathway for the identification of ceRNA‐mediated functional subpathways. It provided an effective level of pathway analysis via integrating ceRNAs, differentially expressed (DE) genes and their key regions within the given pathways. We respectively analysed one pulmonary arterial hypertension (PAH) and one myocardial infarction (MI) data sets and demonstrated that ce‐Subpathway could identify many subpathways whose corresponding entire pathways were ignored by those non‐ceRNA‐mediated pathway identification methods. And these pathways have been well reported to be associated with PAH/MI‐related cardiovascular diseases. Further evidence showed reliability of ceRNA interactions and robustness/reproducibility of the ce‐Subpathway strategy by several data sets of different cancers, including breast cancer, oesophageal cancer and colon cancer. Survival analysis was finally applied to illustrate the clinical application value of the ceRNA‐mediated functional subpathways using another data sets of pancreatic cancer. Comprehensive analyses have shown the power of a joint ceRNAs/DE genes and subpathway strategy based on their topologies.     

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.     

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.     

LINC00958 facilitates cervical cancer cell proliferation and metastasis by sponging miR-625-5p to upregulate LRRC8E expression

Accepted as a malignant tumor worldwide, cervical cancer (CC) has attracted much attention for its high incidence and mortality rates. Previous studies have elucidated the critical regulatory function that long noncoding RNAs (lncRNAs) exert on the tumorigenesis and progression of diverse tumors. Although multiple investigations have depicted that LINC00958 has a great impact on the complex biological process of many cancers, knowledge concerning the regulatory role of LINC00958 in CC remains limited and needs to be further explored. In our study, LINC00958 expression was evidently overexpressed in CC tissues and cells. Besides this, LINC00958 negatively regulated miR-625-5p expression and was verified to bind with miR-625-5p in CC. Subsequently, it was testified by a series of experiments that LINC00958 promotes CC cell proliferation and metastasis by sponging miR-625-5p. Furthermore, the leucine-rich repeat containing the eight family member E (LRRC8E) could bind with miR-625-5p, and its expression was negatively modulated by miR-625-5p, whereas positively regulated by LINC00958 in CC. Final rescue assays verified the effects of LINC0095/LRRC8E interaction and miR-625-5p/LRRC8E interaction on CC cell proliferation and metastasis. Collectively, LINC00958 facilitates CC cell proliferation and metastasis via the miR-625-5p/LRRC8E axis.     

Identification of messenger and long noncoding RNAs associated with gallbladder cancer via gene expression profile analysis

The present study aimed to investigate the long noncoding RNAs (lncRNAs) and messenger RNAs (mRNAs) involved in the progression of gallbladder cancer and explore the potential physiopathologic mechanisms of gallbladder cancer in terms of competing endogenous RNAs (ceRNAs). The original lncRNA and mRNA expression profile data (nine gallbladder cancer tissues samples and nine normal gallbladder samples) in GSE76633 was downloaded from the Gene Expression Omnibus database. Differentially expressed mRNAs and lncRNAs between gallbladder cancer tissue and normal control were selected and the pathways in which they are involved were analyzed using bioinformatics analyses. MicroRNAs (miRNAs) were also predicted based on the differentially expressed mRNAs. Finally, the co-expression relation between lncRNA and mRNA was analyzed and the ceRNA network was constructed by combining the lncRNA-miRNA, miRNA-mRNA, and lncRNA-mRNA pairs. Overall, 373 significantly differentially expressed mRNAs and 47 lncRNAs were identified between cancer and normal tissue samples. The upregulated genes were significantly enriched in the extracellular matrix (ECM)-receptor interaction pathway, while the downregulated genes were involved in the complement and coagulation cascades. Altogether, 128 co-expression relations between lncRNA and mRNA were obtained. In addition, 196 miRNA-mRNA regulatory relations and 145 miRNA-lncRNA relation pairs were predicted. Finally, the lncRNA-miRNA-gene ceRNA network was constructed by combining the three types of relation pairs, such as XLOC_011309-miR-548c-3p-SPOCK1 and XLOC_012588-miR-765-CEACAM6. mRNAs and lncRNAs may be involved in gallbladder cancer progression via ECM-receptor interaction pathways and the complement and coagulation cascades. Moreover, ceRNAs such as XLOC_011309-miR-548c-3p-SPOCK1 and XLOC_012588-miR-765-CEACAM6 can also be implicated in the pathogenesis of gallbladder cancer.     

The mRNA related ceRNA–ceRNA landscape and significance across 20 major cancer types

Cross-talk between competitive endogenous RNAs (ceRNAs) through shared miRNAs represents a novel layer of gene regulation that plays important roles in the physiology and development of cancers. However, a global view of their system-level properties across various types of cancers is still unknown. Here, we constructed the mRNA related ceRNA–ceRNA interaction landscape across 20 cancer types by systematically analyzing molecular profiles of 5203 tumors and miRNA regulations. Our study highlights the conserved features shared by pan-cancer and higher similarity within similar origin cell type. Moreover, a core ceRNA network was identified. Function analysis identified a common theme of cancer hallmarks, however they exhibit phenotype-specific connectivity patterns. Besides, we found a marked rewiring in the ceRNA program between various cancers, and further revealed conserved and rewired network ceRNA hubs in each cancer, which were tensely competitive interactions to constitute conserved and cancer-specific modules. By providing mechanistic linkage between known cancer miRNAs, their mediated ceRNA–ceRNA interactions, and the associations with known cancer hallmarks, the inferred cancer ceRNA–ceRNA interaction landscape will serve as a powerful public resource for further biological discoveries of tumorigenesis.     

lncRNA作为竞争性内源RNA在非小细胞肺癌中的作用

长链非编码RNA(long non-coding RNA,lncRNA)参与肿瘤的多种生理、病理进程.研究表明,lncRNA可通过与微小RNA (microRNA, mi RNA)反应元件相互作用,并与其他RNA分子形成竞争性内源RNA (competing endogenous RNA,ceRNA)的调控网络,参与基因的表达调控.lncRNA以ceRNA方式参与非小细胞肺癌(non-small cell lung cancer,NSCLC)的发生发展过程,为揭示NSCLC的分子机理开拓了新的思路,也为NSCLC的治疗提供新的靶点.本文在课题组前期发现NSCLC相关ceRNA基础上,主要讨论lncRNA作为ceRNA在NSCLC中高表达、低表达及治疗相关方面的作用.     

Long noncoding RNA LINC00941 promotes pancreatic cancer progression by competitively binding miR-335-5p to regulate ROCK1-mediated LIMK1/Cofilin-1 signaling

An accumulation of evidence indicates that long noncoding RNAs are involved in the tumorigenesis and progression of pancreatic cancer (PC). In this study, we investigated the functions and molecular mechanism of action of LINC00941 in PC. Quantitative PCR was used to examine the expression of LINC00941 and miR-335-5p in PC tissues and cell lines, and to investigate the correlation between LINC00941 expression and clinicopathological features. Plasmid vectors or lentiviruses were used to manipulate the expression of LINC00941, miR-335-5p, and ROCK1 in PC cell lines. Gain or loss-of-function assays and mechanistic assays were employed to verify the roles of LINC00941, miR-335-5p, and ROCK1 in PC cell growth and metastasis, both in vivo and in vitro. LINC00941 and ROCK1 were found to be highly expressed in PC, while miR-335-5p exhibited low expression. High LINC00941 expression was strongly associated with larger tumor size, lymph node metastasis, and poor prognosis. Functional experiments revealed that LINC00941 silencing significantly suppressed PC cell growth, metastasis and epithelial–mesenchymal transition. LINC00941 functioned as a molecular sponge for miR-335-5p, and a competitive endogenous RNA (ceRNA) for ROCK1, promoting ROCK1 upregulation, and LIMK1/Cofilin-1 pathway activation. Our observations lead us to conclude that LINC00941 functions as an oncogene in PC progression, behaving as a ceRNA for miR-335-5p binding. LINC00941 may therefore have potential utility as a diagnostic and treatment target in this disease.Subject terms: Pancreatic cancer, Long non-coding RNAs     

Dysregulated ceRNA network modulated by copy number variation-driven lncRNAs in breast cancer: A comprehensive analysis

Breast cancer is a malignancy harmful to physical and mental health in women, with quite high mortality. Copy number variations (CNVs) are vital factors affecting the progression of breast cancer. Detecting CNVs in breast cancer to predict the prognosis of patients has become a promising approach to accurate treatment in recent years. The differential analysis was performed on CNVs of long noncoding RNAs (lncRNAs) as well as the expression of lncRNAs, microRNAs (miRNAs) and mRNAs in normal tissue and breast tumor tissue based on The Cancer Genome Atlas (TCGA) database. The CNV-driven lncRNAs were identified by the Kruskal–Wallis test. Meanwhile, a competitive endogenous RNA (ceRNA) network regulated by CNV-driven lncRNA was constructed. As the Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses revealed, the mRNAs in the dysregulated ceRNA network were mainly enriched in the biological functions and signaling pathways, including the Focal Adhesion-PI3K-Akt–mTOR-signaling pathway, the neuronal system, metapathway biotransformation Phase I and II and blood circulation, etc. The relationship between the CNVs of five lncRNAs and their gene expression in the ceRNA network was analyzed via a chi-square test, which confirmed that except for LINC00243, the expression of four lncRNAs was notably correlated with the CNVs. The survival analysis revealed that only the copy number gain of LINC00536 was evidently related to the poor prognosis of patients. The CIBERSORT algorithm showed that five lncRNAs were correlated with the abundance of immune cell infiltration and immune checkpoints. In a word, by analyzing CNV-driven lncRNAs and the ceRNA network regulated by these lncRNAs, this study explored the mechanism of breast cancer and provided novel insights into new biomarkers.     

Circular RNA-associated ceRNA network involved in HIF-1 signalling in triple-negative breast cancer: circ_0047303 as a potential key regulator

The aggressive and highly metastatic nature of triple-negative breast cancer (TNBC) causes patients to suffer from the poor outcome. HIF-1 signalling pathway is a prominent pathway that contributes to angiogenesis and metastasis progression in tumours. On the contrary, the undeniable importance of circular RNAs (circRNAs) as multifunctional non-coding RNAs (ncRNAs) has been identified in breast cancer. These ncRNAs owing to their high stability and specificity have been becoming a hotspot in cancer researches. circRNAs act as competing endogenous RNAs (ceRNAs) and compete with mRNAs for shared miRNAs, thus modulate gene expression. Since the most dysregulated biological functions in TNBC are associated with cellular invasion, understanding the molecular pathogenesis of these processes is a crucial step towards the development of new treatment approaches. The purpose of this study is to undermine the circRNA-associated ceRNA network involved in HIF-1 signalling in TNBC using an integrative bioinformatics approach. In the next step, the novel circ_0047303-mediated ceRNA regulatory axes have been extracted and validated across TNBC samples. We show that circ_0047303 has the highest degree in the circRNA-associated ceRNA network and shows a significant up-expression in TNBC. Moreover, our results suggest that circ_0047303 could mediate the upregulation of key angiogenesis-related genes, including HIF-1, EIF4E2 and VEGFA in TNBC through sponging the tumour-suppressive miRNAs. The circ_0047303 could be a promising molecular biomarker and/or therapeutic target for TNBC.     

Screening of significant biomarkers related with prognosis of liver cancer by lncRNA-associated ceRNAs analysis

This study aimed to identify significant biomarkers related to the prognosis of liver cancer using long noncoding RNA (lncRNA)-associated competing endogenous RNAs (ceRNAs) analysis. Differentially expressed mRNA and lncRNAs between liver cancer and paracancerous tissues were screened, and the functions of these mRNAs were predicted by gene ontology and pathway enrichment analyses. A ceRNA network consisting of differentially expressed mRNAs and lncRNAs was constructed. LncRNA FENDRR and lncRNA HAND2-AS1 were hub nodes in the ceRNA network. A risk score assessment model consisting of eight genes (PDE2A, ESR1, FBLN5, ALDH8A1, AKR1D1, EHHADH, ADRA1A, and GNE) associated with prognosis were developed. Multivariate Cox regression suggested that both pathologic_T and risk group could be regarded as independent prognostic factors. Furthermore, a nomogram model consisting of pathologic_T and risk group showed a good prediction ability for predicting the survival rate of liver cancer patients. The nomogram model consisting of pathologic_T and a risk score assessment model could be regarded as an independent factor for predicting prognosis of liver cancer.     

Comprehensive analysis of lncRNA-associated ceRNA network in colorectal cancer

Colorectal cancer (CRC) is one of the most common malignancies and morbidity and mortality are increasing rapidly. Increasing evidence showed the close correlation between aberrant expression of certain RNAs and the occurrence and development of CRC. However, comprehensive analyses of differentially expressed profiles of linRNA in CRC based on large sample size have been lacking. In the present study, based on RNA-seq data obtained from the TCGA (The Cancer Genome Atlas) database, we identified 1176 lncRNAs, 245 miRNAs and 2083 mRNAs whichaberrantly expressed in the colorectal cancer tissues compared with the adjacent non-tumorous tissues. A Kaplan-Meier curve analysis was used to study the overall survival rate of the three RNA-related CRC patients. After constructing the ceRNA network, we performed the KEGG enrichment pathway analysis on ceRNA-related differentially expressed mRNAs and found that these mRNAs were remarkably enriched in the pathways associated with CRC. Combining the differentially expressed lncRNAs with clinical pathological variables of CRC patients, we also found that LINC00400 and LINC00355 not only contribute to the regulation of ceRNA network, but also show significantchanges in its expression in multiple CRC pathological stages, indicating that LINC00400 and LINC00355 can be considered as promising therapeutic targets for CRC.     

Integrated analysis of a competing endogenous RNA network reveals key long noncoding RNAs as potential prognostic biomarkers for hepatocellular carcinoma

Growing evidence has revealed that long noncoding RNAs (lncRNAs) have an important impact on tumorigenesis and tumor progression via a mechanism involving competing endogenous RNAs (ceRNAs). However, their use in predicting the survival of a patient with hepatocellular carcinoma (HCC) remains unclear. The aim of this study was to develop a novel lncRNA expression–based risk score system to accurately predict the survival of patients with HCC. In our study, using expression profiles downloaded from The Cancer Genome Atlas database, the differentially expressed messenger RNAs (mRNAs), lncRNAs, and microRNAs (miRNAs) were explored in patients with HCC and normal liver tissues, and then a ceRNA network constructed. A risk score system was established between lncRNA expression of the ceRNA network and overall survival (OS) or recurrence-free survival (RFS); it was further analyzed for associations with the clinical features of patients with HCC. In HCC, 473 differentially expressed lncRNAs, 63 differentially expressed miRNAs, and 1417 differentially expressed mRNAs were detected. The ceRNA network comprised 41 lncRNA nodes, 12 miRNA nodes, 24 mRNA nodes, and 172 edges. The lncRNA expression–based risk score system for OS was constructed based on six lncRNAs (MYLK-AS1, AL359878.1, PART1, TSPEAR-AS1, C10orf91, and LINC00501), while the risk score system for RFS was based on four lncRNAs (WARS2-IT1, AL359878.1, AL357060.1, and PART1). Univariate and multivariate Cox analyses showed the risk score systems for OS or RFS were significant independent factors adjusted for clinical factors. Receiver operating characteristic curve analysis showed the area under the curve for the risk score system was 0.704 for OS, and 0.71 for RFS. Our result revealed a lncRNA expression–based risk score system for OS or RFS can effectively predict the survival of patients with HCC and aid in good clinical decision-making.