Characterization of the genetic determinants of context-specific DNA methylation in primary monocytes.

To better understand inter-individual variation in sensitivity of DNA methylation (DNAm) to immune activity, we characterized effects of inflammatory stimuli on primary monocyte DNAm (n = 190). We find that monocyte DNAm is site-dependently sensitive to lipopolysaccharide (LPS), with LPS-induced demethylation occurring following hydroxymethylation. We identify 7,359 high-confidence immune-modulated CpGs (imCpGs) that differ in genomic localization and transcription factor usage according to whether they represent a gain or loss in DNAm. Demethylated imCpGs are profoundly enriched for enhancers and colocalize to genes enriched for disease associations, especially cancer. DNAm is age associated, and we find that 24-h LPS exposure triggers approximately 6 months of gain in epigenetic age, directly linking epigenetic aging with innate immune activity. By integrating LPS-induced changes in DNAm with genetic variation, we identify 234 imCpGs under local genetic control. Exploring shared causal loci between LPS-induced DNAm responses and human disease traits highlights examples of disease-associated loci that modulate imCpG formation.

Cross-population applications of genomics to understand the risk of multifactorial traits involving inflammation and immunity.

The interplay between genetic and environmental factors plays a significant role in interindividual variation in immune and inflammatory responses. The availability of high-throughput low-cost genotyping and next-generation sequencing has revolutionized our ability to identify human genetic variation and understand how this varies within and between populations, and the relationship with disease. In this review, we explore the potential of genomics for patient benefit, specifically in the diagnosis, prognosis and treatment of inflammatory and immune-related diseases. We summarize the knowledge arising from genetic and functional genomic approaches, and the opportunity for personalized medicine. The review covers applications in infectious diseases, rare immunodeficiencies and autoimmune diseases, illustrating advances in diagnosis and understanding risk including use of polygenic risk scores. We further explore the application for patient stratification and drug target prioritization. The review highlights a key challenge to the field arising from the lack of sufficient representation of genetically diverse populations in genomic studies. This currently limits the clinical utility of genetic-based diagnostic and risk-based applications in non-Caucasian populations. We highlight current genome projects, initiatives and biobanks from diverse populations and how this is being used to improve healthcare globally by improving our understanding of genetic susceptibility to diseases and regional pathogens such as malaria and tuberculosis. Future directions and opportunities for personalized medicine and wider application of genomics in health care are described, for the benefit of individual patients and populations worldwide.

Patient stratification using plasma cytokines and their regulators in sepsis: relationship to outcomes, treatment effect and leucocyte transcriptomic subphenotypes.

RATIONALE: Heterogeneity of the host response within sepsis, acute respiratory distress syndrome (ARDS) and more widely critical illness, limits discovery and targeting of immunomodulatory therapies. Clustering approaches using clinical and circulating biomarkers have defined hyper-inflammatory and hypo-inflammatory subphenotypes in ARDS associated with differential treatment response. It is unknown if similar subphenotypes exist in sepsis populations where leucocyte transcriptomic-defined subphenotypes have been reported. OBJECTIVES: We investigated whether inflammatory clusters based on cytokine protein abundance were seen in sepsis, and the relationships with previously described transcriptomic subphenotypes. METHODS: Hierarchical cluster and latent class analysis were applied to an observational study (UK Genomic Advances in Sepsis (GAinS)) (n=124 patients) and two clinical trial datasets (VANISH, n=155 and LeoPARDS, n=484) in which the plasma protein abundance of 65, 21, 11 circulating cytokines, cytokine receptors and regulators were quantified. Clinical features, outcomes, response to trial treatments and assignment to transcriptomic subphenotypes were compared between inflammatory clusters. MEASUREMENTS AND MAIN RESULTS: We identified two (UK GAinS, VANISH) or three (LeoPARDS) inflammatory clusters. A group with high levels of pro-inflammatory and anti-inflammatory cytokines was seen that was associated with worse organ dysfunction and survival. No interaction between inflammatory clusters and trial treatment response was found. We found variable overlap of inflammatory clusters and leucocyte transcriptomic subphenotypes. CONCLUSIONS: These findings demonstrate that differences in response at the level of cytokine biology show clustering related to severity, but not treatment response, and may provide complementary information to transcriptomic sepsis subphenotypes. TRIAL REGISTRATION NUMBER: ISRCTN20769191, ISRCTN12776039.

Reframing sepsis immunobiology for translation: towards informative subtyping and targeted immunomodulatory therapies.

