Skip to main content

Comprehensive analysis of IgA nephropathy expression profiles: identification of potential biomarkers and therapeutic agents

Abstract

Background

IgA nephropathy (IgAN) is a kidney disease recognized by the presence of IgA antibody depositions in kidneys. The underlying mechanisms of this complicated disease are remained to be explored and still, there is an urgent need for the discovery of noninvasive biomarkers for its diagnosis. In this investigation, an integrative approach was applied to mRNA and miRNA expression profiles in PBMCs to discover a gene signature and novel potential targets/biomarkers in IgAN.

Methods

Datasets were selected from gene expression omnibus database. After quality control checking, two datasets were analyzed by Limma to identify differentially expressed genes/miRNAs (DEGs and DEmiRs). Following identification of DEmiR-target genes and data integration, intersecting mRNAs were subjected to different bioinformatic analyses. The intersecting mRNAs, DEmiRs, related transcription factors (from TRRUST database), and long-non coding RNAs (from LncTarD database) were used for the construction of a multilayer regulatory network via Cytoscape.

Result

“GSE25590” (miRNA) and “GSE73953” (mRNA) datasets were analyzed and after integration, 628 intersecting mRNAs were identified. The mRNAs were mainly associated with “Innate immune system”, “Apoptosis”, as well as “NGF signaling” pathways. A multilayer regulatory network was constructed and several hub-DEGs (Tp53, STAT3, Jun, etc.), DEmiRs (miR-124, let-7b, etc.), TFs (NF-kB, etc.), and lncRNAs (HOTAIR, etc.) were introduced as potential factors in the pathogenesis of IgAN.

Conclusion

Integration of two different expression datasets and construction of a multilayer regulatory network not only provided a deeper insight into the pathogenesis of IgAN, but also introduced several key molecules as potential therapeutic target/non-invasive biomarkers.

Peer Review reports

Highlights

  • Innate immune system, apoptosis, and NGF signaling were the most linked pathways to the pathogenesis of IgAN.

  • A multilayer regulatory network was constructed using 3 regulatory elements: long non-coding RNAs, differentially expressed microRNAs and transcription factors.

  • Tp53, STAT3, and Jun were among top DEGs in the constructed multilayer regulatory network.

  • miR-124 was introduced as a non-invasive biomarker for detection of IgAN disease.

Background

IgA nephropathy (IgAN) is known as the most common primary glomerular disorder and one of the main causes of kidney failure worldwide [1,2,3]. Up to 40% of IgAN patients will finally progress to end-stage renal disease (ESRD), a condition with huge stress for the patients and a massive economic burden for the governments [4]. Besides hematuria and proteinuria, renal biopsy for checking the presence of mesangial IgA kidney deposits is still the gold standard for IgAN diagnosis [5]. However, due to the invasive nature of this procedure, its application is not favorable [6, 7]. Despite a huge number of investigations on the IgAN pathogenesis, the underlying mechanisms of this complicated disease are yet to be fully determined [8, 9]. Therefore, discovery of the disease-related pathways and key regulatory agents with a therapeutic/biomarker potential is of utmost necessity to not only shed a light on the disease pathogenicity, but also provide a tool for a non-invasive diagnosis/efficient treatment of IgAN [10].

Today, a wave of researches is dealing with clarifying the role of genes in IgAN pathogenicity. So far, some genes and miRNAs have been identified as key factors in the progression of IgAN [9, 11, 12]. As a matter of fact, due to the complex nature of IgAN, one can assume that no single gene or circuit, but a remarkable number of genes and regulatory networks are involved in the pathogenesis of this disease. Expression profiling with high-throughput techniques has become an extensively applied technology to identify disease associated genes/miRNAs and to identify novel biomarkers in complex diseases like IgAN [10, 13, 14]. In order to translate the produced raw data, they need to be analyzed via various available bioinformatics tools. Moreover, to take the benefits of bigger sample size, researchers could integrate two or more similar datasets in order to achieve more reliable results.

In the IgAN pathogenicity, kidney is thought to be an innocent bystander and the PBMCs as the places for IgA1 production and post-translational modification, are the primary constituents in seeking for the underlying mechanisms of the disease. Accordingly, a huge number of experiments have shown the principal role of PBMCs in IgAN pathogenicity [15, 16]. Moreover, due to the lack of a precise and none-invasive strategy for diagnosis of this disease, identifying potential biomarkers in blood is of great interest [17].

The objective of this investigation was to re-analyze and integrate the existing expression profiles coming from peripheral blood mononuclear cells (PBMCs) of IgAN patients. Upon integration of datasets and identifying intersecting mRNAs and different enrichment analyses, a multilayer regulatory network comprising of transcription factors (TFs), miRNAs, and long-non coding RNAs (lncRNAs) will be constructed to catch a holistic view over the interactions and involvement of hub molecules in the PBMCs of IgAN patients. In this study, we tried to identify a specific molecular signature, novel therapeutic targets, potential biomarkers, and core involved pathways in the pathogenesis of IgAN. Different steps of this experiment are demonstrated in Fig. 1.

Fig. 1
figure1

The flowchart representing different steps of the study

Methods

Dataset selection

