Genetic inactivation of poliovirus infectivity by increasing the frequencies of CpG and UpA dinucleotides within and across synonymous capsid region codons

Dengue virus preferentially uses human and mosquito non-optimal codons

Codon optimality refers to the effect that codon composition has on messenger RNA (mRNA) stability and translation level and implies that synonymous codons are not silent from a regulatory point of view. Here, we investigated the adaptation of virus genomes to the host optimality code using mosquito-borne dengue virus (DENV) as a model. We demonstrated that codon optimality exists in mosquito cells and showed that DENV preferentially uses nonoptimal (destabilizing) codons and avoids codons that are defined as optimal (stabilizing) in either human or mosquito cells. Human genes enriched in the codons preferentially and frequently used by DENV are upregulated during infection, and so is the tRNA decoding the nonoptimal and DENV preferentially used codon for arginine. We found that adaptation during single-host passaging in human or mosquito cells results in the selection of synonymous mutations towards DENV's preferred nonoptimal codons that increase virus fitness. Finally, our analyses revealed that hundreds of viruses preferentially use nonoptimal codons, with those infecting a single host displaying an even stronger bias, suggesting that host-pathogen interaction shapes virus-synonymous codon choice.

Application of MPBT Assay for Multiplex Determination of Infectious Titers and for Selection of the Optimal Formulation for the Trivalent Novel Oral Poliovirus Vaccine

Recently, a multiplex PCR-based titration (MPBT) assay was developed for simultaneous determination of infectious titers of all three Sabin strains of the oral poliovirus vaccine (OPV) to replace the conventional CCID50 assay, which is both time-consuming and laborious. The MPBT assay was shown to be reproducible, robust and sensitive. The conventional and MPBT assays showed similar results and sensitivity. The MPBT assay can be completed in two to three days, instead of ten days for the conventional assay. To prevent attenuated vaccine strains of poliovirus from reversion to virulence, a novel, genetically stable OPV (nOPV) was developed by modifying the genomes of conventional Sabin strains used in OPV. In this work, we evaluated the MPBT assay as a rapid screening tool to support trivalent nOPV (tnOPV) formulation development by simultaneous titration of the three nOPV strains to confirm stability as needed, for the selection of the lead tnOPV formulation candidate. We first assessed the ability of the MPBT assay to discriminate a 0.5 log10 titer difference by titrating the two tnOPV samples (undiluted and threefold-diluted) on the same plate. Once the assay was shown to be discriminating, we then tested different formulations of tnOPV drug products (DPs) that were subjected to different exposure times at 37 °C (untreated group and treated groups: 2 and 7 days at 37 °C), and to three freeze and thaw (FT) cycles. Final confirmation of the down selected formulation candidates was achieved by performing the conventional CCID50 assay, comparing the stability of untreated and treated groups and FT stability testing on the top three candidates. The results showed that the MPBT assay generates similar titers as the conventional assay. By testing two trivalent samples in the same plate, the assay can differentiate a 0.5 log10 difference between the titers of the tested nOPV samples. Also, the assay was able to detect the gradual degradation of nOPV viruses with different formulation compositions and under different time/temperature conditions and freeze/thaw cycles. We found that there were three tnOPV formulations which met the stability criteria of less than 0.5 log10 loss after 2 days' exposure to 37 ℃ and after three FT cycles, maintaining the potency of all three serotypes in these formulations. The ability of the MPBT assay to titrate two tnOPV lots (six viruses) in the same plate makes it cheaper and gives it a higher throughput for rapid screening. The assay detected the gradual degradation of the tnOPV and was successful in the selection of optimal formulations for the tnOPV. The results demonstrated that the MPBT method can be used as a stability indicating assay to assess the thermal stability of the nOPV. It can be used for rapid virus titer determination during the vaccine manufacturing process, and in clinical trials. The MPBT assay can be automated and applied for other viruses, including those with no cytopathic effect.

Codon usage bias of human papillomavirus type 33 and 58: A comprehensive analysis

Cervical cancer is closely linked to specific strains of human papillomavirus (HPV), notably HPV-33 and HPV-58, which exhibit a significant prevalence among women in China. Nevertheless, the codon usage bias in HPV-33 and HPV-58 is not well comprehended. The objective of this research is to analyze the codon usage patterns HPV-33 and HPV-58, pinpoint the primary factors that influence codon preference. The overall preference for codon usage in two HPV genotypes is not significant. Both HPV genotypes exhibit a preference for codons that end with A/U. The GC3 content for HPV-33 is 25.43% ± 0.35%, and for HPV-58, it is 29.44% ± 0.57%. Out of the 26 favored codons in HPV-33 and HPV-58 (relative synonymous codon usage (RSCU) > 1), 25 conclude with A/U. Principal component analysis (PCA) shows a tight clustering of the entire genome sequences of HPV-33 and HPV-58, suggesting a similarity in their RSCU preferences. Moreover, an examination of dinucleotide abundance indicated that translation selection influenced the development of a distinctive dinucleotide usage pattern in HPV-33 and HPV-58. Additionally, a combined analysis involving an effective number of codons plot, parity rule 2, and neutrality analysis demonstrated that, for HPV-33 and HPV-58, the primary determinant influencing codon usage preference is natural selection. HPV-33 and HPV-58 exhibit a restricted set of favored codons in common with humans, potentially mitigating competition for translation resources. Our discoveries could provide valuable perspectives on the evolutionary patterns and codon usage preferences of HPV-33 and HPV-58 viruses, contributing to the development and application of relevant HPV subtype vaccines.

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.

Antiviral Activity of Zinc Finger Antiviral Protein (ZAP) in Different Virus Families

The CCCH-type zinc finger antiviral protein (ZAP) in humans, specifically isoforms ZAP-L and ZAP-S, is a crucial component of the cell's intrinsic immune response. ZAP acts as a post-transcriptional RNA restriction factor, exhibiting its activity during infections caused by retroviruses and alphaviruses. Its function involves binding to CpG (cytosine-phosphate-guanine) dinucleotide sequences present in viral RNA, thereby directing it towards degradation. Since vertebrate cells have a suppressed frequency of CpG dinucleotides, ZAP is capable of distinguishing foreign genetic elements. The expression of ZAP leads to the reduction of viral replication and impedes the assembly of new virus particles. However, the specific mechanisms underlying these effects have yet to be fully understood. Several questions regarding ZAP's mechanism of action remain unanswered, including the impact of CpG dinucleotide quantity on ZAP's activity, whether this sequence is solely required for the binding between ZAP and viral RNA, and whether the recruitment of cofactors is dependent on cell type, among others. This review aims to integrate the findings from studies that elucidate ZAP's antiviral role in various viral infections, discuss gaps that need to be filled through further studies, and shed light on new potential targets for therapeutic intervention.

Dinucleotide biases in RNA viruses that infect vertebrates or invertebrates

IMPORTANCE

