HLA-DQ allele competition in narcolepsy: a comment on Tafti et al. DQB1 locus alone explains most of the risk and protection in narcolepsy with cataplexy in Europe

HLA-DQB1*05 subtypes and not DRB1*10:01 mediates risk in anti-IgLON5 disease

Anti-IgLON5 disease is a rare and likely underdiagnosed subtype of autoimmune encephalitis. The disease displays a heterogeneous phenotype that includes sleep, movement and bulbar-associated dysfunction. The presence of IgLON5-antibodies in CSF/serum, together with a strong association with HLA-DRB1*10:01∼DQB1*05:01, supports an autoimmune basis. In this study, a multicentric human leukocyte antigen (HLA) study of 87 anti-IgLON5 patients revealed a stronger association with HLA-DQ than HLA-DR. Specifically, we identified a predisposing rank-wise association with HLA-DQA1*01:05∼DQB1*05:01, HLA-DQA1*01:01∼DQB1*05:01 and HLA-DQA1*01:04∼DQB1*05:03 in 85% of patients. HLA sequences and binding cores for these three DQ heterodimers were similar, unlike those of linked DRB1 alleles, supporting a causal link to HLA-DQ. This association was further reflected in an increasingly later age of onset across each genotype group, with a delay of up to 11 years, while HLA-DQ-dosage dependent effects were also suggested by reduced risk in the presence of non-predisposing DQ1 alleles. The functional relevance of the observed HLA-DQ molecules was studied with competition binding assays. These proof-of-concept experiments revealed preferential binding of IgLON5 in a post-translationally modified, but not native, state to all three risk-associated HLA-DQ receptors. Further, a deamidated peptide from the Ig2-domain of IgLON5 activated T cells in two patients, compared with one control carrying HLA-DQA1*01:05∼DQB1*05:01. Taken together, these data support a HLA-DQ-mediated T-cell response to IgLON5 as a potentially key step in the initiation of autoimmunity in this disease.

Correlation analysis between HLA-DQA1*0102/DQB1*0602 genotypes and narcolepsy patients in China

Background and objective

At present, the etiology of narcolepsy is not fully understood, and it is generally believed to be an autoimmune reaction caused by interactions between environmental and genetic factors. Human leukocyte antigen (HLA) class II genes are strongly associated with this gene, especially HLA-DQB1*0602/DQA1*0102. In this study, we mainly analyzed the correlation between different genotypes of HLA-DQB1*0602/DQA1*0102 and clinical manifestations in Chinese patients with narcolepsy.

Experimental method

Narcolepsy patients who were treated at the Department of Neurology, The First Affiliated Hospital of Shandong First Medical University from January 2021 to September 2023 were selected. General information, sleep monitoring data, cerebrospinal fluid (CSF) orexin levels, and human leukocyte antigen gene typing data were collected. The statistical analysis was performed using SPSS 26.0, and the graphs were drawn using GraphPad Prism 9.5.

Main results

A total of 78 patients were included in this study. The DQA1 and DQB1 gene loci were detected in 54 patients, and only the DQB1 gene locus was detected in 24 narcoleptic patients. The most common allele at the HLA-DQB1 locus was *0602 (89.7%), and the most common genotype at this locus was *0602*0301 (19.2%), followed by *0602*0602 (17.9%). The most common phenotype of the HLA-DQA1 locus is *0102 (92.6%), and the most common genotype of this locus is *0102*0102 (27.8%), followed by *0102*0505 (14.8%). There were significant differences (p < 0.05) between HLA-DQB1*0602-positive and HLA-DQB1*0602-negative patients in terms of orexin-A levels, presence or absence of cataplexy, UNS, PSG sleep latency, REM sleep latency, N1 sleep percentage, oxygen depletion index, and average REM latency on the MSLT. The HLA-DQA1*0102-positive and HLA-DQA1*0102-negative patients showed significant differences (p < 0.05) in disease course, presence or absence of sudden onset, PSG REM sleep latency, N1 sleep percentage, and average REM latency on the MSLT. There were significant differences in the average REM latency of the MSLT between HLA-DQB1*0602/DQA1*0102 homozygous and heterozygous patients p < 0.05, and no differences were found in the baseline data, orexin-A levels, scale scores, or other sleep parameters.

Conclusion

