Evidence of escape of SARS-CoV-2 variant B.1.351 from natural and vaccine-induced sera.

The race to produce vaccines against severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) began when the first sequence was published, and this forms the basis for vaccines currently deployed globally. Independent lineages of SARS-CoV-2 have recently been reported: UK, B.1.1.7; South Africa, B.1.351; and Brazil, P.1. These variants have multiple changes in the immunodominant spike protein that facilitates viral cell entry via the angiotensin-converting enzyme-2 (ACE2) receptor. Mutations in the receptor recognition site on the spike are of great concern for their potential for immune escape. Here, we describe a structure-function analysis of B.1.351 using a large cohort of convalescent and vaccinee serum samples. The receptor-binding domain mutations provide tighter ACE2 binding and widespread escape from monoclonal antibody neutralization largely driven by E484K, although K417N and N501Y act together against some important antibody classes. In a number of cases, it would appear that convalescent and some vaccine serum offers limited protection against this variant.

The antigenic anatomy of SARS-CoV-2 receptor binding domain.

Antibodies are crucial to immune protection against SARS-CoV-2, with some in emergency use as therapeutics. Here, we identify 377 human monoclonal antibodies (mAbs) recognizing the virus spike and focus mainly on 80 that bind the receptor binding domain (RBD). We devise a competition data-driven method to map RBD binding sites. We find that although antibody binding sites are widely dispersed, neutralizing antibody binding is focused, with nearly all highly inhibitory mAbs (IC50 < 0.1 µg/mL) blocking receptor interaction, except for one that binds a unique epitope in the N-terminal domain. Many of these neutralizing mAbs use public V-genes and are close to germline. We dissect the structural basis of recognition for this large panel of antibodies through X-ray crystallography and cryoelectron microscopy of 19 Fab-antigen structures. We find novel binding modes for some potently inhibitory antibodies and demonstrate that strongly neutralizing mAbs protect, prophylactically or therapeutically, in animal models.

Reduced neutralization of SARS-CoV-2 B.1.1.7 variant by convalescent and vaccine sera.

SARS-CoV-2 has caused over 2 million deaths in little over a year. Vaccines are being deployed at scale, aiming to generate responses against the virus spike. The scale of the pandemic and error-prone virus replication is leading to the appearance of mutant viruses and potentially escape from antibody responses. Variant B.1.1.7, now dominant in the UK, with increased transmission, harbors 9 amino acid changes in the spike, including N501Y in the ACE2 interacting surface. We examine the ability of B.1.1.7 to evade antibody responses elicited by natural SARS-CoV-2 infection or vaccination. We map the impact of N501Y by structure/function analysis of a large panel of well-characterized monoclonal antibodies. B.1.1.7 is harder to neutralize than parental virus, compromising neutralization by some members of a major class of public antibodies through light-chain contacts with residue 501. However, widespread escape from monoclonal antibodies or antibody responses generated by natural infection or vaccination was not observed.

Antibody evasion by the P.1 strain of SARS-CoV-2.

Terminating the SARS-CoV-2 pandemic relies upon pan-global vaccination. Current vaccines elicit neutralizing antibody responses to the virus spike derived from early isolates. However, new strains have emerged with multiple mutations, including P.1 from Brazil, B.1.351 from South Africa, and B.1.1.7 from the UK (12, 10, and 9 changes in the spike, respectively). All have mutations in the ACE2 binding site, with P.1 and B.1.351 having a virtually identical triplet (E484K, K417N/T, and N501Y), which we show confer similar increased affinity for ACE2. We show that, surprisingly, P.1 is significantly less resistant to naturally acquired or vaccine-induced antibody responses than B.1.351, suggesting that changes outside the receptor-binding domain (RBD) impact neutralization. Monoclonal antibody (mAb) 222 neutralizes all three variants despite interacting with two of the ACE2-binding site mutations. We explain this through structural analysis and use the 222 light chain to largely restore neutralization potency to a major class of public antibodies.

Reduced neutralization of SARS-CoV-2 B.1.617 by vaccine and convalescent serum.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has undergone progressive change, with variants conferring advantage rapidly becoming dominant lineages, e.g., B.1.617. With apparent increased transmissibility, variant B.1.617.2 has contributed to the current wave of infection ravaging the Indian subcontinent and has been designated a variant of concern in the United Kingdom. Here we study the ability of monoclonal antibodies and convalescent and vaccine sera to neutralize B.1.617.1 and B.1.617.2, complement this with structural analyses of Fab/receptor binding domain (RBD) complexes, and map the antigenic space of current variants. Neutralization of both viruses is reduced compared with ancestral Wuhan-related strains, but there is no evidence of widespread antibody escape as seen with B.1.351. However, B.1.351 and P.1 sera showed markedly more reduction in neutralization of B.1.617.2, suggesting that individuals infected previously by these variants may be more susceptible to reinfection by B.1.617.2. This observation provides important new insights for immunization policy with future variant vaccines in non-immune populations.

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.

