Comparative Transcriptome Analysis Reveals Relationship among mRNAs, lncRNAs, and circRNAs of Slow Transit Constipation.

Slow transit constipation (STC) is characterized by persistent, infrequent, or incomplete defecation. Systematic analyses of mRNA, lncRNA, and circRNA expression profiling in STC provide insights to understand the molecular mechanisms of STC pathogenesis. The present study is aimed at observing the interaction of mRNAs, lncRNAs, and circRNAs by RNA sequencing in vivo of STC. A rat model of STC was induced by loperamide. The expression profiles of both mRNAs and miRNAs were performed by RNA sequencing. Enrichment analyses of anomalous expressed mRNAs, lncRNAs, and circRNAs were performed in order to identify the related biological functions and pathologic pathways through the Gene Ontology (GO) database and Kyoto Encyclopedia of Genes and Genomes (KEGG) database. In total, 26435 mRNAs, 5703 lncRNAs, and 7708 circRNAs differentially expressed were identified between the two groups. The analyses of GO and KEGG show that (1) upregulated genes were enriched in a positive regulation of GTPase activity, cell migration, and protein binding and lipid binding and (2) GO annotations revealed that most trans-target mRNAs are involved in the regulation process of immune signal together with the proliferation and differentiation of immune cells. Additionally, the protein-protein interaction (PPI) network of differentially expressed (DE) mRNAs was constructed. Interestingly, all of the core lncRNAs and their coexpression mRNAs in this network are downregulated. Moreover, downregulated circRNAs have a set of target mRNAs related to immunoreaction, which was consistent with the overall tendency. Our investigation enriches the STC transcriptome database and provides a preliminary exploration of novel candidate genes and avenues expression profiles in vivo. The dysregulation of mRNAs, lncRNAs, and circRNAs might contribute to the pathological processes during STC.

Yan S ，Yue Y ，Sun M ，Chen Y ，Wang X ，Qian H ... -  《-》

Analysis of differentially expressed lncRNAs and mRNAs associated with slow‑transit constipation.

Slow transit constipation (STC) is a refractory gastrointestinal disease, accounting for approximately 13 ∼ 37 % of chronic constipation. However, the molecular mechanism of STC remains poorly understood. Herein, this study aims to identify the key mRNAs and lncRNAs associated with STC. To this end, we performed high-throughput RNA sequencing to identify differentially expressed (DE) mRNAs and lncRNAs in the whole-layer sigmoid intestinal tissues from 4 STC patients and 4 non-STC patients. The identified DE lncRNAs and mRNAs were validated through quantitative real-time PCR. Weighted gene co-expression network analysis (WGCNA) and Pearson correlation analysis were conducted to determine the significantly correlated DE mRNA-lncRNA pairs. A total of 1420 DE lncRNAs and 1634 DE mRNAs were identified. Kyoto Encyclopedia of Genes and Genomes analysis of DE mRNAs indicated that these DE mRNAs might be associated with systemic lupus erythematosus, alcoholism, intestinal immune network for IgA production, inflammatory bowel disease, NF-kappa B signaling pathway. WGCNA and Pearson correlation analyses jointly identified 16,577 significantly correlated DE mRNA-lncRNA pairs. Furthermore, lncRNAs LINC00641, LINC02268, LINC03013 were identified as hub lncRNAs. The protein-protein interaction (PPI) network of proteins encoded by DE mRNAs was established, and PPI-based analysis revealed that Interleukin 2(IL2), CD80 molecule (CD80), interleukin-17A (IL-17A) might play significant roles in the development of STC. This study analyzes the expression profiles of lncRNAs and mRNAs associated with STC. Our findings will contribute to further understanding of the molecular mechanism of STC and provide potential diagnostic or therapeutic biomarkers for STC.

Miao Y ，Xie X ，Zhang Y ，Ma X ，Zhu X ，Li R ，Bi J ，Duan R ，Ai X ... -  《-》

Identification of potential biomarkers and pathways associated with carotid atherosclerotic plaques in type 2 diabetes mellitus: A transcriptomics study.

Type 2 diabetes mellitus (T2DM) affects the formation of carotid atherosclerotic plaques (CAPs) and patients are prone to plaque instability. It is crucial to clarify transcriptomics profiles and identify biomarkers related to the progression of T2DM complicated by CAPs. Ten human CAP samples were obtained, and whole transcriptome sequencing (RNA-seq) was performed. Samples were divided into two groups: diabetes mellitus (DM) versus non-DM groups and unstable versus stable groups. The Limma package in R was used to identify lncRNAs, circRNAs, and mRNAs. Gene Ontology (GO) annotation and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses, protein-protein interaction (PPI) network creation, and module generation were performed for differentially expressed mRNAs. Cytoscape was used to create a transcription factor (TF)-mRNA regulatory network, lncRNA/circRNA-mRNA co-expression network, and a competitive endogenous RNA (ceRNA) network. The GSE118481 dataset and RT-qPCR were used to verify potential mRNAs.The regulatory network was constructed based on the verified core genes and the relationships were extracted from the above network. In total, 180 differentially expressed lncRNAs, 343 circRNAs, and 1092 mRNAs were identified in the DM versus non-DM group; 240 differentially expressed lncRNAs, 390 circRNAs, and 677 mRNAs were identified in the unstable versus stable group. Five circRNAs, 14 lncRNAs, and 171 mRNAs that were common among all four groups changed in the same direction. GO/KEGG functional enrichment analysis showed that 171 mRNAs were mainly related to biological processes, such as immune responses, inflammatory responses, and cell adhesion. Five circRNAs, 14 lncRNAs, 46 miRNAs, and 54 mRNAs in the ceRNA network formed a regulatory relationship. C22orf34-hsa-miR-6785-5p-RAB37, hsacirc_013887-hsa-miR-6785-5p/hsa-miR-4763-5p/hsa-miR-30b-3p-RAB37, MIR4435-1HG-hsa-miR-30b-3p-RAB37, and GAS5-hsa-miR-30b-3p-RAB37 may be potential RNA regulatory pathways. Seven upregulated mRNAs were verified using the GSE118481 dataset and RT-qPCR. The regulatory network included seven mRNAs, five circRNAs, six lncRNAs, and 14 TFs. We propose five circRNAs (hsacirc_028744, hsacirc_037219, hsacirc_006308, hsacirc_013887, and hsacirc_045622), six lncRNAs (EPB41L4A-AS1, LINC00969, GAS5, MIR4435-1HG, MIR503HG, and SNHG16), and seven mRNAs (RAB37, CCR7, CD3D, TRAT1, VWF, ICAM2, and TMEM244) as potential biomarkers related to the progression of T2DM complicated with CAP. The constructed ceRNA network has important implications for potential RNA regulatory pathways.

