Historical Contingencies Modulate the Adaptability of Rice Yellow Mottle Virus

The rymv1-2 and rymv1-3 alleles of the RYMV1 resistance to Rice yellow mottle virus (RYMV), coded by an eIF(iso)4G1 gene, occur in a few cultivars of the Asiatic (Oryza sativa) and African (O. glaberrima) rice species, respectively. The most salient feature of the resistance breaking (RB) process is the converse genetic barrier to rymv1-2 and rymv1-3 resistance breakdown. This specificity is modulated by the amino acid (glutamic acid vs. threonine) at codon 49 of the Viral Protein genome-linked (VPg), a position which is adjacent to the virulence codons 48 and 52. Isolates with a glutamic acid (E) do not overcome rymv1-3 whereas those with a threonine (T) rarely overcome rymv1-2. We found that isolates with T49 had a strong selective advantage over isolates with E49 in O. glaberrima susceptible cultivars. This explains the fixation of the mutation T49 during RYMV evolution and accounts for the diversifying selection estimated at codon 49. Better adapted to O. glaberrima, isolates with T49 are also more prone than isolates with E49 to fix rymv1-3 RB mutations at codon 52 in resistant O. glaberrima cultivars. However, subsequent genetic constraints impaired the ability of isolates with T49 to fix rymv1-2 RB mutations at codons 48 and 52 in resistant O. sativa cultivars. The origin and role of the amino acid at codon 49 of the VPg exemplifies the importance of historical contingencies in the ability of RYMV to overcome RYMV1 resistance.     

A general model of codon bias due to GC mutational bias

Background

In spite of extensive research on the effect of mutation and selection on codon usage, a general model of codon usage bias due to mutational bias has been lacking. Because most amino acids allow synonymous GC content changing substitutions in the third codon position, the overall GC bias of a genome or genomic region is highly correlated with GC3, a measure of third position GC content. For individual amino acids as well, G/C ending codons usage generally increases with increasing GC bias and decreases with increasing AT bias. Arginine and leucine, amino acids that allow GC-changing synonymous substitutions in the first and third codon positions, have codons which may be expected to show different usage patterns.

Principal Findings

In analyzing codon usage bias in hundreds of prokaryotic and plant genomes and in human genes, we find that two G-ending codons, AGG (arginine) and TTG (leucine), unlike all other G/C-ending codons, show overall usage that decreases with increasing GC bias, contrary to the usual expectation that G/C-ending codon usage should increase with increasing genomic GC bias. Moreover, the usage of some codons appears nonlinear, even nonmonotone, as a function of GC bias. To explain these observations, we propose a continuous-time Markov chain model of GC-biased synonymous substitution. This model correctly predicts the qualitative usage patterns of all codons, including nonlinear codon usage in isoleucine, arginine and leucine. The model accounts for 72%, 64% and 52% of the observed variability of codon usage in prokaryotes, plants and human respectively. When codons are grouped based on common GC content, 87%, 80% and 68% of the variation in usage is explained for prokaryotes, plants and human respectively.

Conclusions

The model clarifies the sometimes-counterintuitive effects that GC mutational bias can have on codon usage, quantifies the influence of GC mutational bias and provides a natural null model relative to which other influences on codon bias may be measured.     

Inferring the evolutionary history of rice yellow mottle virus from genomic, phylogenetic, and phylogeographic studies

Fourteen isolates of Rice yellow mottle virus (RYMV) were selected as representative of the genetic variability of the virus in Africa from a total set of 320 isolates serologically typed or partially sequenced. The 14 isolates were fully sequenced and analyzed together with two previously reported sequences. RYMV had a genomic organization similar to that of Cocksfoot mottle sobemovirus. The average nucleotide diversity among the 16 isolates of RYMV was 7%, and the maximum diversity between any two isolates was 10%. A strong conservative selection was apparent on both synonymous and nonsynonymous substitutions, through the amino acid replacement pattern, on the genome size, and through the limited number of indel events. Furthermore, there was a lack of positive selection on single amino acid sites and no evidence of recombination events. RYMV diversity had a pronounced and characteristic geographic structure. The branching order of the clades correlated with the geographic origin of the isolates along an east-to-west transect across Africa, and there was a marked decrease in nucleotide diversity moving westward across the continent. The insertion-deletion polymorphism was related to virus phylogeny. There was a partial phylogenetic incongruence between the coat protein gene and the rest of the genome. Overall, our results support the hypothesis that RYMV originated in East Africa and then dispersed and differentiated gradually from the east to the west of the continent.     

