The slowing rate of CpG depletion in SARS-CoV-2 genomes is consistent with adaptations to the human host

Selection on synonymous sites: the unwanted transcript hypothesis

Although translational selection to favour codons that match the most abundant tRNAs is not readily observed in humans, there is nonetheless selection in humans on synonymous mutations. We hypothesize that much of this synonymous site selection can be explained in terms of protection against unwanted RNAs - spurious transcripts, mis-spliced forms or RNAs derived from transposable elements or viruses. We propose not only that selection on synonymous sites functions to reduce the rate of creation of unwanted transcripts (for example, through selection on exonic splice enhancers and cryptic splice sites) but also that high-GC content (but low-CpG content), together with intron presence and position, is both particular to functional native mRNAs and used to recognize transcripts as native. In support of this hypothesis, transcription, nuclear export, liquid phase condensation and RNA degradation have all recently been shown to promote GC-rich transcripts and suppress AU/CpG-rich ones. With such 'traps' being set against AU/CpG-rich transcripts, the codon usage of native genes has, in turn, evolved to avoid such suppression. That parallel filters against AU/CpG-rich transcripts also affect the endosomal import of RNAs further supports the unwanted transcript hypothesis of synonymous site selection and explains the similar design rules that have enabled the successful use of transgenes and RNA vaccines.

Mapping of long stretches of highly conserved sequences in over 6 million SARS-CoV-2 genomes

We identified 11 conserved stretches in over 6.3 million SARS-CoV-2 genomes including all the major variants of concerns. Each conserved stretch is ≥100 nucleotides in length with ≥99.9% conservation at each nucleotide position. Interestingly, six of the eight conserved stretches in ORF1ab overlapped significantly with well-folded experimentally verified RNA secondary structures. Furthermore, two of the conserved stretches were mapped to regions within the S2-subunit that undergo dynamic structural rearrangements during viral fusion. In addition, the conserved stretches were significantly depleted for zinc-finger antiviral protein (ZAP) binding sites, which facilitated the recognition and degradation of viral RNA. These highly conserved stretches in the SARS-CoV-2 genome were poorly conserved at the nucleotide level among closely related β-coronaviruses, thus representing ideal targets for highly specific and discriminatory diagnostic assays. Our findings highlight the role of structural constraints at both RNA and protein levels that contribute to the sequence conservation of specific genomic regions in SARS-CoV-2.

Analysis of SARS-CoV-2 genome evolutionary patterns

The spread of SARS-CoV-2 virus accompanied by public availability of abundant sequence data provides a window for the determination of viral evolutionary patterns. In this study, SARS-CoV-2 genome sequences were collected from seven countries in the period January 2020-December 2022. The sequences were classified into three phases, namely, pre-vaccination, post-vaccination, and recent period. Comparison was performed between these phases based on parameters like mutation rates, selection pressure (dN/dS ratio), and transition to transversion ratios (Ti/Tv). Similar comparisons were performed among SARS-CoV-2 variants. Statistical significance was tested using Graphpad unpaired t-test. The analysis showed an increase in the percent genomic mutation rates post-vaccination and in recent periods across all countries from the pre-vaccination sequences. Mutation rates were highest in NSP3, S, N, and NSP12b before and increased further after vaccination. NSP4 showed the largest change in mutation rates after vaccination. The dN/dS ratios showed purifying selection that shifted toward neutral selection after vaccination. N, ORF8, ORF3a, and ORF10 were under highest positive selection before vaccination. Shift toward neutral selection was driven by E, NSP3, and ORF7a in the after vaccination set. In recent sequences, the largest dN/dS change was observed in E, NSP1, and NSP13. The Ti/Tv ratios decreased with time. C→U and G→U were the most frequent transitions and transversions. However, U→G was the most frequent transversion in recent period. The Omicron variant had the highest genomic mutation rates, while Delta showed the highest dN/dS ratio. Protein-wise dN/dS ratio was also seen to vary across the different variants.IMPORTANCETo the best of our knowledge, there exists no other large-scale study of the genomic and protein-wise mutation patterns during the time course of evolution in different countries. Analyzing the SARS-CoV-2 evolutionary patterns in view of the varying spatial, temporal, and biological signals is important for diagnostics, therapeutics, and pharmacovigilance of SARS-CoV-2.

