scholarly journals Comprehensive analysis of circRNA expression profiles and circRNA-associated competing endogenous RNA networks in IgA nephropathy

PeerJ ◽  
2020 ◽  
Vol 8 ◽  
pp. e10395
Author(s):  
Haiyang Liu ◽  
Di Liu ◽  
Yexin Liu ◽  
Ming Xia ◽  
Yan Li ◽  
...  

Background Immunoglobulin A nephropathy (IgAN) is immune-mediated primary glomerulonephritis, which is the most common reason leading to renal failure worldwide. The exact pathogenesis of IgAN is not well defined. Accumulating evidence indicates that circular RNAs (circRNAs) play crucial roles in the immune disease by involving in the competing endogenous RNA (ceRNA) network mechanism. At present, the studies of the circRNA profiles and circRNA-associated ceRNA networks in the IgAN are still scarce. This study aimed to elucidate the potential roles of circRNA-associated ceRNA networks of peripheral blood mononuclear cells (PBMCs) in IgAN patients Method CircRNA sequencing was used to identify the differential expressed circRNAs (DEcircRNAs) of PBMCs in IgAN and healthy controls; limma packages from data sets GSE25590 and GSE73953 in the Gene Expression Omnibus (GEO) database, were used to identify differentially expressed micro RNAs (miRNAs) and message RNAs (mRNAs). A circRNA-miRNA-mRNA ceRNA network was constructed to further investigate the mechanisms of IgAN. Then, GO analysis and KEGG enrichment analyses were used to annotate the genes involved in the circRNA-associated ceRNA network. Further, Protein-protein interaction (PPI) networks were established to screen potential hub genes, by using Search Tool for the Retrieval of Interacting Genes/Proteins (STRING). Last, a quantitative real-time polymerase chain reaction (qRT-PCR) was applied to verify the hub genes in the ceRNA network. Result A total of 145 circRNAs, 22 miRNAs, and 1,117 mRNAs were differentially expressed in IgAN compared with controls (P < 0.05). A ceRNA network was constructed which contained 16 DEcircRNAs, 72 differential expressed mRNAs (DEmRNAs) and 11 differential expressed miRNAs (DEmiRNAs). KEGG pathway enrichment analysis illustrated the underlying biological functions of the ceRNA-associated genes, such as Nitrogen compound metabolic process, COPII-coated ER to Golgi transport vesicle, CAMP response element protein binding process (P < 0.01); meanwhile, Hepatitis B, GnRH signaling, and Prion disease were the most significant enrichment GO terms (P < 0.01). PPI network based on STRING analysis identified 4 potentially hub genes. Finally, Ankyrin repeat and SOCS box containing 16 (ASB16), SEC24 homolog C, COPII coat complex component (SEC24C) were confirmed by qRT-PCR (P < 0.05) and were identified as the hub genes of the ceRNA network in our study. Conclusion Our study identified a novel circRNA-mediated ceRNA regulatory network mechanisms in the pathogenesis of IgAN.

PeerJ ◽  
2019 ◽  
Vol 7 ◽  
pp. e7135 ◽  
Author(s):  
Fangcao Lei ◽  
Han Zhang ◽  
Xiaoli Xie