Different genotypes of HLA-DQA1*0102/DQB1*0602 are associated with symptoms of cataplexy in Chinese narcoleptic patients. Homozygous individuals have a shorter mean REM latency in the MSLT, greater genetic susceptibility, and relatively more severe sleepiness.

The immunopathogenesis of narcolepsy type 1

Narcolepsy type 1 (NT1) is a chronic sleep disorder resulting from the loss of a small population of hypothalamic neurons that produce wake-promoting hypocretin (HCRT; also known as orexin) peptides. An immune-mediated pathology for NT1 has long been suspected given its exceptionally tight association with the MHC class II allele HLA-DQB1*06:02, as well as recent genetic evidence showing associations with polymorphisms of T cell receptor genes and other immune-relevant loci and the increased incidence of NT1 that has been observed after vaccination with the influenza vaccine Pandemrix. The search for both self-antigens and foreign antigens recognized by the pathogenic T cell response in NT1 is ongoing. Increased T cell reactivity against HCRT has been consistently reported in patients with NT1, but data demonstrating a primary role for T cells in neuronal destruction are currently lacking. Animal models are providing clues regarding the roles of autoreactive CD4+ and CD8+ T cells in the disease. Elucidation of the pathogenesis of NT1 will allow for the development of targeted immunotherapies at disease onset and could serve as a model for other immune-mediated neurological diseases.

Narcolepsy risk loci outline role of T cell autoimmunity and infectious triggers in narcolepsy

Narcolepsy type 1 (NT1) is caused by a loss of hypocretin/orexin transmission. Risk factors include pandemic 2009 H1N1 influenza A infection and immunization with Pandemrix®. Here, we dissect disease mechanisms and interactions with environmental triggers in a multi-ethnic sample of 6,073 cases and 84,856 controls. We fine-mapped GWAS signals within HLA (DQ0602, DQB1*03:01 and DPB1*04:02) and discovered seven novel associations (CD207, NAB1, IKZF4-ERBB3, CTSC, DENND1B, SIRPG, PRF1). Significant signals at TRA and DQB1*06:02 loci were found in 245 vaccination-related cases, who also shared polygenic risk. T cell receptor associations in NT1 modulated TRAJ*24, TRAJ*28 and TRBV*4-2 chain-usage. Partitioned heritability and immune cell enrichment analyses found genetic signals to be driven by dendritic and helper T cells. Lastly comorbidity analysis using data from FinnGen, suggests shared effects between NT1 and other autoimmune diseases. NT1 genetic variants shape autoimmunity and response to environmental triggers, including influenza A infection and immunization with Pandemrix®.

Characterization of T cell receptors reactive to HCRTNH2, pHA273-287, and NP17-31 in control and narcolepsy patients

Narcolepsy type 1 (NT1), a disorder caused by hypocretin/orexin (HCRT) cell loss, is associated with human leukocyte antigen (HLA)-DQ0602 (98%) and T cell receptor (TCR) polymorphisms. Increased CD4⁺ T cell reactivity to HCRT, especially DQ0602-presented amidated C-terminal HCRT (HCRT_NH2), has been reported, and homology with pHA_273-287 flu antigens from pandemic 2009 H1N1, an established trigger of the disease, suggests molecular mimicry. In this work, we extended DQ0602 tetramer and dextramer data to 77 cases and 44 controls, replicating our prior finding and testing 709 TCRs in Jurkat 76 T cells for functional activation. We found that fewer TCRs isolated with HCRT_NH2 (∼11%) versus pHA_273-287 or NP_17-31 antigens (∼50%) were activated by their ligand. Single-cell characterization did not reveal phenotype differences in influenza versus HCRT_NH2-reactive T cells, and analysis of TCR CDR3αβ sequences showed TCR clustering by responses to antigens but no cross-peptide class reactivity. Our results do not support the existence of molecular mimicry between HCRT and pHA_273-287 or NP_17-31.

HLA-DQ heterodimers in hematopoietic cell transplantation