Akin to a molecular signature, dinucleotide composition can be exploited by the zinc-finger antiviral protein (ZAP) to restrict CpG-rich (and UpA-rich) RNA viruses. ZAP evolved in tetrapods, and it is not encoded by invertebrates and fish. Because a systematic analysis is missing, we analyzed the genomes of RNA viruses that infect vertebrates or invertebrates. We show that vertebrate single-stranded (ss) RNA(+) viruses and, to a lesser extent, double-stranded RNA viruses tend to have stronger CpG bias than invertebrate viruses. Conversely, ssRNA(-) viruses have similar dinucleotide composition whether they infect vertebrates or invertebrates. Analysis of ssRNA(+) viruses that infect mammals, reptiles, and fish indicated that ZAP is unlikely to be a major driver of CpG depletion. We also show that, compared to other coronaviruses, the genome of SARS-CoV-2 is not homogeneously CpG-depleted. Our study provides new insights into virus evolution and strategies for recoding RNA virus genomes.

parazitCUB: An R package to streamline the process of investigating the adaptations of parasites' codon usage bias

Examining the intricate association between parasites and their hosts, particularly at the codon level, assumes paramount importance in comprehending evolutionary processes and forecasting the characteristics of novel parasites. While diverse metrics and statistical analyses are available to explore codon usage bias (CUB), there presently exists no dedicated tool for examining the co-adaptation of codon usage between parasites and hosts. Therefore, we introduce the parazitCUB R package to address this challenge in a scalable and efficient manner, as it is capable of handling extensive datasets and simultaneously analyzing of multiple parasites with optimized performance. parazitCUB enables the elucidation of parasite-host interactions and the evolutionary patterns of parasites through the implementation of various indices, cluster analysis, multivariate analysis, and data visualization techniques. The tool can be accessed at the following location: https://github.com/AliYoussef96/parazitCUB.

Evolution and diversity of nucleotide and dinucleotide composition in poxviruses

Poxviruses (family Poxviridae) have long dsDNA genomes and infect a wide range of hosts, including insects, birds, reptiles and mammals. These viruses have substantial incidence, prevalence and disease burden in humans and in other animals. Nucleotide and dinucleotide composition, mostly CpG and TpA, have been largely studied in viral genomes because of their evolutionary and functional implications. We analysed here the nucleotide and dinucleotide composition, as well as codon usage bias, of a set of representative poxvirus genomes, with a very diverse host spectrum. After correcting for overall nucleotide composition, entomopoxviruses displayed low overall GC content, no enrichment in TpA and large variation in CpG enrichment, while chordopoxviruses showed large variation in nucleotide composition, no obvious depletion in CpG and a weak trend for TpA depletion in GC-rich genomes. Overall, intergenome variation in dinucleotide composition in poxviruses is largely accounted for by variation in overall genomic GC levels. Nonetheless, using vaccinia virus as a model, we found that genes expressed at the earliest times in infection are more CpG-depleted than genes expressed at later stages. This observation has parallels in betahepesviruses (also large dsDNA viruses) and suggests an antiviral role for the innate immune system (e.g. via the zinc-finger antiviral protein ZAP) in the early phases of poxvirus infection. We also analysed codon usage bias in poxviruses and we observed that it is mostly determined by genomic GC content, and that stratification after host taxonomy does not contribute to explaining codon usage bias diversity. By analysis of within-species diversity, we show that genomic GC content is the result of mutational biases. Poxvirus genomes that encode a DNA ligase are significantly AT-richer than those that do not, suggesting that DNA repair systems shape mutation biases. Our data shed light on the evolution of poxviruses and inform strategies for their genetic manipulation for therapeutic purposes.

Evaluation of the immunization effectiveness of bOPV booster immunization and IPV revaccination

To provide a basis for further optimization of the polio sequential immunization schedule, this study evaluated the effectiveness of booster immunization with one dose of bivalent oral poliovirus vaccine (bOPV) at 48 months of age after different primary polio immunization schedules. At 48 months of age, one dose of bOPV was administered, and their poliovirus types 1-3 (PV1, PV2, and PV3, respectively)-specific neutralizing antibody levels were determined. Participants found to be negative for any type of PV-specific neutralizing antibody at 24, 36, or 48 months of age were re-vaccinated with inactivated polio vaccine (IPV). The 439 subjects who received a bOPV booster immunization at the age of 48 months had lower PV2-specific antibody levels compared with those who received IPV. One dose of IPV during basic polio immunization induced the lowest PV2-specific antibody levels. On the basis of our findings, to ensure that no less than 70% of the vaccinated have protection efficiency, we recommend the following: if basic immunization was conducted with 1IPV + 2bOPV (especially Sabin strain-based IPV), a booster immunization with IPV is recommended at 36 months of age, whereas if basic immunization was conducted with 2IPV + 1bOPV, a booster immunization with IPV is recommended at 48 months of age. A sequential immunization schedule of 2IPV + 1bOPV + 1IPV can not only maintain high levels of antibody against PV1 and PV3 but also increases immunity to PV2 and induces early intestinal mucosal immunity, with relatively good safety. Thus, this may be the best sequential immunization schedule for polio in countries or regions at high risk for polio.

Transcriptomic, proteomic, and functional consequences of codon usage bias in human cells during heterologous gene expression

Differences in codon frequency between genomes, genes, or positions along a gene, modulate transcription and translation efficiency, leading to phenotypic and functional differences. Here, we present a multiscale analysis of the effects of synonymous codon recoding during heterologous gene expression in human cells, quantifying the phenotypic consequences of codon usage bias at different molecular and cellular levels, with an emphasis on translation elongation. Six synonymous versions of an antibiotic resistance gene were generated, fused to a fluorescent reporter, and independently expressed in HEK293 cells. Multiscale phenotype was analyzed by means of quantitative transcriptome and proteome assessment, as proxies for gene expression; cellular fluorescence, as a proxy for single-cell level expression; and real-time cell proliferation in absence or presence of antibiotic, as a proxy for the cell fitness. We show that differences in codon usage bias strongly impact the molecular and cellular phenotype: (i) they result in large differences in mRNA levels and protein levels, leading to differences of over 15 times in translation efficiency; (ii) they introduce unpredicted splicing events; (iii) they lead to reproducible phenotypic heterogeneity; and (iv) they lead to a trade-off between the benefit of antibiotic resistance and the burden of heterologous expression. In human cells in culture, codon usage bias modulates gene expression by modifying mRNA availability and suitability for translation, leading to differences in protein levels and eventually eliciting functional phenotypic changes.

A Synthetic Biology Approach for Vaccine Candidate Design against Delta Strain of SARS-CoV-2 Revealed Disruption of Favored Codon Pair as a Better Strategy over Using Rare Codons