Background Pulpitis is a common inflammatory disease that affects dental pulp. It is important to understand the molecular signals of inflammation and repair associated with this process. Increasing evidence has revealed that long noncoding RNAs (lncRNAs), via competitively sponging microRNAs (miRNAs), can act as competing endogenous RNAs (ceRNAs) to regulate inflammation and reparative responses. The aim of this study was to elucidate the potential roles of lncRNA, miRNA and messenger RNA (mRNA) ceRNA networks in pulpitis tissues compared to normal control tissues. Methods The oligo and limma packages were used to identify differentially expressed lncRNAs and mRNAs (DElncRNAs and DEmRNAs, respectively) based on expression profiles in two datasets, GSE92681 and GSE77459, from the Gene Expression Omnibus (GEO) database. Differentially expressed genes (DEGs) were further analyzed by Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses. Protein–protein interaction (PPI) networks and modules were established to screen hub genes using the Search Tool for the Retrieval of Interacting Genes/Proteins (STRING) and the Molecular Complex Detection (MCODE) plugin for Cytoscape, respectively. Furthermore, an lncRNA-miRNA-mRNA-hub genes regulatory network was constructed to investigate mechanisms related to the progression and prognosis of pulpitis. Then, quantitative real-time polymerase chain reaction (qRT-PCR) was applied to verify critical lncRNAs that may significantly affect the pathogenesis in inflamed and normal human dental pulp. Results A total of 644 upregulated and 264 downregulated differentially expressed genes (DEGs) in pulpitis samples were identified from the GSE77459 dataset, while 8 up- and 19 downregulated probes associated with lncRNA were identified from the GSE92681 dataset. Protein–protein interaction (PPI) based on STRING analysis revealed a network of DEGs containing 4,929 edges and 623 nodes. Upon combined analysis of the constructed PPI network and the MCODE results, 10 hub genes, including IL6, IL8, PTPRC, IL1B, TLR2, ITGAM, CCL2, PIK3CG, ICAM1, and PIK3CD, were detected in the network. Next, a ceRNA regulatory relationship consisting of one lncRNA (PVT1), one miRNA (hsa-miR-455-5p) and two mRNAs (SOCS3 and PLXNC1) was established. Then, we constructed the network in which the regulatory relationship between ceRNA and hub genes was summarized. Finally, our qRT-PCR results confirmed significantly higher levels of PVT1 transcript in inflamed pulp than in normal pulp tissues (p = 0.03). Conclusion Our study identified a novel lncRNA-mediated ceRNA regulatory mechanisms in the pathogenesis of pulpitis.


2021 ◽  
Author(s):  
Li Guoquan ◽  
Du Junwei ◽  
He Qi ◽  
Fu Xinghao ◽  
Ji Feihong ◽  
...  

Abstract BackgroundHashimoto's thyroiditis (HT), also known as chronic lymphocytic thyroiditis, is a common autoimmune disease, which mainly occurs in women. The early manifestation was hyperthyroidism, however, hypothyroidism may occur if HT was not controlled for a long time. Numerous studies have shown that multiple factors, including genetic, environmental, and autoimmune factors, were involved in the pathogenesis of the disease, but the exact mechanisms were not yet clear. The aim of this study was to identify differentially expressed genes (DEGs) by comprehensive analysis and to provide specific insights into HT. MethodsTwo gene expression profiles (GSE6339, GSE138198) about HT were downloaded from the Gene Expression Omnibus (GEO) database. The DEGs were assessed between the HT and normal groups using the GEO2R. The DEGs were then sent to the Gene Ontology (GO) and the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis. The hub genes were discovered using Cytoscape and CytoHubba. Finally, NetworkAnalyst was utilized to create the hub genes' targeted microRNAs (miRNAs). ResultsA total of 62 DEGs were discovered, including 60 up-regulated and 2 down-regulated DEGs. The signaling pathways were mainly engaged in cytokine interaction and cytotoxicity, and the DEGs were mostly enriched in immunological and inflammatory responses. IL2RA, CXCL9, IL10RA, CCL3, CCL4, CCL2, STAT1, CD4, CSF1R, and ITGAX were chosen as hub genes based on the results of the protein-protein interaction (PPI) network and CytoHubba. Five miRNAs, including mir-24-3p, mir-223-3p, mir-155-5p, mir-34a-5p, mir-26b-5p, and mir-6499-3p, were suggested as likely important miRNAs in HT. ConclusionsThese hub genes, pathways and miRNAs contribute to a better understanding of the pathophysiology of HT and offer potential treatment options for HT.


2020 ◽  
Author(s):  
Yanjie Han ◽  
Xinxin Li ◽  
Jiliang Yan ◽  
Chunyan Ma ◽  
Xin Wang ◽  
...  

