key: cord-0913269-d4pluush authors: Deng, Xu; Gao, Junying; Zhao, Fei title: Identification of differentially expressed genes and pathways in kidney of ANCA-associated vasculitis by integrated bioinformatics analysis date: 2022-02-16 journal: Renal failure DOI: 10.1080/0886022x.2022.2030755 sha: 3364bd381641884d35422432ef2e17653643236f doc_id: 913269 cord_uid: d4pluush The antineutrophil cytoplasmic antibody (ANCA)-associated vasculitis (AAV) is a systematic of relatively rare autoimmune diseases with unknown cause. Kidney involvement is one of the most common clinical manifestations, and the degree of renal damage is closely associated with the development and prognosis of AAV. In this study, we utilized the Robust Rank Aggreg (RRA) method in R to integrate GSE104948, GSE104954, GSE108109, GSE108112, and GSE108113 profile datasets loaded from Gene Expression Omnibus (GEO) database and identified a set of differentially expressed genes (DEGs) in kidney between AAV patients and living donors. Then, the results of gene ontology (GO) functional annotation showed that immunity and metabolism involved process of AAV both in glomerulus and tubulointerstitial. The Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis revealed that following pathways, such as complement and coagulation cascades pathway; Staphylococcus aureus infection; disease-COVID-19; and systemic lupus erythematosus (SLE) pathway play a crucial role in AAV. Next, the results analyzed by protein–protein interaction (PPI) network and Cytoscape software exhibited the hub genes ALB, TYROBP, and CYBB existed in both glomerular and tubulointerstitial compartments datasets. Finally, KEGG analysis using genes of two most important modules also further validated complement and coagulation cascades pathway and S. aureus infection existed both in glomerulus and tubulointerstitial compartments datasets. In conclusion, this study identified key genes and pathways involved in kidney of AAV, which was benefit to further uncover the mechanisms underlying the development and progress of AAV, biomarkers, and potential therapeutic targets as well. The antineutrophil cytoplasmic antibody (ANCA)-associated vasculitis (AAV) is a series of relatively rare autoimmune diseases with etiology unknown, mainly injury small vessels [1] . AAV includes granulomatosis with polyangiitis (GPA), microscopic polyangiitis (MPA), and eosinophilic granulomatosis with polyangiitis (EGPA) [2] . The serological marker for AAV is ANCA [3] . There exist two major ANCA types, namely, directed against proteinase 3 (PR3) and myeloperoxidase (MPO). PR3-ANCA is mainly associated to GPA and MPO-ANCA to MPA. But substantial overlap exists and discordant ANCA type (e.g., MPO-positive GPA) is sometimes observed in certain populations, particularly in Asia [4] . There is no different influence of AAV between males and females [5] [6] [7] [8] . The prevalence of AAV is estimated from dozens to hundreds per million [9, 10] . Global figures of AAV prevalence display substantial geographic variations [11] . During the past decades, both the incidence and prevalence of AAV have grown probably ascribe to improvements in disease classification, ANCA testing, and clinical recognition [7, [12] [13] [14] . Although AAV affects small vessels in various organs, kidney involvement is one of the most common characteristics. Up to now, increasing studies point out that renal involvement is an important factor affecting the prognosis of patients with AAV both in adult and children. Before the employ of high-dose glucocorticoids (GC) and cyclophosphamide (CYC), the mortality rate of patients with severe AAV was up to 80% one year after diagnosis [15] . Recently, the estimated five-year survival is 74-91% and 45-76% for GPA and MPA, respectively [16] . However, infections accompanied by immunosuppressive therapy also attribute to rather high mortality rate [15] . Thus, further revealing the pathogenesis of AAV and finding out the biomarkers and new therapeutic targets have become urgent emergency. In this study, five microarray datasets download from the Gene Expression Omnibus (GEO) database were analyzed to further pick out the differentially expressed genes (DEGs) in AAV. After normalized all the datasets, the Robust Rank Aggreg (RRA) R package was applied to integrate DEGs. Subsequently, captured DEGs were analyzed by R package 'clusterprofiler' for the gene ontology (GO) functional annotation and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis. Then, a protein-protein interaction (PPI) network of the DEGs was constructed using search tool for the retrieval of interacting genes/proteins (STRING) database, hub genes, and modules were identified by plug-in of Cytoscape software. Our research for analysis of AAV DEGs datasets may serve as reliable molecular markers and effective therapeutic targets, and also help to understand the pathophysiological mechanisms in AAV. Microarray data and kidney tissue processing By searching in GEO database (http://www.ncbi.nlm.nih. gov/geo/), GSE104948, GSE104954, GSE108109, GSE108112 and GSE108113 profile datasets which met our requirements were downloaded. The datasets of GSE104948 and GSE104954 were based on a GeneChip Human Genome U133 Plus 2.0 Array platform (Affymetrix; Thermo Fisher Scientific, Inc., Waltham, MA), and datasets of GSE108109, GSE108112, and GSE108113 were based on a GeneChip Human Gene 2.1 ST Array platform (Affymetrix; Thermo Fisher Scientific, Inc., Waltham, MA). The samples of GSE104948 and GSE108109 came from glomerular compartments, samples of GSE104954 and GSE108112 from tubulointerstitial compartments. GSE108113 dataset was composed by GSE108109 and GSE108112. The detail of relevant datasets is shown in Table 1 . Human kidney tissue was treated as previously illustrated [17] . In brief, obtained kidney tissue was stored in RNAlater (ThermoFisher) and microdissected into tubulointerstitial and glomerular compartments [18] . Transcriptional data were employed to evaluate reliability of microdissection by reaching to 16-64 folds enrichment of glomerular selective or tubulointerstitial selective transcripts in each independent compartment [18] . The R package 'limma' was applied to normalize each expression matrix respectively, and then further screen DEGs in each dataset. These genes met the requirement, which were with p-value < 0.05 and logFC >1 were considered upregulated DEGs in each dataset. On the contrary, these genes, which were with p-value <0.05 and logFC <À1 were identified as downregulated DEGs in each dataset. The Venn diagram which exhibited the relationship among these DEGs was visualized applying VennDiagram R package. The captured DEGs were further integrated applying the RRA R package according to glomerular compartments DEGs, tubulointerstitial compartments DEGs, glomerular and tubulointerstitial compartments DEGs. Then, the R package 'VennDiagram' was employed to exhibit the relationship among these integrated DEGs. GO enrichment and KEGG pathway analysis was achieved by the R package 'clusterprofiler' [19] . Data with adjusted p < 0.05 were considered statistically significant. The online database, STRING (https://string-db.org/), was employed to get the interactions between different proteins [20] . Then the obtained data from STRING were further imported into Cytoscape software to acquire PPI network [21] . The cytoHubba plug-in of Cytoscape software was used to pick out the top 10 hub genes based on the value of degree. The MCODE plug-in of Cytoscape software was applied to analyze PPI network and acquire the top two modules. The ANCA-associated vasculitis datasets GSE104948, GSE104954, GSE108109, and GSE108112 were normalized, and the results are displayed in Supplementary Figure 1 . The DEGs were screened using the R package 'limma' (adjusted p < 0.05 and jlogFCj > 1). Firstly, 392 differential genes were contained in the GSE104948 dataset, including 165 upregulated and 227 downregulated genes. Next, there were 241 differential genes in the GSE104954 dataset, including 76 upregulated genes and 165 downregulated genes. In addition, 1519 differential genes, 941 upregulated genes, and 578 downregulated genes were contained in the GSE108109 dataset. Finally, the GSE108112 dataset contained 438 differential genes, including 271 upregulated genes and 167 downregulated genes. The DEGs of four datasets were exhibited in Figure 1 , and the top 75 genes' heatmap in each dataset is shown in Figure 2 . The integrated DEGs were screened utilizing the RRA package (corrected p < 0.05 and jlogFC j>1). About 143 integrated DEGs from glomerular compartments datasets (GSE104948 and GSE108109), consisting of 73 upregulated genes and 70 downregulated genes, were identified by the RRA method. What's more, 81 DEGs from integrated tubulointerstitial compartments datasets (GSE104954 and GSE108112) contained 22 upregulated genes and 59 downregulated genes. In addition, 273 integrated DEGs from glomerular and tubulointerstitial compartments datasets (GSE104948, GSE104954, GSE108109, GSE108112) were consisted by 146 upregulated genes and 127 downregulated genes. As shown in Figure 3 , the top 20 upregulated and downregulated genes in each integrated DEGs were showed on heatmap, respectively. Venn diagram of each integrated DEGs is exhibited in Supplementary Figure 3 (B). As it is well accepted now, the GO terms would be divided into three parts: cell component (CC), molecular function (MF), and biological process (BP). The GO results of each integrated DEGs are exhibited in Figure 4 and Supplementary Figure 4 . These significant genes from glomerular compartments were mainly enriched in neutrophil degranulation (gene ontology: BP), apical part of cell (gene ontology: CC), and organic acid blinding (gene ontology: MF). These significant genes from tubulointerstitial compartments were mainly enriched KEGG pathway annotation of the integrated DEGs was using the clusterProfile R package, and the results of the analysis are displayed in Figure 5 . The upregulated DEGs from glomerular compartments datasets (GSE104948 and GSE108109) were mainly enriched in these pathways, such as phenylalanine metabolism; tyrosine metabolism and glycine, serine, and threonine metabolism ( Figure 5(A) ). The downregulated DEGs from glomerular compartments datasets were mainly enriched in Staphylococcus aureus infection; complement and coagulation cascades; pertussis ( Figure 5(A) ). The upregulated DEGs from tubulointerstitial compartments datasets (GSE104954 and GSE108112) were enriched in renin-angiotensin system; parathyroid hormone synthesis, secretion, and action; MAPK signaling pathway, etc. (Figure 5(B) ). The downregulated DEGs from tubulointerstitial compartments datasets were mainly enriched in S. aureus infection; pertussis; complement and coagulation cascades ( Figure 5(B) ). The upregulated DEGs from glomerular and tubulointerstitial compartments datasets (GSE104948, GSE104954, GSE108109, GSE108112), were mainly enriched in these pathways, for example, phenylalanine metabolism; glycine, serine and threonine metabolism; glycolysis/gluconeogenesis ( Figure 5 (C)). The downregulated DEGs from glomerular and tubulointerstitial compartments datasets were enriched in pertussis; complement and coagulation cascades; viral protein interaction with cytokine and cytokine receptor, etc. ( Figure 5 (C)). The PPI networks from glomerular compartments DEGs, tubulointerstitial compartments DEGs, glomerulular and tubulointerstitial compartments DEGs were separately displayed in Figure 6 . The top 10 essential hub genes screened according value of degree are exhibited in Figure 7 . And the top two crucial modules are displayed in Figure 8 . The genes from top two modules of glomerular compartments mainly enriched in S. aureus infection; complement and coagulation cascades; phagosome. Genes from top two modules of tubulointerstitial compartments enriched in S. aureus infection; complement and coagulation cascades; pertussis. The genes from top two modules of glomerular and tubulointerstitial compartments enriched in following pathways, listing as chemokine signaling pathway; cytokine-cytokine receptor interaction; complement and coagulation cascades, etc. All these results further exhibited that infection and immunity played very vital role in AAV. The occurrence of AAV is relatively low and the disease is recognized as a comparatively rare autoimmune disease. However, it is not rare to nephrologists, rheumatologists, clinical immunologists, as well as other physicians [22] . AAV can attack people in all ages, including children, but the elderly at 65-74 years are more susceptible [23] . In spite of adequate immunosuppressive, the prognosis of AAV patient is still poor [24] . Previous researches revealed that the five-year survival rate of patient with AAV was about 70%, and infection was the main cause of death [25] . Recently, with increasing knowledge of AAV pathogenesis, the management and diagnosis of AAV have greatly improved; however, challenges still yet to be totally resolved. Therefore, further exploring the potential mechanisms related to the malignant biological process of AAV would help to detect reliable biomarkers for the early diagnosis, predicting prognosis, early monitoring recurrence, and new therapeutic intervention targets as well. AAV can undermine small vessels in various systems and organs, such as lungs, skin, joints, and kidneys [26] . Kidney involvement is one of the most common clinical manifestation of AAV that could happen over half of the AAV patients, and dialysis-dependency at disease onset is high [27, 28] . A cohort study included 502 patients with AAV diagnosed between 1985 and 2007 in America detected kidney involvement reached to 97%, ESKD was 32% and 27% died [29] . Consistently, a descriptive study recruited renal-limited AAV patients(N ¼ 16) between 2005 and 2016 from 8 centers worldwide, detected 38% progressing to ESRD, seven deaths, five within 6 months of diagnosis and two another at 16 and 21 months of follow-up [30] . In China, a review from Zhao et.al systematically summarized the progress of AAV in Chinese patients [31] . To some degree different from western countries, it showed AAV was not uncommon in China, the MPA was the compelling preponderance, MPO-ANCA was much more often than PR3-ANCA, and secondary infection (rather than AAV activity) was the main cause of death [31] . What's more, research in children of AAV in China found AAV usually occurred in teenage girls, the most commonly affected organ was also the kidney, and 58.8% of them progressed to ESRD with a median duration of three months [32] . Thus, exploring the pathophysiological mechanism of kidney in AAV will not only benefit to raise awareness occurrence and development of AAV, but also contribute to improve the treatment and prognosis. Considering function of glomerulus focused on filtration and tubular mainly on reabsorption, we analyzed datasets according to the sample collected from glomerular or/and tubulointerstitial, respectively. In this study, the GSE104948, GSE104954, GSE108109, GSE108112, and GSE108113 datasets were analyzed employing RRA method. The captured integrated DEGs were then experienced GO analysis. The genes from glomerular compartment were mainly enriched in neutrophil degranulation (gene ontology: BP), apical part of cell (gene ontology: CC), and organic acid binding (ontology: MF). The genes from tubulointerstitial compartment were mainly enriched in humoral immune response (gene ontology: BP), secretory granule lumen (gene ontology: CC), and endopeptidase activity (ontology: MF). And genes from glomerular and tubulointerstitial were mainly enriched in neutrophil degranulation (gene ontology: BP), collagen-containing (gene ontology: CC), and glycosaminoglycan binding (gene ontology: MF). These results suggested that immunity and metabolism involved in process of AAV both in glomerulus and tubulointerstitial. The KEGG pathway analysis revealed that the integrated upregulated genes in glomerular compartment were mainly enriched in the following five pathways: phenylalanine metabolism; tyrosine metabolism; glycine, serine and threonine metabolism; cholesterol metabolism; arginine and proline metabolism. The downregulated genes in glomerular compartment focused on S. aureus infection; complement and coagulation cascades; pertussis; neutrophil extracellular trap formation and phagosome. The integrated upregulated genes in tubulointerstitial compartment were enriched in these pathways, such as renin-angiotensin system; parathyroid hormone secretion and action; MAPK signaling pathway; colorectal cancer and PD-L1 expression and PD-1 checkpoint pathway in cancer and so on. The downregulated genes in tubulointerstitial compartment were mainly in following pathways: S. aureus infection; pertussis; complement and coagulation cascades; systemic lupus erythematosus; and phagosome, etc. The upregulated genes from glomerular and tubulointerstitial were enriched in phenylalanine metabolism; glycine, serine and threonine metabolism; glycolysis/gluconeogenesis, etc. The downregulated genes from glomerular and tubulointerstitial compartment focused on pertussis; complement and coagulation cascades; viral protein interaction with cytokine and cytokine receptor; S. aureus infection; amebiasis; ECM-receptor interaction, etc. It was important to discover complement and coagulation cascades pathway played a crucial role in AAV for the reason that it existed both in glomerular compartment and tubulointerstitial compartment datasets. The complement system, which included dozens of membrane-bound and circulating proteins, had a crucial role both in humoral and innate immunity. Activation of the complement system was mainly comprised by three pathways: the classical, the alternative pathways, and the lectin. Activation of any of three pathways resulted in formation of C3 convertase, and followed by activation of C5, finally resulted in initiating the assembly of the membrane attack complex (MAC), also called C5b-9 [33] . The complement system was used to being ignored by AAV pathophysiology [34] , for lacking indirect signs of complement consumption in plasma of patients with AAV by routine tests. Until recently, the significance of complement had attracted growing attention [35, 36] . A meta-analysis found plasma MAC, C5a, and factor B in patients with active AAV were increased compared with controls and patients in remission (excluding factor B) [37] . Another study also demonstrated AAV patients with hypo-complementemic and low serum C3 were independent predictors for ESRD [38, 39] . Among these patients with obtainable kidney biopsy, 56% had C3 deposition by immunofluorescence and exhibited more commonly a mixed histological type, low level of serum C3, serious infection [38] . Avacopan, a new oral competitive inhibitor that interfered the binding of C5a to C5aR1 (CD88), revealed to have least similar efficacy contrasted to high-dose corticosteroids in patients with renal involved active AAV [40] . S. aureus infection was also a crucial pathway in AAV. Study had found an increased nasal staphylococci carriage, which associated with relapse of Wegener's granulomatosis [41] . A systematic review exhibited median time between infection and vasculitis development was three months' cases of infection-induced MPO-ANCA, 21.7% patients expired during follow-up, 52.2% regressed after the resolution of infection in 12/23 [42] . It suggested that sufficient anti-infective therapy was very important in MPO-AAV before immunosuppressive treatment, or else might lead to catastrophic consequence [43] . Meanwhile, our data analysis showed coronavirus disease-COVID-19 and systemic lupus erythematosus (SLE) pathway also involved in AAV, which was consistent with previous reports. Up to now, growing studies had reported AAV after COVID-19 infection, which not only included adult [44] [45] [46] , but also children [47] . Meanwhile, some studies also reported there was an overlap between AAV and SLE [48, 49] . In light of these above reports, testing ANCA was vital for patients with COVID-19 infection or SLE. Then, we separately constructed PPI network with integrated genes from different renal part. Ten hub genes from glomerular compartment were ALB, ITGB2, TYROBP, CYBB, CD86, CD53, PLEK, C3AR1, FCER1G, FN1. The top 10 hub genes from tubulointerstitial compartment listed as ALB, EGF, C3, CYBB, IL10RA, TYROBP, TIMP1, CCL5, C1QB, CTSS. The top 10 hub genes from glomerular and tubulointerstitial were ALB, ITGB2, TLR2, PTPRC, FN1, TYROBP, CYBB, C3, KNG1, EGF. Obviously, the hub genes ALB, TYROBP, and CYBB existed in both the hub genes of glomerular or/and tubulointerstitial datasets. Highly consistent with our research, the research from Jae-Seung Moon et al. found that C-reactive protein (CRP) to serum albumin (ALB) ratio identified as an independent predictor of allcause of mortality in patients with AAV [50] . Thus, ALB might be as a biomarker in AAV. TYROBP was an adapter protein which non-covalently associated with activating receptors found on the surface of a variety of immune cells to mediate signaling and cell activation following ligand binding by the receptors [51] , mediate immune cell activation [52, 53] . Despite the relationship between TYROBP and AAV was not yet reported. TYROBP might be associate with pathophysiological process of AAV for the reason that immunity imbalance was an essential factor in AAV. CYBB was a critical component of the membrane-bound oxidase of phagocytes that generated superoxide, as it was the beta chain of Cytochrome b (-245). CYBB deficiency was one of five described biochemical defects associated with chronic granulomatous disease (CGD) [54] . But the relationship between CYBB and AAV was still not yet reported. Considering these, exploring the relationship between hub genes and AAV may help to elaborate the occurrence and development mechanism of AAV. Meanwhile, some hub genes only significantly different expressed in glomerulus or tubulointerstitial compartment dataset. Thus, it hinted the change expression of these genes may be beneficial to predict damage existed in glomerulus or tubulointerstitial. But it needed more deepen research. To identify tissue specific of hub genes, we searched hub genes expression of normal tissues in the GTEx database (https://www.gtexportal.org/). In hub genes from glomerular datasets, ALB expressed relatively high in kidney-cortex compared with other tissues (Supplementary Figure 5(A) ). In hub genes from tubulointerstitial datasets, ALB expressed relatively high in kidney-cortex compared with other tissues and EGF expressed relatively high in kidney-medulia (Supplementary Figure 5(B) ). In hub genes from glomerular and tubulointerstitial datasets, it exhibited ALB and KNG1 expressed relatively high in kidney-cortex, EGF and KNG1 expressed relatively high in kidney-medulia (Supplementary Figure 5(C) ). These results showed that ALB, EGF and KNG1 might be potential crucial genes participated in the pathophysiological process in kidney of AAV. There is increasing evidence that represented the importance of genetic contribution to AAV. A genomewide association study (GWAS) about AAV from European revealed both major-histocompatibility-complex (MHC) and non-MHC associated with AAV and also that GPA and MPO were genetically distinct [55] . It is reported that MPO-ANCA was associated with HLA-DQ, while PR3-ANCA was associated with HLA-DP, the genes encoding a1-antitrypsin (SERPINA1) and proteinase 3(PRTN3) [55] . In addition, eosinophilic granulomatosis with polyangiitis (EGPA) patients with MPO-ANCA positive was associated with HLA-DQ, while there was no association between HLA-DQ and ANCA-negative EGPA [56] . Variants at the GPA33 and IL5/IRF1 loci were only associated with the ANCA-negative EGPA patients [56] . In Chinese patients, studies had also suggested the relationship between AAV and different HLA loci. HLA-DRB1 Ã 1454 was less prevalent in AAV than normal controls, DRB1 Ã 1101 was more frequent in MPA patients, DRB1 Ã 1202 was more prevalent in PR3-ANCApositive granulomatosis patients than normal controls [57] . HLA-DPB1 variant rs3117242 was found to be associated with GPA in the Han Chinese population [58] . All of these variants might affect expression of certain genes. Due to downloaded datasets not including clinic details and GWAS results, these integrated DEGs were specific to AAV patients or not, and which variants of genes existed in these samples, were still yet to be known. To solve these puzzles, we have gathered samples and clinic information of AAV patients and other autoimmune diseases in our center, intending to unveil this tricky puzzle by combining methods, such as microarray analysis, GWAS, expression quantitative trait loci (eQTL) analysis and bidirectional Mendelian randomization study, etc. In summary, this research systematically explored potential mechanism underlying the kidney of AAV via integrated bioinformatics analysis. Meanwhile, this study also screened certain crucial genes and significant pathways which might benefit to hunt new biomarkers or therapeutic targets in kidney of AAV. Additionally, considering the different location and function of glomerular and tubulointerstitial in kidney, we analyzed glomerular and tubulointerstitial datasets separately. It would be more specifical to explore the pathophysiological mechanism of kidney in AAV. In the future, animal and clinic molecular biology experiments will be further required to validate the findings of this research. No potential conflict of interest was reported by the author(s) Antineutrophil cytoplasm antibody-associated vasculitis: recent developments ANCAassociated vasculitis -clinical utility of using ANCA specificity to classify patients Anti-neutrophil cytoplasmic autoantibodies with specificity for myeloperoxidase in patients with systemic vasculitis and idiopathic necrotizing and crescentic glomerulonephritis Characteristics of chinese patients with wegener's granulomatosis with antimyeloperoxidase autoantibodies Prevalence of wegener's granulomatosis, microscopic polyangiitis, polyarteritis nodosa and Churg-Strauss syndrome within a defined population in Southern Sweden Classification, epidemiology and clinical subgrouping of antineutrophil cytoplasmic antibody (ANCA)-associated vasculitis The epidemiology of antineutrophil cytoplasmic autoantibodyassociated vasculitis in Olmsted county, Minnesota: a twenty-year US population-based study The frequency of ANCAassociated vasculitis in a national database of hospitalized patients in China Epidemiology of primary systemic vasculitis in the Australian Capital Territory and South-Eastern New South Wales Giant cell arteritis is more prevalent in urban than in rural populations: results of an epidemiological study of primary systemic vasculitides in Germany Epidemiology of vasculitides: differences between Japan, Europe and North America Increasing incidence of Wegener's granulomatosis in Sweden Incidence of ANCA-associated vasculitis in a UK mixed ethnicity population An update on the epidemiology of ANCA-associated vasculitis Anti-neutrophil cytoplasmic antibody-associated vasculitis: prevalence, treatment, and outcomes Damage in the anca-associated vasculitides: long-term data from the European Vasculitis study group (EUVAS) therapeutic trials Quantitative gene expression analysis in renal biopsies: a novel protocol for a high-throughput multicenter application Metabolic pathways and immunometabolism in rare kidney diseases clusterProfiler: an R package for comparing biological themes among gene clusters The STRING database in 2021: customizable protein-protein networks, and functional characterization of useruploaded gene/measurement sets Cytoscape: a software environment for integrated models of biomolecular interaction networks ANCA-associated vasculitis Epidemiology of systemic vasculitis: a ten-year study in the United Kingdom Pathogenesis of rapidly progressive glomerulonephritis: what do we learn? Analyzing fatal cases of chinese patients with primary antineutrophil cytoplasmic antibodies-associated renal vasculitis: a 10-year retrospective study Outcome of ANCA-associated renal vasculitis: a 5-year retrospective study Clinical analysis of ANCA-associated renal vasculitis patients with chronic dialysis Classification of antineutrophil cytoplasmic autoantibody vasculitides: the role of antineutrophil cytoplasmic autoantibody specificity for myeloperoxidase or proteinase 3 in disease recognition and prognosis Patient outcomes in renal-limited antineutrophil cytoplasmic antibody vasculitis with inactive histology The prevalence and management of anti-neutrophil cytoplasmic antibody-associated vasculitis in China Clinicopathological analysis of 34 cases of primary antineutrophil cytoplasmic antibody-associated vasculitis in Chinese children Molecules great and small: the complement system With complements from ANCA mice Complement alternative pathway in ANCA-associated vasculitis: Two decades from bench to bedside The rise of complement in ANCAassociated vasculitis: from marginal player to target of modern therapy The alternative complement pathway in ANCA-associated vasculitis: further evidence and a meta-analysis Hypocomplementemia is associated with more severe renal disease and worse renal outcomes in patients with ANCA-associated vasculitis: a retrospective cohort study Clinical characteristics and prognostic risk factors of anti-neutrophil cytoplasmic antibody (ANCA)-associated vasculitides (AAV) Avacopan in the treatment of ANCA-associated vasculitis The relation between Staphylococcus aureus and Wegener's granulomatosis: current knowledge and future directions Infection-induced myeloperoxidase specific antineutrophil cytoplasmic antibody (MPO-ANCA) associated vasculitis: a systematic review Case 2-2020: a 64-year-old man with fever and respiratory failure ANCA-associated vasculitis after COVID-19 ANCA glomerulonephritis post moderna COVID-19 vaccination De novo ANCAassociated vasculitis with glomerulonephritis in COVID-19 Acute ANCA vasculitis and asymptomatic COVID-19 Systemic lupus erythematosus and antineutrophil cytoplasmic antibody-associated vasculitis overlap syndrome in a 77-year-old man: a case report A case of systemic lupus erythematosus/antineutrophil cytoplasmic antibody-associated vasculitis overlap syndrome with dissociated pathological and immunological findings C-Reactive protein to serum albumin ratio is an independent predictor of all-cause mortality in patients with ANCA-Associated vasculitis Cutting edge: signal-regulatory protein beta 1 is a DAP12-associated activating receptor expressed in myeloid cells DAP12 impacts trafficking and surface stability of killer immunoglobulin-like receptors on natural killer cells Cutting edge: inflammatory responses can be triggered by TREM-1, a novel receptor expressed on neutrophils and monocytes Germline CYBB mutations that selectively affect macrophages in kindreds with X-linked predisposition to tuberculous mycobacterial disease Genetically distinct subsets within ANCA-associated vasculitis Genome-wide association study of eosinophilic granulomatosis with polyangiitis reveals genomic loci stratified by ANCA status The association of HLA-DRB1 alleles with antineutrophil cytoplasmic antibody-associated systemic vasculitis in Chinese patients HLA-DPB1 variant rs3117242 is associated with anti-neutrophil cytoplasmic antibody-associated vasculitides in a Han Chinese population