Construction and analysis of the abnormal lncRNA–miRNA–mRNA network in hypoxic pulmonary hypertension

The incidence of hypoxic pulmonary hypertension (HPH) is increasing. Accumulating evidence suggests that long noncoding RNAs (lncRNAs) play an important role in HPH, but the functions and mechanism have yet to be fully elucidated. In the present study, we established a HPH rat model with 8 h of hypoxia exposure (10% O₂) per day for 21 days. High-throughput sequencing identified 60 differentially expressed (DE) lncRNAs, 20 DE miRNAs and 695 DE mRNAs in rat lung tissue. qRT-PCR verified the accuracy of the results. The DE mRNAs were significantly enriched in immune response, inflammatory response, leukocyte migration, cell cycle, cellular response to interleukin-1, IL-17 signalling pathway, cytokine–cytokine receptor interaction and Toll-like receptor signalling pathway. According to the theory of competing endogenous RNA (ceRNA) networks, lncRNA–miRNA–mRNA network was constructed by Cytoscape software, 16 miRNAs and 144 mRNAs. The results suggested that seven DE lncRNAs (Ly6l, AABR07038849.2, AABR07069008.2, AABR07064873.1, AABR07001382.1, AABR07068161.1 and AABR07060341.2) may serve as molecular sponges of the corresponding miRNAs and play a major role in HPH.     

Genome-wide identification of a competing endogenous RNA network in cholangiocarcinoma

Cholangiocarcinoma (CCA) is the second widespread liver tumor with relatively poor survival. Increasing evidence in recent studies showed long noncoding RNAs (lncRNAs) exert a crucial impact on the development and progression of CCA based on the mechanism of competing endogenous RNAs (ceRNAs). However, functional roles and regulatory mechanisms of lncRNA-regulated ceRNA in CCA, are only partially understood. The expression profile of messenger RNAs (mRNAs), lncRNAs, and microRNAs (miRNAs) downloaded from The Cancer Genome Atlas were comprehensively investigated. Differential expression of these three types of RNA between CCA and corresponding precancerous tissues were screened out for further analysis. On the basis of interactive information generated from miRDB, miRTarBase, TargetScan, and miRcode public databases, we then constructed an mRNA-miRNA-lncRNA regulatory network. Kyoto Encyclopedia of Genes and Genomes and Gene Ontology analyses were conducted to identify the biological function of the ceRNA network involved in CCA. As a result, 2883 mRNAs, 136 miRNAs, and 993 lncRNAs were screened out as differentially expressed RNAs in CCA. In addition, a ceRNA network in CCA was constructed, composing of 50 up and 27 downregulated lncRNAs, 14 up and 7 downregulated miRNAs, 29 up and 25 downregulated mRNAs. Finally, gene set enrichment and pathway analysis indicated our CCA-specific ceRNA network was related with cancer-related pathway and molecular function. In conclusion, our research identified a novel lncRNA-related ceRNA network in CCA, which might act as a potential therapeutic target for patients with CCA.     

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.     

Construction of competitive endogenous RNA network reveals regulatory role of long non‐coding RNAs in type 2 diabetes mellitus

Increasing epidemic of type 2 diabetes mellitus (T2DM) and its comorbidities makes it urgent to understand the pathogenesis and regulatory mechanism. However, little is known about the regulatory role of lncRNAs in diabetes. Here, we constructed a T2DM‐related competitive endogenous RNA (ceRNA) network (DMCN) to explore biological function of lncRNAs during the development of diabetes mellitus. This network contained 351 nodes including 98 mRNAs, 86 microRNAs and 167 lncRNAs. Functional analysis showed that the mRNAs in DMCN were annotated into some diabetes‐related pathways. Furthermore, mTOR‐centred subnetwork was extracted and ncRNA‐involved mTOR pathway was established. Finally, we validated that NEAT1 was potentially communicated with mTOR signalling target protein mLST8 via the association with miR‐181b. These findings provide significant insight into lncRNA regulatory network in T2DM.     

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.     

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 LINC00339 promotes laryngeal squamous cell carcinoma cell proliferation and invasion via sponging miR-145