Yu T ，Xu B ，Bao M ，Gao Y ，Zhang Q ，Zhang X ，Liu R ... -  《Frontiers in Endocrinology》

Identification of Key lncRNAs, circRNAs, and mRNAs in Osteoarthritis via Bioinformatics Analysis.

Osteoarthritis (OA) is a common degenerative joint disorder that adversely affects the quality of life of patients. Identification of novel diagnostic biomarkers is pivotal for the early detection and prevention of OA. Dataset GSE185059 was selected from Gene Expression Omnibus database to obtain differentially expressed lncRNAs (DE-lncRNAs), mRNAs (DE-mRNAs), and circRNAs (DE-circRNAs) between OA and normal samples. The Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses as well as protein-protein interaction (PPI) network construction of DE-mRNAs were conducted. Hub genes were identified from PPI networks and validated by RT-qPCR. starBase database was utilized for predicting miRNAs binding with hub genes, selected DE-lncRNAs and DE-circRNAs, respectively. The competing endogenous RNA (ceRNA) networks were constructed. A total of 818 DE-mRNAs, 191 DE-lncRNAs, and 2053 DE-circRNAs were identified. The DE-mRNAs were significantly enriched in several inflammation-related GO terms and KEGG pathways such as positive regulation of cell-cell adhesion, TNF-alpha signaling pathway and NF-kappa B signaling pathway. Thirteen hub genes were identified, which were CFTR, GART, SMAD2, NCK1, TJP1, UBE2D1, EFTUD2, PRKACB, IL10, SNRPG, CHD4, RPS24, and SRSF6. OA-related DE-lncRNA/circRNA-miRNA-hub gene networks were constructed. We identified 13 hub genes and constructed the ceRNA networks related to OA, providing a theoretical basis for further research.

Zhang W ，Wei C ，Wang L 《-》

Construction and functional enrichment analysis of the competitive endogenous RNA regulatory network for nonarteritic anterior ischemic optic neuropathy based on high-throughput sequencing.

Based on the transcriptome high-throughput sequencing of nonarteritic anterior ischemic optic neuropathy (NAION), this study constructed a competitive endogenous RNA (ceRNA) network, enriched and analyzed it, screened long noncoding (lnc)RNAs and circular (circ)RNAs that may participate in the competitive endogenous mechanism in NAION, and inferred its function. Four milliliters of peripheral blood from NAION patients and the control group was extracted from clinical samples, transcriptome high-throughput sequencing was performed, and the sequencing data were visualized. Based on the principle of the ceRNA network, the lncRNA-micro (mi)RNA-messenger (m)RNA interaction axis and circRNA-miRNA‒mRNA interaction axes were constructed. The differentially expressed genes (DEGs) in the interaction axis were enriched and analyzed by Gene Ontology (GO) and the Kyoto Encyclopedia of Genes and Genomes (KEGG), and the functions and signalling pathways of lncRNAs and circRNAs in the interaction network were speculated. Fifty-one circRNAs were differentially expressed in the sequencing data: 25 were upregulated, and 26 were downregulated. For 996 differentially expressed lncRNAs, 317 were upregulated and 679 were downregulated, and for 1161 differentially expressed mRNAs, 698 were upregulated and 463 were downregulated. Thirty-three differentially expressed miRNAs, upregulated miRNA 18 and downregulated miRNA 15 were identified. After screening, 13 coexpressed mRNAs, 15 lncRNAs, and 3 miRNAs were finally constructed in the lncRNA-miRNA-mRNA network, and the circRNA-miRNA-mRNA network was constructed by 159 mRNAs, 26 miRNAs, and 34 circRNAs. In the lncRNA network, GO enrichment analysis obtained 182 biological processes, 12 cell components and 38 molecular functions, which are related mainly to the regulation of the activity of proteins and enzymes such as cyclic nucleotide-dependent protein kinase activity and magnesium ion-dependent protein serine/threonine phosphatase activity. KEGG analysis involved mainly the forkhead box O (FoxO) signalling pathway, apelin signalling pathway, etc. In the circRNA enrichment results, 353 biological processes, 52 cell components, and 45 molecular functions were obtained, involving mainly calcium channel regulation, neutrophil activation, mRNA transport, and metabolism. KEGG mainly involved the Wnt signalling pathway, apelin signalling pathway, Hippo signalling pathway, oxytocin signalling pathway, etc. This paper provides a new idea for lncRNAs and circRNAs to mediate the occurrence and development of NAION through the ceRNA mechanism. This study lays a foundation for the in-depth study of the pathogenesis of NAION.

Zhang M ，Wang X 《-》