HLA-DQ heterodimers increase the susceptibility to autoimmune diseases, but their role in hematopoietic cell transplantation is unknown. We tested the hypothesis that outcome after HLA-matched and HLA-DQ-mismatched hematopoietic cell transplantation is influenced by HLA-DQ heterodimers. Heterodimers were defined in 5164 HLA-matched and 520 HLA-DQ-mismatched patients and their transplant donors according to well-established crystallographic criteria. Group 1 (G1) heterodimers are any DQA1*02/03/04/05/06α paired with any DQB1*02/03/04β. Group 2 (G2) heterodimers are DQA1*01α paired with any DQB1*05/06β. Multivariable models identified significantly higher relapse risk in G1G2 and G2G2 compared with G1G1 HLA-matched patients with malignant disease; risk increased with an increasing number of G2 molecules. In HLA-DQ-mismatched transplantation for malignant diseases, matching or mismatching for G2 increased relapse risk. G2 lowered disease-free survival after both HLA-matched and HLA-DQ-mismatched transplantation. A paradigm based on HLA-DQ heterodimers provides a functional definition of the hematopoietic cell transplantation barrier and a means to lower risks for future patients.

Human Leukocyte Antigen Association Study Reveals DRB1*04:02 Effects Additional to DRB1*07:01 in Anti-LGI1 Encephalitis

BACKGROUND AND OBJECTIVES

To study human leukocyte antigen (HLA) allele associations in anti-leucine-rich glioma-inactivated 1 (LGI1) encephalitis.

METHODS

A multiethnic cohort of 269 patients with anti-LGI1 encephalitis and 1,359 controls was included. Four-digit HLA sequencing and genome wide association single-nucleotide polymorphism typing imputation (0.99 concordance) were used for HLA typing. Significance of primary and secondary associations was tested using χ2, Fisher exact tests, or logistic regression with the control of population stratification covariates when applicable.

RESULTS

DRB1*07:01 and DQA1*02:01, 2 alleles in strong linkage disequilibrium, were associated with the disease (90% vs 24%, OR = 27.8, p < 10e-50) across ethnicity independent of variation at DRB3 and DQB1, 2 flanking HLA loci. DRB1*07:01 homozygosity was associated with a doubling of risk (OR = 2.1, p = 0.010), suggesting causality. DRB1*07:01 negative subjects were younger (p = 0.003) and more frequently female (p = 0.015). Three patients with malignant thymomas did not carry DRB1*07:01, whereas patients with other tumors had high DRB1*07:01 frequency, suggesting that the presence of tumors other than thymomas may be coincidental and not causal. In both DRB1*07:01 heterozygous individuals and DRB1*07:01 negative subjects, DRB1*04:02 was associated with anti-LGI1 encephalitis, indicating an independent effect of this allele (OR = 6.85, p = 4.57 × 10-6 and OR = 8.93, p = 2.50 × 10-3, respectively). DRB1*04:02 was also independently associated with younger age at onset (β = -6.68, p = 9.78 × 10-3). Major histocompatibility complex peptide-binding predictions using LGI1-derived peptides revealed divergent binding propensities for DRB1*04:02 and DRB1*07:01 alleles, suggesting independent pathogenic mechanisms.

DISCUSSION

In addition to the established primary DRB1*07:01 association in anti-LGI1 encephalitis, we observe a secondary effect of DRB1*04:02 with lower age at onset. Our study provides evidence for secondary effects within HLA locus that correlate with clinical phenotypes in anti-LGI1 encephalitis.

Association of Human Leukocyte Antigens Class II Variants with Susceptibility to Hidradenitis Suppurativa in a Caucasian Spanish Population

Hidradenitis suppurativa (HS) is a chronic inflammatory cutaneous disease of the hair follicle typically presenting recurrent, painful, and inflamed lesions on the inverse areas of the body. Although its pathogenesis remains unknown, the immune system appears to play a potential role. To date, two previous studies have not found any association between the Human Leukocyte Antigen system (HLA) and HS. In this study we analyzed the HLA-A, -B, -C; and DRB1, -DQA1, and –DQB1 allele distribution in 106 HS patients and 262 healthy controls from a Caucasian population in Cantabria (northern Spain). HLA-A*29 and B*50 were significantly more common in HS patients and A*30 and B*37 in controls, but these associations disappeared after statistical correction. DRB1*07, DQA1*02, and DQB1*02 were significantly more common in controls (p 0.026, p 0.0012, and p 0.0005, respectively) and the HLA allele DQB1*03:01 was significantly more common in HS patients (p 0.00007) after the Bonferroni correction. The DRB1*07~DQA1*02~DQB1*02 haplotype was significantly more common in controls (p < 0.0005). This is the first study showing an association between HLA-class II and HS. Our results suggest that HLA-II alleles (DRB1*07, DQA1*02, DQB1*02, and DQB1*03:01) and the DRB1*07~DQA1*02~DQB1*02 haplotype could influence resistance or susceptibility to HS.

