Comprehensive analysis of pseudogene HSPB1P1 and its potential roles in hepatocellular carcinoma

The incidence and mortality rate of hepatocellular carcinoma (HCC) nowadays is still at high levels. The regulatory roles of pseudogene in cancers have been gradually recognized in recent years. However, comprehensive investigation of abnormally expressed pseudogene and related mechanisms in HCC remains lacking. GSE124535 dataset was used to identify differentially expressed pseudogenes in HCC tissues compared with normal tissues. Prognostic value of these differentially expressed pseudogenes was analyzed at GEPIA. StarBase used to analyze microRNAs (miRNAs) can bind with pseudogene, while the targets for these miRNAs were analyzed at miRTarBase. Protein–protein interaction (PPI) network was then established for miRNA targets, after that hub genes were selected. Expression correlation of pseudogene and hub genes was analyzed at StarBase. In total, 16 upregulated and 17 downregulated pseudogenes were identified. Pseudogene HSPB1P1 was identified abnormally expressed in 20 types of human cancers and could be used as an indicator for poorer overall survival of patients with HCC. Functional analyses showed that HSPB1P1 was strongly correlated with signaling pathways related to cancer progression. Further studied revealed that HSPB1P1 could direct regulate the EZH2 expression in HCC. In summary, our study indicated that HSPB1P1 was a predictor for poorer overall survival of patients with HCC and may be potential therapeutic target against HCC.     

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.     

Identification of lncRNA–miRNA–mRNA regulatory network associated with epithelial ovarian cancer cisplatin-resistant

To construct a long noncoding RNA (lncRNA)–microRNA (miRNA)–messenger RNA (mRNA) regulatory network related to epithelial ovarian cancer (EOC) cisplatin-resistant, differentially expressed genes (DEGs), differentially expressed lncRNAs (DELs), and differentially expressed miRNAs (DEMs) between MDAH and TOV-112D cells lines were identified. Gene Ontology (GO) analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis were conducted to analyze the biological functions of DEGs. Downstream mRNAs or upstream lncRNAs for miRNAs were analyzed at miRTarBase 7.0 or DIANA-LncBase V2, respectively. A total of 485 significant DEGs, 85 DELs, and 5 DEMs were identified. Protein–protein interaction (PPI) network of DEGs contrains 81 nodes and 141 edges was constructed, and 25 hub genes related to EOC cisplatin-resistant were identified. Subsequently, a lncRNA–miRNA–mRNA regulatory network contains 4 lncRNAs, 4 miRNAs, and 35 mRNAs was established. Taken together, our study provided evidence concerning the alteration genes involved in EOC cisplatin-resistant, which will help to unravel the mechanisms underlying drug resistant.     

Cancer Specific Long Noncoding RNAs Show Differential Expression Patterns and Competing Endogenous RNA Potential in Hepatocellular Carcinoma

Long noncoding RNAs (lncRNAs) regulate gene expression by acting with microRNAs (miRNAs). However, the roles of cancer specific lncRNA and its related competitive endogenous RNAs (ceRNA) network in hepatocellular cell carcinoma (HCC) are not fully understood. The lncRNA profiles in 372 HCC patients, including 372 tumor and 48 adjacent non-tumor liver tissues, from The Cancer Genome Atlas (TCGA) and NCBI GEO omnibus (GSE65485) were analyzed. Cancer specific lncRNAs (or HCC related lncRNAs) were identified and correlated with clinical features. Based on bioinformatics generated from miRcode, starBase, and miRTarBase, we constructed an lncRNA-miRNA-mRNA network (ceRNA network) in HCC. We found 177 cancer specific lncRNAs in HCC (fold change ≥ 1.5, P < 0.01), 41 of them were also discriminatively expressed with gender, race, tumor grade, AJCC tumor stage, and AJCC TNM staging system. Six lncRNAs (CECR7, LINC00346, MAPKAPK5-AS1, LOC338651, FLJ90757, and LOC283663) were found to be significantly associated with overall survival (OS, log-rank P < 0.05). Collectively, our results showed the lncRNA expression patterns and a complex ceRNA network in HCC, and identified a complex cancer specific ceRNA network, which includes 14 lncRNAs and 17 miRNAs in HCC.     