Gene Expression Omnibus (GEO) database, as a freely accessible data archive, disseminating functional genomic datasets (https://www.ncbi.nlm.nih.gov/gds) was used for searching and achieving datasets [18]. “IgAN” AND “blood” AND/OR “miRNA” keywords were used in searching for related datasets. Principle component analysis (PCA) using R software (Version 1.2.5033), was used for sample quality control and clustering.

Identification of differentially expressed miRNAs (DEmiRs) and differentially expressed genes (DEGs)

Identification of DEGs and DEmiRs were performed using NetworkAnalyst website (http://www.networkanalyst.ca). After performing normalization procedures (variance stabilizing normalization followed by quantile normalization), Limma was utilized for identification of DEmiRs and DEGs in the selected datasets. False discovery rate (FDR) < 0.05 cutoff was applied to select the most significant DEmiRs and DEGs. Volcano plots for the DEGs and DEmiRs were built using the R software.

Identification of intersecting mRNAs

In this step, at first, the DEmiR-target genes were identified using Mienturnet online platform (http://userver.bio.uniroma1.it/apps/,mienturnet/). In this web tool, MirTarBase option was checked to achieve the most reliable validated DEmiR-target genes. After that, a web-based Venn-diagram tool (http://bioinformatics.psb.ugent.be/webtools/Venn/), was applied to identify the intersecting mRNAs between the DEmiR-target genes and the DEGs from the mRNA dataset.

Enrichment analysis

Gene ontology analysis (molecular function, biological process, and cellular component) was performed using EnrichR, a web-based tool for gene set enrichment analysis)https://amp.pharm.mssm.edu/Enrichr/). Pathway analysis was performed using the same tool based on Reactome database. FunRich (http://www.funrich.org/) was used to represent the enrichment analysis results [19].

Construction of a multilayer regulatory network and identification of hub nodes

The STRING (version 11) with selecting the highest interaction score (highest confidence: 0.900) was used for the construction of an interactive network among DEGs. TRRUST (version 2.0) database was utilized for the prediction of the DEGs-related transcript factors (TFs) [20, 21].

DEGs and DEmiRs related LncRNAs were selected among all the curated lncRNAs in the LncTarD (www.biocc.hrbmu.edu.cn/LncTarD). After merging and construction of a multilayer network using CytoScape (vesrion: 3.2.0) [22], topmost DEGs, DEmiRs, TFs, and lncRNAs were identified separately using cytoHubba plugin [17] in CytoScape.

Results

Dataset selection, extraction, and quality control

Based on our considerations in searching for transcription profiles of IgAN blood samples (PBMCs), two mRNA (GSE58539 and GSE73953) and one miRNA (GSE25590) datasets, were obtained from GEO. The mRNA dataset “GSE73953” included 15 samples from patients with IgAN and two pooled control samples from 16 healthy individuals. The used microarray platform was GPL4133 Agilent-014850 Whole Human Genome Microarray 4x44K G4112F. The miRNA dataset “GSE25590” included 7 samples from healthy individuals and 7 samples from IgAN patients. The samples were assessed using GPL7731 Agilent-019118 Human miRNA Microarray 2.0 G4470B platform. To explore the similarity of the samples in each dataset, sample-level quality control was performed using Principal Component Analysis (PCA). PCA, which is known as a standard technique aiming to reduce the dataset dimensionality, can also be a good tool for representation of the dataset quality [23]. In a so-called good quality dataset, the case and control samples are clustered, separately. Here, due to an inappropriate sample clustering pattern, “GSE58539” dataset was removed from the study and two “GSE73953” and “GSE25590” datasets were selected for further analysis.

Another main advantage of PCA clustering is to find possible sample outliers. Removing the sample outliers will give a better sample clustering (case vs. control) also more efficient data analysis. In the case of miRNA data set “GSE25590”, after PCA analysis and removing one healthy sample (as the outlier sample), more significant DEmiRs were identified in this dataset. The other dataset (mRNA dataset; GSE73953) had no outlier and the ‘healthy vs. patients’ samples were separated on PC1 in the PCA plot (Fig. 2a and b).

Fig. 2
figure2

PCA clustering pattern of the healthy and IgAN patient’s samples in mRNA “GSE73953” (a) and miRNA “GSE25590” (b) datasets. Volcano plots representing the DEGs in the mRNA (c), and DEmiRs in the miRNA datasets (d)

623 Intersecting mRNAs were identified after analysis of mRNA and miRNA datasets

All DEGs and DEmiRs were identified after sample grouping, normalization and analyzing the expression profiles using limma. Considering the FDR, 1390 significantly dysregulated genes (664 down-regulated and 729 up-regulated DEGs), were selected from the analyzed mRNA dataset “GSE73953” (Fig. 2c). In case of DEmiRs, 204 significantly dysregulated miRNAs (149 up-regulated and 55 down-regulated DEmiRs), were identified (Fig. 2d). The topmost differentially expressed mRNAs and miRNAs according to log2FC and adjusted p-value are listed in Table 1.

Table 1 List of top 20 up- and down-regulated DEGs and DEmiRs in the analyzed datasets

After identification of 7116 DEmiR targets for all the up- and down-regulated DEmiRs, and considering DEGs from the mRNA dataset, 628 intersecting mRNAs were selected as the most reliable differentially expressed mRNAs in the PBMCs of IgAN patients (Fig. 3a). These intersecting mRNAs, then were subjected to further enrichment analyses and network constructions.

Fig. 3
figure3

a: Venn diagram representing the intersecting mRNAs among the identified DEGs and DEmiR-gene targets. b: Topmost enriched biological process, c: molecular function, d: cellular component GO terms and e: Reactome pathways for the intersecting mRNAs

Enrichment analysis revealed the involvement of different pathways in the pathogenesis of IgAN

The topmost enriched pathways and GO terms for all the 628 intersecting mRNAs are shown in Fig. 3b-e. According to the GO enrichment analyses, top enriched biological process GO terms included “Positive regulation of transcription from RNA polymerase II”, “Signal transduction”, “Apoptosis process” and “Protein phosphorylation” (Fig. 3b). Also, “RNA binding”, “Metal ion binding”, “Protein kinase binding” and “Transcription factor binding” were among the most enriched molecular function GO terms (Fig. 3c) for the selected mRNAs. Likewise, the hub DEGs were mainly enriched in “Cytosol” and “Nucleus” and “Membrane” compartments (Fig. 3d). According to the Reactome pathway enrichment analysis, intersecting mRNAs were mainly associated with the “Innate immune system”, “Apoptosis”, “Toll-receptor cascades”, “NGF and EGFR signaling”, “Interferon signaling”, “Notch signaling” as well as “MyD88 dependent cascade initiated on endosome” (Fig. 3e).

The constructed multilayer regulatory network highlighted the key regulatory elements

To evaluate three levels of regulatory factors affecting the expression of the DEGs (intersecting mRNAs), a multilayer regulatory network including DEGs interactive network and their related TFs, DEmiRs and lncRNAs was constructed and analyzed (Fig. 4a). In a stepwise manner, after the construction of a regulatory network using DEGs and DEmiRs, their related TFs and lncRNAs were identified separately and inserted into the network as other regulatory layers.

Fig. 4
figure4

Multilayer regulatory network (a), and its derived sub-networks showing top regulatory molecules and their targets (b-d). a: Multilayer regulatory network comprising of the intersecting mRNAs (DEGs), DEmiRs, and their related TFs and lncRNAs. The constituents of this multilayer network are including 46 transcription factors (green circle), 132 DEmiRs (yellow circle), 217 LncRNAs (magenta circle), and 598 DEGs (blue circle). The DEmiRs and DEGs with no interactions are omitted from the network. b: Top miRNA molecule (miR-124) (in case of degree centrality and log2FC in the DEmiR dataset), c: top lncRNA (HOTAIR), (D): top transcription factor (NF-κB) and their interactions. *DEGs, DEmiRs, lncRNAs and TFs are colored in blue, yellow, violet and green, respectively. The multilayer regulatory network is accessible at network data exchange (NDEx) server by the below link: [https://public.ndexbio.org/#/network/a1d9669e-7daf-11eb-9e72-0ac135e8bacf?accesskey=c229bf52538f314edd2617ce1f6c6e6cd81ebcbef1437c6d7c81167c3590d973]

In this multilayer network which comprised of 993 nodes and 4407 edges, the topmost nodes from each layer were identified according to degree centrality, independently. Among all DEGs, the top ones, including 30 DEGs with high degree scores in the network were identified as hub-DEGs. Likewise, 3 TFs, 6 miRNAs and 10 lncRNAs were identified based on their scores (degree score) in the multilayer network (Table 2). The selection of hub molecules was based on their degree centrality in the network and 5% of molecules in each group with highest score was introduced as hub molecules.

Table 2 List of hub-DEGs, and topmost DEmiRs, TFs and lncRNAs in the constructed multilayer regulatory network

Topmost DEmiRs in the network were comprising of miR-16, miR-92a, let-7a, miR-17, miR-124, and mir-93a. miR-124-3p with log2FC: 4.7 in the DEmiR dataset (a member of the top 10 DEmiRs based on the fold change), could be considered as a non-invasive biomarker for detection of IgAN (Fig. 4b). Likewise, let-7b as another hub-DEmiR was identified as a top DEmiR according to the adjusted p-value in the DEmiR dataset.

Due to the importance of gene expression regulation at the transcription level, we were interested to find other regulatory elements like lncRNAs and TFs which could affect the expression of the identified DEGs. LncRNAs are shown to play a remarkable role in transcriptional and post-transcriptional regulations. These molecules also can act as decoys for miRNA binding and form complex regulatory networks that regulate the abundance of miRNA molecules [24]. Among all the curated lncRNAs in the LncTarD database, 217 lncRNA was identified to have interactions with the DEmiRs, DEGs, and the predicted TFs in the regulatory network. In the case of centrality score, HOTAIR was identified as the top lncRNA interacting with 9 DEmiRs, 10 DEGs and 10 TFs in the regulatory network (Fig. 4c). MEG3 and H19 were other top ranked lncRNAs interacting with 24 and 21 DEmiRs/DEGs/TFs in the constructed regulatory network. TFs as other core elements in regulatory machinery also were considered for construction of the multilayer regulatory network. Among all the predicted TFs, NF-κB (nuclear factor-kappa B1) was identified as the top TF interacting with 39 DEGs in the network (Fig. 4d). Other TFs like MYC and EZH2 also were identified as top TFs interacting with both DEGs and lncRNAs.

Discussion

Re-analyzing and integration of omics-based expression profiles could be an advantageous tactic to catch a holistic view of the involving genes and miRNAs in the IgAN pathogenesis. Until now, several re-analyzes have been performed on the IgAN related expression profiles, mostly on kidney tissue samples of the patients [25,26,27,28]. As far as we know, there is just one experiment conducting a meta-analysis study seeking DEGs in the PBMCs of IgAN patients. However, the mentioned experiment is suffering from a lack of dataset quality control checking and not considering other regulatory elements like miRNAs, or lncRNAs in introducing hub-DEGs [29]. In the present study, following quality control checking, we performed re-analysis and integration of two expression datasets (mRNA and miRNA) related to PBMCs of IgAN patients. By applying a systems biology approach, we identified several key genes, TFs, miRNAs, and their related lncRNAs which may have a potential role in the pathogenesis of IgAN. According to the results of pathway enrichment analysis, the intersecting mRNAs were mainly enriched in apoptosis, innate immunity, and NGF signaling pathways.

P53 and XIAP (X-linked inhibitor of apoptosis protein), were among the top DEGs with the highest degree centralities in the regulatory network. These two proteins are considered as the main participants in controlling the apoptosis and proliferation of cells. However, the activity of these two regulatory constituents is unlike; While p53 activity is in the favor of apoptosis and preventing the cell proliferation, XIAP action is to inhibit caspase 3, 7, and 9 and stop the apoptosis [30,31,32]. Considering the well-known roles of p53 and XIAP in controlling apoptosis [33], these proteins might take a part in controlling the PBMCs proliferation in IgAN patients. Besides, similar to the previous results [34], our findings also pointed to the apoptotic phenotype of PBMCs in IgAN patients. Moreover, based on some evidence, P53 can suppress autoimmunity (e.g., systemic lupus erythematosus) by suppressing T cell activity while the precise mechanisms are still unknown [35].

Other than apoptosis, “innate immunity”, “toll-like receptor cascades” and “interferon signaling” were other enriched pathways for the identified DEGs in the PBMCs of IgAN patients. Such enrichments were predictable due to the autoimmune nature of IgAN disease. Similar to previous findings, the involvement of these pathways could be an indication for hyper-activation and enhancement of antigen processing pathways in PBMCs of IgAN patients [36, 37]. Dysregulation of innate immunity has also been comprehensively suggested in IgAN patients [38].

Toll-like receptors (TLRs) are known as the main components of the innate immune system. The upregulation and higher activity of these receptors in IgAN samples have been shown by previous experiments [39, 40]. Following ligand binding, TLRs trigger several cascades that finally result in the activation of immune system. Through the production of inflammatory cytokines, these receptors have been shown to induce glomerular damages in IgAN patients [41]. Moreover, interaction of some TLRs (e.g., TLR1, TLR2) could result in the activation of Nuclear actor Kappa B Subunit 1 (NF-κB), as well as various immune cells like keratinocytes, dendritic cells, mast cells, B cells, and NK cells [42]. Therefore, these receptors could be considered as potential targets of more investigations regarding the IgAN pathogenesis.

Nerve growth factor (NGF) signaling was another main enriched pathway for the DEGs in PBMCs of IgAN patients. NGF as a prototypical example of neurotrophic factors is mainly involved in the growth, proliferation and differentiation of neurons [43]. However, the role of NGF in the proliferation of immune cells, stimulation of IgM, IgA and IgG production in lymphocytes, as well as regulation of inflammatory mediators like IL-1 beta, TNF-alpha, and IL-6 have been shown previously [44, 45]. Moreover, some findings have revealed the mediatory role of NGF in autoimmune diseases like multiple sclerosis, rheumatoid arthritis, and systemic lupus erythematosus [43]. Other findings also disclosed an increased serum level of NGF, as well as higher NGF receptor expression in PBMCs of glomerulonephritis patients [46]. The findings of the present study also could be a shred of evidence for the contribution of NGF signaling in the pathological process of IgAN disease.

Generally, the involvement of all above-mentioned pathways along with other enriched pathways like notch signaling and EGFR pathways could point to an aberrant proliferation rate and distortion of PBMCs in the IgAN patients. Involvement of the enriched pathways may provide new explanations about the underlying immunopathological mechanisms of IgAN.

In the next step, a multilayer regulatory network was constructed to identify different hub molecules as the main potential participants in the pathogenesis of the disease. The concept of constructing a multilayer network was to achieve a holistic aspect of regulatory interactions in the PBMCs of IgAN patients. Similar experiments in this area of research only consider DEGs to construct an interactive regulatory network, while skipping other regulatory factors like TFs, miRNAs, and lncRNAs. Such studies suffer from a narrow view and ignore the complex regulatory nature of the cell. Accordingly, a multilayer regulatory network is thought to be more realistic and closer to the complex nature of cellular regulatory machinery. The key regulatory constituents in the multilayer network are listed in Table 2. Here, we introduced top DEGs with the highest degree scores in the network as potential biomarker/therapeutic targets. Moreover, in case of TFs, NF-κB was identified the top TFs in regulating the expression of DEGs in IgAN patients. The trace of NF-κB in IgAN pathogenicity has been previously shown, where toll-like receptors trigger a cascade of intracellular messages that finally leads to the activation of this transcription regulator [47]. The consequence of such activation is the cellular release of different cytokines and chemokines, which may participate in aberrant galactosylation of IgA1 [48].

MiRNAs as other main regulatory elements also could play a major role in immunological and pathological aspects of IgAN. Along with their regulatory roles, higher stability in the biological milieu makes these small RNAs as potent non-invasive biomarkers for the early detection of diseases. In case of DEmiRs in the constructed regulatory network, most of the top miRNAs including miR-16, miR-92, miR-17, and let-7b have been shown to play a direct role in controlling the cell cycle and cell proliferation [49,50,51,52]. miR-124-3p as another hub-DEmiR in the regulatory network was also among the top 10 DEmiRs in the DEmiR dataset (log2FC: 4.7). According to previous experiments, miR-124 has a crucial role in the innate immunity system and fine-tuning the toll-like receptor responses via targeting signal transducer and activator of transcription 3 (STAT3) [47]. The results of our analysis also pointed to the negative regulatory role of miR-124 in inflammatory responses and introduced this miRNA as a therapeutic target, as well as a non-invasive biomarker for detection of IgAN disease.

Back to the multilayer regulatory network and talking about regulatory roles of miRNAs, let-7b was introduced as another top DEmiR able to target several DEGs and thought to play a pivotal role in the IgAN process/pathogenesis. Up-regulation of this miRNA in PBMCs of IgAN patients and its role in the glycosylation process of IgA with regulation of N-acetylgalactosaminyltransferase 2 (GALNT2) enzyme was previously shown [53, 54]. Overexpression of let-7b in PBMCs of IgAN patients might be a reason for the reduced expression levels of GALNT2 enzyme (log2FC: − 4.43 in the DEG dataset), and consequently aberrant galactosylation of IgA1 in patients. As a result, based on this analysis and previous findings, this miRNA could be considered as a potential therapeutic target for IgAN.

Another layer of the constructed network was specified to lncRNAs as regulatory molecules showing outstanding potential in pre- and post-transcriptional regulations [24]. Despite the regulatory role of lncRNAs, their biomarker potential in different types of kidney diseases like membrane nephropathy [55], diabetic nephropathy [56] and IgAN nephropathy [57] have been proposed by several research teams. In this investigation, after importing different curated lncRNAs in the constructed network, HOTAIR was introduced as the topmost lncRNA in regulating the DEGs and DEmiRs in IgAN pathogenicity.

HOTAIR as an anti-sense lncRNA is involved in epigenetic silencing of various genes by recruiting the PRC2 complex, trimethylation of H3K27 across the HOXD locus and finally repressing the transcription [58]. According to recent findings, transcriptional repression of these RNA molecules in immune cells could lead to NF-κB activation and consequently inflammatory responses [59]. Likewise, the potential role of HOTAIR in provoking inflammation was observed in a study on rheumatoid arthritis patients [60]. Here, we also proposed the probable role of this lncRNA in IgAN disease and introduced this RNA molecule as a key target of further analysis in IgAN disease.

Conclusion

In conclusion, the results of this survey add a piece of evidence to the possible involvement of some pathways like NGF signaling and toll-like receptor pathways in the pathogenesis of IgAN. The constructed multilayer regulatory network introduced several hub-DEGs (Tp53, STAT3, Jun, etc.), DEmiRs (miR-124, let-7b, etc.), TFs (NF-kB, etc.), and lncRNAs (HOTAIR, etc.) as potential factors in the pathogenesis of IgAN. MicroRNA-124 and HOTAIR as two RNA molecules were introduced as potential non-invasive biomarkers in IgAN disease. All in all, the selected hub-molecules from the constructed multi-layer regulatory network could be reliable choices of further investigations and validations aiming to diagnose, treat and clarify the hidden pathological aspects of IgAN, this enigmatic disorder.

Availability of data and materials

The datasets analyzed by the current study are available in the GEO repository, [https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE73953, https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE25590].

References

  1. 1.

    Lee H, Hwang JH, Paik JH, Ryu HJ, Kim DK, Chin HJ, et al. Long-term prognosis of clinically early IgA nephropathy is not always favorable. BMC Nephrol. 2014;15(1):94.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  2. 2.

    Suzuki H, Yasutake J, Makita Y, Tanbo Y, Yamasaki K, Sofue T, et al. IgA nephropathy and IgA vasculitis with nephritis have a shared feature involving galactose-deficient IgA1-oriented pathogenesis. Kidney Int. 2018;93(3):700–5.

    CAS  PubMed  Article  Google Scholar 

  3. 3.

    Sun S, Zhao A, Li R, Du R, He L, Sun W, et al. CD147 renal expression as a biomarker for progressive IgAN. J Nephrol. 2015;28(3):307–14.

    CAS  PubMed  Article  Google Scholar 

  4. 4.

    Xie J, Kiryluk K, Wang W, Wang Z, Guo S, Shen P, et al. Predicting progression of IgA nephropathy: new clinical progression risk score. PLoS One. 2012;7(6):e38904.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  5. 5.

    Suzuki H. Biomarkers for IgA nephropathy on the basis of multi-hit pathogenesis. Clin Exp Nephrol. 2019;23(1):26–31.

    CAS  PubMed  Article  Google Scholar 

  6. 6.

    Guo N, Zhou Q, Huang X, Yu J, Han Q, Nong B, et al. Identification of differentially expressed circulating exosomal lncRNAs in IgA nephropathy patients. BMC Immunol. 2020;21(1):1–9.

    Article  Google Scholar 

  7. 7.

    Neprasova M, Maixnerova D, Novak J, Reily C, Julian BA, Boron J, Novotny P, Suchanek M, Tesar V, Kacer P. Toward noninvasive diagnosis of IgA nephropathy: a pilot urinary metabolomic and proteomic study. Dis Markers. 2016;2016(3650909):9.

  8. 8.

    Selewski DT, Ambruzs JM, Appel GB, Bomback AS, Matar RB, Cai Y, et al. Clinical characteristics and treatment patterns of children and adults with IgA nephropathy or IgA vasculitis: findings from the cureGN study. Kidney Int Rep. 2018;3(6):1373–84.

    PubMed  PubMed Central  Article  Google Scholar 

  9. 9.

    Zhai YL, Zhu L, Shi SF, Liu LJ, Lv JC, Zhang H. Increased APRIL Expression Induces IgA1 Aberrant Glycosylation in IgA Nephropathy. Medicine. 2016;95(11):e3099.

  10. 10.

    Maixnerova D, Reily C, Bian Q, Neprasova M, Novak J, Tesar V. Markers for the progression of IgA nephropathy. J Nephrol. 2016;29(4):535–41.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  11. 11.

    Chun J, Chung H, Wang X, Barry R, Taheri ZM, Platnich JM, et al. NLRP3 localizes to the tubular epithelium in human kidney and correlates with outcome in IgA nephropathy. Sci Rep. 2016;6:24667.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  12. 12.

    Sheng X, Zuo X, Liu X, Zhou Y, Sun X. Crosstalk between TLR4 and Notch1 signaling in the IgA nephropathy during inflammatory response. Int Urol Nephrol. 2018;50(4):779–85.

    CAS  PubMed  Article  Google Scholar 

  13. 13.

    Li H-Z, Xu X-H, Lu H-D. Identification of key genes and construction of CircRNA–miRNA–mRNA regulatory networks in osteoarthritis. Electron J Biotechnol. 2019;37:47–55.

    CAS  Article  Google Scholar 

  14. 14.

    Xu B-y, Meng S-J, Shi S-F, Liu L-J, Lv J-C, Zhu L, et al. MicroRNA-21-5p participates in IgA nephropathy by driving T helper cell polarization. J Nephrol. 2020;33(3):551–60.

    CAS  PubMed  Article  Google Scholar 

  15. 15.

    Zheng N, Fan J, Wang B, Wang D, Feng P, Yang Q, et al. Expression profile of BAFF in peripheral blood from patients of IgA nephropathy: correlation with clinical features and streptococcus pyogenes infection. Mol Med Rep. 2017;15(4):1925–35.

    CAS  PubMed  Article  Google Scholar 

  16. 16.

    Ruszkowski J, Lisowska KA, Pindel M, Heleniak Z, Dębska-Ślizień A, Witkowski JM. T cells in IgA nephropathy: role in pathogenesis, clinical significance and potential therapeutic target. Clin Exp Nephrol. 2019;23(3):291–303.

    CAS  PubMed  Article  Google Scholar 

  17. 17.

    Chin C-H, Chen S-H, Wu H-H, Ho C-W, Ko M-T, Lin C-Y. cytoHubba: identifying hub objects and sub-networks from complex interactome. BMC Syst Biol. 2014;8 Suppl 4(Suppl 4):S11.

    PubMed  Article  Google Scholar 

  18. 18.

    Clough E, Barrett T. The gene expression omnibus database. In: Statistical Genomics. Springer; 2016. p. 93–110.

    Google Scholar 

  19. 19.

    Pathan M, Keerthikumar S, Ang CS, Gangoda L, Quek CY, Williamson NA, et al. FunRich: an open access standalone functional enrichment and interaction network analysis tool. Proteomics. 2015;15(15):2597–601.

    CAS  PubMed  Article  Google Scholar 

  20. 20.

    Szklarczyk D, Franceschini A, Kuhn M, Simonovic M, Roth A, Minguez P, et al. The STRING database in 2011: functional interaction networks of proteins, globally integrated and scored. Nucleic Acid Res. 2010;39(suppl_1):D561–8.

    PubMed  PubMed Central  Google Scholar 

  21. 21.

    Han H, Cho J-W, Lee S, Yun A, Kim H, Bae D, et al. TRRUST v2: an expanded reference database of human and mouse transcriptional regulatory interactions. Nucleic Acids Res. 2018;46(D1):D380–6.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  22. 22.

    Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, Ramage D, et al. Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 2003;13(11):2498–504.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  23. 23.

    Yeung KY, Ruzzo WL. Principal component analysis for clustering gene expression data. Bioinformatics. 2001;17(9):763–74.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  24. 24.

    Fang Y, Fullwood MJ. Roles, functions, and mechanisms of long non-coding RNAs in cancer. Genom Proteomics Bioinformatics. 2016;14(1):42–54.

    Article  Google Scholar 

  25. 25.

    Qian W, Xiaoyi W, Zi Y. Screening and bioinformatics analysis of IgA nephropathy gene based on GEO databases. Biomed Res Int. 2019;2019(8794013):7.

  26. 26.

    Chen X, Sun M. Identification of key genes, pathways and potential therapeutic agents for IgA nephropathy using an integrated bioinformatics analysis. J Renin Angiotensin Aldosterone Syst. 2020;21(2):1470320320919635.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  27. 27.

    Cui Y, Liu S, Cui W, Gao D, Zhou W, Luo P. Identification of potential biomarkers and therapeutic targets for human IgA nephropathy and hypertensive nephropathy by bioinformatics analysis. Mol Med Rep. 2017;16(3):3087–94.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  28. 28.

    Miraji MK, Cheng Y, Ge S, Xu G. Identification of primary genes in glomeruli compartment of immunoglobulin a nephropathy by bioinformatic analysis. PeerJ. 2019;7:e7067.

    PubMed  PubMed Central  Article  Google Scholar 

  29. 29.

    Liu Y, Liu X, Jia J, Zheng J, Yan T. Comprehensive analysis of aberrantly expressed profiles of mRNA and its relationship with serum galactose-deficient IgA1 level in IgA nephropathy. J Transl Med. 2019;17(1):320.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  30. 30.

    Vousden KH, Lane DP. p53 in health and disease. Nat Rev Mol Cell Biol. 2007;8(4):275–83.

    CAS  Article  Google Scholar 

  31. 31.

    Hafner A, Bulyk ML, Jambhekar A, Lahav G. The multiple mechanisms that regulate p53 activity and cell fate. Nat Rev Mol Cell Biol. 2019;20(4):199–210.

    CAS  PubMed  Article  Google Scholar 

  32. 32.

    Qin S, Yang C, Zhang B, Li X, Sun X, Li G, et al. XIAP inhibits mature Smac-induced apoptosis by degrading it through ubiquitination in NSCLC. Int J Oncol. 2016;49(4):1289–96.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  33. 33.

    Zilfou JT, Lowe SW. Tumor suppressive functions of p53. Cold Spring Harb Perspect Biol. 2009;1(5):a001883.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  34. 34.

    Cox SN, Serino G, Sallustio F, Blasi A, Rossini M, Pesce F, et al. Altered monocyte expression and expansion of non-classical monocyte subset in IgA nephropathy patients. Nephrol Dialysis Transplant. 2015;30(7):1122–32.

    CAS  Article  Google Scholar 

  35. 35.

    Kawashima H, Takatori H, Suzuki K, Iwata A, Yokota M, Suto A, et al. Tumor suppressor p53 inhibits systemic autoimmune diseases by inducing regulatory T cells. J Immunol. 2013;191(7):3614–23.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  36. 36.

    Cox SN, Sallustio F, Serino G, Pontrelli P, Verrienti R, Pesce F, et al. Altered modulation of WNT–β-catenin and PI3K/Akt pathways in IgA nephropathy. Kidney Int. 2010;78(4):396–407.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  37. 37.

    Coppo R, Camilla R, Alfarano A, Balegno S, Mancuso D, Peruzzi L, et al. Upregulation of the immunoproteasome in peripheral blood mononuclear cells of patients with IgA nephropathy. Kidney Int. 2009;75(5):536–41.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  38. 38.

    Cox SN, Sallustio F, Serino G, Loverre A, Pesce F, Gigante M, et al. Activated innate immunity and the involvement of CX3CR1–fractalkine in promoting hematuria in patients with IgA nephropathy. Kidney Int. 2012;82(5):548–60.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  39. 39.

    Suzuki H, Suzuki Y, Narita I, Aizawa M, Kihara M, Yamanaka T, et al. Toll-like receptor 9 affects severity of IgA nephropathy. J Am Soc Nephrol. 2008;19(12):2384–95.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  40. 40.

    Kajiyama T, Suzuki Y, Kihara M, Suzuki H, Horikoshi S, Tomino Y. Different pathological roles of toll-like receptor 9 on mucosal B cells and dendritic cells in murine IgA nephropathy. Clin Dev Immunol. 2011;2011.

  41. 41.

    Lim BJ, Lee D, Hong SW, Jeong HJ. Toll-like receptor 4 signaling is involved in IgA-stimulated mesangial cell activation. Yonsei Med J. 2011;52(4):610.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  42. 42.

    Anders H-J, Banas B, Schlöndorff D. Signaling danger: toll-like receptors and their potential roles in kidney disease. J Am Soc Nephrol. 2004;15(4):854–67.

    CAS  PubMed  Article  Google Scholar 

  43. 43.

    Skaper SD. Nerve growth factor: a neuroimmune crosstalk mediator for all seasons. Immunology. 2017;151(1):1–15.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  44. 44.

    Kimata H, Yoshida A, Ishioka C, Kusunoki T, Hosoi S, Mikawa H. Nerve growth factor specifically induces human IgG4 production. Eur J Immunol. 1991;21(1):137–41.

    CAS  PubMed  Article  Google Scholar 

  45. 45.

    Otten U, Ehrhard P, Peck R. Nerve growth factor induces growth and differentiation of human B lymphocytes. Proc Natl Acad Sci. 1989;86(24):10059–63.

    CAS  PubMed  Article  Google Scholar 

  46. 46.

    Antonucci MT, Bonofiglio R, Papalia T, Caruso F, Caroleo MC, Mancuso D, Aloe L: Nerve growth factor and its monocyte receptors are affected in kidney disease. Nephron Clin Pract 2009, 111(1):c21-c28.

  47. 47.

    Vergano L, Loiacono E, Albera R, Coppo R, Camilla R, Peruzzi L, et al. Can tonsillectomy modify the innate and adaptive immunity pathways involved in IgA nephropathy? J Nephrol. 2015;28(1):51–8.

    CAS  PubMed  Article  Google Scholar 

  48. 48.

    Gan T, Li Y, Zhou X-J, Zhang H. Immunoproteasome in IgA nephropathy: state-of-art and future perspectives. Int J Biol Sci. 2020;16(14):2518.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  49. 49.

    Wang X, Li Z, Bai J, Song W, Zhang F. miR-17-5p regulates the proliferation and apoptosis of human trabecular meshwork cells by targeting phosphatase and tensin homolog. Mol Med Rep. 2019;19(4):3132–8.

    CAS  PubMed  PubMed Central  Google Scholar 

  50. 50.

    Zhu S, Zhang L, Zhao Z, Fu W, Fu K, Liu G, et al. MicroRNA-92a-3p inhibits the cell proliferation, migration and invasion of Wilms tumor by targeting NOTCH1. Oncol Rep. 2018;40(2):571–8.

    CAS  PubMed  PubMed Central  Google Scholar 

  51. 51.

    Zhan X-H, Xu Q-Y, Tian R, Yan H, Zhang M, Wu J, et al. MicroRNA16 regulates glioma cell proliferation, apoptosis and invasion by targeting Wip1-ATM-p53 feedback loop. Oncotarget. 2017;8(33):54788.

    PubMed  PubMed Central  Article  Google Scholar 

  52. 52.

    Mizuno R, Kawada K, Sakai Y. The molecular basis and therapeutic potential of Let-7 microRNAs against colorectal cancer. Can J Gastroenterol Hepatol. 2018;2018.

  53. 53.

    Serino G, Sallustio F, Curci C, Cox SN, Pesce F, De Palma G, et al. Role of let-7b in the regulation of N-acetylgalactosaminyltransferase 2 in IgA nephropathy. Nephrol Dialysis Transplant. 2015;30(7):1132–9.

    CAS  Article  Google Scholar 

  54. 54.

    Odani H, Hiki Y, Takahashi M, Nishimoto A, Yasuda Y, Iwase H, et al. Direct evidence for decreased Sialylation and Galactosylation of human serum IgA1 fc O-glycosylated hinge peptides in IgA nephropathy by mass spectrometry. Biochem Biophys Res Commun. 2000;271(1):268–74.

    CAS  PubMed  Article  Google Scholar 

  55. 55.

    Huang Y-S, Hsieh H-Y, Shih H-M, Sytwu H-K, Wu C-C. Urinary Xist is a potential biomarker for membranous nephropathy. Biochem Biophys Res Commun. 2014;452(3):415–21.

    CAS  PubMed  Article  Google Scholar 

  56. 56.

    Alvarez ML, Khosroheidari M, Eddy E, Kiefer J. Role of microRNA 1207-5P and its host gene, the long non-coding RNA Pvt1, as mediators of extracellular matrix accumulation in the kidney: implications for diabetic nephropathy. PLoS One. 2013;8(10):e77468.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  57. 57.

    Zuo N, Li Y, Liu N, Wang L. Differentially expressed long non-coding RNAs and mRNAs in patients with IgA nephropathy. Mol Med Rep. 2017;16(5):7724–30.

    CAS  PubMed  Article  Google Scholar 

  58. 58.

    Hajjari M, Salavaty A. HOTAIR: an oncogenic long non-coding RNA in different cancers. Cancer Biol Med. 2015;12(1):1.

    CAS  PubMed  PubMed Central  Google Scholar 

  59. 59.

    Mandal SS, Obaid M, Chini A. Long non-coding RNAs regulate NF-κB activation, inflammation, and glucose metabolism. FASEB J. 2020;34(S1):1–1.

    Google Scholar 

  60. 60.

    Song J, Kim D, Han J, Kim Y, Lee M, Jin E-J. PBMC and exosome-derived Hotair is a critical regulator and potent marker for rheumatoid arthritis. Clin Exp Med. 2015;15(1):121–6.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

Download references

Acknowledgments

The authors want to thank the members of the Regenerative Medicine Research Center of Isfahan University of Medical Sciences for their kind support.

Funding

This work is partly financed by Isfahan University of Medical Sciences (Grant no: 199182).

Author information

Affiliations

Authors

Contributions

Study concept and design by Amir roointan, Alieh Gholaminejad and Yousof Gheisari. Acquisition of data by Amir roointan, Alieh Gholaminejad. Analysis and interpretation of data by Amir roointan, Sedigheh Jalali, Alieh Gholaminejad and Yousof Gheisari. Drafting of the manuscript by Amir Roointan, Sedigheh Jalali, Alieh Gholaminejad and Yousof Gheisari. Critical revision of the manuscript for important intellectual content by Amir Roointan, Alieh Gholaminejad and Yousof Gheisari. Study supervision by Amir roointan. The author(s) read and approved the final manuscript.

Author information

Alieh Gholaminejad, Yousof Gheisari, and Amir Roointan from “Regenerative medicine research center, Isfahan University of Medical Sciences, Isfahan, Iran” and Sedigheh Jalali from “Department of Pediatrics, The University of Melbourne, Melbourne, Australia” were contributed in this study as follow:

Corresponding author

Correspondence to Amir Roointan.

Ethics declarations

Ethics approval and consent to participate

The Ethics Committee of Isfahan University of Medical Sciences approved the study and waived the need of informed written/verbal consent for this study.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Gholaminejad, A., Gheisari, Y., Jalali, S. et al. Comprehensive analysis of IgA nephropathy expression profiles: identification of potential biomarkers and therapeutic agents. BMC Nephrol 22, 137 (2021). https://doi.org/10.1186/s12882-021-02356-4

Download citation

Keywords

  • IgA nephropathy
  • Computational biology
  • Gene expression
  • Gene regulatory network
  • Biomarkers
\