Genetic Variants in Antineutrophil Cytoplasmic Antibody-Associated Vasculitis: A Bayesian Approach and Systematic Review

A number of genome-wide association studies (GWASs) and meta-analyses of genetic variants have been performed in antineutrophil cytoplasmic antibody (ANCA)-associated vasculitis. We reinterpreted previous studies using false-positive report probability (FPRP) and Bayesian false discovery probability (BFDP). This study searched publications in PubMed and Excerpta Medica Database (EMBASE) up to February 2018. Identification of noteworthy associations were analyzed using FPRP and BFDP, and data (i.e., odds ratio (OR), 95% confidence interval (CI), p-value) related to significant associations were separately extracted. Using filtered gene variants, gene ontology (GO) enrichment analysis and protein–protein interaction (PPI) networks were performed. Overall, 241 articles were identified, and 7 were selected for analysis. Single nucleotide polymorphisms (SNPs) discovered by GWASs were shown to be noteworthy, whereas only 27% of significant results from meta-analyses of observational studies were noteworthy. Eighty-five percent of SNPs with borderline p-values (5.0 × 10−8 < p < 0.05) in GWASs were found to be noteworthy. No overlapping SNPs were found between PR3-ANCA and MPO-ANCA vasculitis. GO analysis revealed immune-related GO terms, including “antigen processing and presentation of peptide or polysaccharide antigen via major histocompatibility complex (MHC) class II”, “interferon-gamma-mediated (IFN-γ) signaling pathway”. By using FPRP and BFDP, network analysis of noteworthy genetic variants discovered genetic risk factors associated with the IFN-γ pathway as novel mechanisms potentially implicated in the complex pathogenesis of ANCA-associated vasculitis.


Introduction
Antineutrophil cytoplasmic antibody (ANCA)-associated vasculitis (AAV) is a group of autoimmune diseases characterized by the inflammation of small vessels. Clinical overlap among the different subtypes of AAV, granulomatosis with polyangiitis (GPA), microscopic polyangiitis (MPA) and eosinophilic granulomatosis with polyangiitis (EGPA) can result in diagnostic difficulties [1]. The clinical manifestations of the disease types can vary widely; most patients have signs of a lung, kidney, and/or ear, nose, and throat (ENT) involvement. ENT involvement is especially common in patients with GPA and EGPA [1]. A positive test for ANCA aids in the initial diagnosis of AAV. However, overlaps between disease phenotype and ANCA serotype limit the ability of the ANCA serotype to discriminate between the different disease phenotypes [2]. In general, a perinuclear ANCA (p-ANCA)-target antigen MPO-is found in 60-80% of MPA patients. On the other hand, cytoplasmic ANCA (c-ANCA)-target antigen proteinase 3 (PR3)-is frequently found in patients with severe GPA (approximately 90% of patients) and to a lesser degree in those with non-severe disease (approximately 50% of patients) [1].
Even though the exact pathogenesis of AAV is unclear, it is generally believed that AAV has a genetic background. Several genome-wide association studies (GWASs) and meta-analyses of observational studies employing single nucleotide polymorphisms (SNPs) revealed a number of genetic variants to be associated with AAV [3][4][5][6]. Up to date, there have been three GWASs and meta-analyses with replication cohorts: two by the US Vasculitis Clinical Research Consortium (VCRC) [4,6] and one by the European Vasculitis Genetic Consortium (EVGC) [3]. Moreover, Rahmattulla and co-workers performed a meta-analysis with inclusion of EVGC GWAS data and other meta-analyses of observational studies [5].
Because the prior probabilities of genetic associations are low, the number of false-positive associations that are generated by chance alone is high. Not accounting for these low probabilities in the statistical analysis leads to an increased likelihood of finding false-positive associations. Since concerns have been raised about true positives of the discovered genetic variants in AAV, the aim of this study was to investigate whether these genetic variants are false-positives or are truly associated with AAV by applying a Bayesian approach. Moreover, we discussed and re-analyzed the filtered data based on the integration of the available results for genetic variants.

