Convergent Evolution in SARS-CoV-2 Spike Creates a Variant Soup from Which New COVID-19 Waves Emerge

SARS-CoV-2 evolution in the Omicron era

Since SARS-CoV-2 BA.5 (Omicron) emerged and spread in 2022, Omicron lineages have markedly diversified. Here we review the evolutionary trajectories and processes that underpin the emergence of these lineages, and identify the most prevalent sublineages. We discuss the potential origins of second-generation BA.2 lineages. Simple and complex recombination, antigenic drift and convergent evolution have enabled SARS-CoV-2 to accumulate mutations that alter its antigenicity. We also discuss the potential evolutionary trajectories of SARS-CoV-2 in the future.

Treatment-Emergent Cilgavimab Resistance Was Uncommon in Vaccinated Omicron BA.4/5 Outpatients

Mutations in the SARS-CoV-2 Spike glycoprotein can affect monoclonal antibody efficacy. Recent findings report the occurrence of resistant mutations in immunocompromised patients after tixagevimab/cilgavimab treatment. More recently, the Food and Drug Agency revoked the authorization for tixagevimab/cilgavimab, while this monoclonal antibody cocktail is currently recommended by the European Medical Agency. We retrospectively reviewed 22 immunocompetent patients at high risk for disease progression who received intramuscular tixagevimab/cilgavimab as early COVID-19 treatment and presented a prolonged high viral load. Complete SARS-CoV-2 genome sequences were obtained for a deep investigation of mutation frequencies in Spike protein before and during treatment. At seven days, only one patient showed evidence of treatment-emergent cilgavimab resistance. Quasispecies analysis revealed two different deletions on the Spike protein (S:del138-144 or S:del141-145) in combination with the resistance S:K444N mutation. The structural and dynamic impact of the two quasispecies was characterized by using molecular dynamics simulations, showing the conservation of the principal functional movements in the mutated systems and their capabilities to alter the structure and dynamics of the RBD, responsible for the interaction with the ACE2 human receptor. Our study underlines the importance of prompting an early virological investigation to prevent drug resistance or clinical failures in immunocompetent patients.

Impact of COVID-19 on Dental Practices in El Salvador and Mexico: A Comprehensive Survey Analysis

BACKGROUND AND OBJECTIVES

This study aimed to identify the relationship between prevention measures and protective barriers in dental practice in El Salvador and Mexico during the COVID-19 pandemic in 2020 and 2021.

MATERIALS AND METHODS

A longitudinal study was conducted from June 2020 to December 2021, involving 1,719 dentists divided into four groups based on location and year. A 20-question survey in Spanish was utilized and validated with a Cronbach's alpha value of 0.84.

RESULTS

The use of phone triage (OR = 1.3), thermometers (OR = 1.4), physical distancing (OR = 1.7), and face shields (OR = 2.6) was significantly associated with dental practice in both countries during the pandemic.

CONCLUSIONS

During 2020 and 2021, dental care in El Salvador and Mexico was significantly linked to COVID-19 preventive measures. Phone triage, thermometers, distancing, and face shields positively correlated with dental services. National health agencies should promote the use of minimum preventive measures in dental care, preparing for potential reinfections or new pandemics from emerging virus variants.

Unraveling the Dynamics of Omicron (BA.1, BA.2, and BA.5) Waves and Emergence of the Deltacton Variant: Genomic Epidemiology of the SARS-CoV-2 Epidemic in Cyprus (Oct 2021-Oct 2022)

Commencing in December 2019 with the emergence of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), three years of the coronavirus disease 2019 (COVID-19) pandemic have transpired. The virus has consistently demonstrated a tendency for evolutionary adaptation, resulting in mutations that impact both immune evasion and transmissibility. This ongoing process has led to successive waves of infections. This study offers a comprehensive assessment spanning genetic, phylogenetic, phylodynamic, and phylogeographic dimensions, focused on the trajectory of the SARS-CoV-2 epidemic in Cyprus. Based on a dataset comprising 4700 viral genomic sequences obtained from affected individuals between October 2021 and October 2022, our analysis is presented. Over this timeframe, a total of 167 distinct lineages and sublineages emerged, including variants such as Delta and Omicron (1, 2, and 5). Notably, during the fifth wave of infections, Omicron subvariants 1 and 2 gained prominence, followed by the ascendancy of Omicron 5 in the subsequent sixth wave. Additionally, during the fifth wave (December 2021-January 2022), a unique set of Delta sequences with genetic mutations associated with Omicron variant 1, dubbed "Deltacron", was identified. The emergence of this phenomenon initially evoked skepticism, characterized by concerns primarily centered around contamination or coinfection as plausible etiological contributors. These hypotheses were predominantly disseminated through unsubstantiated assertions within the realms of social and mass media, lacking concurrent scientific evidence to validate their claims. Nevertheless, the exhaustive molecular analyses presented in this study have demonstrated that such occurrences would likely lead to a frameshift mutation-a genetic aberration conspicuously absent in our provided sequences. This substantiates the accuracy of our initial assertion while refuting contamination or coinfection as potential etiologies. Comparable observations on a global scale dispelled doubt, eventually leading to the recognition of Delta-Omicron variants by the scientific community and their subsequent monitoring by the World Health Organization (WHO). As our investigation delved deeper into the intricate dynamics of the SARS-CoV-2 epidemic in Cyprus, a discernible pattern emerged, highlighting the major role of international connections in shaping the virus's local trajectory. Notably, the United States and the United Kingdom were the central conduits governing the entry and exit of the virus to and from Cyprus. Moreover, notable migratory routes included nations such as Greece, South Korea, France, Germany, Brazil, Spain, Australia, Denmark, Sweden, and Italy. These empirical findings underscore that the spread of SARS-CoV-2 within Cyprus was markedly influenced by the influx of new, highly transmissible variants, triggering successive waves of infection. This investigation elucidates the emergence of new waves of infection subsequent to the advent of highly contagious and transmissible viral variants, notably characterized by an abundance of mutations localized within the spike protein. Notably, this discovery decisively contradicts the hitherto hypothesis of seasonal fluctuations in the virus's epidemiological dynamics. This study emphasizes the importance of meticulously examining molecular genetics alongside virus migration patterns within a specific region. Past experiences also emphasize the substantial evolutionary potential of viruses such as SARS-CoV-2, underscoring the need for sustained vigilance. However, as the pandemic's dynamics continue to evolve, a balanced approach between caution and resilience becomes paramount. This ethos encourages an approach founded on informed prudence and self-preservation, guided by public health authorities, rather than enduring apprehension. Such an approach empowers societies to adapt and progress, fostering a poised confidence rooted in well-founded adaptation.

