Modulation of poliovirus replicative fitness in HeLa cells by deoptimization of synonymous codon usage in the capsid region

Comprehensive Analysis of Codon Usage on Porcine Astrovirus

Porcine astrovirus (PAstV), associated with mild diarrhea and neurological disease, is transmitted in pig farms worldwide. The purpose of this study is to elucidate the main factors affecting codon usage to PAstVs. Phylogenetic analysis showed that the subtype PAstV-5 sat at the bottom of phylogenetic tree, followed by PAstV-3, PAstV-1, PAstV-2, and PAstV-4, indicating that the five existing subtypes (PAstV1-PAstV5) may be formed by multiple differentiations of PAstV ancestors. A codon usage bias was found in the PAstVs-2,3,4,5 from the analyses of effective number of codons (ENC) and relative synonymous codon usage (RSCU). Nucleotides A/U are more frequently used than nucleotides C/G in the genome CDSs of the PAstVs-3,4,5. Codon usage patterns of PAstV-5 are dominated by mutation pressure and natural selection, while natural selection is the main evolutionary force that affects the codon usage pattern of PAstVs-2,3,4. The analyses of codon adaptation index (CAI), relative codon deoptimization index (RCDI), and similarity index (SiD) showed the codon usage similarities between the PAstV and animals might contribute to the broad host range and the cross-species transmission of astrovirus. Our results provide insight into understanding the PAstV evolution and codon usage patterns.

Viral adaption of staphylococcal phage: A genome-based analysis of the selective preference based on codon usage Bias

Given the high therapeutic value of the staphylococcal phage, the genome co-evolution of the phage and the host has gained great attention. Though the genome-wide AT richness in staphylococcal phages has been well-studied with nucleotide usage bias, here we proved that host factor, lifestyle and taxonomy are also important factors in understanding the phage nucleotide usages bias using information entropy formula. Such correlation is especially prominent when it comes to the synonymous codon usages of staphylococcal phages, despite the overall scattered codon usage pattern represented by principal component analysis. This strong relationship is explained by nucleotide skew which testified that the usage biases of nucleotide at different codon positions are acting on synonymous codons. Therefore, our study reveals a hidden relationship of genome evolution with host limitation and phagic phenotype, providing new insight into phage genome evolution at genetic level.

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.

Analysis of codon usage bias in potato virus Y non-recombinant strains

Potato virus Y (PVY) is a member of the genus Potyvirus, family Potyviridae, is considered one of the most devastating pest affecting economically important crops, such as potato, tobacco, tomato and pepper, representing a serious threat due to high incidence and worldwide distribution. Its economic significance as well as it biological and molecular complexities have aroused great attention, thus several studies have explore it genetic characteristics. However, little is known about PVY codon usage. To shed light on the relation of codon usage among viruses and their hosts is extremely important to understand virus survival, fitness and evolution. In this study, we performed a comprehensive analysis of codon usage and composition of PVY non-recombinant strains (PVY^N-NA, PVY^Eu-N, PVY^O, PVY^O5, PVY^C) based on 130 complete open reading frame sequences extracted from public databases. Furthermore, similarities between the synonymous codon usage of PVY and its main hosts were investigated. The results obtained in the current study suggest that the overall codon usage among PVY genotypes is similar and slightly biased. PVY codon usage is strongly influenced by mutational bias, but also by G + C compositional constraint and dinucleotide composition. Furthermore, similarities among codon usage preferences between PVY strains and analyzed hosts were observed.

Long-term context-dependent genetic adaptation of the viral genetic cloud

MOTIVATION

RNA viruses generate a cloud of genetic variants within each host. This cloud contains high-frequency genotypes, and many rare variants. The dynamics of these variants is crucial to understand viral evolution and their effect on their host.

RESULTS

We use an experimental evolution system to show that the genetic cloud surrounding the Coxsackie virus master sequence slowly, but steadily, evolves over hundreds of generations. This movement is determined by strong context-dependent mutations, where the frequency and type of mutations are affected by neighboring positions, even in silent mutations. This context-dependent mutation pattern serves as a spearhead for the viral population's movement within the adaptive landscape and affects which new dominant variants will emerge. The non-local mutation patterns affect the mutated dinucleotide distribution, and eventually lead to a non-uniform dinucleotide distribution in the main viral sequence. We tested these results on other RNA viruses with similar conclusions.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

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.

The rabbit papillomavirus model: a valuable tool to study viral-host interactions

Cottontail rabbit papillomavirus (CRPV) was the first DNA virus shown to be tumorigenic. The virus has since been renamed and is officially known as Sylvilagus floridanus papillomavirus 1 (SfPV1). Since its inception as a surrogate preclinical model for high-risk human papillomavirus (HPV) infections, the SfPV1/rabbit model has been widely used to study viral-host interactions and has played a pivotal role in the successful development of three prophylactic virus-like particle vaccines. In this review, we will focus on the use of the model to gain a better understanding of viral pathogenesis, gene function and host immune responses to viral infections. We will discuss the application of the model in HPV-associated vaccine testing, in therapeutic vaccine development (using our novel HLA-A2.1 transgenic rabbits) and in the development and validation of novel anti-viral and anti-tumour compounds. Our goal is to demonstrate the role the SfPV1/rabbit model has played, and continues to play, in helping to unravel the intricacies of papillomavirus infections and to develop tools to thwart the disease. This article is part of the theme issue 'Silent cancer agents: multi-disciplinary modelling of human DNA oncoviruses'.

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.

A Lassa Fever Live-Attenuated Vaccine Based on Codon Deoptimization of the Viral Glycoprotein Gene

Lassa virus (LASV) is endemic in Western Africa and is estimated to infect hundreds of thousands of individuals annually. A considerable number of these infections result in Lassa fever (LF), which is associated with significant morbidity and a case-fatality rate as high as 69% among hospitalized confirmed patients. U.S. Food and Drug Administration-approved LF vaccines are not available. Current antiviral treatment is limited to off-label use of a nucleoside analogue, ribavirin, that is only partially effective and associated with significant side effects. We generated and characterized a recombinant LASV expressing a codon-deoptimized (CD) glycoprotein precursor gene (GPC), rLASV-GPC/CD. Comparison of growth kinetics and peak titers showed that rLASV-GPC/CD is slightly attenuated in cell culture compared to wild-type (WT) recombinant LASV (rLASV-WT). However, rLASV-GPC/CD is highly attenuated in strain 13 and Hartley guinea pigs, as reflected by the absence of detectable clinical signs in animals inoculated with rLASV-GPC/CD. Importantly, a single subcutaneous dose of rLASV-GPC/CD provides complete protection against an otherwise lethal exposure to LASV. Our results demonstrate the feasibility of implementing a CD approach for developing a safe and effective LASV live-attenuated vaccine candidate. Moreover, rLASV-GPC/CD might provide investigators with a tool to safely study LASV outside maximum (biosafety level 4) containment, which could accelerate the elucidation of basic aspects of the molecular and cell biology of LASV and the development of novel LASV medical countermeasures.IMPORTANCE Lassa virus (LASV) infects several hundred thousand people in Western Africa, resulting in many lethal Lassa fever (LF) cases. Licensed LF vaccines are not available, and anti-LF therapy is limited to off-label use of the nucleoside analog ribavirin with uncertain efficacy. We describe the generation of a novel live-attenuated LASV vaccine candidate. This vaccine candidate is based on mutating wild-type (WT) LASV in a key region of the viral genome, the glycoprotein precursor (GPC) gene. These mutations do not change the encoded GPC but interfere with its production in host cells. This mutated LASV (rLASV-GPC/CD) behaves like WT LASV (rLASV-WT) in cell culture, but in contrast to rLASV-WT, does not cause disease in inoculated guinea pigs. Guinea pigs immunized with rLASV-GPC/CD were protected against an otherwise lethal exposure to WT LASV. Our results support the testing of this candidate vaccine in nonhuman primate models ofLF.

Emerging viruses and current strategies for vaccine intervention

