Orchestrating a biomarker panel with lncRNAs and mRNAs for predicting survival in pancreatic ductal adenocarcinoma

The low survival of patients with pancreatic ductal adenocarcinoma (PDAC) makes the treatment of this disease one of the most challenging task in modern medicine. Here, by mining a large‐scale cancer genome atlas data set of pancreatic cancer tissues, we identified 21 long noncoding RNAs (lncRNAs) that significantly associated with overall survival in patients with PDAC (P < .01). Further analysis revealed that 8 lncRNAs turned out to be independently correlated with patients’ overall survival, and the risk score could be calculated based on their expression. To obtain a better predicting power, we integrated lncRNA data with a total of 410 differently expressed messenger RNAs (mRNAs) screened from PDAC and normal tissues in gene expression omnibus (GEO) database. The integration resulted in a much better panel including 8 lncRNAs (RP3.470B24.5, CTA.941F9.9, RP11.557H15.3, LINC00960, AP000479.1, LINC00635, LINC00636, and AC073133.1) and 8 mRNAs (DHRS9, ONECUT1, OR8D4, MT1M, TCN1, MMP9, DPYSL3, and TTN) to predict prognosis. A functional evaluation showed that these lncRNAs might play roles in pancreatic secretion, cell adhesion, and proteolysis. Using normal and pancreatic cancer cell lines, we confirmed that a majority of identified lncRNAs and mRNAs showed altered expressions in pancreatic cancer cells. Especially, LINC01589, LINC00960, TCN1, and MT1M showed a profoundly increased expression in pancreatic cancer cells, which suggests their potentially important role in pancreatic cancer. The results of our work indicate that lncRNAs have vital roles in PADC and provide new insights to integrate multiple kinds of markers in clinical practices.     

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.     

Identification of LINC01615 as potential metastasis-related long noncoding RNA in hepatocellular carcinoma

Hepatocellular carcinoma is one of the most prevalent and fatal cancers. Studying the long noncoding RNA (lncRNA) alterations in hepatocellular carcinoma may lead to new therapeutic strategies. We checked whether there were correlations between The Cancer Genome Atlas expression profiles of the differentially expressed lncRNAs and their DNA methylation status or the copy number variations for hepatocellular carcinoma. We obtained 41 lncRNAs that were differentially expressed between tumor and normal samples, and their DNA methylation status was negatively correlated with the expression levels. We identified five lncRNAs that were recurrently amplified or deleted in tumor samples, but none of them were associated with the messenger RNA (mRNA) expression levels. To obtain the biological function of these lncRNAs, the coexpressed mRNAs in the hepatocellular carcinoma were figured out. A total of 10 lncRNAs were highly correlated with at least one gene. Six out of the ten lncRNAs were already known to be related with cancer previously. LINC01615 had 72 coexpressed genes, and we carried out the gene ontology (GO) term enrichment for these protein-coding genes. The results suggested that these lncRNAs were associated with extracellular matrix organization. To summarize, we identified 41 potentially cancer-related lncRNAs. In particular, we proposed that LINC01615 potentially affected the extracellular matrix and had further impacts on the metastasis of hepatocellular carcinoma.     

Identification of short-lived long non-coding RNAs as surrogate indicators for chemical stress response

Abiotic and biotic stressors in human cells are often a result of sudden and/or frequent changes in environmental factors. The molecular response to stress involves elaborate modulation of gene expression and is of homeostatic, ecological, and evolutionary importance. Although attention has primarily focused on signaling pathways and protein networks, long non-coding RNAs (ncRNAs) are increasingly involved in the molecular mechanisms associated with responses to cellular stresses. We identified six novel short-lived long ncRNAs (MIR22HG, GABPB-AS1, LINC00152, IDI2-AS1, SNHG15, and FLJ33630) that responded to chemical stressors (cisplatin, cycloheximide, and mercury (II) oxide) in HeLa Tet-off cells. Our results indicate that short-lived long ncRNAs respond to general and specific chemical stressors. The expression levels of the short-lived long ncRNAs were elevated because of prolonged decay rates in response to chemical stressors and interruption of RNA degradation pathways. We propose that these long ncRNAs have the potential to be surrogate indicators of cellular stress responses.     

