Studies have demonstrated that at least 20% of individuals infected with SARS-CoV-2 remain asymptomatic^1–4. Although most global efforts have focused on severe illness in COVID-19, examining asymptomatic infection provides a unique opportunity to consider early immunological features that promote rapid viral clearance. Here, postulating that variation in the human leukocyte antigen (HLA) loci may underly processes mediating asymptomatic infection, we enrolled 29,947 individuals, for whom high-resolution HLA genotyping data were available, in a smartphone-based study designed to track COVID-19 symptoms and outcomes. Our discovery cohort (n = 1,428) comprised unvaccinated individuals who reported a positive test result for SARS-CoV-2. We tested for association of five HLA loci with disease course and identified a strong association between HLA-B*15:01 and asymptomatic infection, observed in two independent cohorts. Suggesting that this genetic association is due to pre-existing T cell immunity, we show that T cells from pre-pandemic samples from individuals carrying HLA-B*15:01 were reactive to the immunodominant SARS-CoV-2 S-derived peptide NQKLIANQF. The majority of the reactive T cells displayed a memory phenotype, were highly polyfunctional and were cross-reactive to a peptide derived from seasonal coronaviruses. The crystal structure of HLA-B*15:01–peptide complexes demonstrates that the peptides NQKLIANQF and NQKLIANAF (from OC43-CoV and HKU1-CoV) share a similar ability to be stabilized and presented by HLA-B*15:01. Finally, we show that the structural similarity of the peptides underpins T cell cross-reactivity of high-affinity public T cell receptors, providing the molecular basis for HLA-B*15:01-mediated pre-existing immunity.

BACKGROUND: Vaccine-elicited adaptive immunity is a prerequisite for control of SARS-CoV-2 infection. Multiple sclerosis (MS) disease-modifying therapies (DMTs) differentially target humoral and cellular immunity. A comprehensive comparison of the effects of MS DMTs on SARS-CoV-2 vaccine–specific immunity is needed, including quantitative and functional B and T cell responses. METHODS: Spike-specific Ab and T cell responses were measured before and following SARS-CoV-2 vaccination in a cohort of 80 study participants, including healthy controls and patients with MS in 6 DMT groups: untreated and treated with glatiramer acetate (GA), dimethyl fumarate (DMF), natalizumab (NTZ), sphingosine-1-phosphate (S1P) receptor modulators, and anti-CD20 mAbs. Anti–spike-Ab responses were assessed by Luminex assay, VirScan, and pseudovirus neutralization. Spike-specific CD4 ⁺ and CD8 ⁺ T cell responses were characterized by activation-induced marker and cytokine expression and tetramer. RESULTS: Anti-spike IgG levels were similar between healthy control participants and patients with untreated MS and those receiving GA, DMF, or NTZ but were reduced in anti-CD20 mAb– and S1P-treated patients. Anti-spike seropositivity in anti-CD20 mAb–treated patients was correlated with CD19 ⁺ B cell levels and inversely correlated with cumulative treatment duration. Spike epitope reactivity and pseudovirus neutralization were reduced in anti-CD20 mAb– and S1P-treated patients. Spike-specific CD4 ⁺ and CD8 ⁺ T cell reactivity remained robust across all groups, except in S1P-treated patients, in whom postvaccine CD4 ⁺ T cell responses were attenuated. CONCLUSION: These findings from a large cohort of patients with MS exposed to a wide spectrum of MS immunotherapies have important implications for treatment-specific COVID-19 clinical guidelines. FUNDING: NIH grants 1K08NS107619, K08NS096117, R01AI159260, R01NS092835, R01AI131624, and R21NS108159; NMSS grants TA-1903-33713 and RG1701-26628; Westridge Foundation; Chan Zuckerberg Biohub; Maisin Foundation.