Database Search and Selection
The eligible studies were selected according to the standardized reporting protocol of systematic reviews and meta-analyses PRISMA (Preferred Reporting Items for Systematic reviews and Meta-Analyses) checklist (Table S1) [7]. A systematic search was performed in PubMed and EMBASE to retrieve studies published before 20 February 2018. Antineutrophil cytoplasmic antibody (ANCA)-associated vasculitis, polymorphisms, meta, genetic, variant and genome-wide association study (GWAS) were the terms used in the systematic search. At the end of the systematic search, 241 articles were identified, 13 were screened in detail, and 7 were selected for analysis [4][5][6][8][9][10][11]. (Figure 1 and Table S2).  The inclusion criteria were (1) a genetic meta-analysis or GWAS providing information of odds ratios (OR), 95% confidence intervals (CI) in AAV (2) studies written in English, and (3) most recent larger meta-analyses (if smaller numbers have been analyzed before). The exclusion criteria were (1) studies not related to genetic polymorphisms or AAV, (2) articles not providing accurate data sets (i.e., review articles), and (3) older meta-analyses with overlapping genetic variants.
Data elements collected from meta-analyses included, when available, author, publication year, gene variant, single nucleotide polymorphism ID (rsID), genotype comparison, OR, 95% CI, minor allelic frequency, ethnicity of study population, number of cases and controls, publication bias, and heterogeneity. Results of various subgroup analyses (i.e., Caucasian vs. mixed-population) were also extracted.

Statistical Method
Data related to statistically significant associations (i.e., OR, 95% CI, p-value) and heterogeneity parameters (i.e., p-value and I 2 ) were separately extracted. Statistically significant associations were selected if the reported p-value was <5 × 10 −8 for results of meta-analyses in GWASs for both the discovery and replication cohorts, and <0.05 for meta-analysis of general observational studies (i.e., case-control studies).
To assess and identify noteworthy associations, false-positive report probability (FPRP) and Bayesian false discovery probability (BFDP) were applied [12,13]. FPRP is defined as "the probability of no true association between a gene variant and disease (null hypothesis)" for a statistically significant association [12], and detailed calculation is presented in the supplementary file (Equation S1). FPRP is calculated with the observed p-value for the association, the statistical power of the test, and the prior probability that a molecular association is real. In this review, we calculated FPRP at two levels of assumed prior probabilities, which were 10 −3 and 10 −6 . The interpretation for the noteworthiness of significant associations using FPRP is that the FPRP value of <0.2 is noteworthy [12].
BFDP is another Bayesian statistical method for detecting the true association between a gene variant and disease [13]. Main differences of BFDP compared to FPRP is that BFDP is independent from a statistical power and its approximation is based on a logistic regression model instead of a standard normal distribution [13]. The interpretation for the noteworthiness of significant associations using BFDP is that the BFDP value of <0.8 is noteworthy [13].

Outcomes
This study presents both FPRP and BFDP methods because the genetic epidemiologists and clinicians use BFDP less frequently than FPRP [12], although BFDP is a more recently developed method with a more substantial justification for its use [13]. By summarizing both results of FPRP and BFDP, we provide the readers options for interpreting noteworthiness.

Construction of PPI (Protein-Protein Interaction) Network
The STRING 9.1 network database is one of the largest databases of direct (physical) protein-protein interactions and indirect (functional) interactions constructed from various data sources including genomic context predictions, high-throughput experiments, co-expression, and known databases [14]. The STRING database covers 9.6 million proteins from more than 2031 organisms. In our study, we used STRING database to identify the PPIs associated with genes mapping to AAV, GPA, PR3-ANCA, and MPO-ANCA SNPs.

