A comprehensive evolutionary and epidemiological characterization of insertion and deletion mutations in SARS-CoV-2 genomes

Evolution of the SARS-CoV-2 Omicron Variants: Genetic Impact on Viral Fitness

Over the last three years, the pandemic of COVID-19 has had a significant impact on people's lives and the global economy. The incessant emergence of variant strains has compounded the challenges associated with the management of COVID-19. As the predominant variant from late 2021 to the present, Omicron and its sublineages, through continuous evolution, have demonstrated iterative viral fitness. The comprehensive elucidation of the biological implications that catalyzed this evolution remains incomplete. In accordance with extant research evidence, we provide a comprehensive review of subvariants of Omicron, delineating alterations in immune evasion, cellular infectivity, and the cross-species transmission potential. This review seeks to clarify the underpinnings of biology within the evolution of SARS-CoV-2, thereby providing a foundation for strategic considerations in the post-pandemic era of COVID-19.

Identifying featured indels associated with SARS-CoV-2 fitness

As an RNA virus, severe acute respiratory coronavirus 2 (SARS-CoV-2) is known for frequent substitution mutations, and substitutions in important genome regions are often associated with viral fitness. However, whether indel mutations are related to viral fitness is generally ignored. Here we developed a computational methodology to investigate indels linked to fitness occurring in over 9 million SARS-CoV-2 genomes. Remarkably, by analyzing 31,642,404 deletion records and 1,981,308 insertion records, our pipeline identified 26,765 deletion types and 21,054 insertion types and discovered 65 indel types with a significant association with Pango lineages. We proposed the concept of featured indels representing the population of specific Pango lineages and variants as substitution mutations and termed these 65 indels as featured indels. The selective pressure of all indel types is assessed using the Bayesian model to explore the importance of indels. Our results exhibited higher selective pressure of indels like substitution mutations, which are important for assessing viral fitness and consistent with previous studies in vitro. Evaluation of the growth rate of each viral lineage indicated that indels play key roles in SARS-CoV-2 evolution and deserve more attention as substitution mutations. IMPORTANCE The fitness of indels in pathogen genome evolution has rarely been studied. We developed a computational methodology to investigate the severe acute respiratory coronavirus 2 genomes and analyze over 33 million records of indels systematically, ultimately proposing the concept of featured indels that can represent specific Pango lineages and identifying 65 featured indels. Machine learning model based on Bayesian inference and viral lineage growth rate evaluation suggests that these featured indels exhibit selection pressure comparable to replacement mutations. In conclusion, indels are not negligible for evaluating viral fitness.

Tracking the evolution of the SARS-CoV-2 Delta variant of concern: analysis of genetic diversity and selection across the whole viral genome

Background

Since 2019, severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has diversified extensively, producing five highly virulent lineages designated as variants of concern (VOCs). The Delta VOC emerged in India with increased transmission, immune evasion, and mortality, causing a massive global case surge in 2021. This study aims to understand how the Delta VOC evolved by characterizing mutation patterns in the viral population before and after its emergence. Furthermore, we aim to identify the influence of positive and negative selection on VOC evolution and understand the prevalence of different mutation types in the viral genome.

Methods

Three groups of whole viral genomes were retrieved from GISAID, sourced from India, with collection periods as follows: Group A-during the initial appearance of SARS-CoV-2; Group B-just before the emergence of the Delta variant; Group C-after the establishment of the Delta variant in India. Mutations in >1% of each group were identified with BioEdit to reveal differences in mutation quantity and type. Sites under positive or negative selection were identified with FUBAR. The results were compared to determine how mutations correspond with selective pressures and how viral mutation profiles changed to reflect genetic diversity before and after VOC emergence.

Results

