T cell reactivity to the SARS-CoV-2 Omicron variant is preserved in most but not all individuals.

The SARS-CoV-2 Omicron variant (B.1.1.529) contains mutations that mediate escape from antibody responses, although the extent to which these substitutions in spike and non-spike proteins affect T cell recognition is unknown. In this study, we show that T cell responses in individuals with prior infection, vaccination, both prior infection and vaccination, and boosted vaccination are largely preserved to Omicron spike and non-spike proteins. However, we also identify a subset of individuals (∼21%) with a >50% reduction in T cell reactivity to the Omicron spike. Evaluation of functional CD4+ and CD8+ memory T cell responses confirmed these findings and revealed that reduced recognition to Omicron spike is primarily observed within the CD8+ T cell compartment potentially due to escape from HLA binding. Booster vaccination enhanced T cell responses to Omicron spike. In contrast to neutralizing immunity, these findings suggest preservation of T cell responses to the Omicron variant, although with reduced reactivity in some individuals.

mRNA-based COVID-19 vaccine boosters induce neutralizing immunity against SARS-CoV-2 Omicron variant.

Recent surveillance has revealed the emergence of the SARS-CoV-2 Omicron variant (BA.1/B.1.1.529) harboring up to 36 mutations in spike protein, the target of neutralizing antibodies. Given its potential to escape vaccine-induced humoral immunity, we measured the neutralization potency of sera from 88 mRNA-1273, 111 BNT162b, and 40 Ad26.COV2.S vaccine recipients against wild-type, Delta, and Omicron SARS-CoV-2 pseudoviruses. We included individuals that received their primary series recently (<3 months), distantly (6-12 months), or an additional "booster" dose, while accounting for prior SARS-CoV-2 infection. Remarkably, neutralization of Omicron was undetectable in most vaccinees. However, individuals boosted with mRNA vaccines exhibited potent neutralization of Omicron, only 4-6-fold lower than wild type, suggesting enhanced cross-reactivity of neutralizing antibody responses. In addition, we find that Omicron pseudovirus infects more efficiently than other variants tested. Overall, this study highlights the importance of additional mRNA doses to broaden neutralizing antibody responses against highly divergent SARS-CoV-2 variants.

Structure-guided T cell vaccine design for SARS-CoV-2 variants and sarbecoviruses.

The emergence of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) variants that escape convalescent and vaccine-induced antibody responses has renewed focus on the development of broadly protective T-cell-based vaccines. Here, we apply structure-based network analysis and assessments of HLA class I peptide stability to define mutationally constrained CD8+ T cell epitopes across the SARS-CoV-2 proteome. Highly networked residues are conserved temporally among circulating variants and sarbecoviruses and disproportionately impair spike pseudotyped lentivirus infectivity when mutated. Evaluation of HLA class I stabilizing activity for 18 globally prevalent alleles identifies CD8+ T cell epitopes within highly networked regions with limited mutational frequencies in circulating SARS-CoV-2 variants and deep-sequenced primary isolates. Moreover, these epitopes elicit demonstrable CD8+ T cell reactivity in convalescent individuals but reduced recognition in recipients of mRNA-based vaccines. These data thereby elucidate key mutationally constrained regions and immunogenic epitopes in the SARS-CoV-2 proteome for a global T-cell-based vaccine against emerging variants and SARS-like coronaviruses.

Multiple SARS-CoV-2 variants escape neutralization by vaccine-induced humoral immunity.

Vaccination elicits immune responses capable of potently neutralizing SARS-CoV-2. However, ongoing surveillance has revealed the emergence of variants harboring mutations in spike, the main target of neutralizing antibodies. To understand the impact of these variants, we evaluated the neutralization potency of 99 individuals that received one or two doses of either BNT162b2 or mRNA-1273 vaccines against pseudoviruses representing 10 globally circulating strains of SARS-CoV-2. Five of the 10 pseudoviruses, harboring receptor-binding domain mutations, including K417N/T, E484K, and N501Y, were highly resistant to neutralization. Cross-neutralization of B.1.351 variants was comparable to SARS-CoV and bat-derived WIV1-CoV, suggesting that a relatively small number of mutations can mediate potent escape from vaccine responses. While the clinical impact of neutralization resistance remains uncertain, these results highlight the potential for variants to escape from neutralizing humoral immunity and emphasize the need to develop broadly protective interventions against the evolving pandemic.

