The emergence and ongoing convergent evolution of the SARS-CoV-2 N501Y lineages

Identifying Markers of Emerging SARS-CoV-2 Variants in Patients With Secondary Immunodeficiency

Since the end of 2019, the world has been challenged by the coronavirus disease 2019 (COVID-19) pandemic. With COVID-19 cases rising globally, severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) continues to evolve, resulting in the emergence of variants of interest (VOI) and of concern (VOC). Of the hundreds of millions infected, immunodeficient patients are one of the vulnerable cohorts that are most susceptible to this virus. These individuals include those with preexisting health conditions and/or those undergoing immunosuppressive treatment (secondary immunodeficiency). In these cases, several researchers have reported chronic infections in the presence of anti-COVID-19 treatments that may potentially lead to the evolution of the virus within the host. Such variations occurred in a variety of viral proteins, including key structural ones involved in pathogenesis such as spike proteins. Tracking and comparing such mutations with those arisen in the general population may provide information about functional sites within the SARS-CoV-2 genome. In this study, we reviewed the current literature regarding the specific features of SARS-CoV-2 evolution in immunocompromised patients and identified recurrent de novo amino acid changes in virus isolates of these patients that can potentially play an important role in SARS-CoV-2 pathogenesis and evolution.

Detection of SARS-CoV-2 intra-host recombination during superinfection with Alpha and Epsilon variants in New York City

Recombination is an evolutionary process by which many pathogens generate diversity and acquire novel functions. Although a common occurrence during coronavirus replication, detection of recombination is only feasible when genetically distinct viruses contemporaneously infect the same host. Here, we identify an instance of SARS-CoV-2 superinfection, whereby an individual was infected with two distinct viral variants: Alpha (B.1.1.7) and Epsilon (B.1.429). This superinfection was first noted when an Alpha genome sequence failed to exhibit the classic S gene target failure behavior used to track this variant. Full genome sequencing from four independent extracts reveals that Alpha variant alleles comprise around 75% of the genomes, whereas the Epsilon variant alleles comprise around 20% of the sample. Further investigation reveals the presence of numerous recombinant haplotypes spanning the genome, specifically in the spike, nucleocapsid, and ORF 8 coding regions. These findings support the potential for recombination to reshape SARS-CoV-2 genetic diversity.

Mutability Patterns Across the Spike Glycoprotein Reveal the Diverging and Lineage-specific Evolutionary Space of SARS-CoV-2

Mutations in the spike glycoprotein of SARS-CoV-2 allow the virus to probe the sequence space in search of higher-fitness states. New sublineages of SARS-CoV-2 variants-of-concern (VOCs) continuously emerge with such mutations. Interestingly, the sites of mutation in these sublineages vary between the VOCs. Whether such differences reflect the random nature of mutation appearance or distinct evolutionary spaces of spike in the VOCs is unclear. Here we show that each position of spike has a lineage-specific likelihood for mutations to appear and dominate descendent sublineages. This likelihood can be accurately estimated from the lineage-specific mutational profile of spike at a protein-wide level. The mutability environment of each position, including adjacent sites on the protein structure and neighboring sites on the network of comutability, accurately forecast changes in descendent sublineages. Mapping of imminent changes within the VOCs can contribute to the design of immunogens and therapeutics that address future forms of SARS-CoV-2.

SARS-CoV-2 Variants of Concern Hijack IFITM2 for Efficient Replication in Human Lung Cells

It has recently been shown that an early SARS-CoV-2 isolate (NL-02-2020) hijacks interferon-induced transmembrane proteins (IFITMs) for efficient replication in human lung cells, cardiomyocytes, and gut organoids. To date, several "variants of concern" (VOCs) showing increased infectivity and resistance to neutralization have emerged and globally replaced the early viral strains. Here, we determined whether the five current SARS-CoV-2 VOCs (Alpha, Beta, Gamma, Delta, and Omicron) maintained the dependency on IFITM proteins for efficient replication. We found that depletion of IFITM2 strongly reduces viral RNA production by all VOCs in the human epithelial lung cancer cell line Calu-3. Silencing of IFITM1 had modest effects, while knockdown of IFITM3 resulted in an intermediate phenotype. Strikingly, depletion of IFITM2 generally reduced infectious virus production by more than 4 orders of magnitude. In addition, an antibody directed against the N terminus of IFITM2 inhibited SARS-CoV-2 VOC replication in induced pluripotent stem cell (iPSC)-derived alveolar epithelial type II cells, thought to represent major viral target cells in the lung. In conclusion, endogenously expressed IFITM proteins (especially IFITM2) are critical cofactors for efficient replication of genuine SARS-CoV-2 VOCs, including the currently dominant Omicron variant. IMPORTANCE Recent data indicate that SARS-CoV-2 requires endogenously expressed IFITM proteins for efficient infection. However, the results were obtained with an early SARS-CoV-2 isolate. Thus, it remained to be determined whether IFITMs are also important cofactors for infection of emerging SARS-CoV-2 VOCs that outcompeted the original strains in the meantime. This includes the Omicron VOC, which currently dominates the pandemic. Here, we show that depletion of endogenous IFITM2 expression almost entirely prevents productive infection of Alpha, Beta, Gamma, Delta, and Omicron SARS-CoV-2 VOCs in human lung cells. In addition, an antibody targeting the N terminus of IFITM2 inhibited SARS-CoV-2 VOC replication in iPSC-derived alveolar epithelial type II cells. Our results show that SARS-CoV-2 VOCs, including the currently dominant Omicron variant, are strongly dependent on IFITM2 for efficient replication, suggesting a key proviral role of IFITMs in viral transmission and pathogenicity.

National Scale Real-Time Surveillance of SARS-CoV-2 Variants Dynamics by Wastewater Monitoring in Israel

In this report, we describe a national-scale monitoring of the SARS-CoV-2 (SC-2) variant dynamics in Israel, using multiple-time sampling of 13 wastewater treatment plants. We used a combination of inclusive and selective quantitative PCR assays that specifically identify variants A19/A20 or B.1.1.7 and tested each sample for the presence and relative viral RNA load of each variant. We show that between December 2020 and March 2021, a complete shift in the SC-2 variant circulation was observed, where the B.1.1.7 replaced the A19 in all examined test points. We further show that the normalized viral load (NVL) values and the average new cases per week reached a peak in January 2021 and then decreased gradually in almost all test points, in parallel with the progression of the national vaccination campaign, during February-March 2021. This study demonstrates the importance of monitoring SC-2 variant by using a combination of inclusive and selective PCR tests on a national scale through wastewater sampling, which is far more amendable for high-throughput monitoring compared with sequencing. This approach may be useful for real-time dynamics surveillance of current and future variants, such as the Omicron (BA.1, BA.2) and other variants.

Genetic Diversity and Evolutionary Convergence of Cryptic SARS-CoV-2 Lineages Detected Via Wastewater Sequencing

Wastewater-based epidemiology (WBE) is an effective way of tracking the appearance and spread of SARS-COV-2 lineages through communities. Beginning in early 2021, we implemented a targeted approach to amplify and sequence the receptor binding domain (RBD) of SARS-COV-2 to characterize viral lineages present in sewersheds. Over the course of 2021, we reproducibly detected multiple SARS-COV-2 RBD lineages that have never been observed in patient samples in 9 sewersheds located in 3 states in the USA. These cryptic lineages contained between 4 to 24 amino acid substitutions in the RBD and were observed intermittently in the sewersheds in which they were found for as long as 14 months. Many of the amino acid substitutions in these lineages occurred at residues also mutated in the Omicron variant of concern (VOC), often with the same substitution. One of the sewersheds contained a lineage that appeared to be derived from the Alpha VOC, but the majority of the lineages appeared to be derived from pre-VOC SARS-COV-2 lineages. Specifically, several of the cryptic lineages from New York City appeared to be derived from a common ancestor that most likely diverged in early 2020. While the source of these cryptic lineages has not been resolved, it seems increasingly likely that they were derived from immunocompromised patients or animal reservoirs. Our findings demonstrate that SARS-COV-2 genetic diversity is greater than what is commonly observed through routine SARS-CoV-2 surveillance. Wastewater sampling may more fully capture SARS-CoV-2 genetic diversity than patient sampling and could reveal new VOCs before they emerge in the wider human population.

