Comprehensive analysis of a ceRNA network reveals potential prognostic cytoplasmic lncRNAs involved in HCC progression

The aberrant expression of long noncoding RNAs (lncRNAs) has drawn increasing attention in the field of hepatocellular carcinoma (HCC) biology. In the present study, we obtained the expression profiles of lncRNAs, microRNAs (miRNAs), and messenger RNAs (mRNAs) in 371 HCC tissues and 50 normal tissues from The Cancer Genome Atlas (TCGA) and identified hepatocarcinogenesis-specific differentially expressed genes (DEGs, log fold change ≥ 2, FDR < 0.01), including 753 lncRNAs, 97 miRNAs, and 1,535 mRNAs. Because the specific functions of lncRNAs are closely related to their intracellular localizations and because the cytoplasm is the main location for competitive endogenous RNA (ceRNA) action, we analyzed not only the interactions among these DEGs but also the distributions of lncRNAs (cytoplasmic, nuclear or both). Then, an HCC-associated deregulated ceRNA network consisting of 37 lncRNAs, 10 miRNAs, and 26 mRNAs was constructed after excluding those lncRNAs located only in the nucleus. Survival analysis of this network demonstrated that 15 lncRNAs, 3 miRNAs, and 16 mRNAs were significantly correlated with the overall survival of HCC patients (p < 0.01). Through multivariate Cox regression and lasso analysis, a risk score system based on 13 lncRNAs was constructed, which showed good discrimination and predictive ability for HCC patient survival time. This ceRNA network-construction approach, based on lncRNA distribution, not only narrowed the scope of target lncRNAs but also provided specific candidate molecular biomarkers for evaluating the prognosis of HCC, which will help expand our understanding of the ceRNA mechanisms involved in the early development of HCC.     

Long noncoding RNA H19 promotes vascular remodeling by sponging let-7a to upregulate the expression of cyclin D1

Vascular remodeling is mainly caused by excessive proliferation of vascular smooth muscle cells (VSMCs). Noncoding RNAs (ncRNAs) have emerged as important regulators in diverse pathological processes. Previous work has shown the functions and mechanisms of long noncoding RNA H19 (LncRNA H19) on VSMCs. As long noncoding RNAs (lncRNAs) are complex in their mechanisms of action, the aim of the study is to identify if there are any other molecular mechanisms of LncRNA H19 on VSMCs. In vivo studies demonstrated that cyclin D1 was overexpressed in neointima of balloon-injured artery. In vitro studies identified that the overexpression of LncRNA H19 promoted VSMCs proliferation and cyclin D1 upregulation. On the contrary, cellular proliferation and expression of cyclin D1 were inhibited in VSMCs after infection with let-7a. Furthermore, luciferase reporter assays and RNA pull-down assays were used to explore the regulatory mechanism, we found that LncRNA H19 functioned as a competing endogenous RNA (ceRNA) by sponging let-7a to promote the expression of the target gene cyclin D1. In conclusion, LncRNA H19 positively regulated cyclin D1 expression through directly binding to let-7a in VSMCs. Our findings provide new insight into the mechanism of LncRNA H19 in VSMCs proliferation and vascular remodeling, and further indicate the implications of LncRNA H19 in the diagnosis and treatment of vascular proliferative diseases.     

Aberrantly expressed long noncoding RNAs in human intervertebral disc degeneration: a microarray related study

Introduction

In addition to the well-known short noncoding RNAs such as microRNAs (miRNAs), increasing evidence suggests that long noncoding RNAs (lncRNAs) act as key regulators in a wide aspect of biologic processes. Dysregulated expression of lncRNAs has been demonstrated being implicated in a variety of human diseases. However, little is known regarding the role of lncRNAs with regards to intervertebral disc degeneration (IDD). In the present study we aimed to determine whether lncRNAs are differentially expressed in IDD.

Methods