HLA class I signal peptide polymorphism determines the level of CD94/NKG2-HLA-E-mediated regulation of effector cell responses.

Human leukocyte antigen (HLA)-E binds epitopes derived from HLA-A, HLA-B, HLA-C and HLA-G signal peptides (SPs) and serves as a ligand for CD94/NKG2A and CD94/NKG2C receptors expressed on natural killer and T cell subsets. We show that among 16 common classical HLA class I SP variants, only 6 can be efficiently processed to generate epitopes that enable CD94/NKG2 engagement, which we term 'functional SPs'. The single functional HLA-B SP, known as HLA-B/-21M, induced high HLA-E expression, but conferred the lowest receptor recognition. Consequently, HLA-B/-21M SP competes with other SPs for providing epitope to HLA-E and reduces overall recognition of target cells by CD94/NKG2A, calling for reassessment of previous disease models involving HLA-B/-21M. Genetic population data indicate a positive correlation between frequencies of functional SPs in humans and corresponding cytomegalovirus mimics, suggesting a means for viral escape from host responses. The systematic, quantitative approach described herein will facilitate development of prediction algorithms for accurately measuring the impact of CD94/NKG2-HLA-E interactions in disease resistance/susceptibility.

Author Correction: CCR5AS lncRNA variation differentially regulates CCR5, influencing HIV disease outcome.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

CCR5AS lncRNA variation differentially regulates CCR5, influencing HIV disease outcome.

Multiple genome-wide studies have identified associations between outcome of human immunodeficiency virus (HIV) infection and polymorphisms in and around the gene encoding the HIV co-receptor CCR5, but the functional basis for the strongest of these associations, rs1015164A/G, is unknown. We found that rs1015164 marks variation in an activating transcription factor 1 binding site that controls expression of the antisense long noncoding RNA (lncRNA) CCR5AS. Knockdown or enhancement of CCR5AS expression resulted in a corresponding change in CCR5 expression on CD4+ T cells. CCR5AS interfered with interactions between the RNA-binding protein Raly and the CCR5 3' untranslated region, protecting CCR5 messenger RNA from Raly-mediated degradation. Reduction in CCR5 expression through inhibition of CCR5AS diminished infection of CD4+ T cells with CCR5-tropic HIV in vitro. These data represent a rare determination of the functional importance of a genome-wide disease association where expression of a lncRNA affects HIV infection and disease progression.

Mechanistic model for booster doses effectiveness in healthy, cancer, and immunosuppressed patients infected with SARS-CoV-2.

SARS-CoV-2 vaccines are effective at limiting disease severity, but effectiveness is lower among patients with cancer or immunosuppression. Effectiveness wanes with time and varies by vaccine type. Moreover, previously prescribed vaccines were based on the ancestral SARS-CoV-2 spike-protein that emerging variants may evade. Here, we describe a mechanistic mathematical model for vaccination-induced immunity. We validate it with available clinical data and use it to simulate the effectiveness of vaccines against viral variants with lower antigenicity, increased virulence, or enhanced cell binding for various vaccine platforms. The analysis includes the omicron variant as well as hypothetical future variants with even greater immune evasion of vaccine-induced antibodies and addresses the potential benefits of the new bivalent vaccines. We further account for concurrent cancer or underlying immunosuppression. The model confirms enhanced immunogenicity following booster vaccination in immunosuppressed patients but predicts ongoing booster requirements for these individuals to maintain protection. We further studied the impact of variants on immunosuppressed individuals as a function of the interval between multiple booster doses. Our model suggests possible strategies for future vaccinations and suggests tailored strategies for high-risk groups.