An extreme codon preference strategy: codon reassignment

We argue that in animal mitochondria codon reassignments, such as those for AGA and AGG from arginine to serine or of AUA from isoleucine to methionine, are the result of an interplay between biased mutational forces and selective ones. In particular, there is a marked tendency for animal mitochondria to have very small genomes and to minimize their investment in components required for gene expression. These tendencies are expressed as a reduction in the diversity of tRNA isoacceptor species. In our view, the pressure to simplify tRNA populations, together with mutational bias against certain codons, will account for the codon reassignments observed in animal mitochondria. A parallel to the major codon bias in microorganisms, which likewise tends to reduce the diversity of the tRNA isoacceptor populations under fast growth conditions, may be drawn. Therefore, we suggest that codon reassignments are usefully viewed as an extreme form of codon bias.     

Mutations in pea seedborne mosaic virus genome-linked protein VPg after pathotype-specific virulence in Pisum sativum

Pisum sativum plant introduction (PI) line 269818 is resistant to potyvirus pea seedborne mosaic virus (PSbMV) isolates, categorized as pathotype P1, and is susceptible to pathotype P4 isolates. This difference in infectivity is determined by the viral genome-linked protein (VPg) cistron. Mutational analysis of VPg of PSbMV isolates DPD1 and NY representing pathotypes P1 and P4 revealed that codon changes affecting amino acids 105 to 117 in the central region of VPg influenced virulence on PI 269818. In contrast, infectivity on pea cultivar Dark Skinned Perfection, which is susceptible to both pathotypes, was not affected by the mutations. Mutants overcoming resistance in PI 269818 were analyzed for changes in the VPg coding region upon passage through PI 269818 and Dark Skinned Perfection. Adaptive changes were observed only upon passage through PI 269818 and only at codons from amino acid 105 to 117. Expression of DPD1 VPg in PI 269818 did not affect infection by NY, which suggests that VPg from DPD1 is not an elicitor of a general resistance response. The results are compatible with the hypothesis that viral amplification depends upon the interaction between VPg and a host factor.     

The consequences of base-pair substitution mutations in AT- and GC-rich bacteria

The likely consequences, in terms of premature stop codons, detectable missense mutants, silent missense mutants, and degenerate codon changes, have been determined for all 12 individual base substitution changes. This has been done for the full, 61 sense codon, genetic code and also for the much more limited codon availabilities of AT- or GC-rich DNA. The specificities and outcomes of individual base substitutions are likely to be rather different at AT- or GC-rich extremes, and also from the situation at an intermediate DNA base-ratio where all 61 sense codons are available. In particular, at DNA base-ratio extremes many mutations will be to non-utilized codons, which may well act as nonsense mutants. These in turn will give novel classes of suppressor-containing revertants. Even in bacteria with intermediate DNA base-ratios, particular codons for a given amino acid may be favoured, over alternatives, because their use maximizes, or minimizes, the mutational consequences of one, or more, base substitution changes.     

Mutations in potato virus Y genome-linked protein determine virulence toward recessive resistances in Capsicum annuum and Lycopersicon hirsutum

The recessive resistance genes pot-1 and pvr2 in Lycopersicon hirsutum and Capsicum annuum, respectively, control Potato virus Y (PVY) accumulation in the inoculated leaves. Infectious cDNA molecules from two PVY isolates differing in their virulence toward these resistances were obtained using two different strategies. Chimeras constructed with these cDNA clones showed that a single nucleotide change corresponding to an amino acid substitution (Arg119His) in the central part of the viral protein genome-linked (VPg) was involved in virulence toward the pot-1 resistance. On the other hand, 15 nucleotide changes corresponding to five putative amino acid differences in the same region of the VPg affected virulence toward the pvr2(1) and pvr2(2) resistances. Substitution models identified six and five codons within the central and C terminal parts of the VPg for PVY and for the related potyvirus Potato virus A, respectively, which undergo positive selection. This suggests that the role of the VPg-encoding region is determined by the protein and not by the viral RNA apart from its protein-encoding capacity.     

Novel in-frame two codon translational hop during synthesis of bovine placental lactogen in a recombinant strain of Escherichia coli.