SARS-CoV-2 and innate immunity: the good, the bad, and the "goldilocks"

An ancient conflict between hosts and pathogens has driven the innate and adaptive arms of immunity. Knowledge about this interplay can not only help us identify biological mechanisms but also reveal pathogen vulnerabilities that can be leveraged therapeutically. The humoral response to SARS-CoV-2 infection has been the focus of intense research, and the role of the innate immune system has received significantly less attention. Here, we review current knowledge of the innate immune response to SARS-CoV-2 infection and the various means SARS-CoV-2 employs to evade innate defense systems. We also consider the role of innate immunity in SARS-CoV-2 vaccines and in the phenomenon of long COVID.

Optimal entropic properties of SARS-CoV-2 RNA sequences

The reaction of the scientific community against the COVID-19 pandemic has generated a huge (approx. 106 entries) dataset of genome sequences collected worldwide and spanning a relatively short time window. These unprecedented conditions together with the certain identification of the reference viral genome sequence allow for an original statistical study of mutations in the virus genome. In this paper, we compute the Shannon entropy of every sequence in the dataset as well as the relative entropy and the mutual information between the reference sequence and the mutated ones. These functions, originally developed in information theory, measure the information content of a sequence and allows us to study the random character of mutation mechanism in terms of its entropy and information gain or loss. We show that this approach allows us to set in new format known features of the SARS-CoV-2 mutation mechanism like the CT bias, but also to discover new optimal entropic properties of the mutation process in the sense that the virus mutation mechanism track closely theoretically computable lower bounds for the entropy decrease and the information transfer.

Long-term dynamic shifts in genomic base content and evolutionary trajectories of SARS-CoV-2 variants

The rapid spread and remarkable mutations of SARS-CoV-2 variants, particularly Omicron, necessitate an understanding of their evolutionary characteristics. In this study, we analyzed representative high-quality whole-genome sequences of 2008 SARS-CoV-2 variants to explore long-term dynamic changes in genomic base (especially GC) content and variations during viral evolution. Our results demonstrated a highly negative correlation between GC content and variant emergence time (r = -0.765, p < 2.22e-16). Major gene partitions (S, N, ORF1ab) displayed similar trends. Omicron exhibited a significantly lower GC content than non-Omicron variants (p < 2.22e-16). Notably, we observed a robust negative correlation between C and T content (r = -0.778, p < 2.22e-16) and between G and A content (r = -0.773, p < 2.22e-16). Among all strains, Omicron showed the greatest base variation, with C->T mutations being the most frequent (median [interquartile range [IQR]]: 29 (27, 31), 37.67%), succeeded by G->A mutations (11 (9, 13), 14.63%). Over a 3-year span, an annual decline rate of 0.12% in SARS-CoV-2 GC content was observed and could become more pronounced in future emerging variants. These findings provided insights into the evolutionary trajectory of SARS-CoV-2, underscoring the significance of continuous genomic surveillance for effective prediction of and response to future variants.

