Evolutionary Basis of Codon Usage and Nucleotide Composition Bias in Vertebrate DNA Viruses

Understanding the extent and causes of biases in codon usage and nucleotide composition is essential to the study of viral evolution, particularly the interplay between viruses and host cells or immune responses. To understand the common features and differences among viruses we analyzed the genomic characteristics of a representative collection of all sequenced vertebrate-infecting DNA viruses. This revealed that patterns of codon usage bias are strongly correlated with overall genomic GC content, suggesting that genome-wide mutational pressure, rather than natural selection for specific coding triplets, is the main determinant of codon usage. Further, we observed a striking difference in CpG content between DNA viruses with large and small genomes. While the majority of large genome viruses show the expected frequency of CpG, most small genome viruses had CpG contents far below expected values. The exceptions to this generalization, the large gammaherpesviruses and iridoviruses and the small dependoviruses, have sufficiently different life-cycle characteristics that they may help reveal some of the factors shaping the evolution of CpG usage in viruses. Electronic Supplementary Material Electronic Supplementary material is available for this article at and accessible for authorised users. [Reviewing Editor: Dr. Nicolas Galtier]     

Emergence,evolution, and vaccine production approaches of SARS-CoV-2 virus: Benefits of getting vaccinated and common questions

The emergence of coronavirus disease 2019 (COVID-19) pandemic in Wuhan city, China at the end of 2019 made it urgent to identify the origin of the causal pathogen and its molecular evolution, to appropriately design an effective vaccine. This study analyzes the evolutionary background of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2 or SARS-2) in accordance with its close relative SARS-CoV (SARS-1), which was emerged in 2002. A comparative genomic and proteomic study was conducted on SARS-2, SARS-1, and Middle East respiratory syndrome coronavirus (MERS), which was emerged in 2012. In silico analysis inferred the genetic variability among the tested viruses. The SARS-1 genome harbored 11 genes encoding 12 proteins, while SARS-2 genome contained only 10 genes encoding for 10 proteins. MERS genome contained 11 genes encoding 11 proteins. The analysis also revealed a slight variation in the whole genome size of SARS-2 comparing to its siblings resulting from sequential insertions and deletions (indels) throughout the viral genome particularly ORF1AB, spike, ORF10 and ORF8. The effective indels were observed in the gene encoding the spike protein that is responsible for viral attachment to the angiotensin-converting enzyme 2 (ACE2) cell receptor and initiating infection. These indels are responsible for the newly emerging COVID-19 variants αCoV, βCoV, γCoV and δCoV. Nowadays, few effective COVID-19 vaccines developed based on spike (S) glycoprotein were approved and become available worldwide. Currently available vaccines can relatively prevent the spread of COVID-19 and suppress the disease. The traditional (killed or attenuated virus vaccine and antibody-based vaccine) and innovated vaccine production technologies (RNA- and DNA-based vaccines and viral vectors) are summarized in this review. We finally highlight the most common questions related to COVID-19 disease and the benefits of getting vaccinated.     

Genome composition and genetic characterization of SARS-CoV-2

SARS-CoV-2 is a type of Betacoronaviruses responsible for COVID-19 pandemic disease, with more than 1.745 million fatalities globally as of December-2020. Genetically, it is considered the second largest genome of all RNA viruses with a 5′ cap and 3′ poly-A tail. Phylogenetic analyses of coronaviruses reveal that SARS-CoV-2 is genetically closely related to the Bat-SARS Like-Corona virus (Bat-SL-Cov) with 96% whole-genome identity. SARS-CoV-2 genome consists of 15 ORFs coded into 29 proteins. At the 5′ terminal of the genome, we have ORF1ab and ORF1a, which encode the 1ab and 1a polypeptides that are proteolytically cleaved into 16 different nonstructural proteins (NSPs). The 3′ terminal of the genome represents four structural (spike, envelope, matrix, and nucleocapsid) and nine accessory (3a, 3b, 6, 7a, 7b, 8b, 9a, 9b, and orf10) proteins. As the number of COVID-19 patients increases dramatically worldwide, there is an urgent need to find a quick and sensitive diagnostic tool for controlling the outbreak of SARS-CoV-2 in the community. Today, molecular testing methods utilizing viral genetic material (e.g., PCR) represent the crucial diagnostic tool for the SARS-CoV-2 virus despite its low sensitivity in the early stage of viral infection. This review summarizes the genome composition and genetic characterization of the SARS-CoV-2.     