Computation of Noteworthy Variants
The re-analysis of genetic variants was mainly driven by two categories of studies: one meta-analysis of observational studies and two GWASs with replication cohorts (combined analysis). All statistically significant SNPs and variants (p < 0.05 for observational studies and p < 5.0 × 10 −8 for GWASs) reported in the meta-analyses were included in this study. In addition, variants found in the GWAS meta-analysis with 5.0 × 10 −8 < p < 0.05, which were rejected as insignificant in their interpretation, were re-analyzed. Whether a variant was noteworthy was determined based on satisfaction of the condition in at least one of the values (FPRP < 0.2 and BFDP < 0.8). The reported results were based on the various subtypes (GPA and MPA). In this study we followed the classification system used in the original studies.
In PR3-ANCA, 5/24 (21%) comparisons were noteworthy in FPRP while 10/24 (42%) were found noteworthy in BFDP. Variants marked with 'NA' in the tables were not assessable with FPRP due to the mathematical error of calculating the inverse of the cumulative distribution, but BFDP was still computable. All the noteworthy variants in PR3-ANCA with FPRP estimation were also noteworthy in BFDP estimation ( Table 2 and Table S4). No noteworthy comparison was observed in both MPA and MPO-ANCA vasculitis (Tables 1 and 2 and Tables S3 and S4).
Secondly, 22/28 (79%) SNPs in the meta-analyses of GWASs were noteworthy in total. Without exception, statistically significant SNPs discovered by GWAS meta-analyses (p < 5.0 × 10 −8 ) were also noteworthy in FPRP estimation at the prior probabilities of either 10 -3 or 10 -6 with the statistical power to detect an OR of 1.2 or 1.5 and in BFDP estimation at the probability of 10 −3 or 10 −6 except where mathematical errors did not allow FPRP. However, among the 14 SNPs with p-values ranging between 0.05 and 5.0 × 10 −8 , which were considered as statistically non-significant in GWAS meta-analyses, 8 SNPs were reported noteworthy in both FPRP and BFDP estimation ( Table 4 and Table S5).
We could re-analyze one GWAS of patients with GPA, while the remaining two GWASs were unable to be re-assessed due to the absence of 95% CI data. With our re-analyses, all the significant SNPs (p < 5.0 × 10 8 ) observed in the GWASs were noteworthy in FPRP and BFDP estimation. However, among the borderline SNPs rejected by the GWAS with p > 5.0 × 10 −8 , 53/62 (85%) of the SNPs were found noteworthy in our re-analysis (Table 4 and Table S6).

Gene Network Analysis
FPRP and BFDP computation reported all the SNPs with genome-wide significance discovered by GWASs as noteworthy, which may indicate that SNPs identified to be significant within GWASs are highly credible due to a conservative statistical standard and larger sample size. Based on this result, the SNPs' related protein-protein interaction (PPI) network of AAV subgroups and a Venn diagram of SNPs were constructed with the noteworthy genetic variants sorted by FPRP and BFDP and all the meaningful GWAS SNPs. The Venn diagram showed that no SNPs overlap between PR3-ANCA and MPO-ANCA vasculitis (Figure 2A and Figure S1). The PPI networks among genes mapping to PR3-ANCA and MPO-ANCA vasculitis-associated SNPs showed closer connections ( Figure 2B,C). Interestingly, the genes associated with interferon-gamma (IFN-γ)-mediated signaling pathway are highly enriched in both PR3-ANCA and MPO-ANCA vasculitis-associated SNPs (e.g., yellow nodes in the PPI networks). In addition, the MPA subtype was not included in the diagram, as no noteworthy SNPs have been reported in our analysis. The PPI network of AAV showed a variety of network connections. The major genes included in the network were major histocompatibility complex (MHC) class I and II genes (i.e., human leukocyte antigen (HLA)-A, -B, -C, -G, -F, -DPA1, -DPB1, -DQB1, -DRA). Other subtypes are shown in Figure 3 and Figure S2. Table 1. Meta-analysis results of observational studies; gene variants with statistical significance (p-value < 0.05), found to be noteworthy by false-positive report probability (FPRP) or Bayesian false discovery probability (BFDP) for each clinical diagnosis antineutrophil cytoplasmic antibody (ANCA)-associated vasculitis (AAV), granulomatosis with polyangiitis (GPA), and microscopic polyangiitis (MPA)).  FPRP: false-positive report probability; BFDP: Bayesian false discovery probability; OR: odds ratio; 95% CI: 95% confidence interval; No.: number; NA: not applicable.        Table S7. The yellow nodes represent the genes associated with the interferongamma-mediated signaling pathway.    Table S7. The yellow nodes represent the genes associated with the interferongamma-mediated signaling pathway.  Table S7. The yellow nodes represent the genes associated with the interferon-gamma-mediated signaling pathway.