GWAS and meta-analysis identifies 49 genetic variants underlying critical COVID-19.

Critical illness in COVID-19 is an extreme and clinically homogeneous disease phenotype that we have previously shown1 to be highly efficient for discovery of genetic associations2. Despite the advanced stage of illness at presentation, we have shown that host genetics in patients who are critically ill with COVID-19 can identify immunomodulatory therapies with strong beneficial effects in this group3. Here we analyse 24,202 cases of COVID-19 with critical illness comprising a combination of microarray genotype and whole-genome sequencing data from cases of critical illness in the international GenOMICC (11,440 cases) study, combined with other studies recruiting hospitalized patients with a strong focus on severe and critical disease: ISARIC4C (676 cases) and the SCOURGE consortium (5,934 cases). To put these results in the context of existing work, we conduct a meta-analysis of the new GenOMICC genome-wide association study (GWAS) results with previously published data. We find 49 genome-wide significant associations, of which 16 have not been reported previously. To investigate the therapeutic implications of these findings, we infer the structural consequences of protein-coding variants, and combine our GWAS results with gene expression data using a monocyte transcriptome-wide association study (TWAS) model, as well as gene and protein expression using Mendelian randomization. We identify potentially druggable targets in multiple systems, including inflammatory signalling (JAK1), monocyte-macrophage activation and endothelial permeability (PDE4A), immunometabolism (SLC2A5 and AK5), and host factors required for viral entry and replication (TMPRSS2 and RAB2A).

Whole-genome sequencing reveals host factors underlying critical COVID-19.

Critical COVID-19 is caused by immune-mediated inflammatory lung injury. Host genetic variation influences the development of illness requiring critical care1 or hospitalization2-4 after infection with SARS-CoV-2. The GenOMICC (Genetics of Mortality in Critical Care) study enables the comparison of genomes from individuals who are critically ill with those of population controls to find underlying disease mechanisms. Here we use whole-genome sequencing in 7,491 critically ill individuals compared with 48,400 controls to discover and replicate 23 independent variants that significantly predispose to critical COVID-19. We identify 16 new independent associations, including variants within genes that are involved in interferon signalling (IL10RB and PLSCR1), leucocyte differentiation (BCL11A) and blood-type antigen secretor status (FUT2). Using transcriptome-wide association and colocalization to infer the effect of gene expression on disease severity, we find evidence that implicates multiple genes-including reduced expression of a membrane flippase (ATP11A), and increased expression of a mucin (MUC1)-in critical disease. Mendelian randomization provides evidence in support of causal roles for myeloid cell adhesion molecules (SELE, ICAM5 and CD209) and the coagulation factor F8, all of which are potentially druggable targets. Our results are broadly consistent with a multi-component model of COVID-19 pathophysiology, in which at least two distinct mechanisms can predispose to life-threatening disease: failure to control viral replication; or an enhanced tendency towards pulmonary inflammation and intravascular coagulation. We show that comparison between cases of critical illness and population controls is highly efficient for the detection of therapeutically relevant mechanisms of disease.

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.

Host genetics and infectious disease: new tools, insights and translational opportunities.

Understanding how human genetics influence infectious disease susceptibility offers the opportunity for new insights into pathogenesis, potential drug targets, risk stratification, response to therapy and vaccination. As new infectious diseases continue to emerge, together with growing levels of antimicrobial resistance and an increasing awareness of substantial differences between populations in genetic associations, the need for such work is expanding. In this Review, we illustrate how our understanding of the host-pathogen relationship is advancing through holistic approaches, describing current strategies to investigate the role of host genetic variation in established and emerging infections, including COVID-19, the need for wider application to diverse global populations mirroring the burden of disease, the impact of pathogen and vector genetic diversity and a broad array of immune and inflammation phenotypes that can be mapped as traits in health and disease. Insights from study of inborn errors of immunity and multi-omics profiling together with developments in analytical methods are further advancing our knowledge of this important area.

Genetic mechanisms of critical illness in COVID-19.