Signatures in SARS-CoV-2 spike protein conferring escape to neutralizing antibodies

Understanding SARS-CoV-2 evolution and host immunity is critical to control COVID-19 pandemics. At the core is an arms-race between SARS-CoV-2 antibody and angiotensin-converting enzyme 2 (ACE2) recognition, a function of the viral protein spike. Mutations in spike impacting antibody and/or ACE2 binding are appearing worldwide, imposing the need to monitor SARS-CoV2 evolution and dynamics in the population. Determining signatures in SARS-CoV-2 that render the virus resistant to neutralizing antibodies is critical. We engineered 25 spike-pseudotyped lentiviruses containing individual and combined mutations in the spike protein, including all defining mutations in the variants of concern, to identify the effect of single and synergic amino acid substitutions in promoting immune escape. We confirmed that E484K evades antibody neutralization elicited by infection or vaccination, a capacity augmented when complemented by K417N and N501Y mutations. In silico analysis provided an explanation for E484K immune evasion. E484 frequently engages in interactions with antibodies but not with ACE2. Importantly, we identified a novel amino acid of concern, S494, which shares a similar pattern. Using the already circulating mutation S494P, we found that it reduces antibody neutralization of convalescent and post-immunization sera, particularly when combined with E484K and with mutations able to increase binding to ACE2, such as N501Y. Our analysis of synergic mutations provides a signature for hotspots for immune evasion and for targets of therapies, vaccines and diagnostics.     

H1N1流感病毒聚合酶片段的CpG抑制及其对密码子偏爱性的影响

H1N1流感病毒的聚合酶具有RNA复制、转录等功能,并且对流感病毒片段包装、子代繁殖及宿主适应性等有着重要作用。通过分析人、猪及禽类H1N1流感病毒聚合酶片段的二核苷酸频率及同义密码子的偏爱性,发现不同宿主中,流感病毒聚合酶片段的CpG频率最低,且均被强烈抑制;通过三类宿主间的比较发现,人流感病毒抑制最为强烈,且其CpG频率随年份呈下降趋势,但2009年毒株的CpG频率突然上升。比较同义密码子使用频率发现,含有CpG的同义密码子相对使用频率均小于1,证明CpG抑制作用是影响同义密码子偏爱性的一个重要因素。以上结果暗示,CpG抑制对H1N1流感病毒的进化及跨宿主传播可能有重要影响。     

SARS-CoV-2 nucleocapsid protein forms condensates with viral genomic RNA

The Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) infection causes Coronavirus Disease 2019 (COVID-19), a pandemic that seriously threatens global health. SARS-CoV-2 propagates by packaging its RNA genome into membrane enclosures in host cells. The packaging of the viral genome into the nascent virion is mediated by the nucleocapsid (N) protein, but the underlying mechanism remains unclear. Here, we show that the N protein forms biomolecular condensates with viral genomic RNA both in vitro and in mammalian cells. While the N protein forms spherical assemblies with homopolymeric RNA substrates that do not form base pairing interactions, it forms asymmetric condensates with viral RNA strands. Cross-linking mass spectrometry (CLMS) identified a region that drives interactions between N proteins in condensates, and deletion of this region disrupts phase separation. We also identified small molecules that alter the size and shape of N protein condensates and inhibit the proliferation of SARS-CoV-2 in infected cells. These results suggest that the N protein may utilize biomolecular condensation to package the SARS-CoV-2 RNA genome into a viral particle.

The packaging of the SARS-CoV-2 genome is mediated by the nucleocapsid (N) protein; this study shows that the N protein forms liquid condensates with viral genomic RNA and identifies small molecules that alter these condensates.     

Dynamics of severe acute respiratory syndrome coronavirus 2 genome variants in the feces during convalescence