Author Summary

During the COVID-19 pandemic, wastewater-based epidemiology has become an effective public health tool. Because many infected individuals shed SARS-CoV-2 in feces, wastewater has been monitored to reveal infection trends in the sewersheds from which the samples were derived. Here we report novel SARS-CoV-2 lineages in wastewater samples obtained from 3 different states in the USA. These lineages appeared in specific sewersheds intermittently over periods of up to 14 months, but generally have not been detected beyond the sewersheds in which they were initially found. Many of these lineages may have diverged in early 2020. Although these lineages share considerable overlap with each other, they have never been observed in patients anywhere in the world. While the wastewater lineages have similarities with lineages observed in long-term infections of immunocompromised patients, animal reservoirs cannot be ruled out as a potential source.

Increased Frequency of Indels in Hypervariable Regions of SARS-CoV-2 Proteins—A Possible Signature of Adaptive Selection

Most attention in the surveillance of evolving SARS-CoV-2 genome has been centered on nucleotide substitutions in the spike glycoprotein. We show that, as the pandemic extends into its second year, the numbers and ratio of genomes with in-frame insertions and deletions (indels) increases significantly, especially among the variants of concern (VOCs). Monitoring of the SARS-CoV-2 genome evolution shows that co-occurrence (i.e., highly correlated presence) of indels, especially deletions on spike N-terminal domain and non-structural protein 6 (NSP6) is a shared feature in several VOCs such as Alpha, Beta, Delta, and Omicron. Indels distribution is correlated with spike mutations associated with immune escape and growth in the number of genomes with indels coincides with the increasing population resistance due to vaccination and previous infections. Indels occur most frequently in the spike, but also in other proteins, especially those involved in interactions with the host immune system. We also showed that indels concentrate in regions of individual SARS-CoV-2 proteins known as hypervariable regions (HVRs) that are mostly located in specific loop regions. Structural analysis suggests that indels remodel viral proteins’ surfaces at common epitopes and interaction interfaces, affecting the virus’ interactions with host proteins. We hypothesize that the increased frequency of indels, the non-random distribution of them and their independent co-occurrence in several VOCs is another mechanism of response to elevated global population immunity.

The Evolution and Biology of SARS-CoV-2 Variants

Our understanding of the still unfolding severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) pandemic would have been extremely limited without the study of the genetics and evolution of this new human coronavirus. Large-scale genome-sequencing efforts have provided close to real-time tracking of the global spread and diversification of SARS-CoV-2 since its entry into the human population in late 2019. These data have underpinned analysis of its origins, epidemiology, and adaptations to the human population: principally immune evasion and increasing transmissibility. SARS-CoV-2, despite being a new human pathogen, was highly capable of human-to-human transmission. During its rapid spread in humans, SARS-CoV-2 has evolved independent new forms, the so-called "variants of concern," that are better optimized for human-to-human transmission. The most important adaptation of the bat coronavirus progenitor of both SARS-CoV-1 and SARS-CoV-2 for human infection (and other mammals) is the use of the angiotensin-converting enzyme 2 (ACE2) receptor. Relaxed structural constraints provide plasticity to SARS-related coronavirus spike protein permitting it to accommodate significant amino acid replacements of antigenic consequence without compromising the ability to bind to ACE2. Although the bulk of research has justifiably concentrated on the viral spike protein as the main determinant of antigenic evolution and changes in transmissibility, there is accumulating evidence for the contribution of other regions of the viral proteome to virus-host interaction. Whereas levels of community transmission of recombinants compromising genetically distinct variants are at present low, when divergent variants cocirculate, recombination between SARS-CoV-2 clades is being detected, increasing the risk that viruses with new properties emerge. Applying computational and machine learning methods to genome sequence data sets to generate experimentally verifiable predictions will serve as an early warning system for novel variant surveillance and will be important in future vaccine planning. Omicron, the latest SARS-CoV-2 variant of concern, has focused attention on step change antigenic events, "shift," as opposed to incremental "drift" changes in antigenicity. Both an increase in transmissibility and antigenic shift in Omicron led to it readily causing infections in the fully vaccinated and/or previously infected. Omicron's virulence, while reduced relative to the variant of concern it replaced, Delta, is very much premised on the past immune exposure of individuals with a clear signal that boosted vaccination protects from severe disease. Currently, SARS-CoV-2 has proven itself to be a dangerous new human respiratory pathogen with an unpredictable evolutionary capacity, leading to a risk of future variants too great not to ensure all regions of the world are screened by viral genome sequencing, protected through available and affordable vaccines, and have non-punitive strategies in place for detecting and responding to novel variants of concern.

Perspectives: SARS-CoV-2 Spike Convergent Evolution as a Guide to Explore Adaptive Advantage

Viruses rapidly co-evolve with their hosts. The 9 million sequenced SARS-CoV-2 genomes by March 2022 provide a detailed account of viral evolution, showing that all amino acids have been mutated many times. However, only a few became prominent in the viral population. Here, we investigated the emergence of the same mutations in unrelated parallel lineages and the extent of such convergent evolution on the molecular level in the spike (S) protein. We found that during the first phase of the pandemic (until mid 2021, before mass vaccination) 31 mutations evolved independently ≥3-times within separated lineages. These included all the key mutations in SARS-CoV-2 variants of concern (VOC) at that time, indicating their fundamental adaptive advantage. The omicron added many more mutations not frequently seen before, which can be attributed to the synergistic nature of these mutations, which is more difficult to evolve. The great majority (24/31) of S-protein mutations under convergent evolution tightly cluster in three functional domains; N-terminal domain, receptor-binding domain, and Furin cleavage site. Furthermore, among the S-protein receptor-binding motif mutations, ACE2 affinity-improving substitutions are favoured. Next, we determined the mutation space in the S protein that has been covered by SARS-CoV-2. We found that all amino acids that are reachable by single nucleotide changes have been probed multiple times in early 2021. The substitutions requiring two nucleotide changes have recently (late 2021) gained momentum and their numbers are increasing rapidly. These provide a large mutation landscape for SARS-CoV-2 future evolution, on which research should focus now.

Deep phylogenetic-based clustering analysis uncovers new and shared mutations in SARS-CoV-2 variants as a result of directional and convergent evolution

Nearly two decades after the last epidemic caused by a severe acute respiratory syndrome coronavirus (SARS-CoV), newly emerged SARS-CoV-2 quickly spread in 2020 and precipitated an ongoing global public health crisis. Both the continuous accumulation of point mutations, owed to the naturally imposed genomic plasticity of SARS-CoV-2 evolutionary processes, as well as viral spread over time, allow this RNA virus to gain new genetic identities, spawn novel variants and enhance its potential for immune evasion. Here, through an in-depth phylogenetic clustering analysis of upwards of 200,000 whole-genome sequences, we reveal the presence of previously unreported and hitherto unidentified mutations and recombination breakpoints in Variants of Concern (VOC) and Variants of Interest (VOI) from Brazil, India (Beta, Eta and Kappa) and the USA (Beta, Eta and Lambda). Additionally, we identify sites with shared mutations under directional evolution in the SARS-CoV-2 Spike-encoding protein of VOC and VOI, tracing a heretofore-undescribed correlation with viral spread in South America, India and the USA. Our evidence-based analysis provides well-supported evidence of similar pathways of evolution for such mutations in all SARS-CoV-2 variants and sub-lineages. This raises two pivotal points: (i) the co-circulation of variants and sub-lineages in close evolutionary environments, which sheds light onto their trajectories into convergent and directional evolution, and (ii) a linear perspective into the prospective vaccine efficacy against different SARS-CoV-2 strains.

In Silico Analysis Predicts a Limited Impact of SARS-CoV-2 Variants on CD8 T Cell Recognition