A recombinant Escherichia coli strain was constructed for the overexpression of bovine placental lactogen (bPL), using a bPL structural gene containing 9 of the rare arginine codons AGA and AGG. When high level bPL synthesis was induced in this strain, cell growth was inhibited and bPL accumulated to less than 10% of total cell protein. In addition, about 2% of the recombinant bPL produced from this strain exhibited an altered trypsin digestion pattern. Amino acid residues 74 through 109 normally produce 2 tryptic peptides, but the altered form of bPL lacked these two peptides and instead had a new peptide which was missing arginine residue 86 and one of the two flanking leucine residues. The codon for arginine residue 86 was AGG and the codons for the flanking leucine residues 85 and 87 were TTG. When 5 of the 9 AGA and AGG codons in the bPL structural gene were changed to more preferred arginine codons, cell growth was not inhibited and bPL accumulated to about 30% of total cell protein. When bPL was purified from this modified strain, which included changing the arginine codon at position 86 from AGG to CGT, none of the altered form of bPL was produced. These observations are consistent with a model in which translational pausing occurs at the arginine residue 86 AGG codon because the corresponding arginyl-tRNA species is reduced by the high level of bPL synthesis, and a translational hop occurs from the leucine residue 85 TTG codon to the leucine residue 87 TTG codon. This observation represents the first report of an error in protein synthesis due to an in-frame translational hop within an open reading frame.     

Molecular characteristics and evolutionary analysis of a very virulent infectious bursal disease virus

Infectious bursal disease virus(IBDV) poses a significant threat to the poultry industry. Viral protein 2(VP2), the major structural protein of IBDV, has been subjected to frequent mutations that have imparted tremendous genetic diversity to the virus. To determine how amino acid mutations may affect the virulence of IBDV, we built a structural model of VP2 of a very virulent strain of IBDV identified in China, vv IBDV Gx, and performed a molecular dynamics simulation of the interaction between virulence sites. The study showed that the amino acid substitutions that distinguish vv IBDV from attenuated IBDV(H253Q and T284A) favor a hydrophobic and flexible conformation of ?-barrel loops in VP2, which could promote interactions between the virus and potential IBDV-specific receptors. Population sequence analysis revealed that the IBDV strains prevalent in East Asia show a significant signal of positive selection at virulence sites 253 and 284. In addition, a signal of co-evolution between sites 253 and 284 was identified. These results suggest that changes in the virulence of IBDV may result from both the interaction and the co-evolution of multiple amino acid substitutions at virulence sites.     

Synonymous Codon Usage,GC<Subscript>3</Subscript>, and Evolutionary Patterns Across Plastomes of Three Pooid Model Species: Emerging Grass Genome Models for Monocots

We have analyzed factors affecting the codon usage pattern of the chloroplasts genomes of representative species of pooid grass family. Correspondence analysis of relative synonymous codon usages (RSCU) showed that genes on secondary axis were correlated with their GC_3S values (all r > 0.3, p < 0.05), indicating mutational bias as an important selective force that shaped the variation in the codon usage among chloroplast genes. The Nc-plot showed that although a majority of the points with low-Nc values were lying below the expected curve, a few genes lied on the expected curve. Nc plot clearly showed that mutational bias plays a major role in codon biology across the monocot plastomes. The hydrophobicity and aromaticity of encoded proteins of each species were found to be other factors of codon usage variation. In the view of above light, besides natural selection, several other factors also likely to be involved in determining the selective constraints on codon bias in plastomes of pooid grass genomes. In addition, five codons (B. distachyon), seven codons (H. vulgare), and four codons (T. aestivum) were identified as optimal codons of the three grass chloroplasts. To identify genes evolving under positive selection, rates of nonsynonymous substitutions (Ka) and synonymous substitutions (Ks) were computed for all groups of orthologous gene pairs.     

Translationally optimal codons associate with aggregation-prone sites in proteins

We analyze the relationship between codon usage bias and residue aggregation propensity in the genomes of four model organisms, Escherichia coli, yeast, fly, and mouse, as well as the archaeon Halobacterium species NRC-1. Using the Mantel-Haenszel procedure, we find that translationally optimal codons associate with aggregation-prone residues. Our results are qualitatively and quantitatively similar to those of an earlier study where we found an association between translationally optimal codons and buried residues. We also combine the aggregation-propensity data with solvent-accessibility data. Although the resulting data set is small, and hence statistical power low, results indicate that the association between optimal codons and aggregation-prone residues exists both at buried and at exposed sites. By comparing codon usage at different combinations of sites (exposed, aggregation-prone sites versus buried, non-aggregation-prone sites; buried, aggregation-prone sites versus exposed, non-aggregation-prone sites), we find that aggregation propensity and solvent accessibility seem to have independent effects of (on average) comparable magnitude on codon usage. Finally, in fly, we assess whether optimal codons associate with sites at which amino acid substitutions lead to an increase in aggregation propensity, and find only a very weak effect. These results suggest that optimal codons may be required to reduce the frequency of translation errors at aggregation-prone sites that coincide with certain functional sites, such as protein-protein interfaces. Alternatively, optimal codons may be required for rapid translation of aggregation-prone regions.     