The diversification of SARS-CoV-2 Omicron variants and evaluation of their detection with molecular and rapid antigen assays

BACKGROUND

The SARS-CoV-2 pandemic saw the rapid rise, global spread, and diversification of the omicron variant in 2022. Given the overwhelming dominance of this variant globally and its diverse lineages, there is an urgent need to ensure that diagnostic assays are capable of detecting widely circulating omicron sub-lineages.

STUDY DESIGN

Remnant clinical VTM samples from SARS-CoV-2 PCR confirmed infections (n = 733) collected in Wisconsin (n = 94), New York (n = 267), and South Carolina (n = 372) throughout 2022 were sequenced, classified, and tested with m2000 RealTime SARS-CoV-2, Alinity m SARS-CoV-2, ID NOW COVID-19 v2.0, BinaxNOW COVID-19 Ag Card, and Panbio COVID-19 Rapid Test Device assays.

RESULTS

Sequences and lineage classifications were obtained for n = 641/733 (87.4%) samples and included delta (n = 6) and representatives from all major SARS-CoV-2 omicron variants circulating in 2022 (BA.1, BA.2, BA.3, BA.4, BA.5, BE, BF, BQ.1, and XBB). Panels of diverse omicron lineages were tested by molecular assays RealTime (n = 624), Alinity m (n = 80), and ID NOW v2.0 (n = 88) with results showing 100% detection for all samples. BinaxNOW and Panbio had sensitivities of 494/533 (92.7%) and 416/469 (88.7%), respectively for specimens with >4 log10 copies/test, consistent with expected performance for frozen specimens. Furthermore, BinaxNOW demonstrated SARS-CoV-2 detection in clinical samples 1-4 days, and up to 18 days post-symptom onset in BA.1 infected patients with >4 log10 copies/test.

CONCLUSIONS

This data highlights the rise and diversification of SARS-CoV-2 omicron variants over the course of 2022 and demonstrate that each of the 5 tested assays can detect the breadth of omicron variants circulating globally.

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

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

Markov State Models and Perturbation-Based Approaches Reveal Distinct Dynamic Signatures and Hidden Allosteric Pockets in the Emerging SARS-Cov-2 Spike Omicron Variant Complexes with the Host Receptor: The Interplay of Dynamics and Convergent Evolution Modulates Allostery and Functional Mechanisms

The new generation of SARS-CoV-2 Omicron variants displayed a significant growth advantage and increased viral fitness by acquiring convergent mutations, suggesting that the immune pressure can promote convergent evolution leading to the sudden acceleration of SARS-CoV-2 evolution. In the current study, we combined structural modeling, microsecond molecular dynamics simulations, and Markov state models to characterize conformational landscapes and identify specific dynamic signatures of the SARS-CoV-2 spike complexes with the host receptor ACE2 for the recently emerged highly transmissible XBB.1, XBB.1.5, BQ.1, and BQ.1.1 Omicron variants. Microsecond simulations and Markovian modeling provided a detailed characterization of the functional conformational states and revealed the increased thermodynamic stabilization of the XBB.1.5 subvariant, which can be contrasted to more dynamic BQ.1 and BQ.1.1 subvariants. Despite considerable structural similarities, Omicron mutations can induce unique dynamic signatures and specific distributions of the conformational states. The results suggested that variant-specific changes of the conformational mobility in the functional interfacial loops of the receptor-binding domain in the SARS-CoV-2 spike protein can be fine-tuned through crosstalk between convergent mutations which could provide an evolutionary path for modulation of immune escape. By combining atomistic simulations and Markovian modeling analysis with perturbation-based approaches, we determined important complementary roles of convergent mutation sites as effectors and receivers of allosteric signaling involved in modulation of conformational plasticity and regulation of allosteric communications. This study also revealed hidden allosteric pockets and suggested that convergent mutation sites could control evolution and distribution of allosteric pockets through modulation of conformational plasticity in the flexible adaptable regions.

