Selection Analysis Identifies Clusters of Unusual Mutational Changes in Omicron Lineage BA.1 That Likely Impact Spike Function

Predicting human and viral protein variants affecting COVID-19 susceptibility and repurposing therapeutics

The COVID-19 disease is an ongoing global health concern. Although vaccination provides some protection, people are still susceptible to re-infection. Ostensibly, certain populations or clinical groups may be more vulnerable. Factors causing these differences are unclear and whilst socioeconomic and cultural differences are likely to be important, human genetic factors could influence susceptibility. Experimental studies indicate SARS-CoV-2 uses innate immune suppression as a strategy to speed-up entry and replication into the host cell. Therefore, it is necessary to understand the impact of variants in immunity-associated human proteins on susceptibility to COVID-19. In this work, we analysed missense coding variants in several SARS-CoV-2 proteins and their human protein interactors that could enhance binding affinity to SARS-CoV-2. We curated a dataset of 19 SARS-CoV-2: human protein 3D-complexes, from the experimentally determined structures in the Protein Data Bank and models built using AlphaFold2-multimer, and analysed the impact of missense variants occurring in the protein-protein interface region. We analysed 468 missense variants from human proteins and 212 variants from SARS-CoV-2 proteins and computationally predicted their impacts on binding affinities for the human viral protein complexes. We predicted a total of 26 affinity-enhancing variants from 13 human proteins implicated in increased binding affinity to SARS-CoV-2. These include key-immunity associated genes (TOMM70, ISG15, IFIH1, IFIT2, RPS3, PALS1, NUP98, AXL, ARF6, TRIMM, TRIM25) as well as important spike receptors (KREMEN1, AXL and ACE2). We report both common (e.g., Y13N in IFIH1) and rare variants in these proteins and discuss their likely structural and functional impact, using information on known and predicted functional sites. Potential mechanisms associated with immune suppression implicated by these variants are discussed. Occurrence of certain predicted affinity-enhancing variants should be monitored as they could lead to increased susceptibility and reduced immune response to SARS-CoV-2 infection in individuals/populations carrying them. Our analyses aid in understanding the potential impact of genetic variation in immunity-associated proteins on COVID-19 susceptibility and help guide drug-repurposing strategies.

AOC: Analysis of Orthologous Collections - an application for the characterization of natural selection in protein-coding sequences

Motivation

Modern molecular sequence analysis increasingly relies on automated and robust software tools for interpretation, annotation, and biological insight. The Analysis of Orthologous Collections (AOC) application automates the identification of genomic sites and species/lineages influenced by natural selection in coding sequence analysis. AOC quantifies different types of selection: negative, diversifying or directional positive, or differential selection between groups of branches. We include all steps necessary to go from unaligned homologous sequences to complete results and interactive visualizations that are designed to aid in the useful interpretation and contextualization.

Results

We are motivated by a desire to make evolutionary analyses as simple as possible, and to close the disparity in the literature between genes which draw a significant amount of interest and those that are largely overlooked and underexplored. We believe that such underappreciated and understudied genetic datasets can hold rich biological information and offer substantial insights into the diverse patterns and processes of evolution, especially if domain experts are able to perform the analyses themselves.

Availability and implementation

A Snakemake [Mölder et al., 2021] application implementation is publicly available on GitHub at https://github.com/aglucaci/AnalysisOfOrthologousCollections and is accompanied by software documentation and a tutorial.

Persistence of SARS-CoV-2 infection and viral intra- and inter-host evolution in COVID-19 hospitalized patients

Severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) persistence in COVID-19 patients could play a key role in the emergence of variants of concern. The rapid intra-host evolution of SARS-CoV-2 may result in an increased transmissibility, immune and therapeutic escape which could be a direct consequence of COVID-19 epidemic currents. In this context, a longitudinal retrospective study on eight consecutive COVID-19 patients with persistent SARS-CoV-2 infection, from January 2022 to March 2023, was conducted. To characterize the intra- and inter-host viral evolution, whole genome sequencing and phylogenetic analysis were performed on nasopharyngeal samples collected at different time points. Phylogenetic reconstruction revealed an accelerated SARS-CoV-2 intra-host evolution and emergence of antigenically divergent variants. The Bayesian inference and principal coordinate analysis analysis showed a host-based genomic structuring among antigenically divergent variants, that might reflect the positive effect of containment practices, within the critical hospital area. All longitudinal antigenically divergent isolates shared a wide range of amino acidic (aa) changes, particularly in the Spike (S) glycoprotein, that increased viral transmissibility (K417N, S477N, N501Y and Q498R), enhanced infectivity (R346T, S373P, R408S, T478K, Q498R, Y505H, D614G, H655Y, N679K and P681H), caused host immune escape (S371L, S375F, T376A, K417N, and K444T/R) and displayed partial or complete resistance to treatments (G339D, R346K/T, S371F/L, S375F, T376A, D405N, N440K, G446S, N460K, E484A, F486V, Q493R, G496S and Q498R). These results suggest that multiple novel variants which emerge in the patient during persistent infection, might spread to another individual and continue to evolve. A pro-active genomic surveillance of persistent SARS-CoV-2 infected patients is recommended to identify genetically divergent lineages before their diffusion.

Mapping the Evolutionary Space of SARS-CoV-2 Variants to Anticipate Emergence of Subvariants Resistant to COVID-19 Therapeutics

New sublineages of SARS-CoV-2 variants-of-concern (VOCs) continuously emerge with mutations in the spike glycoprotein. In most cases, the sublineage-defining mutations vary between the VOCs. It is unclear whether these differences reflect lineage-specific likelihoods for mutations at each spike position or the stochastic nature of their appearance. Here we show that SARS-CoV-2 lineages have distinct evolutionary spaces (a probabilistic definition of the sequence states that can be occupied by expanding virus subpopulations). This space can be accurately inferred from the patterns of amino acid variability at the whole-protein level. Robust networks of co-variable sites identify the highest-likelihood mutations in new VOC sublineages and predict remarkably well the emergence of subvariants with resistance mutations to COVID-19 therapeutics. Our studies reveal the contribution of low frequency variant patterns at heterologous sites across the protein to accurate prediction of the changes at each position of interest.

Alum/CpG adjuvant promotes immunogenicity of inactivated SARS-CoV-2 Omicron vaccine through enhanced humoral and cellular immunity

The SARS-CoV-2 Omicron variant, which was classified as a variant of concern (VOC) by the World Health Organization on 26 November 2021, has attracted worldwide attention for its high transmissibility and immune evasion ability. The existing COVID-19 vaccine has been shown to be less effective in preventing Omicron variant infection and symptomatic infection, which brings new challenges to vaccine development and application. Here, we evaluated the immunogenicity and safety of an Omicron variant COVID-19 inactivated vaccine containing aluminum and CpG adjuvants in a variety of animal models. The results showed that the vaccine candidate could induce high levels of neutralizing antibodies against the Omicron variant virus and binding antibodies, and significantly promoted cellular immune response. Meanwhile, the vaccine candidate was safe. Therefore, it provided more foundation for the development of aluminum and CpG as a combination adjuvant in human vaccines.

Correlates of Breakthrough SARS-CoV-2 Infections in People with HIV: Results from the CIHR CTN 328 Study

COVID-19 breakthrough infection (BTI) can occur despite vaccination. Using a multi-centre, prospective, observational Canadian cohort of people with HIV (PWH) receiving ≥2 COVID-19 vaccines, we compared the SARS-CoV-2 spike (S) and receptor-binding domain (RBD)-specific IgG levels 3 and 6 months post second dose, as well as 1 month post third dose, in PWH with and without BTI. BTI was defined as positivity based on self-report measures (data up to last study visit) or IgG data (up to 1 month post dose 3). The self-report measures were based on their symptoms and either a positive PCR or rapid antigen test. The analysis was restricted to persons without previous COVID-19 infection. Persons without BTI remained COVID-19-naïve until ≥3 months following the third dose. Of 289 participants, 92 developed BTI (31.5 infections per 100 person-years). The median days between last vaccination and BTI was 128 (IQR 67, 176), with the most cases occurring between the third and fourth dose (n = 59), corresponding to the Omicron wave. In analyses adjusted for age, sex, race, multimorbidity, hypertension, chronic kidney disease, diabetes and obesity, a lower IgG S/RBD (log10 BAU/mL) at 1 month post dose 3 was significantly associated with BTI, suggesting that a lower IgG level at this time point may predict BTI in this cohort of PWH.

SARS-CoV-2 specific antibody response after an mRNA vaccine as the third dose: Homologous versus heterologous boost

