Repeated Omicron exposures override ancestral SARS-CoV-2 immune imprinting.

The continuing emergence of SARS-CoV-2 variants highlights the need to update COVID-19 vaccine compositions. However, immune imprinting induced by vaccination based on the ancestral (hereafter referred to as WT) strain would compromise the antibody response to Omicron-based boosters1-5. Vaccination strategies to counter immune imprinting are critically needed. Here we investigated the degree and dynamics of immune imprinting in mouse models and human cohorts, especially focusing on the role of repeated Omicron stimulation. In mice, the efficacy of single Omicron boosting is heavily limited when using variants that are antigenically distinct from WT-such as the XBB variant-and this concerning situation could be mitigated by a second Omicron booster. Similarly, in humans, repeated Omicron infections could alleviate WT vaccination-induced immune imprinting and generate broad neutralization responses in both plasma and nasal mucosa. Notably, deep mutational scanning-based epitope characterization of 781 receptor-binding domain (RBD)-targeting monoclonal antibodies isolated from repeated Omicron infection revealed that double Omicron exposure could induce a large proportion of matured Omicron-specific antibodies that have distinct RBD epitopes to WT-induced antibodies. Consequently, immune imprinting was largely mitigated, and the bias towards non-neutralizing epitopes observed in single Omicron exposures was restored. On the basis of the deep mutational scanning profiles, we identified evolution hotspots of XBB.1.5 RBD and demonstrated that these mutations could further boost the immune-evasion capability of XBB.1.5 while maintaining high ACE2-binding affinity. Our findings suggest that the WT component should be abandoned when updating COVID-19 vaccines, and individuals without prior Omicron exposure should receive two updated vaccine boosters.

Imprinted SARS-CoV-2 humoral immunity induces convergent Omicron RBD evolution.

Continuous evolution of Omicron has led to a rapid and simultaneous emergence of numerous variants that display growth advantages over BA.5 (ref. 1). Despite their divergent evolutionary courses, mutations on their receptor-binding domain (RBD) converge on several hotspots. The driving force and destination of such sudden convergent evolution and its effect on humoral immunity remain unclear. Here we demonstrate that these convergent mutations can cause evasion of neutralizing antibody drugs and convalescent plasma, including those from BA.5 breakthrough infection, while maintaining sufficient ACE2-binding capability. BQ.1.1.10 (BQ.1.1 + Y144del), BA.4.6.3, XBB and CH.1.1 are the most antibody-evasive strains tested. To delineate the origin of the convergent evolution, we determined the escape mutation profiles and neutralization activity of monoclonal antibodies isolated from individuals who had BA.2 and BA.5 breakthrough infections2,3. Owing to humoral immune imprinting, BA.2 and especially BA.5 breakthrough infection reduced the diversity of the neutralizing antibody binding sites and increased proportions of non-neutralizing antibody clones, which, in turn, focused humoral immune pressure and promoted convergent evolution in the RBD. Moreover, we show that the convergent RBD mutations could be accurately inferred by deep mutational scanning profiles4,5, and the evolution trends of BA.2.75 and BA.5 subvariants could be well foreseen through constructed convergent pseudovirus mutants. These results suggest that current herd immunity and BA.5 vaccine boosters may not efficiently prevent the infection of Omicron convergent variants.

Omicron escapes the majority of existing SARS-CoV-2 neutralizing antibodies.

