Comprehensive analysis of the major histocompatibility complex in systemic sclerosis identifies differential HLA associations by clinical and serological subtypes

Abstract

Objective: The greatest genetic effect reported for systemic sclerosis (SSc) lies in the major histocompatibility complex (MHC) locus. Leveraging the largest SSc genome-wide association study, we aimed to fine-map this region to identify novel human leucocyte antigen (HLA) genetic variants associated with SSc susceptibility and its main clinical and serological subtypes.

Methods: 9095 patients with SSc and 17 584 controls genome-wide genotyped were used to impute and test single-nucleotide polymorphisms (SNPs) across the MHC, classical HLA alleles and their composite amino acid residues. Additionally, patients were stratified according to their clinical and serological status, namely, limited cutaneous systemic sclerosis (lcSSc), diffuse cutaneous systemic sclerosis (dcSSc), anticentromere (ACA), antitopoisomerase (ATA) and anti-RNApolIII autoantibodies (ARA).

Results: Sequential conditional analyses showed nine SNPs, nine classical alleles and seven amino acids that modelled the observed associations with SSc. This confirmed previously reported associations with HLA-DRB1*11:04 and HLA-DPB1*13:01, and revealed a novel association of HLA-B*08:01. Stratified analyses showed specific associations of HLA-DQA1*02:01 with lcSSc, and an exclusive association of HLA-DQA1*05:01 with dcSSc. Similarly, private associations were detected in HLA-DRB1*08:01 and confirmed the previously reported association of HLA-DRB1*07:01 with ACA-positive patients, as opposed to the HLA-DPA1*02:01 and HLA-DQB1*03:01 alleles associated with ATA presentation.

Conclusions: This study confirms the contribution of HLA class II and reveals a novel association of HLA class I with SSc, suggesting novel pathways of disease pathogenesis. Furthermore, we describe specific HLA associations with SSc clinical and serological subtypes that could serve as biomarkers of disease severity and progression.

Keywords: autoantibodies; genetic; immune complex diseases; polymorphism; systemic sclerosis.

Figure 1

Association signals for systemic sclerosis in the human leucocyte antigen region. The −log₁₀ of the meta-analysis p values are plotted against their chromosomal position. The red line represents the genome-wide level of significance (p value=5×10⁻⁰⁸). The size of the diamond indicates the degree of linkage disequilibrium with the strongest association from the meta-analysis (rs1048372).

Figure 2

LD among the independent variants. Circos plot depicting the LD relationship among the SNPs, four-digit classical HLA alleles and HLA amino acid residues independently associated from the sequential conditional analysis. HLA, human leucocyte antigen; LD, linkage disequilibrium; SNP, single-nucleotide polymorphism.