The aim of this study was to evaluate immunogenicity and longevity of the humoral immune response within six months after the homologous (BNT162b2/BNT162b2) or heterologous (BBIBP-CorV/BNT162b2) third dose, and to assess breakthrough infections among vaccinees during the Omicron wave in Serbia. Serum samples were analyzed at four timepoints: five months after the primary series; three weeks, three months, and six months after the boost. IgG antibodies against the receptor-binding domain of the spike protein were detected using enzyme-linked fluorescence assay. Both homologous (n = 55) and heterologous group (n = 36) showed a highly significant increase in antibody concentrations (p < 0.001) three weeks after the boost. A moderate inverse correlation between the age of recipients and the antibody levels at three weeks post-boost was observed in the homologous group (p = 0.02, r = -0.37), while the same correlation was not significant for heterologous group (p = 0.55, r = -0.15). Heterologous group had significantly higher antibody concentrations than homologous group at three weeks (Median 851.4(IQR 766.6-894.1); 784.3(676.9-847.4); p = 0.03) and three months post-boost (766.6(534.8-798.9); 496.8(361.6-664.0); p < 0.001). However, a significant decline in antibody response over time was noted for both strategies. The overall incidence of breakthrough cases was estimated at 36.36% (20/55) for homologous, and 16.67% (6/36) for heterologous group, but none of them required hospitalization. Although observed incidence in the homologous group was more than double when compared to the heterologous group, this difference was not statistically significant, most likely due to the small sample size. In conclusion, waning immunity after inactivated vaccine can be recovered by BNT162b2 heterologous boost regardless of the age of recipients, and both boost strategies induced potent humoral immune response and protection against severe COVID-19 during the Omicron wave. However, as the observed incidence of breakthrough infections was higher in the homologous group, although non-significant, this finding could indicate an advantage of heterologous approach.

An updated review of epidemiological characteristics, immune escape, and therapeutic advances of SARS-CoV-2 Omicron XBB.1.5 and other mutants

The rapid evolution of Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) has led to the emergence of new variants with different genetic profiles, with important implications for public health. The continued emergence of new variants with unique genetic features and potential changes in biological properties poses significant challenges to public health strategies, vaccine development, and therapeutic interventions. Omicron variants have attracted particular attention due to their rapid spread and numerous mutations in key viral proteins. This review aims to provide an updated and comprehensive assessment of the epidemiological characteristics, immune escape potential, and therapeutic advances of the SARS-CoV-2 Omicron XBB.1.5 variant, as well as other variants.

Reconstructing relative transmission rates in Bayesian phylodynamics: Two-fold transmission advantage of Omicron in Berlin, Germany during December 2021

Phylodynamic methods have lately played a key role in understanding the spread of infectious diseases. During the coronavirus disease (COVID-19) pandemic, large scale genomic surveillance has further increased the potential of dynamic inference from viral genomes. With the continual emergence of novel severe acute respiratory syndrome coronavirus type 2 (SARS-CoV-2) variants, explicitly allowing transmission rate differences between simultaneously circulating variants in phylodynamic inference is crucial. In this study, we present and empirically validate an extension to the BEAST2 package birth-death skyline model (BDSKY), BDSKY[Formula: see text], which introduces a scaling factor for the transmission rate between independent, jointly inferred trees. In an extensive simulation study, we show that BDSKY[Formula: see text] robustly infers the relative transmission rates under different epidemic scenarios. Using publicly available genome data of SARS-CoV-2, we apply BDSKY[Formula: see text] to quantify the transmission advantage of the Omicron over the Delta variant in Berlin, Germany. We find the overall transmission rate of Omicron to be scaled by a factor of two with pronounced variation between the individual clusters of each variant. These results quantify the transmission advantage of Omicron over the previously circulating Delta variant, in a crucial period of pre-established non-pharmaceutical interventions. By inferring variant- as well as cluster-specific transmission rate scaling factors, we show the differences in transmission dynamics for each variant. This highlights the importance of incorporating lineage-specific transmission differences in phylodynamic inference.

Raman Fingerprints of SARS-CoV-2 Omicron Subvariants: Molecular Roots of Virological Characteristics and Evolutionary Directions

The latest RNA genomic mutation of SARS-CoV-2 virus, termed the Omicron variant, has generated a stream of highly contagious and antibody-resistant strains, which in turn led to classifying Omicron as a variant of concern. We systematically collected Raman spectra from six Omicron subvariants available in Japan (i.e., BA.1.18, BA.2, BA.4, BA.5, XE, and BA.2.75) and applied machine-learning algorithms to decrypt their structural characteristics at the molecular scale. Unique Raman fingerprints of sulfur-containing amino acid rotamers, RNA purines and pyrimidines, tyrosine phenol ring configurations, and secondary protein structures clearly differentiated the six Omicron subvariants. These spectral characteristics, which were linked to infectiousness, transmissibility, and propensity for immune evasion, revealed evolutionary motifs to be compared with the outputs of genomic studies. The availability of a Raman "metabolomic snapshot", which was then translated into a barcode to enable a prompt subvariant identification, opened the way to rationalize in real-time SARS-CoV-2 activity and variability. As a proof of concept, we applied the Raman barcode procedure to a nasal swab sample retrieved from a SARS-CoV-2 patient and identified its Omicron subvariant by coupling a commercially available magnetic bead technology with our newly developed Raman analyses.

Multiple mutations of SARS-CoV-2 Omicron BA.2 variant orchestrate its virological characteristics

IMPORTANCE

Most studies investigating the characteristics of emerging SARS-CoV-2 variants have been focusing on mutations in the spike proteins that affect viral infectivity, fusogenicity, and pathogenicity. However, few studies have addressed how naturally occurring mutations in the non-spike regions of the SARS-CoV-2 genome impact virological properties. In this study, we proved that multiple SARS-CoV-2 Omicron BA.2 mutations, one in the spike protein and another downstream of the spike gene, orchestrally characterize this variant, shedding light on the importance of Omicron BA.2 mutations out of the spike protein.

Computational studies on the interaction of Omicron subvariants (BA.1, BA.2, and BA.3) with ACE2 and polyphenols

INTRODUCTION

The SARS-CoV-2 Omicron variant BA.2 is spreading widely across the globe. The World Health Organization (WHO) designated BA.2 as a variant of concern due to its high transmission rate and pathogenicity. To elucidate the structural changes caused by mutations, we conducted a comparative analysis of BA.2 with variants BA.1 and BA.3.

OBJECTIVE

In the present study, we aimed to investigate the interactions of the spike glycoprotein receptor-binding domain (SGp RBD) of Omicron variants BA.1, BA.2, and BA.3 with the human receptor hACE2. Further, a library of 233 polyphenols was screened by molecular docking with the SGp RBDs of Omicron variants BA.1, BA.2, and BA.3.

METHODS

Protein-protein and protein-ligand molecular docking simulations were performed with AutoDock Vina and the ClusPro 2.0 server, respectively. The protein-ligand interactions were evaluated by BIOVIA Discovery Studio and ChimeraX 1.4. The molecular dynamics simulations for 100 ns were performed using GROMACS 2021.

RESULTS

Compared to other variants of concern, the structural changes in Omicron caused by mutations at key positions improved its ability to cause infection. Despite multiple mutations, many important polyphenols bind effectively at the RBDs of Omicron variants, with the required pharmacokinetic and ADME features and obeying the Lipinski rule.

CONCLUSION

Even though Omicron variants have multiple mutations and their transmission rate is relatively high, the computed binding affinities of lead polyphenols like epigallocatechin-3-O-gallate (EGCG) and luteolin-7-O-glucuronide (L7G) indicate that traditional medicines and proper immunity booster diets may be useful in the long-term fight against SARS-CoV-2.

Rapid, tunable, and multiplexed detection of RNA using convective array PCR

Detection of RNA targets is typically achieved through RT-qPCR or RNAseq. RT-qPCR is rapid but limited in number and complexity of targets detected, while RNAseq is high-throughput but takes multiple days. We demonstrate simultaneous amplification and detection of 28 distinct RNA targets from a single unsplit purified RNA sample in under 40 minutes using our convective array PCR (caPCR) technology. We integrate tunable strand displacement probes into caPCR to allow detection of RNA species with programmable sequence selectivity for either a single, perfectly matched target sequence or for targets with up to 2 single-nucleotide variants within the probe-binding regions. Tunable probes allow for robust detection of desired RNA species against high homology background sequences and robust detection of RNA species with significant sequence diversity due to community-acquired mutations. As a proof-of-concept, we experimentally demonstrated detection of 7 human coronaviruses and 7 key variants of concern of SARS-CoV-2 in a single assay.

Alternative molecular and genomic strategies to provide a rapid response to alerts concerning the introduction of new emerging SARS-CoV-2 variants: the Omicron alert