Minor
To investigate disease pathways, gene ontology (GO) enrichment analysis was performed. Several GO terms enriched in the SNPs set were identified by GO enrichment analysis for each subtype. Immune-related GO terms, including "antigen processing and presentation of peptide or polysaccharide antigen via MHC class II" (GO:0002504) with a p-value of 4.4 × 10 −17 , "interferon-gamma-mediated (IFN-γ) signaling pathway" (GO:0060333) with a p-value of 2.9 × 10 −12 , and "T cell receptor signaling pathway" (GO:0050852) with a p-value of 5.5 × 10 −10 , were identified as the most prominently enriched AAV-related gene sets (Figure 4 and Figure S3)

Discussion
In this study we investigated whether the genetic variants previously found to be associated with AAV are truly associated with AAV or are false-positives. Former studies have demonstrated that the target antigens of AAV, MPO, and PR3 are better discriminators of the disease than the clinical phenotypes, MPA and GPA [2]. Meta-analyses and GWASs have confirmed a genetic background of AAV. Our analysis found no genetic overlap of MPO-ANCA and PR3-ANCA vasculitis. This further underlines the need to distinguish both subtypes in clinical studies and may lead to a re-classification of AAV based on the serotypes. Our analyses further highlight that AAV is not a multigenic disease with the implication of several genetic variants contributing to one's individual risk to develop vasculitis.
This work used FPRP and BFDP estimation to re-analyze genetic associations from the retrieved research papers. Discovered SNPs of each included study with p-values under 0.05 or either 5.0 × 10 −8 were computed by both methods. Nearly half of the significant SNPs reported in meta-analyses of the observational studies were not noteworthy, whereas those SNPs in GWASs and their metaanalyses were noteworthy with FPRP and BFDP. However, when we applied the Bayesian procedures to SNPs associated with borderline error probabilities (5.0 × 10 −8 < p-value < 0.05) observed in a GWAS including patients with GPA, FPRP and BFDP computation yielded a noteworthiness of 85% of these SNPs.
The exact mechanisms leading to AAV onset are unclear, but a genetic predisposition and selfintolerance to environmental exposures, such as infections, have been proposed [14]. A current pathogenetic model proposes the initiation of an inflammatory cascade with the release of PR3 (proinflammatory cytokine) by an unknown insult, further selecting PR3-specific B lymphocytes producing ANCA and PR3 specific T cells [14]. Subsequently, neutrophils are activated by ANCA priming, leading to the destruction of endothelial cells by firmly attached neutrophils on endothelial cells [15,16]. Our GO enrichment analysis showed that the processing of antigens via MHC class II, the IFN-γ-mediated pathway and the T cell receptor signaling pathway are the main biological processes disturbed in AAV. GO enrichment analysis identified a crucial role of the IFN-γ mediated signaling pathway and showed comparable involvement in both, MPO-ANCA and PR3-ANCA vasculitis. PPI analysis revealed a role of the interferon regulatory factor 5 (IRF5) in the model including patients with GPA and AAV. In systemic lupus erythematosus (SLE), IRF5 genetic variants are associated with an increased risk to develop the disease, alongside elevation of IRF5 expression and IFN production [17]. IFN-γ has been studied in two biomarker studies involving a cohort with AAV and one cohort of patients with GPA. Analysis of samples obtained from patients with AAV recruited to participate in the RAVE trial indicated significant higher levels of IFN-γ during active disease compared to remission (137 samples each) and to healthy controls (68 samples) [18].