In response to the current coronavirus disease 2019 (COVID-19) pandemic, it is crucial to understand the origin, transmission, and evolution of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), which relies on close surveillance of genomic diversity in clinical samples. Although the mutation at the population level had been extensively investigated, how the mutations evolve at the individual level is largely unknown. Eighteen time-series fecal samples were collected from nine patients with COVID-19 during the convalescent phase. The nucleic acids of SARS-CoV-2 were enriched by the hybrid capture method. First, we demonstrated the outstanding performance of the hybrid capture method in detecting intra-host variants. We identified 229 intra-host variants at 182 sites in 18 fecal samples. Among them, nineteen variants presented frequency changes > 0.3 within 1–5 days, reflecting highly dynamic intra-host viral populations. Moreover, the evolution of the viral genome demonstrated that the virus was probably viable in the gastrointestinal tract during the convalescent period. Meanwhile, we also found that the same mutation showed a distinct pattern of frequency changes in different individuals, indicating a strong random drift. In summary, dramatic changes of the SARS-CoV-2 genome were detected in fecal samples during the convalescent period; whether the viral load in feces is sufficient to establish an infection warranted further investigation.     

The Slowing Rate of CpG Depletion in SARS-CoV-2 Genomes Is Consistent with Adaptations to the Human Host

Depletion of CpG dinucleotides in severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) genomes has been linked to virus evolution, host-switching, virus replication, and innate immune responses. Temporal variations, if any, in the rate of CpG depletion during virus evolution in the host remain poorly understood. Here, we analyzed the CpG content of over 1.4 million full-length SARS-CoV-2 genomes representing over 170 million documented infections during the first 17 months of the pandemic. Our findings suggest that the extent of CpG depletion in SARS-CoV-2 genomes is modest. Interestingly, the rate of CpG depletion is highest during early evolution in humans and it gradually tapers off, almost reaching an equilibrium; this is consistent with adaptations to the human host. Furthermore, within the coding regions, CpG depletion occurs predominantly at codon positions 2-3 and 3-1. Loss of ZAP (Zinc-finger antiviral protein)-binding motifs in SARS-CoV-2 genomes is primarily driven by the loss of the terminal CpG within the motifs. Nonetheless, majority of the CpG depletion in SARS-CoV-2 genomes occurs outside ZAP-binding motifs. SARS-CoV-2 genomes selectively lose CpGs-motifs from a U-rich context; this may help avoid immune recognition by TLR7. SARS-CoV-2 alpha-, beta-, and delta-variants of concern have reduced CpG content compared to sequences from the beginning of the pandemic. In sum, we provide evidence that the rate of CpG depletion in virus genomes is not uniform and it greatly varies over time and during adaptations to the host. This work highlights how temporal variations in selection pressures during virus adaption may impact the rate and the extent of CpG depletion in virus genomes.     

Vanishing GC-rich isochores in mammalian genomes

To understand the origin and evolution of isochores-the peculiar spatial distribution of GC content within mammalian genomes-we analyzed the synonymous substitution pattern in coding sequences from closely related species in different mammalian orders. In primate and cetartiodactyls, GC-rich genes are undergoing a large excess of GC --> AT substitutions over AT --> GC substitutions: GC-rich isochores are slowly disappearing from the genome of these two mammalian orders. In rodents, our analyses suggest both a decrease in GC content of GC-rich isochores and an increase in GC-poor isochores, but more data will be necessary to assess the significance of this pattern. These observations question the conclusions of previous works that assumed that base composition was at equilibrium. Analysis of allele frequency in human polymorphism data, however, confirmed that in the GC-rich parts of the genome, GC alleles have a higher probability of fixation than AT alleles. This fixation bias appears not strong enough to overcome the large excess of GC --> AT mutations. Thus, whatever the evolutionary force (neutral or selective) at the origin of GC-rich isochores, this force is no longer effective in mammals. We propose a model based on the biased gene conversion hypothesis that accounts for the origin of GC-rich isochores in the ancestral amniote genome and for their decline in present-day mammals.     

Looking for pathways related to COVID-19: confirmation of pathogenic mechanisms by SARS-CoV-2–host interactome

