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.     

急性髓系白血病mRNA与非编码RNA差异表达谱及CeRNA调控网络分析

利用GEO数据库(gene expression omnibus database)通过生物信息学分析方法探讨急性髓系白血病(acute myelogenous leukemia,AML)的发病机制。检索GEO数据库中AML相关芯片数据集GSE142698、GSE142699和GSE96535。利用GEO2R分析得到差异mRNAs、miRNAs以及差异lncRNAs。利用在线生物信息学分析工具DAVID对差异mRNAs进行GO富集分析和KEGG通路分析。利用miRWalk数据库预测AML相关miRNAs的靶向mRNAs,利用Spongescan数据库预测AML相关miRNAs的靶向lncRNAs,构建lncRNA-miRNA-mRNA竞争性内源RNA （competing endogenous RNA,ceRNA）调控网络。共筛选出29个显著差异mRNAs、70个显著差异miRNAs和20 005个显著差异lncRNAs。GO富集分析和KEGG通路分析显示,差异表达基因主要涉及蛋白磷酸化、细胞分裂、细胞增殖的负调控、基因表达的正向调节、周期蛋白依赖的丝氨酸/苏氨酸激酶活性的调节等生物过程以及细胞周期、细胞衰老、癌症通路、PI3K-Akt通路等信号通路。将miRWalk数据库预测的靶向mRNAs与差异mRNAs取交集,Spongescan数据库预测的靶向lncRNAs与差异lncRNAs取交集,分别确定了25个mRNAs、6个lncRNAs参与AML相关ceRNA调控网络的构建。结果表明,lncRNAs可能作为关键的ceRNA,通过调控miRNA和相关靶基因参与AML的发生与发展,研究结果为AML诊断和治疗的分子生物学研究提供了新的依据。     

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.     

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.     

Transcriptional profiling of exosomes derived from plasma of canine with mammary tumor by RNA-seq analysis

Canine mammary tumor (CMT) are the second most common tumor in dogs. Exosomes can act as biomarkers for the early diagnosis of tumors, and also be involved in the pathogenesis and metastasis mechanism of tumors. The expression profile of exosomal RNA revealed that there were a total of 5547 differentially expressed mRNAs, and 196 differentially expressed lncRNAs. GO and KEGG enrichment analysis found that the differentially expressed mRNAs and lncRNA target genes were associated with metabolic processes, DNA replication, cell proliferation, cell junction, and cell adhesion. In conclusion, this study revealed lncRNA and mRNA expression profiles in exosomes derived from plasma of CMT and further annotated their potential functions. The data obtained in this study will also provide valuable resources for understanding lncRNA information in plasma exosomes of dogs with CMT, and contribute to the study of early diagnostic markers and pathogenesis of CMT.     

Aberrantly expressed long non‐coding RNAs in air pollution‐induced congenital defects

Air pollution has been a serious public health issue over the past few decades particularly in developing countries. Air pollution exposure during pregnancy poses potential threat to offspring as the deleterious substances might pass through placenta to alter foetal development. A growing number of studies have demonstrated that long non‐coding RNAs (lncRNAs) participate in the development of many diseases, including congenital defects. Here, we used RNA sequencing to identify differentially expressed lncRNAs in air pollution‐exposed rat embryos compared with control group. Our data suggested that 554 lncRNAs (216 up‐regulated and 338 down‐regulated) were significantly differentially expressed in the air pollution‐exposed embryos. Moreover, potential cellular functions of these deregulated lncRNAs were predicted via KEGG signal pathway/GO enrichment analyses, which suggested the possible involvements of neurological process, sensory perception of smell and the G‐protein signalling pathway. Furthermore, potential functional network of deregulated lncRNAs and their correlated mRNAs in the development of congenital spinal abnormality was established. Our data suggested that lncRNAs may play a vital role in the pathophysiology of air pollution‐exposed congenital spinal malformation.     

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.     

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.     

LncRNA‐mRNA expression profiles and functional networks in osteoclast differentiation

Human osteoclasts are differentiated from CD14⁺ monocytes and are responsible for bone resorption. Long non‐coding RNAs (lncRNAs) have been proved to be significantly involved in multiple biologic processes, especially in cell differentiation. However, the effect of lncRNAs in osteoclast differentiation is less appreciated. In our study, RNA sequencing (RNA‐seq) was used to identify the expression profiles of lncRNAs and mRNAs in osteoclast differentiation. The results demonstrated that expressions of 1117 lncRNAs and 296 mRNAs were significantly altered after osteoclast differentiation. qRT‐PCR assays were performed to confirm the expression profiles, and the results were almost consistent with the RNA‐seq data. GO and KEGG analyses were used to predict the functions of these differentially expressed mRNA and lncRNAs. The Path‐net analysis demonstrated that MAPK pathway, PI3K‐AKT pathway and NF‐kappa B pathway played important roles in osteoclast differentiation. Co‐expression networks and competing endogenous RNA networks indicated that ENSG00000257764.2‐miR‐106a‐5p‐TIMP2 may play a central role in osteoclast differentiation. Our study provides a foundation to further understand the role and underlying mechanism of lncRNAs in osteoclast differentiation, in which many of them could be potential targets for bone metabolic disease.     

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.     