The SARS-CoV-2 B.1.1.529 (Omicron) variant contains 15 mutations of the receptor-binding domain (RBD). How Omicron evades RBD-targeted neutralizing antibodies requires immediate investigation. Here we use high-throughput yeast display screening1,2 to determine the profiles of RBD escaping mutations for 247 human anti-RBD neutralizing antibodies and show that the neutralizing antibodies can be classified by unsupervised clustering into six epitope groups (A-F)-a grouping that is highly concordant with knowledge-based structural classifications3-5. Various single mutations of Omicron can impair neutralizing antibodies of different epitope groups. Specifically, neutralizing antibodies in groups A-D, the epitopes of which overlap with the ACE2-binding motif, are largely escaped by K417N, G446S, E484A and Q493R. Antibodies in group E (for example, S309)6 and group F (for example, CR3022)7, which often exhibit broad sarbecovirus neutralizing activity, are less affected by Omicron, but a subset of neutralizing antibodies are still escaped by G339D, N440K and S371L. Furthermore, Omicron pseudovirus neutralization showed that neutralizing antibodies that sustained single mutations could also be escaped, owing to multiple synergetic mutations on their epitopes. In total, over 85% of the tested neutralizing antibodies were escaped by Omicron. With regard to neutralizing-antibody-based drugs, the neutralization potency of LY-CoV016, LY-CoV555, REGN10933, REGN10987, AZD1061, AZD8895 and BRII-196 was greatly undermined by Omicron, whereas VIR-7831 and DXP-604 still functioned at a reduced efficacy. Together, our data suggest that infection with Omicron would result in considerable humoral immune evasion, and that neutralizing antibodies targeting the sarbecovirus conserved region will remain most effective. Our results inform the development of antibody-based drugs and vaccines against Omicron and future variants.

BA.2.12.1, BA.4 and BA.5 escape antibodies elicited by Omicron infection.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) Omicron sublineages BA.2.12.1, BA.4 and BA.5 exhibit higher transmissibility than the BA.2 lineage1. The receptor binding and immune-evasion capability of these recently emerged variants require immediate investigation. Here, coupled with structural comparisons of the spike proteins, we show that BA.2.12.1, BA.4 and BA.5 (BA.4 and BA.5 are hereafter referred collectively to as BA.4/BA.5) exhibit similar binding affinities to BA.2 for the angiotensin-converting enzyme 2 (ACE2) receptor. Of note, BA.2.12.1 and BA.4/BA.5 display increased evasion of neutralizing antibodies compared with BA.2 against plasma from triple-vaccinated individuals or from individuals who developed a BA.1 infection after vaccination. To delineate the underlying antibody-evasion mechanism, we determined the escape mutation profiles2, epitope distribution3 and Omicron-neutralization efficiency of 1,640 neutralizing antibodies directed against the receptor-binding domain of the viral spike protein, including 614 antibodies isolated from people who had recovered from BA.1 infection. BA.1 infection after vaccination predominantly recalls humoral immune memory directed against ancestral (hereafter referred to as wild-type (WT)) SARS-CoV-2 spike protein. The resulting elicited antibodies could neutralize both WT SARS-CoV-2 and BA.1 and are enriched on epitopes on spike that do not bind ACE2. However, most of these cross-reactive neutralizing antibodies are evaded by spike mutants L452Q, L452R and F486V. BA.1 infection can also induce new clones of BA.1-specific antibodies that potently neutralize BA.1. Nevertheless, these neutralizing antibodies are largely evaded by BA.2 and BA.4/BA.5 owing to D405N and F486V mutations, and react weakly to pre-Omicron variants, exhibiting narrow neutralization breadths. The therapeutic neutralizing antibodies bebtelovimab4 and cilgavimab5 can effectively neutralize BA.2.12.1 and BA.4/BA.5, whereas the S371F, D405N and R408S mutations undermine most broadly sarbecovirus-neutralizing antibodies. Together, our results indicate that Omicron may evolve mutations to evade the humoral immunity elicited by BA.1 infection, suggesting that BA.1-derived vaccine boosters may not achieve broad-spectrum protection against new Omicron variants.

Simultaneous single-cell three-dimensional genome and gene expression profiling uncovers dynamic enhancer connectivity underlying olfactory receptor choice.

The simultaneous measurement of three-dimensional (3D) genome structure and gene expression of individual cells is critical for understanding a genome's structure-function relationship, yet this is challenging for existing methods. Here we present 'Linking mRNA to Chromatin Architecture (LiMCA)', which jointly profiles the 3D genome and transcriptome with exceptional sensitivity and from low-input materials. Combining LiMCA and our high-resolution scATAC-seq assay, METATAC, we successfully characterized chromatin accessibility, as well as paired 3D genome structures and gene expression information, of individual developing olfactory sensory neurons. We expanded the repertoire of known olfactory receptor (OR) enhancers and discovered unexpected rules of their dynamics: OR genes and their enhancers are most accessible during early differentiation. Furthermore, we revealed the dynamic spatial relationship between ORs and enhancers behind stepwise OR expression. These findings offer valuable insights into how 3D connectivity of ORs and enhancers dynamically orchestrate the 'one neuron-one receptor' selection process.

