Time-Series Clustering of lncRNA-mRNA Expression during the Adipogenic Transdifferentiation of Porcine Skeletal Muscle Satellite Cells

Skeletal muscle satellite cells (SMSCs), which are multifunctional muscle-derived stem cells, can differentiate into adipocytes. Long-chain non-coding RNA (lncRNA) has diverse biological functions, including the regulation of gene expression, chromosome silencing, and nuclear transport. However, the regulatory roles and mechanism of lncRNA during adipogenic transdifferentiation in muscle cells have not been thoroughly investigated. Here, porcine SMSCs were isolated, cultured, and induced for adipogenic differentiation. The expressions of lncRNA and mRNA at different time points during transdifferentiation were analysed using RNA-seq analysis. In total, 1005 lncRNAs and 7671 mRNAs showed significant changes in expression at differential differentiation stages. Time-series expression analysis showed that the differentially expressed (DE) lncRNAs and mRNAs were clustered into 5 and 11 different profiles with different changes, respectively. GO, KEGG, and REACTOME enrichment analyses revealed that DE mRNAs with increased expressions during the trans-differentiation were mainly enriched in the pathways for lipid metabolism and fat cell differentiation. The genes with decreased expressions were mainly enriched in the regulation of cell cycle and genetic information processing. In addition, 1883 DE mRNAs were regulated by 193 DE lncRNAs, and these genes were related to the controlling in cell cycle mainly. Notably, three genes in the fatty acid binding protein (FABP) family significantly and continuously increased during trans-differentiation, and 15, 13, and 11 lncRNAs may target FABP3, FABP4, and FABP5 genes by cis- or trans-regulation, respectively. In conclusion, these studies identify a set of new potential regulator for adipogenesis and cell fate and help us in better understanding the molecular mechanisms of trans-differentiation.     

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.     

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.     

Comprehensive analysis of lncRNA–miRNA–mRNA during proliferative phase of rat liver regeneration

This study aims to reveal the regulatory mechanism of lncRNAs–miRNAs–mRNAs network during the proliferative phase of liver regeneration (LR). High-throughput sequencing technology was performed, and a total of 1,738 differentially expressed lncRNAs (DE lncRNAs), 167 known differentially expressed miRNAs (DE miRNAs), and 2,727 differentially expressed mRNAs were identified. Then, the target DE lncRNAs and DE mRNAs regulated by the same miRNAs were screened and a ceRNA regulatory network containing 32 miRNAs, 107 lncRNAs, and 270 mRNAs was constructed. Insulin signaling pathway, pyrimidine metabolism, axon guidance, carbohydrate digestion and absorption, and pyruvate metabolism were significantly enriched in the network. Through literature review and the regulatory relationship between lncRNAs and miRNAs, nine core lncRNAs were identified, which might play important roles during the proliferative phase of rat LR. This study analyzed lncRNA–miRNA–mRNA regulatory network for the first time during the proliferative phase of rat LR, providing clues for exploring the mechanism of LR and the treatment of liver diseases.     

Tissue-specific analysis of Coffea arabica L. transcriptome revealed potential regulatory roles of lncRNAs

Long non-coding RNAs (lncRNAs) play pivot roles in regulating mRNA expression in eukaryotic organisms without coding any proteins. In the current study, a comprehensive analysis of 260 published RNA-Seq datasets collected from different tissues (fruits, leaves, stems, and roots) of Coffea arabica L. was performed to discover potential lncRNAs. A total of 10,564 unique lncRNAs were identified. Our results showed that 77.14% of the lncRNAs were intergenic and 60.39% of them are located within 5 Kbp from the partner gene. In general, all the identified lncRNAs showed shorter lengths, fewer number of exons, and lower expression levels as compared to mRNAs in different studied tissues. Several lncRNAs were determined as differentially expressed (DE) in fruits as compared to leaves, stems, or roots. The functional characterization of the DE lncRNAs revealed their roles in regulating significant biological processes in different tissues of C. arabica. The current study provides a comprehensive analysis and dataset of lncRNAs in C. arabica that could be utilized in further studies concerning the roles of these molecules in plant cells.     

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.     

Differential expression profile of mRNAs,lncRNAs, and circRNAs reveals potential molecular mechanism in breast cancer

In recent years, breast cancer attracts more and more attention because of its high incidence. To explore the molecular functions and mechanisms, we performed RNA sequencing on the tumor tissues and their paired normal tissues from three breast cancer patients. By differential expression analysis, we found 3764 differentially expressed (DE) mRNAs, 5416 DE lncRNAs, and 148 DE circRNAs. Enrichment analysis suggested that the DE lncRNAs and DE circRNAs were enriched in mitochondria and nucleus, which indicated that they may participate in the vital metabolism directly or indirectly, such as fatty acid metabolism. Subsequently, the protein–protein interaction (PPI) network was constructed and we got 8 key proteins, of which the matrix metalloproteinase-9 (MMP9; degree 5) draws our attention. Based on the 38 up-regulated circRNAs and 14 down-regulated circRNAs, we constructed competing endogenous RNA (ceRNA) networks, from which the has-miR-6794-5p has been identified to enriched in the up-regulated network and correlated with the circNFIX directly. At this point, we presented that the circNFIX and MMP9 may play a significant role by regulating fatty acid metabolism in breast cancer.     

