A selective sweep in the Spike gene has driven SARS-CoV-2 human adaptation

Dynamic analysis of SARS-CoV-2 evolution based on different countries

Since late 2019, COVID-19 has significantly impacted the world. Understanding the evolution of SARS-CoV-2 is crucial for protecting against future infectious pathogens. In this study, we conducted a comprehensive chronological analysis of SARS-CoV-2 evolution by examining mutation prevalence from the source countries of VOCs: United Kingdom, India, Brazil, South Africa, plus two countries: United States, Russia, utilizing genomic sequences from GISAID. Our methodological approach involved large-scale genomic sequence alignment using MAFFT, Python-based data processing on a high-performance computing platform, and advanced statistical methods the Maximal Information Coefficient (MIC), and also Long Short-Term Memory (LSTM) models for correlation analysis. Our findings elucidate the dynamics of SARS-CoV-2 evolution, highlighting the virus's changing behaviour over various pandemic stages. Key results include the discovery of three temporal mutation patterns-lineage distinct, long-span, and competitive mutations-with varying levels of impact on the virus. Notably, we observed a convergence of advantageous mutations in the spike protein, especially in the later stages of the pandemic, indicating a substantial evolutionary pressure on the virus. One of the most significant revelations is the predominant role of natural immunity over vaccination-induced immunity in driving these evolutionary changes. This emphasizes the critical need for regular vaccine updates to maintain efficacy against evolving strains. In conclusion, our study not only sheds light on the evolutionary trajectory of SARS-CoV-2 but also underscores the urgency for robust, continuous global data collection and sharing. It highlights the necessity for rapid adaptations in medical countermeasures, including vaccine development, to stay ahead of pathogen evolution. This research provides valuable insights for future pandemic preparedness and response strategies.

Role of N343 glycosylation on the SARS-CoV-2 S RBD structure and co-receptor binding across variants of concern

Glycosylation of the SARS-CoV-2 spike (S) protein represents a key target for viral evolution because it affects both viral evasion and fitness. Successful variations in the glycan shield are difficult to achieve though, as protein glycosylation is also critical to folding and structural stability. Within this framework, the identification of glycosylation sites that are structurally dispensable can provide insight into the evolutionary mechanisms of the shield and inform immune surveillance. In this work, we show through over 45 μs of cumulative sampling from conventional and enhanced molecular dynamics (MD) simulations, how the structure of the immunodominant S receptor binding domain (RBD) is regulated by N-glycosylation at N343 and how this glycan's structural role changes from WHu-1, alpha (B.1.1.7), and beta (B.1.351), to the delta (B.1.617.2), and omicron (BA.1 and BA.2.86) variants. More specifically, we find that the amphipathic nature of the N-glycan is instrumental to preserve the structural integrity of the RBD hydrophobic core and that loss of glycosylation at N343 triggers a specific and consistent conformational change. We show how this change allosterically regulates the conformation of the receptor binding motif (RBM) in the WHu-1, alpha, and beta RBDs, but not in the delta and omicron variants, due to mutations that reinforce the RBD architecture. In support of these findings, we show that the binding of the RBD to monosialylated ganglioside co-receptors is highly dependent on N343 glycosylation in the WHu-1, but not in the delta RBD, and that affinity changes significantly across VoCs. Ultimately, the molecular and functional insight we provide in this work reinforces our understanding of the role of glycosylation in protein structure and function and it also allows us to identify the structural constraints within which the glycosylation site at N343 can become a hotspot for mutations in the SARS-CoV-2 S glycan shield.

Cell type-specific adaptation of the SARS-CoV-2 spike

Severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) can infect various human tissues and cell types, principally via interaction with its cognate receptor angiotensin-converting enzyme-2 (ACE2). However, how the virus evolves in different cellular environments is poorly understood. Here, we used experimental evolution to study the adaptation of the SARS-CoV-2 spike to four human cell lines expressing different levels of key entry factors. After twenty passages of a spike-expressing recombinant vesicular stomatitis virus (VSV), cell-type-specific phenotypic changes were observed and sequencing allowed the identification of sixteen adaptive spike mutations. We used VSV pseudotyping to measure the entry efficiency, ACE2 affinity, spike processing, TMPRSS2 usage, and entry pathway usage of all the mutants, alone or in combination. The fusogenicity of the mutant spikes was assessed with a cell-cell fusion assay. Finally, mutant recombinant VSVs were used to measure the fitness advantage associated with selected mutations. We found that the effects of these mutations varied across cell types, both in terms of viral entry and replicative fitness. Interestingly, two spike mutations (L48S and A372T) that emerged in cells expressing low ACE2 levels increased receptor affinity, syncytia induction, and entry efficiency under low-ACE2 conditions. Our results demonstrate specific adaptation of the SARS-CoV-2 spike to different cell types and have implications for understanding SARS-CoV-2 tissue tropism and evolution.

A Novel Approach Utilizing Domain Adversarial Neural Networks for the Detection and Classification of Selective Sweeps

The identification and classification of selective sweeps are of great significance for improving the understanding of biological evolution and exploring opportunities for precision medicine and genetic improvement. Here, a domain adaptation sweep detection and classification (DASDC) method is presented to balance the alignment of two domains and the classification performance through a domain-adversarial neural network and its adversarial learning modules. DASDC effectively addresses the issue of mismatch between training data and real genomic data in deep learning models, leading to a significant improvement in its generalization capability, prediction robustness, and accuracy. The DASDC method demonstrates improved identification performance compared to existing methods and excels in classification performance, particularly in scenarios where there is a mismatch between application data and training data. The successful implementation of DASDC in real data of three distinct species highlights its potential as a useful tool for identifying crucial functional genes and investigating adaptive evolutionary mechanisms, particularly with the increasing availability of genomic data.

The role of N-glycosylation in spike antigenicity for the SARS-CoV-2 gamma variant

The emergence of SARS-CoV-2 variants alters the efficacy of existing immunity towards the viral spike protein, whether acquired from infection or vaccination. Mutations that impact N-glycosylation of spike may be particularly important in influencing antigenicity, but their consequences are difficult to predict. Here, we compare the glycosylation profiles and antigenicity of recombinant viral spike of ancestral Wu-1 and the Gamma strain, which has two additional N-glycosylation sites due to amino acid substitutions in the N-terminal domain (NTD). We found that a mutation at residue 20 from threonine to asparagine within the NTD caused the loss of NTD-specific antibody COVA2-17 binding. Glycan site-occupancy analyses revealed that the mutation resulted in N-glycosylation switching to the new sequon at N20 from the native N17 site. Site-specific glycosylation profiles demonstrated distinct glycoform differences between Wu-1, Gamma, and selected NTD variant spike proteins, but these did not affect antibody binding. Finally, we evaluated the specificity of spike proteins against convalescent COVID-19 sera and found reduced cross-reactivity against some mutants, but not Gamma spike compared to Wuhan spike. Our results illustrate the impact of viral divergence on spike glycosylation and SARS-CoV-2 antibody binding profiles.

