Phylodynamics and Codon Usage Pattern Analysis of Broad Bean Wilt Virus 2

Broad bean wilt virus 2 (BBWV-2), which belongs to the genus Fabavirus of the family Secoviridae, is an important pathogen that causes damage to broad bean, pepper, yam, spinach and other economically important ornamental and horticultural crops worldwide. Previously, only limited reports have shown the genetic variation of BBWV2. Meanwhile, the detailed evolutionary changes, synonymous codon usage bias and host adaptation of this virus are largely unclear. Here, we performed comprehensive analyses of the phylodynamics, reassortment, composition bias and codon usage pattern of BBWV2 using forty-two complete genome sequences of BBWV-2 isolates together with two other full-length RNA1 sequences and six full-length RNA2 sequences. Both recombination and reassortment had a significant influence on the genomic evolution of BBWV2. Through phylogenetic analysis we detected three and four lineages based on the ORF1 and ORF2 nonrecombinant sequences, respectively. The evolutionary rates of the two BBWV2 ORF coding sequences were 8.895 × 10⁻⁴ and 4.560 × 10⁻⁴ subs/site/year, respectively. We found a relatively conserved and stable genomic composition with a lower codon usage choice in the two BBWV2 protein coding sequences. ENC-plot and neutrality plot analyses showed that natural selection is the key factor shaping the codon usage pattern of BBWV2. Strong correlations between BBWV2 and broad bean and pepper were observed from similarity index (SiD), codon adaptation index (CAI) and relative codon deoptimization index (RCDI) analyses. Our study is the first to evaluate the phylodynamics, codon usage patterns and adaptive evolution of a fabavirus, and our results may be useful for the understanding of the origin of this virus.     

Evolution of codon usage bias in Drosophila

We first review what is known about patterns of codon usage bias in Drosophila and make the following points: (i) Drosophila genes are as biased or more biased than those in microorganisms. (ii) The level of bias of genes and even the particular pattern of codon bias can remain phylogenetically invariant for very long periods of evolution. (iii) However, some genes, even very tightly linked genes, can change very greatly in codon bias across species. (iv) Generally G and especially C are favored at synonymous sites in biased genes. (v) With the exception of aspartic acid, all amino acids contribute significantly and about equally to the codon usage bias of a gene. (vi) While most individual amino acids that can use G or C at synonymous sites display a preference for C, there are exceptions: valine and leucine, which prefer G. (vii) Finally, smaller genes tend to be more biased than longer genes. We then examine possible causes of these patterns and discount mutation bias on three bases: there is little evidence of regional mutation bias in Drosophila, mutation bias is likely toward A+T (the opposite of codon usage bias), and not all amino acids display the preference for the same nucleotide in the wobble position. Two lines of evidence support a selection hypothesis based on tRNA pools: highly biased genes tend to be highly and/or rapidly expressed, and the preferred codons in highly biased genes optimally bind the most abundant isoaccepting tRNAs. Finally, we examine the effect of bias on DNA evolution and confirm that genes with high codon usage bias have lower rates of synonymous substitution between species than do genes with low codon usage bias. Surprisingly, we find that genes with higher codon usage bias display higher levels of intraspecific synonymous polymorphism. This may be due to opposing effects of recombination.     

Synonymous Codon Usage Analysis of Three Narcissus Potyviruses

Narcissus degeneration virus (NDV), narcissus late season yellows virus (NLSYV) and narcissus yellow stripe virus (NYSV), which belong to the genus Potyvirus of the family Potyviridae, cause significant losses in the ornamental value and quality of narcissus. Several previous studies have explored the genetic diversity and evolution rate of narcissus viruses, but the analysis of the synonymous codons of the narcissus viruses is still unclear. Herein, the coat protein (CP) of three viruses is used to analyze the viruses’ phylogeny and codon usage pattern. Phylogenetic analysis showed that NYSV, NDV and NLSYV isolates were divided into five, three and five clusters, respectively, and these clusters seemed to reflect the geographic distribution. The effective number of codon (ENC) values indicated a weak codon usage bias in the CP coding region of the three narcissus viruses. ENC-plot and neutrality analysis showed that the codon usage bias of the three narcissus viruses is all mainly influenced by natural selection compared with the mutation pressure. The three narcissus viruses shared the same best optimal codon (CCA) and the synonymous codon prefers to use codons ending with A/U, compared to C/G. Our study shows the codon analysis of different viruses on the same host for the first time, which indicates the importance of the evolutionary-based design to control these viruses.     