An lncRNA-mRNA microarray analysis of human nucleus pulposus (NP) was employed. Bioinformatics prediction was also applied to delineate the functional roles of the differentially expressed lncRNAs. Several lncRNAs and mRNAs were chosen for quantitative real-time PCR (qRT-PCR) validation.

Results

Microarray data profiling indicated that 116 lncRNAs (67 up and 49 down) and 260 mRNAs were highly differentially expressed with an absolute fold change greater than ten. Moreover, 1,052 lncRNAs and 1,314 mRNAs were differentially expressed in the same direction in at least four of the five degenerative samples with fold change greater than two. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis for the differentially expressed mRNAs indicated a number of pathways, such as extracellular matrix (ECM)-receptor interaction. A coding-noncoding gene co-expression (CNC) network was constructed for the ten most significantly changed lncRNAs. Annotation terms of the coexpressed mRNAs were related to several known degenerative alterations, such as chondrocyte differentiation. Moreover, lncRNAs belonging to a particular subgroup were identified. Functional annotation for the corresponding nearby coding genes showed that these lncRNAs were mainly associated with cell migration and phosphorylation. Interestingly, we found that Fas-associated protein factor-1 (FAF1), which potentiates the Fas-mediated apoptosis and its nearby enhancer-like lncRNA RP11-296A18.3, were highly expressed in the degenerative discs. Subsequent qRT-PCR results confirmed the changes.

Conclusions

This is the first study to demonstrate that aberrantly expressed lncRNAs play a role in the development of IDD. Our study noted that up-regulated RP11-296A18.3 highly likely induced the over-expression of FAF1, which eventually promoted the aberrant apoptosis of disc cells. Such findings further broaden the understanding of the etiology of IDD.

Electronic supplementary material

The online version of this article (doi:10.1186/s13075-014-0465-5) contains supplementary material, which is available to authorized users.     

Expression profile of long non-coding RNAs in cervical spondylotic myelopathy of rats by microarray and bioinformatics analysis

IntroductionIt has been reported that a wide range of long non-coding RNAs (lncRNAs) are implicated in numerous diseases such as tumor, cardiopathy and neurological disorders. Identifying the differentially expressed (DE) profile of lncRNAs in cervical spondylotic myelopathy (CSM) is essential to understand the mechanisms of CSM.MethodsMicroarray assay, quantitative real-time PCR (qRT-PCR) and bioinformatics analysis were employed to reveal the DE profile and potential functions of lncRNAs in CSM.ResultsMicroarray analysis displayed the DE profiles of lncRNAs and mRNAs in rats between the CSM group and the control (CON) group. Thereinto, 1266 DE lncRNAs (738 up-regulation and 528 down-regulation) and 847 mRNAs (487 up-regulation and 360 down-regulation) with >1.1 fold change (FC) were finally identified. Moreover, 17 lncRNAs (13 up-regulation and 4 down-regulation) and 18 mRNAs (13 up-regulation and 5 down-regulation) were found deregulated by >2 FC. Further bioinformatics analysis showed the most remarkable biological processes among up-regulated RNAs contain cellular response to interferon-beta, inflammatory response and innate immune response, which may involve in CSM. Besides, related DE mRNAs of 17 DE lncRNAs in the genome were related to signaling pathway about NOD-like receptor, TNF, and apoptosis. In addition, a co-expression network of lncRNA-mRNA was established for analyzing the biological roles of lncRNAs. Among these, we found a ceRNA network related to CSM. Finally, the expressions of the DE lncRNAs and ceRNA network confirmed by qRT-PCR were in agreement with microarray data.ConclusionsOur study revealed the DE profiles of lncRNAs and mRNAs for CSM. Those dysregulated RNAs may represent potential therapeutic targets of CSM for further study.     

Analysis of long noncoding RNA expression profiles in the whole blood of neonates with hypoxic-ischemic encephalopathy

