Evolutionary basis of codon usage and nucleotide composition bias in vertebrate DNA viruses

Rationalizing the development of live attenuated virus vaccines

The design of vaccines against viral disease has evolved considerably over the past 50 years. Live attenuated viruses (LAVs)-those created by passaging a virus in cultured cells-have proven to be an effective means for preventing many viral diseases, including smallpox, polio, measles, mumps and yellow fever. Even so, empirical attenuation is unreliable in some cases and LAVs pose several safety issues. Although inactivated viruses and subunit vaccines alleviate many of these concerns, they have in general been less efficacious than their LAV counterparts. Advances in molecular virology--creating deleterious gene mutations, altering replication fidelity, deoptimizing codons and exerting control by microRNAs or zinc finger nucleases--are providing new ways of controlling viral replication and virulence and renewing interest in LAV vaccines. Whereas these rationally attenuated viruses may lead to a new generation of safer, more widely applicable LAV vaccines, each approach requires further testing before progression to human testing.

Human alpha and beta papillomaviruses use different synonymous codon profiles

Human papillomaviruses use rare codons relative to their hosts. It has been theorized that this is a mechanism to allow the virus to escape immune surveillance. In the present study, we examined the codings of four major genes of 21 human alpha (mucosatropic) viruses and 16 human beta (cutaneous-tropic) viruses. We compared the codon usage of different genes from a given papillomavirus and also the same genes from different papillomaviruses. Our data showed that codon usage was not always uniform between two genes of a given papillomavirus or between the same genes of papillomaviruses from different genera. We speculate as to why this might be and conclude that codon usage in the papillomaviruses may not only play a role in facilitating escape from immune surveillance but may also underlie some of the unanswered questions in the papillomavirus field.

Categorizing host-dependent RNA viruses by principal component analysis of their codon usage preferences

Viruses have to exploit host transcription and translation mechanisms to replicate in a hostile host cellular environment, and therefore, it is likely that the infected host may impose pressure on viral evolution. In this study, we investigated differences in codon usage preferences among the highly mutable single strain RNA viruses which infect vertebrate or invertebrate hosts, respectively. We incorporate principal component analysis (PCA) and k-mean methods to clustering viruses infected with different type of hosts. The relative synonymous codon usage (RSCU) indices of all genes in 32 RNA viruses were calculated, and the correlation of the RSCU indices among different viruses was analyzed by the PCA. Our results show a positive correlation in codon usage preferences among viruses that target the same host category. Results of k-means clustering analysis further confirmed the statistical significance of this study, demonstrating that viruses infecting vertebrate hosts have different codon usage preferences to those of invertebrate viruses. Based on the analysis of the effective number of codons (ENC) in relation to the GC-content at the synonymous third codon position (GC3s), we further identified that mutational pressure was the dominant evolution driving force in making the different codon usage preferences. This study suggests a new and effective way to characterize host-dependent RNA viruses based on the codon usage pattern.

Analysis of synonymous codon usage in porcine reproductive and respiratory syndrome virus

In this study, we calculated the relative synonymous codon usage (RSCU) values and codon usage bias (CUB) values to implement a comparative analysis of codon usage pattern of open reading frames (ORFs) which belong to the two main genotypes of porcine reproductive and respiratory syndrome virus (PRRSV). By analysis of synonymous codon usage values in each ORF of PRRSV, the optimal codons for most amino acids were all C or G-ended codons except GAU for Asp, CAU for His, UUU for Phe and CCU for Pro. The synonymous codon usage patterns in different ORFs of PRRSV were different and genetically conserved. Among them, ORF1a, ORF4, ORF5 and ORF7 could cluster these strains into the two main serotypes (EU and US). Due to mutational pressure, compositional constraint played an important role in shaping the synonymous codon usage pattern in different ORFs, and the synonymous codon usage diversity in ORFs was correlated with gene function. The degree of CUB for some particular amino acids under strong selection pressure probably served as a potential genetic marker for each ORF in PRRSV. However, gene length and translational selection in nature had no effect on the synonymous codon usage pattern in PRRSV. These conclusions could not only offer an insight into the synonymous codon usage pattern and differentiation of gene function, but also assist in understanding the discrepancy of evolution among ORFs in PRRSV.

Analysis of synonymous codon usage in foot-and-mouth disease virus

