Natural selection in the evolution of SARS-CoV-2 in bats created a generalist virus and highly capable human pathogen

Phylodynamics of a regional SARS-CoV-2 rapid spreading event in Colorado in late 2020

Since the initial reported discovery of SARS-CoV-2 in late 2019, genomic surveillance has been an important tool to understand its transmission and evolution. Here, we sought to describe the underlying regional phylodynamics before and during a rapid spreading event that was documented by surveillance protocols of the United States Air Force Academy (USAFA) in late October-November of 2020. We used replicate long-read sequencing on Colorado SARS-CoV-2 genomes collected July through November 2020 at the University of Colorado Anschutz Medical campus in Aurora and the United States Air Force Academy in Colorado Springs. Replicate sequencing allowed rigorous validation of variation and placement in a phylogenetic relatedness network. We focus on describing the phylodynamics of a lineage that likely originated in the local Colorado Springs community and expanded rapidly over the course of two months in an outbreak within the well-controlled environment of the United States Air Force Academy. Divergence estimates from sampling dates indicate that the SARS-CoV-2 lineage associated with this rapid expansion event originated in late October 2020. These results are in agreement with transmission pathways inferred by the United States Air Force Academy, and provide a window into the evolutionary process and transmission dynamics of a potentially dangerous but ultimately contained variant.

SARS-COV-2/COVID-19: scenario, epidemiology, adaptive mutations, and environmental factors

The coronavirus pandemic of 2019 has already exerted an enormous impact. For over a year, the worldwide pandemic has ravaged the whole globe, with approximately 250 million verified human infection cases and a mortality rate surpassing 4 million. While the genetic makeup of the related pathogen (SARS-CoV-2) was identified, many unknown facets remain a mystery, comprising the virus's origin and evolutionary trend. There were many rumors that SARS-CoV-2 was human-borne and its evolution was predicted many years ago, but scientific investigation proved them wrong and concluded that bats might be the origin of SARS-CoV-2 and pangolins act as intermediary species to transmit the virus from bats to humans. Airborne droplets were found to be the leading cause of human-to-human transmission of this virus, but later studies showed that contaminated surfaces and other environmental factors are also involved in its transmission. The evolution of different SARS-CoV-2 variants worsens the condition and has become a challenge to overcome this pandemic. The emergence of COVID-19 is still a mystery, and scientists are unable to explain the exact origin of SARS-CoV-2. This review sheds light on the possible origin of SARS-CoV-2, its transmission, and the key factors that worsen the situation.

Role and Limits of COVID-19 Vaccines in the Delicate Transition from Pandemic Mitigation to Endemic Control

The recent surge of COVID-19 related to the Omicron variant emergence has thrown a harsh light upon epidemic control in the near future. This should lead the scientific and medical community to question the long-term vaccine strategy for SARS-CoV-2 control. We provide here a critical point of view regarding the virological evolution, epidemiological aspects, and immunological drivers for COVID-19 control, including a vaccination strategy. Overall, we need more innovations in vaccine development to reduce the COVID-19 burden long term. The most adequate answer might be better cooperation between universities, biotech and pharmaceutical companies

The Genetic Characterization of the First Detected Bat Coronaviruses in Poland Revealed SARS-Related Types and Alphacoronaviruses

Bats are a major global reservoir of alphacoronaviruses (alphaCoVs) and betaCoVs. Attempts to discover the causative agents of COVID-19 and SARS have revealed horseshoe bats (Rhinolophidae) to be the most probable source of the virus. We report the first detection of bat coronaviruses (BtCoVs) in insectivorous bats in Poland and highlight SARS-related coronaviruses found in Rhinolophidae bats. The study included 503 (397 oral swabs and 106 fecal) samples collected from 20 bat species. Genetically diverse BtCoVs (n = 20) of the Alpha- and Betacoronavirus genera were found in fecal samples of two bat species. SARS-related CoVs were in 18 out of 58 lesser horseshoe bat (Rhinolophus hipposideros) samples (31%, 95% CI 20.6–43.8), and alphaCoVs were in 2 out of 55 Daubenton’s bat (Myotis daubentonii) samples (3.6%, 95% CI 0.6–12.3). The overall BtCoV prevalence was 4.0% (95% CI 2.6–6.1). High identity was determined for BtCoVs isolated from European M. daubentonii and R. hipposideros bats. The detection of SARS-related and alphaCoVs in Polish bats with high phylogenetic relatedness to reference BtCoVs isolated in different European countries but from the same species confirms their high host restriction. Our data elucidate the molecular epidemiology, prevalence, and geographic distribution of coronaviruses and particularly SARS-related types in the bat population.

SARS-CoV-2 origins: zoonotic Rhinolophus vs contemporary models

The question of the origin of coronavirus spread like wildfire ever since it wreaked havoc among humankind, and ever since the scientific community has worked tirelessly to trace the history of the virus. In this review, we have tried to compile relevant literature pertaining to the different theories of origin of SARS-CoV-2, hopefully without any bias, and we strongly support the zoonotic origin of the infamous SARS-CoV-2 in bats and its transfer to human beings through the most probable evolutionary hosts, pangolins and minks. We also support the contemporary 'Circulation Model' that simply mirrors the concept of evolution to explain the origin of the virus which, the authors believe, is the most rational school of thought. The most recent variant of SARS-CoV-2, Omicron, has been taken as an example to clarify the concept. We recommend the community to refer to this model for further understanding and delving deep into this mystery of the origin of SARS-CoV-2.

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.

Structural aspects of SARS-CoV-2 mutations: Implications to plausible infectivity with ACE-2 using computational modeling approach

Some of the SARS-CoV-2 variants are said to be more infectious than the previous others and are causing panic around the globe. Cases related to Delta plus (δ+) and omicron (ο) variants are on the rise worldwide. This sudden surge warrants an investigation into the reasons for its binding with ACE-2. The present study attempts to find out the structural basis of binding interactions of SARS-CoV-2 mutants based on computational modeling and comparative analysis. In silico strategies including protein-protein docking, mutation analysis, molecular dynamics, and binding energy calculations were used to study the binding of the 'receptor binding domain' (RBD) of the seven 'variants of concern' which include Alpha (α), Beta (β), Gamma (γ), Kappa (κ), Delta (δ), Delta plus (δ+) and omicron (ο) with ACE-2 (human angiotensin-converting enzyme-2) and with antibodies. Among all the variants dealt with in this study, Delta plus and omicron were found to be binding more strongly to ACE-2 than others due to inherent mutations and the consequent change in the hydrophilic and hydrophobic environment of the binding site. Furthermore, molecular dynamic (MD) simulations and subsequent MM/PBSA calculations provided useful structural insights into key residues participating in the interaction. Infectivity of a virus could be dependent on the interplay of evading antibodies and simultaneously attaching strongly with the host receptor. A cross-correlation between mutant spike proteins' binding with ACE-2 and antibodies provides a holistic assessment of the binding nature of these mutants vis-à-vis native virus and offers opportunities for designing potential therapeutics against these new mutants.