Laryngeal squamous cell carcinoma (LSCC) is a very common neoplasm of the head and neck in the world. Long noncoding RNAs play key roles in cell infiltration, fate, apoptosis, and invasion. However, the functional role and expression of LINC00339 remains unclear in LSCC. In this study, we showed that the expression level of LINC00339 was upregulated in LSCC tissues and cell lines. LINC00339 silencing suppressed the proliferation, invasion, and epithelial-mesenchymal transition (EMT) progression of LSCC cells. In addition, we showed that LINC00339 acted as a sponge of miR-145, and LINC00339 silencing promoted the expression of miR-145 in Hep2 cell. Furthermore, the expression of miR-145 was lower in LSCC tissues than in their paired normal samples and the miR-145 expression level was negatively correlated with LINC00339 expression in LSCC tissues. The knockdown of miR-145 promoted the proliferation, invasion, and EMT progression of LSCC cells. Finally, we indicated that LINC00339 silencing inhibited the proliferation, invasion, and EMT progression of LSCC cells by suppressing the miR-145 expression. These data suggested that LINC00339 acted as an oncogene in the development of LSCC, partly by regulating the miR-145 expression.     

Overexpression of ZEB2‐AS1 promotes epithelial‐to‐mesenchymal transition and metastasis by stabilizing ZEB2 mRNA in head neck squamous cell carcinoma

The long noncoding RNAs (lncRNAs) have been increasingly appreciated as key players underlying tumourigenesis and hold great potentials as prognostic biomarkers and therapeutic targets. However, their roles in head neck squamous cell carcinoma (HNSCC) have remained incompletely known. Here, we sought to reveal the oncogenic roles and clinical significance of a tumour‐associated lncRNA, zinc finger E‐box binding homeobox 2 antisense RNA 1 (ZEB2‐AS1), in HNSCC. ZEB2‐AS1 was aberrantly overexpressed in a fraction of HNSCC samples. Its overexpression significantly associated with large tumour size, cervical node metastasis and reduced overall and disease‐free survival. Antisense oligonucleotides (ASO)‐mediated ZEB2‐AS1 depletion markedly inhibited cell proliferation, migration and invasion while triggered apoptosis in HNSCC cells in part via modulating ZEB2 mRNA stability. Enforced overexpression of ZEB2 largely attenuated the phenotypic changes resulted from ZEB2‐AS1 inhibition except the impaired cell proliferation. In addition, ZEB2‐AS1 was required for TGF‐β1‐induced epithelial‐mesenchymal transition (EMT) in vitro. Significantly reduced tumour growth and lung metastasis were observed in ZEB2‐AS1‐depleted cells in HNSCC xenograft animal models. Taken together, our findings reveal that overexpression of ZEB2‐AS1 associates with tumour aggressiveness and unfavourable prognosis by serving as a putative oncogenic lncRNA and a novel prognostic biomarker in HNSCC.     

An 11-lncRNA expression could be potential prognostic biomarkers in head and neck squamous cell carcinoma

The aim of our study is to construct the competing endogenous RNA (ceRNA) network of head and neck squamous cell carcinoma (HNSCC) and identify key long noncoding RNAs (lncRNAs) to predict prognosis. The genes whose expression were differentially in HNSCC and normal tissues were explored by the Cancer Genome Atlas database. The ceRNA network was constructed by the Cytoscape software. The lncRNAs which could estimate the overall survival were explored from Cox proportional hazards regression. There are 1997, 589, and 82 mRNAs, lncRNAs, and miRNAs whose expression were statistically significant different, respectively. Then, the network between miRNA and mRNA or miRNA and lncRNA was constructed by miRcode, miRDB, TargetScan, and miRanda. Five mRNAs, 10 lncRNAs, and 3 miRNAs were associated with overall survival. Then, 11-lncRNAs were found to be prognostic factors. Therefore, our research analyzed the potential signature of novel 11-lncRNA as candidate prognostic biomarker from the ceRNA network for patients with HNSCC.     

Prediction of head and neck squamous cell carcinoma survival based on the expression of 15 lncRNAs

