Identification of Serum Exosome-Derived circRNA-miRNA-TF-mRNA Regulatory Network in Postmenopausal Osteoporosis Using Bioinformatics Analysis and Validation in Peripheral Blood-Derived Mononuclear Cells.
Osteoporosis is one of the most common systemic metabolic bone diseases, especially in postmenopausal women. Circular RNA (circRNA) has been implicated in various human diseases. However, the potential role of circRNAs in postmenopausal osteoporosis (PMOP) remains largely unknown. The study aims to identify potential biomarkers and further understand the mechanism of PMOP by constructing a circRNA-associated ceRNA network.
The PMOP-related datasets GSE161361, GSE64433, and GSE56116 were downloaded from the Gene Expression Omnibus (GEO) database and were used to obtain differentially expressed genes (DEGs). Gene ontology (GO) enrichment analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis were applied to determine possible relevant functions of differentially expressed messenger RNAs (mRNAs). The TRRUST database was used to predict differential transcription factor (TF)-mRNA regulatory pairs. Afterwards, combined CircBank and miRTarBase, circRNA-miRNA as well as miRNA-TF pairs were constructed. Then, a circRNA-miRNA-TF-mRNA network was established. Next, the correlation of mRNAs, TFs, and PMOP was verified by the Comparative Toxicogenomics Database. And expression levels of key genes, including circRNAs, miRNAs, TFs, and mRNAs in the ceRNA network were further validated by quantitative real-time PCR (qRT-PCR). Furthermore, to screen out signaling pathways related to key mRNAs of the ceRNA network, Gene Set Enrichment Analysis (GSEA) was performed.
A total of 1201 DE mRNAs, 44 DE miRNAs, and 1613 DE circRNAs associated with PMOP were obtained. GO function annotation showed DE mRNAs were mainly related to inflammatory responses. KEGG analysis revealed DE mRNAs were mainly enriched in osteoclast differentiation, rheumatoid arthritis, hematopoietic cell lineage, and cytokine-cytokine receptor interaction pathways. We first identified 26 TFs and their target mRNAs. Combining DE miRNAs, miRNA-TF/mRNA pairs were obtained. Combining DE circRNAs, we constructed the ceRNA network contained 6 circRNAs, 4 miRNAs, 4 TFs, and 12 mRNAs. The expression levels of most genes detected by qRT-PCR were generally consistent with the microarray results. Combined with the qRT-PCR validation results, we eventually identified the ceRNA network that contained 4 circRNAs, 3 miRNAs, 3 TFs, and 9 mRNAs. The GSEA revealed that 9 mRNAs participate in many important signaling pathways, such as "olfactory transduction", "T cell receptor signaling pathway", and "neuroactive ligand-receptor interaction". These pathways have been reported to the occurrence and development of PMOP. To sum up, key mRNAs in the ceRNA network may participate in the development of osteoporosis by regulating related signal pathways.
A circRNA-associated ceRNA network containing TFs was established for PMOP. The study may help further explore the molecular mechanisms and may serve as potential biomarkers or therapeutic targets for PMOP.
Dong Q
,Han Z
,Tian L
《Frontiers in Endocrinology》
Identification of circular RNA-associated competing endogenous RNA network in the development of cleft palate.
Circular RNAs (circRNAs) serve as competing endogenous RNAs (ceRNAs) and indirectly regulate gene expression through shared microRNAs (miRNAs). However, the regulatory mechanisms of circRNA as ceRNA associated with the fusion of palatal shelves in palatogenesis are yet unclear. This study aimed to explore the potential mechanism underlying the role of circRNA as ceRNA in cleft palate (CP). First, we systematically analyzed RNA-seq and miRNA-seq data after high-throughput sequencing for embryonic palatal shelf tissues from a mouse CP model induced by maternal exposure to all-trans retinoic acid on embryonic gestation day 14.5 (E14.5). Thirty-nine circRNAs, 18 miRNAs, and 936 messenger RNAs (mRNAs) were significantly dysregulated (log2 [fold change {FC}] > 1; P < 0.05). Thereafter, we constructed a circRNA-associated ceRNA network. Finally, we determined the circRNA_0954-miRNA-881-3p-PRKAR1α ceRNA network as a hub involved in palatogenesis. Gene Ontology analysis revealed that ceRNA-related genes were associated with facial morphogenesis and developmental gene silencing. Kyoto Encyclopedia of Genes and Genomes pathway analysis indicated that ceRNA-related genes are involved in apoptosis (P < 0.05, fold enrichment >1). Quantitative reverse transcription polymerase chain reaction was performed to verify the results of ceRNA analysis. We found that the circRNA-miRNA-mRNA ceRNA network is involved in palatogenesis. The present results imply that circRNA_0954-miRNA-881-3p-PRKAR1α ceRNA network may cause dysfunctional palatal fusion and might facilitate the development of novel epigenetic biomarkers to treat CP in the future.
Shu X
,Cheng L
,Dong Z
,Shu S
... -
《-》
Identification of a Plasma Exosomal lncRNA- and circRNA-Based ceRNA Regulatory Network in Patients With Lung Adenocarcinoma.
Exosomes have been established to be enriched with various long noncoding RNAs (lncRNAs) and circular RNAs (circRNAs) that exert various biological effects. However, the lncRNA- and circRNA-mediated coexpression competing endogenous RNA (ceRNA) regulatory network in exosomes derived from the plasma of patients with lung adenocarcinoma (LUAD) remains elusive.
This study enrolled nine patients with lung adenocarcinoma and three healthy individuals, and the differential expression of messenger RNAs (mRNAs), lncRNAs, and circRNAs was detected using microarray analysis, while microRNAs (miRNAs) were detected through RNA sequencing. Additionally, bioinformatics algorithms were applied to evaluate the lncRNA-miRNA-mRNAs/circRNA-miRNA-mRNA network. Differentially expressed cicRNAs were identified via quantitative reverse transcription polymerase chain reaction (RT-qPCR). A total of 1016 lncRNAs, 1396 circRNAs, 45 miRNAs, and 699 mRNAs were differentially expressed in the plasma exosomes of patients with LUAD compared with healthy controls. Among them, 881 lncRNAs were upregulated and 135 were downregulated, 916 circRNAs were upregulated while 480 were downregulated, 45 miRNAs were upregulated while none were downregulated, and 591 mRNAs were upregulated while 108 were downregulated (p ≤ 0.05, and fold change ≥ 2). Gene Ontology (GO) analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis revealed the biological functions of differentially expressed RNAs. Meanwhile, the RNA networks displayed the regulatory relationship between dysregulated RNAs. Finally, RT-qPCR validated that the expression of circ-0033861, circ-0043273, and circ-0011959 was upregulated in the plasma exosome of patients with LUAD compared to healthy controls (p = 0.0327, p = 0.0002, p = 0.0437, respectively).
This study proposed a newly discovered ncRNA-miRNA-mRNA/circRNA-miRNA-mRNA ceRNA network and identified that the expression of circulating circ-0033861, circ-0043273, and circ-0011959 was up-regulated in the plasma exosomes of patients with LUAD, offering valuable insights for exploring the potential function of exosomal noncoding RNA and identifying potential biomarkers for LUAD.
Zhu W
,Zhang H
,Tang L
,Fang K
,Lin N
,Huang Y
,Zhang Y
,Le H
... -
《-》