The SARS-CoV-2 delta variant (B.1.617.2) appeared for the first time in December 2020 and later spread worldwide. Currently available vaccines are not so efficacious in curbing the viral pathogenesis of the delta strain of COVID; therefore, the development of a safe and effective vaccine is required. In the present study, we envisaged molecular patterns in the structural genes' spike, nucleoprotein, membrane, and envelope of the SARS-CoV-2 delta variant. The study was based on determining compositional features, dinucleotide odds ratio, synonymous codon usage, positive and negative codon contexts, rare codons, and insight into relatedness between the human host isoacceptor tRNA and preferred codons from the structural genes. We found specific patterns, including a significant abundance of T nucleotide over all other three nucleotides. The underrepresentation of GpA, GpG, CpC, and CpG dinucleotides and the overrepresentation of TpT, ApA, CpT, and TpG were observed. A preference towards ACT- (Thr), AAT- (Asn), TTT- (Phe), and TTG- (Leu) initiated codons and aversion towards CGG (Arg), CCG (Pro), and CAC (His) was present in the structural genes of the delta strain. The interaction between the host tRNA pool and preferred codons of the envisaged structural genes revealed that the virus preferred the codons for those suboptimal numbers of isoacceptor tRNA were present. We see this as a strategy adapted by the virus to keep the translation rate low to facilitate the correct folding of viral proteins. The information generated in the study helps design the attenuated vaccine candidate against the SARS-CoV-2 delta variant using a synthetic biology approach. Three strategies were tested: changing TpT to TpA, introducing rare codons, and disrupting favored codons. It found that disrupting favored codons is a better approach to reducing virus fitness and attenuating SARS-CoV-2 delta strain using structural genes.

Evolutionary dynamics of codon usages for peste des petits ruminants virus

Peste des petits ruminants virus (PPRV) is an important agent of contagious, acute and febrile viral diseases in small ruminants, while its evolutionary dynamics related to codon usage are still lacking. Herein, we adopted information entropy, the relative synonymous codon usage values and similarity indexes and codon adaptation index to analyze the viral genetic features for 45 available whole genomes of PPRV. Some universal, lineage-specific, and gene-specific genetic features presented by synonymous codon usages of the six genes of PPRV that encode N, P, M, F, H and L proteins reflected evolutionary plasticity and independence. The high adaptation of PPRV to hosts at codon usages reflected high viral gene expression, but some synonymous codons that are rare in the hosts were selected in high frequencies in the viral genes. Another obvious genetic feature was that the synonymous codons containing CpG dinucleotides had weak tendencies to be selected in viral genes. The synonymous codon usage patterns of PPRV isolated during 2007–2008 and 2013–2014 in China displayed independent evolutionary pathway, although the overall codon usage patterns of these PPRV strains matched the universal codon usage patterns of lineage IV. According to the interplay between nucleotide and synonymous codon usages of the six genes of PPRV, the evolutionary dynamics including mutation pressure and natural selection determined the viral survival and fitness to its host.

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.

Intra-genomic heterogeneity in CpG dinucleotide composition in dengue virus

PURPOSE

Dengue virus is a life-threatening virus and cases of dengue infection have been increasing steadily in the past decades causing millions of deaths every year. So far, there is no vaccine that works effectively on all serotypes. Recently, CpG-recoded vaccines have proved to be effective against few viruses.

METHODS

In this study, evaluation and interpretation of more than 4547 Dengue virus genome sequences were included for analyzing novel CpG dinucleotides rich regions which are shared amid all serotypes. Genomic regions of DENV were synonymously CpG recoded using in silico methods and analyzed for adaptation in both human and Aedes spp. hosts based on CAI scores.

RESULTS

The analysis mirrored that serotypes 1, 3, and 4 shared CpG islands present in common regions. DENV-2 CpG islands showed no similarity with any of the CpG islands present in other serotypes. While DENV-3 sequences were found to possess the maximum number of conserved CpG islands stretches; DENV-2 was found to possess the lowest number. We found that all serotypes (with an exception of serotype 2) have CpG island in their 3' UTR. In silico CpG recoding of DENV genomic regions resulted in ∼ 3 fold increase of CpG dinucleotide frequency and comparative analysis based on CAI scores showed decreased adaptive fitness of CpG recoded DENV inside human host.

CONCLUSION

These CG-dinucleotide-enriched RNA sequences can be targeted by ZAP (zinc-finger antiviral protein) which can differentiate between host mRNA and viral mRNA. Our in silico findings can further be exploited for CpG-recoding of DENV genomes which can evoke cellular and humoral immune responses by recruiting ZAP-induced RNA degradation machinery and hence providing a promising approach for vaccine development.

Compositional biases in RNA viruses: Causes, consequences and applications

If each of the four nucleotides were represented equally in the genomes of viruses and the hosts they infect, each base would occur at a frequency of 25%. However, this is not observed in nature. Similarly, the order of nucleotides is not random (e.g., in the human genome, guanine follows cytosine at a frequency of ~0.0125, or a quarter the number of times predicted by random representation). Codon usage and codon order are also nonrandom. Furthermore, nucleotide and codon biases vary between species. Such biases have various drivers, including cellular proteins that recognize specific patterns in nucleic acids, that once triggered, induce mutations or invoke intrinsic or innate immune responses. In this review we examine the types of compositional biases identified in viral genomes and current understanding of the evolutionary mechanisms underpinning these trends. Finally, we consider the potential for large scale synonymous recoding strategies to engineer RNA virus vaccines, including those with pandemic potential, such as influenza A virus and Severe Acute Respiratory Syndrome Coronavirus Virus 2. This article is categorized under: RNA in Disease and Development > RNA in Disease RNA Evolution and Genomics > Computational Analyses of RNA RNA Interactions with Proteins and Other Molecules > Protein-RNA Recognition.

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

Depletion of CpG dinucleotides in 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 analysed 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-binding motifs in SARS-CoV-2 genomes is primarily driven by the loss of the terminal CpG in 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.

Generation of CpG-Recoded Zika Virus Vaccine Candidates

Experimental increase of cytosine-phosphate-guanine (CpG) dinucleotides in an RNA virus genome impairs infection. Beneficially, this weak infection may lead to robust antiviral host immunity providing a cutting-edge approach for vaccines. For example, we have recently demonstrated that recoded Zika virus variants with the increased CpG content showed considerable attenuated infection phenotypes and protection against lethal challenge in mice. Here, we describe the workflow for the design and generation of CpG-recoded Zika virus vaccine candidates. The workflow can be adapted for other viruses.

Variability in Codon Usage in Coronaviruses Is Mainly Driven by Mutational Bias and Selective Constraints on CpG Dinucleotide

The Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is the third human-emerged virus of the 21st century from the Coronaviridae family, causing the ongoing coronavirus disease 2019 (COVID-19) pandemic. Due to the high zoonotic potential of coronaviruses, it is critical to unravel their evolutionary history of host species breadth, host-switch potential, adaptation and emergence, to identify viruses posing a pandemic risk in humans. We present here a comprehensive analysis of the composition and codon usage bias of the 82 Orthocoronavirinae members, infecting 47 different avian and mammalian hosts. Our results clearly establish that synonymous codon usage varies widely among viruses, is only weakly dependent on their primary host, and is dominated by mutational bias towards AU-enrichment and by CpG avoidance. Indeed, variation in GC3 explains around 34%, while variation in CpG frequency explains around 14% of total variation in codon usage bias. Further insight on the mutational equilibrium within Orthocoronavirinae revealed that most coronavirus genomes are close to their neutral equilibrium, the exception being the three recently infecting human coronaviruses, which lie further away from the mutational equilibrium than their endemic human coronavirus counterparts. Finally, our results suggest that, while replicating in humans, SARS-CoV-2 is slowly becoming AU-richer, likely until attaining a new mutational equilibrium.

Assessment of genetic changes and neurovirulence of shed Sabin and novel type 2 oral polio vaccine viruses