In this study, we calculate the relative synonymous codon usage (RSCU) values and codon usage bias (CUB) values to carry out a comparative analysis of codon usage pattern for open reading frames (ORFs) among 85 samples which belong to all seven serotypes of foot-and-mouth disease virus (FMDV). Although the degree of CUB for ORFs is a relatively slight, there is a significant variation for CUB among different serotypes, which is mainly determined by codon usage pattern depending on RSCU. By comparison with RSCU values for all samples, although RSCU values fail to show the relationship of specific-lineage serotype, there are two main genetic populations existing in FMDV, namely (i) serotypes Asia 1, A, C & O; (ii) serotypes SAT 1, 2 & 3. This interesting characteristic may be formed by the mechanism of RNA virus recombination. The analysis of quantitative & qualitative evaluation based on CUB indicates interesting characteristic of codon usage, which suggests that more FMDV genome diversity may exist in specific-lineage serotypes rather than exist randomly. Furthermore, the relationship between amino acids and codon usage pattern indicates that mutation pressure rather than translational selection in nature is the important determinant of the codon usage bias observed. Our work might give some sight into some characteristics of FMDV ORF and some evolutionary information of this virus.

Viral adaptation to host: a proteome-based analysis of codon usage and amino acid preferences

Viruses differ markedly in their specificity toward host organisms. Here, we test the level of general sequence adaptation that viruses display toward their hosts. We compiled a representative data set of viruses that infect hosts ranging from bacteria to humans. We consider their respective amino acid and codon usages and compare them among the viruses and their hosts. We show that bacteria-infecting viruses are strongly adapted to their specific hosts, but that they differ from other unrelated bacterial hosts. Viruses that infect humans, but not those that infect other mammals or aves, show a strong resemblance to most mammalian and avian hosts, in terms of both amino acid and codon preferences. In groups of viruses that infect humans or other mammals, the highest observed level of adaptation of viral proteins to host codon usages is for those proteins that appear abundantly in the virion. In contrast, proteins that are known to participate in host-specific recognition do not necessarily adapt to their respective hosts. The implication for the potential of viral infectivity is discussed.

Large-scale inference of the point mutational spectrum in human segmental duplications

BACKGROUND

Recent segmental duplications are relatively large (> or = 1 kb) genomic regions of high sequence identity (> or = 90%). They cover approximately 4-5% of the human genome and play important roles in gene evolution and genomic disease. The DNA sequence differences between copies of a segmental duplication represent the result of various mutational events over time, since any two duplication copies originated from the same ancestral DNA sequence. Based on this fact, we have developed a computational scheme for inference of point mutational events in human segmental duplications, which we collectively term duplication-inferred mutations (DIMs). We have characterized these nucleotide substitutions by comparing them with high-quality SNPs from dbSNP, both in terms of sequence context and frequency of substitution types.

RESULTS

Overall, DIMs show a lower ratio of transitions relative to transversions than SNPs, although this ratio approaches that of SNPs when considering DIMs within most recent duplications. Our findings indicate that DIMs and SNPs in general are caused by similar mutational mechanisms, with some deviances at the CpG dinucleotide. Furthermore, we discover a large number of reference SNPs that coincide with computationally inferred DIMs. The latter reflects how sequence variation in duplicated sequences can be misinterpreted as ordinary allelic variation.

CONCLUSION

In summary, we show how DNA sequence analysis of segmental duplications can provide a genome-wide mutational spectrum that mirrors recent genome evolution. The inferred set of nucleotide substitutions represents a valuable complement to SNPs for the analysis of genetic variation and point mutagenesis.

Family Iridoviridae: poor viral relations no longer

Members of the family Iridoviridae infect a diverse array of invertebrate and cold-blooded vertebrate hosts and are currently viewed as emerging pathogens of fish and amphibians. Iridovirid replication is unique and involves both nuclear and cytoplasmic compartments, a circularly permuted, terminally redundant genome that, in the case of vertebrate iridoviruses, is also highly methylated, and the efficient shutoff of host macromolecular synthesis. Although initially neglected largely due to the perceived lack of health, environmental, and economic concerns, members of the genus Ranavirus, and the newly recognized genus Megalocytivirus, are rapidly attracting growing interest due to their involvement in amphibian population declines and their adverse impacts on aquaculture. Herein we describe the molecular and genetic basis of viral replication, pathogenesis, and immunity, and discuss viral ecology with reference to members from each of the invertebrate and vertebrate genera.

