基于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-AS 1)的lncRNA-miRNA-mRNA网络。网络核心基因的KEGG/GO富集和蛋白质互作分析结果显示,lncRNA可能通过miRNA海绵吸附作用,调控胃癌的发生、发展与转移。AC011352.1、AC087636.1、AC093627.1、GAS1RR与胃癌患者的预后相关性具有统计学意义(P<0.05),并可能成为胃癌患者潜在的预后生物标志物。     

Identification and development of long non‐coding RNA‐associated regulatory network in colorectal cancer

Colorectal cancer (CRC) is one of the leading causes of cancer‐associated death globally. Long non‐coding RNAs (lncRNAs) have been identified as micro RNA (miRNA) sponges in a competing endogenous RNA (ceRNA) network and are involved in the regulation of mRNA expression. This study aims to construct a lncRNA‐associated ceRNA network and investigate the prognostic biomarkers in CRC. A total of 38 differentially expressed (DE) lncRNAs, 23 DEmiRNAs and 27 DEmRNAs were identified by analysing the expression profiles of CRC obtained from The Cancer Genome Atlas (TCGA). These RNAs were chosen to develop a ceRNA regulatory network of CRC, which comprised 125 edges. Survival analysis showed that four lncRNAs, six miRNAs and five mRNAs were significantly associated with overall survival. A potential regulatory axis of ADAMTS9‐AS2/miR‐32/PHLPP2 was identified from the network. Experimental validation was performed using clinical samples by quantitative real‐time PCR (qRT‐PCR), which showed that expression of the genes in the axis was associated with clinicopathological features and the correlation among them perfectly conformed to the ‘ceRNA theory’. Overexpression of ADAMTS9‐AS2 in colon cancer cell lines significantly inhibited the miR‐32 expression and promoted PHLPP2 expression, while ADAMTS9‐AS2 knockdown had the opposite effects. The constructed novel ceRNA network may provide a comprehensive understanding of the mechanisms of CRC carcinogenesis. The ADAMTS9‐AS2/miR‐32/PHLPP2 regulatory axis may serve as a potential therapeutic target for CRC.     

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.     

Functional role of a long non-coding RNA LIFR-AS1/miR-29a/TNFAIP3 axis in colorectal cancer resistance to pohotodynamic therapy

Colorectal Cancer (CRC) is one of the most common digestive system malignant tumors. Recently, PDT has been used as a first-line treatment for colon cancer; however, limited curative effect was obtained due to resistance of CRC to PDT. During the past decades, accumulating CRC-related long non-coding RNAs (lncRNAs), microRNAs (miRNAs) and mRNAs have been reported to exert diverse functions through various biological processes; their dysregulation might trigger and/or promote the pathological changes. Herein, we performed microarrays analysis to identify dysregulated lncRNAs, miRNAs and mRNAs in PDT-treated HCT116 cells to figure out the lncRNA-miRNA interactions related to the resistance of CRC to PDT treatment, and the downstream mRNA target, as well as the molecular mechanism. We found a total of 1096 lncRNAs dysregulated in PDT-treated CRC HCT116 cells; among them, LIFR-AS1 negatively interacted with miR-29a, one of the dysregulated miRNAs in PDT-treated CRC cells, to affect the resistance of CRC to PDT. LIFR-AS1 knockdown attenuated, whereas miR-29a inhibition enhanced the cellular effect of PDT on HCT116 cell proliferation and apoptosis. Furthermore, among the dysregulated mRNAs, TNFAIP3 was confirmed to be a direct target of miR-29a and exerted a similar effect to LIFR-AS1 on the cellular effects of PDT. In summary, LIFR-AS1 serves as a competitive endogenous RNA (ceRNA) for miR-29a to inhibit its expression and up-regulate downstream target TNFAIP3 expression, finally modulating the resistance of CRC to PDT. We provide an experimental basis for this lncRNA/miRNA/mRNA network being a promising target in CRC resistance to PDT treatment.     

Genome-wide circular RNA (circRNA) and mRNA profiling identify a circMET-miR-410-3p regulatory motif for cell growth in colorectal cancer