Molecular insights into the adaptive evolution of SARS-CoV-2 spike protein

The COVID-19 pandemic, caused by the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), has substantially damaged the global economy and human health. The spike (S) protein of coronaviruses plays a pivotal role in viral entry by binding to host cell receptors. Additionally, it acts as the primary target for neutralizing antibodies in those infected and is the central focus for currently utilized or researched vaccines. During the virus's adaptation to the human host, the S protein of SARS-CoV-2 has undergone significant evolution. As the COVID-19 pandemic has unfolded, new mutations have arisen and vanished, giving rise to distinctive amino acid profiles within variant of concern strains of SARS-CoV-2. Notably, many of these changes in the S protein have been positively selected, leading to substantial alterations in viral characteristics, such as heightened transmissibility and immune evasion capabilities. This review aims to provide an overview of our current understanding of the structural implications associated with key amino acid changes in the S protein of SARS-CoV-2. These research findings shed light on the intricate and dynamic nature of viral evolution, underscoring the importance of continuous monitoring and analysis of viral genomes. Through these molecular-level investigations, we can attain deeper insights into the virus's adaptive evolution, offering valuable guidance for designing vaccines and developing antiviral drugs to combat the ever-evolving viral threats.

Cellular receptors for mammalian viruses

The interaction of viral surface components with cellular receptors and other entry factors determines key features of viral infection such as host range, tropism and virulence. Despite intensive research, our understanding of these interactions remains limited. Here, we report a systematic analysis of published work on mammalian virus receptors and attachment factors. We build a dataset twice the size of those available to date and specify the role of each factor in virus entry. We identify cellular proteins that are preferentially used as virus receptors, which tend to be plasma membrane proteins with a high propensity to interact with other proteins. Using machine learning, we assign cell surface proteins a score that predicts their ability to function as virus receptors. Our results also reveal common patterns of receptor usage among viruses and suggest that enveloped viruses tend to use a broader repertoire of alternative receptors than non-enveloped viruses, a feature that might confer them with higher interspecies transmissibility.

Chronic shedding of a SARS-CoV-2 Alpha variant in wastewater

BACKGROUND

Central Michigan University (CMU) participated in a state-wide SARS-CoV-2 wastewater monitoring program since 2021. Wastewater samples were collected from on-campus sites and nine off-campus wastewater treatment plants servicing small metropolitan and rural communities. SARS-CoV-2 genome copies were quantified using droplet digital PCR and results were reported to the health department.

RESULTS

One rural, off-campus site consistently produced higher concentrations of SARS-CoV-2 genome copies. Samples from this site were sequenced and contained predominately a derivative of Alpha variant lineage B.1.1.7, detected from fall 2021 through summer 2023. Mutational analysis of reconstructed genes revealed divergence from the Alpha variant lineage sequence over time, including numerous mutations  in the Spike RBD and NTD.

CONCLUSIONS

We discuss the possibility that a chronic SARS-CoV-2 infection accumulated adaptive mutations that promoted long-term infection. This study reveals that small wastewater treatment plants can enhance resolution of rare events and facilitate reconstruction of viral genomes due to the relative lack of contaminating sequences.

Allelic gene conversion softens selective sweeps

The prominence of positive selection, in which beneficial mutations are favored by natural selection and rapidly increase in frequency, is a subject of intense debate. Positive selection can result in selective sweeps, in which the haplotype(s) bearing the adaptive allele "sweep" through the population, thereby removing much of the genetic diversity from the region surrounding the target of selection. Two models of selective sweeps have been proposed: classical sweeps, or "hard sweeps", in which a single copy of the adaptive allele sweeps to fixation, and "soft sweeps", in which multiple distinct copies of the adaptive allele leave descendants after the sweep. Soft sweeps can be the outcome of recurrent mutation to the adaptive allele, or the presence of standing genetic variation consisting of multiple copies of the adaptive allele prior to the onset of selection. Importantly, soft sweeps will be common when populations can rapidly adapt to novel selective pressures, either because of a high mutation rate or because adaptive alleles are already present. The prevalence of soft sweeps is especially controversial, and it has been noted that selection on standing variation or recurrent mutations may not always produce soft sweeps. Here, we show that the inverse is true: selection on single-origin de novo mutations may often result in an outcome that is indistinguishable from a soft sweep. This is made possible by allelic gene conversion, which "softens" hard sweeps by copying the adaptive allele onto multiple genetic backgrounds, a process we refer to as a "pseudo-soft" sweep. We carried out a simulation study examining the impact of gene conversion on sweeps from a single de novo variant in models of human, Drosophila, and Arabidopsis populations. The fraction of simulations in which gene conversion had produced multiple haplotypes with the adaptive allele upon fixation was appreciable. Indeed, under realistic demographic histories and gene conversion rates, even if selection always acts on a single-origin mutation, sweeps involving multiple haplotypes are more likely than hard sweeps in large populations, especially when selection is not extremely strong. Thus, even when the mutation rate is low or there is no standing variation, hard sweeps are expected to be the exception rather than the rule in large populations. These results also imply that the presence of signatures of soft sweeps does not necessarily mean that adaptation has been especially rapid or is not mutation limited.

PARNAS: Objectively Selecting the Most Representative Taxa on a Phylogeny