Analysis of synonymous codon usage and evolution of begomoviruses

Begomoviruses are single-stranded DNA viruses and cause severe diseases in major crop plants worldwide. Based on current genome sequence analyses, we found that synonymous codon usage variations in the protein-coding genes of begomoviruses are mainly influenced by mutation bias. Base composition analysis suggested that the codon usage bias of AV1 and BV1 genes is significant and their expressions are high. Fourteen codons were determined as translational optimal ones according to the comparison of codon usage patterns between highly and lowly expressed genes. Interestingly the codon usages between begomoviruses from the Old and the New Worlds are apparently different, which supports the idea that the bipartite begomoviruses of the New World might originate from bipartite ones of the Old World, whereas the latter evolve from the Old World monopartite begomoviruses.

Analysis of synonymous codon usage in classical swine fever virus

Using the complete genome sequences of 35 classical swine fever viruses (CSFV) representing all three genotypes and all three kinds of virulence, we analyzed synonymous codon usage and the relative dinucleotide abundance in CSFV. The general correlation between base composition and codon usage bias suggests that mutational pressure rather than natural selection is the main factor that determines the codon usage bias in CSFV. Furthermore, we observed that the relative abundance of dinucleotides in CSFV is independent of the overall base composition but is still the result of differential mutational pressure, which also shapes codon usage. In addition, other factors, such as the subgenotypes and aromaticity, also influence the codon usage variation among the genomes of CSFV. This study represents the most comprehensive analysis to date of CSFV codon usage patterns and provides a basic understanding of the mechanisms for codon usage bias.

CRPV genomes with synonymous codon optimizations in the CRPV E7 gene show phenotypic differences in growth and altered immunity upon E7 vaccination

Papillomaviruses use rare codons relative to their hosts. Recent studies have demonstrated that synonymous codon changes in viral genes can lead to increased protein production when the codons are matched to those of cells in which the protein is being expressed. We theorized that the immunogenicity of the virus would be enhanced by matching codons of selected viral genes to those of the host. We report here that synonymous codon changes in the E7 oncogene are tolerated in the context of the cottontail rabbit papillomavirus (CRPV) genome. Papilloma growth rates differ depending upon the changes made indicating that synonymous codons are not necessarily neutral. Immunization with wild type E7 DNA yielded significant protection from subsequent challenge by both wild type and codon-modified genomes. The reduction in growth was most dramatic with the genome containing the greatest number of synonymous codon changes.

Patterns of evolution and host gene mimicry in influenza and other RNA viruses

It is well known that the dinucleotide CpG is under-represented in the genomic DNA of many vertebrates. This is commonly thought to be due to the methylation of cytosine residues in this dinucleotide and the corresponding high rate of deamination of 5-methycytosine, which lowers the frequency of this dinucleotide in DNA. Surprisingly, many single-stranded RNA viruses that replicate in these vertebrate hosts also have a very low presence of CpG dinucleotides in their genomes. Viruses are obligate intracellular parasites and the evolution of a virus is inexorably linked to the nature and fate of its host. One therefore expects that virus and host genomes should have common features. In this work, we compare evolutionary patterns in the genomes of ssRNA viruses and their hosts. In particular, we have analyzed dinucleotide patterns and found that the same patterns are pervasively over- or under-represented in many RNA viruses and their hosts suggesting that many RNA viruses evolve by mimicking some of the features of their host's genes (DNA) and likely also their corresponding mRNAs. When a virus crosses a species barrier into a different host, the pressure to replicate, survive and adapt, leaves a footprint in dinucleotide frequencies. For instance, since human genes seem to be under higher pressure to eliminate CpG dinucleotide motifs than avian genes, this pressure might be reflected in the genomes of human viruses (DNA and RNA viruses) when compared to those of the same viruses replicating in avian hosts. To test this idea we have analyzed the evolution of the influenza virus since 1918. We find that the influenza A virus, which originated from an avian reservoir and has been replicating in humans over many generations, evolves in a direction strongly selected to reduce the frequency of CpG dinucleotides in its genome. Consistent with this observation, we find that the influenza B virus, which has spent much more time in the human population, has adapted to its human host and exhibits an extremely low CpG dinucleotide content. We believe that these observations directly show that the evolution of RNA viral genomes can be shaped by pressures observed in the host genome. As a possible explanation, we suggest that the strong selection pressures acting on these RNA viruses are most likely related to the innate immune response and to nucleotide motifs in the host DNA and RNAs.