Genetic variation that determines TAPBP expression levels associates with the course of malaria in an HLA allotype-dependent manner.

HLA class I (HLA-I) allotypes vary widely in their dependence on tapasin (TAPBP), an integral component of the peptide-loading complex, to present peptides on the cell surface. We identified two single-nucleotide polymorphisms that regulate TAPBP messenger RNA (mRNA) expression in Africans, rs111686073 (G/C) and rs59097151 (A/G), located in an AP-2α transcription factor binding site and a microRNA (miR)-4486 binding site, respectively. rs111686073G and rs59097151A induced significantly higher TAPBP mRNA expression relative to the alternative alleles due to higher affinity for AP-2α and abrogation of miR-4486 binding, respectively. These variants associated with lower Plasmodium falciparum parasite prevalence and lower incidence of clinical malaria specifically among individuals carrying tapasin-dependent HLA-I allotypes, presumably by augmenting peptide loading, whereas tapasin-independent allotypes associated with relative protection, regardless of imputed TAPBP mRNA expression levels. Thus, an attenuated course of malaria may occur through enhanced breadth and/or magnitude of antigen presentation, an important consideration when evaluating vaccine efficacy.

Genome-wide association study of leprosy in Malawi and Mali.

Leprosy is a chronic infection of the skin and peripheral nerves caused by Mycobacterium leprae. Despite recent improvements in disease control, leprosy remains an important cause of infectious disability globally. Large-scale genetic association studies in Chinese, Vietnamese and Indian populations have identified over 30 susceptibility loci for leprosy. There is a significant burden of leprosy in Africa, however it is uncertain whether the findings of published genetic association studies are generalizable to African populations. To address this, we conducted a genome-wide association study (GWAS) of leprosy in Malawian (327 cases, 436 controls) and Malian (247 cases, 368 controls) individuals. In that analysis, we replicated four risk loci previously reported in China, Vietnam and India; MHC Class I and II, LACC1 and SLC29A3. We further identified a novel leprosy susceptibility locus at 10q24 (rs2015583; combined p = 8.81 × 10-9; OR = 0.51 [95% CI 0.40 - 0.64]). Using publicly-available data we characterise regulatory activity at this locus, identifying ACTR1A as a candidate mediator of leprosy risk. This locus shows evidence of recent positive selection and demonstrates pleiotropy with established risk loci for inflammatory bowel disease and childhood-onset asthma. A shared genetic architecture for leprosy and inflammatory bowel disease has been previously described. We expand on this, strengthening the hypothesis that selection pressure driven by leprosy has shaped the evolution of autoimmune and atopic disease in modern populations. More broadly, our data highlights the importance of defining the genetic architecture of disease across genetically diverse populations, and that disease insights derived from GWAS in one population may not translate to all affected populations.

HLA tapasin independence: broader peptide repertoire and HIV control.

Human leukocyte antigen (HLA) class I allotypes vary in their ability to present peptides in the absence of tapasin, an essential component of the peptide loading complex. We quantified tapasin dependence of all allotypes that are common in European and African Americans (n = 97), which revealed a broad continuum of values. Ex vivo examination of cytotoxic T cell responses to the entire HIV-1 proteome from infected subjects indicates that tapasin-dependent allotypes present a more limited set of distinct peptides than do tapasin-independent allotypes, data supported by computational predictions. This suggests that variation in tapasin dependence may impact the strength of the immune responses by altering peptide repertoire size. In support of this model, we observed that individuals carrying HLA class I genotypes characterized by greater tapasin independence progress more slowly to AIDS and maintain lower viral loads, presumably due to increased breadth of peptide presentation. Thus, tapasin dependence level, like HLA zygosity, may serve as a means to restrict or expand breadth of the HLA-I peptide repertoire across humans, ultimately influencing immune responses to pathogens and vaccines.

Comparative Immunogenicity and Effectiveness of mRNA-1273, BNT162b2, and Ad26.COV2.S COVID-19 Vaccines.