During the COVID-19 pandemic, several SARS-CoV-2 variants of concern (VOCs) of particular relevance emerged. Early detection of VOCs entering a country is essential to control spread. The alert triggered by the first suspected case of the Omicron variant in Spain in a traveler arriving from South Africa in November 2021 provided a unique opportunity to evaluate four different methodological strategies tailored to rapid identification of Omicron. The different approaches were designed to respond to the different technical resources available in different settings. First, we used melting probes in RT-PCR to determine the presence of four Omicron signatures (K417N, E484A, P681H, and absence of L452R): three probes showed deviations in temperature (Tm) values relative to the reference codons (E484K-15.8°C, P681H-5.2°C, and L452R-7.2°C) and one maintained the reference value (K417N). The deviation in Tm of P681H suggested the presence of the characteristic Omicron N679K mutation in the probe hybridization region; these data pointed to the presence of Omicron alleles. Second, the presence of 29 of the 33 characteristic single nucleotide polymorphisms (SNPs) in the Omicron variant S-gene was identified by Sanger sequencing of nine amplicons. The final two strategies involved identification of 47 of the 50 non-synonymous and indel mutations attributed to Omicron by rapid nanopore whole genome sequencing (WGS) and by Illumina WGS technology. These strategies enabled us to pre-assign the first Omicron case in Spain with high certainty 2 h after receipt of RNA and to confirm it genomically 3 h later, so that the Public Health authorities could be rapidly notified. IMPORTANCE The study presents different experimental alternatives to identify new variants of concern (VOCs) of SARS-CoV-2 entering a certain population. Early detection of a new VOC is crucial for surveillance and control of spread. The objective is to provide laboratories with tools adapted to their resource capabilities that offer a sufficient level of resolution to rule out, confirm, or pre-assign the presence of a suspected VOC. The study describes four different techniques that were applied simultaneously to the first suspected Omicron case in Spain, highlighting the level of resolution and response time achieved in each case. These techniques are based on the detection of mutations in the S-gene of the virus that can easily adapt to potential emerging variants. The results of the study allow any laboratory to prepare for new alerts of SARS-CoV-2 VOCs.

Natural selection shapes the evolution of SARS-CoV-2 Omicron in Bangladesh

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has evolved to give rise to a highly transmissive and immune-escaping variant of concern, known as Omicron. Many aspects of the evolution of SARS-CoV-2 and the driving forces behind the ongoing Omicron outbreaks remain unclear. Substitution at the receptor-binding domain (RBD) in the spike protein is one of the primary strategies of SARS-CoV-2 Omicron to hinder recognition by the host angiotensin-converting enzyme 2 (ACE2) receptor and avoid antibody-dependent defense activation. Here, we scanned for adaptive evolution within the SARS-CoV-2 Omicron genomes reported from Bangladesh in the public database GISAID (www.gisaid.org; dated 2 April 2023). The ratio of the non-synonymous (Ka) to synonymous (Ks) nucleotide substitution rate, denoted as ω, is an indicator of the selection pressure acting on protein-coding genes. A higher proportion of non-synonymous to synonymous substitutions (Ka/Ks or ω > 1) indicates positive selection, while Ka/Ks or ω near zero indicates purifying selection. An equal amount of non-synonymous and synonymous substitutions (Ka/Ks or ω = 1) refers to neutrally evolving sites. We found evidence of adaptive evolution within the spike (S) gene of SARS-CoV-2 Omicron isolated from Bangladesh. In total, 22 codon sites of the S gene displayed a signature of positive selection. The data also highlighted that the receptor-binding motif within the RBD of the spike glycoprotein is a hotspot of adaptive evolution, where many of the codons had ω > 1. Some of these adaptive sites at the RBD of the spike protein are known to be associated with increased viral fitness. The M gene and ORF6 have also experienced positive selection. These results suggest that although purifying selection is the dominant evolutionary force, positive Darwinian selection also plays a vital role in shaping the evolution of SARS-CoV-2 Omicron in Bangladesh.

The signature of SARS-CoV-2 evolution reflects selective pressures within human guts

In somatic cells, microRNAs (miRNAs) bind to the genomes of RNA viruses and influence their translation and replication. In London and Berlin samples represented in GISAID database, we traced severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) lineages and divided these sequenced in two groups, "Ancestral variants" and "Omicrons," and analyzed them through the prism of the tissue-specific binding between host miRNAs and viral messenger RNAs. We demonstrate a significant number of miRNA-binding sites in the NSP4 region of the SARS-CoV-2 genome, with evidence of evolutionary pressure within this region exerted by human intestinal miRNAs. Notably, in infected cells, NSP4 promotes the formation of double-membrane vesicles, which serve as the scaffolds for replication-transcriptional complexes and protect viral RNA from intracellular destruction. In 3 years of selection, the loss of many miRNA-binding sites in general and those within the NSP4 in particular has shaped the SARS-CoV-2 genomes. With that, the descendants of the BA.2 variants were promoted as dominant strains, which define current momentum of the pandemics.

Deep-learning-enabled protein-protein interaction analysis for prediction of SARS-CoV-2 infectivity and variant evolution

Host-pathogen interactions and pathogen evolution are underpinned by protein-protein interactions between viral and host proteins. An understanding of how viral variants affect protein-protein binding is important for predicting viral-host interactions, such as the emergence of new pathogenic SARS-CoV-2 variants. Here we propose an artificial intelligence-based framework called UniBind, in which proteins are represented as a graph at the residue and atom levels. UniBind integrates protein three-dimensional structure and binding affinity and is capable of multi-task learning for heterogeneous biological data integration. In systematic tests on benchmark datasets and further experimental validation, UniBind effectively and scalably predicted the effects of SARS-CoV-2 spike protein variants on their binding affinities to the human ACE2 receptor, as well as to SARS-CoV-2 neutralizing monoclonal antibodies. Furthermore, in a cross-species analysis, UniBind could be applied to predict host susceptibility to SARS-CoV-2 variants and to predict future viral variant evolutionary trends. This in silico approach has the potential to serve as an early warning system for problematic emerging SARS-CoV-2 variants, as well as to facilitate research on protein-protein interactions in general.

The effects of amino acid substitution of spike protein and genomic recombination on the evolution of SARS-CoV-2

Over three years' pandemic of 2019 novel coronavirus disease (COVID-19), multiple variants and novel subvariants have emerged successively, outcompeted earlier variants and become predominant. The sequential emergence of variants reflects the evolutionary process of mutation-selection-adaption of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). Amino acid substitution/insertion/deletion in the spike protein causes altered viral antigenicity, transmissibility, and pathogenicity of SARS-CoV-2. Early in the pandemic, D614G mutation conferred virus with advantages over previous variants and increased transmissibility, and it also laid a conservative background for subsequent substantial mutations. The role of genomic recombination in the evolution of SARS-CoV-2 raised increasing concern with the occurrence of novel recombinants such as Deltacron, XBB.1.5, XBB.1.9.1, and XBB.1.16 in the late phase of pandemic. Co-circulation of different variants and co-infection in immunocompromised patients accelerate the emergence of recombinants. Surveillance for SARS-CoV-2 genomic variations, particularly spike protein mutation and recombination, is essential to identify ongoing changes in the viral genome and antigenic epitopes and thus leads to the development of new vaccine strategies and interventions.

The recency and geographical origins of the bat viruses ancestral to SARS-CoV and SARS-CoV-2

The emergence of SARS-CoV in 2002 and SARS-CoV-2 in 2019 has led to increased sampling of related sarbecoviruses circulating primarily in horseshoe bats. These viruses undergo frequent recombination and exhibit spatial structuring across Asia. Employing recombination-aware phylogenetic inference on bat sarbecoviruses, we find that the closest-inferred bat virus ancestors of SARS-CoV and SARS-CoV-2 existed just ~1-3 years prior to their emergence in humans. Phylogeographic analyses examining the movement of related sarbecoviruses demonstrate that they traveled at similar rates to their horseshoe bat hosts and have been circulating for thousands of years in Asia. The closest-inferred bat virus ancestor of SARS-CoV likely circulated in western China, and that of SARS-CoV-2 likely circulated in a region comprising southwest China and northern Laos, both a substantial distance from where they emerged. This distance and recency indicate that the direct ancestors of SARS-CoV and SARS-CoV-2 could not have reached their respective sites of emergence via the bat reservoir alone. Our recombination-aware dating and phylogeographic analyses reveal a more accurate inference of evolutionary history than performing only whole-genome or single gene analyses. These results can guide future sampling efforts and demonstrate that viral genomic fragments extremely closely related to SARS-CoV and SARS-CoV-2 were circulating in horseshoe bats, confirming their importance as the reservoir species for SARS viruses.

