Host Immune Response Driving SARS-CoV-2 Evolution

RNA editing sites and triplet usage in exomes of bat RNA virus genomes of the family Paramyxoviridae

Bats contain a diverse spectrum of viral species in their bodies. The RNA virus family Paramyxoviridae tends to infect several vertebrate species, which are accountable for a variety of devastating infections in both humans and animals. Viruses of this kind include measles, mumps, and Hendra. Some synonymous codons are favoured over others in mRNAs during gene-to-protein synthesis process. Such phenomenon is termed as codon usage bias (CUB). Our research emphasized many aspects that shape the CUB of genes in the Paramyxoviridae family found in bats. Here, the nitrogenous base A occurred the most. AT was found to be abundant in the coding sequences of the Paramyxoviridae family. RSCU data revealed that A or T ending codons occurred more frequently than predicted. Furthermore, 3 overrepresented codons (CAT, AGA, and GCA) and 7 underrepresented codons (CCG, TCG, CGC, CGG, CGT, GCG and ACG) were detected in the viral genomes. Correspondence analysis, neutrality plot, and parity plots highlight the combined impact of mutational pressure and natural selection on CUB. The neutrality plot of GC12 against GC3 yielded a regression coefficient value of 0.366, indicating that natural selection had a significant (63.4 %) impact. Moreover, RNA editing analysis was done, which revealed the highest frequency of C to T mutations. The results of our research revealed the pattern of codon usage and RNA editing sites in Paramyxoviridae genomes.

Global lineage evolution pattern of sars-cov-2 in Africa, America, Europe, and Asia: A comparative analysis of variant clusters and their relevance across continents

Objective

The objective of this study is to provide a comparative analysis of variant clusters and their relevance across Africa, America, Europe, and Asia, in order to understand the evolutionary patterns of the virus across different regions and to inform the development of targeted interventions and genomic surveillance eforts.

Methods

The study analyzed the global lineage evolution pattern of 74, 075 severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) genomes from 32 countries across four continents, focusing on variant clusters and their relevance across regions. Variants were weighted according to their hierarchical level. The correlation between variants was visualized through Dimensionality reduction analysis and Pairwise Pearson's correlation. We presented a reconstructed phylogenetic tree based on correlation analysis and variant weights.

Results

The analysis revealed that each continent had distinct variant clusters and different evolutionary patterns. The Americas had two clustered variants before lineage divergence and a downstream confluence lineage, Europe had bifurcation into two global lineages with an early occurrence of certain cluster while Asia had a downstream confluence of two large lineages diverging by two distinct clusters. Based on the cluster patterns of shared variants of the SARS-CoV-2 virus, Africa demonstrated a relatively clear distinction among three distinct regions.

Conclusions

The study provides insights into the evolutionary patterns of SARS-CoV-2 and highlights the importance of international collaboration in tracking and responding to emerging variants. The study found that the global pandemic was driven by Omicron variants that evolved with significant differences between countries and regions, and with different patterns across continents.

Intra-host mutation rate of acute SARS-CoV-2 infection during the initial pandemic wave

SARS-CoV-2 mutation is minimized through a proofreading function encoded by NSP-14. Most estimates of the SARS-CoV-2 mutation rate are derived from population based sequence data. Our understanding of SARS-CoV-2 evolution might be enhanced through analysis of intra-host viral mutation rates in specific populations. Viral genome analysis was performed between paired samples and mutations quantified at allele frequencies (AF) ≥ 0.25, ≥ 0.5 and ≥ 0.75. Mutation rate was determined employing F81 and JC69 evolution models and compared between isolates with (ΔNSP-14) and without (wtNSP-14) non-synonymous mutations in NSP-14 and by patient comorbidity. Forty paired samples with median interval of 13 days [IQR 8.5-20] were analyzed. The estimated mutation rate by F81 modeling was 93.6 (95%CI 90.8-96.4], 40.7 (95%CI 38.9-42.6) and 34.7 (95%CI 33.0-36.4) substitutions/genome/year at AF ≥ 0.25, ≥ 0.5, ≥ 0.75 respectively. Mutation rate in ΔNSP-14 were significantly elevated at AF ≥ 0.25 vs wtNSP-14. Patients with immune comorbidities had higher mutation rate at all allele frequencies. Intra-host SARS-CoV-2 mutation rates are substantially higher than those reported through population analysis. Virus strains with altered NSP-14 have accelerated mutation rate at low AF. Immunosuppressed patients have elevated mutation rate at all AF. Understanding intra-host virus evolution will aid in current and future pandemic modeling.

A tale of two waves: Delineating diverse genomic and transmission landscapes driving the COVID-19 pandemic in Pune, India

BACKGROUND

Modern response to pandemics, critical for effective public health measures, is shaped by the availability and integration of diverse epidemiological outbreak data. Tracking variants of concern (VOC) is integral to understanding the evolution of SARS-CoV-2 in space and time, both at the local level and global context. This potentially generates actionable information when integrated with epidemiological outbreak data.

METHODS

A city-wide network of researchers, clinicians, and pathology diagnostic laboratories was formed for genome surveillance of COVID-19 in Pune, India. The genomic landscapes of 10,496 sequenced samples of SARS-CoV-2 driving peaks of infection in Pune between December-2020 to March-2022, were determined. As a modern response to the pandemic, a "band of five" outbreak data analytics approach was used. This integrated the genomic data (Band 1) of the virus through molecular phylogenetics with key outbreak data including sample collection dates and case numbers (Band 2), demographics like age and gender (Band 3-4), and geospatial mapping (Band 5).

RESULTS

The transmission dynamics of VOCs in 10,496 sequenced samples identified B.1.617.2 (Delta) and BA(x) (Omicron formerly known as B.1.1.529) variants as drivers of the second and third peaks of infection in Pune. Spike Protein mutational profiling during pre and post-Omicron VOCs indicated differential rank ordering of high-frequency mutations in specific domains that increased the charge and binding properties of the protein. Time-resolved phylogenetic analysis of Omicron sub-lineages identified a highly divergent BA.1 from Pune in addition to recombinant X lineages, XZ, XQ, and XM.

CONCLUSIONS

The band of five outbreak data analytics approach, which integrates five different types of data, highlights the importance of a strong surveillance system with high-quality meta-data for understanding the spatiotemporal evolution of the SARS-CoV-2 genome in Pune. These findings have important implications for pandemic preparedness and could be critical tools for understanding and responding to future outbreaks.

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.

Evolution of SARS-CoV-2 Variants: Implications on Immune Escape, Vaccination, Therapeutic and Diagnostic Strategies

The COVID-19 pandemic caused by SARS-CoV-2 is associated with a lower fatality rate than its SARS and MERS counterparts. However, the rapid evolution of SARS-CoV-2 has given rise to multiple variants with varying pathogenicity and transmissibility, such as the Delta and Omicron variants. Individuals with advanced age or underlying comorbidities, including hypertension, diabetes and cardiovascular diseases, are at a higher risk of increased disease severity. Hence, this has resulted in an urgent need for the development of better therapeutic and preventive approaches. This review describes the origin and evolution of human coronaviruses, particularly SARS-CoV-2 and its variants as well as sub-variants. Risk factors that contribute to disease severity and the implications of co-infections are also considered. In addition, various antiviral strategies against COVID-19, including novel and repurposed antiviral drugs targeting viral and host proteins, as well as immunotherapeutic strategies, are discussed. We critically evaluate strategies of current and emerging vaccines against SARS-CoV-2 and their efficacy, including immune evasion by new variants and sub-variants. The impact of SARS-CoV-2 evolution on COVID-19 diagnostic testing is also examined. Collectively, global research and public health authorities, along with all sectors of society, need to better prepare against upcoming variants and future coronavirus outbreaks.

Emerging Dominant SARS-CoV-2 Variants