BACKGROUND: Understanding immunogenicity and effectiveness of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) vaccines is critical to guide rational use. METHODS: We compared the immunogenicity of mRNA-1273, BNT-162b2, and Ad26.COV2.S in healthy ambulatory adults. We performed an inverse-variance meta-analysis of population-level effectiveness from public health reports in > 40 million individuals. RESULTS: A single dose of either mRNA vaccine yielded comparable antibody and neutralization titers to convalescent individuals. Ad26.COV2.S yielded lower antibody concentrations and frequently undetectable neutralization titers. Bulk and cytotoxic T-cell responses were higher in mRNA1273 and BNT162b2 than Ad26.COV2.S recipients. Regardless of vaccine, <50% of vaccinees demonstrated CD8+ T-cell responses. Antibody concentrations and neutralization titers increased comparably after the first dose of either vaccine, and further in recipients of a second dose. Prior infection was associated with high antibody concentrations and neutralization even after a single dose and regardless of vaccine. Neutralization of Beta, Gamma, and Delta strains were poorer regardless of vaccine. In meta-analysis, relative to mRNA1273 the effectiveness of BNT162b2 was lower against infection and hospitalization, and Ad26COV2.S was lower against infection, hospitalization, and death. CONCLUSIONS: Variation in the immunogenicity correlates with variable effectiveness of the 3 vaccines deployed in the United States.

HLA-A*03 and response to immune checkpoint blockade in cancer: an epidemiological biomarker study.

BACKGROUND: Predictive biomarkers could allow more precise use of immune checkpoint inhibitors (ICIs) in treating advanced cancers. Given the central role of HLA molecules in immunity, variation at the HLA loci could differentially affect the response to ICIs. The aim of this epidemiological study was to determine the effect of HLA-A*03 as a biomarker for predicting response to immunotherapy. METHODS: In this epidemiological study, we investigated the clinical outcomes (overall survival, progression free survival, and objective response rate) after treatment for advanced cancer in eight cohorts of patients: three observational cohorts of patients with various types of advanced tumours (the Memorial Sloan Kettering Integrated Mutation Profiling of Actionable Cancer Targets [MSK-IMPACT] cohort, the Dana-Farber Cancer Institute [DFCI] Profile cohort, and The Cancer Genome Atlas) and five clinical trials of patients with advanced bladder cancer (JAVELIN Solid Tumour) or renal cell carcinoma (CheckMate-009, CheckMate-010, CheckMate-025, and JAVELIN Renal 101). In total, these cohorts included 3335 patients treated with various ICI agents (anti-PD-1, anti-PD-L1, and anti-CTLA-4 inhibitors) and 10 917 patients treated with non-ICI cancer-directed therapeutic approaches. We initially modelled the association of HLA amino-acid variation with overall survival in the MSK-IMPACT discovery cohort, followed by a detailed analysis of the association between HLA-A*03 and clinical outcomes in MSK-IMPACT, with replication in the additional cohorts (two further observational cohorts and five clinical trials). FINDINGS: HLA-A*03 was associated in an additive manner with reduced overall survival after ICI treatment in the MSK-IMPACT cohort (HR 1·48 per HLA-A*03 allele [95% CI 1·20-1·82], p=0·00022), the validation DFCI Profile cohort (HR 1·22 per HLA-A*03 allele, 1·05-1·42; p=0·0097), and in the JAVELIN Solid Tumour clinical trial for bladder cancer (HR 1·36 per HLA-A*03 allele, 1·01-1·85; p=0·047). The HLA-A*03 effect was observed across ICI agents and tumour types, but not in patients treated with alternative therapies. Patients with HLA-A*03 had shorter progression-free survival in the pooled patient population from the three CheckMate clinical trials of nivolumab for renal cell carcinoma (HR 1·31, 1·01-1·71; p=0·044), but not in those receiving control (everolimus) therapies. Objective responses were observed in none of eight HLA-A*03 homozygotes in the ICI group (compared with 59 [26·6%] of 222 HLA-A*03 non-carriers and 13 (17·1%) of 76 HLA-A*03 heterozygotes). HLA-A*03 was associated with shorter progression-free survival in patients receiving ICI in the JAVELIN Renal 101 randomised clinical trial for renal cell carcinoma (avelumab plus axitinib; HR 1·59 per HLA-A*03 allele, 1·16-2·16; p=0·0036), but not in those receiving control (sunitinib) therapy. Objective responses were recorded in one (12·5%) of eight HLA-A*03 homozygotes in the ICI group (compared with 162 [63·8%] of 254 HLA-A*03 non-carriers and 40 [55·6%] of 72 HLA-A*03 heterozygotes). HLA-A*03 was associated with impaired outcome in meta-analysis of all 3335 patients treated with ICI at genome-wide significance (p=2·01 × 10-8) with no evidence of heterogeneity in effect (I2 0%, 95% CI 0-0·76) INTERPRETATION: HLA-A*03 is a predictive biomarker of poor response to ICI. Further evaluation of HLA-A*03 is warranted in randomised trials. HLA-A*03 carriage could be considered in decisions to initiate ICI in patients with cancer. FUNDING: National Institutes of Health, Merck KGaA, and Pfizer.