Sepsis is a common and deadly condition. Within the current model of sepsis immunobiology, the framing of dysregulated host immune responses into proinflammatory and immunosuppressive responses for the testing of novel treatments has not resulted in successful immunomodulatory therapies. Thus, the recent focus has been to parse observable heterogeneity into subtypes of sepsis to enable personalised immunomodulation. In this Personal View, we highlight that many fundamental immunological concepts such as resistance, disease tolerance, resilience, resolution, and repair are not incorporated into the current sepsis immunobiology model. The focus for addressing heterogeneity in sepsis should be broadened beyond subtyping to encompass the identification of deterministic molecular networks or dominant mechanisms. We explicitly reframe the dysregulated host immune responses in sepsis as altered homoeostasis with pathological disruption of immune-driven resistance, disease tolerance, resilience, and resolution mechanisms. Our proposal highlights opportunities to identify novel treatment targets and could enable successful immunomodulation in the future.

Optimized protocol for CRISPR knockout of human iPSC-derived macrophages.

Here, we present a protocol for lentiviral delivery of CRISPR-Cas9 to human induced pluripotent stem cell (iPSC)-derived macrophages using co-incubation with VPX virus-like particles (VPX-VLPs). We describe steps for producing polybrene and puromycin kill curves, VPX viral production, and VPX-VLP titration by western blotting. We then detail procedures for iPSC macrophage precursor lentiviral transduction and lentiviral CRISPR-Cas9-based knockout in iPSC-derived macrophages. This protocol uses efficient genome-editing techniques to explore macrophage involvement in immune response, chronic inflammation, neurodegenerative disease, and cancer progression. For complete details on the use and execution of this protocol, please refer to Navarro-Guerrero et al.1.

A common NFKB1 variant detected through antibody analysis in UK Biobank predicts risk of infection and allergy.

Infectious agents contribute significantly to the global burden of diseases through both acute infection and their chronic sequelae. We leveraged the UK Biobank to identify genetic loci that influence humoral immune response to multiple infections. From 45 genome-wide association studies in 9,611 participants from UK Biobank, we identified NFKB1 as a locus associated with quantitative antibody responses to multiple pathogens, including those from the herpes, retro-, and polyoma-virus families. An insertion-deletion variant thought to affect NFKB1 expression (rs28362491), was mapped as the likely causal variant and could play a key role in regulation of the immune response. Using 121 infection- and inflammation-related traits in 487,297 UK Biobank participants, we show that the deletion allele was associated with an increased risk of infection from diverse pathogens but had a protective effect against allergic disease. We propose that altered expression of NFKB1, as a result of the deletion, modulates hematopoietic pathways and likely impacts cell survival, antibody production, and inflammation. Taken together, we show that disruptions to the tightly regulated immune processes may tip the balance between exacerbated immune responses and allergy, or increased risk of infection and impaired resolution of inflammation.

A Transcriptomic Approach to Understand Patient Susceptibility to Pneumonia After Abdominal Surgery.

OBJECTIVE: To describe immune pathways and gene networks altered following major abdominal surgery and to identify transcriptomic patterns associated with postoperative pneumonia. BACKGROUND: Nosocomial infections are a major healthcare challenge, developing in over 20% of patients aged 45 or over undergoing major abdominal surgery, with postoperative pneumonia associated with an almost 5-fold increase in 30-day mortality. METHODS: From a prospective consecutive cohort (n=150) undergoing major abdominal surgery, whole-blood RNA was collected preoperatively and at 3 time-points postoperatively (2-6, 24, and 48 h). Twelve patients diagnosed with postoperative pneumonia and 27 matched patients remaining infection-free were identified for analysis with RNA-sequencing. RESULTS: Compared to preoperative sampling, 3639 genes were upregulated and 5043 downregulated at 2 to 6 hours. Pathway analysis demonstrated innate-immune activation with neutrophil degranulation and Toll-like-receptor signaling upregulation alongside adaptive-immune suppression. Cell-type deconvolution of preoperative RNA-sequencing revealed elevated S100A8/9-high neutrophils alongside reduced naïve CD4 T-cells in those later developing pneumonia. Preoperatively, a gene-signature characteristic of neutrophil degranulation was associated with postoperative pneumonia acquisition ( P =0.00092). A previously reported Sepsis Response Signature (SRSq) score, reflecting neutrophil dysfunction and a more dysregulated host response, at 48 hours postoperatively, differed between patients subsequently developing pneumonia and those remaining infection-free ( P =0.045). Analysis of the novel neutrophil gene-signature and SRSq scores in independent major abdominal surgery and polytrauma cohorts indicated good predictive performance in identifying patients suffering later infection. CONCLUSIONS: Major abdominal surgery acutely upregulates innate-immune pathways while simultaneously suppressing adaptive-immune pathways. This is more prominent in patients developing postoperative pneumonia. Preoperative transcriptomic signatures characteristic of neutrophil degranulation and postoperative SRSq scores may be useful predictors of subsequent pneumonia risk.