Clinical utility of next generation sequencing based HLA typing for disease association and pharmacogenetic testing

HLA associations have been linked to many diseases and are important for risk assessment of drug hypersensitivity reactions. The increasing number of HLA alleles discovered generated a list of ambiguities that cannot be resolved with the current clinical assays, which commonly include sequence-specific oligonucleotide probe (SSOP) genotyping, and real-time PCR with melting curve analysis. HLA typing by next-generation sequencing (NGS) has recently been adopted by clinical laboratories for transplantation testing, as it provides unambiguous and cost-effective HLA typing. The goal of this study was to evaluate the feasibility of using NGS-based HLA-B and DQ genotyping for clinical HLA disease association testing, and provide direct comparison with the currently used clinical tests, including SSOP genotyping, and real-time PCR with melting curve analysis. While the real-time PCR method is easy and inexpensive to perform, ambiguities are rapidly increasing as more and more HLA alleles are discovered. SSOP genotyping identifies the alleles present but limitations include ambiguities and underreporting less common alleles. Our data show that HLA typing by NGS is superior to the existing clinical methods for identifying HLA alleles associated with disease or drug hypersensitivity, and offers a viable approach for high volume clinical diagnostic laboratories.

[Narcolepsy: From the discovery of a wake promoting peptide to autoimmune T cell biology and molecular mimicry with flu epitopes]

Narcolepsy-cataplexy was first described in the late 19th century in Germany and France. Prevalence was established to be 0.05 % and a canine model was discovered in the 1970s. In 1983, a Japanese study found that all patients carried HLA-DR2, suggesting autoimmunity as the cause of the disease. Studies in the canine model established that dopaminergic stimulation underlies anti-narcoleptic action of psychostimulants, while antidepressants were found to suppress cataplexy through adrenergic reuptake inhibition. No HLA association was found in canines. A linkage study initiated in 1988 revealed in hypocretin (orexin) receptor two mutations as the cause of canine narcolepsy in 1999. In 1992, studies on African Americans showed that DQ0602 was a better marker than DR2 across all ethnic groups. In 2000, hypocretin-1/orexin A levels were measured in the cerebrospinal fluid (CSF) and found to be undetectable in most patients, establishing hypocretin deficiency as the cause of narcolepsy. Decreased CSF hypocretin-1 was then found to be secondary to the loss of the 70,000 neurons producing hypocretin in the hypothalamus, suggesting immune destruction of these cells as the cause of the disease. Additional genetic studies, notably genome wide associations (GWAS), found multiple genetic predisposing factors for narcolepsy. These were almost all involved in other autoimmune diseases, although a strong and unique association with T cell receptor (TCR) alpha and beta loci were observed. Nonetheless, all attempts to demonstrate presence of autoantibodies against hypocretin cells in narcolepsy failed, and the presumed autoimmune cause remained unproven. In 2009, association with strep throat infections were found, and narcolepsy onsets were found to occur more frequently in spring and summer, suggesting upper away infections as triggers. Following reports that narcolepsy cases were triggered by vaccinations and infections against influenza A 2009 pH1N1, a new pandemic strain that erupted in 2009, molecular mimicry with influenza A virus was suggested in 2010. This hypothesis was later confirmed by peptide screening showing higher activity of CD4⁺ T cell reactivity to a specific post-translationally amidated segment of hypocretin (HCRT-_NH2) and cross-reactivity of specific TCRs with a pH1N1-specific segment of hemagglutinin that shares homology with HCRT-_NH2. Strikingly, the most frequent TCR recognizing these antigens was found to carry sequences containing TRAJ24 or TRVB4-2, segments modulated by narcolepsy-associated genetic polymorphisms. Cross-reactive CD4⁺ T cells with these cross-reactive TCRs likely subsequently recruit CD8⁺ T cells that are then involved in hypocretin cell destruction. Additional flu mimics are also likely to be discovered since narcolepsy existed prior to 2009. The work that has been conducted over the years on narcolepsy offers a unique perspective on the conduct of research on the etiopathogeny of a specific disease.

HLA high-resolution typing by next-generation sequencing in Pandemrix-induced narcolepsy