Identification and expression analysis of lncRNA in seven organs of Rhinopithecus roxellana

Long non-coding RNA (lncRNA) represents a new direction to identify expression profiles and regulatory mechanisms in various organisms. Here, we report the first dataset of lncRNAs of the golden snub-nosed monkey (GSM), including 12,557 putative lncRNAs identified from seven organs. Compared with mRNA, GSM lncRNA had fewer exons and isoforms, and longer length. LncRNA showed more obvious tissue-specific expression than mRNA. However, for the top ten most abundant genes in each organ, mRNAs expression was more tissue-specific than lncRNAs. By identification of specifically expressed lncRNAs and mRNAs in each organ, it indicates that the expression of SEG-lncRNA (specifically expressed lncRNA) and SEG-mRNA (specifically expressed mRNA) had high correlation. In particular, combined our lncRNA and mRNA data, we identified 92 heart SEG-lncRNAs targeted ten mRNA genes in the oxidative phosphorylation pathway and upregulated the expression of these target genes such as ND4, ATP6, and ATP8. These may contribute to GSM adaption to its high-elevation environment. We also identified 171 liver SEG-lncRNAs, which targeted 27 genes associated with the metabolism of xenobiotics and leaded to high expression of these target genes in liver. These lncRNAs may play important roles in GSM adaptation to a folivory diet.

     

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.     

Reduced LINC00551 expression promotes proliferation and invasion of esophageal squamous cancer by increase in HSP27 phosphorylation

Esophageal squamous cell carcinoma (ESCC) is one of the deadliest cancers, and long noncoding RNAs (lncRNAs) regulate gene expression or activities. This study investigated the role of lncRNA LINC00551 in ESCC development and progression. Three paired ESCC and normal tissues were subjected to next‐generation sequencing and we identified 82 upregulated and 60 downregulated lncRNAs, including LINC00551, which was confirmed to markedly downregulated in 78 ESCC tissues and in the Gene Expression Profiling Interactive Analysis data set. Downregulated LINC00551 expression was associated with lymph node metastasis, advanced TNM stage, and tumor size. Moreover, downregulated LINC00551 expression was also associated with poor progression‐free survival and overall survival of ESCC patients. In vitro and in vivo, LINC00551 overexpression inhibited ESCC cell proliferation and invasion, whereas knockdown of LINC00551 expression promoted ESCC cell proliferation and invasion. RNA pull‐down and mass spectrometry assays identified the potential LINC00551 binding proteins, and HSP27 was a promising LINC00551 targeting proteins after RNA immunoprecipitation assay. At the protein level, LINC00551 bound to and decreased HSP27 phosphorylation, and in turn, downregulated ESCC cell proliferation and invasion. The current study demonstrated the functional significance of LINC00551 in ESCC development, progression, and prognosis. Further study will assess LINC00551 as a novel prognostic marker or therapeutic target for ESCC.     

Systematic identification of lncRNAs and circRNAs-associated ceRNA networks in human lumbar disc degeneration

ABSTRACT

Lumbar disc degeneration (LDD) is a common cause of low back and neck pain. The molecular mechanisms underlying LDD, however, are unclear. Noncoding RNAs have been reported to participate in human diseases. We investigated a series of public datasets (GSE67566, GSE56081 and GSE63492) and identified 568 mRNAs, 55 microRNAs (miRNAs), 765 long noncoding RNAs (lncRNAs), and 586 circular RNAs (circRNAs) that were expressed differently in LDD than in normal discs. We constructed lncRNAs and circRNAs regulated competing endogenous RNAs (ceRNA) networks in LDD. Four lncRNAs, DANCR, CASK-AS1, SCARNA2, and LINC00638), and three circRNAs, hsa_circ_0005139, hsa_circ_0037858, and hsa_circ_0087890, were identified as key regulators of LDD progression. We found that hsa-miR-486-5p regulated the crosstalk among circRNA hsa_circ_0000189, lncRNA DANCR and 6 mRNAs, PYCR2, TOB1, ARHGAP5, RBPJ, CD247, SLC34A1. Gene ontology (GO) analysis demonstrated that these differently expressed lncRNAs and circRNAs were involved in cellular component organization or biogenesis, gene expression and negative regulation of metabolic processes. Our findings provide useful information for exploring new mechanisms for LDD and candidates for therapeutic targets.     

