The Mutational Landscape of SARS-CoV-2

SARS-CoV-2 main protease (M-pro) mutational profiling: An insight into mutation coldspots

SARS-CoV-2 and the COVID-19 pandemic have marked a milestone in the history of scientific research worldwide. To ensure that treatments are successful in the mid-long term, it is crucial to characterize SARS-CoV-2 mutations, as they might lead to viral resistance. Data from >5,700,000 SARS-CoV-2 genomes available at GISAID was used to report SARS-CoV-2 mutations. Given the pivotal role of its main protease (M-pro) in virus replication, a detailed analysis of SARS-CoV-2 M-pro mutations was conducted, with particular attention to mutation-resistant residues or mutation coldspots, defined as those residues that have mutated in five or fewer genomes. 32 mutation coldspots were identified, most of which mediate interprotomer interactions or funneling interaction networks from the substrate-binding site towards the dimerization surface and vice versa. Besides, mutation coldspots were virtually conserved in all main proteases from other CoVs. Our results provide valuable information about key residues to M-pro structure that could be useful in rational target-directed drug design and establish a solid groundwork based on mutation analyses for the inhibition of M-pro dimerization, with a potential applicability to future coronavirus outbreaks.

Dynamic expedition of leading mutations in SARS-CoV-2 spike glycoproteins

The continuous evolution of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), which caused the recent pandemic, has generated countless new variants with varying fitness. Mutations of the spike glycoprotein play a particularly vital role in shaping its evolutionary trajectory, as they have the capability to alter its infectivity and antigenicity. We present a time-resolved statistical method, Dynamic Expedition of Leading Mutations (deLemus), to analyze the evolutionary dynamics of the SARS-CoV-2 spike glycoprotein. The proposed L -index of the deLemus method is effective in quantifying the mutation strength of each amino acid site and outlining evolutionarily significant sites, allowing the comprehensive characterization of the evolutionary mutation pattern of the spike glycoprotein.

Impact of Vaccination on Intra-Host Genetic Diversity of Patients Infected with SARS-CoV-2 Gamma Lineage

The high transmissibility, rapid evolution, and immune escape of SARS-CoV-2 variants can influence the course of infection and, in turn, morbidity and mortality in COVID-19, posing a challenge in controlling transmission rates and contributing to the emergence and spread of new variants. Understanding the factors that shape viral genetic variation is essential for comprehending the evolution and transmission of SARS-CoV-2, especially in vaccinated individuals where immune response plays a role in the progression and spread of this disease. In this context, we evaluated the impact of immunity induced by the CoronaVac vaccine (Butantan/Sinovac) on intra-host genetic diversity, analyzing 118 whole-genome sequences of SARS-CoV-2 from unvaccinated and vaccinated patients infected with the Gamma variant. Vaccination with CoronaVac favors negative selection at the intra-host level in different genomic regions. It prevents greater genetic diversity of SARS-CoV-2, reinforcing the importance of vaccination in reducing the emergence of new mutations and virus transmission.

Refining SARS-CoV-2 intra-host variation by leveraging large-scale sequencing data

Understanding viral genome evolution during host infection is crucial for grasping viral diversity and evolution. Analyzing intra-host single nucleotide variants (iSNVs) offers insights into new lineage emergence, which is important for predicting and mitigating future viral threats. Despite next-generation sequencing's potential, challenges persist, notably sequencing artifacts leading to false iSNVs. We developed a workflow to enhance iSNV detection in large NGS libraries, using over 130 000 SARS-CoV-2 libraries to distinguish mutations from errors. Our approach integrates bioinformatics protocols, stringent quality control, and dimensionality reduction to tackle batch effects and improve mutation detection reliability. Additionally, we pioneer the application of the PHATE visualization approach to genomic data and introduce a methodology that quantifies how related groups of data points are represented within a two-dimensional space, enhancing clustering structure explanation based on genetic similarities. This workflow advances accurate intra-host mutation detection, facilitating a deeper understanding of viral diversity and evolution.

Modulation of biophysical properties of nucleocapsid protein in the mutant spectrum of SARS-CoV-2

Genetic diversity is a hallmark of RNA viruses and the basis for their evolutionary success. Taking advantage of the uniquely large genomic database of SARS-CoV-2, we examine the impact of mutations across the spectrum of viable amino acid sequences on the biophysical phenotypes of the highly expressed and multifunctional nucleocapsid protein. We find variation in the physicochemical parameters of its extended intrinsically disordered regions (IDRs) sufficient to allow local plasticity, but also observe functional constraints that similarly occur in related coronaviruses. In biophysical experiments with several N-protein species carrying mutations associated with major variants, we find that point mutations in the IDRs can have nonlocal impact and modulate thermodynamic stability, secondary structure, protein oligomeric state, particle formation, and liquid-liquid phase separation. In the Omicron variant, distant mutations in different IDRs have compensatory effects in shifting a delicate balance of interactions controlling protein assembly properties, and include the creation of a new protein-protein interaction interface in the N-terminal IDR through the defining P13L mutation. A picture emerges where genetic diversity is accompanied by significant variation in biophysical characteristics of functional N-protein species, in particular in the IDRs.