Sabin-strain oral polio vaccines (OPV) can, in rare instances, cause disease in recipients and susceptible contacts or evolve to become circulating vaccine-derived strains with the potential to cause outbreaks. Two novel type 2 OPV (nOPV2) candidates were designed to stabilize the genome against the rapid reversion that is observed following vaccination with Sabin OPV type 2 (mOPV2). Next-generation sequencing and a modified transgenic mouse neurovirulence test were applied to shed nOPV2 viruses from phase 1 and 2 studies and shed mOPV2 from a phase 4 study. The shed mOPV2 rapidly reverted in the primary attenuation site (domain V) and increased in virulence. In contrast, the shed nOPV2 viruses showed no evidence of reversion in domain V and limited or no increase in neurovirulence in mice. Based on these results and prior published data on safety, immunogenicity, and shedding, the nOPV2 viruses are promising alternatives to mOPV2 for outbreak responses.

Development of a rabies virus-based retrograde tracer with high trans-monosynaptic efficiency by reshuffling glycoprotein

Rabies virus (RV) is the most widely used vector for mapping neural circuits. Previous studies have shown that the RV glycoprotein can be a target to improve the retrograde transsynaptic tracing efficiency. However, the current versions still label only a small portion of all presynaptic neurons. Here, we reshuffled the oG sequence, a chimeric glycoprotein, with positive codon pair bias score (CPBS) based on bioinformatic analysis of mouse codon pair bias, generating ooG, a further optimized glycoprotein. Our experimental data reveal that the ooG has a higher expression level than the oG in vivo, which significantly increases the tracing efficiency by up to 12.6 and 62.1-fold compared to oG and B19G, respectively. The new tool can be used for labeling neural circuits Therefore, the approach reported here provides a convenient, efficient and universal strategy to improve protein expression for various application scenarios such as trans-synaptic tracing efficiency, cell engineering, and vaccine and oncolytic virus designs.

Altering Compositional Properties of Viral Genomes to Design Live-Attenuated Vaccines

Live-attenuated vaccines have been historically used to successfully prevent numerous diseases caused by a broad variety of RNA viruses due to their ability to elicit strong and perdurable immune-protective responses. In recent years, various strategies have been explored to achieve viral attenuation by rational genetic design rather than using classic and empirical approaches, based on successive passages in cell culture. A deeper understanding of evolutionary implications of distinct viral genomic compositional aspects, as well as substantial advances in synthetic biology technologies, have provided a framework to achieve new viral attenuation strategies. Herein, we will discuss different approaches that are currently applied to modify compositional features of viruses in order to develop novel live-attenuated vaccines.

Evolving pharmacovigilance requirements with novel vaccines and vaccine components

This paper explores the pipeline of new and upcoming vaccines as it relates to monitoring their safety. Compared with most currently available vaccines, that are constituted of live attenuated organisms or inactive products, future vaccines will also be based on new technologies. Several products that include such technologies are either already licensed or at an advanced stage of clinical development. Those include viral vectors, genetically attenuated live organisms, nucleic acid vaccines, novel adjuvants, increased number of antigens present in a single vaccine, novel mode of vaccine administration and thermostabilisation. The Global Advisory Committee on Vaccine Safety (GACVS) monitors novel vaccines, from the time they become available for large scale use. GACVS maintains their safety profile as evidence emerges from post-licensure surveillance and observational studies. Vaccines and vaccine formulations produced with novel technologies will have different safety profiles that will require adapting pharmacovigilance approaches. For example, GACVS now considers viral vector templates developed on the model proposed by Brighton Collaboration. The characteristics of those novel products will also have implications for the risk management plans (RMPs). Questions related to the duration of active monitoring for genetic material, presence of adventitious agents more easily detected with enhanced biological screening, or physiological mechanisms of novel adjuvants are all considerations that will belong to the preparation of RMPs. In addition to assessing those novel products and advising experts, GACVS will also consider how to more broadly communicate about risk assessment, so vaccine users can also benefit from the committee’s advice.

The Evolution of Severe Acute Respiratory Syndrome Coronavirus-2 during Pandemic and Adaptation to the Host

Severe Acute Respiratory Syndrome Coronavirus-2 is a zoonotic virus with a possible origin in bats and potential transmission to humans through an intermediate host. When zoonotic viruses jump to a new host, they undergo both mutational and natural selective pressures that result in non-synonymous and synonymous adaptive changes, necessary for efficient replication and rapid spread of diseases in new host species. The nucleotide composition and codon usage pattern of SARS-CoV-2 indicate the presence of a highly conserved, gene-specific codon usage bias. The codon usage pattern of SARS-CoV-2 is mostly antagonistic to human and bat codon usage. SARS-CoV-2 codon usage bias is mainly shaped by the natural selection, while mutational pressure plays a minor role. The time-series analysis of SARS-CoV-2 genome indicates that the virus is slowly evolving. Virus isolates from later stages of the outbreak have more biased codon usage and nucleotide composition than virus isolates from early stages of the outbreak.

Zinc finger proteins in the host-virus interplay: multifaceted functions based on their nucleic acid-binding property

Zinc finger proteins (ZFPs) are a huge family comprised of massive, structurally diverse proteins characterized by zinc ion coordinating. They engage in the host-virus interplay in-depth and occupy a significant portion of the host antiviral arsenal. Nucleic acid-binding is the basic property of certain ZFPs, which draws increasing attention due to their immense influence on viral infections. ZFPs exert multiple roles on the viral replications and host cell transcription profiles by recognizing viral genomes and host mRNAs. Their roles could be either antiviral or proviral and were separately discussed. Our review covers the recent research progress and provides a comprehensive understanding of ZFPs in antiviral immunity based on their DNA/RNA binding property.

Read between the lines: Diversity of non-translational selection pressures on local codon usage

Protein coding genes can contain specific motifs within their nucleotide sequence that function as a signal for various biological pathways. The presence of such sequence motifs within a gene can have beneficial or detrimental effects on the phenotype and fitness of an organism, and this can lead to the enrichment or avoidance of this sequence motif. The degeneracy of the genetic code allows for the existence of alternative synonymous sequences that exclude or include these motifs, while keeping the encoded amino acid sequence intact. This implies that locally, there can be a selective pressure for preferentially using a codon over its synonymous alternative in order to avoid or enrich a specific sequence motif. This selective pressure could -in addition to mutation, drift and selection for translation efficiency and accuracy- contribute to shape the codon usage bias. In this review, we discuss patterns of avoidance of (or enrichment for) the various biological signals contained in specific nucleotide sequence motifs: transcription and translation initiation and termination signals, mRNA maturation signals, and antiviral immune system targets. Experimental data on the phenotypic or fitness effects of synonymous mutations in these sequence motifs confirm that they can be targets of local selection pressures on codon usage. We also formulate the hypothesis that transposable elements could have a similar impact on codon usage through their preferred integration sequences. Overall, selection on codon usage appears to be a combination of a global selection pressure imposed by the translation machinery, and a patchwork of local selection pressures related to biological signals contained in specific sequence motifs.

SARS-CoV-2 hot-spot mutations are significantly enriched within inverted repeats and CpG island loci