Long noncoding RNAs (lncRNAs) have been shown to participate in many biological processes. To investigate the expression profiles of lncRNAs and their potential functions in neonatal hypoxic-ischemic encephalopathy (HIE), we detected the lncRNA and messenger RNA (mRNA) expression in the peripheral blood samples from HIE patients and controls using a microarray. A total of 376 lncRNAs and 126 mRNAs were differentially expressed between the HIE and the non-HIE samples (fold change > 2). Quantitative real-time polymerase chain reaction was used to validate the microarray data. Gene Ontology enrichment and Kyoto Encyclopedia of Genes and Genomes pathway analysis were performed to determine the gene function. Furthermore, the lncRNA-mRNA coexpression network was generated to predict the potential targets of lncRNAs. In conclusion, our study first demonstrated the differential expression profiles of lncRNAs in the whole blood of infants with HIE and may provide a new view of the distinct lncRNA functions in HIE.     

Screening of significant biomarkers related with prognosis of liver cancer by lncRNA-associated ceRNAs analysis

This study aimed to identify significant biomarkers related to the prognosis of liver cancer using long noncoding RNA (lncRNA)-associated competing endogenous RNAs (ceRNAs) analysis. Differentially expressed mRNA and lncRNAs between liver cancer and paracancerous tissues were screened, and the functions of these mRNAs were predicted by gene ontology and pathway enrichment analyses. A ceRNA network consisting of differentially expressed mRNAs and lncRNAs was constructed. LncRNA FENDRR and lncRNA HAND2-AS1 were hub nodes in the ceRNA network. A risk score assessment model consisting of eight genes (PDE2A, ESR1, FBLN5, ALDH8A1, AKR1D1, EHHADH, ADRA1A, and GNE) associated with prognosis were developed. Multivariate Cox regression suggested that both pathologic_T and risk group could be regarded as independent prognostic factors. Furthermore, a nomogram model consisting of pathologic_T and risk group showed a good prediction ability for predicting the survival rate of liver cancer patients. The nomogram model consisting of pathologic_T and a risk score assessment model could be regarded as an independent factor for predicting prognosis of liver cancer.     

Intestinal epithelial cells related lncRNA and mRNA expression profiles in dextran sulphate sodium-induced colitis

Intestinal epithelial barrier damage caused by intestinal epithelial cells (IECs) dysfunction plays a crucial role in the pathogenesis and development of inflammatory bowel disease (IBD). Recently, some studies have suggested the emerging role of long non-coding RNAs (lncRNAs) in IBD. The aim of this study was to reveal lncRNAs and mRNA expression profiles in IECs from a mouse model of colitis and to expand our understanding in the intestinal epithelial barrier regulation. IECs from the colons of wild-type mice and dextran sulphate sodium (DSS)-induced mice were isolated for high-throughput RNA-sequencing. A total of 254 up-regulated and 1013 down-regulated mRNAs and 542 up-regulated and 766 down-regulated lncRNAs were detected in the DSS group compared with the Control group. Four mRNAs and six lncRNAs were validated by real-time quantitative PCR. Function analysis showed that dysregulated mRNAs participated in TLR7 signalling pathway, IL-1 receptor activity, BMP receptor binding and IL-17 signalling pathway. Furthermore, the possibility of indirect interactions between differentially expressed mRNAs and lncRNAs was illustrated by the competing endogenous RNA (ceRNA) network. LncRNA ENSMUST00000128026 was predicted to bind to mmu-miR-6899-3p, regulating Dnmbp expression. LncRNA NONMMUT143162.1 was predicted to competitively bind to mmu-miR-6899-3p, regulating Tnip3 expression. Finally, the protein-protein interaction (PPI) network analysis was constructed with 311 nodes and 563 edges. And the highest connectivity degrees were Mmp9, Fpr2 and Ccl3. These results provide novel insights into the functions of lncRNAs and mRNAs involved in the regulation of the intestinal epithelial barrier.     

Identification of mitochondrial function‐associated lncRNAs in septic mice myocardium