High Seroprevalence of Anti-SARS-CoV-2 Antibodies in Chelsea, Massachusetts.

SARS-CoV-2 antibody testing allows quantitative determination of disease prevalence, which is especially important in high-risk communities. We performed anonymized convenience sampling of 200 currently asymptomatic residents of Chelsea, the epicenter of COVID-19 illness in Massachusetts, by BioMedomics SARS-CoV-2 combined IgM-IgG point-of-care lateral flow immunoassay. The seroprevalence was 31.5% (17.5% IgM+IgG+, 9.0% IgM+IgG-, and 5.0% IgM-IgG+). Of the 200 participants, 50.5% reported no symptoms in the preceding 4 weeks, of which 24.8% (25/101) were seropositive, and 60% of these were IgM+IgG-. These data are the highest seroprevalence rates observed to date and highlight the significant burden of asymptomatic infection.

The eMERGE genotype set of 83,717 subjects imputed to ~40 million variants genome wide and association with the herpes zoster medical record phenotype.

The Electronic Medical Records and Genomics (eMERGE) network is a network of medical centers with electronic medical records linked to existing biorepository samples for genomic discovery and genomic medicine research. The network sought to unify the genetic results from 78 Illumina and Affymetrix genotype array batches from 12 contributing medical centers for joint association analysis of 83,717 human participants. In this report, we describe the imputation of eMERGE results and methods to create the unified imputed merged set of genome-wide variant genotype data. We imputed the data using the Michigan Imputation Server, which provides a missing single-nucleotide variant genotype imputation service using the minimac3 imputation algorithm with the Haplotype Reference Consortium genotype reference set. We describe the quality control and filtering steps used in the generation of this data set and suggest generalizable quality thresholds for imputation and phenotype association studies. To test the merged imputed genotype set, we replicated a previously reported chromosome 6 HLA-B herpes zoster (shingles) association and discovered a novel zoster-associated loci in an epigenetic binding site near the terminus of chromosome 3 (3p29).

HLA-C Level Is Regulated by a Polymorphic Oct1 Binding Site in the HLA-C Promoter Region.

Differential HLA-C levels influence several human diseases, but the mechanisms responsible are incompletely characterized. Using a validated prediction algorithm, we imputed HLA-C cell surface levels in 228 individuals from the 1000 Genomes dataset. We tested 68,726 SNPs within the MHC for association with HLA-C level. The HLA-C promoter region variant, rs2395471, 800 bp upstream of the transcription start site, gave the most significant association with HLA-C levels (p = 4.2 × 10-66). This imputed expression quantitative trait locus, termed impeQTL, was also shown to associate with HLA-C expression in a genome-wide association study of 273 donors in which HLA-C mRNA expression levels were determined by quantitative PCR (qPCR) (p = 1.8 × 10-20) and in two cohorts where HLA-C cell surface levels were determined directly by flow cytometry (n = 369 combined, p < 10-15). rs2395471 is located in an Oct1 transcription factor consensus binding site motif where the A allele is predicted to have higher affinity for Oct1 than the G allele. Mobility shift electrophoresis demonstrated that Oct1 binds to both alleles in vitro, but decreased HLA-C promoter activity was observed in a luciferase reporter assay for rs2395471_G relative to rs2395471_A on a fixed promoter background. The rs2395471 variant accounts for up to 36% of the explained variation of HLA-C level. These data strengthen our understanding of HLA-C transcriptional regulation and provide a basis for understanding the potential consequences of manipulating HLA-C levels therapeutically.