Since the start of the COVID-19 pandemic, mutations have led to the emergence of new SARS-CoV-2 variants, and some of these have become prominent or dominant variants of concern. This natural course of development can have an impact on how protective the previously naturally or vaccine induced immunity is. Therefore, it is crucial to understand whether and how variant specific mutations influence host immunity. To address this, we have investigated how mutations in the recent SARS-CoV-2 variants of interest and concern influence epitope sequence similarity, predicted binding affinity to HLA, and immunogenicity of previously reported SARS-CoV-2 CD8 T cell epitopes. Our data suggests that the vast majority of SARS-CoV-2 CD8 T cell recognized epitopes are not altered by variant specific mutations. Interestingly, for the CD8 T cell epitopes that are altered due to variant specific mutations, our analyses show there is a high degree of sequence similarity between mutated and reference SARS-CoV-2 CD8 T cell epitopes. However, mutated epitopes, primarily derived from the spike protein, in SARS-CoV-2 variants Delta, AY.4.2 and Mu display reduced predicted binding affinity to their restriction element. These findings indicate that the recent SARS-CoV-2 variants of interest and concern have limited ability to escape memory CD8 T cell responses raised by vaccination or prior infection with SARS-CoV-2 early in the pandemic. The overall low impact of the mutations on CD8 T cell cross-recognition is in accordance with the notion that mutations in SARS-CoV-2 are primarily the result of receptor binding affinity and antibody selection pressures exerted on the spike protein, unrelated to T cell immunity.

Quadrivalent mosaic HexaPro-bearing nanoparticle vaccine protects against infection of SARS-CoV-2 variants

Emerging SARS-CoV-2 variants of concern (VOCs) harboring multiple mutations in the spike protein raise concerns on effectiveness of current vaccines that rely on the ancestral spike protein. Here, we design a quadrivalent mosaic nanoparticle vaccine displaying spike proteins from the SARS-CoV-2 prototype and 3 different VOCs. The mosaic nanoparticle elicits equivalent or superior neutralizing antibodies against variant strains in mice and non-human primates with only small reduction in neutralization titers against the ancestral strain. Notably, it provides protection against infection with prototype and B.1.351 strains in mice. These results provide a proof of principle for the development of multivalent vaccines against pandemic and potential pre-emergent SARS-CoV-2 variants.

Emerging SARS-CoV-2 variants with multiple mutations raise concerns on vaccine effectiveness. Here, Kang et al. report that a quadrivalent mosaic nanoparticle vaccine displaying spike proteins from the SARS-CoV-2 prototype and three different VOCs confer protection against SARS-CoV-2 variants in mice.

The International Virus Bioinformatics Meeting 2022

The International Virus Bioinformatics Meeting 2022 took place online, on 23-25 March 2022, and has attracted about 380 participants from all over the world. The goal of the meeting was to provide a meaningful and interactive scientific environment to promote discussion and collaboration and to inspire and suggest new research directions and questions. The participants created a highly interactive scientific environment even without physical face-to-face interactions. This meeting is a focal point to gain an insight into the state-of-the-art of the virus bioinformatics research landscape and to interact with researchers in the forefront as well as aspiring young scientists. The meeting featured eight invited and 18 contributed talks in eight sessions on three days, as well as 52 posters, which were presented during three virtual poster sessions. The main topics were: SARS-CoV-2, viral emergence and surveillance, virus-host interactions, viral sequence analysis, virus identification and annotation, phages, and viral diversity. This report summarizes the main research findings and highlights presented at the meeting.

Distinct mutations and lineages of SARS-CoV-2 virus in the early phase of COVID-19 pandemic and subsequent 1-year global expansion

A novel coronavirus, SARS-CoV-2, has caused over 274 million cases and over 5.3 million deaths worldwide since it occurred in December 2019 in Wuhan, China. Here we conceptualized the temporospatial evolutionary and expansion dynamics of SARS-CoV-2 by taking a series of the cross-sectional view of viral genomes from early outbreak in January 2020 in Wuhan to the early phase of global ignition in early April, and finally to the subsequent global expansion by late December 2020. Based on the phylogenetic analysis of the early patients in Wuhan, Wuhan/WH04/2020 is supposed to be a more appropriate reference genome of SARS-CoV-2, instead of the first sequenced genome Wuhan-Hu-1. By scrutinizing the cases from the very early outbreak, we found a viral genotype from the Seafood Market in Wuhan featured with two concurrent mutations (i.e., M type) had become the overwhelmingly dominant genotype (95.3%) of the pandemic 1 year later. By analyzing 4013 SARS-CoV-2 genomes from different continents by early April, we were able to interrogate the viral genomic composition dynamics of the initial phase of global ignition over a time span of 14 weeks. Eleven major viral genotypes with unique geographic distributions were also identified. WE1 type, a descendant of M and predominantly witnessed in western Europe, consisted of half of all the cases (50.2%) at the time. The mutations of major genotypes at the same hierarchical level were mutually exclusive, which implies that various genotypes bearing the specific mutations were propagated during human-to-human transmission, not by accumulating hot-spot mutations during the replication of individual viral genomes. As the pandemic was unfolding, we also used the same approach to analyze 261 323 SARS-CoV-2 genomes from the world since the outbreak in Wuhan (i.e., including all the publicly available viral genomes) to recapitulate our findings over 1-year time span. By December 25, 2020, 95.3% of global cases were M type and 93.0% of M-type cases were WE1. In fact, at present all the five variants of concern (VOC) are the descendants of WE1 type. This study demonstrates that viral genotypes can be utilized as molecular barcodes in combination with epidemiologic data to monitor the spreading routes of the pandemic and evaluate the effectiveness of control measures. Moreover, the dynamics of viral mutational spectrum in the study may help the early identification of new strains in patients to reduce further spread of infection, guide the development of molecular diagnosis and vaccines against COVID-19, and help assess their accuracy and efficacy in real world at real time.

Mapping Data to Deep Understanding: Making the Most of the Deluge of SARS-CoV-2 Genome Sequences

Next-generation sequencing has been essential to the global response to the COVID-19 pandemic. As of January 2022, nearly 7 million severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) sequences are available to researchers in public databases. Sequence databases are an abundant resource from which to extract biologically relevant and clinically actionable information. As the pandemic has gone on, SARS-CoV-2 has rapidly evolved, involving complex genomic changes that challenge current approaches to classifying SARS-CoV-2 variants. Deep sequence learning could be a potentially powerful way to build complex sequence-to-phenotype models. Unfortunately, while they can be predictive, deep learning typically produces "black box" models that cannot directly provide biological and clinical insight. Researchers should therefore consider implementing emerging methods for visualizing and interpreting deep sequence models. Finally, researchers should address important data limitations, including (i) global sequencing disparities, (ii) insufficient sequence metadata, and (iii) screening artifacts due to poor sequence quality control.

SARS-CoV-2: An Overview of the Genetic Profile and Vaccine Effectiveness of the Five Variants of Concern

With the onset of the COVID-19 pandemic, enormous efforts have been made to understand the genus SARS-CoV-2. Due to the high rate of global transmission, mutations in the viral genome were inevitable. A full understanding of the viral genome and its possible changes represents one of the crucial aspects of pandemic management. Structural protein S plays an important role in the pathogenicity of SARS-CoV-2, mutations occurring at this level leading to viral forms with increased affinity for ACE2 receptors, higher transmissibility and infectivity, resistance to neutralizing antibodies and immune escape, increasing the risk of infection and disease severity. Thus, five variants of concern are currently being discussed, Alpha, Beta, Gamma, Delta and Omicron. In the present review, a comprehensive summary of the following critical aspects regarding SARS-CoV-2 has been made: (i) the genomic characteristics of SARS-CoV-2; (ii) the pathological mechanism of transmission, penetration into the cell and action on specific receptors; (iii) mutations in the SARS-CoV-2 genome; and (iv) possible implications of mutations in diagnosis, treatment, and vaccination.

Rapid and parallel adaptive mutations in spike S1 drive clade success in SARS-CoV-2