During the past decade several notable viruses have suddenly emerged from obscurity or anonymity to become serious global health threats, provoking concern regarding their sustained epidemic transmission in immunologically naive human populations. With each new threat comes the call for rapid vaccine development. Indeed, vaccines are considered a critical component of disease prevention for emerging viral infections because, in many cases, other medical options are limited or non‐existent, or that infections result in such a rapid clinical deterioration that the effectiveness of therapeutics is limited. While classic approaches to vaccine development are still amenable to emerging viruses, the application of molecular techniques in virology has profoundly influenced our understanding of virus biology, and vaccination methods based on replicating, attenuated and non‐replicating virus vector approaches have become useful vaccine platforms. Together with a growing understanding of viral disease emergence, a range of vaccine strategies and international commitment to underpin development, vaccine intervention for new and emerging viruses may become a possibility.

The fit of codon usage of human-isolated avian influenza A viruses to human

Avian influenza A viruses (AIVs) classify into 18 hemagglutinin (HA) and 11 neuraminidase (NA) subtypes. Even though H1N1 and H3N2 subtypes usually circulate among humans leading to infection, occasionally, H5, H6, H7, H9, and H10 that circulate in poultry also infect humans, and especially H5N1 and H7N9. Efficient virus replication is a critical factor that influences infection. Codon usage of a virus must coevolve with its host for efficient viral replication, therefore, we conduct a comprehensive analysis of codon usage bias in human-isolated AIVs to test their adaptation to host expression system. The relative synonymous codon usage (RSCU) pattern, and the codon adaptation index (CAI) are calculated for this purpose. We find that all human-isolated AIVs tend to eliminate GC and CpG compositions, which may prevent activation of the host innate immune system. Although codon usage differs between AIV subtypes, our data support the conclusion that natural selection has played a major role and mutation pressure a minor role in shaping codon usage bias in all AIVs. Our efforts discover that codon usage of genes encoding surface proteins of H5N1, and the polymerase genes of H7N9 has better fit to the human expression system. This may associate with their better replication and infection in human.

Interaction of virus populations with their hosts

Viral population numbers are extremely large compared with those of their host species. Population bottlenecks are frequent during the life cycle of viruses and can reduce viral populations transiently to very few individuals. Viruses have to confront several types of constraints that can be divided into basal, cell-dependent, and organism-dependent constraints. Viruses overcome them exploiting a number of molecular mechanisms, with an important contribution of population numbers and genome variation. The adaptive potential of viruses is reflected in modifications of cell tropism and host range, escape to components of the host immune response, and capacity to alternate among different host species, among other phenotypic changes. Despite a fitness cost of most mutations required to overcome a selective constraint, viruses can find evolutionary pathways that ensure their survival in equilibrium with their hosts.

Genome analysis of coxsackievirus B1 isolates during the consecutive alternating administration course of triple antiviral combination in newborn mice

Background

We developed a new approach for the treatment of enterovirus infections, the consecutive alternating administration (CAA) of a combination of enterovirus inhibitors. On the model of coxsackievirus B1 (CVB1) in mice, two phenomena were observed: absence of drug resistance and increased susceptibility to the antivirals. This study aims to clarify the genetic basis of these phenomena.

Methods

Brain samples from CVB1-infected mice subjected to a CAA course with the combination pleconaril/MDL-860/oxoglaucine were used for viral RNA extraction and next generation sequencing. In parallel, samples from monotherapeutic courses of the three substances included in the combination were studied. Whole genome sequence analysis was carried out on all samples.

Results

Samples of pleconaril monotherapy showed mutations in 5′untranslated region, VP3, 2C, 3C and 2A regions of viral RNA, translated in amino acid substitution of the 2A protein. The MDL-860 course induced changes in CVB1 RNA in the VP3 and 2C regions. The oxoglaucine monotherapy samples showed RNA mutation and amino acid substitution in the VP1 region and nucleotide substitution in the 3D region. In the specimens taken from mice subjected to the CAA course with pleconaril/MDL-860/oxoglaucine, the following RNA mutations were established: 5′ untranslated region, 2A, and 2B, and amino acids substitutions in VP3 and 2A, which differ from those mentioned above. These changes could be the reason for the prevention of drug resistance development and also to be considered as the basis for the phenomenon of increased drug susceptibility.

Conclusions

The results reveal that the high anti-enteroviral efficacy of the CAA course is substantiated by the appearance of specific changes in the viral genome.

vhcub: Virus-host codon usage co-adaptation analysis

Viruses show noticeable evolution to adapt and reproduce within their hosts. Theoretically, patterns and factors that affect the codon usage of viruses should reflect evolutionary changes that allow them to optimize their codon usage to their hosts. Some software tools can analyze the codon usage of organisms; however, their performance has room for improvement, as these tools do not focus on examining the codon usage co-adaptation between viruses and their hosts. This paper describes the vhcub R package, which is a crucial tool used to analyze the co-adaptation of codon usage between a virus and its host, with several implementations of indices and plots. The tool is available from: https://cran.r-project.org/web/packages/vhcub/.

The role of ZAP and OAS3/RNAseL pathways in the attenuation of an RNA virus with elevated frequencies of CpG and UpA dinucleotides

Zinc finger antiviral protein (ZAP) is a powerful restriction factor for viruses with elevated CpG dinucleotide frequencies. We report that ZAP similarly mediates antiviral restriction against echovirus 7 (E7) mutants with elevated frequencies of UpA dinucleotides. Attenuation of both CpG- and UpA-high viruses and replicon mutants was reversed in ZAP k/o cell lines, and restored by plasmid-derived reconstitution of expression in k/o cells. In pull-down assays, ZAP bound to viral RNA transcripts with either CpG- and UpA-high sequences inserted in the R2 region. We found no evidence that attenuation of CpG- or UpA-high mutants was mediated through either translation inhibition or accelerated RNA degradation. Reversal of the attenuation of CpG-high, and UpA-high E7 viruses and replicons was also achieved through knockout of RNAseL and oligodenylate synthetase 3 (OAS3), but not OAS1. WT levels of replication of CpG- and UpA-high mutants were observed in OAS3 k/o cells despite abundant expression of ZAP, indicative of synergy or complementation of these hitherto unconnected pathways. The dependence on expression of ZAP, OAS3 and RNAseL for CpG/UpA-mediated attenuation and the variable and often low level expression of these pathway proteins in certain cell types, such as those of the central nervous system, has implications for the use of CpG-elevated mutants as attenuated live vaccines against neurotropic viruses.

Sequential two-step chromatographic purification of infectious poliovirus using ceramic fluoroapatite and ceramic hydroxyapatite columns

Infectious virus purification techniques are important for vaccine development and gene therapy applications. However, the standardized one-step purification technique using ceramic hydroxyapatite (CHAp) has proven unsuitable for poliovirus. Therefore, we designed a sequential two-step chromatographic technique for purification of the infectious Sabin type 2 vaccine strain of poliovirus from the cell culture supernatant. In the first step, we removed protein contaminants from the Sabin type 2 virus fraction by pH gradient elution on a ceramic fluoroapatite column. In the second step, we removed double-stranded DNA derived from host cells by diluting the virus fraction, directly loading it on a CHAp column, and purifying it using a phosphate gradient with 1 M sodium chloride. This process achieved removal rates of more than 99.95% and 99.99% for proteins and double-stranded DNA, respectively, and was highly reproducible and scalable. Furthermore, it is likely that it will be applicable to other virus species.

Chikungunya Virus Vaccine Candidates with Decreased Mutational Robustness Are Attenuated In Vivo and Have Compromised Transmissibility

