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

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.