Host-mediated lung inflammation is present1, and drives mortality2, in the critical illness caused by coronavirus disease 2019 (COVID-19). Host genetic variants associated with critical illness may identify mechanistic targets for therapeutic development3. Here we report the results of the GenOMICC (Genetics Of Mortality In Critical Care) genome-wide association study in 2,244 critically ill patients with COVID-19 from 208 UK intensive care units. We have identified and replicated the following new genome-wide significant associations: on chromosome 12q24.13 (rs10735079, P = 1.65 × 10-8) in a gene cluster that encodes antiviral restriction enzyme activators (OAS1, OAS2 and OAS3); on chromosome 19p13.2 (rs74956615, P = 2.3 × 10-8) near the gene that encodes tyrosine kinase 2 (TYK2); on chromosome 19p13.3 (rs2109069, P = 3.98 ×  10-12) within the gene that encodes dipeptidyl peptidase 9 (DPP9); and on chromosome 21q22.1 (rs2236757, P = 4.99 × 10-8) in the interferon receptor gene IFNAR2. We identified potential targets for repurposing of licensed medications: using Mendelian randomization, we found evidence that low expression of IFNAR2, or high expression of TYK2, are associated with life-threatening disease; and transcriptome-wide association in lung tissue revealed that high expression of the monocyte-macrophage chemotactic receptor CCR2 is associated with severe COVID-19. Our results identify robust genetic signals relating to key host antiviral defence mechanisms and mediators of inflammatory organ damage in COVID-19. Both mechanisms may be amenable to targeted treatment with existing drugs. However, large-scale randomized clinical trials will be essential before any change to clinical practice.

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.

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.

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.

Priority index: database of genetic targets in immune-mediated disease.

We describe a comprehensive and unique database 'Priority index' (Pi; http://pi.well.ox.ac.uk) of prioritized genes encoding potential therapeutic targets that encompasses all major immune-mediated diseases. We provide targets at the gene level, each receiving a 5-star rating supported by: genomic evidence arising from disease genome-wide associations and functional immunogenomics, annotation evidence using ontologies restricted to genes with genomic evidence, and network evidence from protein interactions. Target genes often act together in related molecular pathways. The underlying Pi approach is unique in identifying a network of highly rated genes that mediate pathway crosstalk. In the Pi website, disease-centric pages are specially designed to enable the users to browse a complete list of prioritized genes and also a manageable list of nodal genes at the pathway crosstalk level; both switchable by clicks. Moreover, target genes are cross-referenced and supported using additional information, particularly regarding tractability, including druggable pockets viewed in 3D within protein structures. Target genes highly rated across diseases suggest drug repurposing opportunity, while genes in a particular disease reveal disease-specific targeting potential. To facilitate the ease of such utility, cross-disease comparisons involving multiple diseases are also supported. This facility, together with the faceted search, enhances integrative mining of the Pi resource to accelerate early-stage therapeutic target identification and validation leveraging human genetics.

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.

Fine-Scale Genetic Structure in the United Arab Emirates Reflects Endogamous and Consanguineous Culture, Population History, and Geography.

The indigenous population of the United Arab Emirates (UAE) has a unique demographic and cultural history. Its tradition of endogamy and consanguinity is expected to produce genetic homogeneity and partitioning of gene pools while population movements and intercontinental trade are likely to have contributed to genetic diversity. Emiratis and neighboring populations of the Middle East have been underrepresented in the population genetics literature with few studies covering the broader genetic history of the Arabian Peninsula. Here, we genotyped 1,198 individuals from the seven Emirates using 1.7 million markers and by employing haplotype-based algorithms and admixture analyses, we reveal the fine-scale genetic structure of the Emirati population. Shared ancestry and gene flow with neighboring populations display their unique geographic position while increased intra- versus inter-Emirati kinship and sharing of uniparental haplogroups, reflect the endogamous and consanguineous cultural traditions of the Emirates and their tribes.

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.

Identifying collagen VI as a target of fibrotic diseases regulated by CREBBP/EP300.

Fibrotic diseases remain a major cause of morbidity and mortality, yet there are few effective therapies. The underlying pathology of all fibrotic conditions is the activity of myofibroblasts. Using cells from freshly excised disease tissue from patients with Dupuytren's disease (DD), a localized fibrotic disorder of the palm, we sought to identify new therapeutic targets for fibrotic disease. We hypothesized that the persistent activity of myofibroblasts in fibrotic diseases might involve epigenetic modifications. Using a validated genetics-led target prioritization algorithm (Pi) of genome wide association studies (GWAS) data and a broad screen of epigenetic inhibitors, we found that the acetyltransferase CREBBP/EP300 is a major regulator of contractility and extracellular matrix production via control of H3K27 acetylation at the profibrotic genes, ACTA2 and COL1A1 Genomic analysis revealed that EP300 is highly enriched at enhancers associated with genes involved in multiple profibrotic pathways, and broad transcriptomic and proteomic profiling of CREBBP/EP300 inhibition by the chemical probe SGC-CBP30 identified collagen VI (Col VI) as a prominent downstream regulator of myofibroblast activity. Targeted Col VI knockdown results in significant decrease in profibrotic functions, including myofibroblast contractile force, extracellular matrix (ECM) production, chemotaxis, and wound healing. Further evidence for Col VI as a major determinant of fibrosis is its abundant expression within Dupuytren's nodules and also in the fibrotic foci of idiopathic pulmonary fibrosis (IPF). Thus, Col VI may represent a tractable therapeutic target across a range of fibrotic disorders.