Modeling of ACE2 and antibodies bound to SARS-CoV-2 provides insights into infectivity and immune evasion

Given the COVID-19 pandemic, there is interest in understanding ligand-receptor features and targeted antibody-binding attributes against emerging SARS-CoV-2 variants. Here, we developed a large-scale structure-based pipeline for analysis of protein-protein interactions regulating SARS-CoV-2 immune evasion. First, we generated computed structural models of the Spike protein of 3 SARS-CoV-2 variants (B.1.1.529, BA.2.12.1, and BA.5) bound either to a native receptor (ACE2) or to a large panel of targeted ligands (n = 282), which included neutralizing or therapeutic monoclonal antibodies. Moreover, by using the Barnes classification, we noted an overall loss of interfacial interactions (with gain of new interactions in certain cases) at the receptor-binding domain (RBD) mediated by substituted residues for neutralizing complexes in classes 1 and 2, whereas less destabilization was observed for classes 3 and 4. Finally, an experimental validation of predicted weakened therapeutic antibody binding was performed in a cell-based assay. Compared with the original Omicron variant (B.1.1.529), derivative variants featured progressive destabilization of antibody-RBD interfaces mediated by a larger set of substituted residues, thereby providing a molecular basis for immune evasion. This approach and findings provide a framework for rapidly and efficiently generating structural models for SARS-CoV-2 variants bound to ligands of mechanistic and therapeutic value.

Evolutionary Shortcuts via Multinucleotide Substitutions and Their Impact on Natural Selection Analyses

Inference and interpretation of evolutionary processes, in particular of the types and targets of natural selection affecting coding sequences, are critically influenced by the assumptions built into statistical models and tests. If certain aspects of the substitution process (even when they are not of direct interest) are presumed absent or are modeled with too crude of a simplification, estimates of key model parameters can become biased, often systematically, and lead to poor statistical performance. Previous work established that failing to accommodate multinucleotide (or multihit, MH) substitutions strongly biases dN/dS-based inference towards false-positive inferences of diversifying episodic selection, as does failing to model variation in the rate of synonymous substitution (SRV) among sites. Here, we develop an integrated analytical framework and software tools to simultaneously incorporate these sources of evolutionary complexity into selection analyses. We found that both MH and SRV are ubiquitous in empirical alignments, and incorporating them has a strong effect on whether or not positive selection is detected (1.4-fold reduction) and on the distributions of inferred evolutionary rates. With simulation studies, we show that this effect is not attributable to reduced statistical power caused by using a more complex model. After a detailed examination of 21 benchmark alignments and a new high-resolution analysis showing which parts of the alignment provide support for positive selection, we show that MH substitutions occurring along shorter branches in the tree explain a significant fraction of discrepant results in selection detection. Our results add to the growing body of literature which examines decades-old modeling assumptions (including MH) and finds them to be problematic for comparative genomic data analysis. Because multinucleotide substitutions have a significant impact on natural selection detection even at the level of an entire gene, we recommend that selection analyses of this type consider their inclusion as a matter of routine. To facilitate this procedure, we developed, implemented, and benchmarked a simple and well-performing model testing selection detection framework able to screen an alignment for positive selection with two biologically important confounding processes: site-to-site synonymous rate variation, and multinucleotide instantaneous substitutions.

An Early SARS-CoV-2 Omicron Outbreak in a Dormitory in Saint Petersburg, Russia

The Omicron variant of SARS-CoV-2 rapidly spread worldwide in late 2021-early 2022, displacing the previously prevalent Delta variant. Before 16 December 2021, community transmission had already been observed in tens of countries globally. However, in Russia, the majority of reported cases at that time had been sporadic and associated with travel. Here, we report an Omicron outbreak at a student dormitory in Saint Petersburg between 16-29 December 2021, which was the earliest known instance of a large-scale community transmission in Russia. Out of the 465 sampled residents of the dormitory, 180 (38.7%) tested PCR-positive. Among the 118 residents for whom the variant had been tested by whole-genome sequencing, 111 (94.1%) were found to carry the Omicron variant. Among these 111 residents, 60 (54.1%) were vaccinated or had reported a previous infection of COVID-19. Phylogenetic analysis confirmed that the outbreak was caused by a single introduction of the BA.1.1 sub-lineage of the Omicron variant. The dormitory-derived clade constituted a significant proportion of BA.1.1 samples in Saint Petersburg and has spread to other regions of Russia and even to other countries. The rapid spread of the Omicron variant in a population with preexisting immunity to previous variants underlines its propensity for immune evasion.

Web Resources for SARS-CoV-2 Genomic Database, Annotation, Analysis and Variant Tracking

The SARS-CoV-2 genomic data continue to grow, providing valuable information for researchers and public health officials. Genomic analysis of these data sheds light on the transmission and evolution of the virus. To aid in SARS-CoV-2 genomic analysis, many web resources have been developed to store, collate, analyze, and visualize the genomic data. This review summarizes web resources used for the SARS-CoV-2 genomic epidemiology, covering data management and sharing, genomic annotation, analysis, and variant tracking. The challenges and further expectations for these web resources are also discussed. Finally, we highlight the importance and need for continued development and improvement of related web resources to effectively track the spread and understand the evolution of the virus.

Convergent evolution of SARS-CoV-2 Omicron subvariants leading to the emergence of BQ.1.1 variant

In late 2022, various Omicron subvariants emerged and cocirculated worldwide. These variants convergently acquired amino acid substitutions at critical residues in the spike protein, including residues R346, K444, L452, N460, and F486. Here, we characterize the convergent evolution of Omicron subvariants and the properties of one recent lineage of concern, BQ.1.1. Our phylogenetic analysis suggests that these five substitutions are recurrently acquired, particularly in younger Omicron lineages. Epidemic dynamics modelling suggests that the five substitutions increase viral fitness, and a large proportion of the fitness variation within Omicron lineages can be explained by these substitutions. Compared to BA.5, BQ.1.1 evades breakthrough BA.2 and BA.5 infection sera more efficiently, as demonstrated by neutralization assays. The pathogenicity of BQ.1.1 in hamsters is lower than that of BA.5. Our multiscale investigations illuminate the evolutionary rules governing the convergent evolution for known Omicron lineages as of 2022.

Genomic Surveillance of SARS-CoV-2 Variants in the Dominican Republic and Emergence of a Local Lineage

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is an RNA virus that evolves over time, leading to new variants. In the current study, we assessed the genomic epidemiology of SARS-CoV-2 in the Dominican Republic. A total of 1149 SARS-CoV-2 complete genome nucleotide sequences from samples collected between March 2020 and mid-February 2022 in the Dominican Republic were obtained from the Global Initiative on Sharing All Influenza Data (GISAID) database. Phylogenetic relationships and evolution rates were analyzed using the maximum likelihood method and the Bayesian Markov chain Monte Carlo (MCMC) approach. The genotyping details (lineages) were obtained using the Pangolin web application. In addition, the web tools Coronapp, and Genome Detective Viral Tools, among others, were used to monitor epidemiological characteristics. Our results show that the most frequent non-synonymous mutation over the study period was D614G. Of the 1149 samples, 870 (75.74%) were classified into 8 relevant variants according to Pangolin/Scorpio. The first Variants Being Monitored (VBM) were detected in December 2020. Meanwhile, in 2021, the variants of concern Delta and Omicron were identified. The mean mutation rate was estimated to be 1.5523 × 10^-3 (95% HPD: 1.2358 × 10^-3, 1.8635 × 10^-3) nucleotide substitutions per site. We also report the emergence of an autochthonous SARS-CoV-2 lineage, B.1.575.2, that circulated from October 2021 to January 2022, in co-circulation with the variants of concern Delta and Omicron. The impact of B.1.575.2 in the Dominican Republic was minimal, but it then expanded rapidly in Spain. A better understanding of viral evolution and genomic surveillance data will help to inform strategies to mitigate the impact on public health.

Coarse-Grained Molecular Simulations and Ensemble-Based Mutational Profiling of Protein Stability in the Different Functional Forms of the SARS-CoV-2 Spike Trimers: Balancing Stability and Adaptability in BA.1, BA.2 and BA.2.75 Variants