The present study aimed to analyze long noncoding RNA (lncRNA) and messenger RNA (mRNA) expression profiles in septic mice heart and to identify potential lncRNAs and mRNAs that be responsible for cardiac mitochondrial dysfunction during sepsis. Mice were treated with 10 mg/kg of lipopolysaccharides to induce sepsis. LncRNAs and mRNAs expression were evaluated by using lncRNA and mRNA microarray or real‐time polymerase chain reaction technique. LncRNA‐mRNA coexpression network assay, Gene Ontology (GO) analysis, Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis were performed. The results showed that 1275 lncRNAs were differentially expressed in septic myocardium compared with those in the control group. A total of 2769 mRNAs were dysregulated in septic mice heart, most of which are mainly related to the process of inflammation, mitochondrial metabolism, oxidative stress, and apoptosis. Coexpression network analysis showed that 14 lncRNAs were highly correlated with 11 mitochondria‐related differentially expressed mRNA. Among all lncRNAs and their cis‐acting mRNAs, 41 lncRNAs‐mRNA pairs (such as NONMMUG004378 and Apaf1 gene) were enriched in GO terms and KEGG pathways. In summary, we gained some specific lncRNAs and their potential target mRNAs that might be involved in mitochondrial dysfunction in septic myocardium. These findings provide a panoramic view of lncRNA and might allow developing new treatment strategies for sepsis.     

Different expression profile of mRNA and long noncoding RNA in autoimmune thyroid diseases patients

Increasing evidence has found that long noncoding RNAs (lncRNAs) and message RNAs (mRNAs) play an important role in the progress of autoimmune thyroid diseases (AITD). So, in this study, the different expressed of lncRNA and mRNA was screened by microarray analysis and quantitative real-time polymerase chain reaction (PCR). To further investigate the relationship among the differentially expressed genes, the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis and Gene ontology (GO) were used for biofunctions and signaling pathways analysis, respectively. Finally, the interaction relationship between lncRNA and mRNAs was analysis with lncRNA-mRNA coexpression network. The result found that the abnormal expression of lncRNAs and mRNAs were 1615 and 1913, respectively. The altered genes included CD40LG, IFNG, CTLA4, FAS, STAT1, STAT3, and STAT4. These were enriched in presentation and antigen processing, Th1 and Th2 cell differentiation, natural killer cell-mediated cytotoxicity, B cell receptor signaling pathway, Th17 cell differentiation, and cell adhesion molecules (CAMs), all of which had been suggested to be associated with immunopathogenic mechanisms and AITD-induced pathophysiologic changes. A coexpression network profile was contained with 126 network nodes and 477 connections which were based on seven mRNAs and 119 interacted lncRNAs. The outcomes of differentially expressed lncRNAs and their coreralated mRNAs in our study revealed that lncRNAs involved in immunopathogenic mechanisms may play a crital role in the pathogenesis of AITD.     

Alteration of LncRNA expression in mice placentae after frozen embryo transfer is associated with increased fetal weight

The birthweight after frozen embryo transfer (FET) was significantly higher compared with fresh embryo transfer (fresh ET), while the mechanism remains unclear. In this study, we transferred vitrified-warmed or fresh mice blastocysts into pseudopregnant recipients (n = 11 each group) produced by natural mating to avoid the influence of superovulation. The fetal weight, placental weight, placental efficiency and placental architecture were studied at E18.5. Placental RNA-Seq analysis was used to identify candidate different lncRNAs and mRNAs between the FET group and the fresh ET group. We found that the fetal weight was increased in the FET group, with increased placental efficiency and the proportion of placental function related labyrinth zone area. 554 lncRNAs and 1012 mRNAs were differentially expressed. KEGG and GO enrichment analyses showed these differentially expressed lncRNAs and their targeted mRNAs might be related to placental morphogenesis. Furthermore, the most differentially expressed 15 lncRNAs and 15 mRNAs were validated by qRT-PCR, we found the LncRNA embryonic stem cells expressed 1 (Lncenc1) was significantly decreased, and Gjb5, which played an important role in labyrinth zone development, was increased. Gjb5 protein increase was further confirmed by Western blot. Lncenc1 and Gjb5 had 48 predicted co-targeted miRNAs, while the correlation analysis of Lncenc1 and Gjb5 mRNA showed a significant inverse correlation. The results showed that FET treatment might enhance the placental function to increase mouse fetal weight via the network diagram of Lncenc1-miRNA-Gjb5.     

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.     