Convergent evolution of SARS-CoV-2 XBB lineages on receptor-binding domain 455-456 synergistically enhances antibody evasion and ACE2 binding.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) XBB lineages have achieved dominance worldwide and keep on evolving. Convergent evolution of XBB lineages on the receptor-binding domain (RBD) L455F and F456L is observed, resulting in variants with substantial growth advantages, such as EG.5, FL.1.5.1, XBB.1.5.70, and HK.3. Here, we show that neutralizing antibody (NAb) evasion drives the convergent evolution of F456L, while the epistatic shift caused by F456L enables the subsequent convergence of L455F through ACE2 binding enhancement and further immune evasion. L455F and F456L evade RBD-targeting Class 1 public NAbs, reducing the neutralization efficacy of XBB breakthrough infection (BTI) and reinfection convalescent plasma. Importantly, L455F single substitution significantly dampens receptor binding; however, the combination of L455F and F456L forms an adjacent residue flipping, which leads to enhanced NAbs resistance and ACE2 binding affinity. The perturbed receptor-binding mode leads to the exceptional ACE2 binding and NAb evasion, as revealed by structural analyses. Our results indicate the evolution flexibility contributed by epistasis cannot be underestimated, and the evolution potential of SARS-CoV-2 RBD remains high.

Highly sensitive single-cell chromatin accessibility assay and transcriptome coassay with METATAC.

Recent advances in single-cell assay for transposase accessible chromatin using sequencing (scATAC-seq) and its coassays have transformed the field of single-cell epigenomics and transcriptomics. However, the low detection efficiency of current methods has limited our understanding of the true complexity of chromatin accessibility and its relationship with gene expression in single cells. Here, we report a high-sensitivity scATAC-seq method, termed multiplexed end-tagging amplification of transposase accessible chromatin (METATAC), which detects a large number of accessible sites per cell and is compatible with automation. Our high detectability and statistical framework allowed precise linking of enhancers to promoters without merging single cells. We systematically investigated allele-specific accessibility in the mouse cerebral cortex, revealing allele-specific accessibility of promotors of certain imprinted genes but biallelic accessibility of their enhancers. Finally, we combined METATAC with our high-sensitivity single-cell RNA sequencing (scRNA-seq) method, multiple annealing and looping based amplification cycles for digital transcriptomics (MALBAC-DT), to develop a joint ATAC-RNA assay, termed METATAC and MALBAC-DT coassay by sequencing (M2C-seq). M2C-seq achieved significant improvements for both ATAC and RNA compared with previous methods, with consistent performance across cell lines and early mouse embryos.

Antigenicity assessment of SARS-CoV-2 saltation variant BA.2.87.1.

The recent emergence of a SARS-CoV-2 saltation variant, BA.2.87.1, which features 65 spike mutations relative to BA.2, has attracted worldwide attention. In this study, we elucidate the antigenic characteristics and immune evasion capability of BA.2.87.1. Our findings reveal that BA.2.87.1 is more susceptible to XBB-induced humoral immunity compared to JN.1. Notably, BA.2.87.1 lacks critical escaping mutations in the receptor binding domain (RBD) thus allowing various classes of neutralizing antibodies (NAbs) that were escaped by XBB or BA.2.86 subvariants to neutralize BA.2.87.1, although the deletions in the N-terminal domain (NTD), specifically 15-23del and 136-146del, compensate for the resistance to humoral immunity. Interestingly, several neutralizing antibody drugs have been found to restore their efficacy against BA.2.87.1, including SA58, REGN-10933 and COV2-2196. Hence, our results suggest that BA.2.87.1 may not become widespread until it acquires multiple RBD mutations to achieve sufficient immune evasion comparable to that of JN.1.