Abstract Background: Melanoma is the most deadly tumor in skin tumors and is prone to distant metastases. The incidence of melanoma has increased rapidly in the past few decades, and current trends indicate that this growth is continuing. This study was aimed to explore the molecular mechanisms of melanoma pathogenesis and discover underlying pathways and genes associated with melanoma.Methods: We used high-throughput expression data to study differential expression profiles of related genes in melanoma. The differentially expressed genes (DEGs) of melanoma in GSE15605, GSE46517, GSE7553 and the Cancer Genome Atlas (TCGA) datasets were analyzed. Differentially expressed genes (DEGs) were identified by paired t-test. Then the DEGs were performed cluster and principal component analyses and protein–protein interaction (PPI) network construction. After that, we analyzed the differential genes through bioinformatics and got hub genes. Finally, the expression of hub genes was confirmed in the TCGA databases and collected patient tissue samples.Results: Total 144 up-regulated DEGs and 16 down-regulated DEGs were identified. A total of 17 gene ontology analysis (GO) terms and 11 pathways were closely related to melanoma. Pathway of pathways in cancer was enriched in 8 DEGs, such as junction plakoglobin (JUP) and epidermal growth factor receptor (EGFR). In the PPI networks, 9 hub genes were obtained, such as loricrin (LOR), filaggrin (FLG), keratin 5 (KRT5), corneodesmosin (CDSN), desmoglein 1 (DSG1), desmoglein 3 (DSG3), keratin 1 (KRT1), involucrin (IVL) and EGFR. The pathway of pathways in cancer and its enriched DEGs may play important roles in the process of melanoma. The hub genes of DEGs may become promising melanoma candidate genes. Five key genes FLG, DSG1, DSG3, IVL and EGFR were identified in the TCGA database and melanoma tissues.Conclusions: The results suggested that FLG, DSG1, DSG3, IVL and EGFR might play important roles and potentially be valuable in the prognosis and treatment of melanoma.


2020 ◽  
Vol 2020 ◽  
pp. 1-12
Author(s):  
Xuefeng Gu ◽  
Dongyang Jiang ◽  
Yue Yang ◽  
Peng Zhang ◽  
Guoqing Wan ◽  
...  

Background. Moyamoya disease (MMD) is a rare cerebrovascular disease characterized by chronic progressive stenosis or occlusion of the bilateral internal carotid artery (ICA), the anterior cerebral artery (ACA), and the middle cerebral artery (MCA). MMD is secondary to the formation of an abnormal vascular network at the base of the skull. However, the etiology and pathogenesis of MMD remain poorly understood. Methods. A competing endogenous RNA (ceRNA) network was constructed by analyzing sample-matched messenger RNA (mRNA), long non-coding RNA (lncRNA), and microRNA (miRNA) expression profiles from MMD patients and control samples. Then, a protein-protein interaction (PPI) network was constructed to identify crucial genes associated with MMD. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes pathway (KEGG) enrichment analyses were employed with the DAVID database to investigate the underlying functions of differentially expressed mRNAs (DEmRNAs) involved in the ceRNA network. CMap was used to identify potential small drug molecules. Results. A total of 94 miRNAs, 3649 lncRNAs, and 2294 mRNAs were differentially expressed between MMD patients and control samples. A synergistic ceRNA lncRNA-miRNA-mRNA regulatory network was constructed. Core regulatory miRNAs (miR-107 and miR-423-5p) and key mRNAs (STAT5B, FOSL2, CEBPB, and CXCL16) involved in the ceRNA network were identified. GO and KEGG analyses indicated that the DEmRNAs were involved in the regulation of the immune system and inflammation in MMD. Finally, two potential small molecule drugs, CAY-10415 and indirubin, were identified by CMap as candidate drugs for treating MMD. Conclusions. The present study used bioinformatics analysis of candidate RNAs to identify a series of clearly altered miRNAs, lncRNAs, and mRNAs involved in MMD. Furthermore, a ceRNA lncRNA-miRNA-mRNA regulatory network was constructed, which provides insights into the novel molecular pathogenesis of MMD, thus giving promising clues for clinical therapy.


2021 ◽  
Author(s):  
Jiefang Zhou ◽  
Xiaowei Ji ◽  
Xiuwei Shen ◽  
Kefeng Yan ◽  
Peng Huang ◽  
...  