Construction of lncRNA-associated ceRNA networks to identify prognostic lncRNA biomarkers for glioblastoma

Long noncoding RNAs (lncRNAs) serve as competitive endogenous RNAs (ceRNAs) that play significant regulatory roles in the pathogenesis of tumors. However, the role of lncRNAs, especially the lncRNA-related ceRNA regulatory network, in glioblastoma (GBM) has not been fully elucidated. The goal of the current study was to construct lncRNA-microRNA-mRNA-related ceRNA networks for further investigation of their mechanism of action in GBM. We downloaded data from The Cancer Genome Atlas (TCGA) and the Gene Expression Omnibus (GEO) databases and identified differential lncRNAs, microRNAs (miRNAs), and messenger RNAs (mRNAs) associated with GBM. A ceRNA network was constructed and analyzed to examine the relationship between lncRNAs and patients’ overall survival. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGGs) were used to analyze the related mRNAs to indirectly explain the mechanism of action of lncRNAs. The potential effective drugs for the treatment of GBM were identified using the connectivity map (CMap). After integrated analysis, we obtained a total of 210 differentially expressed lncRNAs, 90 differentially expressed miRNAs, and 2508 differentially expressed mRNAs (DEmRNAs) from the TCGA and GEO databases. Using these differential genes, we constructed a lncRNA-associated ceRNA network. Six lncRNAs in the ceRNA network were associated with the overall survival of patients with GBM. Through KEGG analysis, it was found that the DEmRNAs involved in the network are related to cancer-associated pathways, for instance, mitogen-activated protein kinase and Ras signaling pathways. CMap analysis revealed four small-molecule compounds that could be used as drugs for the treatment of GBM. In this study, a multi-database joint analysis was used to construct a lncRNA-related ceRNA network to help identify the regulatory functions of lncRNAs in the pathogenesis of GBM.     

Screening and bioinformatics analysis of mRNA,long non-coding RNA and circular RNA expression profiles in mucoepidermoid carcinoma of salivary gland

Mucoepidermoid carcinoma (MEC) of salivary gland is a disease characterized by high rate of diatant metastasis, and associated with poor outcomes. However, the molecular mechanisms underlying the MEC remain poorly understand. Here, we simultaneously detected, for the first time, the expression profiles of mRNAs, lncRNAs, and circRNAs in four pairs of MEC and matched non-carcinoma tissues by microarrays. A total of 3612?mRNA, 3091 lncRNAs, and 284 circRNAs were altered during the pathogenesis of MEC. The functions of these differentially expressed RNAs were predicted by Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis. Co-expression networks of lncRNA-mRNA and circRNA-miRNA were conducted to uncovered the hidden ceRNA mechanisms. Moreover, NONHSAT154433.1 that associated with ADAM12 and hsa_circ_0012342 were further screened and confirmed using qRT-PCR analysis. In conclusion, this study provides a systematic perspective on the potential function of non-coding RNAs (ncRNAs) in the molecular mechanisms of MEC. Among these, NONHSAT154433.1 and hsa_circ_0012342 might be served as potential prognostic biomarkers and therapeutic target of MEC.     

小鼠正常黑素细胞与黑色素瘤细胞（B16）mRNA表达差异分析

黑色素瘤是一种高侵袭性的恶性皮肤肿瘤,转移率高、预后差。研究黑素瘤细胞生物学特性对黑素瘤的治疗和控制具有重要的意义。本研究以C57BL/6J小鼠的正常黑色素细胞及B16黑色素瘤细胞为研究对象,采用二代测序技术分析两种细胞间的转录组表达差异,筛选差异基因,为后续黑色素瘤的形成机制研究提供理论依据。采用差异倍数及错误率分析测序数据,鉴定出1 436个新的mRNA和4 086个差异表达的已知mRNA。GO数据库和KEGG数据库分析显示,差异表达的mRNAs参与了149个调控途径, 主要集中在疾病调控、细胞周期调节和环境信息调控方面。qRT-PCR及Western印迹检测发现,调节细胞增殖、迁移的Pdgf-B、Integrin-β1和Integrin-β5以及调节黑色素颗粒增加的Mitf、Tyr、Tyrp1和Tyrp2在B16细胞中的表达量显著高于在正常黑色素细胞中的表达。本研究获得的差异基因为后续黑色素瘤的研究提供了新的候选基因。