Accurate and reliable forecasting of emerging dominant severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) variants enables policymakers and vaccine makers to get prepared for future waves of infections. The last three waves of SARS-CoV-2 infections caused by dominant variants, Omicron (BA.1), BA.2, and BA.4/BA.5, were accurately foretold by our artificial intelligence (AI) models built with biophysics, genotyping of viral genomes, experimental data, algebraic topology, and deep learning. On the basis of newly available experimental data, we analyzed the impacts of all possible viral spike (S) protein receptor-binding domain (RBD) mutations on the SARS-CoV-2 infectivity. Our analysis sheds light on viral evolutionary mechanisms, i.e., natural selection through infectivity strengthening and antibody resistance. We forecast that BP.1, BL*, BA.2.75*, BQ.1*, and particularly BN.1* have a high potential to become the new dominant variants to drive the next surge. Our key projection about these variants dominance made on Oct. 18, 2022 (see arXiv:2210.09485) became reality in late November 2022.

The clinical characteristics of pediatric patients infected by SARS-CoV-2 Omicron variant and whole viral genome sequencing analysis

Pediatric population was generally less affected clinically by SARS-CoV-2 infection. Few pediatric cases of COVID-19 have been reported compared to those reported in infected adults. However, a rapid increase in the hospitalization rate of SARS-CoV-2 infected pediatric patients was observed during Omicron variant dominated COVID-19 outbreak. In this study, we analyzed the B.1.1.529 (Omicron) genome sequences collected from pediatric patients by whole viral genome amplicon sequencing using Illumina next generation sequencing platform, followed by phylogenetic analysis. The demographic, epidemiologic and clinical data of these pediatric patients are also reported in this study. Fever, cough, running nose, sore throat and vomiting were the more commonly reported symptoms in children infected by Omicron variant. A novel frameshift mutation was found in the ORF1b region (NSP12) of the genome of Omicron variant. Seven mutations were identified in the target regions of the WHO listed SARS-CoV-2 primers and probes. On protein level, eighty-three amino acid substitutions and fifteen amino acid deletions were identified. Our results indicate that asymptomatic infection and transmission among children infected by Omicron subvariants BA.2.2 and BA.2.10.1 are not common. Omicron may have different pathogenesis in pediatric population.

Mutations in SARS-CoV-2: Insights on structure, variants, vaccines, and biomedical interventions

COVID-19 is a worldwide pandemic caused by SARS-coronavirus-2 (SARS-CoV-2). Less than a year after the emergence of the Covid-19 pandemic, many vaccines have arrived on the market with innovative technologies in the field of vaccinology. Based on the use of messenger RNA (mRNA) encoding the Spike SARS-Cov-2 protein or on the use of recombinant adenovirus vectors enabling the gene encoding the Spike protein to be introduced into our cells, these strategies make it possible to envisage the vaccination in a new light with tools that are more scalable than the vaccine strategies used so far. Faced with the appearance of new variants, which will gradually take precedence over the strain at the origin of the pandemic, these new strategies will allow a much faster update of vaccines to fight against these new variants, some of which may escape neutralization by vaccine antibodies. However, only a vaccination policy based on rapid and massive vaccination of the population but requiring a supply of sufficient doses could make it possible to combat the emergence of these variants. Indeed, the greater the number of infected individuals, the faster the virus multiplies, with an increased risk of the emergence of variants in these RNA viruses. This review will discuss SARS-CoV-2 pathophysiology and evolution approaches in altered transmission platforms and emphasize the different mutations and how they influence the virus characteristics. Also, this article summarizes the common vaccines and the implication of the mutations and genetic variety of SARS-CoV-2 on the COVID-19 biomedical arbitrations.

Advances in Next-Generation Coronavirus Vaccines in Response to Future Virus Evolution

Coronavirus disease (COVID-19), caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), has spread to more than 230 countries and territories worldwide since its outbreak in late 2019. In less than three years, infection by SARS-CoV-2 has resulted in over 600 million cases of COVID-19 and over 6.4 million deaths. Vaccines have been developed with unimaginable speed, and 11 have already been approved by the World Health Organization and given Emergency Use Listing. The administration of several first-generation SARS-CoV-2 vaccines has successfully decelerated the spread of COVID-19 but not stopped it completely. In the ongoing fight against viruses, genetic mutations frequently occur in the viral genome, resulting in a decrease in vaccine-induced antibody neutralization and widespread breakthrough infection. Facing the evolution and uncertainty of SARS-CoV-2 in the future, and the possibility of the spillover of other coronaviruses to humans, the need for vaccines with a broad spectrum of antiviral variants against multiple coronaviruses is recognized. It is imperative to develop a universal coronavirus or pan-coronavirus vaccine or drug to combat the ongoing COVID-19 pandemic as well as to prevent the next coronavirus pandemic. In this review, in addition to summarizing the protective effect of approved vaccines, we systematically summarize current work on the development of vaccines aimed at suppressing multiple SARS-CoV-2 variants of concern as well as multiple coronaviruses.

Prediction of Recurrent Mutations in SARS-CoV-2 Using Artificial Neural Networks

Predicting SARS-CoV-2 mutations is difficult, but predicting recurrent mutations driven by the host, such as those caused by host deaminases, is feasible. We used machine learning to predict which positions from the SARS-CoV-2 genome will hold a recurrent mutation and which mutations will be the most recurrent. We used data from April 2021 that we separated into three sets: a training set, a validation set, and an independent test set. For the test set, we obtained a specificity value of 0.69, a sensitivity value of 0.79, and an Area Under the Curve (AUC) of 0.8, showing that the prediction of recurrent SARS-CoV-2 mutations is feasible. Subsequently, we compared our predictions with updated data from January 2022, showing that some of the false positives in our prediction model become true positives later on. The most important variables detected by the model's Shapley Additive exPlanation (SHAP) are the nucleotide that mutates and RNA reactivity. This is consistent with the SARS-CoV-2 mutational bias pattern and the preference of some host deaminases for specific sequences and RNA secondary structures. We extend our investigation by analyzing the mutations from the variants of concern Alpha, Beta, Delta, Gamma, and Omicron. Finally, we analyzed amino acid changes by looking at the predicted recurrent mutations in the M-pro and spike proteins.

In-depth investigation of the point mutation pattern of HIV-1

Mutations may produce highly transmissible and damaging HIV variants, which increase the genetic diversity, and pose a challenge to develop vaccines. Therefore, it is of great significance to understand how mutations drive the virulence of HIV. Based on the 11897 reliable genomes of HIV-1 retrieved from HIV sequence Database, we analyze the 12 types of point mutation (A>C, A>G, A>T, C>A, C>G, C>T, G>A, G>C, G>T, T>A, T>C, T>G) from multiple statistical perspectives for the first time. The global/geographical location/subtype/k-mer analysis results report that A>G, G>A, C>T and T>C account for nearly 64% among all SNPs, which suggest that APOBEC-editing and ADAR-editing may play an important role in HIV-1 infectivity. Time analysis shows that most genomes with abnormal mutation numbers comes from African countries. Finally, we use natural vector method to check the k-mer distribution changing patterns in the genome, and find that there is an important substitution pattern between nucleotides A and G, and 2-mer CG may have a significant impact on viral infectivity. This paper provides an insight into the single mutation of HIV-1 by using the latest data in the HIV sequence Database.

Emerging dominant SARS-CoV-2 variants

Accurate and reliable forecasting of emerging dominant severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) variants enables policymakers and vaccine makers to get prepared for future waves of infections. The last three waves of SARS-CoV-2 infections caused by dominant variants Omicron (BA.1), BA.2, and BA.4/BA.5 were accurately foretold by our artificial intelligence (AI) models built with biophysics, genotyping of viral genomes, experimental data, algebraic topology, and deep learning. Based on newly available experimental data, we analyzed the impacts of all possible viral spike (S) protein receptor-binding domain (RBD) mutations on the SARS-CoV-2 infectivity. Our analysis sheds light on viral evolutionary mechanisms, i.e., natural selection through infectivity strengthening and antibody resistance. We forecast that BA.2.10.4, BA.2.75, BQ.1.1, and particularly, BA.2.75+R346T, have high potential to become new dominant variants to drive the next surge.

The roles of APOBEC-mediated RNA editing in SARS-CoV-2 mutations, replication and fitness