MHC class I expression levels influence the strength of immune responses and represent another variable in determining outcome to disease beyond peptide binding alone. Identification of the HLA loci that vary in allelic expression levels and delineating the mechanism responsible for expression variation may provide the opportunity to modify their expression therapeutically. We have examined the expression levels of allelic lineages at the HLA-A locus in a sample of 216 European Americans using a real-time polymerase chain reaction assay, which amplifies all HLA-A lineages specifically with equal efficiency, and observed a gradient of expression that associates with HLA-A allelic lineage (R = 0.6, P = 5 × 10⁻²⁵). DNA methylation of the HLA-A gene appears to contribute to the variation in HLA-A mRNA expression levels, as a significant inverse correlation was observed between HLA-A mRNA expression levels in untreated cells and the degree to which expression is increased after treatment of the cells with a DNA methyltransferase inhibitor (R = 0.6, P = 2.8 × 10⁻⁶). Further, deep-sequencing and immunoprecipitation assays revealed allelic lineage-specific methylation patterns within the HLA-A promoter region where increased DNA methylation levels correlated significantly with reduced HLA-A expression levels (R = 0.89, P = 3.7 × 10⁻⁹). These data demonstrate HLA-A allelic lineage-specific variation in expression levels, and DNA methylation as a likely factor in contributing to this variation.

Polymorphisms located within the MHC have been linked to many disease outcomes by mechanisms not yet fully understood in most cases. Variants located within untranslated regions of HLA genes are involved in allele-specific expression and may therefore underlie some of these disease associations. We determined sequences extending nearly 2 kb upstream of the transcription start site for 68 alleles from 57 major lineages of classical HLA class I genes. The nucleotide diversity within this promoter segment roughly follows that seen within the coding regions, with HLA-B showing the highest (∼1.9%), followed by HLA-A (∼1.8%), and HLA-C showing the lowest diversity (∼0.9%). Despite its greater diversity, HLA-B mRNA expression levels determined in 178 European Americans do not vary in an allele- or lineage-specific manner, unlike the differential expression levels of HLA-A or HLA-C reported previously. Close proximity of promoter sequences in phylogenetic trees is roughly reflected by similarity of expression pattern for most HLA-A and -C loci. Although promoter sequence divergence might impact promoter activity, we observed no clear link between the phylogenetic structures as represented by pairwise nucleotide differences in the promoter regions with estimated differences in mRNA expression levels for the classical class I loci. Further, no pair of class I loci showed coordinated expression levels, suggesting that distinct mechanisms across loci determine their expression level under nonstimulated conditions. These data serve as a foundation for more in-depth analysis of the functional consequences of promoter region variation within the classical HLA class I loci.

The KIR genes and their HLA class I ligands have thus far not been investigated in pemphigus foliaceus (PF) and related autoimmune diseases, such as pemphigus vulgaris. We genotyped 233 patients and 204 controls for KIR by PCR-SSP. HLA typing was performed by LABType SSO reagent kits. We estimated the odds ratio, 95% confidence interval and performed logistic regression analyses to test the hypothesis that KIR genes and their known ligands influence susceptibility to PF. We found significant negative association between activating genes and PF. The activating KIR genes may have an overlapping effect in the PF susceptibility and the presence of more than three activating genes was protective (OR = 0.49, p = 0.003). A strong protective association was found for higher ratios activating/inhibitory KIR (OR = 0.44, p = 0.001). KIR3DS1 and HLA-Bw4 were negatively associated to PF either isolated or combined, but higher significance was found for the presence of both together (OR = 0.34, p<10⁻³) suggesting that the activating function is the major factor to interfere in the PF pathogenesis. HLA-Bw4 (80I and 80T) was decreased in patients. There is evidence that HLA-Bw4(80T) may also be important as KIR3DS1 ligand, being the association of this pair (OR = 0.07, p = 0.001) stronger than KIR3DS1-Bw4(80I) (OR = 0.31, p = 0.002). Higher levels of activating KIR signals appeared protective to PF. The activating KIR genes have been commonly reported to increase the risk for autoimmunity, but particularities of endemic PF, like the well documented influence the environmental exposure in the pathogenesis of this disease, may be the reason why activated NK cells probably protect against pemphigus foliaceus.