Towards increased accuracy and reproducibility in SARS-CoV-2 next generation sequence analysis for public health surveillance

During the COVID-19 pandemic, SARS-CoV-2 surveillance efforts integrated genome sequencing of clinical samples to identify emergent viral variants and to support rapid experimental examination of genome-informed vaccine and therapeutic designs. Given the broad range of methods applied to generate new viral genomes, it is critical that consensus and variant calling tools yield consistent results across disparate pipelines. Here we examine the impact of sequencing technologies (Illumina and Oxford Nanopore) and 7 different downstream bioinformatic protocols on SARS-CoV-2 variant calling as part of the NIH Accelerating COVID-19 Therapeutic Interventions and Vaccines (ACTIV) Tracking Resistance and Coronavirus Evolution (TRACE) initiative, a public-private partnership established to address the COVID-19 outbreak. Our results indicate that bioinformatic workflows can yield consensus genomes with different single nucleotide polymorphisms, insertions, and/or deletions even when using the same raw sequence input datasets. We introduce the use of a specific suite of parameters and protocols that greatly improves the agreement among pipelines developed by diverse organizations. Such consistency among bioinformatic pipelines is fundamental to SARS-CoV-2 and future pathogen surveillance efforts. The application of analysis standards is necessary to more accurately document phylogenomic trends and support data-driven public health responses.

Imprinted antibody responses against SARS-CoV-2 Omicron sublineages

SARS-CoV-2 Omicron sublineages carry distinct spike mutations and represent an antigenic shift resulting in escape from antibodies induced by previous infection or vaccination. We show that hybrid immunity or vaccine boosters result in potent plasma neutralizing activity against Omicron BA.1 and BA.2 and that breakthrough infections, but not vaccination-only, induce neutralizing activity in the nasal mucosa. Consistent with immunological imprinting, most antibodies derived from memory B cells or plasma cells of Omicron breakthrough cases cross-react with the Wuhan-Hu-1, BA.1 and BA.2 receptor-binding domains whereas Omicron primary infections elicit B cells of narrow specificity. While most clinical antibodies have reduced neutralization of Omicron, we identified an ultrapotent pan-variant antibody, that is unaffected by any Omicron lineage spike mutations and is a strong candidate for clinical development.

Broadly neutralizing antibodies overcome SARS-CoV-2 Omicron antigenic shift

The recently emerged SARS-CoV-2 Omicron variant harbors 37 amino acid substitutions in the spike (S) protein, 15 of which are in the receptor-binding domain (RBD), thereby raising concerns about the effectiveness of available vaccines and antibody therapeutics. Here, we show that the Omicron RBD binds to human ACE2 with enhanced affinity relative to the Wuhan-Hu-1 RBD and acquires binding to mouse ACE2. Severe reductions of plasma neutralizing activity were observed against Omicron compared to the ancestral pseudovirus for vaccinated and convalescent individuals. Most (26 out of 29) receptor-binding motif (RBM)-directed monoclonal antibodies (mAbs) lost in vitro neutralizing activity against Omicron, with only three mAbs, including the ACE2-mimicking S2K146 mAb1, retaining unaltered potency. Furthermore, a fraction of broadly neutralizing sarbecovirus mAbs recognizing antigenic sites outside the RBM, including sotrovimab2, S2X2593 and S2H974, neutralized Omicron. The magnitude of Omicron-mediated immune evasion and the acquisition of binding to mouse ACE2 mark a major SARS-CoV-2 mutational shift. Broadly neutralizing sarbecovirus mAbs recognizing epitopes conserved among SARS-CoV-2 variants and other sarbecoviruses may prove key to controlling the ongoing pandemic and future zoonotic spillovers.

Antibody therapy reverses biological signatures of COVID-19 progression