Chikungunya virus (CHIKV) is a reemerged arbovirus, a member of the Togaviridae family. It circulates through mosquito vectors mainly of the Aedes family and a mammalian host. CHIKV causes chikungunya fever, a mild to severe disease characterized by arthralgia, with some fatal outcomes described. In the past years, several outbreaks mainly caused by enhanced adaptation of the virus to the vector and ineffective control of the contacts between infected mosquito populations and the human host have been reported. Vaccines represent the best solution for the control of insect-borne viruses, including CHIKV, but are often unavailable. We designed live attenuated CHIKVs by applying a rational genomic design based on multiple replacements of synonymous codons. In doing so, the virus mutational robustness (capacity to maintain phenotype despite introduction of mutations to genotype) is decreased, driving the viral population toward deleterious evolutionary trajectories. When the candidate viruses were tested in the insect and mammalian hosts, we observed overall strong attenuation in both and greatly diminished signs of disease. Moreover, we found that the vaccine candidates elicited protective immunity related to the production of neutralizing antibodies after a single dose. During an experimental transmission cycle between mosquitoes and naive mice, vaccine candidates could be transmitted by mosquito bite, leading to asymptomatic infection in mice with compromised dissemination. Using deep-sequencing technology, we observed an increase in detrimental (stop) codons, which confirmed the effectiveness of this genomic design. Because the approach involves hundreds of synonymous modifications to the genome, the reversion risk is significantly reduced, rendering the viruses promising vaccine candidates.IMPORTANCE Chikungunya fever is a debilitating disease that causes severe pain to the joints, which can compromise the patient's lifestyle for several months and even in some grave cases lead to death. The etiological agent is chikungunya virus, an alphavirus transmitted by mosquito bite. Currently, there are no approved vaccines or treatments against the disease. In our research, we developed novel live attenuated vaccine candidates against chikungunya virus by applying an innovative genomic design. When tested in the insect and mammalian host, the vaccine candidates did not cause disease, elicited strong protection against further infection, and had low risk of reversion to pathogenic phenotypes.

Evolutionary and Structural Overview of Human Picornavirus Capsid Antibody Evasion

Picornaviruses constitute one of the most relevant viral groups according to their impact on human and animal health. Etiologic agents of a broad spectrum of illnesses with a clinical presentation that ranges from asymptomatic to fatal disease, they have been the cause of uncountable epidemics throughout history. Picornaviruses are small naked RNA-positive single-stranded viruses that include some of the most important pillars in the development of virology, comprising poliovirus, rhinovirus, and hepatitis A virus. Picornavirus infectious particles use the fecal-oral or respiratory routes as primary modes of transmission. In this regard, successful viral spread relies on the capability of viral capsids to (i) shelter the viral genome, (ii) display molecular determinants for cell receptor recognition, (iii) facilitate efficient genome delivery, and (iv) escape from the immune system. Importantly, picornaviruses display a substantial amount of genetic variability driven by both mutation and recombination. Therefore, the outcome of their replication results in the emergence of a genetically diverse cloud of individuals presenting phenotypic variance. The host humoral response against the capsid protein represents the most active immune pressure and primary weapon to control the infection. Since the preservation of the capsid function is deeply rooted in the virus evolutionary dynamics, here we review the current structural evidence focused on capsid antibody evasion mechanisms from that perspective.

Foot-and-Mouth Disease Virus: Immunobiology, Advances in Vaccines and Vaccination Strategies Addressing Vaccine Failures-An Indian Perspective

A mass vaccination campaign in India seeks to control and eventually eradicate foot-and-mouth disease (FMD). Biosanitary measures along with FMD monitoring are being conducted along with vaccination. The implementation of the FMD control program has drastically reduced the incidence of FMD. However, cases are still reported, even in regions where vaccination is carried out regularly. Control of FMD outbreaks is difficult when the virus remains in circulation in the vaccinated population. Various FMD risk factors have been identified that are responsible for FMD in vaccinated areas. The factors are discussed along with strategies to address these challenges. The current chemically inactivated trivalent vaccine formulation containing strains of serotype O, A, and Asia 1 has limitations including thermolability and induction of only short-term immunity. Advantages and disadvantages of several new-generation alternate vaccine formulations are discussed. It is unfeasible to study every incidence of FMD in vaccinated animals/areas in such a big country as India with its huge livestock population. However, at the same time, it is absolutely necessary to identify the precise reason for vaccination failure. Failure to vaccinate is one reason for the occurrence of FMD in vaccinated areas. FMD epidemiology, emerging and re-emerging virus strains, and serological status over the past 10 years are discussed to understand the impact of vaccination and incidences of vaccination failure in India. Other factors that are important in vaccination failure that we discuss include disrupted herd immunity, health status of animals, FMD carrier status, and FMD prevalence in other species. Recommendations to boost the search of alternate vaccine formulation, strengthen the veterinary infrastructure, bolster the real-time monitoring of FMD, as well as a detailed investigation and documentation of every case of vaccination failure are provided with the goal of refining the control program.

Schlafen 11 Restricts Flavivirus Replication

Schlafen 11 (Slfn11) is an interferon-stimulated gene that controls the synthesis of proteins by regulating tRNA abundance. Likely through this mechanism, Slfn11 has previously been shown to impair human immunodeficiency virus type 1 (HIV-1) infection and the expression of codon-biased open reading frames. Because replication of positive-sense single-stranded RNA [(+)ssRNA] viruses requires the immediate translation of the incoming viral genome, whereas negative-sense single-stranded RNA [(-)ssRNA] viruses carry at infection an RNA replicase that makes multiple translation-competent copies of the incoming viral genome, we reasoned that (+)ssRNA viruses will be more sensitive to the effect of Slfn11 on protein synthesis than (-)ssRNA viruses. To evaluate this hypothesis, we tested the effects of Slfn11 on the replication of a panel of ssRNA viruses in the human glioblastoma cell line A172, which naturally expresses Slfn11. Depletion of Slfn11 significantly increased the replication of (+)ssRNA viruses from the Flavivirus genus, including West Nile virus (WNV), dengue virus (DENV), and Zika virus (ZIKV), but had no significant effect on the replication of the (-)ssRNA viruses vesicular stomatitis virus (VSV) (Rhabdoviridae family) and Rift Valley fever virus (RVFV) (Phenuiviridae family). Quantification of the ratio of genome-containing viral particles to PFU indicated that Slfn11 impairs WNV infectivity. Intriguingly, Slfn11 prevented WNV-induced downregulation of a subset of tRNAs implicated in the translation of 11.8% of the viral polyprotein. Low-abundance tRNAs might promote optimal protein folding and enhance viral infectivity, as previously reported. In summary, this study demonstrates that Slfn11 restricts flavivirus replication by impairing viral infectivity.IMPORTANCE We provide evidence that the cellular protein Schlafen 11 (Slfn11) impairs replication of flaviviruses, including West Nile virus (WNV), dengue virus (DENV), and Zika virus (ZIKV). However, replication of single-stranded negative RNA viruses was not affected. Specifically, Slfn11 decreases the infectivity of WNV potentially by preventing virus-induced modifications of the host tRNA repertoire that could lead to enhanced viral protein folding. Furthermore, we demonstrate that Slfn11 is not the limiting factor of this novel broad antiviral pathway.

Effects of Arbovirus Multi-Host Life Cycles on Dinucleotide and Codon Usage Patterns

Arthropod-borne viruses (arboviruses) of vertebrates including dengue, zika, chikungunya, Rift Valley fever, and blue tongue viruses cause extensive morbidity and mortality in humans, agricultural animals, and wildlife across the globe. As obligate intercellular pathogens, arboviruses must be well adapted to the cellular and molecular environment of both their arthropod (invertebrate) and vertebrate hosts, which are vastly different due to hundreds of millions of years of separate evolution. Here we discuss the comparative pressures on arbovirus RNA genomes as a result of a dual host life cycle, focusing on pressures that do not alter amino acids. We summarize what is currently known about arboviral genetic composition, such as dinucleotide and codon usage, and how cyclical infection of vertebrate and invertebrate hosts results in different genetic profiles compared with single-host viruses. To serve as a comparison, we compile what is known about arthropod tRNA, dinucleotide, and codon usages and compare this with vertebrates. Additionally, we discuss the potential roles of genetic robustness in arboviral evolution and how it may vary from other viruses. Overall, both arthropod and vertebrate hosts influence the resulting genetic composition of arboviruses, but a great deal remains to be investigated.

Enterovirus A71 Containing Codon-Deoptimized VP1 and High-Fidelity Polymerase as Next-Generation Vaccine Candidate