SARS-CoV-2 and the Missing Link of Intermediate Hosts in Viral Emergence - What We Can Learn From Other Betacoronaviruses

The emergence of SARS-CoV-2 in 2019 has resulted in a global pandemic with devastating human health and economic consequences. The development of multiple vaccines, antivirals and supportive care modalities have aided in our efforts to gain control of the pandemic. However, the emergence of multiple variants of concern and spillover into numerous nonhuman animal species could protract the pandemic. Further, these events also increase the difficulty in simultaneously monitoring viral evolution across multiple species and predicting future spillback potential into the human population. Here, we provide historic context regarding the roles of reservoir and intermediate hosts in coronavirus circulation and discuss current knowledge of these for SARS-CoV-2. Increased understanding of SARS-CoV-2 zoonoses are fundamental for efforts to control the global health and economic impacts of COVID-19.

Evidence for yeast artificial synthesis in SARS-CoV-2 and SARS-CoV-1 genomic sequences

Background: Knowledge about the origin of SARS-CoV-2 is necessary for both a biological and epidemiological understanding of the COVID-19 pandemic. Evidence suggests that a proximal evolutionary ancestor of SARS-CoV-2 belongs to the bat coronavirus family. However, as further evidence for a direct zoonosis remains limited, alternative modes of SARS-CoV-2 biogenesis should also be considered.    Results: Here we show that the genomes of SARS-CoV-2 and SARS-CoV-1 significantly diverge from other SARS-like coronaviruses through short chromosomal sequences from the yeast S. cerevisiae at focal positions that are known to be critical for host cell invasion, virus replication, and host immune response. For SARS-CoV-1, we identify two sites: one at the start of the RNA dependent RNA polymerase gene, and the other at the start of the spike protein’s receptor binding domain; for SARS-CoV-2, one at the start of the viral replicase domain, and the other toward the end of the spike gene past its domain junction. At this junction, we detect a highly specific stretch of yeast origin covering the critical furin cleavage site insert PRRA, which has not been seen in other lineage b betacoronaviruses. As yeast is not a natural host for this virus family, we propose an artificial synthesis model for viral constructs in yeast cells based on co-transformation of virus DNA plasmids carrying yeast selectable genetic markers followed by intra-chromosomal homologous recombination through gene conversion. Highly differential yeast sequence patterns congruent with chromosomes harboring specific auxotrophic markers further support yeast artificial synthesis. Conclusions: These results provide evidence that the genomes of SARS-CoV-1 and SARS-CoV-2 contain sequence information that points to their artificial synthesis in genetically modified yeast cells. Our data specifically allow the identification of the yeast S. cerevisiae as a potential recombination donor for the critical furin cleavage site in SARS-CoV-2.

The Impact of the Highly Virulent SARS-CoV-2 Gamma Variant on Young Adults in the State of São Paulo: Was It Inevitable?

Background

The coronavirus disease 2019 (COVID-19) pandemic had and is still having a tremendous impact on people all over the world, but it has been particularly harsh in South America. Nine out of 13 South American countries are among the 50 countries with the highest COVID-19 death rates. The gamma severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) variant that emerged by the end of 2020 in the Brazilian Amazon quickly spread throughout the country causing the harsh COVID-19 second wave. This variant displayed high viral loads, high transmissibility, and increased virulence as compared to previous variants.

Aims

The aim of this retrospective study is to revisit and analyse the epidemiology of the COVID-19 second wave in the state of São Paulo, the most populous Brazilian state. In addition to examining the possible factors that led to the emergence and propagation of the gamma variant, measures that could have prevented its spread and that of other highly virulent variants were also investigated.

Materials and methods

Data from São Paulo's official sources on morbidity, mortality, age distribution, and testing prior to and during the COVID-19 second wave (February - June 2021) and data regarding the distribution of SARS-CoV-2 variants in the country were parsed, analyzed, and compared to the period that anteceded the eruption of the second COVID-19 wave.

Results

In the state of São Paulo, the toll of the COVID-19 second wave surpassed that of the first 11 months of the pandemic (from March 2020 to January 2021), as 56% of the deaths occurred in the five months of the second wave between February and June 2021. The mean age of COVID-19 victims, which was already below life expectancy in the state dropped even further in the pandemic's second wave, reaching an average of 60 years of age. The years of life lost per death per month doubled and the case-fatality rate (CFR) of young adults (20-39 years old) more than trebled during this period. A number of hypotheses have been raised that might explain the emergence and spread of the gamma variant and the measures that could have been taken to prevent it and minimise its impact on the population.

Conclusions

Over 142,000 people died as a result of the SARS-CoV-2 gamma variant sweep in São Paulo in the first semester of 2021. Due to its high viral load, the gamma variant displayed high transmissibility and a high degree of virulence resulting in increased case fatality rates across most age tiers. Notably, this second wave was marked by a very significant increase in deaths among young adults. This increase was at least partially due to a deterioration in general health provoked by non-pharmaceutical interventions. In hindsight, a safer and more effective measure might have been to allow the free spread of the virus among the young and healthy in the first wave, thus conferring immunity against more virulent variants that emerged later on.

The structural role of SARS-CoV-2 genetic background in the emergence and success of spike mutations: The case of the spike A222V mutation

The S:A222V point mutation, within the G clade, was characteristic of the 20E (EU1) SARS-CoV-2 variant identified in Spain in early summer 2020. This mutation has since reappeared in the Delta subvariant AY.4.2, raising questions about its specific effect on viral infection. We report combined serological, functional, structural and computational studies characterizing the impact of this mutation. Our results reveal that S:A222V promotes an increased RBD opening and slightly increases ACE2 binding as compared to the parent S:D614G clade. Finally, S:A222V does not reduce sera neutralization capacity, suggesting it does not affect vaccine effectiveness.

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.