The SARS-CoV-2 pandemic has resulted in numerous virus variants, some of which have altered receptor-binding or antigenic phenotypes. Here, we quantify the degree to which adaptive evolution is driving the accumulation of mutations across the genome. We correlate clade growth with mutation accumulation, compare rates of nonsynonymous to synonymous divergence, assess temporal clustering of mutations, and evaluate the evolutionary success of individual mutations. We find that spike S1 is the focus of adaptive evolution but also identify positively selected mutations in other proteins (notably Nsp6) that are sculpting the evolutionary trajectory of SARS-CoV-2. Adaptive changes in S1 accumulated rapidly, resulting in a remarkably high ratio of nonsynonymous to synonymous divergence that is 2.5× greater than that observed in influenza hemagglutinin HA1 at the beginning of the 2009 H1N1 pandemic. These findings uncover a high degree of adaptation in S1 and suggest that SARS-CoV-2 might undergo antigenic drift.

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

Among the 30 nonsynonymous nucleotide substitutions in the Omicron S-gene are 13 that have only rarely been seen in other SARS-CoV-2 sequences. These mutations cluster within three functionally important regions of the S-gene at sites that will likely impact (1) interactions between subunits of the Spike trimer and the predisposition of subunits to shift from down to up configurations, (2) interactions of Spike with ACE2 receptors, and (3) the priming of Spike for membrane fusion. We show here that, based on both the rarity of these 13 mutations in intrapatient sequencing reads and patterns of selection at the codon sites where the mutations occur in SARS-CoV-2 and related sarbecoviruses, prior to the emergence of Omicron the mutations would have been predicted to decrease the fitness of any virus within which they occurred. We further propose that the mutations in each of the three clusters therefore cooperatively interact to both mitigate their individual fitness costs, and, in combination with other mutations, adaptively alter the function of Spike. Given the evident epidemic growth advantages of Omicron overall previously known SARS-CoV-2 lineages, it is crucial to determine both how such complex and highly adaptive mutation constellations were assembled within the Omicron S-gene, and why, despite unprecedented global genomic surveillance efforts, the early stages of this assembly process went completely undetected.

Advances in Pathogenesis, Progression, Potential Targets and Targeted Therapeutic Strategies in SARS-CoV-2-Induced COVID-19

As the new year of 2020 approaches, an acute respiratory disease quietly caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), also known as coronavirus disease 2019 (COVID-19) was reported in Wuhan, China. Subsequently, COVID-19 broke out on a global scale and formed a global public health emergency. To date, the destruction that has lasted for more than two years has not stopped and has caused the virus to continuously evolve new mutant strains. SARS-CoV-2 infection has been shown to cause multiple complications and lead to severe disability and death, which has dealt a heavy blow to global development, not only in the medical field but also in social security, economic development, global cooperation and communication. To date, studies on the epidemiology, pathogenic mechanism and pathological characteristics of SARS-CoV-2-induced COVID-19, as well as target confirmation, drug screening, and clinical intervention have achieved remarkable effects. With the continuous efforts of the WHO, governments of various countries, and scientific research and medical personnel, the public's awareness of COVID-19 is gradually deepening, a variety of prevention methods and detection methods have been implemented, and multiple vaccines and drugs have been developed and urgently marketed. However, these do not appear to have completely stopped the pandemic and ravages of this virus. Meanwhile, research on SARS-CoV-2-induced COVID-19 has also seen some twists and controversies, such as potential drugs and the role of vaccines. In view of the fact that research on SARS-CoV-2 and COVID-19 has been extensive and in depth, this review will systematically update the current understanding of the epidemiology, transmission mechanism, pathological features, potential targets, promising drugs and ongoing clinical trials, which will provide important references and new directions for SARS-CoV-2 and COVID-19 research.

Clinical Characteristics, Transmissibility, Pathogenicity, Susceptible Populations, and Re-infectivity of Prominent COVID-19 Variants

In addition to the rapid, global spread of SARS-CoV-2, new and comparatively more contagious variants are of considerable concern. These emerging mutations have become a threat to the global public health, creating COVID-19 surges in different countries. However, information on these emerging variants is limited and scattered. In this review, we discuss new variants that have emerged worldwide and identify several variants of concern, such as B.1.1.7, B.1.351, P.1, B.1.617.2 and B.1.1.529, and their basic characteristics. Other significant variants such as C.37, B.1.621, B.1.525, B.1.526, AZ.5, C.1.2, and B.1.617.1 are also discussed. This review highlights the clinical characteristics of these variants, including transmissibility, pathogenicity, susceptible population, and re-infectivity. It provides the latest information on the recent variants of SARS-CoV-2. The summary of this information will help researchers formulate reasonable strategies to curb the ongoing COVID-19 pandemic.

Comparative Analysis of SARS-CoV-2 Variants of Concern, Including Omicron, Highlights Their Common and Distinctive Amino Acid Substitution Patterns, Especially at the Spike ORF

In order to gain a deeper understanding of the recently emerged and highly divergent Omicron variant of concern (VoC), a study of amino acid substitution (AAS) patterns was performed and compared with those of the other four successful variants of concern (Alpha, Beta, Gamma, Delta) and one closely related variant of interest (VoI-Lambda). The Spike ORF consistently emerges as an AAS hotspot in all six lineages, but in Omicron this enrichment is significantly higher. The progenitors of each of these VoC/VoI lineages underwent positive selection in the Spike ORF. However, once they were established, their Spike ORFs have been undergoing purifying selection, despite the application of global vaccination schemes from 2021 onwards. Our analyses reject the hypothesis that the heavily mutated receptor binding domain (RBD) of the Omicron Spike was introduced via recombination from another closely related Sarbecovirus. Thus, successive point mutations appear as the most parsimonious scenario. Intriguingly, in each of the six lineages, we observed a significant number of AAS wherein the new residue is not present at any homologous site among the other known Sarbecoviruses. Such AAS should be further investigated as potential adaptations to the human host. By studying the phylogenetic distribution of AAS shared between the six lineages, we observed that the Omicron (BA.1) lineage had the highest number (8/10) of recurrent mutations.

Tracking SARS-CoV-2 mutations and variants through the COG-UK-Mutation Explorer