Viruses are obligate intracellular parasites that use different strategies to sequester host cell machinery and avoid the host immune system. In this paper we explore the genomes of viruses that encode their genetic information in single-stranded RNA, a different material than the one used by their hosts (double-stranded DNA). It is interesting to observe that these viruses share some of the host's characteristics. For instance, one of the most underrepresented motifs in the DNA of vertebrates is the dinucleotide CpG. This is commonly thought to be due to methylation and deamination of cytosine residues in this dinucleotide. Surprisingly, the same CpG suppression is observed in vertebrate RNA viruses but not in RNA phages. We show that RNA viruses present similar dinucleotide pressures as their host genes. We find that the influenza A virus, which originated from an avian reservoir and replicated in humans over many generations, evolves to reduce the frequency of CpG dinucleotides mimicking the human genes. Influenza B, which has been in humans longer, exhibits an extremely low CpG dinucleotide content. These observations suggest that the evolution of RNA viruses is shaped by pressures observed in the host genome.

Genome-scale relationships between cytosine methylation and dinucleotide abundances in animals

In mammalian genomes CpGs occur at one-fifth their expected frequency. This is accepted as resulting from cytosine methylation and deamination of 5-methylcytosine leading to TpG and CpA dinucleotides. The corollary that a CpG deficit should correlate with TpG excess has not hitherto been systematically tested at a genomic level. I analyzed genome sequences (human, chimpanzee, mouse, pufferfish, zebrafish, sea squirt, fruitfly, mosquito, and nematode) to do this and generally to assess the hypothesis that CpG deficit, TpG excess, and other data are accountable in terms of 5-methylcytosine mutation. In all methylated genomes local CpG deficit decreases with higher G + C content. Local TpG surplus, while positively associated with G + C level in mammalian genomes but negatively associated with G + C in nonmammalian methylated genomes, is always explicable in terms of the CpG trend under the methylation model. Covariance of dinucleotide abundances with G + C demonstrates that correlation analyses should control for G + C. Doing this reveals a strong negative correlation between local CpG and TpG abundances in methylated genomes, in accord with the methylation hypothesis. CpG deficit also correlates with CpT excess in mammals, which may reflect enhanced cytosine mutation in the context 5'-YCG-3'. Analyses with repeat-masked sequences show that the results are not attributable to repetitive elements.

Presence and role of cytosine methylation in DNA viruses of animals

Nucleotide composition varies greatly among DNA viruses of animals, yet the evolutionary pressures and biological mechanisms driving these patterns are unclear. One of the most striking discrepancies lies in the frequency of CpG (the dinucleotide CG, linked by a phosphate group), which is underrepresented in most small DNA viruses (those with genomes below 10 kb) but not in larger DNA viruses. Cytosine methylation might be partially responsible, but research on this topic has focused on a few virus groups. For several viruses that integrate their genome into the host genome, the methylation status during this stage has been studied extensively, and the relationship between methylation and viral-induced tumor formation has been examined carefully. However, for actively replicating viruses—particularly small DNA viruses—the methylation status of CpG motifs is rarely known and the effects on the viral life cycle are obscure. In vertebrate host genomes, most cytosines at CpG sites are methylated, which in vertebrates acts to regulate gene expression and facilitates the recognition of unmethylated, potentially pathogen-associated DNA. Here we briefly introduce cytosine methylation before reviewing what is currently known about CpG methylation in DNA viruses.

Analysis of synonymous codon usage in 11 human bocavirus isolates

Human Bocavirus (HBoV) is a novel virus which can cause respiratory tract disease in infants or children. In this study, the codon usage bias and the base composition variations in the available 11 complete HBoV genome sequences have been investigated. Although, there is a significant variation in codon usage bias among different HBoV genes, codon usage bias in HBoV is a little slight, which is mainly determined by the base compositions on the third codon position and the effective number of codons (ENC) value. The results of correspondence analysis (COA) and Spearman's rank correlation analysis reveals that the G + C compositional constraint is the main factor that determines the codon usage bias in HBoV and the gene's function also contributes to the codon usage in this virus. Moreover, it was found that the hydrophobicity of each protein and the gene length are also critical in affecting these viruses’ codon usage, although they were less important than that of the mutational bias and the genes’ function. At last, the relative synonymous codon usage (RSCU) of 44 genes from these 11 HBoV isolates is analyzed using a hierarchical cluster method. The result suggests that genes with same function yet from different isolates are classified into the same lineage and it does not depend on geographical location. These conclusions not only can offer an insight into the codon usage patterns and gene classification of HBoV, but also may help in increasing the efficiency of gene delivery/expression systems.