Single cell spatial analysis reveals inflammatory foci of immature neutrophil and CD8 T cells in COVID-19 lungs.

Single cell spatial interrogation of the immune-structural interactions in COVID -19 lungs is challenging, mainly because of the marked cellular infiltrate and architecturally distorted microstructure. To address this, we develop a suite of mathematical tools to search for statistically significant co-locations amongst immune and structural cells identified using 37-plex imaging mass cytometry. This unbiased method reveals a cellular map interleaved with an inflammatory network of immature neutrophils, cytotoxic CD8 T cells, megakaryocytes and monocytes co-located with regenerating alveolar progenitors and endothelium. Of note, a highly active cluster of immature neutrophils and CD8 T cells, is found spatially linked with alveolar progenitor cells, and temporally with the diffuse alveolar damage stage. These findings offer further insights into how immune cells interact in the lungs of severe COVID-19 disease. We provide our pipeline [Spatial Omics Oxford Pipeline (SpOOx)] and visual-analytical tool, Multi-Dimensional Viewer (MDV) software, as a resource for spatial analysis.

Extending the 'host response' paradigm from sepsis to cardiogenic shock: evidence, limitations and opportunities.

Recent clinical and research efforts in cardiogenic shock (CS) have largely focussed on the restoration of the low cardiac output state that is the conditio sine qua non of the clinical syndrome. This approach has failed to translate into improved outcomes, and mortality has remained static at 30-50%. There is an unmet need to better delineate the pathobiology of CS to understand the observed heterogeneity of presentation and treatment effect and to identify novel therapeutic targets. Despite data in other critical illness syndromes, specifically sepsis, the role of dysregulated inflammation and immunity is hitherto poorly described in CS. High-dimensional molecular profiling, particularly through leukocyte transcriptomics, may afford opportunity to better characterise subgroups of patients with shared mechanisms of immune dysregulation. In this state-of-the-art review, we outline the rationale for considering molecular subtypes of CS. We describe how high-dimensional molecular technologies can be used to identify these subtypes, and whether they share biological features with sepsis and other critical illness states. Finally, we propose how the identification of molecular subtypes of patients may enrich future clinical trial design and identification of novel therapies for CS.

Uncovering heterogeneity in sepsis: a comparative analysis of subphenotypes.