The number of mutations increased progressively in Groups A-C, with Group C reporting a 2.2- and 1.9-fold increase from Groups A and B, respectively. Among all the observed mutations, Group C had the highest percentage of deletions (22.7%; vs. 4.2% and 2.6% in Groups A and B, respectively), and most mutations altered the final amino acid code, such as non-synonymous substitutions and deletions. Conversely, Group B had the most synonymous substitutions that are effectively silent. The number of sites experiencing positive selection increased in Groups A-C, but Group B had 2.4- and 2.6 times more sites under negative selection compared to Groups A and C, respectively.

Conclusion

Our findings demonstrated that viral genetic diversity continuously increased during and after the emergence of the Delta VOC. Despite this, Group B reports heightened negative selection, which potentially preserves important gene regions during evolution. Group C contains an unprecedented quantity of mutations and positively selected sites, providing strong evidence of active viral adaptation in the population.

High-resolution mapping reveals the mechanism and contribution of genome insertions and deletions to RNA virus evolution

RNA viruses rapidly adapt to selective conditions due to the high intrinsic mutation rates of their RNA-dependent RNA polymerases (RdRps). Insertions and deletions (indels) in viral genomes are major contributors to both deleterious mutational load and evolutionary novelty, but remain understudied. To characterize the mechanistic details of their formation and evolutionary dynamics during infection, we developed a hybrid experimental-bioinformatic approach. This approach, called MultiMatch, extracts insertions and deletions from ultradeep sequencing experiments, including those occurring at extremely low frequencies, allowing us to map their genomic distribution and quantify the rates at which they occur. Mapping indel mutations in adapting poliovirus and dengue virus populations, we determine the rates of indel generation and identify mechanistic and functional constraints shaping indel diversity. Using poliovirus RdRp variants of distinct fidelity and genome recombination rates, we demonstrate tradeoffs between fidelity and Indel generation. Additionally, we show that maintaining translation frame and viral RNA structures constrain the Indel landscape and that, due to these significant fitness effects, Indels exert a significant deleterious load on adapting viral populations. Conversely, we uncover positively selected Indels that modulate RNA structure, generate protein variants, and produce defective interfering genomes in viral populations. Together, our analyses establish the kinetic and mechanistic tradeoffs between misincorporation, recombination, and Indel rates and reveal functional principles defining the central role of Indels in virus evolution, emergence, and the regulation of viral infection.

Clinical characteristics and host immunity responses of SARS-CoV-2 Omicron variant BA.2 with deletion of ORF7a, ORF7b and ORF8

BACKGROUND

The pathogenicity and virulence of the Omicron strain have weakened significantly pathogenesis of Omicron variants. Accumulating data indicated accessory proteins play crucial roles in host immune evasion and virus pathogenesis of SARS-CoV-2. Therefore, the impact of simultaneous deletion of accessory protein ORF7a, ORF7b and ORF8 on the clinical characteristics and specific immunity in Omicron breakthrough infected patients (BIPs) need to be verified.

METHODS

Herein, plasma cytokines were identified using a commercial Multi-cytokine detection kit. Enzyme-linked immunosorbent assay and pseudovirus neutralization assays were utilized to determine the titers of SARS-CoV-2 specific binding antibodies and neutralizing antibodies, respectively. In addition, an enzyme-linked immunospot assay was used to quantify SARS-CoV-2 specific T cells and memory B cells.

RESULTS

A local COVID-19 outbreak was caused by the Omicron BA.2 variant, which featured a deletion of 871 base pairs (∆871 BA.2), resulting in the removal of ORF7a, ORF7b, and ORF8. We found that hospitalized patients with ∆871 BA.2 had significantly shorter hospital stays than those with wild-type (WT) BA.2. Plasma cytokine levels in both ∆871 BA.2 and WT BA.2 patients were within the normal range of reference, and there was no notable difference in the titers of SARS-CoV-2 ancestor or Omicron-specific binding IgG antibodies, neutralizing antibody titers, effector T cells, and memory B cells frequencies between ∆871 BA.2 and WT BA.2 infected adult patients. However, antibody titers in ∆871 BA.2 infected adolescents were higher than in adults.

CONCLUSIONS