Molecular and Clinical Epidemiology of SARS-CoV-2 Infection among Vaccinated and Unvaccinated Individuals in a Large Healthcare Organization from New Jersey

New Jersey was among the first states impacted by the COVID-19 pandemic, with one of the highest overall death rates in the nation. Nevertheless, relatively few reports have been published focusing specifically on New Jersey. Here we report on molecular, clinical, and epidemiologic observations, from the largest healthcare network in the state, in a cohort of vaccinated and unvaccinated individuals with laboratory-confirmed SARS-CoV-2 infection. We conducted molecular surveillance of SARS-CoV-2-positive nasopharyngeal swabs collected in nine hospitals from December 2020 through June 2022, using both whole genome sequencing (WGS) and a real-time RT-PCR screening assay targeting spike protein mutations found in variants of concern (VOCs) within our region. De-identified clinical data were obtained retrospectively, including demographics, COVID-19 vaccination status, ICU admission, ventilator support, mortality, and medical history. Statistical analyses were performed to identify associations between SARS-CoV-2 variants, vaccination status, clinical outcomes, and medical risk factors. A total of 5007 SARS-CoV-2-positive nasopharyngeal swabs were successfully screened and/or sequenced. Variant screening identified three predominant VOCs, including Alpha (n = 714), Delta (n = 1877), and Omicron (n = 1802). Omicron isolates were further sub-typed as BA.1 (n = 899), BA.2 (n = 853), or BA.4/BA.5 (n = 50); the remaining 614 isolates were classified as "Other". Approximately 31.5% (1577/5007) of the samples were associated with vaccine breakthrough infections, which increased in frequency following the emergence of Delta and Omicron. Severe clinical outcomes included ICU admission (336/5007 = 6.7%), ventilator support (236/5007 = 4.7%), and mortality (430/5007 = 8.6%), with increasing age being the most significant contributor to each (p < 0.001). Unvaccinated individuals accounted for 79.7% (268/336) of ICU admissions, 78.3% (185/236) of ventilator cases, and 74.4% (320/430) of deaths. Highly significant (p < 0.001) increases in mortality were observed in individuals with cardiovascular disease, hypertension, cancer, diabetes, and hyperlipidemia, but not with obesity, thyroid disease, or respiratory disease. Significant differences (p < 0.001) in clinical outcomes were also noted between SARS-CoV-2 variants, including Delta, Omicron BA.1, and Omicron BA.2. Vaccination was associated with significantly improved clinical outcomes in our study, despite an increase in breakthrough infections associated with waning immunity, greater antigenic variability, or both. Underlying comorbidities contributed significantly to mortality in both vaccinated and unvaccinated individuals, with increasing risk based on the total number of comorbidities. Real-time RT-PCR-based screening facilitated timely identification of predominant variants using a minimal number of spike protein mutations, with faster turnaround time and reduced cost compared to WGS. Continued evolution of SARS-CoV-2 variants will likely require ongoing surveillance for new VOCs, with real-time assessment of clinical impact.

ACE2 receptor polymorphism in humans and animals increases the risk of the emergence of SARS-CoV-2 variants during repeated intra- and inter-species host-switching of the virus

Like other coronaviruses, SARS-CoV-2 has ability to spread through human-to-human transmission and to circulate from humans to animals and from animals to humans. A high frequency of SARS-CoV-2 mutations has been observed in the viruses isolated from both humans and animals, suggesting a genetic fitness under positive selection in both ecological niches. The most documented positive selection force driving SARS-CoV-2 mutations is the host-specific immune response. However, after electrostatic interactions with lipid rafts, the first contact between the virus and host proteins is the viral spike-cellular receptor binding. Therefore, it is likely that the first level of selection pressure impacting viral fitness relates to the virus's affinity for its receptor, the angiotensin I converting enzyme 2 (ACE2). Although sufficiently conserved in a huge number of species to support binding of the viral spike with enough affinity to initiate fusion, ACE2 is highly polymorphic both among species and within a species. Here, we provide evidence suggesting that when the viral spike-ACE2 receptor interaction is not optimal, due to host-switching, mutations can be selected to improve the affinity of the spike for the ACE2 expressed by the new host. Notably, SARS-CoV-2 is mutation-prone in the spike receptor binding domain (RBD), allowing a better fit for ACE2 orthologs in animals. It is possibly that this may also be true for rare human alleles of ACE2 when the virus is spreading to billions of people. In this study, we present evidence that human subjects expressing the rare E₃₂₉G allele of ACE2 with higher allele frequencies in European populations exhibit a improved affinity for the SARS-CoV-2 spike N₅₀₁Y variant of the virus. This may suggest that this viral N₅₀₁Y variant emerged in the human population after SARS-CoV-2 had infected a human carrying the rare E₃₂₉G allele of ACE2. In addition, this viral evolution could impact viral replication as well as the ability of the adaptive humoral response to control infection with RBD-specific neutralizing antibodies. In a shifting landscape, this ACE2-driven genetic drift of SARS-CoV-2 which we have named the 'boomerang effect', could complicate the challenge of preventing COVID with a SARS-CoV-2 spike-derived vaccine.