Enterovirus A71 (EV-A71) is a major pathogen that causes hand-foot-and-mouth disease (HFMD), which occasionally results in severe neurological complications. In this study, we developed four EV-A71 (rgEV-A71) strains by reverse genetics procedures as possible vaccine candidates. The four rgEV-A71 viruses contained various codon-deoptimized VP1 capsid proteins (VP1-CD) and showed replication rates and antigenicity similar to that of the wild-type virus, while a fifth virus, rg4643C4VP-CD, was unable to form plaques but was still able to be examined by median tissue culture infectious dose (TCID₅₀) titers, which were similar to those of the others, indicating the effect of CD on plaque formation. However, the genome stability showed that there were some mutations which appeared during just one passage of the VP1-CD viruses. Thus, we further constructed VP1-CD rgEV-A71 containing high-fidelity determinants in 3D polymerase (CD-HF), and the number of mutations in CD-HF rgEV-A71 was shown to have decreased. The CD-HF viruses showed less virulence than the parental strain in a mouse infection model. After 14 days postimmunization, antibody titers had increased in mice infected with CD-HF viruses. The mouse antisera showed similar neutralizing antibody titers against various CD-HF viruses and different genotypes of EV-A71. The study demonstrates the proof of concept that VP1 codon deoptimization combined with high-fidelity 3D polymerase decreased EV-A71 mutations and virulence in mice but retained their antigenicity, indicating it is a good candidate for next-generation EV-A71 vaccine development.IMPORTANCE EV-A71 can cause severe neurological diseases with fatality in infants and young children, but there are still no effective drugs to date. Here, we developed a novel vaccine strategy with the combination of CD and HF substitutions to generate the genetically stable reverse genetics virus. We found that CD combined with HF polymerase decreased the virulence but maintained the antigenicity of the virus. This work demonstrated the simultaneous introduction of CD genome sequences and HF substitutions as a potential new strategy to develop attenuated vaccine seed virus. Our work provides insight into the development of a low-virulence candidate vaccine virus through a series of genetic editing of virus sequences while maintaining its antigenicity and genome stability, which will provide an additional strategy for next-generation vaccine development of EV-A71.

Attenuation of Human Respiratory Viruses by Synonymous Genome Recoding

Using computer algorithms and commercial DNA synthesis, one or more ORFs of a microbial pathogen such as a virus can be recoded and deoptimized by several strategies that may involve the introduction of up to thousands of nucleotide (nt) changes without affecting amino acid (aa) coding. The synonymous recoding strategies that have been applied to RNA viruses include: deoptimization of codon or codon-pair usage, which may reduce protein expression among other effects; increased content of immunomodulatory CpG and UpA RNA, which increase immune responses and thereby restrict viral replication; and substitution of serine and leucine codons with synonymous codons for which single-nt substitutions can yield nonsense codons, thus limiting evolutionary potential. This can reduce pathogen fitness and create potential live-attenuated vaccines that may have improved properties. The combined approach of genome recoding, synthetic biology, and reverse genetics offers several advantages for the generation of attenuated RNA viruses. First, synonymous recoding involves many mutations, which should reduce the rate and magnitude of de-attenuation. Second, increasing the amount of recoding can provide increased attenuation. Third, because there are no changes at the aa level, all of the relevant epitopes should be expressed. Fourth, attenuation frequently does not compromise immunogenicity, suggesting that the recoded viruses have increased immunogenicity per infectious particle. Synonymous deoptimization approaches have been applied to two important human viral pathogens, namely respiratory syncytial virus (RSV) and influenza A virus (IAV). This manuscript will briefly review the use of these different methods of synonymous recoding to generate attenuated RSV and IAV strains. It also will review the characterization of these vaccine candidates in vitro and in animal models, and describe several surprising findings with respect to phenotypic and genetic instability of some of these candidates.

Optimizing the dynamics of protein expression

Heterologously expressed genes require adaptation to the host organism to ensure adequate levels of protein synthesis, which is typically approached by replacing codons by the target organism’s preferred codons. In view of frequently encountered suboptimal outcomes we introduce the codon-specific elongation model (COSEM) as an alternative concept. COSEM simulates ribosome dynamics during mRNA translation and informs about protein synthesis rates per mRNA in an organism- and context-dependent way. Protein synthesis rates from COSEM are integrated with further relevant covariates such as translation accuracy into a protein expression score that we use for codon optimization. The scoring algorithm further enables fine-tuning of protein expression including deoptimization and is implemented in the software OCTOPOS. The protein expression score produces competitive predictions on proteomic data from prokaryotic, eukaryotic, and human expression systems. In addition, we optimized and tested heterologous expression of manA and ova genes in Salmonella enterica serovar Typhimurium. Superiority over standard methodology was demonstrated by a threefold increase in protein yield compared to wildtype and commercially optimized sequences.

Analysis of Nipah Virus Codon Usage and Adaptation to Hosts

A recent outbreak of Nipah virus (NiV) in India has caused 17 deaths in people living in districts of Kerala state. Its zoonotic nature, as well as high rate of human-to-human transmission, has led researchers worldwide to work toward understanding the different aspects of the NiV. We performed a codon usage analysis, based on publicly available nucleotide sequences of NiV and its host adaptation, along with other members of the Henipavirus genus in ten hosts. NiV genome encodes nine open reading frames; and overall, no significant bias in codon usage was observed. Aromaticity of proteins had no impact on codon usage. An analysis of preferred codons used by NiV and the tRNA pool in human cells indicated that NiV prefers codons from a suboptimal anticodon tRNA pool. We observed that codon usage by NiV is mainly constrained by compositional and selection pressures, not by mutational forces. Parameters that define NiV and host relatedness in terms of codon usage were analyzed, with a codon adaptation index (CAI), relative codon deoptimization index (RCDI), and similarity index calculations; which indicated that, of all hosts analyzed, NiV was best adapted to African green monkeys. A comparative analysis based on the relative codon deoptimization index (RCDI) for host adaptation of NiV, Hendra virus (HeV), Cedar virus (CedV), and Hendra like Mojiang virus (MojV) revealed that except for dogs and ferrets, all evaluated hosts were more susceptible to HeV than NiV.

Genetic Evolution and Molecular Selection of the HE Gene of Influenza C Virus

Influenza C virus (ICV) was first identified in humans and swine, but recently also in cattle, indicating a wider host range and potential threat to both the livestock industry and public health than was originally anticipated. The ICV hemagglutinin-esterase (HE) glycoprotein has multiple functions in the viral replication cycle and is the major determinant of antigenicity. Here, we developed a comparative approach integrating genetics, molecular selection analysis, and structural biology to identify the codon usage and adaptive evolution of ICV. We show that ICV can be classified into six lineages, consistent with previous studies. The HE gene has a low codon usage bias, which may facilitate ICV replication by reducing competition during evolution. Natural selection, dinucleotide composition, and mutation pressure shape the codon usage patterns of the ICV HE gene, with natural selection being the most important factor. Codon adaptation index (CAI) and relative codon deoptimization index (RCDI) analysis revealed that the greatest adaption of ICV was to humans, followed by cattle and swine. Additionally, similarity index (SiD) analysis revealed that swine exerted a stronger evolutionary pressure on ICV than humans, which is considered the primary reservoir. Furthermore, a similar tendency was also observed in the M gene. Of note, we found HE residues 176, 194, and 198 to be under positive selection, which may be the result of escape from antibody responses. Our study provides useful information on the genetic evolution of ICV from a new perspective that can help devise prevention and control strategies.

Large-Scale Synonymous Substitutions in Cucumber Mosaic Virus RNA 3 Facilitate Amino Acid Mutations in the Coat Protein