SARS-CoV-2 is an intensively investigated virus from the order Nidovirales (Coronaviridae family) that causes COVID-19 disease in humans. Through enormous scientific effort, thousands of viral strains have been sequenced to date, thereby creating a strong background for deep bioinformatics studies of the SARS-CoV-2 genome. In this study, we inspected high-frequency mutations of SARS-CoV-2 and carried out systematic analyses of their overlay with inverted repeat (IR) loci and CpG islands. The main conclusion of our study is that SARS-CoV-2 hot-spot mutations are significantly enriched within both IRs and CpG island loci. This points to their role in genomic instability and may predict further mutational drive of the SARS-CoV-2 genome. Moreover, CpG islands are strongly enriched upstream from viral ORFs and thus could play important roles in transcription and the viral life cycle. We hypothesize that hypermethylation of these loci will decrease the transcription of viral ORFs and could therefore limit the progression of the disease.

Impact of Synonymous Genome Recoding on the HIV Life Cycle

Synonymous mutations within protein coding regions introduce changes in DNA or messenger (m) RNA, without mutating the encoded proteins. Synonymous recoding of virus genomes has facilitated the identification of previously unknown virus biological features. Moreover, large-scale synonymous recoding of the genome of human immunodeficiency virus type 1 (HIV-1) has elucidated new antiviral mechanisms within the innate immune response, and has improved our knowledge of new functional virus genome structures, the relevance of codon usage for the temporal regulation of viral gene expression, and HIV-1 mutational robustness and adaptability. Continuous improvements in our understanding of the impacts of synonymous substitutions on virus phenotype - coupled with the decreased cost of chemically synthesizing DNA and improved methods for assembling DNA fragments - have enhanced our ability to identify potential HIV-1 and host factors and other aspects involved in the infection process. In this review, we address how silent mutagenesis impacts HIV-1 phenotype and replication capacity. We also discuss the general potential of synonymous recoding of the HIV-1 genome to elucidate unknown aspects of the virus life cycle, and to identify new therapeutic targets.

Host-virus interaction and viral evasion

With each infectious pandemic or outbreak, the medical community feels the need to revisit basic concepts of immunology to understand and overcome the difficult times brought about by these infections. Regarding viruses, they have historically been responsible for many deaths, and such a peculiarity occurs because they are known to be obligate intracellular parasites that depend upon the host's cell machinery for their replication. Successful infection with the production of essential viral components requires constant viral evolution as a strategy to manipulate the cellular environment, including host internal factors, the host's nonspecific and adaptive immune responses to viruses, the metabolic and energetic state of the infected cell, and changes in the intracellular redox environment during the viral infection cycle. Based on this knowledge, it is fundamental to develop new therapeutic strategies for controlling viral dissemination, by means of antiviral therapies, vaccines, or antioxidants, or by targeting the inhibition or activation of cell signaling pathways or metabolic pathways that are altered during infection. The rapid recovery of altered cellular homeostasis during viral infection is still a major challenge. Here, we review the strategies by which viruses evade the host's immune response and potential tools used to develop more specific antiviral therapies to cure, control, or prevent viral diseases.

Codon Usage and Adenovirus Fitness: Implications for Vaccine Development

Vaccination is the most effective method to date to prevent viral diseases. It intends to mimic a naturally occurring infection while avoiding the disease, exposing our bodies to viral antigens to trigger an immune response that will protect us from future infections. Among different strategies for vaccine development, recombinant vaccines are one of the most efficient ones. Recombinant vaccines use safe viral vectors as vehicles and incorporate a transgenic antigen of the pathogen against which we intend to generate an immune response. These vaccines can be based on replication-deficient viruses or replication-competent viruses. While the most effective strategy involves replication-competent viruses, they must be attenuated to prevent any health hazard while guaranteeing a strong humoral and cellular immune response. Several attenuation strategies for adenoviral-based vaccine development have been contemplated over time. In this paper, we will review them and discuss novel approaches based on the principle that protein synthesis from individual genes can be modulated by codon usage bias manipulation. We will summarize vaccine approaches that consider recoding of viral proteins to produce adenoviral attenuation and recoding of the transgene antigens for both viral attenuation and efficient viral epitope expression.

Biased Mutation and Selection in RNA Viruses

RNA viruses are responsible for some of the worst pandemics known to mankind, including outbreaks of Influenza, Ebola, and COVID-19. One major challenge in tackling RNA viruses is the fact they are extremely genetically diverse. Nevertheless, they share common features that include their dependence on host cells for replication, and high mutation rates. We set out to search for shared evolutionary characteristics that may aid in gaining a broader understanding of RNA virus evolution, and constructed a phylogeny-based data set spanning thousands of sequences from diverse single-stranded RNA viruses of animals. Strikingly, we found that the vast majority of these viruses have a skewed nucleotide composition, manifested as adenine rich (A-rich) coding sequences. In order to test whether A-richness is driven by selection or by biased mutation processes, we harnessed the effects of incomplete purifying selection at the tips of virus phylogenies. Our results revealed consistent mutational biases toward U rather than A in genomes of all viruses. In +ssRNA viruses, we found that this bias is compensated by selection against U and selection for A, which leads to A-rich genomes. In -ssRNA viruses, the genomic mutational bias toward U on the negative strand manifests as A-rich coding sequences, on the positive strand. We investigated possible reasons for the advantage of A-rich sequences including weakened RNA secondary structures, codon usage bias, and selection for a particular amino acid composition, and conclude that host immune pressures may have led to similar biases in coding sequence composition across very divergent RNA viruses.

Updated Characterization of Outbreak Response Strategies for 2019-2029: Impacts of Using a Novel Type 2 Oral Poliovirus Vaccine Strain

Delays in achieving the global eradication of wild poliovirus transmission continue to postpone subsequent cessation of all oral poliovirus vaccine (OPV) use. Countries must stop OPV use to end all cases of poliomyelitis, including vaccine-associated paralytic polio (VAPP) and cases caused by vaccine-derived polioviruses (VDPVs). The Global Polio Eradication Initiative (GPEI) coordinated global cessation of all type 2 OPV (OPV2) use in routine immunization in 2016 but did not successfully end the transmission of type 2 VDPVs (VDPV2s), and consequently continues to use type 2 OPV (OPV2) for outbreak response activities. Using an updated global poliovirus transmission and OPV evolution model, we characterize outbreak response options for 2019-2029 related to responding to VDPV2 outbreaks with a genetically stabilized novel OPV (nOPV2) strain or with the currently licensed monovalent OPV2 (mOPV2). Given uncertainties about the properties of nOPV2, we model different assumptions that appear consistent with the evidence on nOPV2 to date. Using nOPV2 to respond to detected cases may reduce the expected VDPV and VAPP cases and the risk of needing to restart OPV2 use in routine immunization compared to mOPV2 use for outbreak response. The actual properties, availability, and use of nOPV2 will determine its effects on type 2 poliovirus transmission in populations. Even with optimal nOPV2 performance, countries and the GPEI would still likely need to restart OPV2 use in routine immunization in OPV-using countries if operational improvements in outbreak response to stop the transmission of cVDPV2s are not implemented effectively.

Use of Synonymous Deoptimization to Derive Modified Live Attenuated Strains of Foot and Mouth Disease Virus

