基于TCGA和GEO数据库的胃癌lncRNA筛选与ceRNA网络构建

胃癌（gastric cancer,GC）是我国最常见的恶性肿瘤之一,严重危害人类健康。胃癌发病机制复杂,缺乏特异性预后生物标志物。长链非编码RNA（long non?coding RNA,lncRNA）可作为竞争性内源RNA（competing endogenous RNA,ceRNA）,影响microRNA（miRNA）与mRNA的结合,从而影响胃癌的发生、发展。基于TCGA和GEO数据库的转录组数据,筛选GC中差异表达的lncRNAs,并构建基于6条lncRNAs（HAGLROS、TMEM92?AS1、LINC01745、HOXC?AS3、SEMA3B?AS1、FEZF1?AS1）的lncRNA?miRNA?mRNA网络。网络核心基因的KEGG/GO富集和蛋白质互作分析结果显示,lncRNA可能通过miRNA海绵吸附作用,调控胃癌的发生、发展与转移。AC011352.1、AC087636.1、AC093627.1、GAS1RR与胃癌患者的预后相关性具有统计学意义（P<0.05）,并可能成为胃癌患者潜在的预后生物标志物。     

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中高表达、低表达及治疗相关方面的作用.     

LncRNA作为ceRNA调控肿瘤的发生发展

长链非编码RNA(lncRNA)是长度大于200 bp,不编码蛋白质的内源性RNA分子.近年来的研究表明,lncRNA可以作为一种竞争性内源RNA(competing endogenous RNA,ceRNA)吸附miRNA,参与靶基因的表达调控,从而在肿瘤的发生发展中发挥重要的作用.本文从lncRNA作为ceRNA发挥生物学功能这一角度,概述了相关lncRNA在肿瘤发生发展中的作用及机制.揭秘lncRNA与miRNA在肿瘤发生中的相互作用,将为肿瘤的诊断和治疗提供新思路.     

非编码RNA作为ceRNA在人癌症中的功能及机制

非编码RNA(non-coding RNA, ncRNA),即无编码蛋白质潜力的RNA,包括短非编码RNA,例如微RNA(microRNA, miRNA),长非编码RNA(long coding RNA, lncRNA)和环状RNA(circular RNA, circRNA)等,已被证明可以在RNA水平上调节细胞中多种生理及病理过程。miRNA是约18～22个核苷酸大小的内源性非编码RNA分子,可以通过MRE与靶基因的3′非翻译区(untranslated region,UTR)中的互补序列结合,抑制蛋白质编码基因的表达,并导致mRNA转录物的更新、转换或降解。2011年,Salmena等提出,非编码RNA与具有编码蛋白质能力的RNA(mRNA)之间存在一种被称为竞争性内源RNA(competing endogenous RNA, ceRNA)的相互作用机制假说,即含有miRNA反应元件(MRE)的ncRNA通过与miRNA结合,解除miRNA对靶基因的抑制作用。这一假说对基因表达调控的传统认知进行了补充,在RNA水平上将多种ncRNA的作用补充到经典的miRNA调控mRNA翻译的过程,将其延伸为ceRNA-miRNA-mRNA的网络调控模式。近年来研究发现,ceRNA机制广泛存在于胃癌、结肠癌和膀胱癌等各类癌症中,并且在肿瘤的基因调控及肿瘤细胞的增殖、侵袭、转移、凋亡、细胞周期等生物过程中发挥作用。本文将介绍ceRNA及其网络的机制与分子基础,并且结合两类非编码RNA——lncRNA及circRNA作为ceRNA分别在人类不同癌症类型中的近期研究进展作一综述,讨论该机制在肿瘤中的角色及作用,以期拓宽对肿瘤发生发展机制的视野,为癌症治疗提供新思路。     

长非编码RNA作为ceRNA分子在乳腺癌中的调控机制

乳腺癌是全球女性最常见的异质性恶性肿瘤，分为多种亚型，不同亚型的临床治疗方式和效果不同。多项研究表明长非编码RNA(lncRNA)在乳腺癌的进展中起着关键作用，其中由lncRNA-miRNA-mRNA相互作用形成的ceRNA网络是RNA分子间的一种新调节机制。LncRNA通过该网络广泛参与乳腺癌的增殖、迁移、凋亡等多种生物学过程，可作为诊断和治疗的分子靶点。该文重点阐述了lncRNA介导的ceRNA调控机制，系统总结了lncRNA相关ceRNA网络在调控不同乳腺癌亚型的增殖、迁移和耐药性中的研究进展，旨在为乳腺癌的精准治疗和lncRNA研究提供新见解。     