Understanding who is at risk of progression to severe COVID-19 is key to effective treatment. We studied correlates of disease severity in the COMET-ICE clinical trial that randomized 1:1 to placebo or to sotrovimab, a monoclonal antibody for the treatment of SARS-CoV-2 infection. Several laboratory parameters identified study participants at greater risk of severe disease, including a high neutrophil-lymphocyte ratio (NLR), a negative SARS-CoV-2 serologic test and whole blood transcriptome profiles. Sotrovimab treatment in these groups was associated with normalization of NLR and the transcriptomic profile, and with a decrease of viral RNA in nasopharyngeal samples. Transcriptomics provided the most sensitive detection of participants who would go on to be hospitalized or die. To facilitate timely measurement, we identified a 10-gene signature with similar predictive accuracy. In summary, we identified markers of risk for disease progression and demonstrated that normalization of these parameters occurs with antibody treatment of established infection.

Predicting the mutational drivers of future SARS-CoV-2 variants of concern

SARS-CoV-2 evolution threatens vaccine- and natural infection-derived immunity, and the efficacy of therapeutic antibodies. Herein we sought to predict Spike amino acid changes that could contribute to future variants of concern. We tested the importance of features comprising epidemiology, evolution, immunology, and neural network-based protein sequence modeling. This resulted in identification of the primary biological drivers of SARS-CoV-2 intra-pandemic evolution. We found evidence that resistance to population-level host immunity has increasingly shaped SARS-CoV-2 evolution over time. We identified with high accuracy mutations that will spread, at up to four months in advance, across different phases of the pandemic. Behavior of the model was consistent with a plausible causal structure wherein epidemiological variables integrate the effects of diverse and shifting drivers of viral fitness. We applied our model to forecast mutations that will spread in the future, and characterize how these mutations affect the binding of therapeutic antibodies. These findings demonstrate that it is possible to forecast the driver mutations that could appear in emerging SARS-CoV-2 variants of concern. This modeling approach may be applied to any pathogen with genomic surveillance data, and so may address other rapidly evolving pathogens such as influenza, and unknown future pandemic viruses.

A human antibody that broadly neutralizes betacoronaviruses protects against SARS-CoV-2 by blocking the fusion machinery

The repeated spillovers of {beta}-coronaviruses in humans along with the rapid emergence of SARS-CoV-2 escape variants highlight the need to develop broad coronavirus therapeutics and vaccines. Five monoclonal antibodies (mAbs) were isolated from COVID-19 convalescent individuals and found to cross-react with multiple {beta}-coronavirus spike (S) glycoproteins by targeting the stem helix. One of these mAbs, S2P6, cross-reacts with more than twenty human and animal {beta}-coronavirus S glycoproteins and broadly neutralizes SARS-CoV-2 and pseudotyped viruses from the sarbecovirus, merbecovirus and embecovirus subgenera. Structural and functional studies delineate the molecular basis of S2P6 cross-reactivity and broad neutralization and indicate that this mAb blocks viral entry by inhibiting membrane fusion. S2P6 protects hamsters challenged with SARS-CoV-2 (including the B.1.351 variant of concern) through direct viral neutralization and Fc-mediated effector functions. Serological and B cell repertoire analyses indicate that antibodies targeting the stem helix are found in some convalescent donors and vaccinees but are predominantly of narrow specificity. Germline reversion of the identified cross-reactive mAbs revealed that their unmutated ancestors are specific for the endemic OC43 or HKU1 viruses and acquired enhanced affinity and breadth through somatic mutations. These data demonstrate that conserved epitopes in the coronavirus fusion machinery can be targeted by protective antibodies and provide a framework for structure-guided design of pan-{beta}-coronavirus vaccines eliciting broad protection.

Antibodies to the SARS-CoV-2 receptor-binding domain that maximize breadth and resistance to viral escape