Foot-and-mouth disease (FMD) is one of the most economically important viral diseases that can affect livestock. In the last 70 years, use of an inactivated whole antigen vaccine has contributed to the eradication of disease from many developed nations. However, recent outbreaks in Europe and Eastern Asia demonstrated that infection can spread as wildfire causing economic and social devastation. Therefore, it is essential to develop new control strategies that could confer early protection and rapidly stop disease spread. Live attenuated vaccines (LAV) are one of the best choices to obtain a strong early and long-lasting protection against viral diseases. In proof of concept studies, we previously demonstrated that “synonymous codon deoptimization” could be applied to the P1 capsid coding region of the viral genome to derive attenuated FMDV serotype A12 strains. Here, we demonstrate that a similar approach can be extended to the highly conserved non-structural P2 and P3 coding regions, providing a backbone for multiple serotype FMDV LAV development. Engineered codon deoptimized P2, P3 or P2, and P3 combined regions were included into the A₂₄Cruzeiro infectious clone optimized for vaccine production, resulting in viable progeny that exhibited different degrees of attenuation in cell culture, in mice, and in the natural host (swine). Derived strains were thoroughly characterized in vitro and in vivo. Our work demonstrates that overall, the entire FMDV genome tolerates codon deoptimization, highlighting the potential of using this technology to derive novel improved LAV candidates.

Coronavirus genomes carry the signatures of their habitats

Coronaviruses such as SARS-CoV-2 regularly infect host tissues that express antiviral proteins (AVPs) in abundance. Understanding how they evolve to adapt or evade host immune responses is important in the effort to control the spread of infection. Two AVPs that may shape viral genomes are the zinc finger antiviral protein (ZAP) and the apolipoprotein B mRNA editing enzyme-catalytic polypeptide-like 3 (APOBEC3). The former binds to CpG dinucleotides to facilitate the degradation of viral transcripts while the latter frequently deaminates C into U residues which could generate notable viral sequence variations. We tested the hypothesis that both APOBEC3 and ZAP impose selective pressures that shape the genome of an infecting coronavirus. Our investigation considered a comprehensive number of publicly available genomes for seven coronaviruses (SARS-CoV-2, SARS-CoV, and MERS infecting Homo sapiens, Bovine CoV infecting Bos taurus, MHV infecting Mus musculus, HEV infecting Sus scrofa, and CRCoV infecting Canis lupus familiaris). We show that coronaviruses that regularly infect tissues with abundant AVPs have CpG-deficient and U-rich genomes; whereas those that do not infect tissues with abundant AVPs do not share these sequence hallmarks. Among the coronaviruses surveyed herein, CpG is most deficient in SARS-CoV-2 and a temporal analysis showed a marked increase in C to U mutations over four months of SARS-CoV-2 genome evolution. Furthermore, the preferred motifs in which these C to U mutations occur are the same as those subjected to APOBEC3 editing in HIV-1. These results suggest that both ZAP and APOBEC3 shape the SARS-CoV-2 genome: ZAP imposes a strong CpG avoidance, and APOBEC3 constantly edits C to U. Evolutionary pressures exerted by host immune systems onto viral genomes may motivate novel strategies for SARS-CoV-2 vaccine development.

Use of a small DNA virus model to investigate mechanisms of CpG dinucleotide-induced attenuation of virus replication

Suppression of the CpG dinucleotide is widespread in RNA viruses infecting vertebrates and plants, and in the genomes of retroviruses and small mammalian DNA viruses. The functional basis for CpG suppression in the latter was investigated through the construction of mutants of the parvovirus, minute virus of mice (MVM) with increased CpG or TpA dinucleotides in the VP gene. CpG-high mutants displayed extraordinary attenuation in A9 cells compared to wild-type MVM (>six logs), while TpA elevation showed no replication effect. Attenuation was independent of Toll-like receptor 9 and STING-mediated DNA recognition pathways and unrelated to effects on translation efficiency. While translation from codon-optimized VP RNA was enhanced in a cell-free assay, MVM containing this sequence was highly attenuated. Further mutational analysis indicated that this arose through its increased numbers of CpG dinucleotides (7→70) and separately from its increased G+C content (42.3→57.4 %), which independently attenuated replication. CpG-high viruses showed impaired NS mRNA expression by qPCR and reduced NS and particularly VP protein expression detected by immunofluorescence and replication in A549 cells, effects reversed in zinc antiviral protein (ZAP) knockout cells, even though nuclear relocalization of VP remained defective. The demonstrated functional basis for CpG suppression in MVM and potentially other small DNA viruses and the observed intolerance of CpGs in coding sequences, even after codon optimization, has implications for the use of small DNA virus vectors in gene therapy and immunization.

ΦX174 Attenuation by Whole-Genome Codon Deoptimization

Natural selection acting on synonymous mutations in protein-coding genes influences genome composition and evolution. In viruses, introducing synonymous mutations in genes encoding structural proteins can drastically reduce viral growth, providing a means to generate potent, live-attenuated vaccine candidates. However, an improved understanding of what compositional features are under selection and how combinations of synonymous mutations affect viral growth is needed to predictably attenuate viruses and make them resistant to reversion. We systematically recoded all nonoverlapping genes of the bacteriophage ΦX174 with codons rarely used in its Escherichia coli host. The fitness of recombinant viruses decreases as additional deoptimizing mutations are made to the genome, although not always linearly, and not consistently across genes. Combining deoptimizing mutations may reduce viral fitness more or less than expected from the effect size of the constituent mutations and we point out difficulties in untangling correlated compositional features. We test our model by optimizing the same genes and find that the relationship between codon usage and fitness does not hold for optimization, suggesting that wild-type ΦX174 is at a fitness optimum. This work highlights the need to better understand how selection acts on patterns of synonymous codon usage across the genome and provides a convenient system to investigate the genetic determinants of virulence.

Extreme Genomic CpG Deficiency in SARS-CoV-2 and Evasion of Host Antiviral Defense

Wild mammalian species, including bats, constitute the natural reservoir of betacoronavirus (including SARS, MERS, and the deadly SARS-CoV-2). Different hosts or host tissues provide different cellular environments, especially different antiviral and RNA modification activities that can alter RNA modification signatures observed in the viral RNA genome. The zinc finger antiviral protein (ZAP) binds specifically to CpG dinucleotides and recruits other proteins to degrade a variety of viral RNA genomes. Many mammalian RNA viruses have evolved CpG deficiency. Increasing CpG dinucleotides in these low-CpG viral genomes in the presence of ZAP consistently leads to decreased viral replication and virulence. Because ZAP exhibits tissue-specific expression, viruses infecting different tissues are expected to have different CpG signatures, suggesting a means to identify viral tissue-switching events. The author shows that SARS-CoV-2 has the most extreme CpG deficiency in all known betacoronavirus genomes. This suggests that SARS-CoV-2 may have evolved in a new host (or new host tissue) with high ZAP expression. A survey of CpG deficiency in viral genomes identified a virulent canine coronavirus (alphacoronavirus) as possessing the most extreme CpG deficiency, comparable with that observed in SARS-CoV-2. This suggests that the canine tissue infected by the canine coronavirus may provide a cellular environment strongly selecting against CpG. Thus, viral surveys focused on decreasing CpG in viral RNA genomes may provide important clues about the selective environments and viral defenses in the original hosts.