非编码RNA与基因表达调控

近年来,随着对基因组的深入研究,发现真核生物中存在许多形态和功能各异的非编码RNA分子,这类RNA分子并不表达蛋白质,但它们在基因转录水平、转录后水平及翻译水平起了重要的调控作用。具有调控作用的RNA分子种类非常丰富,如长链非编码RNA(long non-coding RNA,lncRNA)、miRNA、PIWI相互作用RNA(PIWI-interacting RNA,piRNA)、内源性小干扰RNA(endogenous small interfering RNA,endo-siRNA)、竞争性内源RNA(competitive endogenous RNA,ceRNA)等,它们使基因表达过程更为丰富、严谨和有序。本文综述几类典型的非编码RNA对基因表达的调节作用,以助于理解细胞中RNA分子调节网络的功能和机制。     

中华蜜蜂6日龄幼虫响应蜜蜂球囊菌侵染的长链非编码RNA应答研究

[目的] 本研究旨在探究长链非编码RNA（long non-coding RNA,lncRNA）在中华蜜蜂（Apis cerana cerana,简称中蜂）6日龄幼虫应答蜜蜂球囊菌（Ascosphaera apis,简称球囊菌）侵染过程中的差异表达谱及调控作用。[方法] 利用链特异性cDNA建库的RNA-seq技术对未被侵染及球囊菌侵染的中蜂6日龄幼虫肠道（AcCK和AcT）进行深度测序。通过相关生物信息学软件分析lncRNA的结构特征和表达谱。筛选并分析差异表达lncRNA（differentially expressed lncRNA,DElncRNA）的顺式（cis）作用及竞争性内源RNA（competing endogenous RNA,ceRNA）调控网络。采用RT-qPCR验证测序数据及DElncRNA差异变化趋势的可靠性。[结果] AcCK和AcT共预测出642个已知lncRNA和487个新lncRNA。与蛋白编码基因相比,上述中蜂lncRNA外显子数更少、长度更短且表达量更低。43个antisense lncRNA与40个正义链mRNA之间互补配对。AcCK和AcT比较组包含367个上调lncRNA和268个下调lncRNA。有194个DElncRNA潜在调控461个上下游基因,并涉及细胞进程、代谢进程和催化活性等38个功能条目以及氨基酸代谢、内吞作用和MAPK等191条通路。此外,有180个DElncRNA可靶向结合50个DEmiRNA,进而调控6365个mRNA;三者之间形成较为复杂的ceRNA调控网络。[结论] 中蜂6日龄幼虫肠道的部分lncRNA可作为antisense lncRNA参与应答球囊菌侵染;部分DElncRNA可通过cis作用调节物质代谢和免疫途径相关的上下游基因,从而介导宿主的侵染应答;TCONS_00010661和TCONS_00003104等DElncRNA可通过ceRNA网络调控Jak-STAT和氧化磷酸化等通路及富集基因,进而参与宿主的侵染应答。     

lncRNA与miRNA相互调控作用及其与肿瘤的关系研究

长链非编码RNA (long non-coding RNA,lncRNA)是一类转录本长度大于200 bp的非编码RNA,可作为人类基因组中一类重要的调控分子通过多种方式发挥其生物学功能.近年来的研究表明,lncRNA也可以作为一种竞争性内源性RNA (competing endogenous RNA, ceRNA) 与miRNA相互作用,参与靶基因的表达调控,并在肿瘤的发生发展中发挥重要的作用.本综述在简要介绍lncRNA功能研究现状和主要研究方法的基础上,进一步分析了lncRNA与miRNA之间的互相调控关系及其在肿瘤发生发展中的作用,以便为后续的研究提供新的思路.     

竞争性内源RNA的生物学功能及其调控

最新研究表明,RNA之间可以通过竞争结合共同的microRNA反应元件(microRNA response element, MRE)实现相互调节,这种调控模式构成竞争性内源RNA(Competing endogenous RNA, ceRNA)。已发现的ceRNA包括蛋白编码mRNA和非编码RNA,其中后者包括假基因转录物、长链非编码RNA(Long non-coding RNA, lncRNA)、环状RNA(Circular RNA, circRNA)等。文章主要从ceRNA分类的角度,阐述各类ceRNA构成的调控网络发挥的生物学功能在病理和生理相关过程中的作用,以及可能影响ceRNA调控有效性的因素。     

利用TCGA数据库构建肾透明细胞癌相关miRNA预后模型