Recent evidence suggests that long noncoding RNAs (lncRNAs) are essential regulators of many cancer-related processes, including cancer cell proliferation, invasion, and migration. There is thus a reason to believe that the detection of lncRNAs may be useful as a diagnostic and prognostic strategy for cancer detection, however, at present no effective genome-wide tests are available for clinical use, constraining the use of such a strategy. In this study, we performed a comprehensive assessment of lncRNAs expressed in samples in the head and neck squamous cell carcinoma (HNSCC) cohort available in The Cancer Genome Atlas database. A risk score (RS) model was constructed based on the expression data of these 15 lncRNAs in the validation data set of HNSCC patients and was subsequently validated in validation data set and the entire data set. We were able to stratify patients into high- and low-risk categories, using our lncRNA expression panel to determine an RS, with significant differences in overall survival (OS) between these two groups in our test set (median survival, 1.863 vs. 5.484 years; log-rank test, p < 0.001). We were able to confirm the predictive value of our 15-lncRNA signature using both a validation data set and a full data set, finding our signature to be reproducible and effective as a means of predicting HNSCC patient OS. Through the multivariate Cox regression and stratified analyses, we were further able to confirm that the predictive value of this RS was independent of other predictive factors such as clinicopathological parameters. The Gene set enrichment analysis revealed potential functional roles for these 15 lncRNAs in tumor progression. Our findings indicate that an RS established based on a panel of lncRNA expression signatures can effectively predict OS and facilitate patient stratification in HNSCC.     

CERNA2: A predictor for clinical progression and poor prognosis in cervical carcinoma

Competing long noncoding RNA 2 (lncRNA 2) for microRNA let-7b (CERNA2) has emerged as an important regulator of tumorigenesis and cancer progression but the clinical value and regulatory function of CERNA2 is yet to be investigated in cervical carcinoma. In our study, we found the CERNA2 expression was obviously increased in cervical carcinoma tissues compared with adjacent normal cervical tissues. In addition, we observed that metastatic lymph nodes exhibited high levels of CERNA2 expression in contrast to primary cervical carcinoma tissues. Furthermore, high CERNA2 expression was associated with advanced clinical stage, lymph node metastasis, distant metastasis poor histological grade, and short overall survival in cervical carcinoma patients. Moreover, high CERNA2 expression acted as an independent unfavorable predictor for overall survival in cervical carcinoma patients. The cell migration and invasion assays in vitro suggested that knockdown of CERNA2 remarkably inhibited cell migration and invasion in cervical carcinoma. In conclusion, CERNA2 functions as an oncogenic lncRNA and may be as a potential therapeutic target in cervical carcinoma.     

Long noncoding RNA HOXD-AS1 induces epithelial-mesenchymal transition in breast cancer by acting as a competing endogenous RNA of miR-421

Breast cancer (BCa) is the most common malignant tumor in females. Long noncoding RNAs (lncRNAs) are deregulated in many types of human cancers, including BCa. The purpose of the present study was to examine the expression profile and biological role of HOXD cluster antisense RNA 1 (HOXD-AS1) in BCa. Our results revealed that HOXD-AS1 was upregulated in BCa tissues and cell lines, and high HOXD-AS1 expression was correlated with aggressive clinicopathological characteristics of BCa patients. Further gain-of-function and loss-of-function analysis showed that HOXD-AS1 overexpression promoted, whereas HOXD-AS1 knockdown inhibited BCa cell proliferation, cell cycle progression, migration, and invasion, indicating that HOXD-AS1 may function as a novel oncogene in BCa. Mechanistically, HOXD-AS1 could activate epithelial-mesenchymal transition (EMT) in BCa cells. We further proved that HOXD-AS1 might serve as a competing endogenous RNA of miR-421 in BCa cells, and miR-421 was downregulated and negatively correlated with HOXD-AS1 expression in BCa tissues. Besides, we confirmed that SOX4, a master regulator of EMT, was a direct target gene of miR-421. Further, rescue experiments suggested that miR-421 overexpression partly abrogated the oncogenic role of HOXD-AS1 in BCa cells. Therefore, we shed light on that HOXD-AS1/miR-421/SOX4 axis may be considered as a novel therapeutic target for the treatment of BCa patients.