The use of next-generation sequencing technology has enabled phylogenetic studies with hundreds of thousands of taxa. Such large-scale phylogenies have become a critical component in genomic epidemiology in pathogens such as SARS-CoV-2 and influenza A virus. However, detailed phenotypic characterization of pathogens or generating a computationally tractable dataset for detailed phylogenetic analyses requires objective subsampling of taxa. To address this need, we propose parnas, an objective and flexible algorithm to sample and select taxa that best represent observed diversity by solving a generalized k-medoids problem on a phylogenetic tree. parnas solves this problem efficiently and exactly by novel optimizations and adapting algorithms from operations research. For more nuanced selections, taxa can be weighted with metadata or genetic sequence parameters, and the pool of potential representatives can be user-constrained. Motivated by influenza A virus genomic surveillance and vaccine design, parnas can be applied to identify representative taxa that optimally cover the diversity in a phylogeny within a specified distance radius. We demonstrated that parnas is more efficient and flexible than existing approaches. To demonstrate its utility, we applied parnas to 1) quantify SARS-CoV-2 genetic diversity over time, 2) select representative influenza A virus in swine genes derived from over 5 years of genomic surveillance data, and 3) identify gaps in H3N2 human influenza A virus vaccine coverage. We suggest that our method, through the objective selection of representatives in a phylogeny, provides criteria for quantifying genetic diversity that has application in the the rational design of multivalent vaccines and genomic epidemiology. PARNAS is available at https://github.com/flu-crew/parnas.

Tensor Decomposition-based Feature Extraction and Classification to Detect Natural Selection from Genomic Data

Inferences of adaptive events are important for learning about traits, such as human digestion of lactose after infancy and the rapid spread of viral variants. Early efforts toward identifying footprints of natural selection from genomic data involved development of summary statistic and likelihood methods. However, such techniques are grounded in simple patterns or theoretical models that limit the complexity of settings they can explore. Due to the renaissance in artificial intelligence, machine learning methods have taken center stage in recent efforts to detect natural selection, with strategies such as convolutional neural networks applied to images of haplotypes. Yet, limitations of such techniques include estimation of large numbers of model parameters under nonconvex settings and feature identification without regard to location within an image. An alternative approach is to use tensor decomposition to extract features from multidimensional data although preserving the latent structure of the data, and to feed these features to machine learning models. Here, we adopt this framework and present a novel approach termed T-REx, which extracts features from images of haplotypes across sampled individuals using tensor decomposition, and then makes predictions from these features using classical machine learning methods. As a proof of concept, we explore the performance of T-REx on simulated neutral and selective sweep scenarios and find that it has high power and accuracy to discriminate sweeps from neutrality, robustness to common technical hurdles, and easy visualization of feature importance. Therefore, T-REx is a powerful addition to the toolkit for detecting adaptive processes from genomic data.

Remdesivir increases mtDNA copy number causing mild alterations to oxidative phosphorylation

SARS-CoV-2 causes the severe respiratory disease COVID-19. Remdesivir (RDV) was the first fast-tracked FDA approved treatment drug for COVID-19. RDV acts as an antiviral ribonucleoside (adenosine) analogue that becomes active once it accumulates intracellularly. It then diffuses into the host cell and terminates viral RNA transcription. Previous studies have shown that certain nucleoside analogues unintentionally inhibit mitochondrial RNA or DNA polymerases or cause mutational changes to mitochondrial DNA (mtDNA). These past findings on the mitochondrial toxicity of ribonucleoside analogues motivated us to investigate what effects RDV may have on mitochondrial function. Using in vitro and in vivo rodent models treated with RDV, we observed increases in mtDNA copy number in Mv1Lu cells (35.26% increase ± 11.33%) and liver (100.27% increase ± 32.73%) upon treatment. However, these increases only resulted in mild changes to mitochondrial function. Surprisingly, skeletal muscle and heart were extremely resistant to RDV treatment, tissues that have preferentially been affected by other nucleoside analogues. Although our data suggest that RDV does not greatly impact mitochondrial function, these data are insightful for the treatment of RDV for individuals with mitochondrial disease.

Mechanism and evolution of human ACE2 binding by SARS-CoV-2 spike

Spike glycoprotein of SARS-CoV-2 mediates viral entry into host cells by facilitating virus attachment and membrane fusion. ACE2 is the main receptor of SARS-CoV-2 and its interaction with spike has shaped the virus' emergence from an animal reservoir and subsequent evolution in the human host. Many structural studies on the spike:ACE2 interaction have provided insights into mechanisms driving viral evolution during the on-going pandemic. This review describes the molecular basis of spike binding to ACE2, outlines mechanisms that have optimised this interaction during viral evolution, and suggests directions for future research.

Host susceptibility and structural and immunological insight of S proteins of two SARS-CoV-2 closely related bat coronaviruses

The bat coronaviruses (CoV) BANAL-20-52 and BANAL-20-236 are two newly identified severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) closely related coronaviruses (SC2r-CoV) and the genome of BANAL-20-52 shares the highest homology with SARS-CoV-2. However, the risk of their potential zoonotic transmission has not been fully evaluated. Here, we determined their potential host susceptibility among 13 different bat species and 26 different animal species, and found that both might have extensive host ranges, indicating high zoonotic transmission potential. We also determined the cryo-EM structures of BANAL-20-52 and BANAL-20-236 S proteins at pH 5.5 and the complex of BANAL-20-236 S1 and Rhinolophus affinis ACE2, and found that both trimeric S proteins adopt all three receptor binding domains (RBDs) in "closed" conformation and are more compact than SARS-CoV-2. Strikingly, the unique sugar moiety at N370 of bat SC2r-CoVs acts like a "bolt" and crosses over two neighboring subunits, facilitating the S proteins in the locked conformation and underpinning the architecture stability. Removal of the glycosylation at N370 by a T372A substitution substantially enhances virus infectivity but becomes highly sensitive to trypsin digestion at pH 5.5, a condition roughly mimicking the insectivorous bat's stomach digestion. In contrast, WT S proteins of SC2r-CoVs showed considerable resistance to trypsin digestion at pH 5.5, indicating that the highly conserved T372 in bat CoVs might result from the selective advantages in stability during the fecal-oral transmission over A372. Moreover, the results of cross-immunogenicity among S proteins of SARS-CoV-2, BANAL-20-52, and BANAL-20-236 showed that A372 pseudoviruses are more sensitive to anti-S sera than T372, indicating that immune evasion might also play a role in the natural selection of T372 over A372 during evolution. Finally, residues 493 and 498 of the S protein affect host susceptibility, and residue 498 also influences the immunogenicity of the S protein. Together, our findings aid a better understanding of the molecular basis of CoV entry, selective evolution, and immunogenicity and highlight the importance of surveillance of susceptible hosts of these viruses to prevent potential outbreaks.

Lectin Fingerprinting Distinguishes Antibody Neutralization in SARS-CoV-2