During COVID-19 pandemic, mutations of SARS-CoV-2 produce new strains that can be more infectious or evade vaccines. Viral RNA mutations can arise from misincorporation by RNA-polymerases and modification by host factors. Analysis of SARS-CoV-2 sequence from patients showed a strong bias toward C-to-U mutation, suggesting a potential mutational role by host APOBEC cytosine deaminases that possess broad anti-viral activity. We report the first experimental evidence demonstrating that APOBEC3A, APOBEC1, and APOBEC3G can edit on specific sites of SARS-CoV-2 RNA to produce C-to-U mutations. However, SARS-CoV-2 replication and viral progeny production in Caco-2 cells are not inhibited by the expression of these APOBECs. Instead, expression of wild-type APOBEC3 greatly promotes viral replication/propagation, suggesting that SARS-CoV-2 utilizes the APOBEC-mediated mutations for fitness and evolution. Unlike the random mutations, this study suggests the predictability of all possible viral genome mutations by these APOBECs based on the UC/AC motifs and the viral genomic RNA structure.

(Epi)transcriptomics in cardiovascular and neurological complications of COVID-19

Although systemic inflammation and pulmonary complications increase the mortality rate in COVID-19, a broad spectrum of cardiovascular and neurological complications can also contribute to significant morbidity and mortality. The molecular mechanisms underlying cardiovascular and neurological complications during and after SARS-CoV-2 infection are incompletely understood. Recently reported perturbations of the epitranscriptome of COVID-19 patients indicate that mechanisms including those derived from RNA modifications and non-coding RNAs may play a contributing role in the pathogenesis of COVID-19. In this review paper, we gathered recently published studies investigating (epi)transcriptomic fluctuations upon SARS-CoV-2 infection, focusing on the brain-heart axis since neurological and cardiovascular events and their sequelae are of utmost prevalence and importance in this disease.

Phylogenetic and amino acid signature analysis of the SARS-CoV-2s lineages circulating in Tunisia

Since the beginning of the Coronavirus disease-2019 pandemic, there has been a growing interest in exploring SARS-CoV-2 genetic variation to understand the origin and spread of the pandemic, improve diagnostic methods and develop the appropriate vaccines. The objective of this study was to identify the SARS-CoV-2s lineages circulating in Tunisia and to explore their amino acid signature in order to follow their genome dynamics. Whole genome sequencing and genetic analyses of fifty-eight SARS-CoV-2 samples collected during one-year between March 2020 and March 2021 from the National Influenza Center were performed using three sampling strategies.. Multiple lineage introductions were noted during the initial phase of the pandemic, including B.4, B.1.1, B.1.428.2, B.1.540 and B.1.1.189. Subsequently, lineages B1.160 (24.2%) and B1.177 (22.4%) were dominant throughout the year. The Alpha variant (B.1.1.7 lineage) was identified in February 2021 and firstly observed in the center of our country. In addition, A clear diversity of lineages was observed in the North of the country. A total of 335 mutations including 10 deletions were found. The SARS-CoV-2 proteins ORF1ab, Spike, ORF3a, and Nucleocapsid were observed as mutation hotspots with a mutation frequency exceeding 20%. The 2 most frequent mutations, D614G in S protein and P314L in Nsp12 appeared simultaneously and are often associated with increased viral infectivity. Interestingly, deletions in coding regions causing consequent deletions of amino acids and frame shifts were identified in NSP3, NSP6, S, E, ORF7a, ORF8 and N proteins.

These findings contribute to define the COVID-19 outbreak in Tunisia. Despite the country's limited resources, surveillance of SARS-CoV-2 genomic variation should be continued to control the occurrence of new variants.

Methodology-Centered Review of Molecular Modeling, Simulation, and Prediction of SARS-CoV-2

Despite tremendous efforts in the past two years, our understanding of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), virus-host interactions, immune response, virulence, transmission, and evolution is still very limited. This limitation calls for further in-depth investigation. Computational studies have become an indispensable component in combating coronavirus disease 2019 (COVID-19) due to their low cost, their efficiency, and the fact that they are free from safety and ethical constraints. Additionally, the mechanism that governs the global evolution and transmission of SARS-CoV-2 cannot be revealed from individual experiments and was discovered by integrating genotyping of massive viral sequences, biophysical modeling of protein-protein interactions, deep mutational data, deep learning, and advanced mathematics. There exists a tsunami of literature on the molecular modeling, simulations, and predictions of SARS-CoV-2 and related developments of drugs, vaccines, antibodies, and diagnostics. To provide readers with a quick update about this literature, we present a comprehensive and systematic methodology-centered review. Aspects such as molecular biophysics, bioinformatics, cheminformatics, machine learning, and mathematics are discussed. This review will be beneficial to researchers who are looking for ways to contribute to SARS-CoV-2 studies and those who are interested in the status of the field.

SARS-CoV-2: An Overview of the Genetic Profile and Vaccine Effectiveness of the Five Variants of Concern

With the onset of the COVID-19 pandemic, enormous efforts have been made to understand the genus SARS-CoV-2. Due to the high rate of global transmission, mutations in the viral genome were inevitable. A full understanding of the viral genome and its possible changes represents one of the crucial aspects of pandemic management. Structural protein S plays an important role in the pathogenicity of SARS-CoV-2, mutations occurring at this level leading to viral forms with increased affinity for ACE2 receptors, higher transmissibility and infectivity, resistance to neutralizing antibodies and immune escape, increasing the risk of infection and disease severity. Thus, five variants of concern are currently being discussed, Alpha, Beta, Gamma, Delta and Omicron. In the present review, a comprehensive summary of the following critical aspects regarding SARS-CoV-2 has been made: (i) the genomic characteristics of SARS-CoV-2; (ii) the pathological mechanism of transmission, penetration into the cell and action on specific receptors; (iii) mutations in the SARS-CoV-2 genome; and (iv) possible implications of mutations in diagnosis, treatment, and vaccination.

Tracing the SARS-CoV-2 infection on the ocular surface: Overview and preliminary corneoscleral transcriptome sequencing

COVID-19's impact on the ocular surface has already been recognized, however the molecular mechanisms induced by the infection on the ocular surface are still unclear. The aim of this paper is to provide a first overview of the transcriptional perturbations caused by SARS-CoV-2 on the ocular surface by analyzing gene expression profile of corneoscleral ring samples from post-mortem SARS-CoV-2 positive donors (PD). The presence of SARS-CoV-2 on the ocular surface, in tears and corneal tissues has rarely been detected in infected individuals in both the presence and the absence of ocular manifestations. In this preliminary study, 6 human corneoscleral tissues of 3 PD and two tissues from a negative donor (CTRL) were obtained at the local eye bank. The presence of genomic and sub-genomic SARS-CoV-2 RNAs was assessed by qRT-PCR, while transcriptome analysis (RNA-sequencing) was performed by Illumina. Principal Component Analysis (PCA), search for differentially expressed genes (DEGs) and Gene Ontology (GO)-enrichment analysis were performed. Three samples from PD were found positive for SARS-CoV-2 genomic RNA, although the absence of sub-genomic RNAs indicated an inactive virus. PCA analysis grouped 3 different clusters, one including CTRL, and the other two including, respectively, PD with undetected SARS-CoV-2 (PD-SARS-neg) and PD with detected SARS-CoV-2 (PD-SARS-pos). The DEGs in common with the 2 PD clusters included several genes associable to the interferon pathway, such as ADAMTS4, RSAD2, MMP1, IL6, ISG15 and proinflammatory cytokines. Among the down-regulated genes we found AQP5. GO analysis revealed 77 GO terms over-represented in PD-SARS-neg vs. CTRL, and 17 GO terms in PD-SARS-pos vs. CTRL. The presence of SARS-CoV-2 RNA and RNA-sequencing reads in ocular surface tissues supports the possibility that the eye acts as an entry route. The modulation of early responsive genes, together with several ISGs suggests a potential protective responsiveness of the ocular tissues to SARS-CoV-2.

Virus-specific editing identification approach reveals the landscape of A-to-I editing and its impacts on SARS-CoV-2 characteristics and evolution

