Summary
Background and objectives The HLA system plays a central role in the distinction between self antigens and non-self antigens. This study aimed to investigate the association between HLA genes and the outcomes of patients with ANCA-associated vasculitis (AAV).
Design, setting, participants, & measurements This study recruited 152 consecutive Chinese patients with AAV. The predictive value of the HLA alleles for renal outcome, response to treatment, and all-cause mortality were analyzed.
Results The proportion of patients with treatment failure was significantly higher in DRB1*0405-positive patients than in DRB1*0405-negative patients (41.7% versus 12.9%; P=0.008; corrected P=0.02). After adjusting for the other potential predictors, DRB1*0405 was still an independent predictor for the poor response to treatment (hazard ratio [HR], 5.91; 95% confidence interval [95% CI], 1.23–28.52; P=0.03). Renal survival was significantly worse in patients with DRB1*0405 than those without DRB1*0405 (P<0.001; corrected P<0.001). After adjusting for the other potential predictors, DRB1*0405 was still an independent predictor for ESRD (HR, 5.50; 95% CI, 2.18–13.88; P<0.001). The probability of all-cause mortality in patients with DPB1*0402 was significantly higher than those without DPB1*0402 (P=0.02; corrected P=0.04). After adjusting for the other potential predictors, DPB1*0402 was still an independent predictor for all-cause mortality (HR, 2.52; 95% CI, 1.21–5.28; P=0.01).
Conclusions In AAV patients, DRB1*0405 might be an independent risk factor for the poor response to treatment and the deterioration of renal function, whereas DPB1*0402 might be an independent risk factor for all-cause mortality.
Introduction
ANCA-associated vasculitis (AAV) comprises granulomatosis with polyangiitis (GPA), microscopic polyangiitis (MPA), Churg-Strauss syndrome, and renal-limited vasculitis. The etiology of AAV remains unclear. As discussed in a review by Chen et al., environmental factors and certain drugs were assumed to contribute to the development of AAV (1). In addition to the factors above, the genetic background of AAV has long been of concern because some familial cases of AAV have been reported (2–7).
Among the genetic factors, HLA genes are considered a set of suspected genes in the development of AAV (8–14). A number of studies have shown an association of HLA genes, including both class I and class II, with the disease susceptibility of AAV (8–14). There are five loci of MHC class II genes in HLA, including HLA-DR, HLA-DP, HLA-DQ HLA-DO, and HLA-DM. Among the HLA class II genes, HLA-DR, HLA-DP, and HLA-DQ alleles were found responsible for both GPA and MPA (9,10,14). However, to the best of our knowledge, no association has been found between AAV and the HLA-DO and HLA-DM alleles. Previous studies mainly focused on the association between the frequency of HLA genes and disease susceptibility. Whether there is any association between HLA genes and outcomes of AAV patients is far from clear. This study aimed to investigate the association of some HLA class II genes with the outcomes of patients with AAV, including the response to immunosuppressive treatment, renal outcome, and mortality.
Materials and Methods
Patients and Controls
This study recruited 152 consecutive Chinese patients with AAV, diagnosed at Peking University First Hospital from 1996 to 2009. All of the patients met the criteria of the Chapel Hill Consensus Conference definition of AAV (15). Patients with secondary vasculitis or with other comorbid renal diseases were excluded. This study complied with the Declaration of Helsinki and was approved by the ethics committee of our hospital. Informed consent was obtained from each patient.
Detection of Serum ANCA
ANCA tests were performed by both indirect immunofluorescence assay and antigen-specific ELISA. Standard indirect immunofluorescence assays were performed according to the manufacturer’s instructions (Euroimmun, Lübeck, Germany). Ethanol-fixed human PMNs were used to detect ANCA and monkey liver sections were used to exclude antinuclear antibodies. Cytoplasmic ANCA and perinuclear ANCA were distinguished according to staining patterns by two experienced technicians. Two highly purified known ANCA antigens, proteinase 3 (PR3) and myeloperoxidase (MPO), purified as previously reported, were used as solid-phase ligands in ELISA (16).
DNA Samples
Peripheral blood samples (10 ml) from patients with AAV and normal controls were collected in the presence of EDTA. Genomic DNA was obtained from peripheral blood leukocytes with a salting-out procedure (17).
Sequence-Based Typing
Typing of HLA-DRB1, HLA-DPB1, HLA-DQB1, and HLA-DQA1 was performed by bidirectional sequencing of exon 2 using SeCore sequencing kits (Invitrogen, Brown Deer, WI), as previously described (14,18).
Treatment
Treatment protocols were previously described (19–21). The induction therapy included corticosteroids in combination with cyclophosphamide. Oral prednisone was prescribed at an initial dosage of 1 mg/kg per day for 4–6 weeks, with reducing doses over time to 12.5–15 mg by 3 months. We administered 0.5–1.0 g/m2 of cyclophosphamide administered intravenously every month or orally daily (2 mg/kg per day). A 25% dose reduction of cyclophosphamide was made for those participants aged >65 years, those who developed leukocytopenia (<4000 cells/mm3), and those with renal insufficiency. Patients with ARF or pulmonary hemorrhage received pulses of intravenous methylprednisolone (7–15 mg/kg per day) before the above standard induction therapy. Patients with severe pulmonary hemorrhage received additional plasma exchanges. For maintenance therapy, daily oral azathioprine was given.
Definition of the Response to Immunosuppressive Treatment
The response to immunosuppressive treatment was defined as previously described (22). The vasculitic disease activity was measured by the Birmingham Vasculitis Activity Score (BVAS) (23). Remission was defined as the absence of disease activity attributable to active disease qualified by the need for ongoing stable maintenance immunosuppressive therapy (complete remission), or a 50% reduction of disease activity score and absence of new manifestations (partial remission). Those who failed to satisfy the above criteria for remission (complete or partial) were defined as having treatment failure.
Statistical Analyses
Chi-squared tests, one-way ANOVAs, and t tests were performed as appropriate. Kaplan–Meier survival analysis was used to assess patient and renal survival. Multivariate analyses of the treatment response were performed using the logistic regression model. Multivariate analyses of renal and patient survival were performed with the Cox regression model. Results were expressed as hazard ratios (HRs) with 95% confidence intervals (95% CIs). The following variables were assessed as potential predictors of treatment response, ESRD, and all-cause mortality: age, sex, ANCA specificities in ELISA, initial serum creatinine, HLA genes, plasma exchanges, pulses of intravenous methylprednisolone, and the way in which cyclophosphamide was used in the induction therapy (monthly intravenous dosage versus daily oral dosage). The Benjamini–Hochberg false discovery rate correction was applied to correct the P value (24,25).
Statistical power was estimated using Power and Sample Size Estimation Software (http://biostat.mc.vanderbilt.edu/PowerSampleSize).
Differences were considered significant if the P value was <0.05. Statistical analyses were performed using the SPSS statistical software package (version 11.0; SPSS Inc, Chicago, IL).
Results
Clinical Data
Among the 152 patients with AAV, 75 were male and 77 were female, with a median age of 61 years (range, 14–86) at diagnosis. According to the Chapel Hill Consensus Conference definition (15), 44 patients (28.9%) were classified as having GPA, 94 patients (61.8%) were classified as having MPA, and 14 patients (9.3%) were classified as having renal-limited vasculitis. Our results showed that 130 patients were MPO-ANCA positive, 19 were PR3-ANCA positive, and the remaining 3 were both MPO-ANCA and PR3-ANCA positive. Among the 107 patients with MPA, 105 (98.1%) were MPO-ANCA positive, and the other 2 (1.9%) were both MPO-ANCA and PR3-ANCA positive. Among the 44 patients with GPA, 25 (56.8%) were MPO-ANCA positive and 19 (43.2%) were PR3-ANCA positive. The BVAS level at diagnosis was 20 (interquartile range [IQR], 16–24].
HLA-DRB1, HLA-DPB1, HLA-DQB1, and HLA-DQA1 Alleles
A total of 56 HLA-DRB1, HLA-DPB1, HLA-DQB1, and HLA-DQA1 alleles were typed. The frequencies of these typed alleles are shown in Table 1.
Frequencies of all HLA-DRB1 and HLA-DPB1 alleles
Association between HLA Genes and the Response to Immunosuppressive Treatment
At the sixth month after initiation of the induction therapy, 58 of 152 patients (38.2%) achieved complete remission, 80 of 152 patients (52.6%) achieved partial remission, and 14 of 152 (9.2%) patients had treatment failure. Of the 12 patients who were DRB1*0405 positive, 5 (41.7%) had treatment failure and 7 (58.3%) achieved remission. Of the 140 patients who were DRB1*0405 negative, 122 (87.1%) achieved remission and 18 (12.9%) had treatment failure. The proportion of patients with treatment failure was significantly higher in the DRB1*0405 positive group than in the DRB1* 0405 negative group (41.7% versus 12.9%; chi-squared test=7.14; P=0.008; corrected P=0.02). Power calculations indicated that this strategy would provide >70% power to detect HLA loci with a hazard risk >1.5.
Table 2 provides the clinical data of patients with and without DRB1*0405. The clinical characteristics were roughly comparable between these two groups.
Comparison of clinical data between DRB1*0405 positive and DRB1*0405 negative, DPB1*0402 positive, and DPB1*0402 negative patients
To further investigate the role of HLA-DRB1*0405 in predicting the treatment response, variables including age, sex, ANCA specificities in ELISA, initial serum creatinine, HLA genes, plasma exchanges, pulses of intravenous methylprednisolone, and the way in which cyclophosphamide was used in the induction therapy, in addition to DRB1*0405, were assessed as potential predictors of the treatment response in the multivariate analysis (Table 3). After adjusting for the other potential predictors described above, DRB1*0405 was still an independent predictor for the poor treatment response (HR, 5.91; 95% CI, 1.23–28.52; P=0.03).
Predictors of the treatment response in the multivariate analyses
Association between HLA Genes and ESRD
The median follow-up time was 17 months (IQR, 1–47]. During follow-up, 32 patients progressed to ESRD. As shown in Figure 1, the renal survival was significantly worse in patients with DRB1*0405 than those without DRB1*0405 (P<0.001; corrected P<0.001). The cumulative proportion of renal survival at 1 year was significantly lower in patients with DRB1*0405 than those without DRB1*0405 (52.4% versus 89.8%; P<0.01). This was also the case for the cumulative proportions of renal survival at 2 years (52.4% versus 86.3%; P<0.01) and 3 years (37.4% versus 80.7%; P<0.01). There was no significant association between other HLA alleles and ESRD.
Renal survival of patients with and without DRB1*0405.
To further investigate the role of HLA-DRB1*0405 in predicting ESRD, variables including age, sex, ANCA specificities, initial serum creatinine, plasma exchanges, pulses of intravenous methylprednisolone, and the way in which cyclophosphamide was used in induction therapy, in addition to DRB1*0405, were assessed as potential predictors of ESRD in the multivariate analysis (Table 4). After adjusting for the other variables described above, DRB1*0405 was still an independent predictor for ESRD (HR, 5.50; 95% CI, 2.18–13.88; P<0.001).
Predictors of ESRD in the multivariate analyses
Association between HLA Genes and All-Cause Mortality
The median follow-up time was 24 months (IQR, 2–48). During follow-up, 34 of 152 (22.4%) patients died. As shown in Figure 2, the probability of all-cause mortality in patients with DPB1*0402 was significantly higher than those without DPB1*0402 (P=0.02; corrected P=0.04). The cumulative proportion of patient survival at 1 year was significantly lower in patients with DPB1*0402 than those without DPB1*0402 (79.6% versus 89.3%). This was also the case at 2 years (68.0% versus 85.0%) and 3 years (61.0% versus 83.4%). Causes of death of these patients are shown in Table 5. There was no significant association between other HLA genes and all-cause mortality.
Survival of patients with and without DPB1*0402.
Causes of death in patients with and without DPB1*0402
Table 2 shows the clinical data of patients with and without DPB1*0402. The clinical characteristics were roughly comparable between these two groups.
The clinical data of patients with and without DPB1*0402 are presented in Table 6. To further investigate the role of HLA-DPB1*0402 in predicting all-cause mortality, variables including age, sex, ANCA specificities, initial serum creatinine, plasma exchanges, pulses of intravenous methylprednisolone, and the way in which cyclophosphamide was used in induction therapy, in addition to DPB1*0402, were assessed as potential predictors of all-cause mortality in the multivariate analysis (Table 6). After adjusting for the other variables described above, DPB1*0402 was still an independent predictor for all-cause mortality (HR, 2.52; 95% CI, 1.21–5.28; P=0.01).
Predictors of all-cause mortality in the multivariate analysis
Discussion
This study investigated the association of some HLA class II genes, including HLA-DRB1, HLA-DPB1, HLA-DQA1, and HLA-DQB1, with the outcomes of patients with AAV. Of all of these HLA genes, DRB1*0405 might be an independent predictor for the deterioration of renal function. In our study, the treatment response was defined by BVAS, which reflects disease activity of multi-systems. Our study showed that DRB1*0405 was also associated with the treatment response. In the multivariate analysis, DRB1*0405 was still an independent predictor for the poor treatment response. This suggests that DRB1*0405 had a great effect not only on renal function progression but also on the systemic disease activity. The DRB1*0405 allele has been reported to be associated with various autoimmune disorders including rheumatoid arthritis and Vogt-Koyanagi-Harada disease (26–28). Gencik et al. showed that TNFα 2a/DRB1*04 individuals were found to be over-represented in a subgroup of patients with GPA in ESRD as a result of renal vasculitis. ANCA antigens are mainly located in azurophilic granules of polymorphonuclear granulocytes; after being primed by TNFα, there are more ANCA antigens expressed on the membrane of polymorphonuclear granulocytes (29,30). Over-representation of TNFα 2a/DRB1*04 might contribute to the severe courses of GPA (11). We speculated that the mechanism of HLA-DRB1*0405 contributing to the outcome of AAV is possibly through the over-representation of TNFα, resulting in poor prognosis of renal function.
HLA-DPB1*0402 has been reported to be associated with some other autoimmune diseases such as type 1 diabetes (31,32). The association with genes of the HLA class II region is explained by their role in the presentation of peptides derived from exogenous protein to CD4+ T helper cells. In the case of type 1 diabetes, this mechanism may result in an inappropriate T cell immune response against self-antigens on the pancreatic β cell. In this study, the probability of all-cause mortality in patients with DPB1*0402 was significantly higher than those without DPB1*0402, and DPB1*0402 was an independent predictor for all-cause mortality. To the best of our knowledge, this is the first report to show that a certain HLA allele is associated with mortality in AAV patients. The exact mechanism of such an association requires further study.
For the other HLA genes in this study, we conducted further investigations to calculate the power value in the Kaplan–Meier survival analysis. Some of the power values were <0.6 (Supplemental Material), which means there might be a possibility of false negative results for these alleles. Future studies with larger samples are needed. The limitation of our research was the relatively short duration of follow-up. Moreover, it is necessary to validate or replicate our association results using an independent cohort.
In conclusion, we found some associations between certain HLA genes and outcomes of AAV patients. DRB1*0405 might be an independent risk factor for the poor treatment response and the deterioration of renal function, whereas DPB1*0402 might be an independent risk factor for all-cause mortality.
Disclosures
None.
Acknowledgments
We are very grateful to Professor Xue-ying Li for assistance with medical statistics.
This study was supported Chinese 973 Project Grant 2012CB517702 and National Natural Science Fund Grants 30972733 and 81021004.
Footnotes
Published online ahead of print. Publication date available at www.cjasn.org.
This article contains supplemental material online at http://cjasn.asnjournals.org/lookup/suppl/doi:10.2215/CJN.13071211/-/DCSupplemental.
- Received December 21, 2011.
- Accepted April 25, 2012.
- Copyright © 2012 by the American Society of Nephrology