Evolutionary and functional studies have suggested that the emergence of Omicron variants can be determined by multiple fitness tradeoffs including immune escape, binding affinity, conformational plasticity, protein stability, and allosteric modulation. In this study, we embarked on a systematic comparative analysis of the conformational dynamics, electrostatics, protein stability, and allostery in the different functional states of spike trimers for BA.1, BA.2, and BA.2.75 variants. Using efficient and accurate coarse-grained simulations and atomistic reconstruction of the ensembles, we examined the conformational dynamics of the spike trimers that agree with the recent functional studies, suggesting that BA.2.75 trimers are the most stable among these variants. A systematic mutational scanning of the inter-protomer interfaces in the spike trimers revealed a group of conserved structural stability hotspots that play a key role in the modulation of functional dynamics and are also involved in the inter-protomer couplings through local contacts and interaction networks with the Omicron mutational sites. The results of mutational scanning provided evidence that BA.2.75 trimers are more stable than BA.2 and comparable in stability to the BA.1 variant. Using dynamic network modeling of the S Omicron BA.1, BA.2, and BA.2.75 trimers, we showed that the key network mediators of allosteric interactions are associated with the major stability hotspots that are interconnected along potential communication pathways. The network analysis of the BA.1, BA.2, and BA.2.75 trimers suggested that the increased thermodynamic stability of the BA.2.75 variant may be linked with the organization and modularity of the residue interaction network that allows for allosteric communications between structural stability hotspots and Omicron mutational sites. This study provided a plausible rationale for a mechanism in which Omicron mutations may evolve by targeting vulnerable sites of conformational adaptability to elicit immune escape while maintaining their control on balancing protein stability and functional fitness through robust allosteric communications with the stability hotspots.

Evolution of increased positive charge on the SARS-CoV-2 spike protein may be adaptation to human transmission

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) continues to evolve and infect individuals. The exterior surface of the SARS-CoV-2 virion is dominated by the spike protein, and the current work examined spike protein biochemical features that have changed during the 3 years in which SARS-CoV-2 has infected humans. Our analysis identified a striking change in spike protein charge, from -8.3 in the original Lineage A and B viruses to -1.26 in most of the current Omicron viruses. We conclude that in addition to immune selection pressure, the evolution of SARS-CoV-2 has also altered viral spike protein biochemical properties, which may influence virion survival and promote transmission. Future vaccine and therapeutic development should also exploit and target these biochemical properties.

Natural selection differences detected in key protein domains between non-pathogenic and pathogenic feline coronavirus phenotypes

Feline coronaviruses (FCoVs) commonly cause mild enteric infections in felines worldwide (termed feline enteric coronavirus [FECV]), with around 12 per cent developing into deadly feline infectious peritonitis (FIP; feline infectious peritonitis virus [FIPV]). Genomic differences between FECV and FIPV have been reported, yet the putative genotypic basis of the highly pathogenic phenotype remains unclear. Here, we used state-of-the-art molecular evolutionary genetic statistical techniques to identify and compare differences in natural selection pressure between FECV and FIPV sequences, as well as to identify FIPV- and FECV-specific signals of positive selection. We analyzed full-length FCoV protein coding genes thought to contain mutations associated with FIPV (Spike, ORF3abc, and ORF7ab). We identified two sites exhibiting differences in natural selection pressure between FECV and FIPV: one within the S1/S2 furin cleavage site (FCS) and the other within the fusion domain of Spike. We also found fifteen sites subject to positive selection associated with FIPV within Spike, eleven of which have not previously been suggested as possibly relevant to FIP development. These sites fall within Spike protein subdomains that participate in host cell receptor interaction, immune evasion, tropism shifts, host cellular entry, and viral escape. There were fourteen sites (twelve novel sites) within Spike under positive selection associated with the FECV phenotype, almost exclusively within the S1/S2 FCS and adjacent to C domain, along with a signal of relaxed selection in FIPV relative to FECV, suggesting that furin cleavage functionality may not be needed for FIPV. Positive selection inferred in ORF7b was associated with the FECV phenotype and included twenty-four positively selected sites, while ORF7b had signals of relaxed selection in FIPV. We found evidence of positive selection in ORF3c in FCoV-wide analyses, but no specific association with the FIPV or FECV phenotype. We hypothesize that some combination of mutations in FECV may contribute to FIP development, and that it is unlikely to be one singular 'switch' mutational event. This work expands our understanding of the complexities of FIP development and provides insights into how evolutionary forces may alter pathogenesis in coronavirus genomes.

Seasonal effects decouple SARS-CoV-2 haplotypes worldwide

Background: Variants of concern (VOCs) have been replacing each other during the still rampant COVID-19 pandemic. As a result, SARS-CoV-2 populations have evolved increasingly intricate constellations of mutations that often enhance transmissibility, disease severity, and other epidemiological characteristics. The origin and evolution of these constellations remain puzzling. Methods: Here we study the evolution of VOCs at the proteome level by analyzing about 12 million genomic sequences retrieved from GISAID on July 23, 2022. A total 183,276 mutations were identified and filtered with a relevancy heuristic. The prevalence of haplotypes and free-standing mutations was then tracked monthly in various latitude corridors of the world. Results: A chronology of 22 haplotypes defined three phases driven by protein flexibility-rigidity, environmental sensing, and immune escape. A network of haplotypes illustrated the recruitment and coalescence of mutations into major VOC constellations and seasonal effects of decoupling and loss. Protein interaction networks mediated by haplotypes predicted communications impacting the structure and function of proteins, showing the increasingly central role of molecular interactions involving the spike (S), nucleocapsid (N), and membrane (M) proteins. Haplotype markers either affected fusogenic regions while spreading along the sequence of the S-protein or clustered around binding domains. Modeling of protein structure with AlphaFold2 showed that VOC Omicron and one of its haplotypes were major contributors to the distortion of the M-protein endodomain, which behaves as a receptor of other structural proteins during virion assembly. Remarkably, VOC constellations acted cooperatively to balance the more extreme effects of individual haplotypes. Conclusions: Our study uncovers seasonal patterns of emergence and diversification occurring amid a highly dynamic evolutionary landscape of bursts and waves. The mapping of genetically-linked mutations to structures that sense environmental change with powerful ab initio modeling tools demonstrates the potential of deep-learning for COVID-19 predictive intelligence and therapeutic intervention.

Exploring and Learning the Universe of Protein Allostery Using Artificial Intelligence Augmented Biophysical and Computational Approaches

Allosteric mechanisms are commonly employed regulatory tools used by proteins to orchestrate complex biochemical processes and control communications in cells. The quantitative understanding and characterization of allosteric molecular events are among major challenges in modern biology and require integration of innovative computational experimental approaches to obtain atomistic-level knowledge of the allosteric states, interactions, and dynamic conformational landscapes. The growing body of computational and experimental studies empowered by emerging artificial intelligence (AI) technologies has opened up new paradigms for exploring and learning the universe of protein allostery from first principles. In this review we analyze recent developments in high-throughput deep mutational scanning of allosteric protein functions; applications and latest adaptations of Alpha-fold structural prediction methods for studies of protein dynamics and allostery; new frontiers in integrating machine learning and enhanced sampling techniques for characterization of allostery; and recent advances in structural biology approaches for studies of allosteric systems. We also highlight recent computational and experimental studies of the SARS-CoV-2 spike (S) proteins revealing an important and often hidden role of allosteric regulation driving functional conformational changes, binding interactions with the host receptor, and mutational escape mechanisms of S proteins which are critical for viral infection. We conclude with a summary and outlook of future directions suggesting that AI-augmented biophysical and computer simulation approaches are beginning to transform studies of protein allostery toward systematic characterization of allosteric landscapes, hidden allosteric states, and mechanisms which may bring about a new revolution in molecular biology and drug discovery.

Molecular Characterization of AZD7442 (Tixagevimab-Cilgavimab) Neutralization of SARS-CoV-2 Omicron Subvariants