Upon SARS-CoV-2 infection, viral intermediates specifically activate the IFN response through MDA5-mediated sensing and accordingly induce ADAR1 p150 expression, which might lead to viral A-to-I RNA editing. Here, we developed an RNA virus-specific editing identification pipeline, surveyed 7622 RNA-seq data from diverse types of samples infected with SARS-CoV-2, and constructed an atlas of A-to-I RNA editing sites in SARS-CoV-2. We found that A-to-I editing was dynamically regulated, varied between tissue and cell types, and was correlated with the intensity of innate immune response. On average, 91 editing events were deposited at viral dsRNA intermediates per sample. Moreover, editing hotspots were observed, including recoding sites in the spike gene that affect viral infectivity and antigenicity. Finally, we provided evidence that RNA editing accelerated SARS-CoV-2 evolution in humans during the epidemic. Our study highlights the ability of SARS-CoV-2 to hijack components of the host antiviral machinery to edit its genome and fuel its evolution, and also provides a framework and resource for studying viral RNA editing.

Emerging Vaccine-Breakthrough SARS-CoV-2 Variants

The surge of COVID-19 infections has been fueled by new SARS-CoV-2 variants, namely Alpha, Beta, Gamma, Delta, and so forth. The molecular mechanism underlying such surge is elusive due to the existence of 28 554 unique mutations, including 4 653 non-degenerate mutations on the spike protein. Understanding the molecular mechanism of SARS-CoV-2 transmission and evolution is a prerequisite to foresee the trend of emerging vaccine-breakthrough variants and the design of mutation-proof vaccines and monoclonal antibodies. We integrate the genotyping of 1 489 884 SARS-CoV-2 genomes, a library of 130 human antibodies, tens of thousands of mutational data, topological data analysis, and deep learning to reveal SARS-CoV-2 evolution mechanism and forecast emerging vaccine-breakthrough variants. We show that prevailing variants can be quantitatively explained by infectivity-strengthening and vaccine-escape (co-)mutations on the spike protein RBD due to natural selection and/or vaccination-induced evolutionary pressure. We illustrate that infectivity strengthening mutations were the main mechanism for viral evolution, while vaccine-escape mutations become a dominating viral evolutionary mechanism among highly vaccinated populations. We demonstrate that Lambda is as infectious as Delta but is more vaccine-resistant. We analyze emerging vaccine-breakthrough comutations in highly vaccinated countries, including the United Kingdom, the United States, Denmark, and so forth. Finally, we identify sets of comutations that have a high likelihood of massive growth: [A411S, L452R, T478K], [L452R, T478K, N501Y], [V401L, L452R, T478K], [K417N, L452R, T478K], [L452R, T478K, E484K, N501Y], and [P384L, K417N, E484K, N501Y]. We predict they can escape existing vaccines. We foresee an urgent need to develop new virus combating strategies.

Functions and consequences of AID/APOBEC-mediated DNA and RNA deamination

The AID/APOBEC polynucleotide cytidine deaminases have historically been classified as either DNA mutators or RNA editors based on their first identified nucleic acid substrate preference. DNA mutators can generate functional diversity at antibody genes but also cause genomic instability in cancer. RNA editors can generate informational diversity in the transcriptome of innate immune cells, and of cancer cells. Members of both classes can act as antiviral restriction factors. Recent structural work has illuminated differences and similarities between AID/APOBEC enzymes that can catalyse DNA mutation, RNA editing or both, suggesting that the strict functional classification of members of this family should be reconsidered. As many of these enzymes have been employed for targeted genome (or transcriptome) editing, a more holistic understanding will help improve the design of therapeutically relevant programmable base editors.

In this Perspective, Pecori et al. provide an overview of the AID/APOBEC cytidine deaminase family, discussing key structural features, how they contribute to viral and tumour evolution and how they can be harnessed for (potentially therapeutic) base-editing purposes.

Rapid SARS-CoV-2 Intra-Host and Within-Household Emergence of Novel Haplotypes

In February 2020, the municipality of Vo', a small town near Padua (Italy) was quarantined due to the first coronavirus disease 19 (COVID-19)-related death detected in Italy. To investigate the viral prevalence and clinical features, the entire population was swab tested in two sequential surveys. Here we report the analysis of 87 viral genomes, which revealed that the unique ancestor haplotype introduced in Vo' belongs to lineage B, carrying the mutations G11083T and G26144T. The viral sequences allowed us to investigate the viral evolution while being transmitted within and across households and the effectiveness of the non-pharmaceutical interventions implemented in Vo'. We report, for the first time, evidence that novel viral haplotypes can naturally arise intra-host within an interval as short as two weeks, in approximately 30% of the infected individuals, regardless of symptom severity or immune system deficiencies. Moreover, both phylogenetic and minimum spanning network analyses converge on the hypothesis that the viral sequences evolved from a unique common ancestor haplotype that was carried by an index case. The lockdown extinguished both the viral spread and the emergence of new variants.

Spatio-temporal dynamics of intra-host variability in SARS-CoV-2 genomes

During the course of the COVID-19 pandemic, large-scale genome sequencing of SARS-CoV-2 has been useful in tracking its spread and in identifying variants of concern (VOC). Viral and host factors could contribute to variability within a host that can be captured in next-generation sequencing reads as intra-host single nucleotide variations (iSNVs). Analysing 1347 samples collected till June 2020, we recorded 16 410 iSNV sites throughout the SARS-CoV-2 genome. We found ∼42% of the iSNV sites to be reported as SNVs by 30 September 2020 in consensus sequences submitted to GISAID, which increased to ∼80% by 30th June 2021. Following this, analysis of another set of 1774 samples sequenced in India between November 2020 and May 2021 revealed that majority of the Delta (B.1.617.2) and Kappa (B.1.617.1) lineage-defining variations appeared as iSNVs before getting fixed in the population. Besides, mutations in RdRp as well as RNA-editing by APOBEC and ADAR deaminases seem to contribute to the differential prevalence of iSNVs in hosts. We also observe hyper-variability at functionally critical residues in Spike protein that could alter the antigenicity and may contribute to immune escape. Thus, tracking and functional annotation of iSNVs in ongoing genome surveillance programs could be important for early identification of potential variants of concern and actionable interventions.

Identification of genotypic variants and its proteomic mutations of Brazilian SARS-CoV-2 isolates

The second wave of COVID-19 caused by severe acute respiratory syndrome virus (SARS-CoV-2) is rapidly spreading over the world. Mechanisms behind the flee from current antivirals are still unclear due to the continuous occurrence of SARS-CoV-2 genetic variants. Brazil is the world's second-most COVID-19 affected country. In the present study, we identified the genomic and proteomic variants of Brazilian SARS-CoV-2 isolates. We identified 16 different genotypic variants were found among the 27 isolates. The genotypes of three isolates such as Bra/1236/2021 (G15), Bra/MASP2C844R2/2020 (G11), and Bra/RJ-DCVN5/2020 (G9) have a unique mutant in NSP4 (S184N), 2′O-Mutase (R216N), membrane protein (A2V) and Envelope protein (V5A). A mutation in RdRp of SARS-CoV-2, particularly the change of Pro-to Leu-at 323 resulted in the stabilization of the structure in BRA/CD1739-P4/2020. NSP4, NSP5 protein mutants are more virulent in genotype 15 and 16. A fast protein folding rate changes the structural stability and leads to escape for current antivirals. Thus, our findings help researchers to develop the best potent antivirals based on the new mutant of Brazilian isolates.

Evolutionary trend of SARS-CoV-2 inferred by the homopolymeric nucleotide repeats

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is the causative agent of the current global COVID-19 pandemic, in which millions of lives have been lost. Understanding the zoonotic evolution of the coronavirus may provide insights for developing effective vaccines, monitoring the transmission trends, and preventing new zoonotic infections. Homopolymeric nucleotide repeats (HP), the most simple tandem repeats, are a ubiquitous feature of eukaryotic genomes. Yet the HP distributions and roles in coronavirus genome evolution are poorly investigated. In this study, we characterize the HP distributions and trends in the genomes of bat and human coronaviruses and SARS-CoV-2 variants. The results show that the SARS-CoV-2 genome is abundant in HPs, and has augmented HP contents during evolution. Especially, the disparity of HP poly-(A/T) and ploy-(C/G) of coronaviruses increases during the evolution in human hosts. The disparity of HP poly-(A/T) and ploy-(C/G) is correlated to host adaptation and the virulence level of the coronaviruses. Therefore, we propose that the HP disparity can be a quantitative measure for the zoonotic evolution levels of coronaviruses. Peculiarly, the HP disparity measure infers that SARS-CoV-2 Omicron variants have a high disparity of HP poly-(A/T) and ploy-(C/G), suggesting a high adaption to the human hosts.