Enveloped viruses co-opt host glycosylation pathways to decorate their surface proteins. As viruses evolve, emerging strains can modify their glycosylation patterns to influence host interactions and subvert immune recognition. Still, changes in viral glycosylation or their impact on antibody protection cannot be predicted from genomic sequences alone. Using the highly glycosylated SARS-CoV-2 Spike protein as a model system, we present a lectin fingerprinting method that rapidly reports on changes in variant glycosylation state, which are linked to antibody neutralization. In the presence of antibodies or convalescent and vaccinated patient sera, unique lectin fingerprints emerge that distinguish neutralizing versus non-neutralizing antibodies. This information could not be inferred from direct binding interactions between antibodies and the Spike receptor-binding domain (RBD) binding data alone. Comparative glycoproteomics of the Spike RBD of wild-type (Wuhan-Hu-1) and Delta (B.1.617.2) variants reveal O-glycosylation differences as a key determinant of immune recognition differences. These data underscore the interplay between viral glycosylation and immune recognition and reveal lectin fingerprinting to be a rapid, sensitive, and high-throughput assay to distinguish the neutralization potential of antibodies that target critical viral glycoproteins.

The E2 glycoprotein holds key residues for Mayaro virus adaptation to the urban Aedes aegypti mosquito

Adaptation to mosquito vectors suited for transmission in urban settings is a major driver in the emergence of arboviruses. To better anticipate future emergence events, it is crucial to assess their potential to adapt to new vector hosts. In this work, we used two different experimental evolution approaches to study the adaptation process of an emerging alphavirus, Mayaro virus (MAYV), to Ae. aegypti, an urban mosquito vector of many other arboviruses. We identified E2-T179N as a key mutation increasing MAYV replication in insect cells and enhancing transmission after escaping the midgut of live Ae. aegypti. In contrast, this mutation decreased viral replication and binding in human fibroblasts, a primary cellular target of MAYV in humans. We also showed that MAYV E2-T179N generates reduced viremia and displays less severe tissue pathology in vivo in a mouse model. We found evidence in mouse fibroblasts that MAYV E2-T179N is less dependent on the Mxra8 receptor for replication than WT MAYV. Similarly, exogenous expression of human apolipoprotein receptor 2 and Mxra8 enhanced WT MAYV replication compared to MAYV E2-T179N. When this mutation was introduced in the closely related chikungunya virus, which has caused major outbreaks globally in the past two decades, we observed increased replication in both human and insect cells, suggesting E2 position 179 is an important determinant of alphavirus host-adaptation, although in a virus-specific manner. Collectively, these results indicate that adaptation at the T179 residue in MAYV E2 may result in increased vector competence-but coming at the cost of optimal replication in humans-and may represent a first step towards a future emergence event.

Tensor decomposition based feature extraction and classification to detect natural selection from genomic data

Inferences of adaptive events are important for learning about traits, such as human digestion of lactose after infancy and the rapid spread of viral variants. Early efforts toward identifying footprints of natural selection from genomic data involved development of summary statistic and likelihood methods. However, such techniques are grounded in simple patterns or theoretical models that limit the complexity of settings they can explore. Due to the renaissance in artificial intelligence, machine learning methods have taken center stage in recent efforts to detect natural selection, with strategies such as convolutional neural networks applied to images of haplotypes. Yet, limitations of such techniques include estimation of large numbers of model parameters under non-convex settings and feature identification without regard to location within an image. An alternative approach is to use tensor decomposition to extract features from multidimensional data while preserving the latent structure of the data, and to feed these features to machine learning models. Here, we adopt this framework and present a novel approach termed T-REx , which extracts features from images of haplotypes across sampled individuals using tensor decomposition, and then makes predictions from these features using classical machine learning methods. As a proof of concept, we explore the performance of T-REx on simulated neutral and selective sweep scenarios and find that it has high power and accuracy to discriminate sweeps from neutrality, robustness to common technical hurdles, and easy visualization of feature importance. Therefore, T-REx is a powerful addition to the toolkit for detecting adaptive processes from genomic data.

SARS-CoV-2 variant evasion of monoclonal antibodies based on in vitro studies

Monoclonal antibodies (mAbs) offer a treatment option for individuals with severe COVID-19 and are especially important in high-risk individuals where vaccination is not an option. Given the importance of understanding the evolution of resistance to mAbs by SARS-CoV-2, we reviewed the available in vitro neutralization data for mAbs against live variants and viral constructs containing spike mutations of interest. Unfortunately, evasion of mAb-induced protection is being reported with new SARS-CoV-2 variants. The magnitude of neutralization reduction varied greatly among mAb-variant pairs. For example, sotrovimab retained its neutralization capacity against Omicron BA.1 but showed reduced efficacy against BA.2, BA.4 and BA.5, and BA.2.12.1. At present, only bebtelovimab has been reported to retain its efficacy against all SARS-CoV-2 variants considered here. Resistance to mAb neutralization was dominated by the action of epitope single amino acid substitutions in the spike protein. Although not all observed epitope mutations result in increased mAb evasion, amino acid substitutions at non-epitope positions and combinations of mutations also contribute to evasion of neutralization. This Review highlights the implications for the rational design of viral genomic surveillance and factors to consider for the development of novel mAb therapies.

Effect of polymorphism in Rhinolophus affinis ACE2 on entry of SARS-CoV-2 related bat coronaviruses

Bat coronavirus RaTG13 shares about 96.2% nucleotide sequence identity with that of SARS-CoV-2 and uses human and Rhinolophus affinis (Ra) angiotensin-converting enzyme 2 (ACE2) as entry receptors. Whether there are bat species other than R. affinis susceptible to RaTG13 infection remains elusive. Here, we show that, among 18 different bat ACE2s tested, only RaACE2 is highly susceptible to transduction by RaTG13 S pseudovirions, indicating that the bat species harboring RaTG13 might be very limited. RaACE2 has seven polymorphic variants, RA-01 to RA-07, and they show different susceptibilities to RaTG13 S pseudovirions transduction. Sequence and mutagenesis analyses reveal that residues 34, 38, and 83 in RaACE2 might play critical roles in interaction with the RaTG13 S protein. Of note, RaACE2 polymorphisms have minimal effect on S proteins of SARS-CoV-2 and several SARS-CoV-2 related CoVs (SC2r-CoVs) including BANAL-20-52 and BANAL-20-236 in terms of binding, membrane fusion, and pseudovirus entry. Further mutagenesis analyses identify residues 501 and 505 in S proteins critical for the recognition of different RaACE2 variants and pangolin ACE2 (pACE2), indicating that RaTG13 might have not been well adapted to R. affinis bats. While single D501N and H505Y changes in RaTG13 S protein significantly enhance the infectivity and minimize the difference in susceptibility among different RaACE2 variants, an N501D substitution in SARS-CoV-2 S protein displays marked disparity in transduction efficiencies among RaACE2 variants with a significant reduction in infectivity on several RaACE2 variants. Finally, a T372A substitution in RaTG13 S protein not only significantly increases infectivity on all RaACE2 variants, but also markedly enhances entry on several bat ACE2s including R. sinicus YN, R. pearsonii, and R. ferrumeiqunum. However, the T372A mutant is about 4-fold more sensitive to neutralizing sera from mice immunized with BANAL-20-52 S, suggesting that the better immune evasion ability of T372 over A372 might contribute to the natural selective advantage of T372 over A372 among bat CoVs. Together, our study aids a better understanding of coronavirus entry, vaccine design, and evolution.