COG-UK Mutation Explorer (COG-UK-ME, https://sars2.cvr.gla.ac.uk/cog-uk/-last accessed date 16 March 2022) is a web resource that displays knowledge and analyses on SARS-CoV-2 virus genome mutations and variants circulating in the UK, with a focus on the observed amino acid replacements that have an antigenic role in the context of the human humoral and cellular immune response. This analysis is based on more than 2 million genome sequences (as of March 2022) for UK SARS-CoV-2 data held in the CLIMB-COVID centralised data environment. COG-UK-ME curates these data and displays analyses that are cross-referenced to experimental data collated from the primary literature. The aim is to track mutations of immunological importance that are accumulating in current variants of concern and variants of interest that could alter the neutralising activity of monoclonal antibodies (mAbs), convalescent sera, and vaccines. Changes in epitopes recognised by T cells, including those where reduced T cell binding has been demonstrated, are reported. Mutations that have been shown to confer SARS-CoV-2 resistance to antiviral drugs are also included. Using visualisation tools, COG-UK-ME also allows users to identify the emergence of variants carrying mutations that could decrease the neutralising activity of both mAbs present in therapeutic cocktails, e.g. Ronapreve. COG-UK-ME tracks changes in the frequency of combinations of mutations and brings together the curated literature on the impact of those mutations on various functional aspects of the virus and therapeutics. Given the unpredictable nature of SARS-CoV-2 as exemplified by yet another variant of concern, Omicron, continued surveillance of SARS-CoV-2 remains imperative to monitor virus evolution linked to the efficacy of therapeutics.

Rapid epidemic expansion of the SARS-CoV-2 Omicron variant in southern Africa

The SARS-CoV-2 epidemic in southern Africa has been characterized by three distinct waves. The first was associated with a mix of SARS-CoV-2 lineages, while the second and third waves were driven by the Beta (B.1.351) and Delta (B.1.617.2) variants, respectively1-3. In November 2021, genomic surveillance teams in South Africa and Botswana detected a new SARS-CoV-2 variant associated with a rapid resurgence of infections in Gauteng province, South Africa. Within three days of the first genome being uploaded, it was designated a variant of concern (Omicron, B.1.1.529) by the World Health Organization and, within three weeks, had been identified in 87 countries. The Omicron variant is exceptional for carrying over 30 mutations in the spike glycoprotein, which are predicted to influence antibody neutralization and spike function4. Here we describe the genomic profile and early transmission dynamics of Omicron, highlighting the rapid spread in regions with high levels of population immunity.

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

SARS-CoV-2 evolution threatens vaccine- and natural infection-derived immunity as well as the efficacy of therapeutic antibodies. To improve public health preparedness, we sought to predict which existing amino acid mutations in SARS-CoV-2 might contribute to future variants of concern. We tested the predictive value of features comprising epidemiology, evolution, immunology, and neural network-based protein sequence modeling, and identified primary biological drivers of SARS-CoV-2 intra-pandemic evolution. We found evidence that ACE2-mediated transmissibility and resistance to population-level host immunity has waxed and waned as a primary driver of SARS-CoV-2 evolution over time. We retroactively identified with high accuracy (area under the receiver operator characteristic curve, AUROC=0.92-0.97) mutations that will spread, at up to four months in advance, across different phases of the pandemic. The behavior of the model was consistent with a plausible causal structure wherein epidemiological covariates combine the effects of diverse and shifting drivers of viral fitness. We applied our model to forecast mutations that will spread in the future and characterize how these mutations affect the binding of therapeutic antibodies. These findings demonstrate that it is possible to forecast the driver mutations that could appear in emerging SARS-CoV-2 variants of concern. We validate this result against Omicron, showing elevated predictive scores for its component mutations prior to emergence, and rapid score increase across daily forecasts during emergence. This modeling approach may be applied to any rapidly evolving pathogens with sufficiently dense genomic surveillance data, such as influenza, and unknown future pandemic viruses.

Misidentification of the SARS-CoV-2 Mu variant using commercial mutation screening assays

Detection of mutations by multiplex real-time RT-PCR is a widely used method for the screening of SARS-CoV-2 variants, but this method has several limitations. We describe three cases in which a Mu strain containing the mutation K417N was initially misclassified as the Beta variant. We recommend the detection of P681H to distinguish between these two variants. Our experience highlights the importance of keeping track of new variants and mutations in order to adapt the current workflows.

Repurposing Antifungals for Host-Directed Antiviral Therapy?

Because of their epidemic and pandemic potential, emerging viruses are a major threat to global healthcare systems. While vaccination is in general a straightforward approach to prevent viral infections, immunization can also cause escape mutants that hide from immune cell and antibody detection. Thus, other approaches than immunization are critical for the management and control of viral infections. Viruses are prone to mutations leading to the rapid emergence of resistant strains upon treatment with direct antivirals. In contrast to the direct interference with pathogen components, host-directed therapies aim to target host factors that are essential for the pathogenic replication cycle or to improve the host defense mechanisms, thus circumventing resistance. These relatively new approaches are often based on the repurposing of drugs which are already licensed for the treatment of other unrelated diseases. Here, we summarize what is known about the mechanisms and modes of action for a potential use of antifungals as repurposed host-directed anti-infectives for the therapeutic intervention to control viral infections.

The Emergence of SARS-CoV-2 Variants of Concern Is Driven by Acceleration of the Substitution Rate

The ongoing SARS-CoV-2 pandemic has seen an unprecedented amount of rapidly generated genome data. These data have revealed the emergence of lineages with mutations associated to transmissibility and antigenicity, known as variants of concern (VOCs). A striking aspect of VOCs is that many of them involve an unusually large number of defining mutations. Current phylogenetic estimates of the substitution rate of SARS-CoV-2 suggest that its genome accrues around two mutations per month. However, VOCs can have 15 or more defining mutations and it is hypothesized that they emerged over the course of a few months, implying that they must have evolved faster for a period of time. We analyzed genome sequence data from the GISAID database to assess whether the emergence of VOCs can be attributed to changes in the substitution rate of the virus and whether this pattern can be detected at a phylogenetic level using genome data. We fit a range of molecular clock models and assessed their statistical performance. Our analyses indicate that the emergence of VOCs is driven by an episodic increase in the substitution rate of around 4-fold the background phylogenetic rate estimate that may have lasted several weeks or months. These results underscore the importance of monitoring the molecular evolution of the virus as a means of understanding the circumstances under which VOCs may emerge.

Limited spread of a rare spike E484K-harboring SARS-CoV-2 in Marseille, France

We detected SARS-CoV-2 of PANGO lineage R.1 with the spike substitution E484K in three patients. Eleven other sequences in France and 8,831 worldwide were available from GISAID, 92% originating from Japan. The three genome sequences from our institute were phylogenetically closest to another from Guinea-Conakry, where one of the patients had travelled. These viruses did not exhibit any unusual features in cell culture. Spike structural predictions indicated a 1.3-time higher transmissibility index than for the globally spread B.1.1.7 variant but also an affinity loss for gangliosides that might have slowed dissemination. The spread of new SARS-CoV-2 mutants/variants is still not well understood and therefore difficult to predict, and this hinders implementation of effective preventive measures, including adapted vaccines.

Occurrence of a substitution or deletion of SARS-CoV-2 spike amino acid 677 in various lineages in Marseille, France

Great concerns have been raised about SARS-CoV-2 variants over the past six months. At the end of 2020, an increasing incidence of spike substitutions Q677H/P was described in the USA, which involved six independent lineages. We searched for changes to this amino acid in the sequence database of SARS-CoV-2 genomes obtained at the IHU Méditerranée Infection (Marseille, France) from 3634 patients sampled between February 2020 and April 2021. In seven genomes (0.2%), we found a deletion of five amino acids at spike positions 675-679 (QTQTN) including Q677, and in 76 genomes (2.3%) we found a Q677H substitution. The 83 genomes were classified in ten different Pangolin lineages. Genomes with a spike Q677 deletion were obtained from respiratory samples collected in six cases between 28 March 2020 and 12 October 2020 and in one case on 1 February 2021. The Q677H substitution was found in genomes all obtained from respiratory samples collected from 19 January 2021 and were classified in seven different lineages. Most of these genomes (41 cases) were of UK variant. Two others were classified in the B.1.160 Pangolin lineage (Marseille-4 variant) which was first detected in July 2020 in our institute but was devoid of this substitution until 19 January 2021. Also, eight genomes were classified in the A.27/Marseille-501 lineage which was first detected in our institute in January 2021 and which either harboured or did not harbour the Q677H substitution. Thus, the spike Q677H substitution should be considered as another example of convergent evolution, as it is the case of spike substitutions L18F, E484K, L452R, and N501Y which also independently appeared in various lineages.

Increasing the frequency of omicron variant mutations boosts the immune response and may reduce the virus virulence

The coronavirus has posed a serious threat to the world since its discovery in Wuhan in 2019. Beta, gamma, delta, and the final omicron variants have emerged as a result of several mutations in the virion structure. The Australian Omicron S protein variant contains 37 mutations out of a total of 67 mutations. According to preliminary data from South Africa, Omicron variant infection is not associated with any particular symptoms. The purpose of this research was to determine how changes in the structure of the S protein alter the protein's interaction with the ACE2 receptor. The Omicron variant stimulates the immune response more than the wild strain.

Rapid and parallel adaptive mutations in spike S1 drive clade success in SARS-CoV-2

Given the importance of variant SARS-CoV-2 viruses with altered receptor-binding or antigenic phenotypes, we sought to quantify the degree to which adaptive evolution is driving accumulation of mutations in the SARS-CoV-2 genome. Here we assessed adaptive evolution across genes in the SARS-CoV-2 genome by correlating clade growth with mutation accumulation as well as by comparing rates of nonsynonymous to synonymous divergence, clustering of mutations across the SARS-CoV-2 phylogeny and degree of convergent evolution of individual mutations. We find that spike S1 is the focus of adaptive evolution, but also identify positively-selected mutations in other genes that are sculpting the evolutionary trajectory of SARS-CoV-2. Adaptive changes in S1 accumulated rapidly, resulting in a remarkably high ratio of nonsynonymous to synonymous divergence that is 2.5X greater than that observed in HA1 at the beginning of the 2009 H1N1 pandemic.

RASCL: Rapid Assessment Of SARS-CoV-2 Clades Through Molecular Sequence Analysis

An important component of efforts to manage the ongoing COVID19 pandemic is the R apid A ssessment of how natural selection contributes to the emergence and proliferation of potentially dangerous S ARS-CoV-2 lineages and CL ades (RASCL). The RASCL pipeline enables continuous comparative phylogenetics-based selection analyses of rapidly growing clade-focused genome surveillance datasets, such as those produced following the initial detection of potentially dangerous variants. From such datasets RASCL automatically generates down-sampled codon alignments of individual genes/ORFs containing contextualizing background reference sequences, analyzes these with a battery of selection tests, and outputs results as both machine readable JSON files, and interactive notebook-based visualizations.

AVAILABILITY

RASCL is available from a dedicated repository at https://github.com/veg/RASCL and as a Galaxy workflow https://usegalaxy.eu/u/hyphy/w/rascl . Existing clade/variant analysis results are available here: https://observablehq.com/@aglucaci/rascl .

CONTACT

Dr. Sergei L Kosakovsky Pond ( spond@temple.edu ).

SUPPLEMENTARY INFORMATION

N/A.

Selection analysis identifies unusual clustered mutational changes in Omicron lineage BA.1 that likely impact Spike function

Among the 30 non-synonymous nucleotide substitutions in the Omicron S-gene are 13 that have only rarely been seen in other SARS-CoV-2 sequences. These mutations cluster within three functionally important regions of the S-gene at sites that will likely impact (i) interactions between subunits of the Spike trimer and the predisposition of subunits to shift from down to up configurations, (ii) interactions of Spike with ACE2 receptors, and (iii) the priming of Spike for membrane fusion. We show here that, based on both the rarity of these 13 mutations in intrapatient sequencing reads and patterns of selection at the codon sites where the mutations occur in SARS-CoV-2 and related sarbecoviruses, prior to the emergence of Omicron the mutations would have been predicted to decrease the fitness of any genomes within which they occurred. We further propose that the mutations in each of the three clusters therefore cooperatively interact to both mitigate their individual fitness costs, and adaptively alter the function of Spike. Given the evident epidemic growth advantages of Omicron over all previously known SARS-CoV-2 lineages, it is crucial to determine both how such complex and highly adaptive mutation constellations were assembled within the Omicron S-gene, and why, despite unprecedented global genomic surveillance efforts, the early stages of this assembly process went completely undetected.

Heterologous infection and vaccination shapes immunity against SARS-CoV-2 variants

The impact of the initial severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infecting strain on downstream immunity to heterologous variants of concern (VOCs) is unknown. Studying a longitudinal healthcare worker cohort, we found that after three antigen exposures (infection plus two vaccine doses), S1 antibody, memory B cells, and heterologous neutralization of B.1.351, P.1, and B.1.617.2 plateaued, whereas B.1.1.7 neutralization and spike T cell responses increased. Serology using the Wuhan Hu-1 spike receptor binding domain poorly predicted neutralizing immunity against VOCs. Neutralization potency against VOCs changed with heterologous virus encounter and number of antigen exposures. Neutralization potency fell differentially depending on targeted VOCs over the 5 months from the second vaccine dose. Heterologous combinations of spike encountered during infection and vaccination shape subsequent cross-protection against VOC, with implications for future-proof next-generation vaccines.

Detection of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) and its first variants in fourplex real-time quantitative reverse transcription-PCR assays