The incidence of narcolepsy type 1 (NT1) increased in Sweden following the 2009–2010 mass-vaccination with the influenza Pandemrix-vaccine. NT1 has been associated with Human leukocyte antigen (HLA) DQB1*06:02 but full high-resolution HLA-typing of all loci in vaccine-induced NT1 remains to be done. Therefore, here we performed HLA typing by sequencing HLA-DRB3, DRB4, DRB5, DRB1, DQA1, DQB1, DPA1 and DPB1 in 31 vaccine-associated NT1 patients and 66 of their first-degree relatives (FDR), and compared these data to 636 Swedish general population controls (GP). Previously reported disease-related alleles in the HLA-DRB5*01:01:01-DRB1*15:01:01-DQA1*01:02:01-DQB1*06:02:01extended haplotype were increased in NT1 patients (34/62 haplotypes, 54.8%) compared to GP (194/1272 haplotypes, 15.3%, p = 6.17E^-16). Indeed, this extended haplotype was found in 30/31 patients (96.8%) and 178/636 GP (28.0%). In total, 15 alleles, four extended haplotypes, and six genotypes were found to be increased or decreased in frequency among NT1 patients compared to GP. Among subjects with the HLA-DRB5*01:01:01-DRB1*15:01:01-DQA1*01:02-DQB1*06:02 haplotype, a second DRB4*01:03:01-DRB1*04:01:01-DQA1*03:02//*03:03:01-DQB1*03:01:01 haplotype (p = 2.02E^-2), but not homozygosity for DRB1*15:01:01-DQB1*06:02:01 (p = 7.49E^-1) conferred association to NT1. Alleles with increased frequency in DQA1*01:02:01 (p = 1.07E^-2) and DQA1*03:02//*03:03:01 (p = 3.26E^-2), as well as with decreased frequency in DRB3*01:01:02 (p = 8.09E^-3), DRB1*03:01:01 (p = 1.40E^-2), and DQB1*02:01:01 (p = 1.40E^-2) were found among patients compared to their FDR. High-resolution HLA sequencing in Pandemrix-associated NT1 confirmed the strong association with the DQB1*06:02:01-containing haplotype but also revealed an increased association to the not previously reported extended HLA-DRB4*01:03:01-DRB1*04:01:01-DQA1*03:02//*03:03:01-DQB1*03:01:01 haplotype. High-resolution HLA typing should prove useful in dissecting the immunological mechanisms of vaccination-associated NT1.

HLA-DQB1 and HLA-DRB1 Variants Confer Susceptibility to Latent Autoimmune Diabetes in Adults: Relative Predispositional Effects among Allele Groups

Latent autoimmune diabetes in adults (LADA) was recently demonstrated to be the most frequent form of adult-onset autoimmune diabetes mellitus. Case–control studies have investigated the relationship between human leukocyte antigen (HLA)-DQB1 and HLA-DRB1 polymorphisms and LADA risk, but their conclusions are inconsistent. This study aimed to more precisely explore the correlation between these HLA gene variants and LADA development. Eight databases, including PubMed, Embase, and Medline, were systematically searched for relevant studies up to September 15, 2018. We performed this retrospective study using meta-analysis and relative predispositional effect (RPE) methods. The meta-analysis results indicated that DQB1*02 (odds ratio (OR) = 1.685, p_c < 0.005) and DQB1*06 (OR = 0.604, p_c = 0.010) have opposite effects on susceptibility to LADA, while a significant decrease in LADA risk caused by DQB1*05 (OR = 0.764, p_c = 0.100) disappeared upon Bonferroni correction. The RPE method confirmed the roles of DQB1*02 (χ² = 46.475, p < 0.001) and DQB1*06 (χ² = 17.883, p < 0.001) and further suggested protective effects of DQB1*05 (χ² = 16.496, p < 0.001). Additionally, the meta-analysis results showed that DRB1*03 (OR = 2.685, p_c < 0.013), DRB1*04 (OR = 1.954, p_c < 0.013), and DRB1*09 (OR = 1.346, p_c < 0.013) are associated with increased LADA risk, while DRB1*12 (OR = 0.600, p_c < 0.013) and DRB1*13 (OR = 0.583, p_c < 0.013) carriers have a decreased risk of developing LADA. Furthermore, the RPE method revealed that DRB1*03 (χ² = 98.754, p < 0.001), DRB1*04 (χ² = 94.685, p < 0.001), DRB1*09 (χ² = 40.489, p < 0.001), DRB1*01 (χ² = 12.181, p < 0.001), DRB1*07 (χ² = 10.882, p = 0.001), and DRB1*08 (χ² = 5.000, p = 0.025) play protective roles against LADA. LADA showed a close relationship with genetic polymorphisms of HLA-DQB1 and WHLA-DRB1, which could contribute to a better understanding of disease pathogenesis and the identification of predisposing loci in the diagnosis and treatment of LADA.