Drivers of adaptive evolution during chronic SARS-CoV-2 infections

In some immunocompromised patients with chronic severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection, considerable adaptive evolution occurs. Some substitutions found in chronic infections are lineage-defining mutations in variants of concern (VOCs), which has led to the hypothesis that VOCs emerged from chronic infections. In this study, we searched for drivers of VOC-like emergence by consolidating sequencing results from a set of 27 chronic infections. Most substitutions in this set reflected lineage-defining VOC mutations; however, a subset of mutations associated with successful global transmission was absent from chronic infections. We further tested the ability to associate antibody evasion mutations with patient-specific and virus-specific features and found that viral rebound is strongly correlated with the emergence of antibody evasion. We found evidence for dynamic polymorphic viral populations in most patients, suggesting that a compromised immune system selects for antibody evasion in particular niches in a patient’s body. We suggest that a tradeoff exists between antibody evasion and transmissibility and that extensive monitoring of chronic infections is necessary to further understanding of VOC emergence.

Analysis of mutations that arise in chronic SARS-CoV-2 infections shows both overlap and differences with mutations present in pandemic viral variants of concern, highlighting their distinct drivers of evolution.

The past, current and future epidemiological dynamic of SARS-CoV-2

SARS-CoV-2, the agent of the COVID-19 pandemic, emerged in late 2019 in China, and rapidly spread throughout the world to reach all continents. As the virus expanded in its novel human host, viral lineages diversified through the accumulation of around two mutations a month on average. Different viral lineages have replaced each other since the start of the pandemic, with the most successful Alpha, Delta and Omicron variants of concern (VoCs) sequentially sweeping through the world to reach high global prevalence. Neither Alpha nor Delta was characterized by strong immune escape, with their success coming mainly from their higher transmissibility. Omicron is far more prone to immune evasion and spread primarily due to its increased ability to (re-)infect hosts with prior immunity. As host immunity reaches high levels globally through vaccination and prior infection, the epidemic is expected to transition from a pandemic regime to an endemic one where seasonality and waning host immunization are anticipated to become the primary forces shaping future SARS-CoV-2 lineage dynamics. In this review, we consider a body of evidence on the origins, host tropism, epidemiology, genomic and immunogenetic evolution of SARS-CoV-2 including an assessment of other coronaviruses infecting humans. Considering what is known so far, we conclude by delineating scenarios for the future dynamic of SARS-CoV-2, ranging from the good-circulation of a fifth endemic 'common cold' coronavirus of potentially low virulence, the bad-a situation roughly comparable with seasonal flu, and the ugly-extensive diversification into serotypes with long-term high-level endemicity.

Le virus SARS-CoV-2 n’a pas « d’origine »

Depuis le début de la pandémie de COVID-19, la question de l’origine de ce virus fait l’objet d’une vive polémique. C’est la question de « l’origine » qui est biaisée. Darwin a montré qu’il n’y a pas d’origine déterminée à aucune espèce animale ou végétale, simplement un processus évolutif et sélectif. Il en est de même pour les virus, il n’y a pas d’origine, mais un processus évolutif. Les virus circulent d’hôte à hôte, animaux ou humains. Les virus pandémiques circulent déjà chez l’homme et évoluent avant l’apparition d’une maladie. Ce processus évolutif se poursuit et donne naissance à des variants successifs. La solution n’est pas de cibler la maladie ou le possible agent causal, mais plutôt de cibler le processus d’émergence de la maladie lui-même.

A Paradigm Shift in the Combination Changes of SARS-CoV-2 Variants and Increased Spread of Delta Variant (B.1.617.2) across the World

Since September 2020, the SARS-CoV-2 variants have gained their dominance worldwide, especially in Kenya, Italy, France, the UK, Turkey, Indonesia, India, Finland, Ireland, Singapore, Denmark, Germany, and Portugal. In this study, we developed a model on the frequency of delta variants across 28 countries (R²= 0.1497), displaying the inheritance of mutations during the generation of the delta variants with 123,526 haplotypes. The country-wise haplotype network showed the distribution of haplotypes in USA (10,174), Denmark (5,637), India (4,089), Germany (2,350), Netherlands (1,899), Sweden (1,791), Italy (1,720), France (1,293), Ireland (1,257), Belgium (1,207), Singapore (1,193), Portugal (1,184) and Spain (1,133). Our analysis shows the highest haplotype in Europe with 84% and the lowest in Australia with 0.00001%. A model of scatter plot was generated with a regression line which provided the estimated rate of mutation, including 24.048 substitutions yearly. Our study concluded that the high global prevalence of the delta variants is due to a high frequency of infectivity, supporting the paradigm shift of the viral variants.

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.

Transmission of SARS-CoV-2 from humans to animals and potential host adaptation

SARS-CoV-2, the causative agent of the COVID-19 pandemic, can infect a wide range of mammals. Since its spread in humans, secondary host jumps of SARS-CoV-2 from humans to multiple domestic and wild populations of mammals have been documented. Understanding the extent of adaptation to these animal hosts is critical for assessing the threat that the spillback of animal-adapted SARS-CoV-2 into humans poses. We compare the genomic landscapes of SARS-CoV-2 isolated from animal species to that in humans, profiling the mutational biases indicative of potentially different selective pressures in animals. We focus on viral genomes isolated from mink (Neovison vison) and white-tailed deer (Odocoileus virginianus) for which multiple independent outbreaks driven by onward animal-to-animal transmission have been reported. We identify five candidate mutations for animal-specific adaptation in mink (NSP9_G37E, Spike_F486L, Spike_N501T, Spike_Y453F, ORF3a_L219V), and one in deer (NSP3a_L1035F), though they appear to confer a minimal advantage for human-to-human transmission. No considerable changes to the mutation rate or evolutionary trajectory of SARS-CoV-2 has resulted from circulation in mink and deer thus far. Our findings suggest that minimal adaptation was required for onward transmission in mink and deer following human-to-animal spillover, highlighting the 'generalist' nature of SARS-CoV-2 as a mammalian pathogen.

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.

Molecular evolution and structural analyses of the spike glycoprotein from Brazilian SARS-CoV-2 genomes: the impact of selected mutations

