Inclusion body myositis (IBM) is a progressive inflammatory-degenerative muscle disease of older individuals, with some patients producing anti-cytosolic 5′-nucleotidase 1A (NT5C1A, aka cN1A) antibodies. Human Leukocyte Antigens (HLA) is the highest genetic risk factor for developing IBM. In this study, we aimed to further define the contribution of HLA alleles to IBM and the production of anti-cN1A antibodies.

We HLA haplotyped a Western Australian cohort of 113 Caucasian IBM patients and 112 ethnically matched controls using Illumina next-generation sequencing. Allele frequency analysis and amino acid alignments were performed using the Genentech/MiDAS bioinformatics package. Allele frequencies were compared using Fisher's exact test. Age at onset analysis was performed using the ggstatsplot package. All analysis was carried out in RStudio version 1.4.1717.

Our findings validated the independent association of HLA-DRB1*03:01:01 with IBM and attributed the risk to an arginine residue in position 74 within the DRβ1 protein. Conversely, DRB4*01:01:01 and DQA1*01:02:01 were found to have protective effects; the carriers of DRB1*03:01:01 that did not possess these alleles had a fourteenfold increased risk of developing IBM over the general Caucasian population. Furthermore, patients with the abovementioned genotype developed symptoms on average five years earlier than patients without. We did not find any HLA associations with anti-cN1A antibody production.

High-resolution HLA sequencing more precisely characterised the alleles associated with IBM and defined a haplotype linked to earlier disease onset. Identification of the critical amino acid residue by advanced biostatistical analysis of immunogenetics data offers mechanistic insights and future directions into uncovering IBM aetiopathogenesis.