Circular RNA (circRNA) is a non-coding RNA molecule that lacks polyadenylated tails and is highly stable, abundant, and conserved in human cells. CircRNAs can serve as a competing endogenous RNA (ceRNA) to sponge microRNAs (miRNA) and block their effects on target mRNA expression. CircRNAs also have possible relevance to cancer and therefore may be considered as ideal biomarkers for monitoring cancer progression. Of the about 300,000 predicted human circRNAs, only a few have validated biological functions related to cancer. To better understand the ceRNA role of circRNAs in colorectal cancer (CRC), we performed genome-wide circRNA-based RNA-sequencing (RNA-Seq) on nine CRC tumor samples and their paired histologically normal adjacent tissue samples. By profiling the mRNA expression in the same patients, we further explored the expression correlation between circRNAs and mRNAs generated from the same parental gene. Focusing on the concordant differential expression between circRNAs and mRNAs, we substantially reduced the regulatory noise. In total, we identified 694 circRNA-mRNA pairs that were consistently up or downregulated between tumor and normal tissues. These 694 circRNA-mRNA pairs are from 182 protein-coding genes associated with hormone responses and chemotaxis. Of these 182 genes, 43 are downstream targets of three highly conserved miRNAs (miR-410-3p, miR-135a, and miR-30a). Interestingly, these 43 genes are highly mutated in another cohort from eight independent CRC studies, which have significant effects on patient survival time. Focusing on miR-410-3p and its target oncogene MET, we experimentally validated the ceRNA regulatory motif of circMET. Notably, circMET is substantially upregulated in CRC cell lines and could promote cell proliferation and growth. By confirming the regulatory relationship between miR-410-3p and circMET, we propose a new mechanism for the observed sustained activation of MET in CRC. In conclusion, our work identifies a novel regulatory role of circMET and provides a potential diagnostic biomarker for CRC.     

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.     

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

子宫内膜癌(uterine corpus endometrial carcinoma,UCEC)是危害女性健康的癌症之一,但其发生发展机制尚不完全清楚.基因的异常表达在细胞癌变过程中发挥着重要的作用.本研究利用生物信息学方法对UCEC中异常表达的基因进行网络调控分析,为UCEC的机制研究及预后治疗提供理论依据.首先,利...     

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

Competing endogenous RNA network and prognostic nomograms for hepatocellular carcinoma patients who underwent R0 resection

The prognosis of hepatocellular carcinoma (HCC) after R0 resection is unsatisfactory due to the high rate of recurrence. In this study, we investigated the recurrence-related RNAs and the underlying mechanism. The long noncoding RNA (lncRNA), microRNA (miRNA), and messenger RNA (mRNA) expression data and clinical information of 247 patients who underwent R0 resection patients with HCC were obtained from The Cancer Genome Atlas. Comparing the 1-year recurrence group (n = 56) with the nonrecurrence group (n = 60), we detected 34 differentially expressed lncRNAs (DElncRNAs), five DEmiRNAs, and 216 DEmRNAs. Of these, three DElncRNAs, hsa-mir-150-5p, and 11 DEmRNAs were selected for constructing the competing endogenous RNA (ceRNA) network. Next, two nomogram models were constructed based separately on the lncRNAs and mRNAs that were further selected by Cox and least absolute shrinkage and selection operator regression analysis. The two nomogram models that showed a high prediction accuracy for disease-free survival with the concordance indexes at 0.725 and 0.639. Further functional enrichment analysis of DEmRNAs showed that the mRNAs in the ceRNA network and nomogram models were associated with immune pathways. Hence, we constructed a hsa-mir-150-5p-centric ceRNA network and two effective nomogram prognostic models, and the related RNAs may be useful as potential biomarkers for predicting recurrence in patients with HCC.     

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.     

microRNA-214在恶性肿瘤中的表达及相关性的研究进展

微小RNA(microRNA,miRNA)是广泛存在于动植物中的一类不编码蛋白质的短小的单链RNA分子,一般由22个核苷酸组成,它们可以特异性地结合mRNA并通过降解或抑制其翻译而在转录后水平调控基因表达。miRNA的表达及功能可影响许多表观遗传学特征,其功能涉及细胞的发生、生长、发育、分化和凋亡过程,在肿瘤的形成和进展过程中扮演重要角色。microRNA-214(miRNA-214,miR-214)参与肝癌、乳腺癌、宫颈癌、卵巢癌、恶性黑色素瘤、胃癌、胶质瘤、儿童骨肉瘤等恶性肿瘤的发生发展,以及与肿瘤细胞的侵袭及转移密切相关。miRNA-214在不同的肿瘤中表达水平并不相同,miRNA-214在不同肿瘤中的差异表达是通过调控某个或者某些癌基因及抑癌基因而实现其参与肿瘤的发生发展、侵袭及转移的作用。因此,本文主要通过阅读大量国内外文献,总结和概括了miRNA-214参与部分恶性肿瘤发生发展的机制。虽然目前对于miRNA的理论研究已经日渐完善和成熟,但是怎样将这些研究结果应用于临床,怎样能够更准确、更便捷的通过对miRNA的检测达到对疾病的诊断、治疗以及预后评估,想必一定会成为将来研究的热点,我们期待一种新型的恶性肿瘤的分子标志物会使越来越多的肿瘤患者获益。     

The exploration of new therapeutic targets for HPV-negative head and neck squamous cell cancer through the construction of a ceRNA network and immune microenvironment analysis