Dissemination of the Omicron Variant and Its Sub-Lineages among Residents and Travelers in Its First Year of Emergence in Venezuela

The emergence of the SARS-CoV-2 Variant of Concern (VOC), Omicron, has been characterized by an explosive number of cases in almost every part of the world. The dissemination of different sub-lineages and recombinant genomes also led to several posterior waves in many countries. The circulation of this VOC and its major sub-lineages (BA.1 to BA.5) was monitored in community cases and in international travelers returning to Venezuela by a rapid partial sequencing method. The specific sub-lineage assignment was performed by complete genome sequencing. Epidemic waves of SARS-CoV-2 cases were observed among international travelers during 2022, a situation not seen before December 2021. The succession of the Omicron VOC sub-lineages BA.1 to BA.5 occurred sequentially, except for BA.3, which was almost not detected. However, the sub-lineages generally circulated two months earlier in international travelers than in community cases. The diversity of Omicron sub-lineages found in international travelers was related to the one found in the USA, consistent with the most frequent destination of international travel from Venezuela this year. These differences are compatible with the delay observed sometimes in Latin American countries in the circulation of the different lineages of the Omicron VOC. Once the sub-lineages were introduced in the country, community transmission was responsible for generating a characteristic distribution of them, with a predominance of sub-lineages not necessarily similar to the one observed in travelers or neighboring countries.

The Lambda Variant in Argentina: Analyzing the Evolution and Spread of SARS-CoV-2 Lineage C.37

The second wave of COVID-19 occurred in South America in early 2021 and was mainly driven by Gamma and Lambda variants. In this study, we aimed to describe the emergence and local genomic diversity of the SARS-CoV-2 Lambda variant in Argentina, from its initial entry into the country until its detection ceased. Molecular surveillance was conducted on 9356 samples from Argentina between October 2020 and April 2022, and sequencing, phylogenetic, and phylogeographic analyses were performed. Our findings revealed that the Lambda variant was first detected in Argentina in January 2021 and steadily increased in frequency until it peaked in April 2021, with continued detection throughout the year. Phylodynamic analyses showed that at least 18 introductions of the Lambda variant into the country occurred, with nine of them having evidence of onward local transmission. The spatial--temporal reconstruction showed that Argentine clades were associated with Lambda sequences from Latin America and suggested an initial diversification in the Metropolitan Area of Buenos Aires before spreading to other regions in Argentina. Genetic analyses of genome sequences allowed us to describe the mutational patterns of the Argentine Lambda sequences and detect the emergence of rare mutations in an immunocompromised patient. Our study highlights the importance of genomic surveillance in identifying the introduction and geographical distribution of the SARS-CoV-2 Lambda variant, as well as in monitoring the emergence of mutations that could be involved in the evolutionary leaps that characterize variants of concern.

Clinical Characteristics Profile of COVID-19 Patients with Omicron Variant Admitted in a Tertiary Hospital, Central China

Purpose

Omicron, a variant of COVID-19, is becoming a major issue of global concern. Its high transmissibility may bring challenges to the distribution of health care in a large population country like China. Investigating the behavior of the virus in the Chinese population will certainly help to plan for the upcoming surge of Omicron. Therefore, we made a preliminary analysis of the clinical and epidemiological characteristics of suspected cases of Omicron at the early stage of the surge.

Patients and Methods

The study was conducted in Nanyang Central Hospital, a tertiary hospital, from 21st December, 2022 to 8th January, 2023. A total of 210 patients underwent demographic characteristics and clinical symptom collection from their medical records. Moreover, sputum culture was also conducted to explore the types of bacterial or fungal infections.

Results

Our results showed that 5 patients (4.1%) were aged 16-49, 40 patients (32.5%) were aged 50-70, and 78 patients (63.4%) were aged 70 or more in the severe group. The proportion of male patients with severe diseases infected with Omicron is higher than that of female patients and the proportion of severe cases increases with age. The main symptoms of patients infected with Omicron are cough (91, 74.0%), fever (90, 73.2%), and asthma (73, 59.3%). The pathogens Streptococcus pneumoniae (71, 31.0%), Staphylococcus aureus (46, 20.1%), Mycoplasma pneumoniae (26, 11.4%), Klebsiella pneumoniae (18, 7.9%), Acinetobacter baumannii (13, 5.7%), and Haemophilus influenzae were detected in lower respiratory tract.

Conclusion

This study suggests that age >70 is a risk factor for severe COVID-19 and that patients often have bacterial or fungal infections. Our research results may help to provide effective treatment for patients with Omicron infection and also contribute to health economic analysis and research to assist future public health decision-making.

Rapid Detection of SARS-CoV-2 RNA Using Reverse Transcription Recombinase Polymerase Amplification (RT-RPA) with Lateral Flow for N-Protein Gene and Variant-Specific Deletion-Insertion Mutation in S-Protein Gene

