MHC associations with clinical and autoantibody manifestations in European SLE

Abstract

Systemic lupus erythematosus (SLE) is a clinically heterogeneous disease affecting multiple organ systems and characterized by autoantibody formation to nuclear components. Although genetic variation within the major histocompatibility complex (MHC) is associated with SLE, its role in the development of clinical manifestations and autoantibody production is not well defined. We conducted a meta-analysis of four independent European SLE case collections for associations between SLE sub-phenotypes and MHC single-nucleotide polymorphism genotypes, human leukocyte antigen (HLA) alleles and variant HLA amino acids. Of the 11 American College of Rheumatology criteria and 7 autoantibody sub-phenotypes examined, anti-Ro/SSA and anti-La/SSB antibody subsets exhibited the highest number and most statistically significant associations. HLA-DRB1*03:01 was significantly associated with both sub-phenotypes. We found evidence of associations independent of MHC class II variants in the anti-Ro subset alone. Conditional analyses showed that anti-Ro and anti-La subsets are independently associated with HLA-DRB1*0301, and that the HLA-DRB1*03:01 association with SLE is largely but not completely driven by the association of this allele with these sub-phenotypes. Our results provide strong evidence for a multilevel risk model for HLA-DRB1*03:01 in SLE, where the association with anti-Ro and anti-La antibody-positive SLE is much stronger than SLE without these autoantibodies.

Figure 1

HLA-DRB1*03:01 dosage (average number of alleles observed) over levels of disease (a): (healthy controls/anti-Ro( −)/anti-Ro( +)); (b): (healthy controls/anti-La( −)/anti-La( +)); (c) (healthy controls/anti-Ro( −) AND anti-La( −)/anti-Ro( +) AND anti-La( +)/). Average dosage is represented by a square, whereas upper and lower 95% confidence intervals are represented by ‘ −’. Note that dosage ranges from 0 to 2 for each subject and so to convert to allele frequency you must divide by 2. All three plots have been truncated at 1.