Abstract Objectives We identified functional genes and studied the underlying molecular mechanisms of diabetic cardiomyopathy (DCM) using bioinformatics tools. Methods Original gene expression profiles were obtained from the GSE21610 and GSE112556 datasets. We used GEO2R to screen the differentially expressed genes (DEGs). Gene Ontology and Kyoto Encyclopedia of Genes and Genomes pathway enrichment analyses were performed on DEGs. Protein–protein interaction (PPI) networks of DEGs were constructed using STRING and hub genes of signaling pathways were identified using Cytoscape. Aberrant hub gene expression was verified using The Cancer Genome Atlas dataset. Connectivity Map was used to predict the drugs that could treat DCM. Results The DEGs in DCM were mainly enriched in the nuclei and cytoplasm and involved in DCM- and chemokine-related signaling pathways. In the PPI network, 32 nodes were chosen as hub nodes and an RNA interaction network was constructed with 517 interactions. The expression of key genes (JPIK3R1, CCR9, XIST, WDFY3.AS2, hsa-miR-144-5p, and hsa-miR-146b-5p) was significantly different between DCM and normal tissues. Danazol, ikarugamycin, and semustine were identified as therapeutic agents against DCM using CMAP. Conclusion The identified hub genes could be associated with DCM pathogenesis and the above drugs could be used for treating DCM.


2021 ◽  
Author(s):  
jie Yang ◽  
Yan-Nan Tao ◽  
Fang-Xiao Hu ◽  
Yong-Zhi Chen ◽  
Xue-Song Yang ◽  
...  

Abstract Background: Increasing evidences uncover that lncRNAs play an important role in Isolated systolic hypertension (ISH). However, a systematic lncRNA-mRNA regulatory network is still absent in isolated systolic hypertension and atherosclerotic cerebral infarction patients (ISH & ACI).Aim:This research aims to establish a lncRNA-mRNA co-expression network in patients with ISH & ACI, to probe into the potential functions of lncRNA in those patients.Design and Setting:Expression profiles of lncRNA and mRNAs are collected and compared respectively from 8 patients with ISH and 8 patients with ISH & ACI by RNA-seq data.Methods: Differentially expressed lncRNAs and mRNAs were screened out via high-throughput sequencing in the plasma of ISH/ACI patients and control ISH patients. Then, a lncRNA-mRNA interaction network was built using the Pearson correlation coefficient by Cytoscape software. The expression levels of the hub genes and lncRNAs were verified by qRT-PCR in another 10 ISH/ACI patients and 10 control patients. Results: 2768 differentially expressed lncRNAs and 747 differentially expressed mRNAs were identified. 2 hub genes (CD226 and PARVB) and 11 lncRNAs were identified in the lncRNA-mRNA interaction network. qRT-PCR and cell assay results verified that lncRNAs ENST00000590604 and CD226 are highly expressed in patients of ISH & ACI. CD226 was associated with vascular endothelial cells growth and stability through platelet activation and focal adhesion pathway.Conclusion: We established a novel mRNA-lncRNA interaction network. lncRNAs ENST00000590604 and CD226 might be the potential biomarkers of ISH & ACI.


2019 ◽  
Vol 19 (1) ◽  
Author(s):  
Cheng Zhang ◽  
Bingye Zhang ◽  
Di Meng ◽  
Chunlin Ge