Evaluating the effect of SARS-CoV-2 spike mutations with a linear doubly robust learner

Driven by various mutations on the viral Spike protein, diverse variants of SARS-CoV-2 have emerged and prevailed repeatedly, significantly prolonging the pandemic. This phenomenon necessitates the identification of key Spike mutations for fitness enhancement. To address the need, this manuscript formulates a well-defined framework of causal inference methods for evaluating and identifying key Spike mutations to the viral fitness of SARS-CoV-2. In the context of large-scale genomes of SARS-CoV-2, it estimates the statistical contribution of mutations to viral fitness across lineages and therefore identifies important mutations. Further, identified key mutations are validated by computational methods to possess functional effects, including Spike stability, receptor-binding affinity, and potential for immune escape. Based on the effect score of each mutation, individual key fitness-enhancing mutations such as D614G and T478K are identified and studied. From individual mutations to protein domains, this paper recognizes key protein regions on the Spike protein, including the receptor-binding domain and the N-terminal domain. This research even makes further efforts to investigate viral fitness via mutational effect scores, allowing us to compute the fitness score of different SARS-CoV-2 strains and predict their transmission capacity based solely on their viral sequence. This prediction of viral fitness has been validated using BA.2.12.1, which is not used for regression training but well fits the prediction. To the best of our knowledge, this is the first research to apply causal inference models to mutational analysis on large-scale genomes of SARS-CoV-2. Our findings produce innovative and systematic insights into SARS-CoV-2 and promotes functional studies of its key mutations, serving as reliable guidance about mutations of interest.

A New HIV-1 K28E32-Reverse Transcriptase Variant Associated with the Rapid Expansion of CRF07_BC among Men Who Have Sex with Men

HIV-1 CRF07_BC originated among injection drug users (IDUs) in China. After diffusing into men who have sex with men (MSM), CRF07_BC has shown a rapid expansion in this group; however, the mechanism remains unclear. Here, we identified a new K₂₈E₃₂ variant of CRF07_BC that was characterized by five specific mutations (E28K, K32E, E248V, K249Q, and T338S) in reverse transcriptase. This variant was mainly prevalent among MSM, and was overrepresented in transmission clusters, suggesting that it could have driven the rapid expansion of CRF07_BC in MSM, though founder effects cannot be ruled out. It was descended from an evolutionary intermediate accumulating four specific mutations and formed an independent phylogenetic node with an estimated origin time in 2003. The K₂₈E₃₂ variant was demonstrated to have significantly higher in vitro HIV-1 replication ability than the wild type. Mutations E28K and K32E play a critical role in the improvement of in vitro HIV-1 replication ability, reflected by improved reverse transcription activity. The results could allow public health officials to use this marker (especially E28K and K32E mutations in the reverse transcriptase (RT) coding region) to target prevention measures prioritizing MSM population and persons infected with this variant for test and treat initiatives. IMPORTANCE HIV-1 has very high mutation rate that is correlated with the survival and adaption of the virus. The variants with higher transmissibility may be more selective advantage than the strains with higher virulence. Several HIV-1 variants were previously demonstrated to be correlated with higher viral load and lower CD4 T cell count. Here, we first identified a new variant (the K₂₈E₃₂ variant) of HIV-1 CRF07_BC, described its origin and evolutionary dynamics, and demonstrated its higher in vitro HIV-1 replication ability than the wild type. We demonstrated that five RT mutations (especially E28K and K32E) significantly improve in vitro HIV-1 replication ability. The appearance of the new K₂₈E₃₂ variant was associated with the rapidly increasing prevalence of CRF07_BC among MSM.

Contrasting patterns in the early stage of SARS-CoV-2 evolution between humans and minks

One of the unique features of SARS-CoV-2 is its apparent neutral evolution during the early pandemic (before February 2020). This contrasts with the preceding SARS-CoV epidemics, where viruses evolved adaptively. SARS-CoV-2 may exhibit a unique or adaptive feature which deviates from other coronaviruses. Alternatively, the virus may have been cryptically circulating in humans for a sufficient time to have acquired adaptive changes before the onset of the current pandemic. To test the scenarios above, we analyzed the SARS-CoV-2 sequences from minks (Neovision vision) and parental humans. In the early phase of the mink epidemic (April to May 2020), nonsynonymous to synonymous mutation ratio per site in the spike protein is 2.93, indicating a selection process favoring adaptive amino acid changes. Mutations in the spike protein were concentrated within its receptor binding domain and receptor binding motif. An excess of high frequency derived variants produced by genetic hitchhiking was found during the middle (June to July 2020) and late phase I (August to September 2020) of the mink epidemic. In contrast, the site frequency spectra of early SARS-CoV-2 in humans only show an excess of low frequency mutations, consistent with the recent outbreak of the virus. Strong positive selection in the mink SARS-CoV-2 implies the virus may not be pre-adapted to a wide range of hosts and illustrates how a virus evolves to establish a continuous infection in a new host. Therefore, the lack of positive selection signal during the early pandemic in humans deserves further investigation.

The spike glycoprotein of highly pathogenic human coronaviruses: structural insights for understanding infection, evolution and inhibition

Highly pathogenic human coronaviruses (CoV) including SARS‐CoV, MERS‐CoV and SARS‐CoV‐2 have emerged over the past two decades, resulting in infectious disease outbreaks that have greatly affected public health. The CoV surface spike (S) glycoprotein mediates receptor binding and membrane fusion for cell entry, playing critical roles in CoV infection and evolution. The S glycoprotein is also the major target molecule for prophylactic and therapeutic interventions, including neutralizing antibodies and vaccines. In this review, we summarize key studies that have revealed the structural basis of S‐mediated cell entry of SARS‐CoV, MERS‐CoV and SARS‐CoV‐2. Additionally, we discuss the evolution of the S glycoprotein to realize cross‐species transmission from the viewpoint of structural biology. Lastly, we describe the recent progress in developing antibodies, nanobodies and peptide inhibitors that target the SARS‐CoV‐2 S glycoprotein for therapeutic purposes.