An ideal anti-SARS-CoV-2 antibody would resist viral escape1-3, have activity against diverse SARS-related coronaviruses4-7, and be highly protective through viral neutralization8-11 and effector functions12,13. Understanding how these properties relate to each other and vary across epitopes would aid development of antibody therapeutics and guide vaccine design. Here, we comprehensively characterize escape, breadth, and potency across a panel of SARS-CoV-2 antibodies targeting the receptor-binding domain (RBD), including S3094, the parental antibody of the late-stage clinical antibody VIR-7831. We observe a tradeoff between SARS-CoV-2 in vitro neutralization potency and breadth of binding across SARS-related coronaviruses. Nevertheless, we identify several neutralizing antibodies with exceptional breadth and resistance to escape, including a new antibody (S2H97) that binds with high affinity to all SARS-related coronavirus clades via a unique RBD epitope centered on residue E516. S2H97 and other escape-resistant antibodies have high binding affinity and target functionally constrained RBD residues. We find that antibodies targeting the ACE2 receptor binding motif (RBM) typically have poor breadth and are readily escaped by mutations despite high neutralization potency, but we identify one potent RBM antibody (S2E12) with breadth across sarbecoviruses closely related to SARS-CoV-2 and with a high barrier to viral escape. These data highlight functional diversity among antibodies targeting the RBD and identify epitopes and features to prioritize for antibody and vaccine development against the current and potential future pandemics.

Structural basis for broad sarbecovirus neutralization by a human monoclonal antibody

The recent emergence of SARS-CoV-2 variants of concern (VOC) and the recurrent spillovers of coronaviruses in the human population highlight the need for broadly neutralizing antibodies that are not affected by the ongoing antigenic drift and that can prevent or treat future zoonotic infections. Here, we describe a human monoclonal antibody (mAb), designated S2x259, recognizing a highly conserved cryptic receptor-binding domain (RBD) epitope and cross-reacting with spikes from all sarbecovirus clades. S2x259 broadly neutralizes spike-mediated entry of SARS-CoV-2 including the B.1.1.7, B.1.351, P.1 and B.1.427/B.1.429 VOC, as well as a wide spectrum of human and zoonotic sarbecoviruses through inhibition of ACE2 binding to the RBD. Furthermore, deep-mutational scanning and in vitro escape selection experiments demonstrate that S2x259 possesses a remarkably high barrier to the emergence of resistance mutants. We show that prophylactic administration of S2x259 protects Syrian hamsters against challenges with the prototypic SARS-CoV-2 and the B.1.351 variant, suggesting this mAb is a promising candidate for the prevention and treatment of emergent VOC and zoonotic infections. Our data unveil a key antigenic site targeted by broadly-neutralizing antibodies and will guide the design of pan-sarbecovirus vaccines.

SARS-CoV-2 immune evasion by variant B.1.427/B.1.429

SARS-CoV-2 entry is mediated by the spike (S) glycoprotein which contains the receptor-binding domain (RBD) and the N-terminal domain (NTD) as the two main targets of neutralizing antibodies (Abs). A novel variant of concern (VOC) named CAL.20C (B.1.427/B.1.429) was originally detected in California and is currently spreading throughout the US and 29 additional countries. It is unclear whether antibody responses to SARS-CoV-2 infection or to the prototypic Wuhan-1 isolate-based vaccines will be impacted by the three B.1.427/B.1.429 S mutations: S13I, W152C and L452R. Here, we assessed neutralizing Ab responses following natural infection or mRNA vaccination using pseudoviruses expressing the wildtype or the B.1.427/B.1.429 S protein. Plasma from vaccinated or convalescent individuals exhibited neutralizing titers, which were reduced 3-6 fold against the B.1.427/B.1.429 variant relative to wildtype pseudoviruses. The RBD L452R mutation reduced or abolished neutralizing activity of 14 out of 35 RBD-specific monoclonal antibodies (mAbs), including three clinical-stage mAbs. Furthermore, we observed a complete loss of B.1.427/B.1.429 neutralization for a panel of mAbs targeting the N-terminal domain due to a large structural rearrangement of the NTD antigenic supersite involving an S13I-mediated shift of the signal peptide cleavage site. These data warrant closer monitoring of signal peptide variants and their involvement in immune evasion and show that Abs directed to the NTD impose a selection pressure driving SARS-CoV-2 viral evolution through conventional and unconventional escape mechanisms.