Adaptive Evolution as a Driving Force of the Emergence and Re-Emergence of Mosquito-Borne Viral Diseases

Emerging and re-emerging mosquito-borne viral diseases impose a significant burden on global public health. The most common mosquito-borne viruses causing recent epidemics include flaviviruses in the family Flaviviridae, including Dengue virus (DENV), Zika virus (ZIKV), Japanese encephalitis virus (JEV) and West Nile virus (WNV) and Togaviridae viruses, such as chikungunya virus (CHIKV). Several factors may have contributed to the recent re-emergence and spread of mosquito-borne viral diseases. Among these important causes are the evolution of mosquito-borne viruses and the genetic mutations that make them more adaptive and virulent, leading to widespread epidemics. RNA viruses tend to acquire genetic diversity due to error-prone RNA-dependent RNA polymerases, thus promoting high mutation rates that support adaptation to environmental changes or host immunity. In this review, we discuss recent findings on the adaptive evolution of mosquito-borne viruses and their impact on viral infectivity, pathogenicity, vector fitness, transmissibility, epidemic potential and disease emergence.     

Infectious Bronchitis Coronavirus: Genome Evolution in Vaccinated and Non-Vaccinated SPF Chickens

Infectious Bronchitis virus (IBV) continues to cause significant economic losses for the chicken industry despite the use of many live IBV vaccines around the world. Several authors have suggested that vaccine-induced partial protection may contribute to the emergence of new IBV strains. In order to study this hypothesis, three passages of a challenge IBV were made in SPF chickens sham inoculated or vaccinated at day of age using a live vaccine heterologous to the challenge virus. All birds that were challenged with vaccine heterologous virus were positive for viral RNA. NGS analysis of viral RNA in the unvaccinated group showed a rapid selection of seven genetic variants, finally modifying the consensus genome of the viral population. Among them, five were non-synonymous, modifying one position in NSP 8, one in NSP 13, and three in the Spike protein. In the vaccinated group, one genetic variant was selected over the three passages. This synonymous modification was absent from the unvaccinated group. Under these conditions, the genome population of an IBV challenge virus evolved rapidly in both heterologous vaccinated and non-vaccinated birds, while the genetic changes that were selected and the locations of these were very different between the two groups.     

Ongoing Assessment of the Molecular Evolution of Peste Des Petits Ruminants Virus Continues to Question Viral Origins

Understanding the evolution of viral pathogens is critical to being able to define how viruses emerge within different landscapes. Host susceptibility, which is spread between different species and is a contributing factor to the subsequent epidemiology of a disease, is defined by virus detection and subsequent characterization. Peste des petits ruminants virus is a plague of small ruminant species that is a considerable burden to the development of sustainable agriculture across Africa and much of Asia. The virus has also had a significant impact on populations of endangered species in recent years, highlighting its significance as a pathogen of high concern across different regions of the globe. Here, we have re-evaluated the molecular evolution of this virus using novel genetic data to try and further resolve the molecular epidemiology of this disease. Viral isolates are genetically characterized into four lineages (I−IV), and the historic origin of these lineages is of considerable interest to the molecular evolution of the virus. Our re-evaluation of viral emergence using novel genome sequences has demonstrated that lineages I, II and IV likely originated in West Africa, in Senegal (I) and Nigeria (II and IV). Lineage III sequences predicted emergence in either East Africa (Ethiopia) or in the Arabian Peninsula (Oman and/or the United Arab Emirates), with a paucity of data precluding a more refined interpretation. Continual refinements of evolutionary emergence, following the generation of new data, is key to both understanding viral evolution from a historic perspective and informing on the ongoing genetic emergence of this virus.     

Evolution of hepatitis C virus hypervariable region 1 in immunocompetent children born to HCV-infected mothers