PURPOSE: The heterogeneity in sepsis is held responsible, in part, for the lack of precision treatment. Many attempts to identify subtypes of sepsis patients identify those with shared underlying biology or outcomes. To date, though, there has been limited effort to determine overlap across these previously identified subtypes. We aimed to determine the concordance of critically ill patients with sepsis classified by four previously described subtype strategies. METHODS: This secondary analysis of a multicenter prospective observational study included 522 critically ill patients with sepsis assigned to four previously established subtype strategies, primarily based on: (i) clinical data in the electronic health record (α, ß, γ, and δ), (ii) biomarker data (hyper- and hypoinflammatory), and (iii-iv) transcriptomic data (Mars1-Mars4 and SRS1-SRS2). Concordance was studied between different subtype labels, clinical characteristics, biological host response aberrations, as well as combinations of subtypes by sepsis ensembles. RESULTS: All four subtype labels could be adjudicated in this cohort, with the distribution of the clinical subtype varying most from the original cohort. The most common subtypes in each of the four strategies were γ (61%), which is higher compared to the original classification, hypoinflammatory (60%), Mars2 (35%), and SRS2 (54%). There was no clear relationship between any of the subtyping approaches (Cramer's V = 0.086-0.456). Mars2 and SRS1 were most alike in terms of host response biomarkers (p = 0.079-0.424), while other subtype strategies showed no clear relationship. Patients enriched for multiple subtypes revealed that characteristics and outcomes differ dependent on the combination of subtypes made. CONCLUSION: Among critically ill patients with sepsis, subtype strategies using clinical, biomarker, and transcriptomic data do not identify comparable patient populations and are likely to reflect disparate clinical characteristics and underlying biology.

Mapping the epigenomic landscape of human monocytes following innate immune activation reveals context-specific mechanisms driving endotoxin tolerance.

BACKGROUND: Monocytes are key mediators of innate immunity to infection, undergoing profound and dynamic changes in epigenetic state and immune function which are broadly protective but may be dysregulated in disease. Here, we aimed to advance understanding of epigenetic regulation following innate immune activation, acutely and in endotoxin tolerant states. METHODS: We exposed human primary monocytes from healthy donors (n = 6) to interferon-γ or differing combinations of endotoxin (lipopolysaccharide), including acute response (2 h) and two models of endotoxin tolerance: repeated stimulations (6 + 6 h) and prolonged exposure to endotoxin (24 h). Another subset of monocytes was left untreated (naïve). We identified context-specific regulatory elements based on epigenetic signatures for chromatin accessibility (ATAC-seq) and regulatory non-coding RNAs from total RNA sequencing. RESULTS: We present an atlas of differential gene expression for endotoxin and interferon response, identifying widespread context specific changes. Across assayed states, only 24-29% of genes showing differential exon usage are also differential at the gene level. Overall, 19.9% (6,884 of 34,616) of repeatedly observed ATAC peaks were differential in at least one condition, the majority upregulated on stimulation and located in distal regions (64.1% vs 45.9% of non-differential peaks) within which sequences were less conserved than non-differential peaks. We identified enhancer-derived RNA signatures specific to different monocyte states that correlated with chromatin accessibility changes. The endotoxin tolerance models showed distinct chromatin accessibility and transcriptomic signatures, with integrated analysis identifying genes and pathways involved in the inflammatory response, detoxification, metabolism and wound healing. We leveraged eQTL mapping for the same monocyte activation states to link potential enhancers with specific genes, identifying 1,946 unique differential ATAC peaks with 1,340 expression associated genes. We further use this to inform understanding of reported GWAS, for example involving FCHO1 and coronary artery disease. CONCLUSION: This study reports context-specific regulatory elements based on transcriptomic profiling and epigenetic signatures for enhancer-derived RNAs and chromatin accessibility in immune tolerant monocyte states, and demonstrates the informativeness of linking such elements and eQTL to inform future mechanistic studies aimed at defining therapeutic targets of immunosuppression and diseases.

SARS-CoV-2 infection induces a long-lived pro-inflammatory transcriptional profile.

BACKGROUND: The immune response to severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection in COVID-19 patients has been extensively investigated. However, much less is known about the long-term effects of infection in patients and how it could affect the immune system and its capacity to respond to future perturbations. METHODS: Using a targeted single-cell multiomics approach, we have recently identified a prolonged anti-inflammatory gene expression signature in T and NK cells in type 1 diabetes patients treated with low-dose IL-2. Here, we investigated the dynamics of this signature in three independent cohorts of COVID-19 patients: (i) the Oxford COVID-19 Multi-omics Blood Atlas (COMBAT) dataset, a cross-sectional cohort including 77 COVID-19 patients and ten healthy donors; (ii) the INCOV dataset, consisting of 525 samples taken from 209 COVID-19 patients during and after infection; and (iii) a longitudinal dataset consisting of 269 whole-blood samples taken from 139 COVID-19 patients followed for a period of up to 7 months after the onset of symptoms using a bulk transcriptomic approach. RESULTS: We discovered that SARS-CoV-2 infection leads to a prolonged alteration of the gene expression profile of circulating T, B and NK cells and monocytes. Some of the genes affected were the same as those present in the IL-2-induced anti-inflammatory gene expression signature but were regulated in the opposite direction, implying a pro-inflammatory status. The altered transcriptional profile was detected in COVID-19 patients for at least 2 months after the onset of the disease symptoms but was not observed in response to influenza infection or sepsis. Gene network analysis suggested a central role for the transcriptional factor NF-κB in the regulation of the observed transcriptional alterations. CONCLUSIONS: SARS-CoV-2 infection causes a prolonged increase in the pro-inflammatory transcriptional status that could predispose post-acute patients to the development of long-term health consequences, including autoimmune disease, reactivation of other viruses and disruption of the host immune system-microbiome ecosystem.

A genome-wide association study of blood cell morphology identifies cellular proteins implicated in disease aetiology.

Blood cells contain functionally important intracellular structures, such as granules, critical to immunity and thrombosis. Quantitative variation in these structures has not been subjected previously to large-scale genetic analysis. We perform genome-wide association studies of 63 flow-cytometry derived cellular phenotypes-including cell-type specific measures of granularity, nucleic acid content and reactivity-in 41,515 participants in the INTERVAL study. We identify 2172 distinct variant-trait associations, including associations near genes coding for proteins in organelles implicated in inflammatory and thrombotic diseases. By integrating with epigenetic data we show that many intracellular structures are likely to be determined in immature precursor cells. By integrating with proteomic data we identify the transcription factor FOG2 as an early regulator of platelet formation and α-granularity. Finally, we show that colocalisation of our associations with disease risk signals can suggest aetiological cell-types-variants in IL2RA and ITGA4 respectively mirror the known effects of daclizumab in multiple sclerosis and vedolizumab in inflammatory bowel disease.

Large clones of pre-existing T cells drive early immunity against SARS-COV-2 and LCMV infection.

T cell responses precede antibody and may provide early control of infection. We analyzed the clonal basis of this rapid response following SARS-COV-2 infection. We applied T cell receptor (TCR) sequencing to define the trajectories of individual T cell clones immediately. In SARS-COV-2 PCR+ individuals, a wave of TCRs strongly but transiently expand, frequently peaking the same week as the first positive PCR test. These expanding TCR CDR3s were enriched for sequences functionally annotated as SARS-COV-2 specific. Epitopes recognized by the expanding TCRs were highly conserved between SARS-COV-2 strains but not with circulating human coronaviruses. Many expanding CDR3s were present at high frequency in pre-pandemic repertoires. Early response TCRs specific for lymphocytic choriomeningitis virus epitopes were also found at high frequency in the preinfection naive repertoire. High-frequency naive precursors may allow the T cell response to respond rapidly during the crucial early phases of acute viral infection.

Comprehensive epigenomic profiling reveals the extent of disease-specific chromatin states and informs target discovery in ankylosing spondylitis.

Ankylosing spondylitis (AS) is a common, highly heritable inflammatory arthritis characterized by enthesitis of the spine and sacroiliac joints. Genome-wide association studies (GWASs) have revealed more than 100 genetic associations whose functional effects remain largely unresolved. Here, we present a comprehensive transcriptomic and epigenomic map of disease-relevant blood immune cell subsets from AS patients and healthy controls. We find that, while CD14+ monocytes and CD4+ and CD8+ T cells show disease-specific differences at the RNA level, epigenomic differences are only apparent upon multi-omics integration. The latter reveals enrichment at disease-associated loci in monocytes. We link putative functional SNPs to genes using high-resolution Capture-C at 10 loci, including PTGER4 and ETS1, and show how disease-specific functional genomic data can be integrated with GWASs to enhance therapeutic target discovery. This study combines epigenetic and transcriptional analysis with GWASs to identify disease-relevant cell types and gene regulation of likely pathogenic relevance and prioritize drug targets.

Evaluation of T cell responses to naturally processed variant SARS-CoV-2 spike antigens in individuals following infection or vaccination.

Most existing studies characterizing severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2)-specific T cell responses are peptide based. This does not allow evaluation of whether tested peptides are processed and presented canonically. In this study, we use recombinant vaccinia virus (rVACV)-mediated expression of SARS-CoV-2 spike protein and SARS-CoV-2 infection of angiotensin-converting enzyme (ACE)-2-transduced B cell lines to evaluate overall T cell responses in a small cohort of recovered COVID-19 patients and uninfected donors vaccinated with ChAdOx1 nCoV-19. We show that rVACV expression of SARS-CoV-2 antigen can be used as an alternative to SARS-CoV-2 infection to evaluate T cell responses to naturally processed spike antigens. In addition, the rVACV system can be used to evaluate the cross-reactivity of memory T cells to variants of concern (VOCs) and to identify epitope escape mutants. Finally, our data show that both natural infection and vaccination could induce multi-functional T cell responses with overall T cell responses remaining despite the identification of escape mutations.