The early diagnosis of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infections is required to identify and isolate contagious patients to prevent further transmission of SARS-CoV-2. In this study, we present a multitarget real-time TaqMan reverse transcription PCR (rRT-PCR) assay for the quantitative detection of SARS-CoV-2 and some of its circulating variants harboring mutations that give the virus a selective advantage. Seven different primer-probe sets that included probes containing locked nucleic acid (LNA) nucleotides were designed to amplify specific wild-type and mutant sequences in Orf1ab, Envelope (E), Spike (S), and Nucleocapsid (N) genes. Furthermore, a newly developed primer-probe set targeted human β₂-microglobulin (B2M) as a highly sensitive internal control for RT efficacy. All singleplex and fourplex assays detected ≤ 14 copies/reaction of quantified synthetic RNA transcripts, with a linear amplification range of nine logarithmic orders. Primer-probe sets for detection of SARS-CoV-2 exhibited no false-positive amplifications with other common respiratory pathogens, including human coronaviruses NL63, 229E, OC43, and HKU-1. Fourplex assays were evaluated using 160 clinical samples positive for SARS-CoV-2. Results showed that SARS-CoV-2 viral RNA was detected in all samples, including viral strains harboring mutations in the Spike coding sequence that became dominant in the pandemic. Given the emergence of SARS-CoV-2 variants and their rapid spread in some populations, fourplex rRT-PCR assay containing four primer-probe sets represents a reliable approach to allow quicker detection of circulating relevant variants in a single reaction.

An entropy-based study on mutational trajectory of SARS-CoV-2 in India

The pandemic of COVID-19 has been haunting us for almost the past two years. Although, the vaccination drive is in full swing throughout the world, different mutations of the SARS-CoV-2 virus are making it very difficult to put an end to the pandemic. The second wave in India, one of the worst sufferers of this pandemic, can be mainly attributed to the Delta variant i.e. B.1.617.2. Thus, it is very important to analyse and understand the mutational trajectory of SARS-CoV-2 through the study of the 26 virus proteins. In this regard, more than 17,000 protein sequences of Indian SARS-CoV-2 genomes are analysed using entropy-based approach in order to find the monthly mutational trajectory. Furthermore, Hellinger distance is also used to show the difference of the mutation events between the consecutive months for each of the 26 SARS-CoV-2 protein. The results show that the mutation rates and the mutation events of the viral proteins though changing in the initial months, start stabilizing later on for mainly the four structural proteins while the non-structural proteins mostly exhibit a more constant trend. As a consequence, it can be inferred that the evolution of the new mutative configurations will eventually reduce.

The atomic portrait of SARS-CoV-2 as captured by cryo-electron microscopy

Transmission electron microscopy has historically been indispensable for virology research, as it offers unique insight into virus function. In the past decade, as cryo-electron microscopy (cryo-EM) has matured and become more accessible, we have been able to peer into the structure of viruses at the atomic level and understand how they interact with the host cell, with drugs or with antibodies. Perhaps, there was no time in recent history where cryo-EM was more needed, as SARS-CoV-2 has spread around the globe, causing millions of deaths and almost unquantifiable economic devastation. In this concise review, we aim to mark the most important contributions of cryo-EM to understanding the structure and function of SARS-CoV-2 proteins, from surface spikes to the virus core and from virus-receptor interactions to antibody binding.

Intra-Host SARS-CoV-2 Evolution in the Gut of Mucosally-Infected Chlorocebus aethiops (African Green Monkeys)

In recent months, several SARS-CoV-2 variants have emerged that enhance transmissibility and escape host humoral immunity. Hence, the tracking of viral evolutionary trajectories is clearly of great importance. Little is known about SARS-CoV-2 evolution in nonhuman primate models used to test vaccines and therapies and to model human disease. Viral RNA was sequenced from rectal swabs from Chlorocebus aethiops (African green monkeys) after experimental respiratory SARS-CoV-2 infection. Two distinct patterns of viral evolution were identified that were shared between all collected samples. First, mutations in the furin cleavage site that were initially present in the virus as a consequence of VeroE6 cell culture adaptation were not detected in viral RNA recovered in rectal swabs, confirming the necessity of this motif for viral infection in vivo. Three amino acid changes were also identified; ORF 1a S2103F, and spike D215G and H655Y, which were detected in rectal swabs from all sampled animals. These findings are demonstrative of intra-host SARS-CoV-2 evolution and may identify a host-adapted variant of SARS-CoV-2 that would be useful in future primate models involving SARS-CoV-2 infection.

[The different phases of molecular and antigenic evolution of SARS-CoV-2 viruses during the 20 months following its emergence]