Detection and Characterisation of SARS-CoV-2 in Eastern Province of Zambia: A Retrospective Genomic Surveillance Study

Mutations have driven the evolution and development of new variants of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) with potential implications for increased transmissibility, disease severity and vaccine escape among others. Genome sequencing is a technique that allows scientists to read the genetic code of an organism and has become a powerful tool for studying emerging infectious diseases. Here, we conducted a cross-sectional study in selected districts of the Eastern Province of Zambia, from November 2021 to February 2022. We analyzed SARS-CoV-2 samples (n = 76) using high-throughput sequencing. A total of 4097 mutations were identified in 69 SARS-CoV-2 genomes with 47% (1925/4097) of the mutations occurring in the spike protein. We identified 83 unique amino acid mutations in the spike protein of the seven Omicron sublineages (BA.1, BA.1.1, BA.1.14, BA.1.18, BA.1.21, BA.2, BA.2.23 and XT). Of these, 43.4% (36/83) were present in the receptor binding domain, while 14.5% (12/83) were in the receptor binding motif. While we identified a potential recombinant XT strain, the highly transmissible BA.2 sublineage was more predominant (40.8%). We observed the substitution of other variants with the Omicron strain in the Eastern Province. This work shows the importance of pandemic preparedness and the need to monitor disease in the general population.

Upping the ante: enhanced expression of interferon-antagonizing ORF6 and ORF9b proteins by SARS-CoV-2 variants of concern

SARS-CoV-2 exhibits a remarkable capability to subvert the host antiviral innate immune system. This adeptness is orchestrated by viral proteins, which initially attempt to obstruct the activation of the antiviral immune program and then act as a fail-safe mechanism to mitigate the downstream effects of the activated immune response. This dual strategy leads to delayed expression and enfeebled action of type-I and -III interferons at the infection site, enabling the virus to replicate extensively in the lungs and subsequently disseminate to other organs. Throughout the course of the COVID-19 pandemic, SARS-CoV-2 has undergone evolution, giving rise to several variants of concern, some with exceedingly higher transmission and virulence. These improved features have been linked, at least in part, to the heightened expression or activity of specific viral proteins involved in circumventing host defense mechanisms. In this review, we aim to provide a concise summary of two SARS-CoV-2 proteins, ORF6 and ORF9b, which provided selective advantage to certain variants, affecting their biology and pathogenesis.

Role of SARS-CoV-2 mutations in the evolution of the COVID-19 pandemic

OBJECTIVES

The SARS-CoV-2 pandemic and large-scale genomic surveillance provided an exceptional opportunity to analyze mutations that appeared over three years in viral genomes. Here we studied mutations and their epidemic consequences for SARS-CoV-2 genomes from our center.

METHODS

We analyzed 61,397 SARS-CoV-2 genomes we sequenced from respiratory samples for genomic surveillance. Mutations frequencies were calculated using Nextclade, Microsoft Excel, and an in-house Python script.

RESULTS

A total of 22,225 nucleotide mutations were identified, 220 (1.0%) being each at the root of ≥836 genomes, classifying mutations as 'hyperfertile'. Two seeded the European pandemic: P323L in RNA polymerase, associated with an increased mutation rate, and D614G in spike that improved fitness. Most 'hyperfertile' mutations occurred in areas not predicted with increased virulence. Their mean number was 8±6 (0-22) per 1000 nucleotides per gene. They were 3.7-times more frequent in accessory than informational genes (13.8 versus 3.7/1000 nucleotides). Particularly, they were 4.1-times more frequent in ORF8 than in the RNA polymerase gene. Interestingly, stop codons were present in 97 positions, almost only in accessory genes, including ORF8 (21/100 codons).

CONCLUSIONS

most 'hyperfertile' mutations did not predict emergence of a new epidemic, and some were stop codons indicating the existence of so-named 'non-virulence' genes.

Properties and Mechanisms of Deletions, Insertions, and Substitutions in the Evolutionary History of SARS-CoV-2