The effect of large-scale synonymous substitutions in a small icosahedral, single-stranded RNA viral genome on virulence, viral titer, and protein evolution were analyzed. The coat protein (CP) gene of the Fny stain of cucumber mosaic virus (CMV) was modified. We created four CP mutants in which all the codons of nine amino acids in the 5' or 3' half of the CP gene were replaced by either the most frequently or the least frequently used synonymous codons in monocot plants. When the dicot host (Nicotiana benthamiana) was inoculated with these four CP mutants, viral RNA titers in uninoculated symptomatic leaves decreased, while all mutants eventually showed mosaic symptoms similar to those for the wild type. The codon adaptation index of these four CP mutants against dicot genes was similar to those of the wild-type CP gene, indicating that the reduction of viral RNA titer was due to deleterious changes of the secondary structure of RNAs 3 and 4. When two 5' mutants were serially passaged in N. benthamiana, viral RNA titers were rapidly restored but competitive fitness remained decreased. Although no nucleic acid changes were observed in the passaged wild-type CMV, one to three amino acid changes were observed in the synonymously mutated CP of each passaged virus, which were involved in recovery of viral RNA titer of 5' mutants. Thus, we demonstrated that deleterious effects of the large-scale synonymous substitutions in the RNA viral genome facilitated the rapid amino acid mutation(s) in the CP to restore the viral RNA titer.IMPORTANCE Recently, it has been known that synonymous substitutions in RNA virus genes affect viral pathogenicity and competitive fitness by alteration of global or local RNA secondary structure of the viral genome. We confirmed that large-scale synonymous substitutions in the CP gene of CMV resulted in decreased viral RNA titer. Importantly, when viral evolution was stimulated by serial-passage inoculation, viral RNA titer was rapidly restored, concurrent with a few amino acid changes in the CP. This novel finding indicates that the deleterious effects of large-scale nucleic acid mutations on viral RNA secondary structure are readily tolerated by structural changes in the CP, demonstrating a novel part of the adaptive evolution of an RNA viral genome. In addition, our experimental system for serial inoculation of large-scale synonymous mutants could uncover a role for new amino acid residues in the viral protein that have not been observed in the wild-type virus strains.

Evolution of Chikungunya virus in mosquito cells

It has been observed that replication of Chikungunya virus (CHIKV) in C6/36 Aedes albopictus cells has little effect on virus evolution. To characterize evolutionary patterns associated with CHIKV replication in mosquito cells, we performed serial passages of the LR2006 strain in Ae. albopictus cells (75 and 30 passages in C6/36 and U4.4 respectively) and Ae. aegypti cells (100 passages in AA-A20 and in AE) and studied genotypic changes accompanying adaptation during this evolutionary process. Quantitative analysis revealed cell specific patterns: low mutation rates in C6/36 cells except when a CHIKV strain pre-adapted to mammalian was used and typical features of adaptation to cell culture conditions with a high number of fixed mutations in AE and AA-A20 cells probably due to the weak permissiveness of these latter cell lines. Altogether, these results suggested that both cell line and viral strain influence rates of viral evolution. In contrast, characteristics and distribution of mutations were qualitatively very similar in all mosquito cells with a high level of parallel evolution including 4 deletion mutations. Serial passage in mammalian cells of viruses pre-adapted to mosquito cells revealed disappearance of almost all shared mutations suggesting that many of these mutational patterns are vector-specific.

The adaptation of codon usage of +ssRNA viruses to their hosts

Viruses depend on their host's cellular structure to survive. Most of them do not have tRNAs, their translation relies on hosts' tRNA pools. Over the course of evolution, viruses needed to optimally exploit cellular processes of their host. Thus, codon usage of a virus should coevolve with its host to efficiently and rapidly replicate. Some viruses can invade a broad spectrum of hosts (BSTVs), while others can invade a narrow spectrum only (NSTVs). Consequently, we test the hypothesis that similarity of codon usage preference and the degree of matching between BSTVs and their hosts will be lower than that of NSTVs, which only need to coevolve with few hosts. We compare the patterns of codon usage in 255 virus genomes to test this hypothesis. Our results show that NSTVs have a higher degree of matching to their hosts' tRNA pools than BSTVs. Further, analysis of the effective number of codons (ENC) infers that codon usage bias of NSTVs is relatively stronger than that of BSTVs. Thus, codon usage of NSTVs tends to better match their host than that of BSTVs. This supports the hypothesis that viruses adapt to the expression system of their host(s).

Zika Virus Attenuation by Codon Pair Deoptimization Induces Sterilizing Immunity in Mouse Models

Zika virus (ZIKV) infection during the large epidemics in the Americas is related to congenital abnormities or fetal demise. To date, there is no vaccine, antiviral drug, or other modality available to prevent or treat Zika virus infection. Here we designed novel live attenuated ZIKV vaccine candidates using a codon pair deoptimization strategy. Three codon pair-deoptimized ZIKVs (Min E, Min NS1, and Min E+NS1) were de novo synthesized and recovered by reverse genetics and contained large amounts of underrepresented codon pairs in the E gene and/or NS1 gene. The amino acid sequence was 100% unchanged. The codon pair-deoptimized variants had decreased replication fitness in Vero cells (Min NS1 ≫ Min E > Min E+NS1), replicated more efficiently in insect cells than in mammalian cells, and demonstrated diminished virulence in a mouse model. In particular, Min E+NS1, the most restrictive variant, induced sterilizing immunity with a robust neutralizing antibody titer, and a single immunization achieved complete protection against lethal challenge and vertical ZIKV transmission during pregnancy. More importantly, due to the numerous synonymous substitutions in the codon pair-deoptimized strains, reversion to wild-type virulence through gradual nucleotide sequence mutations is unlikely. Our results collectively demonstrate that ZIKV can be effectively attenuated by codon pair deoptimization, highlighting the potential of Min E+NS1 as a safe vaccine candidate to prevent ZIKV infections.IMPORTANCE Due to unprecedented epidemics of Zika virus (ZIKV) across the Americas and the unexpected clinical symptoms, including Guillain-Barré syndrome, microcephaly, and other birth defects in humans, there is an urgent need for ZIKV vaccine development. Here we provided the first attenuated versions of ZIKV with two important genes (E and/or NS1) that were subjected to codon pair deoptimization. Compared to parental ZIKV, the codon pair-deoptimized ZIKVs were mammal attenuated and preferred insect to mammalian cells. Min E+NS1, the most restrictive variant, induced sterilizing immunity with a robust neutralizing antibody titer and achieved complete protection against lethal challenge and vertical virus transmission during pregnancy. More importantly, the massive synonymous mutational approach made it impossible for the variant to revert to wild-type virulence. Our results have proven the feasibility of codon pair deoptimization as a strategy to develop live attenuated vaccine candidates against flaviviruses such as ZIKV, Japanese encephalitis virus, and West Nile virus.

HIV-1 Protease Evolvability Is Affected by Synonymous Nucleotide Recoding

One unexplored aspect of HIV-1 genetic architecture is how codon choice influences population diversity and evolvability. Here we compared the levels of development of HIV-1 resistance to protease inhibitors (PIs) between wild-type (WT) virus and a synthetic virus (MAX) carrying a codon-pair-reengineered protease sequence including 38 (13%) synonymous mutations. The WT and MAX viruses showed indistinguishable replication in MT-4 cells or peripheral blood mononuclear cells (PBMCs). Both viruses were subjected to serial passages in MT-4 cells, with selective pressure from the PIs atazanavir (ATV) and darunavir (DRV). After 32 successive passages, both the WT and MAX viruses developed phenotypic resistance to PIs (50% inhibitory concentrations [IC₅₀s] of 14.6 ± 5.3 and 21.2 ± 9 nM, respectively, for ATV and 5.9 ± 1.0 and 9.3 ± 1.9, respectively, for DRV). Ultradeep sequence clonal analysis revealed that both viruses harbored previously described mutations conferring resistance to ATV and DRV. However, the WT and MAX virus proteases showed different resistance variant repertoires, with the G16E and V77I substitutions observed only in the WT and the L33F, S37P, G48L, Q58E/K, and L89I substitutions detected only in the MAX virus. Remarkably, the G48L and L89I substitutions are rarely found in vivo in PI-treated patients. The MAX virus showed significantly higher nucleotide and amino acid diversity of the propagated viruses with and without PIs (P < 0.0001), suggesting a higher selective pressure for change in this recoded virus. Our results indicate that the HIV-1 protease position in sequence space delineates the evolution of its mutant spectrum. Nevertheless, the investigated synonymously recoded variant showed mutational robustness and evolvability similar to those of the WT virus.IMPORTANCE Large-scale synonymous recoding of virus genomes is a new tool for exploring various aspects of virus biology. Synonymous virus genome recoding can be used to investigate how a virus's position in sequence space defines its mutant spectrum, evolutionary trajectory, and pathogenesis. In this study, we evaluated how synonymous recoding of the human immunodeficiency virus type 1 (HIV-1) protease affects the development of protease inhibitor (PI) resistance. HIV-1 protease is a main target of current antiretroviral therapies. Our present results demonstrate that the wild-type (WT) virus and a virus with recoded protease exhibited different patterns of resistance mutations after PI treatment. Nevertheless, the developed PI resistance phenotypes were indistinguishable between the recoded virus and the WT virus, suggesting that the HIV-1 strain with synonymously recoded protease and the WT virus are equally robust and evolvable.