From its emergence in December 2019 and until the end of the fourth pandemic wave in October 2021, SARS-CoV-2 circulation has been associated with significant molecular evolutions of the virus. These were linked to mutations that have led to new virus linages with replication advantages as a result of increased transmission, or partial immune escape in the context of progressively increasing global immunisation. The pandemic context with large scale epidemics massive outbreaks observed in highly populated areas has favoured this emergence of "variants". During the 20 months period, at least three evolutionary phases have been observed, leading to the situation observed in October 2021. For the first time, an unprecedented worldwide surveillance effort has been conducted to monitor the circulation of the emerging virus, with rapid data sharing. This molecular surveillance system has provided an accurate description of the circulating viruses, and their evolution. The implementation of these tools and skills able to provide SARS-CoV-2 molecular epidemiological data has upgraded the global capacity for surveillance worldwide, and may allow us to be better prepared for a future pandemic episode.

OUP accepted manuscript

The Zika virus (ZIKV) disease caused a public health emergency of international concern that started in February 2016. The overall number of ZIKV-related cases increased until November 2016, after which it declined sharply. While the evaluation of the potential risk and impact of future arbovirus epidemics remains challenging, intensified surveillance efforts along with a scale-up of ZIKV whole-genome sequencing provide an opportunity to understand the patterns of genetic diversity, evolution, and spread of ZIKV. However, a classification system that reflects the true extent of ZIKV genetic variation is lacking. Our objective was to characterize ZIKV genetic diversity and phylodynamics, identify genomic footprints of differentiation patterns, and propose a dynamic classification system that reflects its divergence levels. We analysed a curated dataset of 762 publicly available sequences spanning the full-length coding region of ZIKV from across its geographical span and collected between 1947 and 2021. The definition of genetic groups was based on comprehensive evolutionary dynamics analyses, which included recombination and phylogenetic analyses, within- and between-group pairwise genetic distances comparison, detection of selective pressure, and clustering analyses. Evidence for potential recombination events was detected in a few sequences. However, we argue that these events are likely due to sequencing errors as proposed in previous studies. There was evidence of strong purifying selection, widespread across the genome, as also detected for other arboviruses. A total of 50 sites showed evidence of positive selection, and for a few of these sites, there was amino acid (AA) differentiation between genetic clusters. Two main genetic clusters were defined, ZA and ZB, which correspond to the already characterized ‘African’ and ‘Asian’ genotypes, respectively. Within ZB, two subgroups, ZB.1 and ZB.2, represent the Asiatic and the American (and Oceania) lineages, respectively. ZB.1 is further subdivided into ZB.1.0 (a basal Malaysia sequence sampled in the 1960s and a recent Indian sequence), ZB.1.1 (South-Eastern Asia, Southern Asia, and Micronesia sequences), and ZB.1.2 (very similar sequences from the outbreak in Singapore). ZB.2 is subdivided into ZB.2.0 (basal American sequences and the sequences from French Polynesia, the putative origin of South America introduction), ZB.2.1 (Central America), and ZB.2.2 (Caribbean and North America). This classification system does not use geographical references and is flexible to accommodate potential future lineages. It will be a helpful tool for studies that involve analyses of ZIKV genomic variation and its association with pathogenicity and serve as a starting point for the public health surveillance and response to on-going and future epidemics and to outbreaks that lead to the emergence of new variants.

A Highly Conserved Peptide Vaccine Candidate Activates Both Humoral and Cellular Immunity Against SARS-CoV-2 Variant Strains

Facing the imminent need for vaccine candidates with cross-protection against globally circulating severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) mutants, we present a conserved antigenic peptide RBD9.1 with both T-cell and B-cell epitopes. RBD9.1 can be recognized by coronavirus disease 2019 (COVID-19) convalescent serum, particularly for those with high neutralizing potency. Immunization with RBD9.1 can successfully induce the production of the receptor-binding domain (RBD)-specific antibodies in Balb/c mice. Importantly, the immunized sera exhibit sustained neutralizing efficacy against multiple dominant SARS-CoV-2 variant strains, including B.1.617.2 that carries a point mutation (S^L452R) within the sequence of RBD9.1. Specifically, S^Y451 and S^Y454 are identified as the key amino acids for the binding of the induced RBD-specific antibodies to RBD9.1. Furthermore, we have confirmed that the RBD9.1 antigenic peptide can induce a S^448-456 (NYNYLYRLF)-specific CD8⁺ T-cell response. Both RBD9.1-specific B cells and the S^448-456-specific T cells can still be activated more than 3 months post the last immunization. This study provides a potential vaccine candidate that can generate long-term protective efficacy over SARS-CoV-2 variants, with the unique functional mechanism of activating both humoral and cellular immunity.

Omicron-B.1.1.529 leads to widespread escape from neutralizing antibody responses

On the 24 th November 2021 the sequence of a new SARS CoV-2 viral isolate spreading rapidly in Southern Africa was announced, containing far more mutations in Spike (S) than previously reported variants. Neutralization titres of Omicron by sera from vaccinees and convalescent subjects infected with early pandemic as well as Alpha, Beta, Gamma, Delta are substantially reduced or fail to neutralize. Titres against Omicron are boosted by third vaccine doses and are high in cases both vaccinated and infected by Delta. Mutations in Omicron knock out or substantially reduce neutralization by most of a large panel of potent monoclonal antibodies and antibodies under commercial development. Omicron S has structural changes from earlier viruses, combining mutations conferring tight binding to ACE2 to unleash evolution driven by immune escape, leading to a large number of mutations in the ACE2 binding site which rebalance receptor affinity to that of early pandemic viruses.

Clinical and genomic signatures of rising SARS-CoV-2 Delta breakthrough infections in New York

In 2021, Delta has become the predominant SARS-CoV-2 variant worldwide. While vaccines effectively prevent COVID-19 hospitalization and death, vaccine breakthrough infections increasingly occur. The precise role of clinical and genomic determinants in Delta infections is not known, and whether they contribute to increased rates of breakthrough infections compared to unvaccinated controls. Here, we show a steep and near complete replacement of circulating variants with Delta between May and August 2021 in metropolitan New York. We observed an increase of the Delta sublineage AY.25, its spike mutation S112L, and nsp12 mutation F192V in breakthroughs. Delta infections were associated with younger age and lower hospitalization rates than Alpha. Delta breakthroughs increased significantly with time since vaccination, and, after adjusting for confounders, they rose at similar rates as in unvaccinated individuals. Our data indicate a limited impact of vaccine escape in favor of Delta's increased epidemic growth in times of waning vaccine protection.

Visualizing in deceased COVID-19 patients how SARS-CoV-2 attacks the respiratory and olfactory mucosae but spares the olfactory bulb

Anosmia, the loss of smell, is a common and often the sole symptom of COVID-19. The onset of the sequence of pathobiological events leading to olfactory dysfunction remains obscure. Here, we have developed a postmortem bedside surgical procedure to harvest endoscopically samples of respiratory and olfactory mucosae and whole olfactory bulbs. Our cohort of 85 cases included COVID-19 patients who died a few days after infection with SARS-CoV-2, enabling us to catch the virus while it was still replicating. We found that sustentacular cells are the major target cell type in the olfactory mucosa. We failed to find evidence for infection of olfactory sensory neurons, and the parenchyma of the olfactory bulb is spared as well. Thus, SARS-CoV-2 does not appear to be a neurotropic virus. We postulate that transient insufficient support from sustentacular cells triggers transient olfactory dysfunction in COVID-19. Olfactory sensory neurons would become affected without getting infected.

The impact of viral mutations on recognition by SARS-CoV-2 specific T cells

We identify amino acid variants within dominant SARS-CoV-2 T cell epitopes by interrogating global sequence data. Several variants within nucleocapsid and ORF3a epitopes have arisen independently in multiple lineages and result in loss of recognition by epitope-specific T cells assessed by IFN-γ and cytotoxic killing assays. Complete loss of T cell responsiveness was seen due to Q213K in the A∗01:01-restricted CD8+ ORF3a epitope FTSDYYQLY207-215; due to P13L, P13S, and P13T in the B∗27:05-restricted CD8+ nucleocapsid epitope QRNAPRITF9-17; and due to T362I and P365S in the A∗03:01/A∗11:01-restricted CD8+ nucleocapsid epitope KTFPPTEPK361-369. CD8+ T cell lines unable to recognize variant epitopes have diverse T cell receptor repertoires. These data demonstrate the potential for T cell evasion and highlight the need for ongoing surveillance for variants capable of escaping T cell as well as humoral immunity.