SARS-CoV-2 has accumulated many mutations since its emergence in late 2019. Nucleotide substitutions leading to amino acid replacements constitute the primary material for natural selection. Insertions, deletions, and substitutions appear to be critical for coronavirus's macro- and microevolution. Understanding the molecular mechanisms of mutations in the mutational hotspots (positions, loci with recurrent mutations, and nucleotide context) is important for disentangling roles of mutagenesis and selection. In the SARS-CoV-2 genome, deletions and insertions are frequently associated with repetitive sequences, whereas C>U substitutions are often surrounded by nucleotides resembling the APOBEC mutable motifs. We describe various approaches to mutation spectra analyses, including the context features of RNAs that are likely to be involved in the generation of recurrent mutations. We also discuss the interplay between mutations and natural selection as a complex evolutionary trend. The substantial variability and complexity of pipelines for the reconstruction of mutations and the huge number of genomic sequences are major problems for the analyses of mutations in the SARS-CoV-2 genome. As a solution, we advocate for the development of a centralized database of predicted mutations, which needs to be updated on a regular basis.

Modulation of Biophysical Properties of Nucleocapsid Protein in the Mutant Spectrum of SARS-CoV-2

Genetic diversity is a hallmark of RNA viruses and the basis for their evolutionary success. Taking advantage of the uniquely large genomic database of SARS-CoV-2, we examine the impact of mutations across the spectrum of viable amino acid sequences on the biophysical phenotypes of the highly expressed and multifunctional nucleocapsid protein. We find variation in the physicochemical parameters of its extended intrinsically disordered regions (IDRs) sufficient to allow local plasticity, but also exhibiting functional constraints that similarly occur in related coronaviruses. In biophysical experiments with several N-protein species carrying mutations associated with major variants, we find that point mutations in the IDRs can have nonlocal impact and modulate thermodynamic stability, secondary structure, protein oligomeric state, particle formation, and liquid-liquid phase separation. In the Omicron variant, distant mutations in different IDRs have compensatory effects in shifting a delicate balance of interactions controlling protein assembly properties, and include the creation of a new protein-protein interaction interface in the N-terminal IDR through the defining P13L mutation. A picture emerges where genetic diversity is accompanied by significant variation in biophysical characteristics of functional N-protein species, in particular in the IDRs.

Inhibitory Efficacy of Main Components of Scutellaria baicalensis on the Interaction between Spike Protein of SARS-CoV-2 and Human Angiotensin-Converting Enzyme II

Blocking the interaction between the SARS-CoV-2 spike protein and the human angiotensin-converting enzyme II (hACE2) protein serves as a therapeutic strategy for treating COVID-19. Traditional Chinese medicine (TCM) treatments containing bioactive products could alleviate the symptoms of severe COVID-19. However, the emergence of SARS-CoV-2 variants has complicated the process of developing broad-spectrum drugs. As such, the aim of this study was to explore the efficacy of TCM treatments against SARS-CoV-2 variants through targeting the interaction of the viral spike protein with the hACE2 receptor. Antiviral activity was systematically evaluated using a pseudovirus system. Scutellaria baicalensis (S. baicalensis) was found to be effective against SARS-CoV-2 infection, as it mediated the interaction between the viral spike protein and the hACE2 protein. Moreover, the active molecules of S. baicalensis were identified and analyzed. Baicalein and baicalin, a flavone and a flavone glycoside found in S. baicalensis, respectively, exhibited strong inhibitory activities targeting the viral spike protein and the hACE2 protein, respectively. Under optimized conditions, virus infection was inhibited by 98% via baicalein-treated pseudovirus and baicalin-treated hACE2. In summary, we identified the potential SARS-CoV-2 inhibitors from S. baicalensis that mediate the interaction between the Omicron spike protein and the hACE2 receptor. Future studies on the therapeutic application of baicalein and baicalin against SARS-CoV-2 variants are needed.

Diversity of short linear interaction motifs in SARS-CoV-2 nucleocapsid protein

IMPORTANCE

Short linear motifs (SLiMs) are 3-10 amino acid long binding motifs in intrinsically disordered protein regions (IDRs) that serve as ubiquitous protein-protein interaction modules in eukaryotic cells. Through molecular mimicry, viruses hijack these sequence motifs to control host cellular processes. It is thought that the small size of SLiMs and the high mutation frequencies of viral IDRs allow rapid host adaptation. However, a salient characteristic of RNA viruses, due to high replication errors, is their obligate existence as mutant swarms. Taking advantage of the uniquely large genomic database of SARS-CoV-2, here, we analyze the role of sequence diversity in the presentation of SLiMs, focusing on the highly abundant, multi-functional nucleocapsid protein. We find that motif mimicry is a highly dynamic process that produces an abundance of motifs transiently present in subsets of mutant species. This diversity allows the virus to efficiently explore eukaryotic motifs and evolve the host-virus interface.