The dual function monoclonal antibodies VIR-7831 and VIR-7832 demonstrate potent in vitro and in vivo activity against SARS-CoV-2

Sotrovimab (VIR-7831) and VIR-7832 are dual action monoclonal antibodies (mAbs) targeting the spike glycoprotein of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). Sotrovimab and VIR-7832 were derived from a parent antibody (S309) isolated from memory B cells of a 2003 severe acute respiratory syndrome coronavirus (SARS-CoV) survivor. Both mAbs contain an "LS" mutation in the Fc region to prolong serum half-life. In addition, VIR-7832 encodes an Fc GAALIE mutation that has been shown previously to evoke CD8+ T-cells in the context of an in vivo viral respiratory infection. Sotrovimab and VIR-7832 neutralize wild-type and variant pseudotyped viruses and authentic virus in vitro. In addition, they retain activity against monoclonal antibody resistance mutations conferring reduced susceptibility to previously authorized mAbs. The sotrovimab/VIR-7832 epitope continues to be highly conserved among circulating sequences consistent with the high barrier to resistance observed in vitro. Furthermore, both mAbs can recruit effector mechanisms in vitro that may contribute to clinical efficacy via elimination of infected host cells. In vitro studies with these mAbs demonstrated no enhancement of infection. In a Syrian Golden hamster proof-of concept wildtype SARS-CoV-2 infection model, animals treated with sotrovimab had less weight loss, and significantly decreased total viral load and infectious virus levels in the lung compared to a control mAb. Taken together, these data indicate that sotrovimab and VIR-7832 are key agents in the fight against COVID-19.

N-terminal domain antigenic mapping reveals a site of vulnerability for SARS-CoV-2

SARS-CoV-2 entry into host cells is orchestrated by the spike (S) glycoprotein that contains an immunodominant receptor-binding domain (RBD) targeted by the largest fraction of neutralizing antibodies (Abs) in COVID-19 patient plasma. Little is known about neutralizing Abs binding to epitopes outside the RBD and their contribution to protection. Here, we describe 41 human monoclonal Abs (mAbs) derived from memory B cells, which recognize the SARS-CoV-2 S N-terminal domain (NTD) and show that a subset of them neutralize SARS-CoV-2 ultrapotently. We define an antigenic map of the SARS-CoV-2 NTD and identify a supersite recognized by all known NTD-specific neutralizing mAbs. These mAbs inhibit cell-to-cell fusion, activate effector functions, and protect Syrian hamsters from SARS-CoV-2 challenge. SARS-CoV-2 variants, including the 501Y.V2 and B.1.1.7 lineages, harbor frequent mutations localized in the NTD supersite suggesting ongoing selective pressure and the importance of NTD-specific neutralizing mAbs to protective immunity.

Transfer transcriptomic signatures for infectious diseases

The modulation of the transcriptome is among the earliest responses to infection, and vaccination. However, defining transcriptome signatures of disease is challenging because logistic, technical and cost factors limit the size and representativeness of samples in clinical studies. These limitations lead to poor performance of signatures when applied to new datasets or varying study settings. Using a novel approach, we leverage existing transcriptomic signatures as classifiers in unseen datasets from prospective studies, with the goal of predicting individual outcomes. Machine learning allowed the identification of sets of genes, which we name transfer transcriptomic signatures, that are predictive across diverse datasets and/or species (rhesus to humans) and that are also suggestive of activated pathways and cell type composition. We demonstrate the usefulness of transfer signatures in two use cases: progression of latent to active tuberculosis, and severity of COVID-19 and influenza A H1N1 infection. The broad significance of our work lies in the concept that a small set of archetypal human immunophenotypes, captured by transfer signatures, can explain a larger set of responses to diverse diseases.