Genome‐Wide Association Analysis Reveals Genetic Heterogeneity of Sjögren's Syndrome According to Ancestry

OBJECTIVE: The Sjögren's International Collaborative Clinical Alliance (SICCA) is an international data registry and biorepository derived from a multisite observational study of participants in whom genotyping was performed on the Omni2.5M platform and who had undergone deep phenotyping using common protocol‐directed methods. The aim of this study was to examine the genetic etiology of Sjögren's syndrome (SS) across ancestry and disease subsets. METHODS: We performed genome‐wide association study analyses using SICCA subjects and external controls obtained from dbGaP data sets, one using all participants (1,405 cases, 1,622 SICCA controls, and 3,125 external controls), one using European participants (585, 966, and 580, respectively), and one using Asian participants (460, 224, and 901, respectively) with ancestry adjustments via principal components analyses. We also investigated whether subphenotype distributions differ by ethnicity, and whether this contributes to the heterogeneity of genetic associations. RESULTS: We observed significant associations in established regions of the major histocompatibility complex (MHC), IRF5, and STAT4 (P = 3 × 10(−42), P = 3 × 10(−14), and P = 9 × 10(−10), respectively), and several novel suggestive regions (those with 2 or more associations at P < 1 × 10(−5)). Two regions have been previously implicated in autoimmune disease: KLRG1 (P = 6 × 10(−7) [Asian cluster]) and SH2D2A (P = 2 × 10(−6) [all participants]). We observed striking differences between the associations in Europeans and Asians, with high heterogeneity especially in the MHC; representative single‐nucleotide polymorphisms from established and suggestive regions had highly significant differences in the allele frequencies in the study populations. We showed that SSA/SSB autoantibody production and the labial salivary gland focus score criteria were associated with the first worldwide principal component, indicative of higher non‐European ancestry (P = 4 × 10(−15) and P = 4 × 10(−5), respectively), but that subphenotype differences did not explain most of the ancestry differences in genetic associations. CONCLUSION: Genetic associations with SS differ markedly according to ancestry; however, this is not explained by differences in subphenotypes.