The COVID-19 pandemic caused by SARS-CoV-2 has reached by February 2022 more than 380 million cases and 5.5 million deaths worldwide since its beginning in late 2019, leading to enhanced concern in the scientific community and the general population. One of the most important pieces of this host-pathogen interaction is the spike protein, which binds to the hACE2 cell receptor, mediates the membrane fusion and is the major target of neutralizing antibodies against SARS-CoV-2. The multiple amino acid substitutions observed in this region, specially in RBD have enhanced the hACE2 binding affinity and led to several modifications in the mechanisms of SARS-CoV-2 pathogenesis, improving the viral fitness and/or promoting immune evasion, with potential impact in the vaccine development. In this work, we identified 48 sites under selective pressures, 17 of them with the strongest evidence by the HyPhy tests, including VOC related mutation sites 138, 142, 222, 262, 484, 681, and 845, among others. The coevolutionary analysis identified 28 sites found not to be conditionally independent, such as E484K-N501Y. The molecular dynamics and free energy estimates showed the structural stabilizing effect and the higher impact of E484K for enhanced binding affinity between the spike RBD and hACE2 in P.1 and P.2 lineages (specially with L452V). Structural changes were also identified in the hACE molecule when interacting with B.1.1.7 RDB. Despite some destabilizing substitutions, a stabilizing effect was identified for the majority of the positively selected mutations.Communicated by Ramaswamy H. Sarma.

There is no "origin" to SARS-CoV-2

Since the beginning of the COVID-19 pandemic in 2020 caused by SARS-CoV-2, the question of the origin of this virus has been a highly debated issue. Debates have been, and are still, very disputed and often violent between the two main hypotheses: a natural origin through the "spillover" model or a laboratory-leak origin. Tenants of these two options are building arguments often based on the discrepancies of the other theory. The main problem is that it is the initial question of the origin itself which is biased. Charles Darwin demonstrated in 1859 that all species are appearing through a process of evolution, adaptation and selection. There is no determined origin to any animal or plant species, simply an evolutionary and selective process in which chance and environment play a key role. The very same is true for viruses. There is no determined origin to viruses, simply also an evolutionary and selective process in which chance and environment play a key role. However, in the case of viruses the process is slightly more complex because the "environment" is another living organism. Pandemic viruses already circulate in humans prior to the emergence of a disease. They are simply not capable of triggering an epidemic yet. They must evolve in-host, i.e. in-humans, for that. The evolutionary process which gave rise to SARS-CoV-2 is still ongoing with regular emergence of novel variants more adapted than the previous ones. The real relevant question is how these viruses can emerge as pandemic viruses and what the society can do to prevent the future emergence of pandemic viruses.

The Impact of Evolving SARS-CoV-2 Mutations and Variants on COVID-19 Vaccines

The emergence of several new variants of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) in recent months has raised concerns around the potential impact on ongoing vaccination programs. Data from clinical trials and real-world evidence suggest that current vaccines remain highly effective against the alpha variant (B.1.1.7), while some vaccines have reduced efficacy and effectiveness against symptomatic disease caused by the beta variant (B.1.351) and the delta variant (B.1.617.2); however, effectiveness against severe disease and hospitalization caused by delta remains high. Although data on the effectiveness of the primary regimen against omicron (B.1.1.529) are limited, booster programs using mRNA vaccines have been shown to restore protection against infection and symptomatic disease (regardless of the vaccine used for the primary regimen) and maintain high effectiveness against hospitalization. However, effectiveness against infection and symptomatic disease wanes with time after the booster dose. Studies have demonstrated reductions of varying magnitude in neutralizing activity of vaccine-elicited antibodies against a range of SARS-CoV-2 variants, with the omicron variant in particular exhibiting partial immune escape. However, evidence suggests that T-cell responses are preserved across vaccine platforms, regardless of variant of concern. Nevertheless, various mitigation strategies are under investigation to address the potential for reduced efficacy or effectiveness against current and future SARS-CoV-2 variants, including modification of vaccines for certain variants (including omicron), multivalent vaccine formulations, and different delivery mechanisms.

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.

Translational Control of COVID-19 and Its Therapeutic Implication

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is the causative agent of COVID-19, which has broken out worldwide for more than two years. However, due to limited treatment, new cases of infection are still rising. Therefore, there is an urgent need to understand the basic molecular biology of SARS-CoV-2 to control this virus. SARS-CoV-2 replication and spread depend on the recruitment of host ribosomes to translate viral messenger RNA (mRNA). To ensure the translation of their own mRNAs, the SARS-CoV-2 has developed multiple strategies to globally inhibit the translation of host mRNAs and block the cellular innate immune response. This review provides a comprehensive picture of recent advancements in our understanding of the molecular basis and complexity of SARS-CoV-2 protein translation. Specifically, we summarize how this viral infection inhibits host mRNA translation to better utilize translation elements for translation of its own mRNA. Finally, we discuss the potential of translational components as targets for therapeutic interventions.

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.

SARS-CoV-2: tracing the origin, tracking the evolution

The origin of SARS-CoV-2 is uncertain. Findings support a "bat origin" but results are not highly convincing. Studies found evidence that SARS-CoV-2 was around for many years before the pandemic outbreak. Evidence has been published that the progenitor of SARS-CoV-2 already had the capability to bind strongly to the human ACE2 receptor. This may be an indication that many other animal viruses are capable to jump to humans, having already affinity for a human receptor. This is quite worrying since current ecosystems' collapse brings people to high proximity with animals, increasing probabilities for random viral transitions. On the other hand, future adaptation of SARS-CoV-2 is of great concern. Virus-host interactions are complicated and unfortunately, we still do not have accurate tools for predicting viruses' future evolution. Viral adaptation is a multifactorial process and probably SARS-CoV-2 will not become soon, as we wish, a harmless infection. However, humanity is currently under the largest vaccination program and it's of great interest to see if vaccinations will change the evolutionary game against the virus.

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.

SARS-CoV-2 host prediction based on virus-host genetic features