LINC00657 expedites neuropathic pain development by modulating miR-136/ZEB1 axis in a rat model

Long non-coding RNAs (lncRNAs) are involved in the progression of several diseases. The interactions among lncRNAs, microRNA (miRNAs) or their targeting genes are reported to play crucial roles in the development of diseases. LINC00657 is observed to be upregulated in several cancers. However, the biological role of LINC00657 in neuropathic pain progress is unclear. Hence, in our study, we aimed to investigate the function of LINC00657 in neuropathic pain development. A chronic constriction injury (CCI) rat model was established, and we found that LINC00657 was greatly increased in CCI rats associated with a decrease of miR-136. Inhibition of LINC00657 suppressed neuropathic pain via alleviating mechanical and thermal hyperalgesia. In addition, miR-136 overexpression can also inhibit the neuropathic pain development. MiR-136 was predicted to serve as a miRNA target of LINC00657, and dual-luciferase reporter assay confirmed the correlation between LINC00657 and miR-136. Moreover, we observed that the decrease of LINC00657 was able to inhibit the neuroinflammation of CCI rats by targeting expression of cyclooxygenase-2, tumor necrosis factor-α and interleukin-1β while miR-136 inhibitors reversed this phenomenon. Next, by using bioinformatics analysis, ZEB1 was predicted as a direct target of miR-136, and miR-136 could negatively modulate ZEB1 expression. Besides these, ZEB1 was remarkably increased in the CCI rats. Knockdown of ZEB1 can inhibit neuropathic pain development, while miR-136 inhibitors can reverse it. In conclusion, it was implied that LINC00657 can induce the neuropathic pain development via regulating miR-136/ZEB1 axis.     

The oncogenic potentials and diagnostic significance of long non‐coding RNA LINC00310 in breast cancer

Recent studies have revealed that long non‐coding RNAs (lncRNAs) are involved in different physiological processes and human diseases. However, to date, the function and overall clinical significance of the vast majority of lncRNAs in breast cancer remain largely unexplored. Here, we focused on LINC00310 by interrogating the breast invasive carcinoma data set of the Cancer Genome Atlas (TCGA). The results showed that LINC00310 was increased as breast cancer progressed, and the deregulation of LINC00310 was significantly associated with patients’ survival. Experiments with knockout (KO) approach by CRISPR/Cas9 system and the subsequent rescue experiments revealed that LINC00310 promoted cell proliferation by regulating c‐Myc expression in vitro. Nude mouse xenograft assay demonstrated that LINC00310 KO significantly suppressed tumour growth in vivo. Furthermore, we found that serum LINC00310 expression was significantly up‐regulated in patients with breast cancer, and receiver operating characteristic (ROC) curve analysis indicated that LINC00310 had a powerful capability of distinguishing patients with breast cancer from healthy individuals (the area under curve 0.828). Taken together, these results provide a more intuitive approach to explore the clinical relevance and functional roles of lncRNAs. As a result, lncRNAs, such as LINC00310, may be used in clinical applications as circulating markers for breast cancer.     

Downregulation of LINC00894-002 Contributes to Tamoxifen Resistance by Enhancing the TGF-β Signaling Pathway

Tamoxifen is a widely used personalized medicine for estrogen receptor (ER)-positive breast cancer, but approximately 30% of patients receiving the treatment relapse due to tamoxifen resistance (TamR). Recently, several reports have linked lncRNAs to cancer drug resistance. However, the role of lncRNAs in TamR is unclear. To identify TamR-related lncRNAs, we first used a bioinformatic approach to predict whether they have connection with known TamR-associated genes by starBase v2.0 and divided them into two groups. Group A contains lncRNAs that connect with known TamR genes and group B contains lncRNAs that show no predicted interaction. Among the 12 lncRNAs in group A, 58.3% of them are either up- or downregulated in MCF-7/TamR cells compared to the sensitive cells. In contrast, the expression levels of all group B lncRNAs are not changed in MCF-7/TamR cells. LINC00894-002 exhibits the most sophisticated network pattern and is the most downregulated lncRNA in MCF-7/TamR cells. Moreover, we find that LINC00894-002 is directly upregulated by ERα. Knocking down LINC00894-002 downregulates expression of miR-200a-3p and miR-200b-3p, upregulates the expression of TGF-β2 and ZEB1, and finally contributes to TamR. Herein, we report the first case of an inhibitory lncRNA against TamR through the miR-200-TGF-β2-ZEB1 signaling pathway.     

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.     