SARS-CoV-2 within-host diversity of human hosts and its implications for viral immune evasion

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is continuously evolving, bringing great challenges to the control of the virus. In the present study, we investigated the characteristics of SARS-CoV-2 within-host diversity of human hosts and its implications for immune evasion using about 2,00,000 high-depth next-generation genome sequencing data of SARS-CoV-2. A total of 44% of the samples showed within-host variations (iSNVs), and the average number of iSNVs in the samples with iSNV was 1.90. C-to-U is the dominant substitution pattern for iSNVs. C-to-U/G-to-A and A-to-G/U-to-C preferentially occur in 5'-CG-3' and 5'-AU-3' motifs, respectively. In addition, we found that SARS-CoV-2 within-host variations are under negative selection. About 15.6% iSNVs had an impact on the content of the CpG dinucleotide (CpG) in SARS-CoV-2 genomes. We detected signatures of faster loss of CpG-gaining iSNVs, possibly resulting from zinc-finger antiviral protein-mediated antiviral activities targeting CpG, which could be the major reason for CpG depletion in SARS-CoV-2 consensus genomes. The non-synonymous iSNVs in the S gene can largely alter the S protein's antigenic features, and many of these iSNVs are distributed in the amino-terminal domain (NTD) and receptor-binding domain (RBD). These results suggest that SARS-CoV-2 interacts actively with human hosts and attempts to take different evolutionary strategies to escape human innate and adaptive immunity. These new findings further deepen and widen our understanding of the within-host evolutionary features of SARS-CoV-2. IMPORTANCE Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the causative pathogen of the coronavirus disease 2019, has evolved rapidly since it was discovered. Recent studies have pointed out that some mutations in the SARS-CoV-2 S protein could confer SARS-CoV-2 the ability to evade the human adaptive immune system. In addition, it is observed that the content of the CpG dinucleotide in SARS-CoV-2 genome sequences has decreased over time, reflecting the adaptation to the human host. The significance of our research is revealing the characteristics of SARS-CoV-2 within-host diversity of human hosts, identifying the causes of CpG depletion in SARS-CoV-2 consensus genomes, and exploring the potential impacts of non-synonymous within-host variations in the S gene on immune escape, which could further deepen and widen our understanding of the evolutionary features of SARS-CoV-2.

The how of counter-defense: viral evolution to combat host immunity

Viruses are locked in an evolutionary arms race with their hosts. What ultimately determines viral evolvability, or capacity for adaptive evolution, is their ability to efficiently explore and expand sequence space while under the selective regime imposed by their ecology, which includes innate and adaptive host defenses. Viral genomes have significantly higher evolutionary rates than their host counterparts and should have advantages relative to their slower-evolving hosts. However, functional constraints on virus evolutionary landscapes along with the modularity and mutational tolerance of host defense proteins may help offset the advantage conferred to viruses by high evolutionary rates. Additionally, cellular life forms from all domains of life possess many highly complex defense mechanisms that act as hurdles to viral replication. Consequently, viruses constantly probe sequence space through mutation and genetic exchange and are under pressure to optimize diverse counter-defense strategies.

Selective Depletion of ZAP-Binding CpG Motifs in HCV Evolution

Hepatitis C virus (HCV) is a bloodborne pathogen that can cause chronic liver disease and hepatocellular carcinoma. The loss of CpGs from virus genomes allows escape from restriction by the host zinc-finger antiviral protein (ZAP). The evolution of HCV in the human host has not been explored in the context of CpG depletion. We analysed 2616 full-length HCV genomes from 1977 to 2021. During the four decades of evolution in humans, we found that HCV genomes have become significantly depleted in (a) CpG numbers, (b) CpG O/E ratios (i.e., relative abundance of CpGs), and (c) the number of ZAP-binding motifs. Interestingly, our data suggests that the loss of CpGs in HCV genomes over time is primarily driven by the loss of ZAP-binding motifs; thus suggesting a yet unknown role for ZAP-mediated selection pressures in HCV evolution. The HCV core gene is significantly enriched for the number of CpGs and ZAP-binding motifs. In contrast to the rest of the HCV genome, the loss of CpGs from the core gene does not appear to be driven by ZAP-mediated selection. This work highlights CpG depletion in HCV genomes during their evolution in humans and the role of ZAP-mediated selection in HCV evolution.