The Roles of APOBEC-mediated RNA Editing in SARS-CoV-2 Mutations, Replication and Fitness.. [PMID: 34981048 PMCID: PMC8722585 DOI: 10.1101/2021.12.18.473309]

During COVID-19 pandemic, mutations of SARS-CoV-2 produce new strains that can be more infectious or evade vaccines. Viral RNA mutations can arise from misincorporation by RNA-polymerases and modification by host factors. Analysis of SARS-CoV-2 sequence from patients showed a strong bias toward C-to-U mutation, suggesting a potential mutational role by host APOBEC cytosine deaminases that possess broad anti-viral activity. We report the first experimental evidence demonstrating that APOBEC3A, APOBEC1, and APOBEC3G can edit on specific sites of SARS-CoV-2 RNA to produce C-to-U mutations. However, SARS-CoV-2 replication and viral progeny production in Caco-2 cells are not inhibited by the expression of these APOBECs. Instead, expression of wild-type APOBEC3 greatly promotes viral replication/propagation, suggesting that SARS-CoV-2 utilizes the APOBEC-mediated mutations for fitness and evolution. Unlike the random mutations, this study suggests the predictability of all possible viral genome mutations by these APOBECs based on the UC/AC motifs and the viral genomic RNA structure.

Efficient Editing of SARS-CoV-2 genomic RNA by Host APOBEC deaminases and Its Potential Impacts on the Viral Replication and Emergence of New Strains in COVID-19 Pandemic

Whole genome amplicon sequencing and phylogenetic analysis of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) from lineage B.1.36.27 isolated in Hong Kong

BACKGROUND

The import of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) lineage B.1.36.27 has sparked the fourth wave of COVID-19 outbreak in Hong Kong. This strain has been circulating in Hong Kong since September 2020 but rarely found in other countries (<1%).

RESEARCH DESIGN AND METHODS

A total of 14 SARS-CoV-2 genome sequences collected from patients in Hong Kong between July 2020 and March 2021 were determined by whole viral genome sequencing using Illumina next-generation sequencing platform, followed by phylogenetic analysis.

RESULTS

Of the 14 SARS-CoV-2 genome sequences analyzed, 9 strains belonged to the PANGO lineage B.1.36.27, GISAID clade GH, and Nextclade clade 20A. Compared to the reference genome, 31 nucleotide differences and 11 amino acid differences were identified in the genome of the SARS-CoV-2 from PANGO lineage B.1.36.27.

CONCLUSIONS

We reported the nucleotides and amino acids mutations identified in the SARS-CoV-2 from PANGO lineage B.1.36.27. Our viral genome sequences enriched the understanding of SARS-CoV-2 mutational landscape and improved the repertoire of known SARS-CoV-2 variants for tracking and tracing. From this study, we found no evidence to show that SARS-CoV-2 from lineage B.1.36.27 can compromise existing vaccines and antibody therapies.

Mechanisms of SARS-CoV-2 Evolution Revealing Vaccine-Resistant Mutations in Europe and America

The importance of understanding SARS-CoV-2 evolution cannot be overlooked. Recent studies confirm that natural selection is the dominating mechanism of SARS-CoV-2 evolution, which favors mutations that strengthen viral infectivity. Here, we demonstrate that vaccine-breakthrough or antibody-resistant mutations provide a new mechanism of viral evolution. Specifically, vaccine-resistant mutation Y449S in the spike (S) protein receptor-binding domain, which occurred in co-mutations Y449S and N501Y, has reduced infectivity compared to that of the original SARS-CoV-2 but can disrupt existing antibodies that neutralize the virus. By tracking the evolutionary trajectories of vaccine-resistant mutations in more than 2.2 million SARS-CoV-2 genomes, we reveal that the occurrence and frequency of vaccine-resistant mutations correlate strongly with the vaccination rates in Europe and America. We anticipate that as a complementary transmission pathway, vaccine-breakthrough or antibody-resistant mutations, like those in Omicron, will become a dominating mechanism of SARS-CoV-2 evolution when most of the world's population is either vaccinated or infected. Our study sheds light on SARS-CoV-2 evolution and transmission and enables the design of the next-generation mutation-proof vaccines and antibody drugs.

Growth, reproduction numbers and factors affecting the spread of SARS-CoV-2 novel variants of concern in the UK from October 2020 to July 2021: a modelling analysis

OBJECTIVES

Importations of novel variants of concern (VOC), particularly B.1.617.2, have become the impetus behind recent outbreaks of SARS-CoV-2. Concerns around the impact on vaccine effectiveness, transmissibility and severity are now driving the public health response to these variants. This paper analyses the patterns of growth in hospitalisations and confirmed cases for novel VOCs by age groups, geography and ethnicity in the context of changing behaviour, non-pharmaceutical interventions (NPIs) and the UK vaccination programme. We seek to highlight where strategies have been effective and periods that have facilitated the establishment of new variants.

DESIGN

We have algorithmically linked the most complete testing and hospitalisation data in England to create a data set of confirmed infections and hospitalisations by SARS-CoV-2 genomic variant. We have used these linked data sets to analyse temporal, geographic and demographic distinctions.

SETTING AND PARTICIPANTS

The setting is England from October 2020 to July 2021. Participants included all COVID-19 tests that included RT-PCR CT gene target data or underwent sequencing and hospitalisations that could be linked to these tests.

METHODS

To calculate the instantaneous growth rate for VOCs we have developed a generalised additive model fit to multiple splines and varying day of the week effects. We have further modelled the instantaneous reproduction number Rt for the B.1.1.7 and B.1.617.2 variants and included a doubly interval censored model to temporally adjust the confirmed variant cases.

RESULTS

We observed a clear replacement of the predominant B.1.1.7 by the B.1.617.2 variant without observing sustained exponential growth in other novel variants. Modelled exponential growth of RT PCR gene target triple-positive cases was initially detected in the youngest age groups, although we now observe across all ages a very small doubling time of 10.7 (95% CI 9.1 to 13.2) days and 8 (95% CI 6.9 to 9.1) days for cases and hospitalisations, respectively. We observe that growth in RT PCR gene target triple-positive cases was first detected in the Indian ethnicity group in late February, with a peak of 0.06 (95% CI 0.07 to 0.05) in the instantaneous growth rate, but is now maintained by the white ethnicity groups, observing a doubling time of 6.8 (95% CI 4.9 to 11) days. Rt analysis indicates a reproduction number advantage of 0.45 for B.1.617.2 relative to B.1.1.7, with the Rt value peaking at 1.85 for B.1.617.2.

CONCLUSIONS

Our results illustrate a clear transmission advantage for the B.1.617.2 variant and the growth in hospitalisations illustrates that this variant is able to maintain exponential growth within age groups that are largely doubly vaccinated. There are concerning signs of intermittent growth in the B.1.351 variant, reaching a 28-day doubling time peak in March 2021, although this variant is presently not showing any evidence of a transmission advantage over B.1.617.2. Step 1b of the UK national lockdown easing was sufficient to precipitate exponential growth in B.1.617.2 cases for most regions and younger adult age groups. The final stages of NPI easing appeared to have a negligible impact on the growth of B.1.617.2 with every region experiencing sustained exponential growth from step 2. Nonetheless, early targeted local NPIs appeared to markedly reduced growth of B.1.617.2. Later localised interventions, at a time of higher prevalence and greater geographic dispersion of this variant, appeared to have a negligible impact on growth.