Combinatorial Approaches to Viral Attenuation

Live viral vaccines rely on attenuated viruses that can successfully infect their host but have reduced fitness or virulence. Such attenuated viruses were originally developed through trial and error, typically by adaptation of the wild-type virus to novel conditions. That method was haphazard, with no way of controlling the degree of attenuation or the number of attenuating mutations or preventing evolutionary reversion. Synthetic biology now enables rational design and engineering of viral attenuation, but rational design must be informed by biological principles to achieve stable, quantitative attenuation. This work shows that in a model system for viral attenuation, bacteriophage T7, attenuation can be obtained from rational design principles, and multiple different attenuation approaches can be combined for enhanced overall effect.

Attenuated viruses have numerous applications, in particular in the context of live viral vaccines. However, purposefully designing attenuated viruses remains challenging, in particular if the attenuation is meant to be resistant to rapid evolutionary recovery. Here we develop and analyze a new attenuation method, promoter ablation, using an established viral model, bacteriophage T7. Ablation of promoters of the two most highly expressed T7 proteins (scaffold and capsid) led to major reductions in transcript abundance of the affected genes, with the effect of the double knockout approximately additive of the effects of single knockouts. Fitness reduction was moderate and also approximately additive; fitness recovery on extended adaptation was partial and did not restore the promoters. The fitness effect of promoter knockouts combined with a previously tested codon deoptimization of the capsid gene was less than additive, as anticipated from their competing mechanisms of action. In one design, the engineering created an unintended consequence that led to further attenuation, the effect of which was studied and understood in hindsight. Overall, the mechanisms and effects of genome engineering on attenuation behaved in a predictable manner. Therefore, this work suggests that the rational design of viral attenuation methods is becoming feasible.

IMPORTANCE Live viral vaccines rely on attenuated viruses that can successfully infect their host but have reduced fitness or virulence. Such attenuated viruses were originally developed through trial and error, typically by adaptation of the wild-type virus to novel conditions. That method was haphazard, with no way of controlling the degree of attenuation or the number of attenuating mutations or preventing evolutionary reversion. Synthetic biology now enables rational design and engineering of viral attenuation, but rational design must be informed by biological principles to achieve stable, quantitative attenuation. This work shows that in a model system for viral attenuation, bacteriophage T7, attenuation can be obtained from rational design principles, and multiple different attenuation approaches can be combined for enhanced overall effect.

Exploratory re-encoding of yellow fever virus genome: new insights for the design of live-attenuated viruses

Virus attenuation by genome re-encoding is a pioneering approach for generating effective live-attenuated vaccine candidates. Its core principle is to introduce a large number of synonymous substitutions into the viral genome to produce stable attenuation of the targeted virus. Introduction of large numbers of mutations has also been shown to maintain stability of the attenuated phenotype by lowering the risk of reversion and recombination of re-encoded genomes. Identifying mutations with low fitness cost is pivotal as this increases the number that can be introduced and generates more stable and attenuated viruses. Here, we sought to identify mutations with low deleterious impact on the in vivo replication and virulence of yellow fever virus (YFV). Following comparative bioinformatic analyses of flaviviral genomes, we categorised synonymous transition mutations according to their impact on CpG/UpA composition and secondary RNA structures. We then designed seventeen re-encoded viruses with 100–400 synonymous mutations in the NS2A-to-NS4B coding region of YFV Asibi and Ap7M (hamster-adapted) genomes. Each virus contained a panel of synonymous mutations designed according to the above categorisation criteria. The replication and fitness characteristics of parent and re-encoded viruses were compared in vitro using cell culture competition experiments. In vivo laboratory hamster models were also used to compare relative virulence and immunogenicity characteristics. Most of the re-encoded strains showed no decrease in replicative fitness in vitro. However, they showed reduced virulence and, in some instances, decreased replicative fitness in vivo. Importantly, the most attenuated of the re-encoded strains induced robust, protective immunity in hamsters following challenge with Ap7M, a virulent virus. Overall, the introduction of transitions with no or a marginal increase in the number of CpG/UpA dinucleotides had the mildest impact on YFV replication and virulence in vivo. Thus, this strategy can be incorporated in procedures for the finely tuned creation of substantially re-encoded viral genomes.

PA-X: a key regulator of influenza A virus pathogenicity and host immune responses

PA-X, a fusion protein belonging to influenza A viruses (IAVs), integrating the N-terminal 191 amino acids of PA gene and the ribosomal frame-shifting product that lengthens out to 41 or 61 amino acids. Since its discovery in 2012, multiple functions have been attributed to this small protein, including a process, where wide-spread protein synthesis in infected host cells is shut down (called host shutoff), and viral replication, polymerase activity, viral-induced cell apoptosis, PA nuclear localization, and virulence are modulated. However, many of its proposed functions may be specific to strain, subtype, host, or cell line. In this review, we start by describing the well-defined global host-shutoff ability of PA-X and the potential mechanisms underlying it. We move on to the role played by PA-X in modulating innate and acquired immune responses in the host. We then systematically discuss the role played by PA-X in modulating the virulence of influenza viruses of different subtypes and host origins, and finish with a general overview of the research advances made in identifying the host cell partners that interact with PA-X. To uncover possible clues about the differential effects of PA-X in modulating viral virulence, we focus on systemically evaluating polymorphisms in PA-X from various viral subtypes and hosts, including avian and human H5N1, H5N6, H9N2, and H7N9 viruses. Finally, we conclude with a proposition regarding the possible future research directions for this important protein.

The Evolution of Molecular Compatibility between Bacteriophage ΦX174 and its Host

Viruses rely upon their hosts for biosynthesis of viral RNA, DNA and protein. This dependency frequently engenders strong selection for virus genome compatibility with potential hosts, appropriate gene regulation and expression necessary for a successful infection. While bioinformatic studies have shown strong correlations between codon usage in viral and host genomes, the selective factors by which this compatibility evolves remain a matter of conjecture. Engineered to include codons with a lesser usage and/or tRNA abundance within the host, three different attenuated strains of the bacterial virus ФX174 were created and propagated via serial transfers. Molecular sequence data indicate that biosynthetic compatibility was recovered rapidly. Extensive computational simulations were performed to assess the role of mutational biases as well as selection for translational efficiency in the engineered phage. Using bacteriophage as a model system, we can begin to unravel the evolutionary processes shaping codon compatibility between viruses and their host.

Generation and comparative genomics of synthetic dengue viruses

Background

Synthetic virology is an important multidisciplinary scientific field, with emerging applications in biotechnology and medicine, aiming at developing methods to generate and engineer synthetic viruses. In particular, many of the RNA viruses, including among others the Dengue and Zika, are widespread pathogens of significant importance to human health. The ability to design and synthesize such viruses may contribute to exploring novel approaches for developing vaccines and virus based therapies.

Results