Characterization of Differentially Expressed Genes Involved in Pathways Associated with Gastric Cancer

To explore the patterns of gene expression in gastric cancer, a total of 26 paired gastric cancer and noncancerous tissues from patients were enrolled for gene expression microarray analyses. Limma methods were applied to analyze the data, and genes were considered to be significantly differentially expressed if the False Discovery Rate (FDR) value was < 0.01, P-value was <0.01 and the fold change (FC) was >2. Subsequently, Gene Ontology (GO) categories were used to analyze the main functions of the differentially expressed genes. According to the Kyoto Encyclopedia of Genes and Genomes (KEGG) database, we found pathways significantly associated with the differential genes. Gene-Act network and co-expression network were built respectively based on the relationships among the genes, proteins and compounds in the database. 2371 mRNAs and 350 lncRNAs considered as significantly differentially expressed genes were selected for the further analysis. The GO categories, pathway analyses and the Gene-Act network showed a consistent result that up-regulated genes were responsible for tumorigenesis, migration, angiogenesis and microenvironment formation, while down-regulated genes were involved in metabolism. These results of this study provide some novel findings on coding RNAs, lncRNAs, pathways and the co-expression network in gastric cancer which will be useful to guide further investigation and target therapy for this disease.     

EA15, MIR22, LINC00472 as diagnostic markers for diabetic kidney disease

This study aimed to investigate the molecular mechanisms of diabetic kidney disease (DKD) and to explore new potential therapeutic strategies and biomarkers for DKD. First we analyzed the differentially expressed changes between patients with DKD and the control group using the chip data in Gene Expression Omnibus (GEO) database. Then the gene chip was subjected to be annotated again, so as to screen long noncoding RNAs (lncRNAs) and study expression differences of these lncRNAs in DKD and controlled samples. At last, the function of the differential lncRNAs was analyzed. A total of 252 lncRNAs were identified, and 14 were differentially expressed. In addition, there were 1,629 differentially expressed messenger RNAs (mRNAs) genes, and proliferation and apoptosis adapter protein 15 (PEA15), MIR22, and long intergenic nonprotein coding RNA 472 ( LINC00472) were significantly differentially expressed in DKD samples. Through functional analysis of the encoding genes coexpressed by the three lncRNAs, we found these genes were mainly enriched in type 1 diabetes and autoimmune thyroid disease pathways, whereas in Gene Ontology (GO) function classification, they were also mainly enriched in the immune response, type I interferon signaling pathways, interferon-γ mediated signaling pathways, and so forth. To summary, we identified EA15, MIR22, and LINC00472 may serve as the potential diagnostic markers of DKD.     

Genome-wide differential expression of long noncoding RNAs and mRNAs in ovarian follicles of two different chicken breeds

Bashang long-tail chickens are an indigenous breed with dual purpose in China (meat and eggs) but have low egg laying performance. To improve the low egg laying performance, a genome-wide analysis of mRNAs and long noncoding RNAs (lncRNAs) from Bashang long-tail chickens and Hy-Line brown layers was performed. A total of 16,354 mRNAs and 8691 lncRNAs were obtained from ovarian follicles. Between the breeds, 160 mRNAs and 550 lncRNAs were found to be significantly differentially expressed. Integrated network analysis suggested some differentially expressed genes were involved in ovarian follicular development through oocyte meiosis, progesterone-mediated oocyte maturation, and cell cycle. The impact of lncRNAs on cis and trans target genes, indicating some lncRNAs may play important roles in ovarian follicular development. The current results provided a catalog of chicken ovarian follicular lncRNAs and genes for further study to understand their roles in regulation of egg laying performance.