Intra-genome variability in the dinucleotide composition of SARS-CoV-2

CpG dinucleotides are under-represented in the genomes of single-stranded RNA viruses, and SARS-CoV-2 is no exception to this. Artificial modification of CpG frequency is a valid approach for live attenuated vaccine development; if this is to be applied to SARS-CoV-2, we must first understand the role CpG motifs play in regulating SARS-CoV-2 replication. Accordingly, the CpG composition of the SARS-CoV-2 genome was characterised. CpG suppression among coronaviruses does not differ between virus genera but does vary with host species and primary replication site (a proxy for tissue tropism), supporting the hypothesis that viral CpG content may influence cross-species transmission. Although SARS-CoV-2 exhibits overall strong CpG suppression, this varies considerably across the genome, and the Envelope (E) open reading frame (ORF) and ORF10 demonstrate an absence of CpG suppression. Across the Coronaviridae, E genes display remarkably high variation in CpG composition, with those of SARS and SARS-CoV-2 having much higher CpG content than other coronaviruses isolated from humans. This is an ancestrally derived trait reflecting their bat origins. Conservation of CpG motifs in these regions suggests that they have a functionality which over-rides the need to suppress CpG; an observation relevant to future strategies towards a rationally attenuated SARS-CoV-2 vaccine.

Insights into The Codon Usage Bias of 13 Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) Isolates from Different Geo-locations

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is the causative agent of Coronavirus disease 2019 (COVID-19) which is an infectious disease that spread throughout the world and was declared as a pandemic by the World Health Organization (WHO). In this study, we performed a genome-wide analysis on the codon usage bias (CUB) of 13 SARS-CoV-2 isolates from different geo-locations (countries) in an attempt to characterize it, unravel the main force shaping its pattern, and understand its adaptation to Homo sapiens . Overall results revealed that, SARS-CoV-2 codon usage is slightly biased similarly to other RNA viruses. Nucleotide and dinucleotide compositions displayed a bias toward A/U content in all codon positions and CpU-ended codons preference, respectively. Eight common putative preferred codons were identified, and all of them were A/U-ended (U-ended: 7, A-ended: 1). In addition, natural selection was found to be the main force structuring the codon usage pattern of SARS-CoV-2. However, mutation pressure and other factors such as compositional constraints and hydrophobicity had an undeniable contribution. Two adaptation indices were utilized and indicated that SARS-CoV-2 is moderately adapted to Homo sapiens compared to other human viruses. The outcome of this study may help in understanding the underlying factors involved in the evolution of SARS-CoV-2 and may aid in vaccine design strategies.

The effects of synonymous codon usages on genotypic formation of open reading frames in hepatitis E virus

Hepatitis E virus (HEV) infection has emerged as an important public health issue. As a zoonotic RNA virus, new strains are continuously discovered from human or various animal species. However, the capability of cross-species infection varies largely among different strains. Because the classical nucleotide-based genotyping system provides little functional insight, this study aimed to comprehensively investigate codon usage of the HEV coding regions for better understanding the evolutional orientation, virus-host interaction and cross-species transmission. We observed significant differences of the four nucleotide usages in the three open reading frames, indicating that the evolutional tendency of HEV caused by mutation pressure is modified by the evolutional dynamic related to positive selection. Furthermore, significant differences of nucleotide usages were found among HEV isolated from different host species, suggesting an important role of natural selection related to the host. Analysis of effective number of codons revealed distinct degrees of biased codon usage caused by mutation pressure or the host. Finally, we have mapped the similarity levels of the overall codon usage between the virus and the host to assess the potential of cross-species infection. Thus, this study has provided a novel aspect for better understanding the HEV genetic orientation and the zoonotic nature.

Synonymous genome recoding: a tool to explore microbial biology and new therapeutic strategies

Synthetic genome recoding is a new means of generating designed organisms with altered phenotypes. Synonymous mutations introduced into the protein coding region tolerate modifications in DNA or mRNA without modifying the encoded proteins. Synonymous genome-wide recoding has allowed the synthetic generation of different small-genome viruses with modified phenotypes and biological properties. Recently, a decreased cost of chemically synthesizing DNA and improved methods for assembling DNA fragments (e.g. lambda red recombination and CRISPR-based editing) have enabled the construction of an Escherichia coli variant with a 4-Mb synthetic synonymously recoded genome with a reduced number of sense codons (n = 59) encoding the 20 canonical amino acids. Synonymous genome recoding is increasing our knowledge of microbial interactions with innate immune responses, identifying functional genome structures, and strategically ameliorating cis-inhibitory signaling sequences related to splicing, replication (in eukaryotes), and complex microbe functions, unraveling the relevance of codon usage for the temporal regulation of gene expression and the microbe mutant spectrum and adaptability. New biotechnological and therapeutic applications of this methodology can easily be envisaged. In this review, we discuss how synonymous genome recoding may impact our knowledge of microbial biology and the development of new and better therapeutic methodologies.

Cellular Factors Targeting HIV-1 Transcription and Viral RNA Transcripts

Restriction factors are structurally and functionally diverse cellular proteins that constitute a first line of defense against viral pathogens. Exceptions exist, but typically these proteins are upregulated by interferons (IFNs), target viral components, and are rapidly evolving due to the continuous virus–host arms race. Restriction factors may target HIV replication at essentially each step of the retroviral replication cycle, and the suppression of viral transcription and the degradation of viral RNA transcripts are emerging as major innate immune defense mechanisms. Recent data show that some antiviral factors, such as the tripartite motif-containing protein 22 (TRIM22) and the γ-IFN-inducible protein 16 (IFI16), do not target HIV-1 itself but limit the availability of the cellular transcription factor specificity protein 1 (Sp1), which is critical for effective viral gene expression. In addition, several RNA-interacting cellular factors including RNAse L, the NEDD4-binding protein 1 (N4BP1), and the zinc finger antiviral protein (ZAP) have been identified as important immune effectors against HIV-1 that may be involved in the maintenance of the latent viral reservoirs, representing the major obstacle against viral elimination and cure. Here, we review recent findings on specific cellular antiviral factors targeting HIV-1 transcription or viral RNA transcripts and discuss their potential role in viral latency.

CpG-creating mutations are costly in many human viruses

Mutations can occur throughout the virus genome and may be beneficial, neutral or deleterious. We are interested in mutations that yield a C next to a G, producing CpG sites. CpG sites are rare in eukaryotic and viral genomes. For the eukaryotes, it is thought that CpG sites are rare because they are prone to mutation when methylated. In viruses, we know less about why CpG sites are rare. A previous study in HIV suggested that CpG-creating transition mutations are more costly than similar non-CpG-creating mutations. To determine if this is the case in other viruses, we analyzed the allele frequencies of CpG-creating and non-CpG-creating mutations across various strains, subtypes, and genes of viruses using existing data obtained from Genbank, HIV Databases, and Virus Pathogen Resource. Our results suggest that CpG sites are indeed costly for most viruses. By understanding the cost of CpG sites, we can obtain further insights into the evolution and adaptation of viruses.