Neutrophils and emergency granulopoiesis drive immune suppression and an extreme response endotype during sepsis.

Sepsis arises from diverse and incompletely understood dysregulated host response processes following infection that leads to life-threatening organ dysfunction. Here we showed that neutrophils and emergency granulopoiesis drove a maladaptive response during sepsis. We generated a whole-blood single-cell multiomic atlas (272,993 cells, n = 39 individuals) of the sepsis immune response that identified populations of immunosuppressive mature and immature neutrophils. In co-culture, CD66b+ sepsis neutrophils inhibited proliferation and activation of CD4+ T cells. Single-cell multiomic mapping of circulating hematopoietic stem and progenitor cells (HSPCs) (29,366 cells, n = 27) indicated altered granulopoiesis in patients with sepsis. These features were enriched in a patient subset with poor outcome and a specific sepsis response signature that displayed higher frequencies of IL1R2+ immature neutrophils, epigenetic and transcriptomic signatures of emergency granulopoiesis in HSPCs and STAT3-mediated gene regulation across different infectious etiologies and syndromes. Our findings offer potential therapeutic targets and opportunities for stratified medicine in severe infection.

Chromosome conformation capture approaches to investigate 3D genome architecture in Ankylosing Spondylitis.

Ankylosing Spondylitis (AS) is a chronic inflammatory arthritis of the spine exhibiting a strong genetic background. The mechanistic and functional understanding of the AS-associated genomic loci, identified with Genome Wide Association Studies (GWAS), remains challenging. Chromosome conformation capture (3C) and derivatives are recent techniques which are of great help in elucidating the spatial genome organization and of enormous support in uncover a mechanistic explanation for disease-associated genetic variants. The perturbation of three-dimensional (3D) genome hierarchy may lead to a plethora of human diseases, including rheumatological disorders. Here we illustrate the latest approaches and related findings on the field of genome organization, highlighting how the instability of 3D genome conformation may be among the causes of rheumatological disease phenotypes. We suggest a new perspective on the inclusive potential of a 3C approach to inform GWAS results in rheumatic diseases. 3D genome organization may ultimately lead to a more precise and comprehensive functional interpretation of AS association, which is the starting point for emerging and more specific therapies.