Long Non-coding RNA URHC Regulates Cell Proliferation and Apoptosis via ZAK through the ERK/MAPK Signaling Pathway in Hepatocellular Carcinoma

Long non-coding RNAs (lncRNAs) have previously been implicated in human disease states, especially cancer. Although the aberrant expression of lncRNAs has been observed in cancer, the biological functions and molecular mechanisms underlying aberrantly expressed lncRNAs in hepatocellular carcinoma (HCC) have not been widely established. In the present study, we investigated a novel lncRNA, termed URHC (up-regulated in hepatocellular carcinoma), and evaluated its role in the progression of HCC. Expression profiling using a lncRNA microarray revealed that URHC was highly expressed in 3 HCC cell lines compared to normal hepatocytes. Quantitative real-time polymerase chain reaction (qRT-PCR) analyses confirmed that URHC expression was increased in hepatoma cells and HCC tissues. Moreover, using qRT-PCR, we confirmed that URHC expression was up-regulated in 30 HCC cases (57.7%) and that its higher expression was correlated with poor overall survival. We further demonstrated that URHC inhibition reduced cell proliferation and promoted apoptosis. We hypothesize that URHC may function by regulating the sterile alpha motif and leucine zipper containing kinase AZK (ZAK) gene, which is located near URHC on the same chromosome. We found that ZAK mRNA levels were down-regulated in HCC tissues and the expression levels of ZAK were negatively correlated with those of URHC in the above HCC tissues. Next, we confirmed that URHC down-regulated ZAK, which is involved in URHC-mediated cell proliferation and apoptosis. Furthermore, ERK/MAPK pathway inactivation partially accounted for URHC-ZAK-induced cell growth and apoptosis. Thus, we concluded that high URHC expression can promote cell proliferation and inhibit apoptosis by repressing ZAK expression through inactivation of the ERK/MAPK pathway. These findings may provide a novel mechanism and therapeutic targets for the treatment of HCC.     

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.     

Identification and validation of key genes associated with non-small-cell lung cancer

Non-small-cell lung cancer (NSCLC) is one of the main causes of death induced by cancer globally. However, the molecular aberrations in NSCLC patients remain unclearly. In the present study, four messenger RNA microarray datasets (GSE18842, GSE40275, GSE43458, and GSE102287) were downloaded from the Gene Expression Omnibus (GEO) database. Differentially expressed genes (DEGs) between NSCLC tissues and adjacent lung tissues were obtained from GEO2R and the overlapping DEGs were identified. Moreover, functional and pathway enrichment were performed by Funrich, while the protein–protein interaction (PPI) network construction were obtained from STRING and hub genes were visualized and identified by Cytoscape software. Furthermore, validation, overall survival (OS) and tumor staging analysis of selected hub genes were performed by GEPIA. A total of 367 DEGs (95 upregulated and 272 downregulated) were obtained through gene integration analysis. The PPI network consisted of 94 nodes and 1036 edges in the upregulated DEGs and 272 nodes and 464 edges in the downregulated DEGs, respectively. The PPI network identified 46 upregulated and 27 downregulated hub genes among the DEGs, and six (such as CENPE, NCAPH, MYH11, LRRK2, HSD17B6, and A2M) of that have not been identified to be associated with NSCLC so far. Moreover, the expression differences of the mentioned hub genes were consistent with that in lung adenocarcinoma and lung squamous cell carcinoma in the TCGA database. Further analysis showed that all the six hub genes were associated with tumor staging except MYH11, while only the upregulated DEG CENPE was associated with the worse OS of patients with NSCLC. In conclusion, the current study showed that CENPE, NCAPH, MYH11, LRRK2, HSD17B6, and A2M might be the key genes contributed to tumorigenesis or tumor progression in NSCLC, further functional study is needed to explore the involved mechanisms.