The simultaneous deletion of ORF7a, ORF7b, and ORF8 facilitates the rapid clearance of the BA.2 variant, without impacting cytokine levels or affecting SARS-CoV-2 specific humoral and cellular immunity in Omicron-infected individuals.

Spike substitution T813S increases Sarbecovirus fusogenicity by enhancing the usage of TMPRSS2

SARS-CoV Spike (S) protein shares considerable homology with SARS-CoV-2 S, especially in the conserved S2 subunit (S2). S protein mediates coronavirus receptor binding and membrane fusion, and the latter activity can greatly influence coronavirus infection. We observed that SARS-CoV S is less effective in inducing membrane fusion compared with SARS-CoV-2 S. We identify that S813T mutation is sufficient in S2 interfering with the cleavage of SARS-CoV-2 S by TMPRSS2, reducing spike fusogenicity and pseudoparticle entry. Conversely, the mutation of T813S in SARS-CoV S increased fusion ability and viral replication. Our data suggested that residue 813 in the S was critical for the proteolytic activation, and the change from threonine to Serine at 813 position might be an evolutionary feature adopted by SARS-2-related viruses. This finding deepened the understanding of Spike fusogenicity and could provide a new perspective for exploring Sarbecovirus' evolution.

Tracking of Mutational Signature of SARS-CoV-2 Omicron on Distinct Continents and Little Difference was Found

The Omicron variant is currently ravaging the world, raising serious concern globally. Monitoring genomic variations and determining their influence on biological features are critical for tracing its ongoing transmission and facilitating effective measures. Based on large-scale sequences from different continents, this study found that: (i) The genetic diversity of Omicron is much lower than that of the Delta variant. Still, eight deletions (Del 1-8) and 1 insertion, as well as 130 SNPs, were detected on the Omicron genomes, with two deletions (Del 3 and 4) and 38 SNPs commonly detected on all continents and exhibiting high-occurring frequencies. (ii) Four groups of tightly linked SNPs (linkage I-IV) were detected, among which linkage I, containing 38 SNPs, with 6 located in the RBD, increased its occurring frequency remarkably over time. (iii) The third codons of the Omicron shouldered the most mutation pressures, while the second codons presented the least flexibility. (iv) Four major mutants with amino acid substitutions in the RBD were detected, and further structural analysis suggested that the substitutions did not alter the viral receptor binding ability greatly. It was inferred that though the Omicron genome harbored great changes in antigenicity and remarkable ability to evade immunity, it was immune-pressure selected. This study tracked mutational signatures of Omicron variant and the potential biological significance of the SNPs, and the linkages await further functional verification.

A New Way to Trace SARS-CoV-2 Variants Through Weighted Network Analysis of Frequency Trajectories of Mutations

Early detection of SARS-CoV-2 variants enables timely tracking of clinically important strains in order to inform the public health response. Current subtype-based variant surveillance depending on prior subtype assignment according to lag features and their continuous risk assessment may delay this process. We proposed a weighted network framework to model the frequency trajectories of mutations (FTMs) for SARS-CoV-2 variant tracing, without requiring prior subtype assignment. This framework modularizes the FTMs and conglomerates synchronous FTMs together to represent the variants. It also generates module clusters to unveil the epidemic stages and their contemporaneous variants. Eventually, the module-based variants are assessed by phylogenetic tree through sub-sampling to facilitate communication and control of the epidemic. This process was benchmarked using worldwide GISAID data, which not only demonstrated all the methodology features but also showed the module-based variant identification had highly specific and sensitive mapping with the global phylogenetic tree. When applying this process to regional data like India and South Africa for SARS-CoV-2 variant surveillance, the approach clearly elucidated the national dispersal history of the viral variants and their co-circulation pattern, and provided much earlier warning of Beta (B.1.351), Delta (B.1.617.2), and Omicron (B.1.1.529). In summary, our work showed that the weighted network modeling of FTMs enables us to rapidly and easily track down SARS-CoV-2 variants overcoming prior viral subtyping with lag features, accelerating the understanding and surveillance of COVID-19.