Polymorphism in a lincRNA Associates with a Doubled Risk of Pneumococcal Bacteremia in Kenyan Children.

Bacteremia (bacterial bloodstream infection) is a major cause of illness and death in sub-Saharan Africa but little is known about the role of human genetics in susceptibility. We conducted a genome-wide association study of bacteremia susceptibility in more than 5,000 Kenyan children as part of the Wellcome Trust Case Control Consortium 2 (WTCCC2). Both the blood-culture-proven bacteremia case subjects and healthy infants as controls were recruited from Kilifi, on the east coast of Kenya. Streptococcus pneumoniae is the most common cause of bacteremia in Kilifi and was thus the focus of this study. We identified an association between polymorphisms in a long intergenic non-coding RNA (lincRNA) gene (AC011288.2) and pneumococcal bacteremia and replicated the results in the same population (p combined = 1.69 × 10(-9); OR = 2.47, 95% CI = 1.84-3.31). The susceptibility allele is African specific, derived rather than ancestral, and occurs at low frequency (2.7% in control subjects and 6.4% in case subjects). Our further studies showed AC011288.2 expression only in neutrophils, a cell type that is known to play a major role in pneumococcal clearance. Identification of this novel association will further focus research on the role of lincRNAs in human infectious disease.

HLA Class I-Mediated HIV-1 Control in Vietnamese Infected with HIV-1 Subtype A/E.

HIV-1-specific cytotoxic T cells (CTLs) play an important role in the control of HIV-1 subtype B or C infection. However, the role of CTLs in HIV-1 subtype A/E infection still remains unclear. Here we investigated the association of HLA class I alleles with clinical outcomes in treatment-naive Vietnamese infected with subtype A/E virus. We found that HLA-C*12:02 was significantly associated with lower plasma viral loads (pVL) and higher CD4 counts and that the HLA-A*29:01-B*07:05-C*15:05 haplotype was significantly associated with higher pVL and lower CD4 counts than those for individuals without these respective genotypes. Nine Pol and three Nef mutations were associated with at least one HLA allele in the HLA-A*29:01-B*07:05-C*15:05 haplotype, with a strong negative correlation between the number of HLA-associated Pol mutations and CD4 count as well as a positive correlation with pVL for individuals with these HLA alleles. The results suggest that the accumulation of mutations selected by CTLs restricted by these HLA alleles affects HIV control.IMPORTANCE Most previous studies on HLA association with disease progression after HIV-1 infection have been performed on cohorts infected with HIV-1 subtypes B and C, whereas few such population-based studies have been reported for cohorts infected with the Asian subtype A/E virus. In this study, we analyzed the association of HLA class I alleles with clinical outcomes for 536 HIV-1 subtype A/E-infected Vietnamese individuals. We found that HLA-C*12:02 is protective, while the HLA haplotype HLA-A*29:01-B*07:05-C*15:05 is deleterious. The individuals with HIV-1 mutations associated with at least one of the HLA alleles in the deleterious HLA haplotype had higher plasma viral loads and lower CD4 counts than those of individuals without the mutations, suggesting that viral adaptation and escape from HLA-mediated immune control occurred. The present study identifies a protective allele and a deleterious haplotype for HIV-1 subtype A/E infection which are different from those identified for cohorts infected with HIV-1 subtypes B and C.