Therapeutic anti-severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) monoclonal antibodies (MAbs) provide immunosuppressed and vulnerable populations with prophylactic and treatment interventions against coronavirus disease 2019 (COVID-19). AZD7442 (tixagevimab-cilgavimab) is a combination of extended-half-life neutralizing MAbs that bind to distinct epitopes on the receptor binding domain (RBD) of the SARS-CoV-2 spike protein. The Omicron variant of concern carries mutations at >35 positions in the spike protein and has undergone further genetic diversification since its emergence in November 2021. Here, we characterize the in vitro neutralization activity of AZD7442 toward major viral subvariants circulating worldwide during the first 9 months of the Omicron wave. BA.2 and its derived subvariants showed the highest susceptibility to AZD7442, while BA.1 and BA.1.1 showed a lower susceptibility. BA.4/BA.5 had a susceptibility level intermediate between BA.1 and BA.2. Mutagenesis of parental Omicron subvariant spike proteins was performed to establish a molecular model to describe the underlying determinants of neutralization by AZD7442 and its component MAbs. The concurrent mutation of residues at positions 446 and 493, located in the tixagevimab and cilgavimab binding sites, was sufficient to enhance in vitro susceptibility of BA.1 to AZD7442 and its component MAbs to levels similar to the Wuhan-Hu-1+D614G virus. AZD7442 maintained neutralization activity against all Omicron subvariants tested up to and including BA.5. The evolving nature of the SARS-CoV-2 pandemic warrants continuing real-time molecular surveillance and assessment of in vitro activity of MAbs used in prophylaxis against and the treatment of COVID-19. IMPORTANCE MAbs are key therapeutic options for COVID-19 prophylaxis and treatment in immunosuppressed and vulnerable populations. Due to the emergence of SARS-CoV-2 variants, including Omicron, it is vital to ensure that neutralization is maintained for MAb-based interventions. We studied the in vitro neutralization of AZD7442 (tixagevimab-cilgavimab), a cocktail of two long-acting MAbs targeting the SARS-CoV-2 spike protein, toward Omicron subvariants circulating from November 2021 to July 2022. AZD7442 neutralized major Omicron subvariants up to and including BA.5. The mechanism of action responsible for the lower in vitro susceptibility of BA.1 to AZD7442 was investigated using in vitro mutagenesis and molecular modeling. A combination of mutations at two spike protein positions, namely, 446 and 493, was sufficient to enhance BA.1 susceptibility to AZD7442 to levels similar to the Wuhan-Hu-1+D614G ancestral virus. The evolving nature of the SARS-CoV-2 pandemic warrants continuing real-time global molecular surveillance and mechanistic studies of therapeutic MAbs for COVID-19.

Effects of natural polymorphisms in SARS-CoV-2 RNA-dependent RNA polymerase on its activity and sensitivity to inhibitors in vitro

SARS-CoV-2 RNA-dependent RNA polymerase (RdRp) is the key enzyme required for viral replication and mRNA synthesis. RdRp is one of the most conserved viral proteins and a promising target for antiviral drugs and inhibitors. At the same time, analysis of public databases reveals multiple variants of SARS-CoV-2 genomes with substitutions in the catalytic RdRp subunit nsp12. Structural mapping of these mutations suggests that some of them may affect the interactions of nsp12 with its cofactors nsp7/nsp8 as well as with RNA substrates. We have obtained several mutations of these types and demonstrated that some of them decrease specific activity of RdRp in vitro, possibly by changing RdRp assembly and/or its interactions with RNA. Therefore, natural polymorphisms in RdRp may potentially affect viral replication. Furthermore, we have synthesized a series of polyphenol and diketoacid derivatives based on previously studied inhibitors of hepatitis C virus RdRp and found that several of them can inhibit SARS-CoV-2 RdRp. Tested mutations in RdRp do not have strong effects on the efficiency of inhibition. Further development of more efficient non-nucleoside inhibitors of SARS-CoV-2 RdRp should take into account the existence of multiple polymorphic variants of RdRp.

Insight into SARS-CoV-2 Omicron variant immune escape possibility and variant independent potential therapeutic opportunities

The Omicron, the latest variant of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), was first detected in November 2021 in Botswana, South Africa. Compared to other variants of SARS-CoV-2, the Omicron is the most highly mutated, with 50 mutations throughout the genome, most of which are in the spike (S) protein. These mutations may help the Omicron to evade host immunity against the vaccine. Epidemiological studies suggest that Omicron is highly infectious and spreads rapidly, but causes significantly less severe disease than the wild-type strain and the other variants of SARS-CoV-2. With the increased transmissibility and a higher rate of re-infection, Omicron has now become a dominant variant worldwide and is predicted to be able to evade vaccine-induced immunity. Several clinical studies using plasma samples from individuals receiving two doses of US Food and Drugs Administration (FDA)-approved COVID-19 vaccines have shown reduced humoral immune response against Omicron infection, but T cell-mediated immunity was well preserved. In fact, T cell-mediated immunity protects against severe disease, and thus the disease caused by Omicron remains mild. In this review, I surveyed the current status of Omicron variant mutations and mechanisms of immune response in the context of immune escape from COVID-19 vaccines. I also discuss the potential implications of therapeutic opportunities that are independent of SARS-CoV-2 variants, including Omicron. A better understanding of vaccine-induced immune responses and variant-independent therapeutic interventions that include potent antiviral, antioxidant, and anti-cytokine activities may pave the way to reducing Omicron-related COVID-19 complications, severity, and mortality. Collectively, these insights point to potential research gaps and will aid in the development of new-generation COVID-19 vaccines and antiviral drugs to combat Omicron, its sublineages, or upcoming new variants of SARS-CoV-2.

COVID-19: Challenges of Viral Variants

The COVID-19 pandemic has been accompanied by SARS-CoV-2 evolution and emergence of viral variants that have far exceeded initial expectations. Five major variants of concern (Alpha, Beta, Gamma, Delta, and Omicron) have emerged, each having both unique and overlapping amino acid substitutions that have affected transmissibility, disease severity, and susceptibility to natural or vaccine-induced immune responses and monoclonal antibodies. Several of the more recent variants appear to have evolved properties of immune evasion, particularly in cases of prolonged infection. Tracking of existing variants and surveillance for new variants are critical for an effective pandemic response.

Coordinated evolution at amino acid sites of SARS-CoV-2 spike

SARS-CoV-2 has adapted in a stepwise manner, with multiple beneficial mutations accumulating in a rapid succession at origins of VOCs, and the reasons for this are unclear. Here, we searched for coordinated evolution of amino acid sites in the spike protein of SARS-CoV-2. Specifically, we searched for concordantly evolving site pairs (CSPs) for which changes at one site were rapidly followed by changes at the other site in the same lineage. We detected 46 sites which formed 45 CSP. Sites in CSP were closer to each other in the protein structure than random pairs, indicating that concordant evolution has a functional basis. Notably, site pairs carrying lineage defining mutations of the four VOCs that circulated before May 2021 are enriched in CSPs. For the Alpha VOC, the enrichment is detected even if Alpha sequences are removed from analysis, indicating that VOC origin could have been facilitated by positive epistasis. Additionally, we detected nine discordantly evolving pairs of sites where mutations at one site unexpectedly rarely occurred on the background of a specific allele at another site, for example on the background of wild-type D at site 614 (four pairs) or derived Y at site 501 (three pairs). Our findings hint that positive epistasis between accumulating mutations could have delayed the assembly of advantageous combinations of mutations comprising at least some of the VOCs.

Clinical development of variant-adapted BNT162b2 COVID-19 vaccines: the early Omicron era

INTRODUCTION

The Omicron BA.1 variant of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) and subsequent sub-lineages exhibit partial escape from neutralizing antibodies elicited by vaccines containing or encoding wild-type spike protein. In response to the emergence of Omicron sub-lineages, variant-adapted vaccines that contain or encode for Omicron spike protein components have been developed.

AREAS COVERED

This review presents currently available clinical immunogenicity and safety data on Omicron variant-adapted versions of the BNT162b2 messenger RNA (mRNA) vaccine and summarizes the expected mechanism of action, and rationale for development, of these vaccines. In addition, challenges encountered during development and regulatory approval are discussed.

EXPERT OPINION

Omicron-adapted BNT162b2 vaccines provide a wider breadth and potentially more durable protection against Omicron sub-lineages and antigenically aligned variants when compared with the original vaccine. As SARS-CoV-2 continues to evolve, further vaccine updates may be required. To facilitate this, a globally harmonized regulatory process for the transition to updated vaccines is needed. Next-generation vaccine approaches may provide broader protection against future variants.

mRNA-1273 boost after BNT162b2 vaccination generates comparable SARS-CoV-2-specific functional responses in naïve and COVID-19-recovered individuals

Introduction

COVID-19 vaccines based on mRNA have represented a revolution in the biomedical research field. The initial two-dose vaccination schedule generates potent humoral and cellular responses, with a massive protective effect against severe COVID-19 and death. Months after this vaccination, levels of antibodies against SARS-CoV-2 waned, and this promoted the recommendation of a third vaccination dose.

Methods

We have performed an integral and longitudinal study of the immunological responses triggered by the booster mRNA-1273 vaccination, in a cohort of health workers previously vaccinated with two doses of the BNT162b2 vaccine at University Hospital La Paz located in Madrid, Spain. Circulating humoral responses and SARS-CoV-2-specific cellular reactions, after ex vivo restimulation of both T and B cells (cytokines production, proliferation, class switching), have been analyzed. Importantly, all along these studies, the analyses have been performed comparing naïve and subjects recovered from COVID-19, addressing the influence of a previous infection by SARS-CoV-2. Furthermore, as the injection of the third vaccination dose was contemporary to the rise of the Omicron BA.1 variant of concern, T- and B-cell-mediated cellular responses have been comparatively analyzed in response to this variant.