The genetic diversity of the Coronaviruses gives them different biological abilities, such as infect different cells and/or organisms, a wide spectrum of clinical manifestations, their different routes of dispersion, and viral transmission in a specific host. In recent decades, different Coronaviruses have emerged that are highly adapted for humans and causing serious diseases, leaving their host of unknown origin. The viral genome information is particularly important to enable the recognition of patterns linked to their biological characteristics, such as the specificity in the host-parasite relationship. Here, based on a previously computational tool, the Seq2Hosts, we developed a novel approach which uses new variables obtained from the frequency of spike-Coronaviruses codons, the Relative Synonymous Codon Usage (RSCU) to shed new light on the molecular mechanisms involved in the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) host specificity. By using the RSCU obtained from nucleotide sequences before the SARS-CoV-2 pandemic, we assessed the possibility of know the hosts capable to be infected by these new emerging species, which was first identified infecting humans during 2019 in Wuhan, China. According to the model trained and validated using sequences available before the pandemic, bats are the most likely the natural host to the SARS-CoV-2 infection, as previously suggested in other studies that searched for the host viral origin.

Evidence for yeast artificial synthesis in SARS-CoV-2 and SARS-CoV-1 genomic sequences

Background: Knowledge about the origin of SARS-CoV-2 is necessary for both a biological and epidemiological understanding of the COVID-19 pandemic. Evidence suggests that a proximal evolutionary ancestor of SARS-CoV-2 belongs to the bat coronavirus family. However, as further evidence for a direct zoonosis remains limited, alternative modes of SARS-CoV-2 biogenesis should be considered.    Results: Here we show that the genomes from SARS-CoV-2 and from SARS-CoV-1 are differentially enriched with short chromosomal sequences from the yeast S. cerevisiae at focal positions that are known to be critical for host cell invasion, virus replication, and host immune response. For SARS-CoV-1, we identify two sites: one at the start of the RNA dependent RNA polymerase gene, and the other at the start of the spike protein’s receptor binding domain; for SARS-CoV-2, one at the start of the viral replicase domain, and the other toward the end of the spike gene past its domain junction. At this junction, we detect a highly specific stretch of yeast DNA origin covering the critical furin cleavage site insert PRRA, which has not been seen in other lineage b betacoronaviruses. As yeast is not a natural host for this virus family, we propose a passage model for viral constructs in yeast cells based on co-transformation of virus DNA plasmids carrying yeast selectable genetic markers followed by intra-chromosomal homologous recombination through gene conversion. Highly differential sequence homology data across yeast chromosomes congruent with chromosomes harboring specific auxotrophic markers further support this passage model. Conclusions: These results provide evidence that among SARS-like coronaviruses only the genomes of SARS-CoV-1 and SARS-CoV-2 contain information that points to a synthetic passage in genetically modified yeast cells. Our data specifically allow the identification of the yeast S. cerevisiae as a potential recombination donor for the critical furin cleavage site in SARS-CoV-2.

Transmissibility and pathogenicity of SARS-CoV-2 variants in animal models

As of February 2022, SARS-CoV-2 is still one of the most serious public health threats due to its high mortality rate and rapid spread of novel variants. Since the first outbreak in 2019, general understanding of SARS-CoV-2 has been improved through basic and clinical studies; however, knowledge gaps still exist in our understanding of the emerging novel SARSCoV-2 variants, which impacts the corresponding development of vaccines and therapeutics. Especially, accumulation of mutations in SARS-CoV-2 and rapid spread in populations with previous immunity has resulted in selection of variants that evade the host immune response. This phenomenon threatens to render current SARS-CoV-2 vaccines ineffective for controlling the pandemic. Proper animal models are essential for detailed investigations into the viral etiology, transmission and pathogenesis mechanisms, as well as evaluation of the efficacy of vaccine candidates against recent SARS-CoV-2 variants. Further, the choice of animal model for each research topic is important for researchers to gain better knowledge of recent SARS-CoV-2 variants. Here, we review the advantages and limitations of each animal model, including mice, hamsters, ferrets, and non-human primates, to elucidate variant SARS-CoV-2 etiology and transmission and to evaluate therapeutic and vaccine efficacy.

Analysis of SARS-CoV-2 synonymous codon usage evolution throughout the COVID-19 pandemic

SARS-CoV-2, the seventh coronavirus known to infect humans, can cause severe life-threatening respiratory pathologies. To better understand SARS-CoV-2 evolution, genome-wide analyses have been made, including the general characterization of its codons usage profile. Here we present a bioinformatic analysis of the evolution of SARS-CoV-2 codon usage over time using complete genomes collected since December 2019. Our results show that SARS-CoV-2 codon usage pattern is antagonistic to, and it is getting farther away from that of the human host. Further, a selection of deoptimized codons over time, which was accompanied by a decrease in both the codon adaptation index and the effective number of codons, was observed. All together, these findings suggest that SARS-CoV-2 could be evolving, at least from the perspective of the synonymous codon usage, to become less pathogenic.

Evolutionary Shift from Purifying Selection towards Divergent Selection of SARS-CoV2 Favors its Invasion into Multiple Human Organs

SARS-CoV2 virus is believed to be originated from a closely related bat Coronavirus RaTG13 lineage and uses its key entry-point residues in S1 protein to attach with human ACE2 receptor. SARS-CoV2 could enter human from bat with its poorly developed entry-point residues much before its known appearance with slower mutation rate or recently with efficiently developed entry-point residues with higher mutation rate or through an intermediate host. Temporal analysis of SARS-CoV2 genome shows that its nucleotide substitution rate is as low as 27nt/year with an evolutionary rate of 9×10⁻⁴/site/year, which is well within the range of other RNA virus (10⁻⁴ to 10⁻⁶/site/year). TMRCA of SARS-CoV2 from bat RaTG13 lineage appears to be in between 9 and 14 years. Evolution of a critical entry-point residue Y493Q needs two substitutions with an intermediate virus carrying Y493H (Y>H>Q) but has not been identified in known twenty-nine bat CoV virus. Genetic codon analysis indicates that SARS-CoV2 evolution during propagation in human disobeys neutral evolution as nonsynonymous mutations surpass synonymous mutations with the increase of ω (d_n/d_s). Taken together, genetic data suggests that SARS-CoV2 is originated long time back before its appearance in human in 2019. Increase of ω signifies that SARs-CoV2 evolution is approaching towards diversifying selection from purifying selection predictably for its infection power to evade multiple human organs.

Exploring the Natural Origins of SARS-CoV-2 in the Light of Recombination

