Mesangial cell-derived CircRNAs in chronic glomerulonephritis: RNA sequencing and bioinformatics analysis

Abstract

Background Circular RNAs (circRNAs) have been shown to play critical roles in the initiation and progression of chronic glomerulonephritis (CGN), while their role from mesangial cells in contributing to the pathogenesis of CGN is rarely understood. Our study aims to explore the potential functions of mesangial cell-derived circRNAs using RNA sequencing (RNA-seq) and bioinformatics analysis. Methods Mouse mesangial cells (MMCs) were stimulated by lipopolysaccharide (LPS) to establish an in vitro model of CGN. Pro-inflammatory cytokines and cell cycle stages were detected by Enzyme-linked immunosorbent assay (ELISA) and Flow Cytometry experiment, respectively. Subsequently, differentially expressed circRNAs (DE-circRNAs) were identified by RNA-seq. GEO microarrays were used to identify differentially expressed mRNAs (DE-mRNAs) between CGN and healthy populations. Weighted co-expression network analysis (WGCNA) was utilized to explore clinically significant modules of CGN. CircRNA-associated CeRNA networks were constructed by bioinformatics analysis. The hub mRNAs from CeRNA network were identified using LASSO algorithms. Furthermore, utilizing protein–protein interaction (PPI), gene ontology (GO), pathway enrichment (KEGG), and GSEA analyses to explore the potential biological function of target genes from CeRNA network. In addition, we investigated the relationships between immune cells and hub mRNAs from CeRNA network using CIBERSORT. Results The expression of pro-inflammatory cytokines IL-1β, IL-6, and TNF-α was drastically increased in LPS-induced MMCs. The number of cells decreased significantly in the G1 phase but increased significantly in the S/G2 phase. A total of 6 DE-mRNAs were determined by RNA-seq, including 4 up-regulated circRNAs and 2 down-regulated circRNAs. WGCNA analysis identified 1747 DE-mRNAs of the turquoise module from CGN people in the GEO database. Then, the CeRNA networks, including 6 circRNAs, 38 miRNAs, and 80 mRNAs, were successfully constructed. The results of GO and KEGG analyses revealed that the target mRNAs were mainly enriched in immune, infection, and inflammation-related pathways. Furthermore, three hub mRNAs (BOC, MLST8, and HMGCS2) from the CeRNA network were screened using LASSO algorithms. GSEA analysis revealed that hub mRNAs were implicated in a great deal of immune system responses and inflammatory pathways, including IL-5 production, MAPK signaling pathway, and JAK-STAT signaling pathway. Moreover, according to an evaluation of immune infiltration, hub mRNAs have statistical correlations with neutrophils, plasma cells, monocytes, and follicular helper T cells. Conclusions Our findings provide fundamental and novel insights for further investigations into the role of mesangial cell-derived circRNAs in CGN pathogenesis.