Results

All these analyses indicated that differential responses to vaccination due to a previous SARS-CoV-2 infection were balanced following the boost. The increase in circulating humoral responses due to this booster dropped after 6 months, whereas T-cell-mediated responses were more stable along the time. Finally, all the analyzed immunological features were dampened in response to the Omicron variant of concern, particularly late after the booster vaccination.

Conclusion

This work represents a follow-up longitudinal study for almost 1.5 years, analyzing in an integral manner the immunological responses triggered by the prime-boost mRNA-based vaccination schedule against COVID-19.

The SARS-CoV-2 spike S375F mutation characterizes the Omicron BA.1 variant

Recent studies have revealed the unique virological characteristics of Omicron, particularly those of its spike protein, such as less cleavage efficacy in cells, reduced ACE2 binding affinity, and poor fusogenicity. However, it remains unclear which mutation(s) determine these three virological characteristics of Omicron spike. Here, we show that these characteristics of the Omicron spike protein are determined by its receptor-binding domain. Of interest, molecular phylogenetic analysis revealed that acquisition of the spike S375F mutation was closely associated with the explosive spread of Omicron in the human population. We further elucidated that the F375 residue forms an interprotomer pi-pi interaction with the H505 residue of another protomer in the spike trimer, conferring the attenuated cleavage efficiency and fusogenicity of Omicron spike. Our data shed light on the evolutionary events underlying the emergence of Omicron at the molecular level.

A Glimpse on the Evolution of RNA Viruses: Implications and Lessons from SARS-CoV-2

RNA viruses are characterised by extremely high genetic variability due to fast replication, large population size, low fidelity, and (usually) a lack of proofreading mechanisms of RNA polymerases leading to high mutation rates. Furthermore, viral recombination and reassortment may act as a significant evolutionary force among viruses contributing to greater genetic diversity than obtainable by mutation alone. The above-mentioned properties allow for the rapid evolution of RNA viruses, which may result in difficulties in viral eradication, changes in virulence and pathogenicity, and lead to events such as cross-species transmissions, which are matters of great interest in the light of current severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) pandemics. In this review, we aim to explore the molecular mechanisms of the variability of viral RNA genomes, emphasising the evolutionary trajectory of SARS-CoV-2 and its variants. Furthermore, the causes and consequences of coronavirus variation are explored, along with theories on the origin of human coronaviruses and features of emergent RNA viruses in general. Finally, we summarise the current knowledge on the circulating variants of concern and highlight the many unknowns regarding SARS-CoV-2 pathogenesis.

Contributions of adaptation and purifying selection to SARS-CoV-2 evolution

Continued evolution and adaptation of SARS-CoV-2 has led to more transmissible and immune-evasive variants with profound impacts on the course of the pandemic. Here I analyze the evolution of the virus over 2.5 years since its emergence and estimate the rates of evolution for synonymous and non-synonymous changes separately for evolution within clades-well-defined monophyletic groups with gradual evolution-and for the pandemic overall. The rate of synonymous mutation is found to be around 6 changes per year. Synonymous rates within variants vary little from variant to variant and are compatible with the overall rate of 7 changes per year (or [Formula: see text] per year and codon). In contrast, the rate at which variants accumulate amino acid changes (non-synonymous mutations) was initially around 12-16 changes per year, but in 2021 and 2022 it dropped to 6-9 changes per year. The overall rate of non-synonymous evolution, that is across variants, is estimated to be about 26 amino acid changes per year (or [Formula: see text] per year and codon). This strong acceleration of the overall rate compared to within clade evolution indicates that the evolutionary process that gave rise to the different variants is qualitatively different from that in typical transmission chains and likely dominated by adaptive evolution. I further quantify the spectrum of mutations and purifying selection in different SARS-CoV-2 proteins and show that the massive global sampling of SARS-CoV-2 is sufficient to estimate site-specific fitness costs across the entire genome. Many accessory proteins evolve under limited evolutionary constraints with little short-term purifying selection. About half of the mutations in other proteins are strongly deleterious.

Sierra SARS-CoV-2 sequence and antiviral resistance analysis program

INTRODUCTION

Although most laboratories are capable of employing established protocols to perform full-genome SARS-CoV-2 sequencing, many are unable to assess sequence quality, select appropriate mutation-detection thresholds, or report on the potential clinical significance of mutations in the targets of antiviral therapy METHODS: We describe the technical aspects and benchmark the performance of Sierra SARS-CoV-2, a program designed to perform these functions on user-submitted FASTQ and FASTA sequence files and lists of Spike mutations. Sierra SARS-CoV-2 indicates which sequences contain an unexpectedly large number of unusual mutations and which mutations are associated with reduced susceptibility to clinical stage mAbs, the RdRP inhibitor remdesivir, or the Mpro inhibitor nirmatrelvir RESULTS: To assess the performance of Sierra SARS-CoV-2 on FASTQ files, we applied it to 600 representative FASTQ sequences and compared the results to the COVID-19 EDGE program. To assess its performance on FASTA files, we applied it to nearly one million representative FASTA sequences and compared the results to the GISAID mutation annotation. To assess its performance on mutations lists, we applied it to 13,578 distinct Spike RBD mutation patterns and showed that exactly or partially matching annotations were available for 88% of patterns CONCLUSION: Sierra SARS-CoV-2 leverages previously published data to improve the quality control of submitted viral genomic data and to provide functional annotation on the impact of mutations in the targets of antiviral SARS-CoV-2 therapy. The program can be found at https://covdb.stanford.edu/sierra/sars2/ and its source code at https://github.com/hivdb/sierra-sars2.

Binding and neutralizing abilities of antibodies towards SARS-CoV-2 S2 domain

Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) variants have been reported to be resistant to several neutralizing antibodies (NAbs) targeting Receptor Binding Domain (RBD) and N Terminal Domain (NTD) of spike (S) protein and thus inducing immune escape. However, fewer studies were carried out to investigate the neutralizing ability of S2-specific antibodies. In this research, 10 monoclonal antibodies (mAbs) targeting SARS-CoV-2 S2 subunit were generated from Coronavirus Disease 2019 (COVID-19) convalescent patients by phage display technology and molecular cloning technology. The binding activity of these S2-mAbs toward SARS-CoV-2 S, SARS-CoV-2 S2, SARS-CoV-2 RBD, SARS-CoV-2 NTD, severe acute respiratory syndrome coronavirus (SARS-CoV) S, SARS-CoV S2 and Middle East Respiratory Syndrome Coronavirus (MERS-CoV) S proteins were evaluated by enzyme-linked immunosorbent assay (ELISA). Their neutralizing potency toward SARS-CoV-2 wild-type (WT), B.1.1.7, B.1.351, P.1, B.1.617.2, B.1.1.1 and B.1.621 variants were determined by pseudo-virus-based neutralization assay. Results showed that S2E7-mAb had cross-activity to S or S2 proteins of SARS-CoV-2, SARS-CoV and MERS-CoV, while with limited neutralizing activity to pseudo-viruses of SARS-CoV-2 WT and variants. It is undeniable that the binding and neutralizing activities of the S2-targeting mAbs are significantly weaker than the previously reported antibodies targeting RBD and NTD, but our study may provide some evidences for understanding immune protection and identifying targets for vaccine design based on the conserved S2 subunit.

RASCL: Rapid Assessment of Selection in CLades through molecular sequence analysis

An important unmet need revealed by the COVID-19 pandemic is the near-real-time identification of potentially fitness-altering mutations within rapidly growing SARS-CoV-2 lineages. Although powerful molecular sequence analysis methods are available to detect and characterize patterns of natural selection within modestly sized gene-sequence datasets, the computational complexity of these methods and their sensitivity to sequencing errors render them effectively inapplicable in large-scale genomic surveillance contexts. Motivated by the need to analyze new lineage evolution in near-real time using large numbers of genomes, we developed the Rapid Assessment of Selection within CLades (RASCL) pipeline. RASCL applies state of the art phylogenetic comparative methods to evaluate selective processes acting at individual codon sites and across whole genes. RASCL is scalable and produces automatically updated regular lineage-specific selection analysis reports: even for lineages that include tens or hundreds of thousands of sampled genome sequences. Key to this performance is (i) generation of automatically subsampled high quality datasets of gene/ORF sequences drawn from a selected "query" viral lineage; (ii) contextualization of these query sequences in codon alignments that include high-quality "background" sequences representative of global SARS-CoV-2 diversity; and (iii) the extensive parallelization of a suite of computationally intensive selection analysis tests. Within hours of being deployed to analyze a novel rapidly growing lineage of interest, RASCL will begin yielding JavaScript Object Notation (JSON)-formatted reports that can be either imported into third-party analysis software or explored in standard web-browsers using the premade RASCL interactive data visualization dashboard. By enabling the rapid detection of genome sites evolving under different selective regimes, RASCL is well-suited for near-real-time monitoring of the population-level selective processes that will likely underlie the emergence of future variants of concern in measurably evolving pathogens with extensive genomic surveillance.