Abstract Background The incidence of cholangiocarcinoma (CCA) has risen in recent years, and it has become a significant health burden worldwide. However, the mechanisms underlying tumorigenesis and progression of this disease remain largely unknown. An increasing number of studies have demonstrated crucial biological functions of epigenetic modifications, especially DNA methylation, in CCA. The present study aimed to identify and analyze methylation-regulated differentially expressed genes (MeDEGs) involved in CCA tumorigenesis and progression by bioinformatics analysis. Methods The gene expression profiling dataset (GSE119336) and gene methylation profiling dataset (GSE38860) were obtained from the Gene Expression Omnibus (GEO) database. Differentially expressed genes (DEGs) and differentially methylated genes (DMGs) were identified using the limma packages of R and GEO2R, respectively. The MeDEGs were obtained by overlapping the DEGs and DMGs. Functional enrichment analyses of these genes were then carried out. Protein–protein interaction (PPI) networks were constructed using STRING and visualized in Cytoscape to determine hub genes. Finally, the results were verified based on The Cancer Genome Atlas (TCGA) database. Results We identified 98 hypermethylated, downregulated genes and 93 hypomethylated, upregulated genes after overlapping the DEGs and DMGs. These genes were mainly enriched in the biological processes of the cell cycle, nuclear division, xenobiotic metabolism, drug catabolism, and negative regulation of proteolysis. The top nine hub genes of the PPI network were F2, AHSG, RRM2, AURKB, CCNA2, TOP2A, BIRC5, PLK1, and ASPM. Moreover, the expression and methylation status of the hub genes were significantly altered in TCGA. Conclusions Our study identified novel methylation-regulated differentially expressed genes (MeDEGs) and explored their related pathways and functions in CCA, which may provide novel insights into a further understanding of methylation-mediated regulatory mechanisms in CCA.


2021 ◽  
Author(s):  
De-Bin Liu ◽  
You-Fu He ◽  
Gui-Jian Chen ◽  
Hua Huang ◽  
Xu-Ling Xie ◽  
...  

Abstract Background Aortic dissection (AD) is a rare and lethal disorder with its genetic basis remains largely unknown. Many studies have confirmed that circular RNAs (circRNAs) play important roles in various physiological and pathological processes. However, the roles of circRNAs in AD are still unclear and need further investigation. The present study aimed to elucidate the underlying molecular mechanisms of circRNAs regulation in aortic dissection based on the circRNA-miRNA-mRNA competing endogenous RNA (ceRNA) network. Methods Expression profiles of circRNAs (GSE97745), miRNAs (GSE92427), and mRNAs (GSE52093) were downloaded from Gene Expression Omnibus (GEO) databases, and the differentially expressed RNAs (DERNAs) were subsequently identified in AD by bioinformatics analysis. Further bioinformatics analyses, including circRNA-miRNA-mRNA ceRNA network, Gene Ontology (GO), and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses, were used to predict the potential functions of circRNA-associated ceRNA regulatory network. RNA was isolated from human arterial blood samples after which quantitative real-time PCR (qRT-PCR) was performed to confirm the DERNAs. Results We identified 14 (5 up-regulated and 9 down-regulated) differentially expressed circRNAs (DEcircRNAs), 17 (8 up-regulated and 9 down-regulated) differentially expressed miRNAs (DEmiRNAs) and 527 (297 up-regulated and 230 down-regulated) differentially expressed mRNAs (DEmRNAs) when AD samples were compared with normal ascending aorta samples (adjusted P-value < 0.05 and | log2FC |> 1.0). KEGG pathway analysis indicated that DEmRNAs were related to focal adhesion and extracellular matrix (ECM) receptor interaction signaling pathways. Simultaneously, the present study successfully constructed a ceRNA regulatory network based on 1 circRNAs (hsa_circRNA_082317), 1 miRNAs (hsa-miR-149-3p) and 10 mRNAs (MLEC, ENTPD7, SLC16A3, SLC7A8, TBC1D16, PAQR4, MAPK13, PIK3R2, ITGA5, SERPINA1) in AD. Furthermore, qRT-PCR demonstrated that hsa_circRNA_082317 andα5 integrin (ITGA5) were significantly up-regulated in AD (n = 3), and hsa-miR-149-3p was dramatically down-regulated in AD (n = 3). The expression of hsa-miR-149-3p target mRNA, ITGA5, was positively modulated by hsa_circRNA_082317. Conclusion This is the first study to demonstrate the circRNA-associated ceRNA regulatory network is altered in AD, implying that circRNAs may play important roles in regulating the onset and progression of AD and thus may serve as potential biomarkers for the diagnosis and treatment of AD.


2021 ◽  
Vol 16 (1) ◽  
Author(s):  
Zhanyu Yang ◽  
Delong Liu ◽  
Rui Guan ◽  
Xin Li ◽  
Yiwei Wang ◽  
...  