Pleiotropic genetic effects influencing sleep and neurological disorders

Research evidence increasingly points to the large impact of sleep disturbances on public health. Many aspects of sleep are heritable and genes influencing traits such as timing, EEG characteristics, sleep duration, and response to sleep loss have been identified. Notably, large-scale genome-wide analyses have implicated numerous genes with small effects on sleep timing. Additionally, there has been considerable progress in the identification of genes influencing risk for some neurological sleep disorders. For restless legs syndrome, implicated variants are typically in genes associated with neuronal development. By contrast, genes conferring risk for narcolepsy function in the immune system. Many genetic variants associated with sleep disorders are also implicated in neurological disorders in which sleep abnormalities are common; for example, variation in genes involved in synaptic homoeostasis are implicated in autism spectrum disorder and sleep-wake control. Further investigation into pleiotropic roles of genes influencing both sleep and neurological disorders could lead to new treatment strategies for a variety of sleep disturbances.

Use of Drosophila in the investigation of sleep disorders

Genetic underpinnings for sleep disorders in humans remain poorly identified, investigated and understood. This is due to the inherent complexity of sleep and a disruption of normal sleep parameters in a number of neurological disorders. On the other hand, there have been steady and remarkable developments in the investigation of sleep using model organisms such as Drosophila. These studies have illuminated conserved genetic pathways, neural circuits and intra-cellular signaling modules in the regulation of sleep. Additionally, work in model systems is beginning to clarify the role of the circadian clock and basal sleep need in this process. There have also been initial efforts to directly model sleep disorders in flies in a few instances where a genetic basis has been suspected. Here, we discuss the opportunities and limitations of studying sleep disorders in Drosophila and propose that a greater convergence of basic sleep research in model organisms and human genetics should catalyze better understanding of sleep disorders and generate viable therapeutic options.

HLA-DPB1 and HLA class I confer risk of and protection from narcolepsy

Type 1 narcolepsy, a disorder caused by a lack of hypocretin (orexin), is so strongly associated with human leukocyte antigen (HLA) class II HLA-DQA1(∗)01:02-DQB1(∗)06:02 (DQ0602) that very few non-DQ0602 cases have been reported. A known triggering factor for narcolepsy is pandemic 2009 influenza H1N1, suggesting autoimmunity triggered by upper-airway infections. Additional effects of other HLA-DQ alleles have been reported consistently across multiple ethnic groups. Using over 3,000 case and 10,000 control individuals of European and Chinese background, we examined the effects of other HLA loci. After careful matching of HLA-DR and HLA-DQ in case and control individuals, we found strong protective effects of HLA-DPA1(∗)01:03-DPB1(∗)04:02 (DP0402; odds ratio [OR] = 0.51 [0.38-0.67], p = 1.01 × 10(-6)) and HLA-DPA1(∗)01:03-DPB1(∗)04:01 (DP0401; OR = 0.61 [0.47-0.80], p = 2.07 × 10(-4)) and predisposing effects of HLA-DPB1(∗)05:01 in Asians (OR = 1.76 [1.34-2.31], p = 4.71 × 10(-05)). Similar effects were found by conditional analysis controlling for HLA-DR and HLA-DQ with DP0402 (OR = 0.45 [0.38-0.55] p = 8.99 × 10(-17)) and DP0501 (OR = 1.38 [1.18-1.61], p = 7.11 × 10(-5)). HLA-class-II-independent associations with HLA-A(∗)11:01 (OR = 1.32 [1.13-1.54], p = 4.92 × 10(-4)), HLA-B(∗)35:03 (OR = 1.96 [1.41-2.70], p = 5.14 × 10(-5)), and HLA-B(∗)51:01 (OR = 1.49 [1.25-1.78], p = 1.09 × 10(-5)) were also seen across ethnic groups in the HLA class I region. These effects might reflect modulation of autoimmunity or indirect effects of HLA class I and HLA-DP alleles on response to viral infections such as that of influenza.