Summary.  Hepatitis C virus (HCV) hypervariable region 1 (HVR1) is the most variable region of the viral genome and its heterogeneity reflects the virus-host interplay during chronicity. Paediatric HCV-infected patients develop liver disease with typical clinical features. Here, the evolution of HVR1 and its adjacent regions were ascertained in plasma samples of two HCV-positive children during a 5-year follow-up period. We report an almost complete conservation of the HVR1 amino acid sequence over time, with underlying nucleotide variability both within and outside HVR1, suggesting some kind of constraint on virus evolution, particularly within HVR1. Although overall d _N / d _S rates [rates of nonsynonymous nucleotide substitutions per nonsynonymous site ( d _N ) and synonymous nucleotide substitutions per synonymous site ( d _S )] were <1 in both patients, a high resolution analysis of selection pressures exerted at the codon level revealed few sites subject to selection and an absolute predominance of invariable positions within HVR1. The HVR1 amino acid sequences showed the antigenic properties expected for this region. Taken together, these data suggest peculiar evolutionary dynamics in our patients, which could be attributed to a mechanism of nucleotide invariability along with purifying selection operating on the HVR1. The lack of HVR1 variability may reflect the adaptation of the virus to a particular environment within each patient or a phenomenon of immune tolerance generated in these immunocompetent patients earlier in life.     

Natural Selection Plays an Important Role in Shaping the Codon Usage of Structural Genes of the Viruses Belonging to the Coronaviridae Family

Viruses belonging to the Coronaviridae family have a single-stranded positive-sense RNA with a poly-A tail. The genome has a length of ~29.9 kbps, which encodes for genes that are essential for cell survival and replication. Different evolutionary constraints constantly influence the codon usage bias (CUB) of different genes. A virus optimizes its codon usage to fit the host environment on which it savors. This study is a comprehensive analysis of the CUB for the different genes encoded by viruses of the Coronaviridae family. Different methods including relative synonymous codon usage (RSCU), an Effective number of codons (ENc), parity plot 2, and Neutrality plot, were adopted to analyze the factors responsible for the genetic evolution of the Coronaviridae family. Base composition and RSCU analyses demonstrated the presence of A-ended and U-ended codons being preferred in the 3rd codon position and are suggestive of mutational selection. The lesser ENc value for the spike ‘S’ gene suggests a higher bias in the codon usage of this gene compared to the other structural genes. Parity plot 2 and neutrality plot analyses demonstrate the role and the extent of mutational and natural selection towards the codon usage pattern. It was observed that the structural genes of the Coronaviridae family analyzed in this study were at the least under 84% influence of natural selection, implying a major role of natural selection in shaping the codon usage.     

The atypical codon usage of the plant psbA gene may be the remnant of an ancestral bias

The psbA gene of the chloroplast genome has a codon usage that is unusual for plant chloroplast genes. In the present study the evolutionary status of this codon usage is tested by reconstructing putative ancestral psbA sequences to determine the pattern of change in codon bias during angiosperm divergence. It is shown that the codon biases of the ancestral genes are much stronger than all extant flowering plant psbA genes. This is related to previous work that demonstrated a significant increase in synonymous substitution in psbA relative to other chloroplast genes. It is suggested, based on the two lines of evidence, that the codon bias of this gene currently is not being maintained by selection. Rather, the atypical codon bias simply may be a remnant of an ancestral codon bias that now is being degraded by the mutation bias of the chloroplast genome, in other words, that the psbA gene is not at equilibrium. A model for the evolution of selective pressure on the codon usage of plant chloroplast genes is discussed.     

The Stem-Loop I of Senecavirus A IRES Is Essential for Cap-Independent Translation Activity and Virus Recovery