The differential formation of the LINC-mediated perinuclear actin cap in pluripotent and somatic cells

The actin filament cytoskeleton mediates cell motility and adhesion in somatic cells. However, whether the function and organization of the actin network are fundamentally different in pluripotent stem cells is unknown. Here we show that while conventional actin stress fibers at the basal surface of cells are present before and after onset of differentiation of mouse (mESCs) and human embryonic stem cells (hESCs), actin stress fibers of the actin cap, which wrap around the nucleus, are completely absent from undifferentiated mESCs and hESCs and their formation strongly correlates with differentiation. Similarly, the perinuclear actin cap is absent from human induced pluripotent stem cells (hiPSCs), while it is organized in the parental lung fibroblasts from which these hiPSCs are derived and in a wide range of human somatic cells, including lung, embryonic, and foreskin fibroblasts and endothelial cells. During differentiation, the formation of the actin cap follows the expression and proper localization of nuclear lamin A/C and associated linkers of nucleus and cytoskeleton (LINC) complexes at the nuclear envelope, which physically couple the actin cap to the apical surface of the nucleus. The differentiation of hESCs is accompanied by the progressive formation of a perinuclear actin cap while induced pluripotency is accompanied by the specific elimination of the actin cap, and that, through lamin A/C and LINC complexes, this actin cap is involved in progressively shaping the nucleus of hESCs undergoing differentiation. While, the localization of lamin A/C at the nuclear envelope is required for perinuclear actin cap formation, it is not sufficient to control nuclear shape.     

A novel autophagy-related lncRNA prognostic risk model for breast cancer

Long non-coding RNAs (lncRNAs) are well known as crucial regulators to breast cancer development and are implicated in controlling autophagy. LncRNAs are also emerging as valuable prognostic factors for breast cancer patients. It is critical to identify autophagy-related lncRNAs with prognostic value in breast cancer. In this study, we identified autophagy-related lncRNAs in breast cancer by constructing a co-expression network of autophagy-related mRNAs-lncRNAs from The Cancer Genome Atlas (TCGA). We evaluated the prognostic value of these autophagy-related lncRNAs by univariate and multivariate Cox proportional hazards analyses and eventually obtained a prognostic risk model consisting of 11 autophagy-related lncRNAs (U62317.4, LINC01016, LINC02166, C6orf99, LINC00992, BAIAP2-DT, AC245297.3, AC090912.1, Z68871.1, LINC00578 and LINC01871). The risk model was further validated as a novel independent prognostic factor for breast cancer patients based on the calculated risk score by Kaplan-Meier analysis, univariate and multivariate Cox regression analyses and time-dependent receiver operating characteristic (ROC) curve analysis. Moreover, based on the risk model, the low-risk and high-risk groups displayed different autophagy and oncogenic statues by principal component analysis (PCA) and Gene Set Enrichment Analysis (GSEA) functional annotation. Taken together, these findings suggested that the risk model of the 11 autophagy-related lncRNAs has significant prognostic value for breast cancer and might be autophagy-related therapeutic targets in clinical practice.     

LINC01116 modulates EMT process via binding with AGO1 mRNA in oesophageal squamous cell carcinoma

Recent researches have uncovered that long non-coding RNAs (lncRNAs) are closely correlated with the development of different diseases, while biological functions and hidden molecular mechanisms of antisense lncRNAs in oesophageal squamous cell carcinoma (OSCC) remain unclear. Here, we identified upregulation of LINC01116 in RNA sequencing data, online database, and in OSCC and intraepithelial neoplasia (IEN) specimens. Functionally, LINC01116 facilitates OSCC advancement and metastasis in vitro and vivo. Mechanistically, elevated expression of LINC01116 in OSCC cells other than tumor stroma and cytoplasmic enables it to activate AGO1 expression via complementary binding with AGO1 mRNA to facilitate EMT process of OSCC.     