Rapid molecular testing for severe acute respiratory coronavirus 2 (SARS-CoV-2) variants may contribute to the development of public health measures, particularly in resource-limited areas. Reverse transcription recombinase polymerase amplification using a lateral flow assay (RT-RPA-LF) allows rapid RNA detection without thermal cyclers. In this study, we developed two assays to detect SARS-CoV-2 nucleocapsid (N) gene and Omicron BA.1 spike (S) gene-specific deletion-insertion mutations (del211/ins214). Both tests had a detection limit of 10 copies/µL in vitro and the detection time was approximately 35 min from incubation to detection. The sensitivities of SARS-CoV-2 (N) RT-RPA-LF by viral load categories were 100% for clinical samples with high (>9015.7 copies/µL, cycle quantification (Cq): < 25) and moderate (385.5-9015.7 copies/µL, Cq: 25-29.9) viral load, 83.3% for low (16.5-385.5 copies/µL, Cq: 30-34.9), and 14.3% for very low (<16.5 copies/µL, Cq: 35-40). The sensitivities of the Omicron BA.1 (S) RT-RPA-LF were 94.9%, 78%, 23.8%, and 0%, respectively, and the specificity against non-BA.1 SARS-CoV-2-positive samples was 96%. The assays seemed more sensitive than rapid antigen detection in moderate viral load samples. Although implementation in resource-limited settings requires additional improvements, deletion-insertion mutations were successfully detected by the RT-RPA-LF technique.

Cross-Protection Induced by Highly Conserved Human B, CD4+, and CD8+ T Cell Epitopes-Based Coronavirus Vaccine Against Severe Infection, Disease, and Death Caused by Multiple SARS-CoV-2 Variants of Concern

Background

The Coronavirus disease 2019 (COVID-19) pandemic has created one of the largest global health crises in almost a century. Although the current rate of SARS-CoV-2 infections has decreased significantly; the long-term outlook of COVID-19 remains a serious cause of high death worldwide; with the mortality rate still surpassing even the worst mortality rates recorded for the influenza viruses. The continuous emergence of SARS-CoV-2 variants of concern (VOCs), including multiple heavily mutated Omicron sub-variants, have prolonged the COVID-19 pandemic and outlines the urgent need for a next-generation vaccine that will protect from multiple SARS-CoV-2 VOCs.

Methods

In the present study, we designed a multi-epitope-based Coronavirus vaccine that incorporated B, CD4+, and CD8+ T cell epitopes conserved among all known SARS-CoV-2 VOCs and selectively recognized by CD8+ and CD4+ T-cells from asymptomatic COVID-19 patients irrespective of VOC infection. The safety, immunogenicity, and cross-protective immunity of this pan-Coronavirus vaccine were studied against six VOCs using an innovative triple transgenic h-ACE-2-HLA-A2/DR mouse model.

Results

The Pan-Coronavirus vaccine: (i) is safe; (ii) induces high frequencies of lung-resident functional CD8+ and CD4+ TEM and TRM cells; and (iii) provides robust protection against virus replication and COVID-19-related lung pathology and death caused by six SARS-CoV-2 VOCs: Alpha (B.1.1.7), Beta (B.1.351), Gamma or P1 (B.1.1.28.1), Delta (lineage B.1.617.2) and Omicron (B.1.1.529).

Conclusions

A multi-epitope pan-Coronavirus vaccine bearing conserved human B and T cell epitopes from structural and non-structural SARS-CoV-2 antigens induced cross-protective immunity that cleared the virus, and reduced COVID-19-related lung pathology and death caused by multiple SARS-CoV-2 VOCs.

The Mutational Landscape of SARS-CoV-2

Mutation research is crucial for detecting and treating SARS-CoV-2 and developing vaccines. Using over 5,300,000 sequences from SARS-CoV-2 genomes and custom Python programs, we analyzed the mutational landscape of SARS-CoV-2. Although almost every nucleotide in the SARS-CoV-2 genome has mutated at some time, the substantial differences in the frequency and regularity of mutations warrant further examination. C>U mutations are the most common. They are found in the largest number of variants, pangolin lineages, and countries, which indicates that they are a driving force behind the evolution of SARS-CoV-2. Not all SARS-CoV-2 genes have mutated in the same way. Fewer non-synonymous single nucleotide variations are found in genes that encode proteins with a critical role in virus replication than in genes with ancillary roles. Some genes, such as spike (S) and nucleocapsid (N), show more non-synonymous mutations than others. Although the prevalence of mutations in the target regions of COVID-19 diagnostic RT-qPCR tests is generally low, in some cases, such as for some primers that bind to the N gene, it is significant. Therefore, ongoing monitoring of SARS-CoV-2 mutations is crucial. The SARS-CoV-2 Mutation Portal provides access to a database of SARS-CoV-2 mutations.

Markov State Models and Perturbation-Based Approaches Reveal Distinct Dynamic Signatures and Hidden Allosteric Pockets in the Emerging SARS-Cov-2 Spike Omicron Variants Complexes with the Host Receptor: The Interplay of Dynamics and Convergent Evolution Modulates Allostery and Functional Mechanisms