The lack of an identifiable intermediate host species for the proximal animal ancestor of SARS-CoV-2, and the large geographical distance between Wuhan and where the closest evolutionary related coronaviruses circulating in horseshoe bats (members of the Sarbecovirus subgenus) have been identified, is fueling speculation on the natural origins of SARS-CoV-2. We performed a comprehensive phylogenetic study on SARS-CoV-2 and all the related bat and pangolin sarbecoviruses sampled so far. Determining the likely recombination events reveals a highly reticulate evolutionary history within this group of coronaviruses. Distribution of the inferred recombination events is nonrandom with evidence that Spike, the main target for humoral immunity, is beside a recombination hotspot likely driving antigenic shift events in the ancestry of bat sarbecoviruses. Coupled with the geographic ranges of their hosts and the sampling locations, across southern China, and into Southeast Asia, we confirm that horseshoe bats, Rhinolophus, are the likely reservoir species for the SARS-CoV-2 progenitor. By tracing the recombinant sequence patterns, we conclude that there has been relatively recent geographic movement and cocirculation of these viruses' ancestors, extending across their bat host ranges in China and Southeast Asia over the last 100 years. We confirm that a direct proximal ancestor to SARS-CoV-2 has not yet been sampled, since the closest known relatives collected in Yunnan shared a common ancestor with SARS-CoV-2 approximately 40 years ago. Our analysis highlights the need for dramatically more wildlife sampling to: 1) pinpoint the exact origins of SARS-CoV-2's animal progenitor, 2) the intermediate species that facilitated transmission from bats to humans (if there is one), and 3) survey the extent of the diversity in the related sarbecoviruses' phylogeny that present high risk for future spillovers.

The slowing rate of CpG depletion in SARS-CoV-2 genomes is consistent with adaptations to the human host

Depletion of CpG dinucleotides in SARS-CoV-2 genomes has been linked to virus evolution, host-switching, virus replication, and innate immune responses. Temporal variations, if any, in the rate of CpG depletion during virus evolution in the host remain poorly understood. Here, we analysed the CpG content of over 1.4 million full-length SARS-CoV-2 genomes representing over 170 million documented infections during the first 17 months of the pandemic. Our findings suggest that the extent of CpG depletion in SARS-CoV-2 genomes is modest. Interestingly, the rate of CpG depletion is highest during early evolution in humans and it gradually tapers off almost reaching an equilibrium; this is consistent with adaptations to the human host. Furthermore, within the coding regions, CpG depletion occurs predominantly at codon positions 2-3 and 3-1. Loss of ZAP-binding motifs in SARS-CoV-2 genomes is primarily driven by the loss of the terminal CpG in the motifs. Nonetheless, majority of the CpG depletion in SARS-CoV-2 genomes occurs outside ZAP-binding motifs. SARS-CoV-2 genomes selectively lose CpGs-motifs from a U-rich context; this may help avoid immune recognition by TLR7. SARS-CoV-2 alpha-, beta- and delta-variants of concern have reduced CpG content compared to sequences from the beginning of the pandemic. In sum, we provide evidence that the rate of CpG depletion in virus genomes is not uniform and it greatly varies over time and during adaptations to the host. This work highlights how temporal variations in selection pressures during virus adaption may impact the rate and the extent of CpG depletion in virus genomes.

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.

SARS-CoV-2 competes with host mRNAs for efficient translation by maintaining the mutations favorable for translation initiation

During SARS-CoV-2 proliferation, the translation of viral RNAs is usually the rate-limiting step. Understanding the molecular details of this step is beneficial for uncovering the origin and evolution of SARS-CoV-2 and even for controlling the pandemic. To date, it is unclear how SARS-CoV-2 competes with host mRNAs for ribosome binding and efficient translation. We retrieved the coding sequences of all human genes and SARS-CoV-2 genes. We systematically profiled the GC content and folding energy of each CDS. Considering that some fixed or polymorphic mutations exist in SARS-CoV-2 and human genomes, all algorithms and analyses were applied to both pre-mutate and post-mutate versions. In SARS-CoV-2 but not human, the 5-prime end of CDS had lower GC content and less RNA structure than the 3-prime part, which was favorable for ribosome binding and efficient translation initiation. Globally, the fixed and polymorphic mutations in SARS-CoV-2 had created an even lower GC content at the 5-prime end of CDS. In contrast, no similar patterns were observed for the fixed and polymorphic mutations in human genome. Compared with human RNAs, the SARS-CoV-2 RNAs have less RNA structure in the 5-prime end and thus are more favorable of fast translation initiation. The fixed and polymorphic mutations in SARS-CoV-2 are further amplifying this advantage. This might serve as a strategy for SARS-CoV-2 to adapt to the human host.

Distinct codon usage bias evolutionary patterns between weakly and strongly virulent respiratory viruses

Human respiratory viruses are of vastly different virulence, giving rise to symptoms ranging from common cold to severe pneumonia or even death. Although this most likely impacts molecular evolution of the corresponding viruses, the specific differences in their evolutionary patterns remain largely unknown. By comparing structural and nonstructural genes within respiratory viruses, greater similarities in codon usage bias (CUB) between nonstructural genes and humans were observed in weakly virulent viruses, whereas in strongly virulent viruses, it was structural genes whose CUBs were more similar to that of humans. Further comparisons between genomes of weakly and strongly virulent coronaviruses revealed greater similarities in CUBs between strongly virulent viruses and humans. Finally, using phylogenetic independent contrasts, dissimilation of viral CUB from that of humans was observed in SARS-CoV-2. Our work revealed distinct CUB evolutionary patterns between weakly and strongly virulent viruses, a previously unrecognized interaction between CUB and virulence in respiratory viruses.

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.