LINC00665 induces gastric cancer progression through activating Wnt signaling pathway

Long noncoding RNAs (lncRNAs) are recently recognized as noteworthy regulators of different tumors, counting gastric cancer (GC). Lately, long intergenic noncoding RNA (LINC) 00665 has been verified to display significant parts in several cancers. Be that as it may, its role and mechanism in GC movement still stay uninvestigated. As of now, we observed LINC00665 was obviously GC cells (MKN28, BGC-823, SGC7-901, AGS, HGC-27) in comparison to GES-1 cells, which was identified as human normal gastric epithelial cells. Then, LINC00665 was obviously downregulated in GC cells including AGS and BGC-823 cells. Loss of LINC00665 greatly repressed AGS and BGC-823 cell survival and cell expansion. Moreover, GC cell apoptosis was significantly induced by the loss of LINC00665. For another, we found that the GC cell cycle was also captured in G1 and G2 phases. The experiments on cell migration and invasion indicated that knockdown of LINC00665 restrained GC cell migration and invasion. Modifications in Wnt signaling are closely associated with the development of cancers. Here, we found that Wnt signaling was significantly inactivated by the silence of LINC00665 in GC cells. β-catenin and cyclinD1 were restrained whereas GSK-3β was induced by the inhibition of LINC00665 in GC cells. Furthermore, we confirmed the impact of LINC00665 in vivo using xenograft models. Taken these together, we indicated that LINC00665 could function as a novel biomarker in GC progression.     

RUNX2 stabilization by long non-coding RNAs contributes to hypertrophic changes in human chondrocytes

Background: Chondrocyte hypertrophy has been implicated in endochondral ossification and osteoarthritis (OA). In OA, hypertrophic chondrocytes contribute to the destruction and focal calcification of the joint cartilage. Although studies in this field have remarkably developed the modulation of joint inflammation using gene therapy and regeneration of damaged articular cartilage using cell therapy, studies that can modulate or prevent hypertrophic changes in articular chondrocytes are still lacking.Methods: In vitro hypertrophic differentiation and inflammation assays were conducted using human normal chondrocyte cell lines, TC28a2 cells. Human cartilage tissues and primary articular chondrocytes were obtained from OA patients undergoing total knee arthroplasty. Long non-coding RNAs (lncRNAs), LINC02035 and LOC100130207, were selected through RNA-sequencing analysis using RNAs extracted from TC28a2 cells cultured in hypertrophic medium. The regulatory mechanism was evaluated using western blotting, real-time quantitative polymerase chain reaction, osteocalcin reporter assay, RNA-immunoprecipitation (RNA-IP), RNA-in situ hybridization, and IP.Results: LncRNAs are crucial regulators of various biological processes. In this study, we identified two important lncRNAs, LINC02035 and LOC100130207, which play important roles in hypertrophic changes in normal chondrocytes, through RNA sequencing. Interestingly, the expression level of RUNX2, a master regulator of chondrocyte hypertrophy, was regulated at the post-translational level during hypertrophic differentiation of the normal human chondrocyte cell line, TC28a2. RNA-immunoprecipitation proved the potential interaction between RUNX2 protein and both lncRNAs. Knockdown (KD) of LINC02035 or LOC100130207 promoted ubiquitin-mediated proteasomal degradation of RUNX2 and prevented hypertrophic differentiation of normal chondrocyte cell lines, whereas overexpression of both lncRNAs stabilized RUNX2 protein and generated hypertrophic changes. Furthermore, the KD of the two lncRNAs mitigated the destruction of important cartilage matrix proteins, COL2A1 and ACAN, by hypertrophic differentiation or inflammatory conditions. We also confirmed that the phenotypic changes raised by the two lncRNAs could be rescued by modulating RUNX2 expression. In addition, the KD of these two lncRNAs suppressed hypertrophic changes during chondrogenic differentiation of mesenchymal stem cells.Conclusion: Therefore, this study suggests that LINC02035 and LOC100130207 contribute to hypertrophic changes in normal chondrocytes by regulating RUNX2, suggesting that these two novel lncRNAs could be potential therapeutic targets for delaying or preventing OA development, especially for preventing chondrocyte hypertrophy.