Massive-scale genomic analysis reveals SARS-CoV-2 mutation characteristics and evolutionary trends

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) pandemic resulted in significant societal costs. Hence, an in-depth understanding of SARS-CoV-2 virus mutation and its evolution will help determine the direction of the COVID-19 pandemic. In this study, we identified 296,728 de novo mutations in more than 2,800,000 high-quality SARS-CoV-2 genomes. All possible factors affecting the mutation frequency of SARS-CoV-2 in human hosts were analyzed, including zinc finger antiviral proteins, sequence context, amino acid change, and translation efficiency. As a result, we proposed that when adenine (A) and tyrosine (T) bases are in the context of AM (M stands for adenine or cytosine) or TA motif, A or T base has lower mutation frequency. Furthermore, we hypothesized that translation efficiency can affect the mutation frequency of the third position of the codon by the selection, which explains why SARS-CoV-2 prefers AT3 codons usage. In addition, we found a host-specific asymmetric dinucleotide mutation frequency in the SARS-CoV-2 genome, which provides a new basis for determining the origin of the SARS-CoV-2. Finally, we summarize all possible factors affecting mutation frequency and provide insights into the mutation characteristics and evolutionary trends of SARS-CoV-2.

CpG content in the Zika virus genome affects infection phenotypes in the adult brain and fetal lymph nodes

Increasing the number of CpG dinucleotides in RNA viral genomes, while preserving the original amino acid composition, leads to impaired infection which does not cause disease. Beneficially, impaired infection evokes antiviral host immune responses providing a cutting-edge vaccine approach. For example, we previously showed that CpG-enriched Zika virus variants cause attenuated infection phenotypes and protect against lethal challenge in mice. While CpG recoding is an emerging and promising vaccine approach, little is known about infection phenotypes caused by recoded viruses in vivo, particularly in non-rodent species. Here, we used well-established mouse and porcine models to study infection phenotypes of the CpG-enriched neurotropic and congenital virus—Zika virus, directly in the target tissues—the brain and placenta. Specifically, we used the uttermost challenge and directly injected mice intracerebrally to compare infection phenotypes caused by wild-type and two CpG-recoded Zika variants and model the scenario where vaccine strains breach the blood-brain barrier. Also, we directly injected porcine fetuses to compare in utero infection phenotypes and model the scenario where recoded vaccine strains breach the placental barrier. While overall infection kinetics were comparable between wild-type and recoded virus variants, we found convergent phenotypical differences characterized by reduced pathology in the mouse brain and reduced replication of CpG-enriched variants in fetal lymph nodes. Next, using next-generation sequencing for the whole virus genome, we compared the stability of de novo introduced CpG dinucleotides during prolonged virus infection in the brain and placenta. Most de novo introduced CpG dinucleotides were preserved in sequences of recoded Zika viruses showing the stability of vaccine variants. Altogether, our study emphasized further directions to fine-tune the CpG recoding vaccine approach for better safety and can inform future immunization strategies.

Interaction between Sars-CoV-2 structural proteins and host cellular receptors: From basic mechanisms to clinical perspectives

Severe acute respiratory syndrome coronavirus 2 (Sars-CoV-2) has caused a global pandemic that has affected the lives of billions of individuals. Sars-CoV-2 primarily infects human cells by binding of the viral spike protein to angiotensin-converting enzyme 2 (ACE2). In addition, novel means of viral entry are currently being investigated, including Neuropillin 1, toll-like receptors (TLRs), cluster of differentiation 147 (CD147), and integrin α5β1. Enriched expression of these proteins across metabolic regulatory organs/tissues, including the circulatory system, liver, pancreas, and intestine contributes to major clinical complications among COVID-19 patients, particularly the development of hypertension, myocardial injury, arrhythmia, acute coronary syndrome and increased coagulation in the circulatory system during and post-infection. Pre-existing metabolic disease, such as cardiovascular disease, obesity, diabetes, and non-alcoholic fatty liver disease, is associated with increased risk of hospitalization, persistent post-infection complications and worse outcomes in patients with COVID-19. This review overviews the biological features of Sars-CoV-2, highlights recent findings that delineate the pathological mechanisms of COVID-19 and the consequent clinical diseases.