Synonymous codon usage pattern among the S,M, and L segments in Crimean-congo hemorrhagic fever virus

Crimean-Congo hemorrhagic fever (CCHF) virus is one among the major zoonosis viral diseases that use the Hyalomma ticks as their transmission vector to cause viral infection to the human and mammalian community. The fatality of infectious is high across the world especially in Africa, Asia, Middle East, and Europe. This study regarding codon usage bias of S, M, and L segments of the CCHF virus pertaining to the host Homo sapiens, reveals in-depth information about the evolutionary characteristics of CCHFV. Relative Synonymous Codon Usage (RSCU), Effective number of codons (ENC) were calculated, to determine the codon usage pattern in each segment. Correlation analysis between Codon adaptation index (CAI), GRAVY (Hydrophobicity), AROMO (Aromaticity), and nucleotide composition revealed bias in the codon usage pattern. There was no strong codon bias found among any segments of the CCHF virus, indicating both the factors i.e., natural selection and mutational pressure shapes the codon usage bias.     

Sequences of E/NS1 Gene Junction from Four Dengue-2 Viruses of Northeastern Thailand and Their Evolutionary Relationships with Other Dengue-2 Viruses

We determined the 240-nueleotide sequences of the E/NS1 gene junction of four dengue-2 viruses by the primer extension dideoxy chain termination method. These viruses were isolated from dengue patients with different clinical severities in Nakhon Phanom, Northeastern Thailand in 1993. The results were compared with the 52 published dengue-2 sequences of the same gene region. Sequence divergence of four new isolates varied from 4.17% to 5.42% compared with dengue-2 prototype New Guinea C strain whereas it varied from 5.42% to 6.67% and from 6.67% to 7.09% when compared with Jamaica 1409 strain and PR159/S1 strain, respectively. All nucleotide substitutions were found at the 3rd position of the codons which were silent mutations. All 56 isolates studied were classified into five genotypic groups by constructing the dendrogram. The results indicated that four new isolates from Northeastern Thailand belong to genotype II of dengue virus serotype 2, and were most closely related to prototype New Guinea C strain. We also observed the variation in nucleotide and amino acid sequences among clusters of isolates (Thailand-1980, Malaysia-1989 and Thailand-1993) which were obtained from the dengue patients with different clinical severities. The significance of these genetic differences have been discussed in terms of the possible correlation between genetic variability and virulence.     

Study of Completed Archaeal Genomes and Proteomes： Hypothesis of Strong Mutational AT Pressure Existed in Their Common Predecessor

The number of completely sequenced archaeal genomes has been sufficient for a large-scale bioinformatic study.We have conducted analyses for each coding region from 36 archaeal genomes using the original CGS algorithm by calculating the total GC content(G+C),GC content in first,second and third codon positions as well as in fourfold and twofold degenerated sites from third codon positions,levels of arginine codon usage(Arg2:AGA/G;Arg4:CGX),levels of amino acid usage and the entropy of amino acid content distribution.In archaeal genomes with strong GC pressure,arginine is coded preferably by GC-rich Arg4 codons,whereas in most of archaeal genomes with G+C0.6,arginine is coded preferably by AT-rich Arg2 codons.In the genome of Haloquadratum walsbyi,which is closely related to GC-rich archaea,GC content has decreased mostly in third codon positions,while Arg4Arg2 bias still persists.Proteomes of archaeal species carry characteristic amino acid biases:levels of isoleucine and lysine are elevated,while levels of alanine,histidine,glutamine and cytosine are relatively decreased.Numerous genomic and proteomic biases observed can be explained by the hypothesis of previously existed strong mutational AT pressure in the common predecessor of all archaea.     

Role of nonstructural proteins 3A and 3B in host range and pathogenicity of foot-and-mouth disease virus