Abstract Background Heterotopic ossification (HO) represents pathological lesions that refer to the development of heterotopic bone in extraskeletal tissues around joints. This study investigates the genetic characteristics of bone marrow mesenchymal stem cells (BMSCs) from HO tissues and explores the potential pathways involved in this ailment. Methods Gene expression profiles (GSE94683) were obtained from the Gene Expression Omnibus (GEO), including 9 normal specimens and 7 HO specimens, and differentially expressed genes (DEGs) were identified. Then, protein–protein interaction (PPI) networks and Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses were performed for further analysis. Results In total, 275 DEGs were differentially expressed, of which 153 were upregulated and 122 were downregulated. In the biological process (BP) category, the majority of DEGs, including EFNB3, UNC5C, TMEFF2, PTH2, KIT, FGF13, and WISP3, were intensively enriched in aspects of cell signal transmission, including axon guidance, negative regulation of cell migration, peptidyl-tyrosine phosphorylation, and cell-cell signaling. Moreover, KEGG analysis indicated that the majority of DEGs, including EFNB3, UNC5C, FGF13, MAPK10, DDIT3, KIT, COL4A4, and DKK2, were primarily involved in the mitogen-activated protein kinase (MAPK) signaling pathway, Ras signaling pathway, phosphatidylinositol-3-kinase/protein kinase B (PI3K/Akt) signaling pathway, and Wnt signaling pathway. Ten hub genes were identified, including CX3CL1, CXCL1, ADAMTS3, ADAMTS16, ADAMTSL2, ADAMTSL3, ADAMTSL5, PENK, GPR18, and CALB2. Conclusions This study presented novel insight into the pathogenesis of HO. Ten hub genes and most of the DEGs intensively involved in enrichment analyses may be new candidate targets for the prevention and treatment of HO in the future.


2020 ◽  
Vol 9 (3) ◽  
pp. 90-98 ◽  
Author(s):  
Haitao Chen ◽  
Liaobin Chen

Aims This study aimed to uncover the hub long non-coding RNAs (lncRNAs) differentially expressed in osteoarthritis (OA) cartilage using an integrated analysis of the competing endogenous RNA (ceRNA) network and co-expression network. Methods Expression profiles data of ten OA and ten normal tissues of human knee cartilage were obtained from the Gene Expression Omnibus (GEO) database (GSE114007). The differentially expressed messenger RNAs (DEmRNAs) and lncRNAs (DElncRNAs) were identified using the edgeR package. We integrated human microRNA (miRNA)-lncRNA/mRNA interactions with DElncRNA/DEmRNA expression profiles to construct a ceRNA network. Likewise, lncRNA and mRNA expression profiles were used to build a co-expression network with the WGCNA package. Potential hub lncRNAs were identified based on an integrated analysis of the ceRNA network and co-expression network. StarBase and Multi Experiment Matrix databases were used to verify the lncRNAs. Results We detected 1,212 DEmRNAs and 49 DElncRNAs in OA and normal knee cartilage. A total of 75 dysregulated lncRNA-miRNA interactions and 711 dysregulated miRNA-mRNA interactions were obtained in the ceRNA network, including ten DElncRNAs, 69 miRNAs, and 72 DEmRNAs. Similarly, 1,330 dysregulated lncRNA-mRNA interactions were used to construct the co-expression network, which included ten lncRNAs and 407 mRNAs. We finally identified seven hub lncRNAs, named MIR210HG, HCP5, LINC00313, LINC00654, LINC00839, TBC1D3P1-DHX40P1, and ISM1-AS1. Subsequent enrichment analysis elucidated that these lncRNAs regulated extracellular matrix organization and enriched in osteoclast differentiation, the FoxO signalling pathway, and the tumour necrosis factor (TNF) signalling pathway in the development of OA. Conclusion The integrated analysis of the ceRNA network and co-expression network identified seven hub lncRNAs associated with OA. These lncRNAs may regulate extracellular matrix changes and chondrocyte homeostasis in OA progress. Cite this article: Bone Joint Res. 2020;9(3):90–98.


Sign in / Sign up

Export Citation Format

Share Document