The new generation of SARS-CoV-2 Omicron variants displayed a significant growth advantage and the increased viral fitness by acquiring convergent mutations, suggesting that the immune pressure can promote convergent evolution leading to the sudden acceleration of SARS-CoV-2 evolution. In the current study, we combined structural modeling, extensive microsecond MD simulations and Markov state models to characterize conformational landscapes and identify specific dynamic signatures of the SARS-CoV-2 spike complexes with the host receptor ACE2 for the recently emerged highly transmissible XBB.1, XBB.1.5, BQ.1, and BQ.1.1 Omicron variants. Microsecond simulations and Markovian modeling provided a detailed characterization of the conformational landscapes and revealed the increased thermodynamic stabilization of the XBB.1.5 subvariant which is contrasted to more dynamic BQ.1 and BQ.1.1 subvariants. Despite considerable structural similarities, Omicron mutations can induce unique dynamic signatures and specific distributions of conformational states. The results suggested that variant-specific changes of conformational mobility in the functional interfacial loops of the spike receptor binding domain can be fine-tuned through cross-talk between convergent mutations thereby providing an evolutionary path for modulation of immune escape. By combining atomistic simulations and Markovian modeling analysis with perturbation-based approaches, we determined important complementary roles of convergent mutation sites as effectors and receivers of allosteric signaling involved in modulating conformational plasticity at the binding interface and regulating allosteric responses. This study also characterized the dynamics-induced evolution of allosteric pockets in the Omicron complexes that revealed hidden allosteric pockets and suggested that convergent mutation sites could control evolution and distribution of allosteric pockets through modulation of conformational plasticity in the flexible adaptable regions. Through integrative computational approaches, this investigation provides a systematic analysis and comparison of the effects of Omicron subvariants on conformational dynamics and allosteric signaling in the complexes with the ACE2 receptor.

For Table of Contents Use Only

Molecular Epidemiology of SARS-CoV-2 during Five COVID-19 Waves and the Significance of Low-Frequency Lineages

SARS-CoV-2 lineages and variants of concern (VOC) have gained more efficient transmission and immune evasion properties with time. We describe the circulation of VOCs in South Africa and the potential role of low-frequency lineages on the emergence of future lineages. Whole genome sequencing was performed on SARS-CoV-2 samples from South Africa. Sequences were analysed with Nextstrain pangolin tools and Stanford University Coronavirus Antiviral & Resistance Database. In 2020, 24 lineages were detected, with B.1 (3%; 8/278), B.1.1 (16%; 45/278), B.1.1.348 (3%; 8/278), B.1.1.52 (5%; 13/278), C.1 (13%; 37/278) and C.2 (2%; 6/278) circulating during the first wave. Beta emerged late in 2020, dominating the second wave of infection. B.1 and B.1.1 continued to circulate at low frequencies in 2021 and B.1.1 re-emerged in 2022. Beta was outcompeted by Delta in 2021, which was thereafter outcompeted by Omicron sub-lineages during the 4th and 5th waves in 2022. Several significant mutations identified in VOCs were also detected in low-frequency lineages, including S68F (E protein); I82T (M protein); P13L, R203K and G204R/K (N protein); R126S (ORF3a); P323L (RdRp); and N501Y, E484K, D614G, H655Y and N679K (S protein). Low-frequency variants, together with VOCs circulating, may lead to convergence and the emergence of future lineages that may increase transmissibility, infectivity and escape vaccine-induced or natural host immunity.

Divergent mutations of Delta and Omicron variants: key players behind differential viral attributes across the COVID-19 waves

The third SARS-CoV-2 pandemic wave causing Omicron variant has comparatively higher replication rate and transmissibility than the second wave-causing Delta variant. The exact mechanism behind the differential properties of Delta and Omicron in respect to infectivity and virulence is not properly understood yet. This study reports the analysis of different mutations within the receptor binding domain (RBD) of spike glycoprotein and non-structural protein (nsp) of Delta and Omicron strains. We have used computational studies to evaluate the properties of Delta and Omicron variants in this work. Q498R, Q493R and S375F mutations of RBD showed better docking scores for Omicron compared to Delta variant of SARS-CoV-2, whereas nsp3_L1266I with PARP15 (7OUX), nsp3_L1266I with PARP15 (7OUX), and nsp6_G107 with ISG15 (1Z2M) showed significantly higher docking score. The findings of the present study might be helpful to reveal the probable cause of relatively milder form of COVID-19 disease manifested by Omicron in comparison to Delta variant of SARS-CoV-2 virus.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13337-023-00823-0.

Balancing Functional Tradeoffs between Protein Stability and ACE2 Binding in the SARS-CoV-2 Omicron BA.2, BA.2.75 and XBB Lineages: Dynamics-Based Network Models Reveal Epistatic Effects Modulating Compensatory Dynamic and Energetic Changes

Evolutionary and functional studies suggested that the emergence of the Omicron variants can be determined by multiple fitness trade-offs including the immune escape, binding affinity for ACE2, conformational plasticity, protein stability and allosteric modulation. In this study, we systematically characterize conformational dynamics, structural stability and binding affinities of the SARS-CoV-2 Spike Omicron complexes with the host receptor ACE2 for BA.2, BA.2.75, XBB.1 and XBB.1.5 variants. We combined multiscale molecular simulations and dynamic analysis of allosteric interactions together with the ensemble-based mutational scanning of the protein residues and network modeling of epistatic interactions. This multifaceted computational study characterized molecular mechanisms and identified energetic hotspots that can mediate the predicted increased stability and the enhanced binding affinity of the BA.2.75 and XBB.1.5 complexes. The results suggested a mechanism driven by the stability hotspots and a spatially localized group of the Omicron binding affinity centers, while allowing for functionally beneficial neutral Omicron mutations in other binding interface positions. A network-based community model for the analysis of epistatic contributions in the Omicron complexes is proposed revealing the key role of the binding hotspots R498 and Y501 in mediating community-based epistatic couplings with other Omicron sites and allowing for compensatory dynamics and binding energetic changes. The results also showed that mutations in the convergent evolutionary hotspot F486 can modulate not only local interactions but also rewire the global network of local communities in this region allowing the F486P mutation to restore both the stability and binding affinity of the XBB.1.5 variant which may explain the growth advantages over the XBB.1 variant. The results of this study are consistent with a broad range of functional studies rationalizing functional roles of the Omicron mutation sites that form a coordinated network of hotspots enabling a balance of multiple fitness tradeoffs and shaping up a complex functional landscape of virus transmissibility.