Construction and analysis of the abnormal lncRNA–miRNA–mRNA network in hypoxic pulmonary hypertension

The incidence of hypoxic pulmonary hypertension (HPH) is increasing. Accumulating evidence suggests that long noncoding RNAs (lncRNAs) play an important role in HPH, but the functions and mechanism have yet to be fully elucidated. In the present study, we established a HPH rat model with 8 h of hypoxia exposure (10% O₂) per day for 21 days. High-throughput sequencing identified 60 differentially expressed (DE) lncRNAs, 20 DE miRNAs and 695 DE mRNAs in rat lung tissue. qRT-PCR verified the accuracy of the results. The DE mRNAs were significantly enriched in immune response, inflammatory response, leukocyte migration, cell cycle, cellular response to interleukin-1, IL-17 signalling pathway, cytokine–cytokine receptor interaction and Toll-like receptor signalling pathway. According to the theory of competing endogenous RNA (ceRNA) networks, lncRNA–miRNA–mRNA network was constructed by Cytoscape software, 16 miRNAs and 144 mRNAs. The results suggested that seven DE lncRNAs (Ly6l, AABR07038849.2, AABR07069008.2, AABR07064873.1, AABR07001382.1, AABR07068161.1 and AABR07060341.2) may serve as molecular sponges of the corresponding miRNAs and play a major role in HPH.     

The local immune response in the mammary gland of the sow following infusion of a protein antigen

A study was made of the local immune response in the udder of the sow following infusion of a soluble antigen. Four mammary glands of each of four pregnant sows were infused with ferritin prepartum. Samples of blood, colostrum, and milk were collected during the following lactation; animals were slaughtered and mammary tissue removed for immunohistology. Blood, colostrum, milk, and mammary tissues were similarly collected from nonimmunized (control) sows. Colostral and milk whey from immunized sows contained higher levels of immunoglobulins than whey from control sows. There was an increase in numbers of IgA-containing plasma cells and total lymphoid cells in mammary tissue of immunized sows compared with controls. The results suggested that the local immune response was at least as great in non-infused glands as infused glands of immunized sows.     

Analyses of changes in myocardial long non-coding RNA and mRNA profiles after severe hemorrhagic shock and resuscitation via RNA sequencing in a rat model

Background

Ischemia–reperfusion injury has been proven to induce organ dysfunction and death, although the mechanism is not fully understood. Long non-coding RNAs (lncRNAs) have drawn wide attention with their important roles in the gene expression of some biological processes and diseases, including myocardial ischemia–reperfusion (I/R) injury. In this paper, a total of 26 Sprague–Dawley (SD) rats were randomized into two groups: sham and ischemia–reperfusion (I/R) injury. Hemorrhagic shock was induced by removing 45% of the estimated total blood volume followed by reinfusion of shed blood. High-throughput RNA sequencing was used to analyze differentially expressed (DE) lncRNAs and messenger RNAs (mRNAs) in the heart tissue 4 h after reperfusion. Myocardial function was also evaluated.

Results

After resuscitation, the decline of myocardial function of shocked animals, expressed by cardiac output, ejection fraction, and myocardial performance index (MPI), was significant (p?<?0.05). DE lncRNAs and mRNAs were identified by absolute value of fold change?≥?2 and the false discovery rate ≤?0.001. In rats from the I/R injury group, 851 lncRNAs and 1015 mRNAs were significantly up-regulated while 1533 lncRNAs and 1702 m RNAs were significantly down-regulated when compared to the sham group. Among the DE lncRNAs, we found 12 location-associated with some known apoptosis-related protein-coding genes which were up-regulated or down-regulated accordingly, including STAT3 and Il1r1. Real time PCR assays confirmed that the expression levels of five location-associated lncRNAs (NONRATT006032.2, NONRATT006033.2, NONRATT006034.2, NONRATT006035.2 and NONRATT029969.2) and their location-associated mRNAs (STAT3 and Il1r1) in the rats from the I/R injury group were all significantly up-regulated versus the sham group.

Conclusions

The DE lncRNAs (NONRATT006032.2, NONRATT006033.2, NONRATT006034.2 and NONRATT006035.2) could be compatible with their role in myocardial protection by stimulating their co-located gene (STAT3) after hemorrhagic shock and resuscitation. The final prognosis of I/R injury might be regulated by different genes, which is regarded as a complex network.

     

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.