Here we develop a full multidisciplinary pipeline for generation and analysis of synthetic RNA viruses and specifically apply it to Dengue virus serotype 2 (DENV-2). The major steps of the pipeline include comparative genomics of endogenous and synthetic viral strains. Specifically, we show that although the synthetic DENV-2 viruses were found to have lower nucleotide variability, their phenotype, as reflected in the study of the AG129 mouse model morbidity, RNA levels, and neutralization antibodies, is similar or even more pathogenic in comparison to the wildtype master strain. Additionally, the highly variable positions, identified in the analyzed DENV-2 population, were found to overlap with less conserved homologous positions in Zika virus and other Dengue serotypes. These results may suggest that synthetic DENV-2 could enhance virulence if the correct sequence is selected.

Conclusions

The approach reported in this study can be used to generate and analyze synthetic RNA viruses both on genotypic and on phenotypic level. It could be applied for understanding the functionality and the fitness effects of any set of mutations in viral RNA and for editing RNA viruses for various target applications.

Electronic supplementary material

The online version of this article (10.1186/s12859-018-2132-3) contains supplementary material, which is available to authorized users.

Quasispecies and virus

Quasispecies theory has been instrumental in the understanding of RNA virus population dynamics because it considered for the first time mutation as an integral part of the replication process. The key influences of quasispecies theory on experimental virology have been: (1) to disclose the mutant spectrum nature of viral populations and to evaluate its consequences; (2) to unveil collective properties of genome ensembles that can render a mutant spectrum a unit of selection; and (3) to identify new vulnerability points of pathogenic RNA viruses on three fronts: the need to apply multiple selective constraints (in the form of drug combinations) to minimize selection of treatment-escape variants, to translate the error threshold concept into antiviral designs, and to construct attenuated vaccine viruses through alterations of viral polymerase copying fidelity or through displacements of viral genomes towards unfavorable regions of sequence space. These three major influences on the understanding of viral pathogens preceded extensions of quasispecies to non-viral systems such as bacterial and tumor cell collectivities and prions. These developments are summarized here.

Structural divergence creates new functional features in alphavirus genomes

Alphaviruses are mosquito-borne pathogens that cause human diseases ranging from debilitating arthritis to lethal encephalitis. Studies with Sindbis virus (SINV), which causes fever, rash, and arthralgia in humans, and Venezuelan equine encephalitis virus (VEEV), which causes encephalitis, have identified RNA structural elements that play key roles in replication and pathogenesis. However, a complete genomic structural profile has not been established for these viruses. We used the structural probing technique SHAPE-MaP to identify structured elements within the SINV and VEEV genomes. Our SHAPE-directed structural models recapitulate known RNA structures, while also identifying novel structural elements, including a new functional element in the nsP1 region of SINV whose disruption causes a defect in infectivity. Although RNA structural elements are important for multiple aspects of alphavirus biology, we found the majority of RNA structures were not conserved between SINV and VEEV. Our data suggest that alphavirus RNA genomes are highly divergent structurally despite similar genomic architecture and sequence conservation; still, RNA structural elements are critical to the viral life cycle. These findings reframe traditional assumptions about RNA structure and evolution: rather than structures being conserved, alphaviruses frequently evolve new structures that may shape interactions with host immune systems or co-evolve with viral proteins.

Emergency Services of Viral RNAs: Repair and Remodeling

Reproduction of RNA viruses is typically error-prone due to the infidelity of their replicative machinery and the usual lack of proofreading mechanisms. The error rates may be close to those that kill the virus. Consequently, populations of RNA viruses are represented by heterogeneous sets of genomes with various levels of fitness. This is especially consequential when viruses encounter various bottlenecks and new infections are initiated by a single or few deviating genomes. Nevertheless, RNA viruses are able to maintain their identity by conservation of major functional elements. This conservatism stems from genetic robustness or mutational tolerance, which is largely due to the functional degeneracy of many protein and RNA elements as well as to negative selection. Another relevant mechanism is the capacity to restore fitness after genetic damages, also based on replicative infidelity. Conversely, error-prone replication is a major tool that ensures viral evolvability. The potential for changes in debilitated genomes is much higher in small populations, because in the absence of stronger competitors low-fit genomes have a choice of various trajectories to wander along fitness landscapes. Thus, low-fit populations are inherently unstable, and it may be said that to run ahead it is useful to stumble. In this report, focusing on picornaviruses and also considering data from other RNA viruses, we review the biological relevance and mechanisms of various alterations of viral RNA genomes as well as pathways and mechanisms of rehabilitation after loss of fitness. The relationships among mutational robustness, resilience, and evolvability of viral RNA genomes are discussed.

Codon usage variation of Zika virus: The potential roles of NS2B and NS4A in its global pandemic

A comprehensive demonstration of Zika virus (ZIKV) molecular evolution is essential for understanding its adaptation and expansion in its recent pandemics. Despite several studies on mutations and codon usage in ZIKVs, the variations in codon usage patterns across individual genes and their biological implication remains unclear. Here, we performed a gene-by-gene comparison of the codon usage variation in ZIKVs of the African and Asian lineages. We found that besides the evidence of positive selection (Ka/Ks >1) in the Asian lineage of the ZIKV genome, codon usage patterns were gene-specific and codon usage variation of ZIKV genes, was possibly constrained by their individual functional features, such as transmembrane domains, or antigenicity. In particular, the NS2B and NS4A genes showed distinct codon usage patterns, clearly separating them from the clusters of other genes in the correspondence analysis (CA). In the Asian lineage, the NS2B and NS4A genes showed the highest codon usage bias (ENC values: 51.01 ± 0.72 and 48.89 ± 0.99 respectively), and were subjected to the highest translation selection (ENC_obs/ENC_exp ratio: 0.847 ± 0.0297 and 0.828 ± 0.0233 respectively) in comparison to the African lineages of ZIKV. The CpG frequency of the NS2B showed a gradual ascending trend in the Asian ZIKV lineages, while in NS4A it was constrained along with the expansion of the Asian lineage. Furthermore, between the African and Asian lineages, differentiated and specific over-represented codons were more prominent in the NS2B and NS4A. Together, our study implies that ZIKVs are in the process of evolutionary fine tuning their codon as seen in the recent pandemics, and NS2B and NS4A could have played a potential role in the molecular evolution of the Asian lineage and their establishment.

Diagnostic Assay Development for Poliovirus Eradication

With poliovirus eradication nearing, few pockets of active wild poliovirus (WPV) transmission remain in the world. Intratypic differentiation (ITD) plays a crucial part in laboratory surveillance as the molecular detection method that can identify and distinguish wild and vaccine-like polioviruses isolated from acute flaccid paralysis cases or environmental sources. The need to detect new variants of WPV serotype 1 (WPV1) and the containment of all serotype 2 polioviruses (PV2) in 2015 required changes to the previous version of the method. The ITD version 5.0 is a set of six real-time reverse transcription-PCR (rRT-PCR) assays that serve as accurate diagnostic tools to easily detect and differentiate PV serotypes and genotypes. We describe the creation and properties of quantitation standards, including 16 control RNA transcripts and nine plaque-isolated viruses. All ITD rRT-PCR assays were validated using these standards, and the limits of detection were determined for each assay. We designed and pilot tested two new assays targeting recently circulating WPV1 genotypes and all PV2 viruses. The WPV1 assay had 99.1% specificity and 100% sensitivity, and the PV2 assay had 97.7% specificity and 92% sensitivity. Before proceeding to the next step in the global poliovirus eradication program, we needed to gain a better understanding of the performance of the ITD 5.0 suite of molecular assays and their limits of detection and specificities. The findings and conclusions in this evaluation serve as building blocks for future development work.

Attenuation of Viruses by Large-Scale Recoding of their Genomes: the Selection Is Always Biased

PURPOSE OF REVIEW

This review summarizes the current understanding of virus attenuation by large-scale recoding of viral genomes and discusses what would ultimately be necessary for construction of better and safer live modified virus vaccines.

RECENT FINDINGS

It has been shown that codon and codon pair deoptimization are rapid and robust methods that can be used for the development of attenuated vaccine candidates. The viruses attenuated by large-scale recoding have the added value that they are extremely genetically stable. However, the exact mechanisms that lead to viral attenuation by recoding are yet to be determined.

SUMMARY