Development of a new oral poliovirus vaccine for the eradication end game using codon deoptimization

Enormous progress has been made in global efforts to eradicate poliovirus, using live-attenuated Sabin oral poliovirus vaccine (OPV). However, as the incidence of disease due to wild poliovirus has declined, vaccine-derived poliovirus (VDPV) has emerged in areas of low-vaccine coverage. Coordinated global cessation of routine, type 2 Sabin OPV (OPV2) use has not resulted in fewer VDPV outbreaks, and continued OPV use in outbreak-response campaigns has seeded new emergences in low-coverage areas. The limitations of existing vaccines and current eradication challenges warranted development of more genetically stable OPV strains, most urgently for OPV2. Here, we report using codon deoptimization to further attenuate Sabin OPV2 by changing preferred codons across the capsid to non-preferred, synonymous codons. Additional modifications to the 5′ untranslated region stabilized known virulence determinants. Testing of this codon-deoptimized new OPV2 candidate (nOPV2-CD) in cell and animal models demonstrated that nOPV2-CD is highly attenuated, grows sufficiently for vaccine manufacture, is antigenically indistinguishable from Sabin OPV2, induces neutralizing antibodies as effectively as Sabin OPV2, and unlike Sabin OPV2 is genetically stable and maintains an attenuation phenotype. In-human clinical trials of nOPV2-CD are ongoing, with potential for nOPV strains to serve as critical vaccine tools for achieving and maintaining polio eradication.

CpG Dinucleotides Inhibit HIV-1 Replication through Zinc Finger Antiviral Protein (ZAP)-Dependent and -Independent Mechanisms

CpG dinucleotides are suppressed in the genomes of many vertebrate RNA viruses, including HIV-1. The cellular antiviral protein ZAP (zinc finger antiviral protein) binds CpGs and inhibits HIV-1 replication when CpGs are introduced into the viral genome. However, it is not known if ZAP-mediated restriction is the only mechanism driving CpG suppression. To determine how CpG dinucleotides affect HIV-1 replication, we increased their abundance in multiple regions of the viral genome and analyzed the effect on RNA expression, protein abundance, and infectious-virus production. We found that the antiviral effect of CpGs was not correlated with their abundance. Interestingly, CpGs inserted into some regions of the genome sensitize the virus to ZAP antiviral activity more efficiently than insertions into other regions, and this sensitivity can be modulated by interferon treatment or ZAP overexpression. Furthermore, the sensitivity of the virus to endogenous ZAP was correlated with its sensitivity to the ZAP cofactor KHNYN. Finally, we show that CpGs in some contexts can also inhibit HIV-1 replication by ZAP-independent mechanisms, and one of these is the activation of a cryptic splice site at the expense of a canonical splice site. Overall, we show that the location and sequence context of the CpG in the viral genome determines its antiviral activity.IMPORTANCE Some RNA virus genomes are suppressed in the nucleotide combination of a cytosine followed by a guanosine (CpG), indicating that they are detrimental to the virus. The antiviral protein ZAP binds viral RNA containing CpGs and prevents the virus from multiplying. However, it remains unknown how the number and position of CpGs in viral genomes affect restriction by ZAP and whether CpGs have other antiviral mechanisms. Importantly, manipulating the CpG content in viral genomes could help create new vaccines. HIV-1 shows marked CpG suppression, and by introducing CpGs into its genome, we show that ZAP efficiently targets a specific region of the viral genome, that the number of CpGs does not predict the magnitude of antiviral activity, and that CpGs can inhibit HIV-1 gene expression through a ZAP-independent mechanism. Overall, the position of CpGs in the HIV-1 genome determines the magnitude and mechanism through which they inhibit the virus.

CpG-Recoding in Zika Virus Genome Causes Host-Age-Dependent Attenuation of Infection With Protection Against Lethal Heterologous Challenge in Mice

Experimental increase of CpG dinucleotides in an RNA virus genome impairs infection providing a promising approach for vaccine development. While CpG recoding is an emerging and promising vaccine approach, little is known about infection phenotypes caused by recoded viruses in vivo. For example, infection phenotypes, immunogenicity, and protective efficacy induced by CpG-recoded viruses in different age groups were not studied yet. This is important, because attenuation of infection phenotypes caused by recoded viruses may depend on the population-based expression of cellular components targeting viral CpG dinucleotides. In the present study, we generated several Zika virus (ZIKV) variants with the increasing CpG content and compared infection in neonatal and adult mice. Increasing the CpG content caused host-age-dependent attenuation of infection with considerable attenuation in neonates and high attenuation in adults. Expression of the zinc-finger antiviral protein (ZAP)—the host protein targeting viral CpG dinucleotides—was also age-dependent. Similar to the wild-type virus, ZIKV variants with the increased CpG content evoked robust cellular and humoral immune responses and protection against lethal challenge. Collectively, the host age should be accounted for in future studies on mechanisms targeting viral CpG dinucleotides, development of safe dinucleotide recoding strategies, and applications of CpG-recoded vaccines.

CpG Frequency in the 5' Third of the env Gene Determines Sensitivity of Primary HIV-1 Strains to the Zinc-Finger Antiviral Protein

CpG dinucleotide suppression has been reported to allow HIV-1 to evade inhibition by the zinc-finger antiviral protein (ZAP). Here, we show that primate lentiviruses display marked differences in CpG frequencies across their genome, ranging from 0.44% in simian immunodeficiency virus SIVwrc from Western red colobus to 2.3% in SIVmon infecting mona monkeys. Moreover, functional analyses of a large panel of human and simian immunodeficiency viruses revealed that the magnitude of CpG suppression does not correlate with their susceptibility to ZAP. However, we found that the number of CpG dinucleotides within a region of ∼700 bases at the 5' end of the env gene determines ZAP sensitivity of primary HIV-1 strains but not of HIV-2. Increased numbers of CpGs in this region were associated with reduced env mRNA expression and viral protein production. ZAP sensitivity profiles of chimeric simian-human immunodeficiency viruses (SHIVs) expressing different HIV-1 env genes were highly similar to those of the corresponding HIV-1 strains. The frequency of CpGs in the identified env region correlated with differences in clinical progression rates. Thus, the CpG frequency in a specific part of env, rather than the overall genomic CpG content, governs the susceptibility of HIV-1 to ZAP and might affect viral pathogenicity in vivoIMPORTANCE Evasion of the zinc-finger antiviral protein (ZAP) may drive CpG dinucleotide suppression in HIV-1 and many other viral pathogens but the viral determinants of ZAP sensitivity are poorly defined. Here, we examined CpG suppression and ZAP sensitivity in a large number of primate lentiviruses and demonstrate that their genomic frequency of CpGs varies substantially and does not correlate with ZAP sensitivity. We further show that the number of CpG residues in a defined region at the 5' end of the env gene together with structural features plays a key role in HIV-1 susceptibility to ZAP and correlates with differences in clinical progression rates in HIV-1-infected individuals. Our identification of a specific part of env as a major determinant of HIV-1 susceptibility to ZAP restriction provides a basis for future studies of the underlying inhibitory mechanisms and their potential relevance in the pathogenesis of AIDS.