Natural Selection, Intracellular Bottlenecks of Virus Populations, and Viral Superinfection Exclusion

Natural selection acts on cellular organisms by ensuring the genes responsible for an advantageous phenotype consistently reap the phenotypic advantage. This is possible because reproductive cells of these organisms are almost always haploid, separating the beneficial gene from its rival allele at every generation. How natural selection acts on plus-strand RNA viruses is unclear because these viruses frequently load host cells with numerous genome copies and replicate thousands of progeny genomes in each cell. Recent studies suggest that these viruses encode the Bottleneck, Isolate, Amplify, Select (BIAS) mechanism that blocks all but a few viral genome copies from replication, thus creating the environment in which the bottleneck-escaping viral genome copies are isolated from each other, allowing natural selection to reward beneficial mutations and purge lethal errors. This BIAS mechanism also blocks the genomes of highly homologous superinfecting viruses, thus explaining cellular-level superinfection exclusion. Expected final online publication date for the Annual Review of Virology, Volume 9 is September 2022. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

Principles of SARS-CoV-2 Glycosylation

The structure and post-translational processing of the SARS-CoV-2 spike glycoprotein (S) is intimately associated with the function of the virus and of sterilising vaccines. The surface of the S protein is extensively modified by glycans, and their biosynthesis is driven by both the wider cellular context, and importantly, the underlining protein structure and local glycan density. Comparison of virally derived S protein with both recombinantly derived and adenovirally induced proteins, reveal hotspots of protein-directed glycosylation that drive conserved glycosylation motifs. Molecular dynamics simulations revealed that, while the S surface is extensively shielded by N-glycans, it presents regions vulnerable to neutralising antibodies. Furthermore, glycans have been shown to influence the accessibility of the receptor binding domain and the binding to the cellular receptor. The emerging picture is one of unifying, principles of S protein glycosylation and an intimate role of glycosylation in immunogen structure and efficacy.

SARS-CoV-2-specific T-cell epitope repertoire in convalescent and mRNA-vaccinated individuals

Continuously emerging variants of concern (VOCs) sustain the SARS-CoV-2 pandemic. The SARS-CoV-2 Omicron/B.1.1.529 VOC harbours multiple mutations in the spike protein associated with high infectivity and efficient evasion from humoral immunity induced by previous infection or vaccination. By performing in-depth comparisons of the SARS-CoV-2-specific T-cell epitope repertoire after infection and messenger RNA vaccination, we demonstrate that spike-derived epitopes were not dominantly targeted in convalescent individuals compared to non-spike epitopes. In vaccinees, however, we detected a broader spike-specific T-cell response compared to convalescent individuals. Booster vaccination increased the breadth of the spike-specific T-cell response in convalescent individuals but not in vaccinees with complete initial vaccination. In convalescent individuals and vaccinees, the targeted T-cell epitopes were broadly conserved between wild-type SARS-CoV-2 variant B and Omicron/B.1.1.529. Hence, our data emphasize the relevance of vaccine-induced spike-specific CD8⁺ T-cell responses in combating VOCs including Omicron/B.1.1.529 and support the benefit of boosting convalescent individuals with mRNA vaccines.

A comparison of the repertoire of SARS-CoV-2-specific epitopes targeted by T cells induced by vaccination or natural infection reveals that T cells predominantly target non-spike epitopes in convalescent individuals, while there is a broader spike-specific CD8⁺ T-cell response in vaccinees. Despite differences in T-cell response, the targeted T-cell epitopes were conserved between the wild-type and Omicron variants in both groups.

Structural Dynamics and Molecular Evolution of the SARS-CoV-2 Spike Protein

The ongoing coronavirus disease 2019 (COVID-19) pandemic demonstrates the threat posed by novel coronaviruses to human health. Coronaviruses share a highly conserved cell entry mechanism mediated by the spike protein, the sole product of the S gene. The structural dynamics by which the spike protein orchestrates infection illuminate how antibodies neutralize virions and how S mutations contribute to viral fitness. Here, we review the process by which spike engages its proteinaceous receptor, angiotensin converting enzyme 2 (ACE2), and how host proteases prime and subsequently enable efficient membrane fusion between virions and target cells. We highlight mutations common among severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) variants of concern and discuss implications for cell entry. Ultimately, we provide a model by which sarbecoviruses are activated for fusion competency and offer a framework for understanding the interplay between humoral immunity and the molecular evolution of the SARS-CoV-2 Spike. In particular, we emphasize the relevance of the Canyon Hypothesis (M. G. Rossmann, J Biol Chem 264:14587-14590, 1989) for understanding evolutionary trajectories of viral entry proteins during sustained intraspecies transmission of a novel viral pathogen.

Transient Expression of Glycosylated SARS-CoV-2 Antigens in Nicotiana benthamiana

The current COVID-19 pandemic very dramatically shows that the world lacks preparedness for novel viral diseases. In addition to newly emerging viruses, many known pathogenic viruses such as influenza are constantly evolving, leading to frequent outbreaks with severe diseases and deaths. Hence, infectious viruses are a recurrent burden to our daily life, and powerful strategies to stop the spread of human pathogens and disease progression are of utmost importance. Transient plant-based protein expression is a technology that allows fast and highly flexible manufacturing of recombinant viral proteins and, thus, can contribute to infectious disease detection and prevention. This review highlights recent progress in the transient production of viral glycoproteins in N. benthamiana with a focus on SARS-CoV-2-derived viral antigens.

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.

SARS-CoV-2 spike E156G/Δ157-158 mutations contribute to increased infectivity and immune escape

This study underscores the significance of NTD-specific SARS-CoV-2 spike mutations E156G/Δ157-158 in determining virion infectivity and neutralization by vaccine-elicited antibodies.