Molnupiravir, Nirmatrelvir/Ritonavir, or Sotrovimab for High-Risk COVID-19 Patients Infected by the Omicron Variant: Hospitalization, Mortality, and Time until Negative Swab Test in Real Life

Background. Several drugs which are easy to administer in outpatient settings have been authorized and endorsed for high-risk COVID-19 patients with mild-moderate disease to prevent hospital admission and death, complementing COVID-19 vaccines. However, the evidence on the efficacy of COVID-19 antivirals during the Omicron wave is scanty or conflicting. Methods. This retrospective controlled study investigated the efficacy of Molnupiravir or Nirmatrelvir/Ritonavir (Paxlovid®) or Sotrovimab against standard of care (controls) on three different endpoints among 386 high-risk COVID-19 outpatients: hospital admission at 30 days; death at 30 days; and time between COVID-19 diagnosis and first negative swab test result. Multivariable logistic regression was employed to investigate the determinants of hospitalization due to COVID-19-associated pneumonia, whereas time to first negative swab test result was investigated by means of multinomial logistic analysis as well as Cox regression analysis. Results. Only 11 patients (overall rate of 2.8%) developed severe COVID-19-associated pneumonia requiring admission to hospital: 8 controls (7.2%); 2 patients on Nirmatrelvir/Ritonavir (2.0%); and 1 on Sotrovimab (1.8%). No patient on Molnupiravir was institutionalized. Compared to controls, hospitalization was less likely for patients on Nirmatrelvir/Ritonavir (aOR = 0.16; 95% CI: 0.03; 0.89) or Molnupiravir (omitted estimate); drug efficacy was 84% for Nirmatrelvir/Ritonavir against 100% for Molnupiravir. Only two patients died of COVID-19 (rate of 0.5%), both were controls, one (a woman aged 96 years) was unvaccinated and the other (a woman aged 72 years) had adequate vaccination status. At Cox regression analysis, the negativization rate was significantly higher in patients treated with both antivirals-Nirmatrelvir/Ritonavir (aHR = 1.68; 95% CI: 1.25; 2.26) or Molnupiravir (aHR = 1.45; 95% CI: 1.08; 1.94). However, COVID-19 vaccination with three (aHR = 2.03; 95% CI: 1.51; 2.73) or four (aHR = 2.48; 95% CI: 1.32; 4.68) doses had a slightly stronger effect size on viral clearance. In contrast, the negativization rate reduced significantly in patients who were immune-depressed (aHR = 0.70; 95% CI: 0.52; 0.93) or those with a Charlson index ≥5 (aHR = 0.63; 0.41; 0.95) or those who had started the respective treatment course 3+ days after COVID-19 diagnosis (aOR = 0.56; 95% CI: 0.38; 0.82). Likewise, at internal analysis (excluding patients on standard of care), patients on Molnupiravir (aHR = 1.74; 95% CI: 1.21; 2.50) or Nirmatrelvir/Ritonavir (aHR = 1.96; 95% CI: 1.32; 2.93) were more likely to turn negative earlier than those on Sotrovimab (reference category). Nonetheless, three (aHR = 1.91; 95% CI: 1.33; 2.74) or four (aHR = 2.20; 95% CI: 1.06; 4.59) doses of COVID-19 vaccine were again associated with a faster negativization rate. Again, the negativization rate was significantly lower if treatment started 3+ days after COVID-19 diagnosis (aHR = 0.54; 95% CI: 0.32; 0.92). Conclusions. Molnupiravir, Nirmatrelvir/Ritonavir, and Sotrovimab were all effective in preventing hospital admission and/or mortality attributable to COVID-19. However, hospitalizations also decreased with higher number of doses of COVID-19 vaccines. Although they are effective against severe disease and mortality, the prescription of COVID-19 antivirals should be carefully scrutinized by double opinion, not only to contain health care costs but also to reduce the risk of generating resistant SARS-CoV-2 strains. Only 64.7% of patients were in fact immunized with 3+ doses of COVID-19 vaccines in the present study. High-risk patients should prioritize COVID-19 vaccination, which is a more cost-effective approach than antivirals against severe SARS-CoV-2 pneumonia. Likewise, although both antivirals, especially Nirmatrelvir/Ritonavir, were more likely than standard of care and Sotrovimab to reduce viral shedding time (VST) in high-risk SARS-CoV-2 patients, vaccination had an independent and stronger effect on viral clearance. However, the effect of antivirals or COVID-19 vaccination on VST should be considered a secondary benefit. Indeed, recommending Nirmatrelvir/Ritonavir in order to control VST in high-risk COVID-19 patients is rather questionable since other cheap, large spectrum and harmless nasal disinfectants such as hypertonic saline solutions are available on the market with proven efficacy in containing VST.

