Short sequence motif dynamics in the SARS-CoV-2 genome suggest a role for cytosine deamination in CpG reduction

Differential RNA editing landscapes in host cell versus the SARS-CoV-2 genome

The SARS-CoV-2 pandemic was defined by the emergence of new variants formed through virus mutation originating from random errors not corrected by viral proofreading and/or the host antiviral response introducing mutations into the viral genome. While sequencing information hints at cellular RNA editing pathways playing a role in viral evolution, here, we use an in vitro human cell infection model to assess RNA mutation types in two SARS-CoV-2 strains representing the original and the alpha variants. The variants showed both different cellular responses and mutation patterns with alpha showing higher mutation frequency with most substitutions observed being C-U, indicating an important role for apolipoprotein B mRNA editing catalytic polypeptide-like editing. Knockdown of select APOBEC3s through RNAi increased virus production in the original virus, but not in alpha. Overall, these data suggest a deaminase-independent anti-viral function of APOBECs in SARS-CoV-2 while the C-U editing itself might function to enhance genetic diversity enabling evolutionary adaptation.

Fitness effects of mutations to SARS-CoV-2 proteins

Knowledge of the fitness effects of mutations to SARS-CoV-2 can inform assessment of new variants, design of therapeutics resistant to escape, and understanding of the functions of viral proteins. However, experimentally measuring effects of mutations is challenging: we lack tractable lab assays for many SARS-CoV-2 proteins, and comprehensive deep mutational scanning has been applied to only two SARS-CoV-2 proteins. Here, we develop an approach that leverages millions of publicly available SARS-CoV-2 sequences to estimate effects of mutations. We first calculate how many independent occurrences of each mutation are expected to be observed along the SARS-CoV-2 phylogeny in the absence of selection. We then compare these expected observations to the actual observations to estimate the effect of each mutation. These estimates correlate well with deep mutational scanning measurements. For most genes, synonymous mutations are nearly neutral, stop-codon mutations are deleterious, and amino acid mutations have a range of effects. However, some viral accessory proteins are under little to no selection. We provide interactive visualizations of effects of mutations to all SARS-CoV-2 proteins (https://jbloomlab.github.io/SARS2-mut-fitness/). The framework we describe is applicable to any virus for which the number of available sequences is sufficiently large that many independent occurrences of each neutral mutation are observed.

Fitness effects of mutations to SARS-CoV-2 proteins

Knowledge of the fitness effects of mutations to SARS-CoV-2 can inform assessment of new variants, design of therapeutics resistant to escape, and understanding of the functions of viral proteins. However, experimentally measuring effects of mutations is challenging: we lack tractable lab assays for many SARS-CoV-2 proteins, and comprehensive deep mutational scanning has been applied to only two SARS-CoV-2 proteins. Here we develop an approach that leverages millions of publicly available SARS-CoV-2 sequences to estimate effects of mutations. We first calculate how many independent occurrences of each mutation are expected to be observed along the SARS-CoV-2 phylogeny in the absence of selection. We then compare these expected observations to the actual observations to estimate the effect of each mutation. These estimates correlate well with deep mutational scanning measurements. For most genes, synonymous mutations are nearly neutral, stop-codon mutations are deleterious, and amino-acid mutations have a range of effects. However, some viral accessory proteins are under little to no selection. We provide interactive visualizations of effects of mutations to all SARS-CoV-2 proteins (https://jbloomlab.github.io/SARS2-mut-fitness/). The framework we describe is applicable to any virus for which the number of available sequences is sufficiently large that many independent occurrences of each neutral mutation are observed.

A Synthetic SARS-CoV-2-Derived T-Cell and B-Cell Peptide Cocktail Elicits Full Protection against Lethal Omicron BA.1 Infection in H11-K18-hACE2 Mice

Emerging variants of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) have been developing the capacity for immune evasion and resistance to existing vaccines and drugs. To address this, development of vaccines against coronavirus disease 2019 (COVID-19) has focused on universality, strong T cell immunity, and rapid production. Synthetic peptide vaccines, which are inexpensive and quick to produce, show low toxicity, and can be selected from the conserved SARS-CoV-2 proteome, are promising candidates. In this study, we evaluated the effectiveness of a synthetic peptide cocktail containing three murine CD4+ T-cell epitopes from the SARS-CoV-2 nonspike proteome and one B-cell epitope from the Omicron BA.1 receptor-binding domain (RBD), along with aluminum phosphate (Al) adjuvant and 5' cytosine-phosphate-guanine 3' oligodeoxynucleotide (CpG-ODN) adjuvant in mice. The peptide cocktail induced good Th1-biased T-cell responses and effective neutralizing-antibody titers against the Omicron BA.1 variant. Additionally, H11-K18-hACE2 transgenic mice were fully protected against lethal challenge with the BA.1 strain, with a 100% survival rate and reduced pulmonary viral load and pathological lesions. Subcutaneous administration was found to be the superior route for synthetic peptide vaccine delivery. Our findings demonstrate the effectiveness of the peptide cocktail in mice, suggesting the feasibility of synthetic peptide vaccines for humans. IMPORTANCE Current vaccines based on production of neutralizing antibodies fail to prevent the infection and transmission of SARS-CoV-2 Omicron and its subvariants. Understanding the critical factors and avoiding the disadvantages of vaccine strategies are essential for developing a safe and effective COVID-19 vaccine, which would include a more effective and durable cellular response, minimal effects of viral mutations, rapid production against emerging variants, and good safety. Peptide-based vaccines are an excellent alternative because they are inexpensive, quick to produce, and very safe. In addition, human leukocyte antigen T-cell epitopes could be targeted at robust T-cell immunity and selected in the conserved region of the SARS-CoV-2 variants. Our study showed that a synthetic SARS-CoV-2-derived peptide cocktail induced full protection against lethal infection with Omicron BA.1 in H11-K18-hACE2 mice for the first time. This could have implications for the development of effective COVID-19 peptide vaccines for humans.

Control of Human Anelloviruses by Cytosine to Uracil Genome Editing

Anelloviruses are the most common viruses infecting humans. Every human carries a nonpathogenic personal anellovirus virome (anellome), yet it is unknown which mechanisms contribute to its stability. Here, we assessed the dynamics and impact of a host antiviral defense mechanism-cytidine deaminase activity leading to C to U editing in anelloviruses-on the stability of the anellome. We investigated anellome sequence data obtained from serum samples collected every 6 months from two healthy subjects followed for more than 30 years. The subjects were infected by a total of 64 anellovirus lineages. Minus-stranded C to U editing was observed in lineages belonging to the Alpha-, Beta-, and Gammatorquevirus genera. The edited genomes were present within virus particles, therefore editing must have occurred at the late stages of the virus life cycle. Editing was favored by 5'-TC contexts in the virus genome, indicating that apolipoprotein B mRNA-editing enzyme, catalytic polypeptide-like, catalytic subunit 3 or A3 (APOBEC3) proteins are involved. Within a lineage, mutational dynamics varied over time and few fixations of mutations were detected, indicating that C to U editing is a dead end for a virus genome. We detected an editing coldspot in the GC-rich regions, suggesting that the GC-rich region is crucial for genome packaging, since only packaged virus particles were included in the analysis. Finally, we noticed a lineage-specific reduced concentration after an editing event, yet no clearance. In conclusion, cytidine deaminase activity does not clear anelloviruses, nor does it play a major role in virus evolution, but it does contribute to the stability of the anellome. IMPORTANCE Despite significant attention on anellovirus research, the interaction between the anellovirus virome and the human host remains unknown. We show the dynamics of APOBEC3-mediated cytidine deaminase activity on anelloviruses during a 30-year period of chronic infection and postulate that this antiviral mechanism controls anelloviruses. These results expand our knowledge of anellovirus-host interactions, which may be important for the design of gene therapies.

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.

Commentary on "Poor evidence for host-dependent regular RNA editing in the transcriptome of SARS-CoV-2"

Analysis of the SARS-CoV-2 transcriptome has revealed a background of low-frequency intra-host genetic changes with a strong bias towards transitions. A similar pattern is also observed when inter-host variability is considered. We and others have shown that the cellular RNA editing machinery based on ADAR and APOBEC host-deaminases could be involved in the onset of SARS-CoV-2 genetic variability. Our hypothesis is based both on similarities with other known forms of viral genome editing and on the excess of transition changes, which is difficult to explain with errors during viral replication. Zong et al. criticize our analysis on both conceptual and technical grounds. While ultimate proof of an involvement of host deaminases in viral RNA editing will depend on experimental validation, here, we address the criticism to suggest that viral RNA editing is the most reasonable explanation for the observed intra- and inter-host variability.

Statistical modeling of SARS-CoV-2 substitution processes: predicting the next variant

We build statistical models to describe the substitution process in the SARS-CoV-2 as a function of explanatory factors describing the sequence, its function, and more. These models serve two different purposes: first, to gain knowledge about the evolutionary biology of the virus; and second, to predict future mutations in the virus, in particular, non-synonymous amino acid substitutions creating new variants. We use tens of thousands of publicly available SARS-CoV-2 sequences and consider tens of thousands of candidate models. Through a careful validation process, we confirm that our chosen models are indeed able to predict new amino acid substitutions: candidates ranked high by our model are eight times more likely to occur than random amino acid changes. We also show that named variants were highly ranked by our models before their appearance, emphasizing the value of our models for identifying likely variants and potentially utilizing this knowledge in vaccine design and other aspects of the ongoing battle against COVID-19.

As the virus that causes COVID-19 continues to mutate and spread, new methods are needed to predict new potential variants. Here, the authors identify the best regression models for predicting likely mutation sites in the SARS-CoV-2 genome using a candidate set that considers sequence, gene location, and biological function.

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

The Heterogeneous Landscape and Early Evolution of Pathogen-Associated CpG Dinucleotides in SARS-CoV-2

COVID-19 can lead to acute respiratory syndrome, which can be due to dysregulated immune signaling. We analyze the distribution of CpG dinucleotides, a pathogen-associated molecular pattern, in the SARS-CoV-2 genome. We characterize CpG content by a CpG force that accounts for statistical constraints acting on the genome at the nucleotidic and amino acid levels. The CpG force, as the CpG content, is overall low compared with other pathogenic betacoronaviruses; however, it widely fluctuates along the genome, with a particularly low value, comparable with the circulating seasonal HKU1, in the spike coding region and a greater value, comparable with SARS and MERS, in the highly expressed nucleocapside coding region (N ORF), whose transcripts are relatively abundant in the cytoplasm of infected cells and present in the 3'UTRs of all subgenomic RNA. This dual nature of CpG content could confer to SARS-CoV-2 the ability to avoid triggering pattern recognition receptors upon entry, while eliciting a stronger response during replication. We then investigate the evolution of synonymous mutations since the outbreak of the COVID-19 pandemic, finding a signature of CpG loss in regions with a greater CpG force. Sequence motifs preceding the CpG-loss-associated loci in the N ORF match recently identified binding patterns of the zinc finger antiviral protein. Using a model of the viral gene evolution under human host pressure, we find that synonymous mutations seem driven in the SARS-CoV-2 genome, and particularly in the N ORF, by the viral codon bias, the transition-transversion bias, and the pressure to lower CpG content.

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.

Glimpses into evolutionary trajectories of SARS-CoV-2: emerging variants and potential immune evasion routes

Tweetable abstract An opinion on the coronaviruses' evolution paradoxes, the continuing adaptation of the SARS-CoV-2 in humans following the zoonotic transmission, and clues into escape routes from host immune responses.

Host-directed editing of the SARS-CoV-2 genome

The extensive sequence data generated from SARS-CoV-2 during the 2020 pandemic has facilitated the study of viral genome evolution over a brief period of time. This has highlighted instances of directional mutation pressures exerted on the SARS-CoV-2 genome from host antiviral defense systems. In this brief review we describe three such human defense mechanisms, the apolipoprotein B mRNA editing catalytic polypeptide-like proteins (APOBEC), adenosine deaminase acting on RNA proteins (ADAR), and reactive oxygen species (ROS), and discuss their potential implications on SARS-CoV-2 evolution.