1.
Several studies in ANCA-associated vasculitis have reported statistically significant genetic variants.

2.
Integration and further analysis of the data from several studies using FPRP, BFDP, and GO analysis identified novel biological pathways in ANCA-associated vasculitis, including a role of the IFN-γ pathway.

3.
Our results underline that there is a dichotomy between the respective target antigens, PR3 and MPO.

Discussion
In this study we investigated whether the genetic variants previously found to be associated with AAV are truly associated with AAV or are false-positives. Former studies have demonstrated that the target antigens of AAV, MPO, and PR3 are better discriminators of the disease than the clinical phenotypes, MPA and GPA [2]. Meta-analyses and GWASs have confirmed a genetic background of AAV. Our analysis found no genetic overlap of MPO-ANCA and PR3-ANCA vasculitis. This further underlines the need to distinguish both subtypes in clinical studies and may lead to a re-classification of AAV based on the serotypes. Our analyses further highlight that AAV is not a multigenic disease with the implication of several genetic variants contributing to one's individual risk to develop vasculitis.
This work used FPRP and BFDP estimation to re-analyze genetic associations from the retrieved research papers. Discovered SNPs of each included study with p-values under 0.05 or either 5.0 × 10 −8 were computed by both methods. Nearly half of the significant SNPs reported in meta-analyses of the observational studies were not noteworthy, whereas those SNPs in GWASs and their meta-analyses were noteworthy with FPRP and BFDP. However, when we applied the Bayesian procedures to SNPs associated with borderline error probabilities (5.0 × 10 −8 < p-value < 0.05) observed in a GWAS including patients with GPA, FPRP and BFDP computation yielded a noteworthiness of 85% of these SNPs.
The exact mechanisms leading to AAV onset are unclear, but a genetic predisposition and self-intolerance to environmental exposures, such as infections, have been proposed [14]. A current pathogenetic model proposes the initiation of an inflammatory cascade with the release of PR3 (pro-inflammatory cytokine) by an unknown insult, further selecting PR3-specific B lymphocytes producing ANCA and PR3 specific T cells [14]. Subsequently, neutrophils are activated by ANCA priming, leading to the destruction of endothelial cells by firmly attached neutrophils on endothelial cells [15,16]. Our GO enrichment analysis showed that the processing of antigens via MHC class II, the IFN-γ-mediated pathway and the T cell receptor signaling pathway are the main biological processes disturbed in AAV. GO enrichment analysis identified a crucial role of the IFN-γ mediated signaling pathway and showed comparable involvement in both, MPO-ANCA and PR3-ANCA vasculitis. PPI analysis revealed a role of the interferon regulatory factor 5 (IRF5) in the model including patients with GPA and AAV. In systemic lupus erythematosus (SLE), IRF5 genetic variants are associated with an increased risk to develop the disease, alongside elevation of IRF5 expression and IFN production [17]. IFN-γ has been studied in two biomarker studies involving a cohort with AAV and one cohort of patients with GPA. Analysis of samples obtained from patients with AAV recruited to participate in the RAVE trial indicated significant higher levels of IFN-γ during active disease compared to remission (137 samples each) and to healthy controls (68 samples) [18]. Furthermore, the expression of chemokine (C-X-C motif) ligand (CXCL10 or IFN-γ-induced protein 10) was tested in the same study. CXCL10 expression is dependent on IFN-γ and is implicated in several autoimmune disorders [19]. In patients with AAV, CXCL10 levels did not differ among patients with active disease or in remission, while CXCL10 levels were significantly lower in healthy controls. This argues for a sufficient reduction of IFN-γ once remission is achieved in AAV, but in contrast, the levels of CXCL10 remain elevated. More research is clearly needed to understand the regulation of the IFN-γ pathway in AAV.
IFN-γ increases MHC class II expression in antigen-presenting cells (B cells, macrophages, and dendritic cells) that are engaged by microbial invasions [20,21]. Furthermore, IFN-γ negatively regulates anti-inflammatory IL-10 production and induces pro-inflammatory IL-12 secretion [22]. IL-12 drives helper T(Th) cells to differentiate towards a Th1 phenotype and stimulates natural killer (NK) cells to produce IFN-γ. IFN-γ secreted from NK cells induces monocyte differentiation towards active macrophages in the local inflammatory site [23]. In addition, MPO stimulation leads to a dose-dependent production of IFN-γ by monocytes [24]. MPO induces the formation of neutrophil extracellular traps (NET), which is thought to be a crucial factor in the cascade of inflammation in AAV [25]. The IFN-γ pathway may play a role in the development of local inflammation. Alongside the immune system activation and NET formation as crucial mechanisms in the inflammatory processes, our data confirm that genetic predisposition and the activation of antigen-presenting cells with the selection of autoimmune T cells are important in the pathogenesis of AAV.
Our PPI network shows the mainframe structure of genetic interactions among HLA-associated genes and COL11A2 in GPA, PR3-ANCA vasculitis, and AAV itself. COL11A2 gene codes for the pro-α2 chain of collagen type XI. COL11A2 was demonstrated to be in linkage disequilibrium with the HLA-DP gene in an AAV GWAS [3,26]. Except MHC coding genes, other genetic risk factors associated with AAV were PTPN22, RXRB, CTLA-4, MICA, SERPINA1, PRTN3, CD226, TLR9, IRF5, NOTCH4, AGER, and CFB, which is consistent with previous findings (Figure 3). PTPN22 encodes protein tyrosine phosphatase in lymphoid tissues with an abnormal regulatory CD4 T-cell (Treg) function and increased neutrophil function reported in PTPN22 variant (rs2476601) [27]. RXRB is a gene encoding a family of retinoid X receptors, which form homodimers and heterodimers with retinoic acid, thyroid hormone, and vitamin D receptors [28]. There was a functional difference in the responsiveness to vitamin A and vitamin D in AAV patients, which could not be explained by RXRB polymorphisms alone [29]. CTLA-4 is a gene encoding an inhibitory surface protein on activated T cells that interacts with CD80 or CD86 and competes with CD28 (co-stimulatory molecule) [30]. Abatacept, a selective modulator of the CD80/86-CD28 costimulatory signal was successfully used in patients with non-severe GPA and is currently in a phase III trial (ABROGATE, NCT02108860), which suggests that the CTLA-4 protein might be of importance in AAV pathogenesis [31]. SERPINA1 and PRTN3 are also referred to as crucial genes in the pathogenesis of PR3-ANCA vasculitis [27]. PRTN3 is the gene encoding PR3 that exists on the surface or in the cytoplasm of neutrophils, whereas SERPINA1 encodes α-1 antitrypsin, the major inhibitory molecule of PR3. Previous findings suggest that PRTN3 genetic variants determine the predominant location of PR3 in neutrophils, where it may affect the activity of neutrophils [27]. No noteworthy variants were observed in MPA, but genetic variants related to the HLA-DR and HLA-DQ loci were found to be associated with MPO-ANCA vasculitis ( Figure S2).
FPRP was criticized with respect to its heuristic derivation of the formula, especially with the use of α and 1-β as the probabilities of observing values greater or less than the test quantities under a null and alternative hypothesis [32]. Even though a labored statistical derivation exists, the concept of "noteworthiness" originated from FPRP has influenced recent genetic studies [12]. Moreover, although direct comparison between BFDP and FPRP is not possible, the behavior of both is similar in promoting SNPs' rankings [13]. BFDP produced more noteworthy findings than FPRP because the latter derives smaller posterior null estimates due to the conditioning on tail areas unlike the reliance on point estimates provided by BFDP [13,33].
Our study has some limitations. First, different GWASs included in this study share the same cohorts. Therefore, slight overestimation of some sharing genes found in these studies may exist. However, no effect will be present in the gene network analysis because weighting of a gene corresponded to the number of significant SNPs in each gene. Second, we were not able to analyze variants with incomplete data (those without 95% CIs). Thus, it was not possible to confirm that all genome-wide significant SNPs satisfied the FPRP and BFDP thresholds. In fact, SNPs significant in GWASs were included in GO enrichment analysis, suggesting that these SNPs would be noteworthy due to previous results (Tables 4 and 4 and Tables S5 and S6). Furthermore, the low frequency of patients with MPA in the included studies may have led to the finding of no noteworthy SNPs in our re-analysis. Lastly, linkage disequilibrium and expression changes of SNPs were not considered in this analysis, perhaps treating SNPs and its properties too simply. Further research comprising quantitative expression of loci is needed.
Nevertheless, our study merged current results of genetic associations in AAV and its subtypes. In addition, the investigation of false-positive results in genetic research proves that researchers should pay careful attention when interpreting the positive results reported in previous papers. All GWAS SNPs associated with a borderline p-value are worth further examination with various statistical methods, as we found many of them to be noteworthy with the used Bayesian methods. Moreover, we highlight the importance of the IFN-γ pathway in the pathogenesis of AAV through GO enrichment analysis using GWAS and meta-analysis gene sets. This suggests that further research into the IFN-γ pathway in AAV may lead to the development of novel therapeutic approaches for this complex disease.
Supplementary Materials: The following are available online at http://www.mdpi.com/2077-0383/8/2/266/s1. Figure S1: The Venn diagram for all subtypes and AAV with the noteworthy SNPs from our FPRP and BFDP analyses. Figure S2: Protein-protein interaction network of GPA (A), MPO-ANCA vasculitis (B), PR3-ANCA vasculitis (C), and ANCA-associated vasculitis (D). Genes are weighted according to the number of the SNPs of each gene. The gene set was composed of GWAS discovered SNP genes and noteworthy SNP genes in the meta-analyses of observational studies and GWAS meta-analyses. Borderline SNPs (5.0 × 10 −8 < p < 0.05) from GWASs that are noteworthy were also included in this analysis. Figure S3. Gene ontology enrichment analysis of GPA (A), MPO-ANCA vasculitis (B), PR3-ANCA vasculitis (C), and ANCA-associated vasculitis (D). The gene set was composed of GWAS discovered SNP genes and noteworthy SNP genes in the meta-analyses of observational studies and GWAS meta-analyses. Borderline SNPs (5.0 × 10 −8 < p < 0.05) from GWASs that are noteworthy were also included in this analysis. Table S1: Preferred Reporting Items for Systematic reviews and Meta-Analyses (PRISMA) checklist. Table S2. List of publications included in this re-analysis. Table S3. Meta-analysis results of observational studies, gene variants with statistical significance (p-value < 0.05) for each clinical diagnosis (AAV, GPA, and MPA). Summary of genetic variants with FPRP and BFDP values, and noteworthy genetic variants are shown in bold. Table S4. Meta-analysis results of observational studies; gene variants with statistical significance (p-value < 0.05) for each serologic diagnosis (MPO-ANCA, PR3-ANCA). Summary of genetic variants with FPRP and BFDP values, and noteworthy genetic variants are shown in bold. Table S5. Results of meta analyses (combined analysis) with genome-wide association studies and replication cohort. Noteworthy genetic variants with satisfied FPRP and BFDP values are shown in bold. Table S6. Re-analysis of the SNPs discovered in genome-wide association studies of patients with GPA. Xie et al. only provided GPA SNPs with ORs and 95% CIs. Table S7: These gene sets were computed in gene ontology enrichment analysis and protein-protein interaction network analysis. This gene set include noteworthy variants from meta-analyses and statistically significant GWAS SNPs. Equation S1: Calculation of false positive report probability (FPRP) and Bayesian false discovery probability (BFDP).