A Web Tool to Estimate Baseline Anti-Spike Monoclonal Antibody Efficacy Based on Regional Genomic Surveillance

Drug appropriateness is a pillar of modern evidence-based medicine, but the turnaround times of genomic sequencing are not compatible with the urgent need to deliver treatments against microorganisms. Massive worldwide genomic surveillance has created an unprecedented landscape for exploiting viral sequencing for therapeutic purposes. When it comes to therapeutic antiviral antibodies, using IC50 against specific polymorphisms of the target antigen can be calculated in vitro, and a list of mutations leading to drug resistance (immune escape) can be compiled. The author encountered this type of knowledge (available from the Stanford University Coronavirus Antiviral Resistance Database,) in a publicly accessible repository of SARS-CoV-2 sequences. The author used a custom function of the CoV-Spectrum.org web portal to deliver up-to-date, regional prevalence estimates of baseline efficacy for each authorized anti-spike mAb across all co-circulating SARS-CoV-2 sublineages at a given time point. This publicly accessible tool can inform therapeutic choices that would otherwise be blind.

COVID-19 and Diarylamidines: The Parasitic Connection

As emerging severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) variants (Omicron) continue to outpace and negate combinatorial vaccines and monoclonal antibody therapies targeting the spike protein (S) receptor binding domain (RBD), the appetite for developing similar COVID-19 treatments has significantly diminished, with the attention of the scientific community switching to long COVID treatments. However, treatments that reduce the risk of "post-COVID-19 syndrome" and associated sequelae remain in their infancy, particularly as no established criteria for diagnosis currently exist. Thus, alternative therapies that reduce infection and prevent the broad range of symptoms associated with 'post-COVID-19 syndrome' require investigation. This review begins with an overview of the parasitic-diarylamidine connection, followed by the renin-angiotensin system (RAS) and associated angiotensin-converting enzyme 2 (ACE2) and transmembrane serine protease 2 (TMPRSSR2) involved in SARS-CoV-2 infection. Subsequently, the ability of diarylamidines to inhibit S-protein binding and various membrane serine proteases associated with SARS-CoV-2 and parasitic infections are discussed. Finally, the roles of diarylamidines (primarily DIZE) in vaccine efficacy, epigenetics, and the potential amelioration of long COVID sequelae are highlighted.

The rapid rise of SARS-CoV-2 Omicron subvariants with immune evasion properties: XBB.1.5 and BQ.1.1 subvariants

As the fifth variant of concern of the SARS-CoV-2 virus, the Omicron variant (B.1.1.529) has quickly become the dominant type among the previous circulating variants worldwide. During the Omicron wave, several subvariants have emerged, with some exhibiting greater infectivity and immune evasion, accounting for their fast spread across many countries. Recently, two Omicron subvariants, BQ.1 and XBB lineages, including BQ.1.1, XBB.1, and XBB.1.5, have become a global public health issue given their ability to escape from therapeutic monoclonal antibodies and herd immunity induced by prior coronavirus disease 2019 (COVID-19) vaccines, boosters, and infection. In this respect, XBB.1.5, which has been established to harbor a rare mutation F486P, demonstrates superior transmissibility and immune escape ability compared to other subvariants and has emerged as the dominant strain in several countries. This review provides a comprehensive overview of the epidemiological features, spike mutations, and immune evasion of BQ.1 and XBB lineages. We expounded on the mechanisms underlying mutations and immune escape from neutralizing antibodies from vaccinated or convalescent COVID-19 individuals and therapeutic monoclonal antibodies (mAbs) and proposed strategies for prevention against BQ.1 and XBB sublineages.

Deep Structural Analysis of Myriads of Omicron Sub-Variants Revealed Hotspot for Vaccine Escape Immunity

The emergence of the Omicron variant has reinforced the importance of continued SARS-CoV-2 evolution and its possible impact on vaccine effectiveness. Specifically, mutations in the receptor-binding domain (RBD) are critical to comprehend the flexibility and dynamicity of the viral interaction with the human agniotensin-converting enzyme 2 (hACE2) receptor. To this end, we have applied a string of deep structural and genetic analysis tools to map the substitution patterns in the S protein of major Omicron sub-variants (n = 51) with a primary focus on the RBD mutations. This head-to-head comparison of Omicron sub-variants revealed multiple simultaneous mutations that are attributed to antibody escape, and increased affinity and binding to hACE2. Our deep mapping of the substitution matrix indicated a high level of diversity at the N-terminal and RBD domains compared with other regions of the S protein, highlighting the importance of these two domains in a matched vaccination approach. Structural mapping identified highly variable mutations in the up confirmation of the S protein and at sites that critically define the function of the S protein in the virus pathobiology. These substitutional trends offer support in tracking mutations along the evolutionary trajectories of SAR-CoV-2. Collectively, the findings highlight critical areas of mutations across the major Omicron sub-variants and propose several hotspots in the S proteins of SARS-CoV-2 sub-variants to train the future design and development of COVID-19 vaccines.