The genome of foot-and-mouth disease virus (FMDV) differs from that of other picornaviruses in that it encodes a larger 3A protein (>50% longer than poliovirus 3A), as well as three copies of protein 3B (also known as VPg). Previous studies have shown that a deletion of amino acids 93 to 102 of the 153-codon 3A protein is associated with an inability of a Taiwanese strain of FMDV (O/TAW/97) to cause disease in bovines. Recently, an Asian virus with a second 3A deletion (amino acids 133 to 143) has also been detected (N. J. Knowles et al., J. Virol. 75:1551-1556, 2001). Genetically engineered viruses harboring the amino acids 93 to 102 or 133 to 143 grew well in porcine cells but replicated poorly in bovine cells, whereas a genetically engineered derivative of the O/TAW/97 virus expressing a full-length 3A (strain A12) grew well in both cell types. Interestingly, a virus with a deletion spanning amino acid 93 to 144 also grew well in porcine cells and caused disease in swine. Further, genetically engineered viruses containing only a single copy of VPg were readily recovered with the full-length 3A, the deleted 3A (amino acids 93 to 102), or the "super" deleted forms of 3A (missing amino acids 93 to 144). All of the single-VPg viruses were attenuated in porcine cells and replicated poorly in bovine cells. The single-VPg viruses produced a mild disease in swine, indicating that the VPg copy number is an important determinant of host range and virulence. The association of VPg copy number with increased virulence in vivo may help to explain why all naturally occurring FMDVs have retained three copies of VPg.     

Structure of native and expanded sobemoviruses by electron cryo-microscopy and image reconstruction

Rice yellow mottle virus (RYMV) and southern bean mosaic virus, cowpea strain (SCPMV) are members of the Sobemovirus genus of RNA-containing viruses. We used electron cryo-microscopy (cryo-EM) and icosahedral image analysis to examine the native structures of these two viruses at 25 A resolution. Both viruses have a single tightly packed capsid layer with 180 subunits assembled on a T=3 icosahedral lattice. Distinctive crown-like pentamers emanate from the 12 5-fold axes of symmetry. The exterior face of SCPMV displays deep valleys along the 2-fold axes and protrusions at the quasi-3-fold axes. While having a similar topography, the surface of RYMV is comparatively smooth. Two concentric shells of density reside beneath the capsid layer of RYMV and SCPMV, which we interpret as ordered regions of genomic RNA. In the presence of divalent cations, SCPMV particles swell and fracture, whereas the expanded form of RYMV is stable. We previously proposed that the cell-to-cell movement of RYMV in xylem involves chelation of Ca(2+) from pit membranes of infected cells, thereby stabilizing the capsid shells and allowing a pathway for spread of RYMV through destabilized membranes. In the context of this model, we propose that the expanded form of RYMV is an intermediate in the in vivo assembly of virions.     

Challenges of the genetic code for exploring sequence space in directed protein evolution

Directed protein evolution is the most versatile method for studying protein structure-function relationships, and for tailoring a protein's properties to the needs of industrial applications. In this review, we performed a statistical analysis on the genetic code to study the extent and consequence of the organization of the genetic code on amino acid substitution patterns generated in directed evolution experiments. In detail, we analyzed amino acid substitution patterns caused by (a) a single nucleotide (nt) exchange at each position of all 64 codons, and (b) two subsequent nt exchanges (first and second nt, first and third nt, second and third nt). Additionally, transitions and transversions mutations were compared at the level of amino acid substitution patterns. The latter analysis showed that single nucleotide substitution in a codon generates only 39.5% of the natural diversity on the protein level with 5.2-7 amino acid substitutions per codon. Transversions generate more complex amino acid substitution patterns (increased number and chemically more diverse amino acid substitutions) than transitions. Simultaneous nt exchanges at both first and second nt of a codon generates very diverse amino acid substitution patterns, achieving 83.2% of the natural diversity. The statistical analysis described in this review sets the objectives for novel random mutagenesis methods that address the consequences of the organization of the genetic code. Random mutagenesis methods that favor transversions or introduce consecutive nt exchanges can contribute in this regard.     

Different mutations in the genome-linked protein VPg of potato virus Y confer virulence on the pvr2(3) resistance in pepper

Five different amino acid substitutions in the VPg of Potato virus Y were shown to be independently responsible for virulence toward pvr2(3) resistance gene of pepper. A consequence of these multiple mutations toward virulence involving single nucleotide substitutions is a particularly high frequency of resistance breaking (37% of inoculated plants from the first inoculation) and suggests a potentially low durability of pvr2(3) resistance. These five mutants were observed with significantly different frequencies, one of them being overrepresented. Genetic drift alone could not explain the observed distribution of virulent mutants. More plausible scenarios were obtained by taking into account either the relative substitution rates, the relative fitness of the mutants in pvr2(3) pepper plants, or both.