In the last months, many studies have clearly described several mechanisms of SARS-CoV-2 infection at cell and tissue level, but the mechanisms of interaction between host and SARS-CoV-2, determining the grade of COVID-19 severity, are still unknown. We provide a network analysis on protein–protein interactions (PPI) between viral and host proteins to better identify host biological responses, induced by both whole proteome of SARS-CoV-2 and specific viral proteins. A host-virus interactome was inferred, applying an explorative algorithm (Random Walk with Restart, RWR) triggered by 28 proteins of SARS-CoV-2. The analysis of PPI allowed to estimate the distribution of SARS-CoV-2 proteins in the host cell. Interactome built around one single viral protein allowed to define a different response, underlining as ORF8 and ORF3a modulated cardiovascular diseases and pro-inflammatory pathways, respectively. Finally, the network-based approach highlighted a possible direct action of ORF3a and NS7b to enhancing Bradykinin Storm. This network-based representation of SARS-CoV-2 infection could be a framework for pathogenic evaluation of specific clinical outcomes. We identified possible host responses induced by specific proteins of SARS-CoV-2, underlining the important role of specific viral accessory proteins in pathogenic phenotypes of severe COVID-19 patients.Subject terms: Protein-protein interaction networks, Viral infection     

The SARS-CoV-2 ORF10 is not essential in vitro or in vivo in humans

SARS-CoV-2 genome annotation revealed the presence of 10 open reading frames (ORFs), of which the last one (ORF10) is positioned downstream of the N gene. It is a hypothetical gene, which was speculated to encode a 38 aa protein. This hypothetical protein does not share sequence similarity with any other known protein and cannot be associated with a function. While the role of this ORF10 was proposed, there is growing evidence showing that the ORF10 is not a coding region. Here, we identified SARS-CoV-2 variants in which the ORF10 gene was prematurely terminated. The disease was not attenuated, and the transmissibility between humans was maintained. Also, in vitro, the strains replicated similarly to the related viruses with the intact ORF10. Altogether, based on clinical observation and laboratory analyses, it appears that the ORF10 protein is not essential in humans. This observation further proves that the ORF10 should not be treated as the protein-coding gene, and the genome annotations should be amended.     

The landscape of antibody binding in SARS-CoV-2 infection

The search for potential antibody-based diagnostics, vaccines, and therapeutics for pandemic severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has focused almost exclusively on the spike (S) and nucleocapsid (N) proteins. Coronavirus membrane (M), ORF3a, and ORF8 proteins are humoral immunogens in other coronaviruses (CoVs) but remain largely uninvestigated for SARS-CoV-2. Here, we use ultradense peptide microarray mapping to show that SARS-CoV-2 infection induces robust antibody responses to epitopes throughout the SARS-CoV-2 proteome, particularly in M, in which 1 epitope achieved excellent diagnostic accuracy. We map 79 B cell epitopes throughout the SARS-CoV-2 proteome and demonstrate that antibodies that develop in response to SARS-CoV-2 infection bind homologous peptide sequences in the 6 other known human CoVs. We also confirm reactivity against 4 of our top-ranking epitopes by enzyme-linked immunosorbent assay (ELISA). Illness severity correlated with increased reactivity to 9 SARS-CoV-2 epitopes in S, M, N, and ORF3a in our population. Our results demonstrate previously unknown, highly reactive B cell epitopes throughout the full proteome of SARS-CoV-2 and other CoV proteins.

Profiling of antibody binding from naïve and COVID-19 convalescent human sera to the entire proteome of SARS-CoV-2 and other human, bat and pangolin coronaviruses identifies 79 B cell epitopes throughout the SARS-CoV-2 proteome, finding that the most sensitive and specific binding occurred in the membrane (M) protein, and revealing cross-reactivity patterns.     

In vitro evolution predicts emerging SARS-CoV-2 mutations with high affinity for ACE2 and cross-species binding