利用TCGA数据库中肾透明细胞癌的miRNA与mRNA数据及临床信息,构建由miRNA组成的预后风险评分模型,并筛选与生存预后相关的miRNA-mRNA调控关系对,为研究提供理论依据。下载并整理TCGA[JP+1]数据库中肾透明细胞癌的miRNA与mRNA数据;对数据进行差异分析,将差异表达的miRNA与临床信息进行合并,利用单因素与多因素Cox回归分析,构建预后模型并进行模型评价;对模型中的miRNA进行靶基因预测,结果与差异表达的mRNA进行取交集,构建miRNA-mRNA调控网络;对网络中的mRNA进行生存分析,筛选生存相关的miRNA-mRNA调控关系对。共得到49个差异表达的miRNA与3 613个差异表达的mRNA;预后模型计算公式为:风险值(risk score)=hsa-miR-21-5p表达量×0.603＋hsa-miR-1251-5p表达量×-0.093;调控网络中共纳入31个miRNA-mRNA调控关系对;对mRNA进行生存分析,共得到7个有价值的关系对。所构建预后模型可有效预测肾透明细胞癌患者生存预后情况,筛选到的miRNA-mRNA调控关系对可为相关研究与治疗提供参考。     

Identification of a potential prognostic lncRNA-miRNA-mRNA signature in endometrial cancer based on the competing endogenous RNA network

Endometrial cancer is one of the most common gynecological malignant tumors. The roles of competing endogenous RNAs (ceRNAs) in this disease, however, remain unclear. In this study, we constructed a ceRNA network to reveal the core ceRNAs in endometrial cancer. Differentially expressed genes were summarized from The Cancer Genome Atlas database, whereupon 140 genes were identified for building the network. Further correlation, survival, and enrichment analyses suggested that these genes may help towards elucidating the molecular mechanisms of endometrial cancer. After validation of the findings with the GSE17025 data set, LINC00958, microRNA-761, and DOLPP1 were highlighted as the critical genes in the ceRNA network. Our work suggests that LINC00958 may regulate DOLPP1 by “sponging” miR-761 in endometrial 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.     

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.     

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.     

Comprehensive analysis to identify the RP11–478C19.2/ E2F7 axis as a novel biomarker for treatment decisions in clear cell renal cell carcinoma

Clear cell renal cell carcinoma (ccRCC), accounting for 70–80% of all renal cell carcinomas, is a common malignancy. Survival rates decrease significantly in patients with advanced and metastatic ccRCC. Furthermore, ccRCC is less responsive to radiation and chemotherapy than other cancers. Therefore, targeted therapy and immunotherapy are particularly important for ccRCC management. A growing body of literature recognizes that competitive endogenous RNA (ceRNA) regulatory networks play a crucial role in various cancers. However, the biological functions of the ceRNA network in ccRCC require further investigation. In this study, we built the ceRNA network for ccRCC using the “GDCRNATools” package. After survival analysis, the RP11–478C19.2/hsa-miR-181b-5p, hsa-miR-181a-5p, and hsa-miR-181c-5p/E2F7 axes were obtained for further analysis. Unsupervised clustering was conducted basing on this ceRNA network. The results indicated that the prognosis and immune infiltration levels differed between the two clusters. Furthermore, we conducted correlation analysis, immune infiltration analysis, tumor mutation burden analysis, GSEA analysis, drug sensitivity analysis and pan-cancer analysis of E2F7 to explore its potential role in oncogenesis. Experiments in vitro were performed to confirm the pro-oncogenic impact of E2F7. The results suggest that the RP11–478C19.2/E2F7 axis might be a biomarker for the inclusion of cabozantinib, pazopanib, sunitinib, and immunotherapy in the therapeutic regimen. In summary, we found that the ceRNA-based RP11–478C19.2/E2F7 axis is involved in ccRCC and that it could be a novel biomarker for treatment decisions and a possible therapeutic target to increase the success of targeted therapy and immunotherapy in ccRCC.     

Effects of long non-coding RNA LINC00667 on renal tubular epithelial cell proliferation,apoptosis and renal fibrosis via the miR-19b-3p/LINC00667/CTGF signaling pathway in chronic renal failure