Immune evasion and ACE2 binding affinity contribute to SARS-CoV-2 evolution.

Mutations in the SARS-CoV-2 genome could confer resistance to pre-existing antibodies and/or increased transmissibility. The recently emerged Omicron subvariants exhibit a strong tendency for immune evasion, suggesting adaptive evolution. However, because previous studies have been limited to specific lineages or subsets of mutations, the overall evolutionary trajectory of SARS-CoV-2 and the underlying driving forces are still not fully understood. Here we analysed all open-access SARS-CoV-2 genomes (up to November 2022) and correlated the mutation incidence and fitness changes with the impacts of mutations on immune evasion and ACE2 binding affinity. Our results show that the Omicron lineage had an accelerated mutation rate in the RBD region, while the mutation incidence in other genomic regions did not change dramatically over time. Mutations in the RBD region exhibited a lineage-specific pattern and tended to become more aggregated over time, and the mutation incidence was positively correlated with the strength of antibody pressure. Additionally, mutation incidence was positively correlated with changes in ACE2 binding affinity, but with a lower correlation coefficient than with immune evasion. In contrast, the effect of mutations on fitness was more closely correlated with changes in ACE2 binding affinity than with immune evasion. Our findings suggest that immune evasion and ACE2 binding affinity play significant and diverse roles in the evolution of SARS-CoV-2.

Rational identification of potent and broad sarbecovirus-neutralizing antibody cocktails from SARS convalescents.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) Omicron sublineages have escaped most receptor-binding domain (RBD)-targeting therapeutic neutralizing antibodies (NAbs), which proves that previous NAb drug screening strategies are deficient against the fast-evolving SARS-CoV-2. Better broad NAb drug candidate selection methods are needed. Here, we describe a rational approach for identifying RBD-targeting broad SARS-CoV-2 NAb cocktails. Based on high-throughput epitope determination, we propose that broad NAb drugs should target non-immunodominant RBD epitopes to avoid herd-immunity-directed escape mutations. Also, their interacting antigen residues should focus on sarbecovirus conserved sites and associate with critical viral functions, making the antibody-escaping mutations less likely to appear. Following these criteria, a featured non-competing antibody cocktail, SA55+SA58, is identified from a large collection of broad sarbecovirus NAbs isolated from SARS-CoV-2-vaccinated SARS convalescents. SA55+SA58 potently neutralizes ACE2-utilizing sarbecoviruses, including circulating Omicron variants, and could serve as broad SARS-CoV-2 prophylactics to offer long-term protection, especially for individuals who are immunocompromised or with high-risk comorbidities.

Characterization of the enhanced infectivity and antibody evasion of Omicron BA.2.75.

Recently emerged SARS-CoV-2 Omicron subvariant, BA.2.75, displayed a growth advantage over circulating BA.2.38, BA.2.76, and BA.5 in India. However, the underlying mechanisms for enhanced infectivity, especially compared with BA.5, remain unclear. Here, we show that BA.2.75 exhibits substantially higher affinity for host receptor angiotensin-converting enzyme 2 (ACE2) than BA.5 and other variants. Structural analyses of BA.2.75 spike shows its decreased thermostability and increased frequency of the receptor binding domain (RBD) in the "up" conformation under acidic conditions, suggesting enhanced low-pH-endosomal cell entry. Relative to BA.4/BA.5, BA.2.75 exhibits reduced evasion of humoral immunity from BA.1/BA.2 breakthrough-infection convalescent plasma but greater evasion of Delta breakthrough-infection convalescent plasma. BA.5 breakthrough-infection plasma also exhibits weaker neutralization against BA.2.75 than BA.5, mainly due to BA.2.75's distinct neutralizing antibody (NAb) escape pattern. Antibody therapeutics Evusheld and Bebtelovimab remain effective against BA.2.75. These results suggest BA.2.75 may prevail after BA.4/BA.5, and its increased receptor-binding capability could support further immune-evasive mutations.