Emerging SARS-CoV-2 variants are creating major challenges in the ongoing COVID-19 pandemic. Being able to predict mutations that could arise in SARS-CoV-2 leading to increased transmissibility or immune evasion would be extremely valuable in development of broad-acting therapeutics and vaccines, and prioritising viral monitoring and containment. Here we use in vitro evolution to seek mutations in SARS-CoV-2 receptor binding domain (RBD) that would substantially increase binding to ACE2. We find a double mutation, S477N and Q498H, that increases affinity of RBD for ACE2 by 6.5-fold. This affinity gain is largely driven by the Q498H mutation. We determine the structure of the mutant-RBD:ACE2 complex by cryo-electron microscopy to reveal the mechanism for increased affinity. Addition of Q498H to SARS-CoV-2 RBD variants is found to boost binding affinity of the variants for human ACE2 and confer a new ability to bind rat ACE2 with high affinity. Surprisingly however, in the presence of the common N501Y mutation, Q498H inhibits binding, due to a clash between H498 and Y501 side chains. To achieve an intermolecular bonding network, affinity gain and cross-species binding similar to Q498H alone, RBD variants with the N501Y mutation must acquire instead the related Q498R mutation. Thus, SARS-CoV-2 RBD can access large affinity gains and cross-species binding via two alternative mutational routes involving Q498, with route selection determined by whether a variant already has the N501Y mutation. These mutations are now appearing in emerging SARS-CoV-2 variants where they have the potential to influence human-to-human and cross-species transmission.     

Natural selection on synonymous sites is correlated with gene length and recombination in Drosophila.

Evolutionary analysis of codon bias in Drosophila indicates that synonymous mutations are not neutral, but rather are subject to weak selection at the translation level. Here we show that the effectiveness of natural selection on synonymous sites is strongly correlated with the rate of recombination, in accord with the nearly neutral hypothesis. This correlation, however, is apparent only in genes encoding short proteins. Long coding regions have both a lower codon bias and higher synonymous substitution rates, suggesting that they are affected less efficiently by selection. Therefore, both the length of the coding region and the recombination rate modulate codon bias. In addition, the data indicate that selection coefficients for synonymous mutations must vary by a minimum of one or two orders of magnitude. Two hypotheses are proposed to explain the relationship among the coding region length, the codon bias, and the synonymous divergence and polymorphism levels across the range of recombination rates in Drosophila. The first hypothesis is that selection coefficients on synonymous mutations are inversely related to the total length of the coding region. The second hypothesis proposes that interference among synonymous mutations reduces the efficacy of selection on these mutations. We investigated this second hypothesis by carrying out forward simulations of weakly selected mutations in model populations. These simulations show that even with realistic recombination rates, this interference, which we call the "small-scale" Hill-Robertson effect, can have a moderately strong influence on codon bias.     

Non-synonymous mutations of SARS-CoV-2 leads epitope loss and segregates its variants

The non-synonymous mutations of SARS-CoV-2 isolated from across the world have been identified during the last few months. The surface glycoprotein spike of SARS-CoV-2 forms the most important hotspot for amino acid alterations followed by the ORF1a/ORF1ab poly-proteins. It is evident that the D614G mutation in spike glycoprotein and P4715L in RdRp is the important determinant of SARS-CoV-2 evolution since its emergence. P4715L in RdRp, G251V in ORF3a and S1498F of Nsp3 is associated with the epitope loss that may influence pathogenesis caused by antibody escape variants. The phylogenomics distinguished the ancestral viral samples from China and most part of Asia, isolated since the initial outbreak and the later evolved variants isolated from Europe and Americas. The evolved variants have been found to predominant globally with the loss of epitopes from its proteins. These have implications for SARS-CoV-2 transmission, pathogenesis and immune interventions.     

B.1.617.2 (Delta) Variant of SARS-CoV-2: features,transmission and potential strategies

Coronavirus disease 2019 (COVID-19) caused by the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has become a pandemic. With the continuous evolution of the viral genome, SARS-CoV-2 has evolved many variants. B.1.617.2, also called Delta, is one of the most concerned variants. The Delta variant was first reported in India at the end of 2020 but has spread globally, by now, to 135 countries and is not stand still. Delta shared some mutations with other variants, and owned its special mutations on spike proteins, which may be responsible for its strong transmission and increasing virulence. Under these circumstances, a systematic summary of Delta is necessary. This review will focus on the Delta variant. We will describe all the characteristics of Delta (including biological features and clinical characteristics), analyze potential reasons for its strong transmission, and provide potential protective ways for combating Delta.