Parvoviral host range and cell entry mechanisms

Parvoviruses elaborate rugged nonenveloped icosahedral capsids of approximately 260 A in diameter that comprise just 60 copies of a common core structural polypeptide. While serving as exceptionally durable shells, capable of protecting the single-stranded DNA genome from environmental extremes, the capsid also undergoes sequential conformational changes that allow it to translocate the genome from its initial host cell nucleus all the way into the nucleus of its subsequent host. Lacking a duplex transcription template, the virus must then wait for its host to enter S-phase before it can initiate transcription and usurp the cell's synthetic pathways. Here we review cell entry mechanisms used by parvoviruses. We explore two apparently distinct modes of host cell specificity, first that used by Minute virus of mice, where subtle glycan-specific interactions between host receptors and residues surrounding twofold symmetry axes on the virion surface mediate differentiated cell type target specificity, while the second involves novel protein interactions with the canine transferrin receptor that allow a mutant of the feline leukopenia serotype, Canine parvovirus, to bind to and infect dog cells. We then discuss conformational shifts in the virion that accompany cell entry, causing exposure of a capsid-tethered phospholipase A2 enzymatic core that acts as an endosomolytic agent to mediate virion translocation across the lipid bilayer into the cell cytoplasm. Finally, we discuss virion delivery into the nucleus, and consider the nature of transcriptionally silent DNA species that, escaping detection by the cell, might allow unhampered progress into S-phase and hence unleash the parvoviral Trojan horse.

Cytosine deamination and selection of CpG suppressed clones are the two major independent biological forces that shape codon usage bias in coronaviruses

Using the complete genome sequences of 19 coronavirus genomes, we analyzed the codon usage bias, dinucleotide relative abundance and cytosine deamination in coronavirus genomes. Of the eight codons that contain CpG, six were markedly suppressed. The mean NNU/NNC ratio of the six amino acids using either NNC or NNU as codon is 3.262, suggesting cytosine deamination. Among the 16 dinucleotides, CpG was most markedly suppressed (mean relative abundance 0.509). No correlation was observed between CpG abundance and mean NNU/NNC ratio. Among the 19 coronaviruses, CoV-HKU1 showed the most extreme codon usage bias and extremely high NNU/NNC ratio of 8.835. Cytosine deamination and selection of CpG suppressed clones by the immune system are the two major independent biochemical and biological selective forces that shape codon usage bias in coronavirus genomes. The underlying mechanism for the extreme codon usage bias, cytosine deamination and G + C content in CoV-HKU1 warrants further studies.

Phylogeny and evolution of papillomaviruses based on the E1 and E2 proteins

Papillomaviridae are a family of small double-stranded DNA viruses that infect stratified squamous epithelia in vertebrates. Members of this family are causative agents of malignant tumours, such as cervical cancer while others are associated with benign proliferative lesions. So far, Papillomaviruses (PVs) are classified according to the sequence identity in the capsid gene L1. However, evidence has accumulated indicating a discontinuity in the evolutionary history of the L1 and L2 genes of many PVs, giving rise to differences in the phylogenetic reconstructions of the early and of the late genes. Neither the oncogenes E5, E6 and E7 nor the upstream regulatory region are suitable for phylogenetic inference due to the poor conservation along the Papillomaviridae family. We have analysed here the evolutionary relationships of the PVs with respect to the E1 and E2 proteins, and the results provide both phylogeny and biologic behaviour of the viruses. The hierarchical taxonomic relationships can be structured as an alternative classification system in which mucosal high-risk viruses, mucosal low-risk viruses and viruses associated with cutaneous lesions are grouped separately and do not appear intermingled. Some important trends are also observed: first, evolution of the PVs has not been homogeneous, even in viruses that infect the same host, and second mucosal human PVs have evolved faster than their cutaneous counterparts. The evolutionary analysis based on the E1 and E2 proteins will allow us to better understand the generation of the diversity of the PVs and the development of malignancy associated with these viruses.