Skip to main content

Identification of diagnostic markers and immune cell infiltration characteristics in antineutrophil cytoplasmic antibody-associated vasculitis by weighted gene co-expression network analysis

Abstract

Background

Antineutrophil cytoplasmic antibody-associated vasculitis (AAV) is a group of life-threatening systemic autoimmune diseases. The aim of this study was to determine the relationship between the AAV hub gene and immune cell infiltration, and its value for clinical disease treatment.

Methods

We downloaded the microarray information of 37 AAV patients and 27 controls from Gene Expression Omnibus (GEO). Genes were classified into totally different modules exploitation weighted gene co-expression network analysis (WGCNA). AAV diagnostic indicators were screened and then assessed immune cell infiltration by the least absolute shrinkage and selection operator (LASSO) and CIBERSORT. Finally, Connectivity Map analysis was applied to predict possible AAV glomerulus injury improvement therapies.

Results

WGCNA was developed and differentially expressed genes were classified into 6 modules, the black module was most tightly correlated to AAV. Among them, TIMP1 and FCER1G were most closely related to clinical features. Resting mast cells and monocytes emerged as having the foremost distinguished variations in AAV. C3AR1 and FCER1G were involved in AAV development by immune regulation. Connectivity Map analysis indicated the most significant compound was fisetin.

Conclusions

The present study is that the initial to spot immune cell infiltration with microarray data of glomeruli in AAV, which provides novel proof and clues for additional analysis of the molecular mechanisms.

Introduction

Antineutrophil cytoplasmic antibody (ANCA)-associated vasculitis (AAV) is a group of life-threatening systemic autoimmune diseases characterized by inflammation and destruction of small- and medium-sized blood vessels, including three types: microscopic polyangiitis (MPA), granulomatosis with polyangiitis (GPA), and eosinophilic granulomatosis with polyangiitis (EGPA) [1]. Similar to other autoimmune diseases, the mechanism by that ANCAs cause vasculitis involves ANCA-mediated excessive activation of immune cells consisting of innate and adaptive immune populations [2]. Activation or inactivation of various types of immune cells in AAV contributes to disease initiation and augmentation by regulating the suppression, maintenance, or promotion of immune responses [2, 3]. AAV is usually responsible for grievous kidney failure or pulmonary hemorrhage [1], especially for the kidney, with > 75% of AAV patients have renal involvement, which is called ANCA-associated glomerulonephritis (GN) [4, 5]. In that case, GPA and MPA are the major types in AAV, both leading a pauci-immune necrotizing glomerulonephritis [1]. As an important complication of AAV, ANCA-associated GN directly affects the clinical outcome of the patients [6]. Preventing or ameliorating kidney injury has become an essential part of AAV treatment [7]. Given the functionally distinct cell sorts that comprise the immune reaction, assessing immune infiltration and deciding whether or not variations within the composition of the immune infiltration will improve the event of novel therapy medicine to focus on these cells is vital in AAV and ANCA-associated GN [8]. Although advances in treatment methods and the prognosis of patients have improved over the past decades, there are still a relatively large number of patients who end up in End-Stage Renal Disease (ESRD) present with low eGFR and do not recover renal function [9] and a serious burden of morbidity and mortality in AAV [10]. Hence, there is a requirement for new strategies for exploring AAV.

Gene sequence technologies and bioinformatic analyses have been used in recent years for the identification of disease-related genes which could be used for prognostic biomarks and produced in future as therapeutic targets [11]. Therefore, tracking the biological changes in AAV at the genomic level is a valuable technique [12]. Although several previous studies had focused on genes closely associated with AVV [13, 14], transcriptomics at the cellular level was few reported. To our knowledge, few studies utilized bioinformatic analysis and immune cell infiltration to characterize kidney tissue in the context of AVV. The aim of this study was to contribute to further mechanistic studies of AAV by analyzing the relationship between the hub gene and immune infiltrating cells and their associated immunological processes in AAV kidney tissue.

Materials and methods

Data collection and preprocessing