Two Years of Genomic Surveillance in Belgium during the SARS-CoV-2 Pandemic to Attain Country-Wide Coverage and Monitor the Introduction and Spread of Emerging Variants

An adequate SARS-CoV-2 genomic surveillance strategy has proven to be essential for countries to obtain a thorough understanding of the variants and lineages being imported and successfully established within their borders. During 2020, genomic surveillance in Belgium was not structurally implemented but performed by individual research laboratories that had to acquire the necessary funds themselves to perform this important task. At the start of 2021, a nationwide genomic surveillance consortium was established in Belgium to markedly increase the country's genomic sequencing efforts (both in terms of intensity and representativeness), to perform quality control among participating laboratories, and to enable coordination and collaboration of research projects and publications. We here discuss the genomic surveillance efforts in Belgium before and after the establishment of its genomic sequencing consortium, provide an overview of the specifics of the consortium, and explore more details regarding the scientific studies that have been published as a result of the increased number of Belgian SARS-CoV-2 genomes that have become available.

Cellular interferon-gamma and interleukin-2 responses to SARS-CoV-2 structural proteins are broader and higher in those vaccinated after SARS-CoV-2 infection compared to vaccinees without prior SARS-CoV-2 infection

Class I- and Class II-restricted epitopes have been identified across the SARS-CoV-2 structural proteome. Vaccine-induced and post-infection SARS-CoV-2 T-cell responses are associated with COVID-19 recovery and protection, but the precise role of T-cell responses remains unclear, and how post-infection vaccination ('hybrid immunity') further augments this immunity To accomplish these goals, we studied healthy adult healthcare workers who were (a) uninfected and unvaccinated (n = 12), (b) uninfected and vaccinated with Pfizer-BioNTech BNT162b2 vaccine (2 doses n = 177, one dose n = 1) or Moderna mRNA-1273 vaccine (one dose, n = 1), and (c) previously infected with SARS-CoV-2 and vaccinated (BNT162b2, two doses, n = 6, one dose n = 1; mRNA-1273 two doses, n = 1). Infection status was determined by repeated PCR testing of participants. We used FluoroSpot Interferon-gamma (IFN-γ) and Interleukin-2 (IL-2) assays, using subpools of 15-mer peptides covering the S (10 subpools), N (4 subpools) and M (2 subpools) proteins. Responses were expressed as frequencies (percent positive responders) and magnitudes (spot forming cells/106 cytokine-producing peripheral blood mononuclear cells [PBMCs]). Almost all vaccinated participants with no prior infection exhibited IFN-γ, IL-2 and IFN-γ+IL2 responses to S glycoprotein subpools (89%, 93% and 27%, respectively) mainly directed to the S2 subunit and were more robust than responses to the N or M subpools. However, in previously infected and vaccinated participants IFN-γ, IL-2 and IFN-γ+IL2 responses to S subpools (100%, 100%, 88%) were substantially higher than vaccinated participants with no prior infection and were broader and directed against nine of the 10 S glycoprotein subpools spanning the S1 and S2 subunits, and all the N and M subpools. 50% of uninfected and unvaccinated individuals had IFN-γ but not IL2 or IFN-γ+IL2 responses against one S and one M subpools that were not increased after vaccination of uninfected or SARS-CoV-2-infected participants. Summed IFN-γ, IL-2, and IFN-γ+IL2 responses to S correlated with IgG responses to the S glycoprotein. These studies demonstrated that vaccinations with BNT162b2 or mRNA-1273 results in T cell-specific responses primarily against epitopes in the S2 subunit of the S glycoprotein, and that individuals that are vaccinated after SARS-CoV-2 infection develop broader and greater T cell responses to S1 and S2 subunits as well as the N and M proteins.

Continued Complexity of Mutations in Omicron Sublineages

The latest SARS-CoV-2 variant of concern (VOC), Omicron (B.1.1.529), has diversified into more than 300 sublineages. With an expanding number of newly emerging sublineages, the mutation profile is also becoming complicated. There exist mutually exclusive and revertant mutations in different sublineages. Omicron sublineages share some common mutations with previous VOCs (Alpha, Beta, Gamma, and Delta), indicating an evolutionary relationship between these VOCs. A diverse mutation profile at the spike–antibody interface, flexibility of the regions harboring mutations, mutation types, and coexisting mutations suggest that SARS-CoV-2’s evolution is far from over.

Omicron variant (B.1.1.529) and its sublineages: What do we know so far amid the emergence of recombinant variants of SARS-CoV-2?

Since the start of the COVID-19 pandemic, numerous variants of SARS-CoV-2 have been reported worldwide. The advent of variants of concern (VOCs) raises severe concerns amid the serious containment efforts against COVID-19 that include physical measures, pharmacological repurposing, immunization, and genomic/community surveillance. Omicron variant (B.1.1.529) has been identified as a highly modified, contagious, and crucial variant among the five VOCs of SARS-CoV-2. The increased affinity of the spike protein (S-protein), and host receptor, angiotensin converting enzyme-2 (ACE-2), due to a higher number of mutations in the receptor-binding domain (RBD) of the S-protein has been proposed as the primary reason for the decreased efficacy of majorly available vaccines against the Omicron variant and the increased transmissible nature of the Omicron variant. Because of its significant competitive advantage, the Omicron variant and its sublineages swiftly surpassed other variants to become the dominant circulating lineages in a number of nations. The Omicron variant has been identified as a prevalent strain in the United Kingdom and South Africa. Furthermore, the emergence of recombinant variants through the conjunction of the Omicron variant with other variants or by the mixing of the Omicron variant's sublineages/subvariants poses a major threat to humanity. This raises various issues and hazards regarding the Omicron variant and its sublineages, such as an Omicron variant breakout in susceptible populations among fully vaccinated persons. As a result, understanding the features and genetic implications of this variant is crucial. Hence, we explained in depth the evolution and features of the Omicron variant and analyzed the repercussions of spike mutations on infectiousness, dissemination ability, viral entry mechanism, and immune evasion. We also presented a viewpoint on feasible strategies for precluding and counteracting any future catastrophic emergence and spread of the omicron variant and its sublineages that could result in a detrimental wave of COVID-19 cases.

Epistatic Variations in the Omicron Receptor Binding Domain Can Enhance Host Recognition: An In Silico Assessment and Prediction

The hypermutated receptor binding domain (RBD) of the Omicron (B.1.1.529) lineage exhibits a different binding interface with human angiotensin-converting enzyme 2 (ACE2) relative to that of the wild-type Wuhan Hu-1, yet how the altered interaction will affect viral evolution is largely unknown. Here, we used molecular dynamics simulation to characterize the binding features of the Omicron BA.1/hACE2 complex and used free energy perturbation calculations to assess the ongoing and putative variations. The complex reveals a substantial rearrangement of the interfacial hydrogen-bond network: R493 of RBD forms a dynamic electrostatic interaction with both E35 and D38 of hACE2, which prohibits the hydrogen bonds of R498-D38 and Y449-D38. Whereas most circulating mutations minimally affect RBD binding to hACE2, the charge-altering mutation R493Q attenuates the affinity by abolishing the electrostatic interaction. However, the potential variants H505Y or N417K/R493Q could restore and gain even greater binding affinities than BA.1 as a result of their optimized interaction network and epistasis effects.

Molecular adaptations during viral epidemics

In 1977, the world witnessed both the eradication of smallpox and the beginning of the modern age of genomics. Over the following half-century, 7 epidemic viruses of international concern galvanized virologists across the globe and led to increasingly extensive virus genome sequencing. These sequencing efforts exerted over periods of rapid adaptation of viruses to new hosts, in particular, humans provide insight into the molecular mechanisms underpinning virus evolution. Investment in virus genome sequencing was dramatically increased by the unprecedented support for phylogenomic analyses during the COVID-19 pandemic. In this review, we attempt to piece together comprehensive molecular histories of the adaptation of variola virus, HIV-1 M, SARS, H1N1-SIV, MERS, Ebola, Zika, and SARS-CoV-2 to the human host. Disruption of genes involved in virus-host interaction in animal hosts, recombination including genome segment reassortment, and adaptive mutations leading to amino acid replacements in virus proteins involved in host receptor binding and membrane fusion are identified as the key factors in the evolution of epidemic viruses.