SARS-CoV-2 becoming more infectious as revealed by algebraic topology and deep learning

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) caused by coronavirus disease 2019 (COVID-19) has led to a tremendous human fatality and economic loss. SARS-CoV-2 infectivity is a key reason for the widespread viral transmission, but its rigorous experimental measurement is essentially impossible due to the ongoing genome evolution around the world. We show that artificial intelligence (AI) and algebraic topology (AT) offer an accurate and efficient alternative to the experimental determination of viral infectivity. AI and AT analysis indicates that the on-going mutations make SARS-CoV-2 more infectious.

Cytidine deamination-induced perpetual immunity to SAR-CoV-2 infection is a potential new therapeutic target

As the world is racing to develop perpetual immunity to the SARS-CoV-2 virus. The emergence of new viral strains, together with vaccination and reinfections, are all contributing to a long-term immunity against the deadly virus that has taken over the world since its introduction to humans in late December 2019. The discovery that more than 95 percent of people who recovered from COVID-19 had long-lasting immunity and that asymptomatic people have a different immune response to SARS-CoV-2 than symptomatic people has shifted attention to how our immune system initiates such diverse responses. These findings have provided reason to believe that SARS-CoV-2 days are numbered. Hundreds of research papers have been published on the causes of long-lasting immune responses and variations in the numbers of different immune cell types in COVID 19 survivors, but the main reason of these differences has still not been adequately identified. In this article, we focus on the activation-induced cytidine deaminase (AID), which initiates molecular processes that allow our immune system to generate antibodies against SARS-CoV-2. To establish lasting immunity to SARS-CoV-2, we suggest that AID could be the key to unlocking it.

The evolution of the mechanisms of SARS-CoV-2 evolution revealing vaccine-resistant mutations in Europe and America

The importance of understanding SARS-CoV-2 evolution cannot be overemphasized. Recent studies confirm that natural selection is the dominating mechanism of SARS-CoV-2 evolution, which favors mutations that strengthen viral infectivity. We demonstrate that vaccine-breakthrough or antibody-resistant mutations provide a new mechanism of viral evolution. Specifically, vaccine-resistant mutation Y449S in the spike (S) protein receptor-bonding domain (RBD), which occurred in co-mutation [Y449S, N501Y], has reduced infectivity compared to the original SARS-CoV-2 but can disrupt existing antibodies that neutralize the virus. By tracing the evolutionary trajectories of vaccine-resistant mutations in over 1.9 million SARS-CoV-2 genomes, we reveal that the occurrence and frequency of vaccine-resistant mutations correlate strongly with the vaccination rates in Europe and America. We anticipate that as a complementary transmission pathway, vaccine-resistant mutations will become a dominating mechanism of SARS-CoV-2 evolution when most of the world's population is vaccinated. Our study sheds light on SARS-CoV-2 evolution and transmission and enables the design of the next-generation mutation-proof vaccines and antibody drugs.

The mutational landscape of SARS-CoV-2 variants diversifies T cell targets in an HLA-supertype-dependent manner

The rapid, global dispersion of SARS-CoV-2 has led to the emergence of a diverse range of variants. Here, we describe how the mutational landscape of SARS-CoV-2 has shaped HLA-restricted T cell immunity at the population level during the first year of the pandemic. We analyzed a total of 330,246 high-quality SARS-CoV-2 genome assemblies, sampled across 143 countries and all major continents from December 2019 to December 2020 before mass vaccination or the rise of the Delta variant. We observed that proline residues are preferentially removed from the proteome of prevalent mutants, leading to a predicted global loss of SARS-CoV-2 T cell epitopes in individuals expressing HLA-B alleles of the B7 supertype family; this is largely driven by a dominant C-to-U mutation type at the RNA level. These results indicate that B7-supertype-associated epitopes, including the most immunodominant ones, were more likely to escape CD8⁺ T cell immunosurveillance during the first year of the pandemic.

Review of the mechanisms of SARS-CoV-2 evolution and transmission

The mechanism of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) evolution and transmission is elusive and its understanding, a prerequisite to forecast emerging variants, is of paramount importance. SARS-CoV-2 evolution is driven by the mechanisms at molecular and organism scales and regulated by the transmission pathways at the population scale. In this review, we show that infectivity-based natural selection was discovered as the mechanism for SARS-CoV-2 evolution and transmission in July 2020. In April 2021, we proved beyond all doubt that such a natural selection via infectivity-based transmission pathway remained the sole mechanism for SARS-CoV-2 evolution. However, we reveal that antibody-disruptive co-mutations [Y449S, N501Y] on the spike protein receptor-binding domain (RBD) debuted as a new vaccine-resistant transmission pathway of viral evolution in highly vaccinated populations a few months ago. Over one year ago, we foresaw that mutations on RBD residues, 452 and 501, would "both have high chances to mutate into significantly more infectious COVID-19 strains". Mutations on these residues underpin prevailing SARS-CoV-2 variants Alpha, Beta, Gamma, Delta, Epsilon, Theta, Kappa, Lambda, and Mu at present and are expected to be vital to emerging variants in the future. We anticipate that viral evolution will combine RBD co-mutations at these two sites, creating future variants that are about ten times more infectious than the original SARS-CoV-2. Additionally, two complementary transmission pathways of viral evolution, i.e., infectivity and vaccine resistance will prolong our battle with COVID-19 for years. We predict that RBD co-mutation sets [A411S, L452R, T478K], [L452R, T478K, N501Y], [L452R, T478K, E484K, N501Y], [K417N, L452R, T478K], and [P384L, K417N, E484K, N501Y] will have a high chance to grow into dominating variants due to their high infectivity and/or strong ability to break through current vaccines, calling for the development of new vaccines and antibody therapies.

Emerging vaccine-breakthrough SARS-CoV-2 variants

The recent global surge in COVID-19 infections has been fueled by new SARS-CoV-2 variants, namely Alpha, Beta, Gamma, Delta, etc. The molecular mechanism underlying such surge is elusive due to 4,653 non-degenerate mutations on the spike protein, which is the target of most COVID-19 vaccines. The understanding of the molecular mechanism of transmission and evolution is a prerequisite to foresee the trend of emerging vaccine-breakthrough variants and the design of mutation-proof vaccines and monoclonal antibodies. We integrate the genotyping of 1,489,884 SARS-CoV-2 genomes isolates, 130 human antibodies, tens of thousands of mutational data points, topological data analysis, and deep learning to reveal SARS-CoV-2 evolution mechanism and forecast emerging vaccine-escape variants. We show that infectivity-strengthening and antibody-disruptive co-mutations on the S protein RBD can quantitatively explain the infectivity and virulence of all prevailing variants. We demonstrate that Lambda is as infectious as Delta but is more vaccine-resistant. We analyze emerging vaccine-breakthrough co-mutations in 20 countries, including the United Kingdom, the United States, Denmark, Brazil, and Germany, etc. We envision that natural selection through infectivity will continue to be the main mechanism for viral evolution among unvaccinated populations, while antibody disruptive co-mutations will fuel the future growth of vaccine-breakthrough variants among fully vaccinated populations. Finally, we have identified the co-mutations that have the great likelihood of becoming dominant: [A411S, L452R, T478K], [L452R, T478K, N501Y], [V401L, L452R, T478K], [K417N, L452R, T478K], [L452R, T478K, E484K, N501Y], and [P384L, K417N, E484K, N501Y]. We predict they, particularly the last four, will break through existing vaccines. We foresee an urgent need to develop new vaccines that target these co-mutations.

Epigenetic Mechanisms Underlying COVID-19 Pathogenesis

In 2019, a novel severe acute respiratory syndrome called coronavirus disease 2019 (COVID-19), caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), was reported and was declared a pandemic by the World Health Organization (WHO) in March 2020. With the advancing development of COVID-19 vaccines and their administration globally, it is expected that COVID-19 will converge in the future; however, the situation remains unpredictable because of a series of reports regarding SARS-CoV-2 variants. Currently, there are still few specific effective treatments for COVID-19, as many unanswered questions remain regarding the pathogenic mechanism of COVID-19. Continued elucidation of COVID-19 pathogenic mechanisms is a matter of global importance. In this regard, recent reports have suggested that epigenetics plays an important role; for instance, the expression of angiotensin I converting enzyme 2 (ACE2) receptor, an important factor in human infection with SARS-CoV-2, is epigenetically regulated; further, DNA methylation status is reported to be unique to patients with COVID-19. In this review, we focus on epigenetic mechanisms to provide a new molecular framework for elucidating the pathogenesis of SARS-CoV-2 infection in humans and of COVID-19, along with the possibility of new diagnostic and therapeutic strategies.