The AAV RNA expression data were acquired from the Gene Expression Omnibus (GEO; http://www.ncbi.nlm.nih.gov/geo/). For differential expression analysis, data from GSE108109 and GSE104948 were employed. The GSE 108109 (GPL19983 platform, Affymetrix Human Gene 2.1 ST Array, http://dx.doi.org/10.1136/annrheumdis-2017-212935) contains 21 specimens, including 6 normal controls and 15 specimens of AAV. The GSE 104948 (GPL22945 platform, Affymetrix human genome U133 Plus 2.0 array, http://dx.doi.org/10.1136/annrheumdis-2017-212935) contains 43 specimens, including 21 normal controls and 22 AAV specimens. AAV samples were obtained from kidney tissue specimens from patients with confirmed ANCA-GN diagnosis after kidney puncture pathology, and healthy kidney tissue was obtained from living transplant donors. Among them, the types of AAV are GPA and MPA. A positive ANCA antibody was present in both AAV patients. In order to obtain the matrix of gene expression, R v3.6.1 was used to extract and sort data. The batch effect caused by study heterogeneity was eliminated using the “ComBat module” of the “SVA package,” which utilizes empirical “Bayes” methods. Background adjustments and data normalization were performed with the “limma package” [15].

Identification of differentially expressed genes

We filtered the differentially expressed genes (DEGs) among AAVs and controls using the “limma” R package in the expressing results. The significance analysis of microarrays (SAM) was accustomed choose considerably altered genes with false discovery rate (FDR) < 0.05 and log2 fold change (FC) ≥ 0.5.

Establishment of co-expression network

Based on the R package “WGCNA” the DEG co-expression network was designed [15]. The soft threshold power we selected was eight when 0.8 was the threshold of the correlation and a minimum of 10 genes in the modules were chosen. We have specified 0.2 as the cutting height threshold to combine possibly related modules.

Detection of hub genes

Hub genes with gene significance higher than 0.6 and high module group members (MM) (weighted correlation index > 0.9) were considered in the module trait correlation analysis, demonstrating a significant correlation to certain clinical characteristics.

GSEA-based enriched GO and pathway analysis

Gene Ontology (GO) analysis included biological process (BP), cellular component (CC), and molecular function (MF). Reactome is a primary database for the research of pathways. We selected BP to perform GO analysis, as well as Reactome to perform pathway analysis for the crucial module; it was accomplished in the "cluster Profiler" package via gene set enrichment analysis (GSEA). The cut-off parameters were the adjusted P-value < 0.05.

Assessment of immune cell infiltration

To delete the null values, gene expression datasets had been analyzed. In the "impute" page [16], the missing value was completed with KNN; the format was established in compliance with the recognized CIBERSORT format, followed by the uploading of data to the CIBERSORT website (http://cibersort.stanford.edu/). We used the original CIBERSORT LM22 gene signature file to study datasets, which described 22 immune cell subtypes. Samples that fulfill the CIBERSORT P < 0.05 conditions were included with immune cells profiles.

Analysis of differences among types of immune cell infiltration

We analyzed the main differential expression of different cell immune cell types using the differential analysis between the group AAV and control group. A linear model was used in the “limma” bundle and Bayesian method [17]. The cut-off standard was P-value < 0.05. Spearman's correlating coefficient was used to find a similarity among certain differentially expressed types of immune cells, in order to further see the association between these different immune cells types.

LASSO model and receiver operating characteristic (ROC) curve analysis

A good predictive value and low correlation were extended to the least absolute shrinkage and selection operator (LASSO), which chose the best features for the high-dimensional data [18]. We also extracted the hub gene expression profile for building LASSO model by “glmnet” package (https://CRAN.R-project.org/package=glmnet) to separate AAV from control. Datasets were at random appointed to the training set (70%) and test set (30%). The “pROC” packet [15] was utilized for analyzing ROC curve in training and test set to assess the ability of the LASSO model to classify AAV.

Statistical and CMap analysis

In the case of AAV, the online portal (http://v5.nephroseq.org) was used for the Spearman correlation analysis among hub genes and glomerular filtration rate (GFR) [19, 20] and serum creatinine (SCr) [21, 22]. Connectivity Map (CMap) (https://portals.broadinstitute.org/cmap) is an open source, linking disease, genes, and medicines through similar or opposite gene expression profiles [23]. A CMap analysis was applied to forecast possible AAV glomerulus injury improvement therapies.

Results

Bioinformatic analysis workflows and DEGs related to AAV

The workflows are shown in Fig. 1. After eliminating the batch effect by cross-platform from the microarray data (Fig. 2), 445 DEGs included in AAV were screened in the heatmap (Fig. 3a) by the limma package (adjusted p < 0.05, |logFC|≥ 0.5). 225 genes were downregulated and 220 genes were upregulated (Fig. 3b).

Fig. 1
figure 1

Workflows of our bioinformatic analysis

Fig. 2
figure 2

Principal component analysis (PCA) of the gene expression datasets. The scatter points show the samples on the base of the top two principal components of gene expression profiles (PC1 and PC2) without a and b removing the batch effect. The colors reflect samples in two different datasets

Fig. 3
figure 3

Difference analysis. a Heatmap (top 30 genes) of 445 DEGs screened by “limma package.” Red represent highly expressed genes and blue represent lowly expressed genes in AAV compared with normal controls. b Volcano plot analysis identifies DEGs. Red dots represent upregulated genes and green dots represent downregulated genes in AAV compared with normal controls. DEG differentially expressed gene; AAV ANCA-associated vasculitis

Weighted gene co-expression network analysis (WGCNA) co-expression networks

The soft thresholding power was chosen as 10, while 0.8 was used as the correlation coefficient threshold (Fig. 4a). 6 co-expression modules were built by WGCNA analysis (Fig. 4b). The module comprised the majority of genes was brown, followed by black, green, and blue (Fig. 4b).

Fig. 4
figure 4

WGCNA revealed gene co-expression networks in the glomerulus of AAV. a Analysis of the scale-free fit index for various soft thresholding powers (Left) and analysis of the mean connectivity for various soft thresholding powers (Right); b Clustering dendrogram of DEGs related to the glomerulus of 37 AAV

AAV module trait correlations and Hub Genes identification

The study of module trait correlations found that multiple modules were correlated to AAV (Fig. 5a). The summary of the significance of all genes associated with AAV in each module is shown in Fig. 5b. AAV (correlation coefficient = 0.75, P = 1E-12; Fig. 5a) was positively correlated to the black module, while the brown module was negatively correlated with AAV (correlation coefficient = -0.81, P = 8E-16; Fig. 5a). In the black module, Fig. 5c shows the significance of these genes. In black modules 7 genes (TIMP1, FCER1G, SH3BGRL3, HCLS1, C3AR1, TYROBP, and CD53) were identified as hub genes, according to GS > 0.6 and MM > 0.9 (Fig. 5c, d).

Fig. 5
figure 5

Main findings in the module trait correlations analyses. a Heatmap representing the relationship between modules and AAV (the correlation coefficient and corresponding P-value were present in each cell); b The module significance values of certain co-expression modules related to AAV (the module significance value expressed a description of the gene significance of all genes in each module, and different column colors represent different modules); c AAV gene significance in the black module (In the black module, one dot represents one gene.); d 7 hub genes in the black module were strongly associated with each other

Enrichment analysis of hub genes

In the black module of the merged dataset, the expression of 7 hub genes, including CD53 (Fig. 6a), FCER1G (Fig. 6b), TYROBP (Fig. 6c), C3AR1 (Fig. 6d), SH3BGRL3 (Fig. 6e), HCLS1 (Fig. 6f), and TIMP1 (Fig. 6g), was significantly increased in the AAV patients. Module function enrichment analysis demonstrated the substantial involvement of the black module genes (Fig. 7a) in biological processes associated with immunity mediated by neutrophils. Pathway analysis revealed that the only pathway enriched in the black module was the immune pathway (Fig. 7b).

Fig. 6
figure 6

Expression of hub genes. (a–g) Expression levels of CD53 a, FCER1G b, TYROBP c, C3AR1 d, SH3BGRL3 e, HCLS1 f, and TIMP1 g were significantly increased in AAV patients

Fig. 7
figure 7

Enrichment analysis of hub genes. a The biological processes in module genes increasingly raise the significance of enrichment from blue to red, and point size shows the number of differential genes in the functional pathway. b The pathway of black module enrichment

GSEA-based analysis of all detected genes

GSEA was conducted to detect genes with a statistically significant difference in AAV. Cell activation in the immune process included the most important enriched genes positively associated with the AAV (Fig. 8). Further evidence showed that the immune process was the most significant cause of AAV. Table1 shows the most critical pathways after the black module GO screening. Among the 7 hub genes screened by the black module, the genes involved in the immune effector process were FCER1G, C3AR1, TYROBP, and CD53.

Fig. 8
figure 8

GSEA-based analysis of the main enrichment processes in AAV

Table 1 GO analysis in the black module filtered with meaningful 4 items of BP

Screening and verification of the potential diagnostic marker

We extracted the expression profile for the LASSO model of the hub genes (Fig. 9a). 2 hub genes (TIMP1 and FCER1G) with non-zero regression coefficients had been defined using the LASSO procedure, and the value of lambda. 1se = 0.103313. The index model of the genes was as follows: index = TIMP1 *(0.4335136) + FCER1G* (0.1505385). The ROC (Fig. 9b) study revealed a 2-gene model AUC of 0.98 for the training and 1.0 for the test sample.

Fig. 9
figure 9

A model for predicting and verification of AAV. a LASSO model. b ROC curves analysis of train and test set

Performance of CIBERSORT

Figure 10a shows the ratio of immune cells in 37 renal tissues, among which the most infiltrated cells in AAV were mast cells. It showed the difference in the expression of the immune infiltration cells in the AAV and control groups in Fig. 10b. 11 immunocellular types were differential expression, namely plasma cell T, naive CD4 T cell, activated memory CD4 T cell, regulatory T cell (Tregs), γδT cell, monocyte, M1 macrophage, M2 macrophagic, resting dendritic cell, resting mast cell, and activated mast cell. Tregs, γδT cells, monocytes, M2 macrophages, and resting mast cells were upregulated in ANCA-associated GN tissue. Among them, the increase in resting mast cells and monocytes was the most significant. Plasma cells, naive CD4 T cells, activated memory CD4 T cells, M1 macrophages, resting dendritic cells, and activated mast cells were downregulated. The relation between these differentially expressing immune cell types is shown in Fig. 10c. Resting mast cells were negatively related to with plasma cells and M1 macrophages (r =  − 0.6 and r =  − 0.57, respectively), that indicated that the performance of resting mast cells, plasma cells, and M1 macrophages in AAV was antagonistic. The relationship between plasma memory cells and M1 macrophages was synergistic.

Fig. 10
figure 10

Landscape of immune infiltration in AAV. a Bar charts of 22 immune cell ratios in AAV and normal tissues. b Differential expression of different types of immune cells between AAV (green) and normal tissues (blue). c Correlation matrix of 11 types of immune cell proportions. Data were collated by using R package “tidyverse” (version 1.2.1). R package “ggpubr” (version 0.1.8) was used for T test. Results visualization was performed by using R package “ggplot2” (version 3.1.0). Correlation analysis and visualization were showed by using R package “corrplot” (version 0.84)

Discovery of core genes

The correlation between the core genes of the “immune effector process” and 11 types of the immune cell infiltration process is shown in Fig. 11a. Among 44 genes that showed close association with immune infiltrating cells, 4 hub genes involved in the “immune effector process” in the black module, those associated with differentially immune infiltrating cells were C3AR1 positively correlated with M2 Macrophages, FCER1G negatively correlated with naive CD4 T cells, and no significant correlation was found for the remaining hub genes. Analysis at the human protein atlas website revealed that C3AR1 and FCER1G were significantly highly expressed in C9 Macrophages of renal tissue (Fig. 11b, c).

Fig. 11
figure 11

Correlation analysis between key genes in immune effector process and 11 types of immune infiltrating cells and the expression of C3AR1 and FCER1G in the kidney. a The vertical items are the names of immune cells. The horizontal items indicate the correlation coefficient (different colors represent different immune cell types). b The expression of C3AR1 in the kidney. c The expression of FCER1G in the kidney

Association between hub genes and clinical features of AAV

After online “Nephroseq v5 platform” analysis, the results showed that mRNA expression of TIMP1 (Fig. 12a), C3AR1 (Fig. 12b), CD53 (Fig. 12c), FCER1G (Fig. 12d), HCLS1 (Fig. 12e), SH3BGRL3 (Fig. 12f), and TYROBP (Fig. 12g) in glomerulus reversely correlated with glomerular filtration rate (GFR) and positively correlated with serum creatinine (SCr) in AAV patients. In this case, upregulated hub genes were required to enhance the implementation of AAV.

Fig. 12
figure 12

Correlation between mRNA expression of 7 hub genes in glomerulus and GFR, SCr in AAV. TIMP1 a, C3AR1 b, CD53 c, FCER1G d, HCLS1 e, SH3BGRL3 f, and TYROBP g correlated with GFR and positively correlated with SCr in AAV

Identification of agents to prevent AAV injury

CMap analytics was used to identify small molecular compounds that reversed the expression of 7 hub genes in cell levels to predict potential medications for the improvement of the disease. The most significant compounds are shown in Table 2.

Table 2 Small molecular compounds provided by CMap analysis to reverse expression

Discussion

So far, our study is the first to use WGCNA to probe the gene-network signature of kidney tissue to AAV. Two mRNA microarray datasets downloaded from GEO were enclosed for further analyses. DEGs between renal glomerulus tissues of AAV patients and normal controls were screened to work out the hub genes within the most connected module. The study of GSEA examined possible underlying causes of ANCA glomerulus damage. We examined the variance between AAV kidney tissue and normal kidney tissue in immune cell by using CIBERSORT. Furthermore, the Nephroseq v5 online tool was used for correlation study of hub genes and clinical features of AAV. The CMap analysis offered small molecular compounds to reverse altered expression of hub gene in cell lines.

Most patients with AAV have renal involvement, and both patient survival and risk of ESRD are closely related to renal function at the time of presentation [24]. It has long been recognized that both environmental and genetic factors contribute to the development of AAV. Moreover, different types of AAV have genetic differences. Multiple studies found that anti-proteinase 3 ANCA was correlated with HLA-DP and serpin A1 gene (SERPINA1) and proteinase 3 (PRTN3), but anti-myeloperoxidase ANCA was linked with HLA-DQ, which depended on antigenic specificity of ANCA [25,26,27]. Recent clinical and basic science studies have demonstrated the roles of neutrophils, other immune cells, and humoral factors in the pathogenesis of AAV and improved to several breakthroughs in the grasp and treatment of this disease [28]. Although in-depth efforts have been created, the underlying mechanisms of AAV renal injury stay elusive. The utilization of bioinformatics techniques has allowed us to identify important biomarkers in the development and progression of AAV, which may provide new insight into the further study of this disease.

In this study, we detected 445 DEGs between kidney biopsy specimens of AAV patients and normal controls based on two microarray datasets. WGCNA clusters DEGs in the 6 most relevant gene modules. After analyzing the gene significance between each module and AAV, the interaction of the black module is significantly higher than other modules. Enrichment analyses revealed that major genes in this module joined in various immune responses. As antecedently according, immune responses, like the activation of adaptive immune system triggering immune advanced deposition, complement activation, and self-antigen production, displayed a harmful impact on renal glomerular cells, that were related to impaired kidney performance in AAV [29].

A total of 7 DEGs were identified as hub genes, including TIMP1, FCER1G, SH3BGRL3, HCLS1, C3AR1, TYROBP, and CD53. Among them, TIMP1 and FCER1G were also identified with non-zero regression coefficients as diagnostic markers of AAV by combining LASSO methods. The gene tissue inhibitor matrix Metalloproteinase 1 (TIMP1), which encodes a 931 base pair mRNA and a 207 amino acid protein [30], is a tissue inhibitor of Metalloproteinases family. This protein can inhibit the proteolytic activity of matrix metalloproteinases (MMPs) by building non-covalent complexes and control the stability of the matrix reshaping during extracellular matrix degradation [31] and its production increases to antagonize the elevated MMPs expression [32, 33]. It has been confirmed MMPs will be upregulated in inflammation and physiological remodeling processes and implicated in both the inflammatory and fibrotic phase of crescent formation in active renal ANCA-related vasculitis [34], and the increased renal TIMP1 expression also imaged the inflammatory process in AVV and ANCA-related GN [35]. FCER1G is involved in innate immunity [36], which shifted activation signals from various immunoreceptors [37, 38]. FCER1G was functionally joined to mediate neutrophil activation and was additionally concerned in platelet activation [39], and improved prognosis by affecting the immune-related pathways in the progression of clear cell renal cell carcinoma (ccRCC). Moreover, FCER1G could be a crucial molecule in signal pathways that are wide concerned in an exceedingly sort of immune responses and cell varieties [40]. Combining the results, we speculate that FCER1G may play a vital role in the disease progression of AAV given its function in immune regulation. The above analysis shows that AAV can be developed with TIMP1 and FCER1G and that they can also be used as diagnostic markers for AAV. Except the two genes mentioned above, our research demonstrated that the expression level of SH3BGRL3, HCLS1, C3AR1, TYROBP, and CD53 all showed an increasing trend in the AAV kidney tissue compared with normal tissue. By Nephroseq v5 online tool analysis, the upregulation of the above-mentioned 7 hub genes was significantly positively correlated with the GFR and SCr in the progression of AVV. Therefore, we prefer to consider them as an indicator of disease prediction. Although, at present, there is no direct evidence that these 5 genes play a specific role in kidney disease but it provides direction and hypothesis for future AAV gene studies.

To further explore the association between hub genes and immune cell infiltration in AVV, we used CIBERSORT. We found an increased infiltration of resting mast cells, monocytes, regulatory T cells (Tregs), αβ T‐cell, and M2 macrophages, while a decreased infiltration of plasma cells, naive CD4 T cells, activated memory CD4 T cells, M1 macrophages, resting dendritic cells, and activated mast cells may be related to the occurrence and development of AAV. Mast cells are critical antigen-presenting cells that have arisen as potential main players in general and progressively progressive ways of GN activation, amplification, and paradoxically immunomodulation [41, 42]. Confirmed through in vivo experiments that the earlier vasculitis at 48 h was thought to be mast cell dependent [43], and MCs are the largest populations of leukocytes in AAV infiltration of injured kidneys and the most serious tubulointerstitial injury has occurred to patients with higher MC density [44]. The previous experimental studies suggested MCs can improve the development of immunomodulation by interacting with Tregs in AAV [45], which was presented in the form of a strong positive correlation between mast cells and Tregs in our study. In addition, the dominant form of infiltrating cell in the early biopsies of ANCA-associated GN is CD68 + monocytes/macrophages. In this setting, monocytes can be used to orchestrate the immunologic response [46] and triggered in vasculitis [47]. Our results were also consistent with the previously reported skewed balance of CD4 + T cells, including the increase of Tregs and the decrease of naive cells in AAV patients [48,49,50]. Another interesting finding was that neutrophils had always been thought to be closely related to AVV [51]. However, there were no significant differences in the expression of neutrophils in the present study, which may be related to the formation of neutrophil extracellular traps (NETs), resulting in the number of activated neutrophils had increased in the meantime, while the total number of neutrophils remained almost unchanged [52]. Previous evidence in combination with our findings has shown the importance of resting mast cells, monocytes, regulatory T cells, and activated neutrophils in AAVs and should be the highlight of further research. Furthermore, our results demonstrate 11 types of AAV immune cells. The mast cells are closely linked to the infiltration of regulatory T cells and plasma cells. However, there is no relevant research report on it, further scientific evidence is essential for the basic processes of such associations. In an analysis of the relationship between the Hub genes and immune cells, the C3AR1 association with the M2 Macrophages was found to be extremely positive. FCER1G was significantly negatively correlated with naive CD4 T cells. Studies have shown that M2 macrophages are participated in acute renal injury of glomerulonephritis with crescents especially in AAV [53, 54]. Therefore, combine results about CD4 + T cells indicated that C3AR1 and FCER1G may be involved in AAV development by immune regulation. This requires further study in order to clarify the dynamic possible relationship between genes and immune cells.

We uploaded 7 hub genes into the CMap database and matched them with small molecule treatment. As a natural flavonoid, fisetin was confirmed to prevent high fat diet-induced diabetic nephropathy through the suppression of insulin resistance and the inflammation mediated by RIP3 [55], and reduce the inflammation of kidneys and apoptoses by inhibiting the Src-mediated NF-κB p65 and MAPK pathways to guard against the LPS-induced septic AKI mice [56]. As yet, no relevant studies have shown the effect of these compounds on kidney damage in AAV.

Although our study provides new ideas for the diagnosis and treatment of AAV, it still has some limitations. Firstly, as we know, there are three different types of AAV and that they are all genetically different from each other. However, based on the results of our study, no correlation between the above-mentioned genes that distinguish AAV types and immunomodulatory cells was found. For further studies, we will work on an in-depth analysis of the subtle differences between the different types of AAV and verify our current findings and clarify the biological functions of the AAV genes by basic studies. Second, we need to perform in vivo experiments to verify the specific mechanism of action of hub gene with differential immune cell. Third, there are not enough available sequencing datasets for AAV disease model, and we need to further refine the lack of this field.

Conclusion

Taken together, our study identified TIMP1 and FCER1G had the potential to be used as diagnostic markers for AAV. 11 types of immune cells revealed important associations with AAV, resting mast cells and monocytes showed the greatest variations. In this case, C3AR1and FCER1G may be involved in AAV development by immune regulation.

Data availability

Datasets used and/or analyzed in this study are available from the corresponding author on reasonable request.

Abbreviations

AAV:

Antineutrophil cytoplasmic antibody (ANCA)-associated vasculitis

GEO:

Gene expression omnibus

GSEA:

Gene set enrichment analysis

LASSO:

The least absolute logistic regression of the shrinkage and selection operator

CMap:

Connectivity map

DEGs:

Differentially expressed genes

GO:

Gene ontology

WGCNA:

Weighted gene co-expression network analysis

PCA:

Principal component analysis

ROC:

Receiver operating characteristic

MPA:

Microscopic polyangiitis

GPA:

Granulomatosis with polyangiitis

EGPA:

Eosinophilic granulomatosis with polyangiitis

ANCA-GN:

ANCA-associated glomerulonephritis

ESRD:

End-stage renal-disease

SAM:

Significance analysis of microarrays

FDR:

False discovery rate

MM:

Module group members

CC:

Cellular component

MF:

Molecular function

GFR:

Glomerular filtration rate

SCr:

Serum creatinine

SERPINA1:

Serpin A1 gene

PRTN3:

Proteinase 3

TIMP1:

Metalloproteinase 1

MMPs:

Matrix metalloproteinases

ccRCC:

Clear cell renal cell carcinoma

NETs:

Neutrophil extracellular traps

References

  1. Lyons PA, Rayner TF, Trivedi S, Holle JU, Watts RA, Jayne DR. Genetically distinct subsets within ANCA-associated vasculitis. N Engl J Med. 2012;367(3):214–23. https://doi.org/10.1056/NEJMoa1108735.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  2. Nakazawa D, Masuda S, Tomaru U, Ishizu A. Pathogenesis and therapeutic interventions for ANCA-associated vasculitis. Nat Rev Rheumatol. 2019;15(2):91–101. https://doi.org/10.1038/s41584-018-0145-y.

    CAS  Article  PubMed  Google Scholar 

  3. Rahmattulla C, Mooyaart AL, van Hooven D, Schoones JW, Bruijn JA, Dekkers OM, Bajema IM, European Vasculitis Genetics Consortium. Genetic variants in ANCA-associated vasculitis: a meta-analysis. Ann Rheum Dis. 2016;75(9):1687–92. https://doi.org/10.1136/annrheumdis-2015-207601.

    Article  PubMed  Google Scholar 

  4. Chen M, Jayne D, Zhao MH. Complement in ANCA-associated vasculitis: mechanisms and implications for management. Nat Rev Nephrol. 2017;13(6):359–67. https://doi.org/10.1038/nrneph.2017.37.

    CAS  Article  PubMed  Google Scholar 

  5. Geetha D, Jefferson JA. ANCA-associated vasculitis: core curriculum 2020. Am J Kidney Dis. 2020;75(1):124–37. https://doi.org/10.1053/j.ajkd.2019.04.031.

    CAS  Article  PubMed  Google Scholar 

  6. Geetha D, Specks U, Stone JH. Rituximab versus cyclophosphamide for ANCA-associated vasculitis with renal involvement. J Am Soc Nephrol. 2015;26(4):976–85. https://doi.org/10.1681/ASN.2014010046.

    CAS  Article  PubMed  Google Scholar 

  7. Koneth I. Pauciimmune glomerulonephritis ANCA-vaskulitis. Ther Umsch. 2015;72(3):161–9. https://doi.org/10.1024/0040-5930/a000660.

    Article  PubMed  Google Scholar 

  8. Xiong Y, Wang K, Zhou H. Profiles of immune infiltration in colorectal cancer and their clinical significant: a gene expression-based study. Cancer Med. 2018;7(9):4496–508. https://doi.org/10.1002/cam4.1745.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  9. Kronbichler A, Leierer J, Gauckler P. Comorbidities in ANCA-associated vasculitis. Rheumatology. 2020. https://doi.org/10.1093/rheumatology/kez617.

    Article  PubMed  PubMed Central  Google Scholar 

  10. Yates M, Watts R. ANCA-associated vasculitis. Clin Med (Lond). 2017;17(1):60–4. https://doi.org/10.7861/clinmedicine.17-1-60.

    Article  Google Scholar 

  11. Khan MI, Dębski KJ, Dabrowski M. Gene set enrichment analysis and ingenuity pathway analysis of metastatic clear cell renal cell carcinoma cell line. Am J Physiol Renal Physiol. 2016. https://doi.org/10.1152/ajprenal.00138.2016.

    Article  PubMed  Google Scholar 

  12. Xie G, Roshandel D, Sherva R. Association of granulomatosis with polyangiitis (Wegener’s) with HLA-DPB1*04 and SEMA6A gene variants: evidence from genome-wide analysis. Arthritis Rheum. 2013;65(9):2457–68. https://doi.org/10.1002/art.38036.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  13. Blitzer RD, Gil O, Landau EM. Cholinergic stimulation enhances long-term potentiation in the CA1 region of rat hippocampus. Neurosci Lett. 1990;119(2):207–10. https://doi.org/10.1016/0304-3940(90)90835-w.

    CAS  Article  PubMed  Google Scholar 

  14. Zou D, Li R, Huang X. Identification of molecular correlations of RBM8A with autophagy in Alzheimer’s disease. Aging. 2019;11(23):11673–85. https://doi.org/10.18632/aging.102571.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  15. Langfelder P, Horvath S. WGCNA: an R package for weighted correlation network analysis. BMC Bioinformatics. 2008;9:559. https://doi.org/10.1186/1471-2105-9-559.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  16. de Souto MC, Jaskowiak PA, Costa IG. Impact of missing data imputation methods on gene expression clustering and classification. BMC Bioinformatics. 2015;16:64. https://doi.org/10.1186/s12859-015-0494-3.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  17. Anders S, Huber W. Differential expression analysis for sequence count data. Genome Biol. 2010;11(10):R106. https://doi.org/10.1186/gb-2010-11-10-r106.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  18. Algamal ZY, Alhamzawi R, Mohammad Ali HT. Gene selection for microarray gene expression classification using Bayesian Lasso quantile regression. Comput Biol Med. 2018;97:145–52. https://doi.org/10.1016/j.compbiomed.2018.04.018.

    CAS  Article  PubMed  Google Scholar 

  19. Woroniecka KI, Park AS, Mohtat D. Transcriptome analysis of human diabetic kidney disease. Diabetes. 2011;60(9):2354–69. https://doi.org/10.2337/db10-1181.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  20. Provaznik Z, Philipp A, Zeman F. Extracorporeal life support in postcardiotomy cardiogenic shock: a view on scenario, outcome, and risk factors in 261 patients. Thorac Cardiovasc Surg. 2021;69(3):271–8. https://doi.org/10.1055/s-0040-1714746.

    Article  PubMed  Google Scholar 

  21. Kerr JA, Liu RS, Gasser CE. Diet quality trajectories and cardiovascular phenotypes/metabolic syndrome risk by 11–12 years. Int J Obes. 2021. https://doi.org/10.1038/s41366-021-00800-x.

    Article  Google Scholar 

  22. Lay AC, Hale LJ, Stowell-Connolly H. IGFBP-1 expression is reduced in human type 2 diabetic glomeruli and modulates β1-integrin/FAK signalling in human podocytes. Diabetologia. 2021. https://doi.org/10.1007/s00125-021-05427-1.

    Article  PubMed  PubMed Central  Google Scholar 

  23. Filler G, Taheri S, McIntyre C. Chronic kidney disease stage affects small, dense low-density lipoprotein but not glycated low-density lipoprotein in younger chronic kidney disease patients: a cross-sectional study. Clin Kidney J. 2018;11(3):383–8. https://doi.org/10.1093/ckj/sfx115.

    CAS  Article  PubMed  Google Scholar 

  24. Binda V, Moroni G, Messa P. ANCA-associated vasculitis with renal involvement. J Nephrol. 2018;31(2):197–208. https://doi.org/10.1007/s40620-017-0412-z.

    CAS  Article  PubMed  Google Scholar 

  25. Lyons PA, Peters JE, Alberici F, Liley J, Coulson RMR, Astle W. Genome-wide association study of eosinophilic granulomatosis with polyangiitis reveals genomic loci stratified by ANCA status. Nat Commun. 2019;10(1):5120. https://doi.org/10.1038/s41467-019-12515-9.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  26. Qasim A, Patel J. ANCA Positive Vasculitis. In: StatPearls [Internet]. Treasure Island (FL): StatPearls Publishing; 2021. https://www.ncbi.nlm.nih.gov/books/NBK554372/. Accessed 5 Jan 2021.

  27. Li ZY, Ma TT, Chen M, Zhao MH. The prevalence and management of anti-neutrophil cytoplasmic antibody-associated vasculitis in China. Kidney Dis (Basel). 2016;1(4):216–23. https://doi.org/10.1159/000441912.

    Article  Google Scholar 

  28. von Borstel A, Sanders JS, Rutgers A. Cellular immune regulation in the pathogenesis of ANCA-associated vasculitides. Autoimmun Rev. 2018;17(4):413–21. https://doi.org/10.1016/j.autrev.2017.12.002.

    CAS  Article  Google Scholar 

  29. Konstantinov KN, Ulff-Møller CJ, Tzamaloukas AH. Infections and antineutrophil cytoplasmic antibodies: triggering mechanisms. Autoimmun Rev. 2015. https://doi.org/10.1016/j.autrev.2014.10.020.

    Article  PubMed  Google Scholar 

  30. Song G, Xu S, Zhang H. TIMP1 is a prognostic marker for the progression and metastasis of colon cancer through FAK-PI3K/AKT and MAPK pathway. J Exp Clin Cancer Res. 2016;35(1):148. https://doi.org/10.1186/s13046-016-0427-7.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  31. Batra J, Robinson J, Soares AS. Matrix metalloproteinase-10 (MMP-10) interaction with tissue inhibitors of metalloproteinases TIMP-1 and TIMP-2: binding studies and crystal structure. J Biol Chem. 2012;287(19):15935–46. https://doi.org/10.1074/jbc.M112.341156.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  32. Brew K, Nagase H. The tissue inhibitors of metalloproteinases (TIMPs): an ancient family with structural and functional diversity. Biochem Biophys Acta. 2010;1803(1):55–71. https://doi.org/10.1016/j.bbamcr.2010.01.003.

    CAS  Article  PubMed  Google Scholar 

  33. Grünwald B, Schoeps B, Krüger A. Recognizing the Molecular Multifunctionality and Interactome of TIMP-1. Trends Cell Biol. 2019;29(1):6–19. https://doi.org/10.1016/j.tcb.2018.08.006.

    CAS  Article  PubMed  Google Scholar 

  34. Lenz O, Elliot SJ, Stetler-Stevenson WG. Matrix metalloproteinases in renal development and disease. J Am Soc Nephrol. 2000;11(3):574–81.

    CAS  Article  PubMed  Google Scholar 

  35. Sanders JS, van Goor H, Hanemaaijer R. Renal expression of matrix metalloproteinases in human ANCA-associated glomerulonephritis. Nephrol Dialysis Transplant. 2004;19(6):1412–9. https://doi.org/10.1093/ndt/gfh186.

    CAS  Article  Google Scholar 

  36. Küster H, Thompson H, Kinet JP. Characterization and expression of the gene for the human Fc receptor gamma subunit. Definition of a new gene family. J Biol Chem. 1990;265(11):6448–52.

    Article  PubMed  Google Scholar 

  37. Shah S, Gibson AW, Ji C. Regulation of FcRγ function by site-specific serine phosphorylation. J Leukoc Biol. 2017;101(2):421–8. https://doi.org/10.1189/jlb.2AB0516-228R.

    CAS  Article  PubMed  Google Scholar 

  38. Sweet RA, Nickerson KM, Cullen JL. B cell-extrinsic Myd88 and Fcer1g negatively regulate autoreactive and normal B cell immune responses. J Immunol. 2017;199(3):885–93. https://doi.org/10.4049/jimmunol.1600861.

    CAS  Article  PubMed  Google Scholar 

  39. Le Coniat M, Kinet JP, Berger R. The human genes for the alpha and gamma subunits of the mast cell receptor for immunoglobulin E are located on human chromosome band 1q23. Immunogenetics. 1990;32(3):183–6. https://doi.org/10.1007/BF02114971.

    Article  PubMed  Google Scholar 

  40. Liang Y, Zhao M, Liang G. Construction of special reporter to detect DNA methylation regulatory activity in FCER1G gene promoter through patch-methylation. J Cent South University Med Sci. 2013;38(2):120–4. https://doi.org/10.3969/j.issn.1672-7347.2013.02.002.

    CAS  Article  Google Scholar 

  41. Walker ME, Hatfield JK, Brown MA. New insights into the role of mast cells in autoimmunity: evidence for a common mechanism of action? Biochim Biophys Acta. 2012. https://doi.org/10.1016/j.bbadis.2011.02.009.

    Article  PubMed  PubMed Central  Google Scholar 

  42. Grabauskas G, Wu X, Gao J. Prostaglandin E2, produced by mast cells in colon tissues from patients with irritable bowel syndrome, contributes to visceral hypersensitivity in mice. Gastroenterology. 2020;158(8):2195-2207.e6. https://doi.org/10.1053/j.gastro.2020.02.022.

    CAS  Article  PubMed  Google Scholar 

  43. Vinen CS, Turner DR, Oliveira DB. A central role for the mast cell in early phase vasculitis in the Brown Norway rat model of vasculitis: a histological study. Int J Exp Pathol. 2004;85(3):165–74. https://doi.org/10.1111/j.0959-9673.2004.00382.x.

    Article  PubMed  PubMed Central  Google Scholar 

  44. Gan PY, O’Sullivan KM, Ooi JD. Mast cell stabilization ameliorates autoimmune anti-myeloperoxidase glomerulonephritis. J Am Soc Nephrol. 2016;27(5):1321–33. https://doi.org/10.1681/ASN.2014090906.

    CAS  Article  PubMed  Google Scholar 

  45. Gan PY, Summers SA, Ooi JD. Mast cells contribute to peripheral tolerance and attenuate autoimmune vasculitis. J Am Soc Nephrol. 2012;23(12):1955–66. https://doi.org/10.1681/ASN.2012060572.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  46. Brunini F, Page TH, Gallieni M. The role of monocytes in ANCA-associated vasculitides. Autoimmun Rev. 2016;15(11):1046–53. https://doi.org/10.1016/j.autrev.2016.07.031.

    CAS  Article  PubMed  Google Scholar 

  47. Hilhorst M, Shirai T, Berry G. T cell-macrophage interactions and granuloma formation in vasculitis. Front Immunol. 2014;5:432. https://doi.org/10.3389/fimmu.2014.00432.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  48. Marinaki S, Neumann I, Kälsch AI. Abnormalities of CD4 T cell subpopulations in ANCA-associated vasculitis. Clin Exp Immunol. 2005;140(1):181–91. https://doi.org/10.1111/j.1365-2249.2005.02731.x.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  49. Abdulahad WH, Lamprecht P, Kallenberg CG. T-helper cells as new players in ANCA-associated vasculitides. Arthritis Res Ther. 2011;13(4):236. https://doi.org/10.1186/ar3362.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  50. Lilliebladh S, Johansson Å, Pettersson Å. Phenotypic characterization of circulating CD4+ T Cells in ANCA-associated vasculitis. J Immunol Res. 2018;2018:6984563. https://doi.org/10.1155/2018/6984563.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  51. Brogan P, Eleftheriou D. Vasculitis update: pathogenesis and biomarkers. Pediatr Nephrol. 2018;33(2):187–98. https://doi.org/10.1007/s00467-017-3597-4.

    Article  PubMed  Google Scholar 

  52. Fousert E, Toes R, Desai J. Neutrophil extracellular traps (NETs) take the central stage in driving autoimmune responses. Cells. 2020;9(4):915. https://doi.org/10.3390/cells9040915.

    CAS  Article  PubMed Central  Google Scholar 

  53. Zhao L, David MZ, Hyjek E. M2 macrophage infiltrates in the early stages of ANCA-associated pauci-immune necrotizing GN. CJASN. 2015;10(1):54–62. https://doi.org/10.2215/CJN.03230314.

    CAS  Article  PubMed  Google Scholar 

  54. Li J, Yu YF, Liu CH. Significance of M2 macrophages in glomerulonephritis with crescents. Pathol Res Pract. 2017;213(9):1215–20. https://doi.org/10.1016/j.prp.2017.04.011.

    CAS  Article  PubMed  Google Scholar 

  55. Ge C, Xu M, Qin Y. Fisetin supplementation prevents high fat diet-induced diabetic nephropathy by repressing insulin resistance and RIP3-regulated inflammation. Food Funct. 2019;10(5):2970–85. https://doi.org/10.1039/c8fo01653d.

    CAS  Article  PubMed  Google Scholar 

  56. Ren Q, Guo F, Tao S. Flavonoid fisetin alleviates kidney inflammation and apoptosis via inhibiting Src-mediated NF-κB p65 and MAPK signaling pathways in septic AKI mice. Biomed Pharmacother. 2020;122:109772. https://doi.org/10.1016/j.biopha.2019.109772.

    CAS  Article  PubMed  Google Scholar 

Download references

Acknowledgements

This work has benefited from GEO. The authors thank the GEO network for its generous sharing of large amounts of data.

Funding

This research received no external funding.

Author information

Authors and Affiliations

Authors

Contributions

FZ, YZ, and TL contributed to conceptualization; TL, FZ, and MX were involved in methodology. All the authors contributed to the implementation and analyses of the experiments. TL, YZ, FZ, and MX drafted the article, and all the authors contributed to its completion. All the authors read and approved the final manuscript.

Corresponding authors

Correspondence to Yun Zhou or Tong Liu.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare no conflict of interest.

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

Xia, M., Zhao, F., Zhang, Y. et al. Identification of diagnostic markers and immune cell infiltration characteristics in antineutrophil cytoplasmic antibody-associated vasculitis by weighted gene co-expression network analysis. Eur J Med Res 27, 37 (2022). https://doi.org/10.1186/s40001-022-00666-3

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s40001-022-00666-3

Keywords

  • Antineutrophil cytoplasmic antibody-associated vasculitis (AAV)
  • ANCA-associated glomerulonephritis
  • Hub gene
  • CIBERSORT
  • Immune infiltration