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.     

A landscape of Long non-coding RNAs reveals the leading transcriptome alterations in murine aorta during aging

Considerable studies have given convincing evidence of a forefront position for vascular aging in preventing cardiovascular disease. Various functions of Long non-coding RNAs (lncRNAs) are becoming increasingly distinct in aging-related diseases. This study aims at a better insight into the expression profile and mechanisms of lncRNAs in vascular senescence. High-throughput sequencing was used to detect the differential expression (DE) of lncRNAs and mRNAs in the aorta of 96 W and 8 W-old mice, while 1423 lncRNAs and 80 mRNAs were differentially expressed. By performing GO and KEGG enrichment analysis, we found that DE lncRNAs were mainly involved in purine metabolism and cGMP-PKG signaling pathways. In addition, a co-expression functional network of DE lncRNAs and DE mRNAs was constructed, and ENSMUST00000218874 could interact with 41 DE mRNAs, suggesting that it may play an essential role in vascular senescence. This study reveals DE lncRNAs in naturally aging vascular, which may provide new ideas and targets for aging-related cardiovascular diseases.     

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.     

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.     

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 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.     

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.     

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.     

Changing expression profiles of mRNA,lncRNA, circRNA,and miRNA in lung tissue reveal the pathophysiological of bronchopulmonary dysplasia (BPD) in mouse model

New perinatal care technologies have improved the survival rate of preterm neonates, but the prevalence of bronchopulmonary dysplasia (BPD), one of the most intractable problems in neonatal intensive care unit (NICU), remains unchanged. In present study, high-throughput sequencing (HTS) was performed to detect the expression profiles of long noncoding RNAs (lncRNAs), messenger RNAs (mRNAs), circular RNAs (circRNAs), and microRNAs (miRNAs) in hyperoxia-induced BPD mouse model. Significant differentially expressed RNAs were selected and clustered between the BPD group and the control group. The results revealed that expressions of 1778 lncRNAs, 1240 mRNAs, 97 circRNAs, and 201 miRNAs were significantly altered in the BPD group. Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) were performed to predict the potential functions of differentially expressed RNAs. lncRNA-mRNA and circRNA-miRNA coexpression networks were constructed to detect their association with the pathogenesis of BPD. Our study provides a systematic perspective on the potential function of RNAs during BPD.     

Expression profiling of long noncoding RNAs associated with vasculogenic mimicry in osteosarcoma

Osteosarcoma (OS) is the most common highly malignant bone tumor in teens. Vasculogenic mimicry (VM) is defined as de novo extracellular matrix-rich vascular-like networks formed by highly aggressive tumor cells. We previously reported the presence of VM and it is an unfavorable prognostic factor in OS patients. Long noncoding RNAs (lncRNAs) are aberrantly expressed in OS and involved in cancer cell VM. However, lncRNAs in VM formation of OS have not been investigated. We, therefore, profiled the expression of lncRNAs in highly aggressive OS cell line 143B compared with its parental poorly aggressive cell line HOS. The differentially expressed (DE) lncRNAs and messenger RNA (mRNAs) were subjected to constructed lncRNA-mRNA coexpressed network. The top-ranked hub gene lncRNA n340532 knockdown 143B cells were used for in vitro and in vivo VM assays. The annotation of DE lncRNAs was performed according to the coexpressed mRNAs by Gene Ontology and pathway analysis. A total of 1360 DE lncRNAs and 1353 DE mRNAs were screened out. lncRNA MALAT1 and FTX, which have known functions related to VM formation and tumorigenesis were identified in our data. The coexpression network composed of 226 lncRNAs and 118 mRNAs in which lncRNA n340532 had the highest degree number. lncRNA n340532 knockdown reduced VM formation in vitro. The suppression of n340532 also exhibited potent anti-VM and antimetastasis effect in vivo, suggesting its potential role in OS VM and metastasis. Furthermore, n340532 coexpressed with 10 upregulation mRNAs and 3 downregulation mRNAs. The enriched transforming growth factor-β signaling pathway, angiogenesis and so forth were targeted by those coexpressed mRNAs, implying n340532 may facilitate VM formation in OS through these pathways and gene functions. Our findings provide evidence for the potential role of lncRNAs in VM formation of OS that could be used in the clinic for anti-VM therapy in OS.     

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.