While the advantages of large-scale recoding (speed, simplicity, potency, and universal applicability) have been known for more than a decade, this approach has been only inadequately explored and the attention was focused on a limited number of RNA viruses. Attenuation of viruses by large-scale recoding should be explored to combat known and future viral threats.

Eradicating infectious disease using weakly transmissible vaccines

Viral vaccines have had remarkable positive impacts on human health as well as the health of domestic animal populations. Despite impressive vaccine successes, however, many infectious diseases cannot yet be efficiently controlled or eradicated through vaccination, often because it is impossible to vaccinate a sufficient proportion of the population. Recent advances in molecular biology suggest that the centuries-old method of individual-based vaccine delivery may be on the cusp of a major revolution. Specifically, genetic engineering brings to life the possibility of a live, transmissible vaccine. Unfortunately, releasing a highly transmissible vaccine poses substantial evolutionary risks, including reversion to high virulence as has been documented for the oral polio vaccine. An alternative, and far safer approach, is to rely on genetically engineered and weakly transmissible vaccines that have reduced scope for evolutionary reversion. Here, we use mathematical models to evaluate the potential efficacy of such weakly transmissible vaccines. Our results demonstrate that vaccines with even a modest ability to transmit can significantly lower the incidence of infectious disease and facilitate eradication efforts. Consequently, weakly transmissible vaccines could provide an important tool for controlling infectious disease in wild and domestic animal populations and for reducing the risks of emerging infectious disease in humans.

Selecting among three basic fitness landscape models: Additive, multiplicative and stickbreaking

Fitness landscapes map genotypes to organismal fitness. Their topographies depend on how mutational effects interact - epistasis - andare important for understanding evolutionary processes such as speciation, the rate of adaptation, the advantage of recombination, and the predictability versus stochasticity of evolution. The growing amount of data has made it possible to better test landscape models empirically. We argue that this endeavor will benefit from the development and use of meaningful basic models against which to compare more complex models. Here we develop statistical and computational methods for fitting fitness data from mutation combinatorial networks to three simple models: additive, multiplicative and stickbreaking. We employ a Bayesian framework for doing model selection. Using simulations, we demonstrate that our methods work and we explore their statistical performance: bias, error, and the power to discriminate among models. We then illustrate our approach and its flexibility by analyzing several previously published datasets. An R-package that implements our methods is available in the CRAN repository under the name Stickbreaker.

Significant differences in terms of codon usage bias between bacteriophage early and late genes: a comparative genomics analysis

Background

Viruses undergo extensive evolutionary selection for efficient replication which effects, among others, their codon distribution. In the current study, we aimed at understanding the way evolution shapes the codon distribution in early vs. late viral genes in terms of their expression during different stages in the viral replication cycle. To this end we analyzed 14 bacteriophages and 11 human viruses with available information about the expression phases of their genes.

Results

We demonstrated evidence of selection for distinct composition of synonymous codons in early and late viral genes in 50% of the analyzed bacteriophages. Among others, this phenomenon may be related to the time specific adaptation of the viral genes to the translation efficiency factors involved at different bacteriophage developmental stages. Specifically, we showed that the differences in codon composition in different temporal gene groups cannot be explained only by phylogenetic proximities between the analyzed bacteriophages, and can be partially explained by differences in the adaptation to the host tRNA pool, nucleotide bias, GC content and more.

In contrast, no difference in temporal regulation of synonymous codon usage was observed in human viruses, possibly because of a stronger selection pressure due to a larger effective population size in bacteriophages and their bacterial hosts.

Conclusions

The codon distribution in large fractions of bacteriophage genomes tend to be different in early and late genes. This phenomenon seems to be related to various aspects of the viral life cycle, and to various intracellular processes. We believe that the reported results should contribute towards better understanding of viral evolution and may promote the development of relevant procedures in synthetic virology.

Electronic supplementary material

The online version of this article (10.1186/s12864-017-4248-7) contains supplementary material, which is available to authorized users.

Transmissible Viral Vaccines

Genetic engineering now enables the design of live viral vaccines that are potentially transmissible. Some designs merely modify a single viral genome to improve on the age-old method of attenuation whereas other designs create chimeras of viral genomes. Transmission has the benefit of increasing herd immunity above that achieved by direct vaccination alone but also increases the opportunity for vaccine evolution, which typically undermines vaccine utility. Different designs have different epidemiological consequences but also experience different evolution. Approaches that integrate vaccine engineering with an understanding of evolution and epidemiology will reap the greatest benefit from vaccine transmission.

Arresting Evolution

Evolution in the form of selective breeding has long been harnessed as a useful tool by humans. However, rapid evolution can also be a danger to our health and a stumbling block for biotechnology. Unwanted evolution can underlie the emergence of drug and pesticide resistance, cancer, and weeds. It makes live vaccines and engineered cells inherently unreliable and unpredictable, and therefore potentially unsafe. Yet, there are strategies that have been and can possibly be used to stop or slow many types of evolution. We review and classify existing population genetics-inspired methods for arresting evolution. Then, we discuss how genome editing techniques enable a radically new set of approaches to limit evolution.

MicroRNA-134 regulates poliovirus replication by IRES targeting

Global poliovirus eradication efforts include high vaccination coverage with live oral polio vaccine (OPV), surveillance for acute flaccid paralysis, and OPV “mop-up” campaigns. An important objective involves host-directed strategies to reduce PV replication to diminish viral shedding in OPV recipients. In this study, we show that microRNA-134-5p (miR-134) can regulate Sabin-1 replication but not Sabin-2 or Sabin-3 via direct interaction with the PV 5′UTR. Hypochromicity data showed miR-134 binding to Sabin-1 and 3 but not Sabin-2 IRES. Transfection of a miR-134 mimic repressed translation of Sabin-1 5′UTR driven luciferase validating the mechanism of miR-134-mediated repression of Sabin-1. Further, site directed mutagenesis of the miR-134 binding site in Sabin-1 IRES relieved miR-134-mediated repression indicating that these regulatory molecules have an important role in regulating the host gene response to PV. Binding of miR-134 to Sabin-1 IRES caused degradation of the IRES transcript in a miR-134 and sequence specific manner. The miR-134 binding site was found to be highly conserved in wild type PV-1 as well as EV71 strains indicating that miR-134 may regulate function of these IRES sequences in circulation.

CpG and UpA dinucleotides in both coding and non-coding regions of echovirus 7 inhibit replication initiation post-entry

Most vertebrate and plant RNA and small DNA viruses suppress genomic CpG and UpA dinucleotide frequencies, apparently mimicking host mRNA composition. Artificially increasing CpG/UpA dinucleotides attenuates viruses through an entirely unknown mechanism. Using the echovirus 7 (E7) model in several cell types, we show that the restriction in E7 replication in mutants with increased CpG/UpA dinucleotides occurred immediately after viral entry, with incoming virions failing to form replication complexes. Sequences of CpG/UpA-high virus stocks showed no evidence of increased mutational errors that would render them replication defective, these viral RNAs were not differentially sequestered in cytoplasmic stress granules nor did they induce a systemic antiviral state. Importantly, restriction was not mediated through effects on translation efficiency since replicons with high CpG/UpA sequences inserted into a non-coding region were similarly replication defective. Host-cells thus possess intrinsic defence pathways that prevent replication of viruses with increased CpG/UpA frequencies independently of codon usage.

Reduced Protein Expression in a Virus Attenuated by Codon Deoptimization

A general means of viral attenuation involves the extensive recoding of synonymous codons in the viral genome. The mechanistic underpinnings of this approach remain unclear, however. Using quantitative proteomics and RNA sequencing, we explore the molecular basis of attenuation in a strain of bacteriophage T7 whose major capsid gene was engineered to carry 182 suboptimal codons. We do not detect transcriptional effects from recoding. Proteomic observations reveal that translation is halved for the recoded major capsid gene, and a more modest reduction applies to several coexpressed downstream genes. We observe no changes in protein abundances of other coexpressed genes that are encoded upstream. Viral burst size, like capsid protein abundance, is also decreased by half. Together, these observations suggest that, in this virus, reduced translation of an essential polycistronic transcript and diminished virion assembly form the molecular basis of attenuation.