Breakthrough infections by emerging SARS-CoV-2 variants raise significant concerns. Here, we sequence-characterized the spike gene from breakthrough infections that corresponded to B.1.617 sublineage. Delineating the functional impact of spike mutations revealed that N-terminal domain (NTD)-specific E156G/Δ157-158 contributed to increased infectivity and reduced sensitivity to vaccine-induced antibodies. A six-nucleotide deletion (467–472) in the spike-coding region introduced this change in the NTD. We confirmed the presence of E156G/Δ157-158 from cases concurrently screened, in addition to other circulating spike (S1) mutations such as T19R, T95I, L452R, E484Q, and D614G. Notably, E156G/Δ157-158 was present in more than 90% of the sequences reported from the USA and UK in October 2021. The spike-pseudotyped viruses bearing a combination of E156G/Δ157-158 and L452R exhibited higher infectivity and reduced sensitivity to neutralization. Notwithstanding, the post-recovery plasma robustly neutralized viral particles bearing the mutant spike. When the spike harbored E156G/Δ157-158 along with L452R and E484Q, increased cell-to-cell fusion was also observed, suggesting a combinatorial effect of these mutations. Our study underscores the importance of non-RBD changes in determining infectivity and immune escape.

Continuous genomic diversification of long polynucleotide fragments drives the emergence of new SARS-CoV-2 variants of concern

Highly transmissible or immuno-evasive SARS-CoV-2 variants have intermittently emerged, resulting in repeated COVID-19 surges. With over 6 million SARS-CoV-2 genomes sequenced, there is unprecedented data to decipher the evolution of fitter SARS-CoV-2 variants. Much attention has been directed to studying the functional importance of specific mutations in the Spike protein, but there is limited knowledge of genomic signatures shared by dominant variants. Here, we introduce a method to quantify the genome-wide distinctiveness of polynucleotide fragments (3- to 240-mers) that constitute SARS-CoV-2 sequences. Compared to standard phylogenetic metrics and mutational load, the new metric provides improved separation between Variants of Concern (VOCs; Reference = 89, IQR: 65-108; Alpha = 166, IQR: 149-181; Beta 131, IQR: 114-149; Gamma = 164, IQR: 150-178; Delta = 235, IQR: 217-255; and Omicron = 459, IQR: 395-521). Omicron's high genomic distinctiveness may confer an advantage over prior VOCs and the recently emerged and highly mutated B.1.640.2 (IHU) lineage. Evaluation of 883 lineages highlights that genomic distinctiveness has increased over time (R ² = 0.37) and that VOCs score significantly higher than contemporary non-VOC lineages, with Omicron among the most distinctive lineages observed. This study demonstrates the value of characterizing SARS-CoV-2 variants by genome-wide polynucleotide distinctiveness and emphasizes the need to go beyond a narrow set of mutations at known sites on the Spike protein. The consistently higher distinctiveness of each emerging VOC compared to prior VOCs suggests that monitoring of genomic distinctiveness would facilitate rapid assessment of viral fitness.

Induction of Broadly Cross-Reactive Antibodies by Displaying Receptor Binding Domains of SARS-CoV-2 on Virus-like Particles

The impact of the COVID-19 pandemic has been reduced since the application of vaccination programs, mostly shown in the reduction of hospitalized patients. However, the emerging variants, in particular Omicron, have caused a steep increase in the number of infections; this increase is, nevertheless, not matched by an increase in hospitalization. Therefore, a vaccine that induces cross-reactive antibodies against most or all variants is a potential solution for the issue of emerging new variants. Here, we present a vaccine candidate which displays receptor-binding domain (RBD) of SARS-CoV-2 on virus-like particles (VLP) that, in mice, not only induce strong antibody responses against RBD but also bind RBDs from other variants of concern (VOCs). The antibodies induced by wild-type (wt) RBD displayed on immunologically optimized Cucumber mosaic virus incorporated tetanus toxin (CuMVTT) VLPs bind to wt as well as RBDs of VOCs with high avidities, indicating induction of strongly cross-reactive IgG antibodies. Interestingly, similar cross-reactive IgA antibodies were induced in immunized mice. Furthermore, these cross-reactive antibodies demonstrated efficacy in neutralizing wt (Wuhan) as well as SARS-CoV-2 VOCs (Beta, Delta, and Gamma). In summary, RBDs displayed on VLPs are capable of inducing protective cross-reactive IgG and IgA antibodies in mice, indicating that it may be possible to cover emerging VOCs with a single vaccine based on wt RBD.

Evolutionary dynamics of the severe acute respiratory syndrome coronavirus 2 genomes

The coronavirus disease 2019 (COVID-19) pandemic has caused immense losses in human lives and the global economy and posed significant challenges for global public health. As severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the causative agent of COVID-19, has evolved, thousands of single nucleotide variants (SNVs) have been identified across the viral genome. The roles of individual SNVs in the zoonotic origin, evolution, and transmission of SARS-CoV-2 have become the focus of many studies. This review summarizes recent comparative genomic analyses of SARS-CoV-2 and related coronaviruses (SC2r-CoVs) found in non-human animals, including delineation of SARS-CoV-2 lineages based on characteristic SNVs. We also discuss the current understanding of receptor-binding domain (RBD) evolution and characteristic mutations in variants of concern (VOCs) of SARS-CoV-2, as well as possible co-evolution between RBD and its receptor, angiotensin-converting enzyme 2 (ACE2). We propose that the interplay between SARS-CoV-2 and host RNA editing mechanisms might have partially resulted in the bias in nucleotide changes during SARS-CoV-2 evolution. Finally, we outline some current challenges, including difficulty in deciphering the complicated relationship between viral pathogenicity and infectivity of different variants, and monitoring transmission of SARS-CoV-2 between humans and animals as the pandemic progresses.

Structural basis for continued antibody evasion by the SARS-CoV-2 receptor binding domain

Many studies have examined the impact of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) variants on neutralizing antibody activity after they have become dominant strains. Here, we evaluate the consequences of further viral evolution. We demonstrate mechanisms through which the SARS-CoV-2 receptor binding domain (RBD) can tolerate large numbers of simultaneous antibody escape mutations and show that pseudotypes containing up to seven mutations, as opposed to the one to three found in previously studied variants of concern, are more resistant to neutralization by therapeutic antibodies and serum from vaccine recipients. We identify an antibody that binds the RBD core to neutralize pseudotypes for all tested variants but show that the RBD can acquire an N-linked glycan to escape neutralization. Our findings portend continued emergence of escape variants as SARS-CoV-2 adapts to humans.

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.

Fine-tuning the spike: role of the nature and topology of the glycan shield in the structure and dynamics of the SARS-CoV-2 S