Recent progress on the mutations of SARS-CoV-2 spike protein and suggestions for prevention and controlling of the pandemic

Severe acute respiratory syndrome coronavirus type 2 (SARS-CoV-2) infection has caused a global pandemic in the past year, which poses continuing threat to human beings. To date, more than 3561 mutations in the viral spike protein were identified, including 2434 mutations that cause amino acid changes with 343 amino acids located in the viral receptor-binding domain (RBD). Among these mutations, the most representative ones are substitution mutations such as D614G, N501Y, Y453F, N439K/R, P681H, K417N/T, and E484K, and deletion mutations of ΔH69/V70 and Δ242-244, which confer the virus with enhanced infectivity, transmissibility, and resistance to neutralization. In this review, we discussed the recent findings of SARS-CoV-2 for highlighting mutations and variants on virus transmissibility and pathogenicity. Moreover, several suggestions for prevention and controlling the pandemic are also proposed.

Mutational Landscape and Interaction of SARS-CoV-2 with Host Cellular Components

The emergence of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) and its rapid evolution has led to a global health crisis. Increasing mutations across the SARS-CoV-2 genome have severely impacted the development of effective therapeutics and vaccines to combat the virus. However, the new SARS-CoV-2 variants and their evolutionary characteristics are not fully understood. Host cellular components such as the ACE2 receptor, RNA-binding proteins (RBPs), microRNAs, small nuclear RNA (snRNA), 18s rRNA, and the 7SL RNA component of the signal recognition particle (SRP) interact with various structural and non-structural proteins of the SARS-CoV-2. Several of these viral proteins are currently being examined for designing antiviral therapeutics. In this review, we discuss current advances in our understanding of various host cellular components targeted by the virus during SARS-CoV-2 infection. We also summarize the mutations across the SARS-CoV-2 genome that directs the evolution of new viral strains. Considering coronaviruses are rapidly evolving in humans, this enables them to escape therapeutic therapies and vaccine-induced immunity. In order to understand the virus's evolution, it is essential to study its mutational patterns and their impact on host cellular machinery. Finally, we present a comprehensive survey of currently available databases and tools to study viral-host interactions that stand as crucial resources for developing novel therapeutic strategies for combating SARS-CoV-2 infection.

Effect of Train-Induced Wind on the Transmission of COVID-19: A New Insight into Potential Infectious Risks

The unprecedented COVID-19 pandemic has caused a traffic tie-up across the world. In addition to home quarantine orders and travel bans, the social distance guideline of about six feet was enacted to reduce the risk of contagion. However, with recent life gradually returning to normal, the crisis is not over. In this research, a moving train test and a Gaussian puff model were employed to investigate the impact of wind raised by a train running on the transmission and dispersion of SARS-CoV-2 from infected individuals. Our findings suggest that the 2 m social distance guideline may not be enough; under train-induced wind action, human respiratory disease-carrier droplets may travel to unexpected places. However, there are deficiencies in passenger safety guidelines and it is necessary to improve the quantitative research in the relationship between train-induced wind and virus transmission. All these findings could provide a fresh insight to contain the spread of COVID-19 and provide a basis for preventing and controlling the pandemic virus, and probe into strategies for control of the disease in the future.

SARS-CoV-2 variant evolution in the United States: High accumulation of viral mutations over time likely through serial Founder Events and mutational bursts

Since the first case of COVID-19 in December 2019 in Wuhan, China, SARS-CoV-2 has spread worldwide and within a year and a half has caused 3.56 million deaths globally. With dramatically increasing infection numbers, and the arrival of new variants with increased infectivity, tracking the evolution of its genome is crucial for effectively controlling the pandemic and informing vaccine platform development. Our study explores evolution of SARS-CoV-2 in a representative cohort of sequences covering the entire genome in the United States, through all of 2020 and early 2021. Strikingly, we detected many accumulating Single Nucleotide Variations (SNVs) encoding amino acid changes in the SARS-CoV-2 genome, with a pattern indicative of RNA editing enzymes as major mutators of SARS-CoV-2 genomes. We report three major variants through October of 2020. These revealed 14 key mutations that were found in various combinations among 14 distinct predominant signatures. These signatures likely represent evolutionary lineages of SARS-CoV-2 in the U.S. and reveal clues to its evolution such as a mutational burst in the summer of 2020 likely leading to a homegrown new variant, and a trend towards higher mutational load among viral isolates, but with occasional mutation loss. The last quartile of 2020 revealed a concerning accumulation of mostly novel low frequency replacement mutations in the Spike protein, and a hypermutable glutamine residue near the putative furin cleavage site. Finally, end of the year data and 2021 revealed the gradual increase to prevalence of known variants of concern, particularly B.1.1.7, that have acquired additional Spike mutations. Overall, our results suggest that predominant viral genomes are dynamically evolving over time, with periods of mutational bursts and unabated mutation accumulation. This high level of existing variation, even at low frequencies and especially in the Spike-encoding region may become problematic when super-spreader events, akin to serial Founder Events in evolution, drive these rare mutations to prominence.

Possible Targets of Pan-Coronavirus Antiviral Strategies for Emerging or Re-Emerging Coronaviruses

Severe Acute Respiratory Syndrome Coronavirus-2 (SARS-CoV-2), which caused Coronaviruses Disease 2019 (COVID-19) and a worldwide pandemic, is the seventh human coronavirus that has been cross-transmitted from animals to humans. It can be predicted that with continuous contact between humans and animals, more viruses will spread from animals to humans. Therefore, it is imperative to develop universal coronavirus or pan-coronavirus vaccines or drugs against the next coronavirus pandemic. However, a suitable target is critical for developing pan-coronavirus antivirals against emerging or re-emerging coronaviruses. In this review, we discuss the latest progress of possible targets of pan-coronavirus antiviral strategies for emerging or re-emerging coronaviruses, including targets for pan-coronavirus inhibitors and vaccines, which will provide prospects for the current and future research and treatment of the disease.

Vaccine-escape and fast-growing mutations in the United Kingdom, the United States, Singapore, Spain, India, and other COVID-19-devastated countries

Recently, the SARS-CoV-2 variants from the United Kingdom (UK), South Africa, and Brazil have received much attention for their increased infectivity, potentially high virulence, and possible threats to existing vaccines and antibody therapies. The question remains if there are other more infectious variants transmitted around the world. We carry out a large-scale study of 506,768 SARS-CoV-2 genome isolates from patients to identify many other rapidly growing mutations on the spike (S) protein receptor-binding domain (RBD). We reveal that essentially all 100 most observed mutations strengthen the binding between the RBD and the host angiotensin-converting enzyme 2 (ACE2), indicating the virus evolves toward more infectious variants. In particular, we discover new fast-growing RBD mutations N439K, S477N, S477R, and N501T that also enhance the RBD and ACE2 binding. We further unveil that mutation N501Y involved in United Kingdom (UK), South Africa, and Brazil variants may moderately weaken the binding between the RBD and many known antibodies, while mutations E484K and K417N found in South Africa and Brazilian variants, L452R and E484Q found in India variants, can potentially disrupt the binding between the RBD and many known antibodies. Among these RBD mutations, L452R is also now known as part of the California variant B.1.427. Finally, we hypothesize that RBD mutations that can simultaneously make SARS-CoV-2 more infectious and disrupt the existing antibodies, called vaccine escape mutations, will pose an imminent threat to the current crop of vaccines. A list of most likely vaccine escape mutations is given, including S494P, Q493L, K417N, F490S, F486L, R403K, E484K, L452R, K417T, F490L, E484Q, and A475S. Mutation T478K appears to make the Mexico variant B.1.1.222 the most infectious one. Our comprehensive genetic analysis and protein-protein binding study show that the genetic evolution of SARS-CoV-2 on the RBD, which may be regulated by host gene editing, viral proofreading, random genetic drift, and natural selection, gives rise to more infectious variants that will potentially compromise existing vaccines and antibody therapies.