COMBATdb: a database for the COVID-19 Multi-Omics Blood ATlas.

Advances in our understanding of the nature of the immune response to SARS-CoV-2 infection, and how this varies within and between individuals, is important in efforts to develop targeted therapies and precision medicine approaches. Here we present a database for the COvid-19 Multi-omics Blood ATlas (COMBAT) project, COMBATdb (https://db.combat.ox.ac.uk). This enables exploration of multi-modal datasets arising from profiling of patients with different severities of illness admitted to hospital in the first phase of the pandemic in the UK prior to vaccination, compared with community cases, healthy controls, and patients with all-cause sepsis and influenza. These data include whole blood transcriptomics, plasma proteomics, epigenomics, single-cell multi-omics, immune repertoire sequencing, flow and mass cytometry, and cohort metadata. COMBATdb provides access to the processed data in a well-defined framework of samples, cell types and genes/proteins that allows exploration across the assayed modalities, with functionality including browse, search, download, calculation and visualisation via shiny apps. This advances the ability of users to leverage COMBAT datasets to understand the pathogenesis of COVID-19, and the nature of specific and shared features with other infectious diseases.

Human leukocyte antigen alleles associate with COVID-19 vaccine immunogenicity and risk of breakthrough infection.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) vaccine immunogenicity varies between individuals, and immune responses correlate with vaccine efficacy. Using data from 1,076 participants enrolled in ChAdOx1 nCov-19 vaccine efficacy trials in the United Kingdom, we found that inter-individual variation in normalized antibody responses against SARS-CoV-2 spike and its receptor-binding domain (RBD) at 28 days after first vaccination shows genome-wide significant association with major histocompatibility complex (MHC) class II alleles. The most statistically significant association with higher levels of anti-RBD antibody was HLA-DQB1*06 (P = 3.2 × 10-9), which we replicated in 1,677 additional vaccinees. Individuals carrying HLA-DQB1*06 alleles were less likely to experience PCR-confirmed breakthrough infection during the ancestral SARS-CoV-2 virus and subsequent Alpha variant waves compared to non-carriers (hazard ratio = 0.63, 0.42-0.93, P = 0.02). We identified a distinct spike-derived peptide that is predicted to bind differentially to HLA-DQB1*06 compared to other similar alleles, and we found evidence of increased spike-specific memory B cell responses in HLA-DQB1*06 carriers at 84 days after first vaccination. Our results demonstrate association of HLA type with Coronavirus Disease 2019 (COVID-19) vaccine antibody response and risk of breakthrough infection, with implications for future vaccine design and implementation.