The substitution spectra of coronavirus genomes

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) pandemic has triggered an unprecedented international effort to sequence complete viral genomes. We leveraged this wealth of information to characterize the substitution spectrum of SARS-CoV-2 and to compare it with those of other human and animal coronaviruses. We show that, once nucleotide composition is taken into account, human and most animal coronaviruses display a mutation spectrum dominated by C to U and G to U substitutions, a feature that is not shared by other positive-sense RNA viruses. However, the proportions of C to U and G to U substitutions tend to decrease as divergence increases, suggesting that, whatever their origin, a proportion of these changes is subsequently eliminated by purifying selection. Analysis of the sequence context of C to U substitutions showed little evidence of apolipoprotein B mRNA editing catalytic polypeptide-like (APOBEC)-mediated editing and such contexts were similar for SARS-CoV-2 and Middle East respiratory syndrome coronavirus sampled from different hosts, despite different repertoires of APOBEC3 proteins in distinct species. Conversely, we found evidence that C to U and G to U changes affect CpG dinucleotides at a frequency higher than expected. Whereas this suggests ongoing selective reduction of CpGs, this effect alone cannot account for the substitution spectra. Finally, we show that, during the first months of SARS-CoV-2 pandemic spread, the frequency of both G to U and C to U substitutions increased. Our data suggest that the substitution spectrum of SARS-CoV-2 is determined by an interplay of factors, including intrinsic biases of the replication process, avoidance of CpG dinucleotides and other constraints exerted by the new host.

Genomic Characterization of SARS-CoV-2 Isolated from Patients with Distinct Disease Outcomes in Mexico

The coronavirus disease 2019 (COVID-19), caused by the severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2), has shown a wide spectrum of clinical manifestations ranging from asymptomatic infections to severe disease and death. Pre-existing medical conditions and age have been mainly linked to the development of severe disease; however, the potential association of viral genetic characteristics with different clinical conditions remains unclear. SARS-CoV-2 variants with increased transmissibility were detected early in the pandemics, and several variants with potential relevance for public health are currently circulating around the world. In this study, we characterized 57 complete SARS-CoV-2 genomes during the exponential growth phase of the early epidemiological curve in Mexico, in April 2020. Patients were categorized under distinct disease severity outcomes: mild disease or ambulatory care, severe disease or hospitalized, and deceased. To reduce bias related to risk factors, the patients were less than 60 years old and with no diagnosed comorbidities A trait-association phylogenomic approach was used to explore genotype–phenotype associations, represented by the co-occurrence of mutations, disease severity outcome categories, and clusters of Mexican sequences. Phylogenetic results revealed a higher genomic diversity compared to the initial viruses detected during the early stage of the local epidemic. We identified a total of 90 single nucleotide variants compared to the Wuhan-Hu-1 genome, including 54 nonsynonymous mutations. We did not find evidence for the co-occurrence of mutations associated with specific disease outcomes. Therefore, in the group of patients studied, disease severity was likely mainly driven by the host genetic background and other demographic factors.

IMPORTANCE The genetic association of severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) with different clinical conditions remains unclear and needs further investigation. In this study, we characterized 57 complete SARS-CoV-2 genomes from patients in Mexico with distinct disease severity outcomes: mild disease or ambulatory care, severe disease or hospitalized, and deceased. To reduce bias related to risk factors the patients were less than 60 years old and with no diagnosed comorbidities. We did not find evidence for the co-occurrence of mutations associated with specific disease outcomes. Therefore, in the group of patients studied, disease severity was likely mainly driven by the host genetic background and other demographic factors.

The Remarkable Evolutionary Plasticity of Coronaviruses by Mutation and Recombination: Insights for the COVID-19 Pandemic and the Future Evolutionary Paths of SARS-CoV-2

Coronaviruses (CoVs) constitute a large and diverse subfamily of positive-sense single-stranded RNA viruses. They are found in many mammals and birds and have great importance for the health of humans and farm animals. The current SARS-CoV-2 pandemic, as well as many previous epidemics in humans that were of zoonotic origin, highlights the importance of studying the evolution of the entire CoV subfamily in order to understand how novel strains emerge and which molecular processes affect their adaptation, transmissibility, host/tissue tropism, and patho non-homologous genicity. In this review, we focus on studies over the last two years that reveal the impact of point mutations, insertions/deletions, and intratypic/intertypic homologous and non-homologous recombination events on the evolution of CoVs. We discuss whether the next generations of CoV vaccines should be directed against other CoV proteins in addition to or instead of spike. Based on the observed patterns of molecular evolution for the entire subfamily, we discuss five scenarios for the future evolutionary path of SARS-CoV-2 and the COVID-19 pandemic. Finally, within this evolutionary context, we discuss the recently emerged Omicron (B.1.1.529) VoC.

A history of coronaviruses

The history of coronaviruses is an account of the discovery of coronaviruses and the diseases they cause. It starts with a report of a new type of upper-respiratory tract disease among chickens in North Dakota, US, in 1931. The causative agent was identified as a virus in 1933. By 1936, the disease and the virus were recognised as unique from other viral diseases. The virus became known as infectious bronchitis virus (IBV), but later officially renamed as Avian coronavirus. A new brain disease of mice (murine encephalomyelitis) was discovered in 1947 at Harvard Medical School in Boston. The virus was called JHM (after Harvard pathologist John Howard Mueller). Three years later a new mouse hepatitis was reported from the National Institute for Medical Research in London. The causative virus was identified as mouse hepatitis virus (MHV), later renamed Murine coronavirus. In 1961, a virus was obtained from a school boy in Epsom, England, who was suffering from common cold. The sample, designated B814, was confirmed as novel virus in 1965. New common cold viruses (assigned 229E) collected from medical students at the University of Chicago were also reported in 1966. Structural analyses of IBV, MHV, B18 and 229E using transmission electron microscopy revealed that they all belong to the same group of viruses. Making a crucial comparison in 1967, June Almeida and David Tyrrell invented the collective name coronavirus, as all those viruses were characterised by solar corona-like projections (called spikes) on their surfaces. Other coronaviruses have been discovered from pigs, dogs, cats, rodents, cows, horses, camels, Beluga whales, birds and bats. As of 2022, 52 species are described. Bats are found to be the richest source of different species of coronaviruses. All coronaviruses originated from a common ancestor about 293 million years ago. Zoonotic species such as Severe acute respiratory syndrome-related coronavirus (SARS-CoV), Middle East respiratory syndrome-related coronavirus (MERS-CoV) and severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), a variant of SARS-CoV, emerged during the past two decades and caused the first pandemics of the 21st century.

Detection of SARS-CoV-2 in two cats during the second wave of the COVID-19 pandemic in France