SARS-CoV-2 variants associated with vaccine breakthrough in the Delaware Valley through summer 2021

The severe acute respiratory coronavirus-2 (SARS-CoV-2) is the cause of the global outbreak of COVID-19. Evidence suggests that the virus is evolving to allow efficient spread through the human population, including vaccinated individuals. Here we report a study of viral variants from surveillance of the Delaware Valley, including the city of Philadelphia, and variants infecting vaccinated subjects. We sequenced and analyzed complete viral genomes from 2621 surveillance samples from March 2020 to September 2021 and compared them to genome sequences from 159 vaccine breakthroughs. In the early spring of 2020, all detected variants were of the B.1 and closely related lineages. A mixture of lineages followed, notably including B.1.243 followed by B.1.1.7 (alpha), with other lineages present at lower levels. Later isolations were dominated by B.1.617.2 (delta) and other delta lineages; delta was the exclusive variant present by the last time sampled. To investigate whether any variants appeared preferentially in vaccine breakthroughs, we devised a model based on Bayesian autoregressive moving average logistic multinomial regression to allow rigorous comparison. This revealed that B.1.617.2 (delta) showed three-fold enrichment in vaccine breakthrough cases (odds ratio of 3; 95% credible interval 0.89-11). Viral point substitutions could also be associated with vaccine breakthroughs, notably the N501Y substitution found in the alpha, beta and gamma variants (odds ratio 2.04; 95% credible interval of 1.25-3.18). This study thus provides a detailed picture of viral evolution in the Delaware Valley and a geographically matched analysis of vaccine breakthroughs; it also introduces a rigorous statistical approach to interrogating enrichment of viral variants.

Beta RBD boost broadens antibody-mediated protection against SARS-CoV-2 variants in animal models

Severe acute respiratory syndrome-coronavirus 2 (SARS-CoV-2) variants of concern (VOCs) with resistance to neutralizing antibodies are threatening to undermine vaccine efficacy. Vaccination and infection have led to widespread humoral immunity against the pandemic founder (Wu-Hu-1). Against this background, it is critical to assess the outcomes of subsequent immunization with variant antigens. It is not yet clear whether heterotypic boosts would be compromised by original antigenic sin, where pre-existing responses to a prior variant dampen responses to a new one, or whether the memory B cell repertoire would bridge the gap between Wu-Hu-1 and VOCs. We show, in macaques immunized with Wu-Hu-1 spike, that a single dose of adjuvanted beta variant receptor binding domain (RBD) protein broadens neutralizing antibody responses to heterologous VOCs. Passive transfer of plasma sampled after Wu-Hu-1 spike immunization only partially protects K18-hACE2 mice from lethal challenge with a beta variant isolate, whereas plasma sampled following heterotypic RBD boost protects completely against disease.

A Strategy to Detect Emerging Non-Delta SARS-CoV-2 Variants with a Monoclonal Antibody Specific for the N501 Spike Residue

Efforts to control SARS-CoV-2 have been challenged by the emergence of variant strains that have important implications for clinical and epidemiological decision making. Four variants of concern (VOCs) have been designated by the Centers for Disease Control and Prevention (CDC), namely, B.1.617.2 (delta), B.1.1.7 (alpha), B.1.351 (beta), and P.1 (gamma), although the last three have been downgraded to variants being monitored (VBMs). VOCs and VBMs have shown increased transmissibility and/or disease severity, resistance to convalescent SARS-CoV-2 immunity and antibody therapeutics, and the potential to evade diagnostic detection. Methods are needed for point-of-care (POC) testing to rapidly identify these variants, protect vulnerable populations, and improve surveillance. Antigen-detection rapid diagnostic tests (Ag-RDTs) are ideal for POC use, but Ag-RDTs that recognize specific variants have not yet been implemented. Here, we describe a mAb (2E8) that is specific for the SARS-CoV-2 spike protein N501 residue. The 2E8 mAb can distinguish the delta VOC from variants with the N501Y meta-signature, which is characterized by convergent mutations that contribute to increased virulence and evasion of host immunity. Among the N501Y-containing mutants formerly designated as VOCs (alpha, beta, and gamma), a previously described mAb, CB6, can distinguish beta from alpha and gamma. When used in a sandwich ELISA, these mAbs sort these important SARS-CoV-2 variants into three diagnostic categories, namely, (1) delta, (2) alpha or gamma, and (3) beta. As delta is currently the predominant variant globally, they will be useful for POC testing to identify N501Y meta-signature variants, protect individuals in high-risk settings, and help detect epidemiological shifts among SARS-CoV-2 variants.

Running With Scissors: Evolutionary Conflicts Between Viral Proteases and the Host Immune System

Many pathogens encode proteases that serve to antagonize the host immune system. In particular, viruses with a positive-sense single-stranded RNA genome [(+)ssRNA], including picornaviruses, flaviviruses, and coronaviruses, encode proteases that are not only required for processing viral polyproteins into functional units but also manipulate crucial host cellular processes through their proteolytic activity. Because these proteases must cleave numerous polyprotein sites as well as diverse host targets, evolution of these viral proteases is expected to be highly constrained. However, despite this strong evolutionary constraint, mounting evidence suggests that viral proteases such as picornavirus 3C, flavivirus NS3, and coronavirus 3CL, are engaged in molecular 'arms races' with their targeted host factors, resulting in host- and virus-specific determinants of protease cleavage. In cases where protease-mediated cleavage results in host immune inactivation, recurrent host gene evolution can result in avoidance of cleavage by viral proteases. In other cases, such as recently described examples in NLRP1 and CARD8, hosts have evolved 'tripwire' sequences that mimic protease cleavage sites and activate an immune response upon cleavage. In both cases, host evolution may be responsible for driving viral protease evolution, helping explain why viral proteases and polyprotein sites are divergent among related viruses despite such strong evolutionary constraint. Importantly, these evolutionary conflicts result in diverse protease-host interactions even within closely related host and viral species, thereby contributing to host range, zoonotic potential, and pathogenicity of viral infection. Such examples highlight the importance of examining viral protease-host interactions through an evolutionary lens.

S-farnesylation is essential for antiviral activity of the long ZAP isoform against RNA viruses with diverse replication strategies

The zinc finger antiviral protein (ZAP) is a broad inhibitor of virus replication. Its best-characterized function is to bind CpG dinucleotides present in viral RNAs and, through the recruitment of TRIM25, KHNYN and other cofactors, target them for degradation or prevent their translation. The long and short isoforms of ZAP (ZAP-L and ZAP-S) have different intracellular localization and it is unclear how this regulates their antiviral activity against viruses with different sites of replication. Using ZAP-sensitive and ZAP-insensitive human immunodeficiency virus type I (HIV-1), which transcribe the viral RNA in the nucleus and assemble virions at the plasma membrane, we show that the catalytically inactive poly-ADP-ribose polymerase (PARP) domain in ZAP-L is essential for CpG-specific viral restriction. Mutation of a crucial cysteine in the C-terminal CaaX box that mediates S-farnesylation and, to a lesser extent, the residues in place of the catalytic site triad within the PARP domain, disrupted the activity of ZAP-L. Addition of the CaaX box to ZAP-S partly restored antiviral activity, explaining why ZAP-S lacks antiviral activity for CpG-enriched HIV-1 despite conservation of the RNA-binding domain. Confocal microscopy confirmed the CaaX motif mediated localization of ZAP-L to vesicular structures and enhanced physical association with intracellular membranes. Importantly, the PARP domain and CaaX box together jointly modulate the interaction between ZAP-L and its cofactors TRIM25 and KHNYN, implying that its proper subcellular localisation is required to establish an antiviral complex. The essential contribution of the PARP domain and CaaX box to ZAP-L antiviral activity was further confirmed by inhibition of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) replication, which replicates in double-membrane vesicles derived from the endoplasmic reticulum. Thus, compartmentalization of ZAP-L on intracellular membranes provides an essential effector function in ZAP-L-mediated antiviral activity against divergent viruses with different subcellular replication sites.