Positive selection as a key player for SARS-CoV-2 pathogenicity: Insights into ORF1ab, S and E genes

The human β-coronavirus SARS-CoV-2 epidemic started in late December 2019 in Wuhan, China. It causes Covid-19 disease which has become pandemic. Each of the five-known human β-coronaviruses has four major structural proteins (E, M, N and S) and 16 non-structural proteins encoded by ORF1a and ORF1b together (ORF1ab) that are involved in virus pathogenicity and infectivity. Here, we performed detailed positive selection analyses for those six genes among the four previously known human β-coronaviruses and within 38 SARS-CoV-2 genomes to assess signatures of adaptive evolution using maximum likelihood approaches. Our results suggest that three genes (E, S and ORF1ab genes) are under strong signatures of positive selection among human β-coronavirus, influencing codons that are located in functional important protein domains. The E protein-coding gene showed signatures of positive selection in two sites, Asp 66 and Ser 68, located inside a putative transmembrane α-helical domain C-terminal part, which is preferentially composed by hydrophilic residues. Such Asp and Ser sites substitutions (hydrophilic residues) increase the stability of the transmembrane domain in SARS-CoV-2. Moreover, substitutions in the spike (S) protein S1 N-terminal domain have been found, all of them were located on the S protein surface, suggesting their importance in viral transmissibility and survival. Furthermore, evidence of strong positive selection was detected in three of the SARS-CoV-2 nonstructural proteins (NSP1, NSP3, NSP16), which are encoded by ORF1ab and play vital roles in suppressing host translation machinery, viral replication and transcription and inhibiting the host immune response. These results are insightful to assess the role of positive selection in the SARS-CoV-2 encoded proteins, which will allow to better understand the virulent pathogenicity of the virus and potentially identifying targets for drug or vaccine strategy design

Mutations of SARS-CoV-2 spike protein: Implications on immune evasion and vaccine-induced immunity

Responsible for more than 4.9 million deaths so far, COVID-19, caused by SARS-CoV-2, is instigating devastating effects on the global health care system whose impacts could be longer for the years to come. Acquiring a comprehensive knowledge of host-virus interaction is critical for designing effective vaccines and/or drugs. Understanding the evolution of the virus and the impact of genetic variability on host immune evasion and vaccine efficacy is helpful to design novel strategies to minimize the effects of the emerging variants of concern (VOC). Most vaccines under development and/or in current use target the spike protein owning to its unique function of host receptor binding, relatively conserved nature, potent immunogenicity in inducing neutralizing antibodies, and being a good target of T cell responses. However, emerging SARS-CoV-2 strains are exhibiting variability on the spike protein which could affect the efficacy of vaccines and antibody-based therapies in addition to enhancing viral immune evasion mechanisms. Currently, the degree to which mutations on the spike protein affect immunity and vaccination, and the ability of the current vaccines to confer protection against the emerging variants attracts much attention. This review discusses the implications of SARS-CoV-2 spike protein mutations on immune evasion and vaccine-induced immunity and forward directions which could contribute to future studies focusing on designing effective vaccines and/or immunotherapies to consider viral evolution. Combining vaccines derived from different regions of the spike protein that boost both the humoral and cellular wings of adaptive immunity could be the best options to cope with the emerging VOC.

Genomic variation, origin tracing, and vaccine development of SARS-CoV-2: A systematic review

COVID-19 has spread globally to over 200 countries with more than 40 million confirmed cases and one million deaths as of November 1, 2020. The SARS-CoV-2 virus, leading to COVID-19, shows extremely high rates of infectivity and replication, and can result in pneumonia, acute respiratory distress, or even mortality. SARS-CoV-2 has been found to continue to rapidly evolve, with several genomic variants emerging in different regions throughout the world. In addition, despite intensive study of the spike protein, its origin, and molecular mechanisms in mediating host invasion are still only partially resolved. Finally, the repertoire of drugs for COVID-19 treatment is still limited, with several candidates still under clinical trial and no effective therapeutic yet reported. Although vaccines based on either DNA/mRNA or protein have been deployed, their efficacy against emerging variants requires ongoing study, with multivalent vaccines supplanting the first-generation vaccines due to their low efficacy against new strains. Here, we provide a systematic review of studies on the epidemiology, immunological pathogenesis, molecular mechanisms, and structural biology, as well as approaches for drug or vaccine development for SARS-CoV-2.

Mutational Asymmetries in the SARS-CoV-2 Genome May Lead to Increased Hydrophobicity of Virus Proteins

The genomic diversity of SARS-CoV-2 has been a focus during the ongoing COVID-19 pandemic. Here, we analyzed the distribution and character of emerging mutations in a data set comprising more than 95,000 virus genomes covering eight major SARS-CoV-2 lineages in the GISAID database, including genotypes arising during COVID-19 therapy. Globally, the C>U transitions and G>U transversions were the most represented mutations, accounting for the majority of single-nucleotide variations. Mutational spectra were not influenced by the time the virus had been circulating in its host or medical treatment. At the amino acid level, we observed about a 2-fold excess of substitutions in favor of hydrophobic amino acids over the reverse. However, most mutations constituting variants of interests of the S-protein (spike) lead to hydrophilic amino acids, counteracting the global trend. The C>U and G>U substitutions altered codons towards increased amino acid hydrophobicity values in more than 80% of cases. The bias is explained by the existing differences in the codon composition for amino acids bearing contrasting biochemical properties. Mutation asymmetries apparently influence the biochemical features of SARS CoV-2 proteins, which may impact protein-protein interactions, fusion of viral and cellular membranes, and virion assembly.

Genome-wide bioinformatic analyses predict key host and viral factors in SARS-CoV-2 pathogenesis

The novel betacoronavirus severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) caused a worldwide pandemic (COVID-19) after emerging in Wuhan, China. Here we analyzed public host and viral RNA sequencing data to better understand how SARS-CoV-2 interacts with human respiratory cells. We identified genes, isoforms and transposable element families that are specifically altered in SARS-CoV-2-infected respiratory cells. Well-known immunoregulatory genes including CSF2, IL32, IL-6 and SERPINA3 were differentially expressed, while immunoregulatory transposable element families were upregulated. We predicted conserved interactions between the SARS-CoV-2 genome and human RNA-binding proteins such as the heterogeneous nuclear ribonucleoprotein A1 (hnRNPA1) and eukaryotic initiation factor 4 (eIF4b). We also identified a viral sequence variant with a statistically significant skew associated with age of infection, that may contribute to intracellular host-pathogen interactions. These findings can help identify host mechanisms that can be targeted by prophylactics and/or therapeutics to reduce the severity of COVID-19.

Mutation Rates and Selection on Synonymous Mutations in SARS-CoV-2

The COVID-19 pandemic has seen an unprecedented response from the sequencing community. Leveraging the sequence data from more than 140,000 SARS-CoV-2 genomes, we study mutation rates and selective pressures affecting the virus. Understanding the processes and effects of mutation and selection has profound implications for the study of viral evolution, for vaccine design, and for the tracking of viral spread. We highlight and address some common genome sequence analysis pitfalls that can lead to inaccurate inference of mutation rates and selection, such as ignoring skews in the genetic code, not accounting for recurrent mutations, and assuming evolutionary equilibrium. We find that two particular mutation rates, G →U and C →U, are similarly elevated and considerably higher than all other mutation rates, causing the majority of mutations in the SARS-CoV-2 genome, and are possibly the result of APOBEC and ROS activity. These mutations also tend to occur many times at the same genome positions along the global SARS-CoV-2 phylogeny (i.e., they are very homoplasic). We observe an effect of genomic context on mutation rates, but the effect of the context is overall limited. Although previous studies have suggested selection acting to decrease U content at synonymous sites, we bring forward evidence suggesting the opposite.