Regulatory variants of APOBEC3 genes potentially associate with COVID-19 severity in populations with African ancestry

Since November 2019, the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), has caused the worldwide pandemic of the coronavirus disease 2019 (COVID-19), the impact of which is huge to the lives of world populations. Many studies suggested that such situation will continue due to the endless mutations in SARS-CoV-2 genome that result in complexity of the efforts for the control of SARS-CoV-2, since the special enrichment of nucleotide substitution C>U in SARS-CoV-2 sequences were discovered mainly due to the editing by human host factors APOBEC3 genes. The observation of SARS-CoV-2 variants Beta (B.1.351) and Omicron (B.1.1.529) firstly spreading in South Africa promoted us to hypothesize that genetic variants of APOBEC3 special in African populations may be attributed to the higher mutation rate of SARS-CoV-2 variants in Africa. Current study was conducted to search for functional variants of APOBEC3 genes associate with COVID-19 hospitalization in African population. By integrating data from the 1000 Genomes Project, Genotype-Tissue Expression (GTEx), and Host Genetics Initiative (HGI) of COVID-19, we identified potential functional SNPs close to APOBEC3 genes that are associated with COVID-19 hospitalization in African but not with other populations. Our study provides new insights on the potential contribution of APOBEC3 genes on the evolution of SARS-CoV-2 mutations in African population, but further replication is needed to confirm our results.

A Computational Workflow to Predict Biological Target Mutations: The Spike Glycoprotein Case Study

The biological target identification process, a pivotal phase in the drug discovery workflow, becomes particularly challenging when mutations affect proteins' mechanisms of action. COVID-19 Spike glycoprotein mutations are known to modify the affinity toward the human angiotensin-converting enzyme ACE2 and several antibodies, compromising their neutralizing effect. Predicting new possible mutations would be an efficient way to develop specific and efficacious drugs, vaccines, and antibodies. In this work, we developed and applied a computational procedure, combining constrained logic programming and careful structural analysis based on the Structural Activity Relationship (SAR) approach, to predict and determine the structure and behavior of new future mutants. "Mutations rules" that would track statistical and functional types of substitutions for each residue or combination of residues were extracted from the GISAID database and used to define constraints for our software, having control of the process step by step. A careful molecular dynamics analysis of the predicted mutated structures was carried out after an energy evaluation of the intermolecular and intramolecular interactions using the HINT (Hydrophatic INTeraction) force field. Our approach successfully predicted, among others, known Spike mutants.

Diversity of Short Linear Interaction Motifs in SARS-CoV-2 Nucleocapsid Protein

Molecular mimicry of short linear interaction motifs has emerged as a key mechanism for viral proteins binding host domains and hijacking host cell processes. Here, we examine the role of RNA-virus sequence diversity in the dynamics of the virus-host interface, by analyzing the uniquely vast sequence record of viable SARS-CoV-2 species with focus on the multi-functional nucleocapsid protein. We observe the abundant presentation of motifs encoding several essential host protein interactions, alongside a majority of possibly non-functional and randomly occurring motif sequences absent in subsets of viable virus species. A large number of motifs emerge ex nihilo through transient mutations relative to the ancestral consensus sequence. The observed mutational landscape implies an accessible motif space that spans at least 25% of known eukaryotic motifs. This reveals motif mimicry as a highly dynamic process with the capacity to broadly explore host motifs, allowing the virus to rapidly evolve the virus-host interface.

The effects of amino acid substitution of spike protein and genomic recombination on the evolution of SARS-CoV-2

Over three years' pandemic of 2019 novel coronavirus disease (COVID-19), multiple variants and novel subvariants have emerged successively, outcompeted earlier variants and become predominant. The sequential emergence of variants reflects the evolutionary process of mutation-selection-adaption of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). Amino acid substitution/insertion/deletion in the spike protein causes altered viral antigenicity, transmissibility, and pathogenicity of SARS-CoV-2. Early in the pandemic, D614G mutation conferred virus with advantages over previous variants and increased transmissibility, and it also laid a conservative background for subsequent substantial mutations. The role of genomic recombination in the evolution of SARS-CoV-2 raised increasing concern with the occurrence of novel recombinants such as Deltacron, XBB.1.5, XBB.1.9.1, and XBB.1.16 in the late phase of pandemic. Co-circulation of different variants and co-infection in immunocompromised patients accelerate the emergence of recombinants. Surveillance for SARS-CoV-2 genomic variations, particularly spike protein mutation and recombination, is essential to identify ongoing changes in the viral genome and antigenic epitopes and thus leads to the development of new vaccine strategies and interventions.