Previous studies have shown that human papillomavirus (HPV)-negative patients with head and neck squamous cell cancer (HNSCC) suffer from an unsatisfactory prognosis. Long noncoding RNAs (lncRNAs) have been verified to participate in many biological processes, including regulating gene expression as competing endogenous RNAs (ceRNAs), while few studies focused the ceRNA network regulation mechanism in patients with HPV-negative HNSCC tumor. Meanwhile, the immune microenvironment may be critical in the development and prognosis of HPV-negative tumors. Our study aimed to further investigate the pathogenesis and potential biomarkers for the diagnosis, therapy and prognosis of HPV-negative HNSCC through a ceRNA network. Comprehensively analyzing the sequencing data of lncRNAs, microRNAs (miRNAs), and messenger RNAs (mRNAs) in The Cancer Genome Atlas HNSCC dataset, we constructed a differentially expressed ceRNA network containing 131 lncRNAs, 35 miRNAs and 162 mRNAs. Then, survival analysis in the network was cited to explore the prognostic biomarkers. Eight mRNAs, nine lncRNAs, and one miRNA were identified to be associated with prognosis. Neuropilin (NRP) binding function, retinoid X receptor (RXR) binding, and the vascular endothelial growth factor (VEGF) signaling pathway were associated with the enrichment analysis, and they also related to the immune microenvironment. Combined with the analysis of the immune microenvironment differences, we obtained new targeted therapies using an RXR agonist, or a combination of the VEGF monoclonal antibody and an NRP antagonist, which may provide a promising future for HPV-negative HNSCC patients.     

H3K27 acetylation-induced lncRNA EIF3J-AS1 improved proliferation and impeded apoptosis of colorectal cancer through miR-3163/YAP1 axis

Long noncoding RNAs (lncRNAs) are found to be aberrantly expressed and pose significant impacts in colorectal cancer (CRC), the most prevalent type malignancy in the gastrointestinal tract. This study aimed to find out the regulation of lncRNA EIF3J antisense RNA 1 (EIF3J-AS1) on CRC progression. Expressions of EIF3J-AS1, microRNA-3163 (miR-3163), and Yes-associated protein 1 (YAP1) in tissues and cells were evaluated by real-time quantitative polymerase chain reaction (RT-qPCR) and Western blot analysis. Association of EIF3J-AS1 with CRC prognosis was analyzed through the online bioinformatics tool GEPIA. The biological function of EIF3J-AS1 in CRC was investigated by Cell Counting Kit-8, colony formation, caspase-3 activity, and TUNEL staining. Competitive endogenous RNA (ceRNA) network of EIF3J-AS1/miR-3163/YAP1 was determined by luciferase reporter and RNA immunoprecipitation assays. Results showed that EIF3J-AS1 was upregulated in CRC tissues and cell lines, indicating poor prognosis of CRC patients. The silence of EIF3J-AS1 led to reduced proliferation and facilitated apoptosis of CRC cells. Mechanistcally, EIF3J-AS1 was upregulated by cAMP-response element-binding protein-binding protein-mediated histone H3 on lysine 27 acetylation (H3K27ac) at the promoter region, and EIF3J-AS1 upregulated YAP1 expression through sponging miR-3163 in CRC cells. In conclusion, we first found that H3K27 acetylation-induced lncRNA EIF3J-AS1 improved proliferation and impeded apoptosis of colorectal cancer through the miR-3163/YAP1 axis, which might potentially provide a novel molecular-targeted strategy for CRC treatment.     

Prediction of a competing endogenous RNA co‐expression network as a prognostic marker in glioblastoma

Due to its high proliferation capacity and rapid intracranial spread, glioblastoma (GBM) has become one of the least curable malignant cancers. Recently, the competing endogenous RNAs (ceRNAs) hypothesis has become a focus in the researches of molecular biological mechanisms of cancer occurrence and progression. However, there is a lack of correlation studies on GBM, as well as a lack of comprehensive analyses of GBM molecular mechanisms based on high‐throughput sequencing and large‐scale sample sizes. We obtained RNA‐seq data from The Cancer Genome Atlas (TCGA) and Genotype‐Tissue Expression (GTEx) databases. Further, differentially expressed mRNAs were identified from normal brain tissue and GBM tissue. The similarities between the mRNA modules with clinical traits were subjected to weighted correlation network analysis (WGCNA). With the mRNAs from clinical‐related modules, a survival model was constructed by univariate and multivariate Cox proportional hazard regression analyses. Thereafter, we carried out Gene Ontology (GO) and the Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses. Finally, we predicted interactions between lncRNAs, miRNAs and mRNAs by TargetScan, miRDB, miRTarBase and starBase. We identified 2 lncRNAs (NORAD, XIST), 5 miRNAs (hsa‐miR‐3613, hsa‐miR‐371, hsa‐miR‐373, hsa‐miR‐32, hsa‐miR‐92) and 2 mRNAs (LYZ, PIK3AP1) for the construction of a ceRNA network, which might act as a prognostic biomarker of GBM. Combined with previous studies and our enrichment analysis results, we hypothesized that this ceRNA network affects immune activities and tumour microenvironment variations. Our research provides novel aspects to study GBM development and treatment.