The global prevalence of chronic renal failure (CRF) has significantly elevated with various reports indicating there to be a 10% worldwide rate. The functions of long non-coding RNAs (lncRNAs) and their deeper association with CRF at present remain poorly understood. Hence, the aim of the present study was to investigate the altered expressions of lncRNA LINC00667 in CRF and its associated effects on renal tubular epithelial cell proliferation, apoptosis and renal fibrosis through the microRNA-19b-3p (miR-19b-3p)/LINC00667/connective tissue growth factor (CTGF) signaling pathway. Initially, verification of the targeting relationship between LINC00667, CTGF and miR-19b-3p was performed, after which evidence was obtained indicating that miR-19b-3p could negatively regulate LINC00667 and CTGF. The expressions of CTGF in both the CRF and normal renal tissues were determined by immunohistochemistry means, with LINC00667 and CTGF determined to be highly expressed, while poor expression levels of miR-19b-3p were detected among the CRF tissues. The expressions of LINC00667, miR-19b-3p, fibrosis- and epithelial-mesenchymal transition (EMT)-related genes were also examined. The successfully established CRF rat models were treated with varying mimics, inhibitors, and siRNA. ELISA was applied to determine the renal function-related factors. Besides, the renal cell proliferation, migration and apoptosis were detected. In response to LINC00667 silencing, the renal tubular epithelial cells displayed increased proliferation and migration accompanied by reduced apoptosis based on upregulated miR-19b-3p, along with inhibited renal fibrosis and EMT detected. Taken together, the key findings of our study demonstrated that decreased lncRNA LINC00667 could promote renal tubular epithelial cell proliferation and ameliorate renal fibrosis in CRF via the miR-19b-3p/LINC00667/CTGF signaling pathway.     

Integrated analysis of ceRNA network and tumor-infiltrating immune cells in esophageal cancer

Background: Esophageal cancer (ESCA) is one of the most commonly diagnosed cancers in the world. Tumor immune microenvironment is closely related to tumor prognosis. The present study aimed at analyzing the competing endogenous RNA (ceRNA) network and tumor-infiltrating immune cells in ESCA.Methods: The expression profiles of mRNAs, lncRNAs, and miRNAs were downloaded from the Cancer Genome Atlas database. A ceRNA network was established based on the differentially expressed RNAs by Cytoscape. CIBERSORT was applied to estimate the proportion of immune cells in ESCA. Prognosis-associated genes and immune cells were applied to establish prognostic models basing on Lasso and multivariate Cox analyses. The survival curves were constructed with Kaplan–Meier method. The predictive efficacy of the prognostic models was evaluated by the receiver operating characteristic (ROC) curves.Results: The differentially expressed mRNAs, lncRNAs, and miRNAs were identified. We constructed the ceRNA network including 23 lncRNAs, 19 miRNAs, and 147 mRNAs. Five key molecules (HMGB3, HOXC8, HSPA1B, KLHL15, and RUNX3) were identified from the ceRNA network and five significant immune cells (plasma cells, T cells follicular helper, monocytes, dendritic cells activated, and neutrophils) were selected via CIBERSORT. The ROC curves based on key genes and significant immune cells all showed good sensitivity (AUC of 3-year survival: 0.739, AUC of 5-year survival: 0.899, AUC of 3-year survival: 0.824, AUC of 5-year survival: 0.876). There was certain correlation between five immune cells and five key molecules.Conclusion: The present study provides an effective bioinformatics basis for exploring the potential biomarkers of ESCA and predicting its prognosis.     

Identification of key pathways and RNAs associated with skeletal muscle atrophy after spinal cord injury

Objective:This study was performed to investigate the potential key molecules involved in the progression of skeletal muscle atrophy after SCI.Methods:Based on {"type":"entrez-geo","attrs":{"text":"GSE21497","term_id":"21497"}}GSE21497 dataset, the DEmRNAs and DElncRNAs were screened after differentially expressed analysis. Then the enrichment analyses were performed on DEmRNAs. Then the PPI network and ceRNA network were constructed. Finally, the DGIdb was utilized to predict drug-gene interactions.Results:A total of 412 DEmRNAs and 21 DElncRNAs were obtained. The DEmRNAs were significantly enriched in MAPK signaling pathway and FoxO signaling pathway. In addition, UBE2D1, JUN, and FBXO32 had higher node degrees in PPI network, and the top 20 genes with high degree were significantly enriched in FoxO signaling pathway and Endometrial cancer. Moreover, FOXO3 was regulated by hsa-miR-1207-5p and hsa-miR-1207-5p was regulated by lncRNA RP11-253E3.3 in ceRNA network. Finally, 37 drug-gene interactions were obtained based on the 26 genes in ceRNA network.Conclusion:UBE2D1, JUN, and FBXO32 are likely to be related to the progression of skeletal muscle atrophy after SCI, and activating of MAPK signaling pathway, Endometrial cancer and FoxO signaling pathway may induce skeletal muscle inflammation, apoptosis, autophagy and atrophy after SCI. Moreover, RP11-253E3.3-hsa-miR-1207-5p-FOXO3 axis may be a promising therapeutic target for skeletal muscle atrophy after SCI.