The dense glycan shield is an essential feature of the SARS-CoV-2 spike (S) architecture, key to immune evasion and to the activation of the prefusion conformation. Recent studies indicate that the occupancy and structures of the SARS-CoV-2 S glycans depend not only on the nature of the host cell, but also on the structural stability of the trimer; a point that raises important questions about the relative competence of different glycoforms. Moreover, the functional role of the glycan shield in the SARS-CoV-2 pathogenesis suggests that the evolution of the sites of glycosylation is potentially intertwined with the evolution of the protein sequence to affect optimal activity. Our results from multi-microsecond molecular dynamics simulations indicate that the type of glycosylation at N234, N165 and N343 greatly affects the stability of the receptor binding domain (RBD) open conformation, and thus its exposure and accessibility. Furthermore, our results suggest that the loss of glycosylation at N370, a newly acquired modification in the SARS-CoV-2 S glycan shield's topology, may have contributed to increase the SARS-CoV-2 infectivity as we find that N-glycosylation at N370 stabilizes the closed RBD conformation by binding a specific cleft on the RBD surface. We discuss how the absence of the N370 glycan in the SARS-CoV-2 S frees the RBD glycan binding cleft, which becomes available to bind cell-surface glycans, and potentially increases host cell surface localization.

Reinventing positive-strand RNA virus reverse genetics

Reverse genetics is the prospective analysis of how genotype determines phenotype. In a typical experiment, a researcher alters a viral genome, then observes the phenotypic outcome. Among RNA viruses, this approach was first applied to positive-strand RNA viruses in the mid-1970s and over nearly 50 years has become a powerful and widely used approach for dissecting the mechanisms of viral replication and pathogenesis. During this time the global health importance of two virus groups, flaviviruses (genus Flavivirus, family Flaviviridae) and betacoronaviruses (genus Betacoronavirus, subfamily Orthocoronavirinae, family Coronaviridae), have dramatically increased, yet these viruses have genomes that are technically challenging to manipulate. As a result, several new techniques have been developed to overcome these challenges. Here I briefly review key historical aspects of positive-strand RNA virus reverse genetics, describe some recent reverse genetic innovations, particularly as applied to flaviviruses and coronaviruses, and discuss their benefits and limitations within the larger context of rigorous genetic analysis.

A defective viral genome strategy elicits broad protective immunity against respiratory viruses

RNA viruses generate defective viral genomes (DVGs) that can interfere with replication of the parental wild-type virus. To examine their therapeutic potential, we created a DVG by deleting the capsid-coding region of poliovirus. Strikingly, intraperitoneal or intranasal administration of this genome, which we termed eTIP1, elicits an antiviral response, inhibits replication, and protects mice from several RNA viruses, including enteroviruses, influenza, and SARS-CoV-2. While eTIP1 replication following intranasal administration is limited to the nasal cavity, its antiviral action extends non-cell-autonomously to the lungs. eTIP1 broad-spectrum antiviral effects are mediated by both local and distal type I interferon responses. Importantly, while a single eTIP1 dose protects animals from SARS-CoV-2 infection, it also stimulates production of SARS-CoV-2 neutralizing antibodies that afford long-lasting protection from SARS-CoV-2 reinfection. Thus, eTIP1 is a safe and effective broad-spectrum antiviral generating short- and long-term protection against SARS-CoV-2 and other respiratory infections in animal models.

Expression pattern and function of SARS-CoV-2 receptor ACE2

Since the outbreak at the end of 2019, SARS-CoV-2 has been spreading around the world for more than one year. Scientists have been intensely conducting research on this newly emerged coronavirus and the disease caused by it. Angiotensin-converting enzyme 2 (ACE2), as a receptor mediating the cellular entry of SARS-CoV-2, has become a hot spot for researchers. Here, we summarized the recent progresses on the function, expression and distribution characteristics of ACE2 in human body and among populations. We further discussed the interaction mechanism of ACE2 and SARS-CoV-2 S protein, focusing on key residues that effect interaction and binding ability of SARS-CoV-2 variants. This will facilitate researchers to better understand SARS-CoV-2 infection and transmission route, adaptation mechanism, and designing treatment strategies.

Delta variant of COVID-19: A simple explanation

Severe acute respiratory syndrome coronavirus 2, the virus that causes coronavirus disease (COVID-19), has undergone numerous mutations since its initial identification, leading to challenges in controlling the pandemic. Till date, several variants of concern have been identified. However, currently, the Delta variant (B.1.617.2) is the most dreaded one owing to its enhanced transmissibility and increased virulence. In addition, this variant can potentially facilitate fusion of the spike protein to cells or inhibit antibodies from binding to it. In this commentary, we have simplified the complexity of the nomenclature of variants related to COVID-19, concentrating on the Delta variant including its transmissibility, response to vaccines, and prevention.

Emerging SARS-CoV-2 variants follow a historical pattern recorded in outgroups infecting non-human hosts

The ability to predict emerging variants of SARS-CoV-2 would be of enormous value, as it would enable proactive design of vaccines in advance of such emergence. We estimated diversity of each site on a multiple sequence alignment (MSA) of the Spike (S) proteins from close relatives of SARS-CoV-2 that infected bat and pangolin before the pandemic. Then we compared the locations of high diversity sites in this MSA and those of mutations found in multiple emerging lineages of human-infecting SARS-CoV-2. This comparison revealed a significant correspondence, which suggests that a limited number of sites in this protein are repeatedly substituted in different lineages of this group of viruses. It follows, therefore, that the sites of future emerging mutations in SARS-CoV-2 can be predicted by analyzing their relatives (outgroups) that have infected non-human hosts. We discuss a possible evolutionary basis for these substitutions and provide a list of frequently substituted sites that potentially include future emerging variants in SARS-CoV-2.

Why are some coronavirus variants more infectious?

Since the start of the pandemic, SARS-CoV-2 has infected almost 200 million human hosts and is set to encounter and gain entry in many more in the coming months. As the coronavirus flourish, the evolutionary pressure selects those variants that can complete the infection cycle faster and reproduce in large numbers compared to others. This increase in infectivity and transmissibility coupled with the immune response from high viral load may cause moderate to severe disease. Whether this leads to enhanced virulence in the prevalent Alpha and Delta variants is still not clear. This review describes the different types of SARS-CoV-2 variants that are now prevalent, their emergence, the mutations responsible for their growth advantages, and how they affect vaccine efficacy and increase chances of reinfection. Finally, we have also summarized the efforts made to recognize and predict the mutations, which can cause immune escape and track their emergence through impactful genomic surveillance.