Senecavirus A (SVA) is a picornavirus that causes vesicular disease in swine and the only member of the Senecavirus genus. Like in all members of Picornaviridae, the 5′ untranslated region (5’UTR) of SVA contains an internal ribosome entry site (IRES) that initiates cap-independent translation. For example, the replacement of the IRES of foot-and-mouth disease virus (FMDV) with its relative bovine rhinitis B virus (BRBV) affects the viral translation efficiency and virulence. Structurally, the IRES from SVA resembles that of hepatitis C virus (HCV), a flavivirus. Given the roles of the IRES in cap-independent translation for picornaviruses, we sought to functionally characterize the IRES of this genus by studying chimeric viruses generated by exchanging the native SVA IRES with that of HCV either entirely or individual domains. First, the results showed that a chimeric SVA virus harboring the IRES from HCV, H-SVA, is viable and replicated normally in rodent-derived BHK-21 cells but displays replication defects in porcine-derived ST cells. In the generation of chimeric viruses in which domain-specific elements from SVA were replaced with those of HCV, we identified an essential role for the stem-loop I element for IRES activity and recombinant virus recovery. Furthermore, a series of stem-loop I mutants allowed us to functionally characterize discrete IRES regions and correlate impaired IRES activities, using reporter systems with our inability to recover recombinant viruses in two different cell types. Interestingly, mutant viruses harboring partially defective IRES were viable. However, no discernable replication differences were observed, relative to the wild-type virus, suggesting the cooperation of additional factors, such as intermolecular viral RNA interactions, act in concert in regulating IRES-dependent translation during infection. Altogether, we found that the stem-loop I of SVA is an essential element for IRES-dependent translation activity and viral replication.     

Constraint of Base Pairing on HDV Genome Evolution

The hepatitis delta virus is a single-stranded circular RNA virus, which is characterized by high self-complementarity. About 70% of the genome sequences can form base-pairs with internal nucleotides. There are many studies on the evolution of the hepatitis delta virus. However, the secondary structure has not been taken into account in these studies. In this study, we developed a method to examine the effect of base pairing as a constraint on the nucleotide substitutions during the evolution of the hepatitis delta virus. The method revealed that the base pairing can reduce the evolutionary rate in the non-coding region of the virus. In addition, it is suggested that the non-coding nucleotides without base pairing may be under some constraint, and that the intensity of the constraint is weaker than that by the base pairing but stronger than that on the synonymous site.     

Evolutionary Patterns of Codon Usage in Major Lineages of Porcine Reproductive and Respiratory Syndrome Virus in China

Porcine reproductive and respiratory syndrome virus (PRRSV) is economically important and characterized by its extensive variation. The codon usage patterns and their influence on viral evolution and host adaptation among different PRRSV strains remain largely unknown. Here, the codon usage of ORF5 genes from lineages 1, 3, 5, and 8, and MLV strains of type 2 PRRSV in China was analyzed. A compositional property analysis of ORF5 genes revealed that nucleotide C is most frequently used at the third position of codons, accompanied by rich GC3s. The effective number of codon (ENC) and codon pair bias (CPB) values indicate that all ORF5 genes have low codon bias and the differences in CPB scores among four lineages are almost not significant. When compared with host codon usage patterns, lineage 1 strains show higher CAI and SiD values, with a high similarity to pig, which might relate to its predominant epidemic propensity in the field. The CAI, RCDI, and SiD values of ORF5 genes from different passages of MLV JXA1R indicate no relation between attenuation and CPB or codon adaptation decrease during serial passage on non-host cells. These findings provide a novel way of understanding the PRRSV’s evolution, related to viral survival, host adaptation, and virulence.     

Identification of two functionally redundant RNA elements in the coding sequence of poliovirus using computer-generated design

Genomes of RNA viruses contain multiple functional RNA elements required for translation or RNA replication. We use unique approaches to identify functional RNA elements in the coding sequence of poliovirus (PV), a plus strand RNA virus. The general method is to recode large segments of the genome using synonymous codons, such that protein sequences, codon use, and codon pair bias are conserved but the nucleic acid sequence is changed. Such recoding does not affect the growth of PV unless it destroys the sequence/structure of a functional RNA element. Using genetic analyses and a method called "signal location search," we detected two unique functionally redundant RNA elements (α and β), each about 75 nt long and separated by 150 nt, in the 3'-terminal coding sequence of RNA polymerase, 3D(pol). The presence of wild type (WT) α or β was sufficient for the optimal growth of PV, but the alteration of both segments in the same virus yielded very low titers and tiny plaques. The nucleotide sequences and predicted RNA structures of α and β have no apparent resemblance to each other. In α, we narrowed down the functional domain to a 48-nt-long, highly conserved segment. The primary determinant of function in β is a stable and highly conserved hairpin. Reporter constructs showed that the α- and β-segments are required for RNA replication. Recoding offers a unique and effective method to search for unknown functional RNA elements in coding sequences of RNA viruses, particularly if the signals are redundant in function.     