Although there are several reports in the literature of SARS-CoV-2 infection in cats, few SARS-CoV-2 sequences from infected cats have been published. In this study, SARS-CoV-2 infection was evaluated in two cats by clinical observation, molecular biology (qPCR and NGS), and serology (microsphere immunoassay and seroneutralization). Following the observation of symptomatic SARS-CoV-2 infection in two cats, infection status was confirmed by RT-qPCR and, in one cat, serological analysis for antibodies against N-protein and S-protein, as well as neutralizing antibodies. Comparative analysis of five SARS-CoV-2 sequence fragments obtained from one of the cats showed that this infection was not with one of the three recently emerged variants of SARS-CoV-2. This study provides additional information on the clinical, molecular, and serological aspects of SARS-CoV-2 infection in cats.

An update on host immunity correlates and prospects of re-infection in COVID-19

Reinfection with SARS-CoV-2 is not frequent yet the incidence rate of it is increasing globally owing to the slow emergence of drift variants that pose a perpetual threat to vaccination strategies and have a greater propensity for disease reoccurrence. Long-term protection against SARS-CoV-2 reinfection relies on the induction of the innate as well as the adaptive immune response endowed with immune memory. However, a multitude of factors including the selection pressure, the waning immunity against SARS-CoV-2 over the first year after infection possibly favors evolution of more infectious immune escape variants, amplifying the risk of reinfection. Additionally, the correlates of immune protection, the novel SARS-CoV-2 variants of concern (VOC), the durability of the adaptive and mucosal immunity remain major challenges for the development of therapeutic and prophylactic interventions. Interestingly, a recent body of evidence indicated that the gastrointestinal (GI) tract is another important target organ for SARS-CoV-2 besides the respiratory system, potentially increasing the likelihood of reinfection by impacting the microbiome and the immune response via the gut-lung axis. In this review, we summarized the latest development in SARS-CoV-2 reinfection, and explored the untapped potential of trained immunity. We also highlighted the immune memory kinetics of the humoral and cell-mediated immune response, genetic drift of the emerging viral variants, and discussed the current challenges in vaccine development. Understanding the dynamics and the quality of immune response by unlocking the power of the innate, humoral and cell-mediated immunity during SARS-CoV-2 reinfection would open newer avenues for drug discovery and vaccine designs.

Molecular and Clinical Aspects of COVID-19 Vaccines and Other Therapeutic Interventions Apropos Emerging Variants of Concern

Coronavirus disease 2019 (COVID-19) has overwhelmed the healthcare and economy of the world, with emerging new variants of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) posing an everlasting threat to humanity. While most COVID-19 vaccines provide adequate protective immunological response against the original SARS-CoV-2 variant, there is a pressing need to understand their biological and clinical responses. Recent evidence suggests that some of the new variants of SARS-CoV-2 evade the protection conferred by the existing vaccines, which may impede the ongoing efforts to expedite the vaccination programs worldwide. These concerns have also highlighted the importance of a pan-COVID-19 vaccine, which is currently in the making. Thus, it is imperative to have a better molecular and clinical understanding of the various COVID-19 vaccines and their immunological trajectory against any emerging variant of concerns (VOCs) in particular to break this vicious cycle. Furthermore, other treatment regimens based on cellular therapies and monoclonal antibodies should be explored systematically as an alternative and readily available option considering the possibility of the emergence of more virulent SARS-CoV-2 mutants. In this review, we shed light on the various molecular mechanisms and clinical responses of COVID-19 vaccines. Importantly, we review the recent findings of their long-term immune protection and efficacy against emerging VOCs. Considering that other targeted and effective treatments will complement vaccine therapy, we provide a comprehensive understanding of the role of cell-based therapies, monoclonal antibodies, and immunomodulatory agents as alternative and readily available treatment modalities against any emerging SARS-CoV-2 variant.

D614G mutation eventuates in all VOI and VOC in SARS-CoV-2: Is it part of the positive selection pioneered by Darwin?

Recently, several emerging variants of SARS-CoV-2 have originated from the Wuhan strain and spread throughout the globe within one and a half years. One mutation, D614G, is very prominent in all VOI and VOC in SARS-CoV-2. This mutation might help to increase the viral fitness in all emerging variants where the mutation is present. With the help of this mutation (D614G), the SARS-CoV-2 variants have gained viral fitness to enhance viral replication and increase transmission. This paper attempts to answer the question of whether the mutation (D614G) occurs due to positive selection or not.

The emergence, genomic diversity and global spread of SARS-CoV-2

Since the first cases of COVID-19 were documented in Wuhan, China in 2019, the world has witnessed a devastating global pandemic, with more than 238 million cases, nearly 5 million fatalities and the daily number of people infected increasing rapidly. Here we describe the currently available data on the emergence of the SARS-CoV-2 virus, the causative agent of COVID-19, outline the early viral spread in Wuhan and its transmission patterns in China and across the rest of the world, and highlight how genomic surveillance, together with other data such as those on human mobility, has helped to trace the spread and genetic variation of the virus and has also comprised a key element for the control of the pandemic. We pay particular attention to characterizing and describing the international spread of the major variants of concern of SARS-CoV-2 that were first identified in late 2020 and demonstrate that virus evolution has entered a new phase. More broadly, we highlight our currently limited understanding of coronavirus diversity in nature, the rapid spread of the virus and its variants in such an increasingly connected world, the reduced protection of vaccines, and the urgent need for coordinated global surveillance using genomic techniques. In summary, we provide important information for the prevention and control of both the ongoing COVID-19 pandemic and any new diseases that will inevitably emerge in the human population in future generations.

The biological and clinical significance of emerging SARS-CoV-2 variants

The past several months have witnessed the emergence of SARS-CoV-2 variants with novel spike protein mutations that are influencing the epidemiological and clinical aspects of the COVID-19 pandemic. These variants can increase rates of virus transmission and/or increase the risk of reinfection and reduce the protection afforded by neutralizing monoclonal antibodies and vaccination. These variants can therefore enable SARS-CoV-2 to continue its spread in the face of rising population immunity while maintaining or increasing its replication fitness. The identification of four rapidly expanding virus lineages since December 2020, designated variants of concern, has ushered in a new stage of the pandemic. The four variants of concern, the Alpha variant (originally identified in the UK), the Beta variant (originally identified in South Africa), the Gamma variant (originally identified in Brazil) and the Delta variant (originally identified in India), share several mutations with one another as well as with an increasing number of other recently identified SARS-CoV-2 variants. Collectively, these SARS-CoV-2 variants complicate the COVID-19 research agenda and necessitate additional avenues of laboratory, epidemiological and clinical research.