Weak selection revealed by the whole-genome comparison of the X chromosome and autosomes of human and chimpanzee

The effect of weak selection driving genome evolution has attracted much attention in the last decade, but the task of measuring the strength of such selection is particularly difficult. A useful approach is to contrast the evolution of X-linked and autosomal genes in two closely related species in a whole-genome analysis. If the fitness effect of mutations is recessive, X-linked genes should evolve more rapidly than autosomal genes when the mutations are advantageous, and they should evolve more slowly than autosomal genes when the mutations are deleterious. We found synonymous substitutions on the X chromosome of human and chimpanzee to be less frequent than those on the autosomes. When calibrated against substitutions in the intergenic regions and pseudogenes to filter out the differences in the mutation rate and ancestral population size between X chromosomes and autosomes, X-linked synonymous substitutions are still 10% less frequent. At least 90% of the synonymous substitutions in human and chimpanzee are estimated to be deleterious, but the fitness effect is weaker than the effect of genetic drift. However, X-linked nonsynonymous substitutions are approximately 30% more frequent than autosomal ones, suggesting the fixation of advantageous mutations that are recessive.     

Comprehensive Analysis of Codon Usage Patterns in Chinese Porcine Circoviruses Based on Their Major Protein-Coding Sequences

Porcine circoviruses (PCVs) are distributed in swine herds worldwide and represent a threat to the health of domestic pigs and the profits of the swine industry. Currently, four PCV species, including PCV-1, PCV-2, PCV-3 and PCV-4, have been identified in China. Considering the ubiquitous characteristic of PCVs, the new emerged PCV-4 and the large scale of swine breeding in China, an overall analysis on codon usage bias for Chinese PCV sequences was performed by using the major proteins coding sequences (ORF1 and ORF2) to better understand the relationship of these viruses with their host. The data from genome nucleotide frequency composition and relative synonymous codon usage (RSCU) analysis revealed an overrepresentation of AT pair and the existence of a certain codon usage bias in all PCVs. However, the values of an effective number of codons (ENC) revealed that the bias was of low magnitude. Principal component analysis, ENC-plot, parity rule two analysis and correlation analysis suggested that natural selection and mutation pressure were both involved in the shaping of the codon usage patterns of PCVs. However, a neutrality plot revealed a stronger effect of natural selection than mutation pressure on codon usage patterns. Good host adaptation was also shown by the codon adaptation index analysis for all these viruses. Interestingly, obtained data suggest that PCV-4 might be more adapted to its host compared to other PCVs. The present study obtained insights into the codon usage pattern of PCVs based on ORF1 and ORF2, which further helps the understanding the molecular evolution of these swine viruses.     

Analysis of Adaptive Evolution in Lyssavirus Genomes Reveals Pervasive Diversifying Selection during Species Diversification

Lyssavirus is a diverse genus of viruses that infect a variety of mammalian hosts, typically causing encephalitis. The evolution of this lineage, particularly the rabies virus, has been a focus of research because of the extensive occurrence of cross-species transmission, and the distinctive geographical patterns present throughout the diversification of these viruses. Although numerous studies have examined pattern-related questions concerning Lyssavirus evolution, analyses of the evolutionary processes acting on Lyssavirus diversification are scarce. To clarify the relevance of positive natural selection in Lyssavirus diversification, we conducted a comprehensive scan for episodic diversifying selection across all lineages and codon sites of the five coding regions in lyssavirus genomes. Although the genomes of these viruses are generally conserved, the glycoprotein (G), RNA-dependent RNA polymerase (L) and polymerase (P) genes were frequently targets of adaptive evolution during the diversification of the genus. Adaptive evolution is particularly manifest in the glycoprotein gene, which was inferred to have experienced the highest density of positively selected codon sites along branches. Substitutions in the L gene were found to be associated with the early diversification of phylogroups. A comparison between the number of positively selected sites inferred along the branches of RABV population branches and Lyssavirus intespecies branches suggested that the occurrence of positive selection was similar on the five coding regions of the genome in both groups.