Contrasts in codon usage of latent versus productive genes of Epstein-Barr virus: data and hypotheses

Nucleotide and codon usage biases involved in the evolution of African swine fever virus: A comparative genomics analysis

Since African swine fever virus (ASFV) replication is closely related to its host's machinery, codon usage of viral genome can be subject to selection pressures. A better understanding of codon usage can give new insights into viral evolution. We implemented information entropy and revealed that the nucleotide usage pattern of ASFV is significantly associated with viral isolation factors (region and time), especially the usages of thymine and cytosine. Despite the domination of adenine and thymine in the viral genome, we found that mutation pressure alters the overall codon usage pattern of ASFV, followed by selective forces from natural selection. Moreover, the nucleotide skew index at the gene level indicates that nucleotide usages influencing synonymous codon bias of ASFV are significantly correlated with viral protein hydropathy. Finally, evolutionary plasticity is proved to contribute to the weakness in synonymous codons with A- or T-end serving as optimal codons of ASFV, suggesting that fine-tuning translation selection plays a role in synonymous codon usages of ASFV for adapting host. Taken together, ASFV is subject to evolutionary dynamics on nucleotide selections and synonymous codon usage, and our detailed analysis offers deeper insights into the genetic characteristics of this newly emerging virus around the world.

Enterovirus A Shows Unique Patterns of Codon Usage Bias in Conventional Versus Unconventional Clade

Enterovirus A (EV-A) species cause hand, foot and mouth disease (HFMD), threatening the health of young children. Understanding the mutual codon usage pattern of the virus and its host(s) has fundamental and applied values. Here, through examining multiple codon usage parameters, we found that the codon usage bias among EV-A strains varies and is clade-specific. EVA76, EVA89, EVA90, EVA91 and EVA92, the unconventional clade of EV-A strains, show unique codon usage pattern relative to the two conventional clades, including EVA71, CVA16, CVA6 and CVA10, etc. Analyses of Effective Number of Codon (ENC), Correspondence Analysis (COA) and Parity Rule 2 (PR2), etc., revealed that the codon usage patterns of EV-A strains are shaped by mutation pressure and natural selection. Based on the neutrality analysis, we determined the dominant role of natural selection in the formation of the codon usage bias of EV-A. In addition, we have determined the codon usage compatibility of potential hosts for EV-A strains using codon adaptation index (CAI), relative codon deoptimization index (RCDI) and similarity index (SiD) analyses, and found that EV-A showed host-specific codon adaptation patterns in different clades. Finally, we confirmed that the unique codon usage pattern of the unconventional clade affected protein expression level in human cell lines. In conclusion, we identified novel characteristics of codon usage bias in distinct EV-A clades associated with their host range, transmission and pathogenicity.

Composition, codon usage pattern, protein properties, and influencing factors in the genomes of members of the family Anelloviridae

The present study was carried out on 62 genome sequences of members of the family Anelloviridae, as there have been no reports of genome analysis of these DNA viruses using a bioinformatics approach. The genes were found to be rich in AC content with low codon usage bias (CUB). Relative synonymous codon usage (RSCU) values identified the preferred codons for each amino acid in the family. The codon AGA was overrepresented, while the codons TCG, TTG, CGG, CGT, ACG, GCG and GAT were underrepresented in all of the genomes. A significant correlation was found between the effective number of codons (ENC) and base constraints, indicating that compositional properties might have influenced the CUB. A highly significant correlation was observed between the overall base content and the base content at the third codon position, indicating that mutations might have affected the CUB. A highly significant positive correlation was observed between GC12 and GC3 (r = 0.904, p < 0.01), which indicated that directional mutation pressure influenced all three codon positions. A neutrality plot revealed that the contribution of mutation and natural selection in determining the CUB was 58.6% and 41.4%, respectively.

Origin and Evolution of H1N1/pdm2009: A Codon Usage Perspective

The H1N1/pdm2009 virus is a new triple-reassortant virus. While Eurasian avian-like and triple-reassortant swine influenza viruses are the direct ancestors of H1N1/pdm2009, the classic swine influenza virus facilitate the spectrum of influenza A diversity in pig population when the reassortant events occurred during 1998 to April 2009. The factors that facilitate the final formation of this gene constellation for H1N1/pdm2009 virus from this complex gene pool remain unknown. Since a novel successful virus should efficiently replicate and transmit in their hosts, in this study, we estimated the adaptability of the codon usage patterns of the pool of genes from these lineages of swine influenza viruses to the human expression system. We found that the MP and NA genes of Eurasian avian-like swine influenza viruses, and the PB2, PB1 and PA genes of triple-reassortant swine influenza viruses were best adapted to the human codon usage pattern. As these genes participated in the development of H1N1/pdm2009, they might help in viral replication and strengthen its competitiveness during its emergence. After its emergence in the human population, a gradual optimization of codon usage patterns between 2009 and 2019 to the human codon usage for the H1N1/pdm2009 genes was detected. This reveals that ongoing adaptive evolution, after its original incursion, occurred to further increase the adaptability of overall gene cassette to human expression system.

Codon usage pattern and its influencing factors in different genomes of hepadnaviruses

Codon usage bias (CUB) arises from the preference for a codon over codons for the same amino acid. The major factors contributing to CUB are evolutionary forces, compositional properties, gene expression, and protein properties. The present analysis was performed to investigate the compositional properties and the extent of CUB across the genomes of members of the family Hepadnaviridae, as previously no work using bioinformatic tools has been reported. The viral genes were found to be AT rich with low CUB. Analysis of relative synonymous codon usage (RSCU) was used to identify overrepresented and underrepresented codons for each amino acid. Correlation analysis of overall nucleotide composition and its composition at the third codon position suggested that mutation pressure might influence the CUB. A highly significant correlation was observed between GC12 and GC3 (r = 0.910, p < 0.01), indicating that directional mutation affected all three codon positions across the genome. Translational selection (P2) and mutational responsive index (MRI) values of genes suggested that mutation plays a more important role than translational selection in members of the family Hepadnaviridae.

Canine parvovirus type 2 (CPV-2) and Feline panleukopenia virus (FPV) codon bias analysis reveals a progressive adaptation to the new niche after the host jump

Based on virus dependence from host cell machinery, their codon usage is expected to show a strong relation with the host one. Even if this association has been stated, especially for bacteria viruses, the linkage is considered to be less consistent for more complex organisms and a codon bias adaptation after host jump has never been proven. Canine parvovirus type 2 (CPV-2) was selected as a model because it represents a well characterized case of host jump, originating from Feline panleukopenia virus (FPV). The current study demonstrates that the adaptation to specific tissue and host codon bias affected CPV-2 evolution. Remarkably, FPV and CPV-2 showed a higher closeness toward the codon bias of the tissues they display the higher tropism for. Moreover, after the host jump, a clear and significant trend was evidenced toward a reduction in the distance between CPV-2 and the dog codon bias over time. This evidence was not confirmed for FPV, suggesting that an equilibrium has been reached during the prolonged virus-host co-evolution. Additionally, the presence of an intermediate pattern displayed by some strains infecting wild species suggests that these could have facilitated the host switch also by acting on codon bias.

Codon Optimization Leads to Functional Impairment of RD114-TR Envelope Glycoprotein

Lentiviral vectors (LVs) are a highly valuable tool for gene transfer currently exploited in basic, applied, and clinical studies. Their optimization is therefore very important for the field of vectorology and gene therapy. A key molecule for LV function is the envelope because it guides cell entry. The most commonly used in transiently produced LVs is the vesicular stomatitis virus glycoprotein (VSV-G) envelope, whose continuous expression is, however, toxic for stable LV producer cells. In contrast, the feline endogenous retroviral RD114-TR envelope is suitable for stable LV manufacturing, being well tolerated by producer cells under constitutive expression. We have previously reported successful, transient and stable production of LVs pseudotyped with RD114-TR for good transduction of T lymphocytes and CD34⁺ cells. To further improve RD114-TR-pseudotyped LV cell entry by increasing envelope expression, we codon-optimized the RD114-TR open reading frame (ORF). Here we show that, despite the RD114-TRco precursor being produced at a higher level than the wild-type counterpart, it is unexpectedly not duly glycosylated, exported to the cytosol, and processed. Correct cleavage of the precursor in the functional surface and transmembrane subunits is prevented in vivo, and, consequently, the unprocessed precursor is incorporated into LVs, making them inactive.

Testing hypotheses for the presence of tRNA genes in mycobacteriophage genomes

The presence of tRNA genes in bacteriophages has been explained on the basis of codon usage (tRNA genes are retained in the phage genome if they correspond to codons more common in the phage than in its host) or amino acid usage (independent of codon, the amino acid corresponding to the retained tRNA gene is more common in the phage genome than in the bacterial host). The existence of a large database of sequenced mycobacteriophages, isolated on the common host Mycobacterium smegmatis, allows us to test the above hypotheses as well as explore other hypotheses for the presence of tRNA genes. Our analyses suggest that amino acid rather than codon usage better explains the presence of tRNA genes in mycobacteriophages. However, closely related phages that differ in the presence of tRNA genes in their genomes are capable of lysing the common bacterial host and do not differ in codon or amino acid usage. This suggests that the benefits of having tRNA genes may be associated with either growth in the host or the ability to infect more hosts (i.e., host range) rather than simply infecting a particular host.

Cancer, Warts, or Asymptomatic Infections: Clinical Presentation Matches Codon Usage Preferences in Human Papillomaviruses

Viruses rely completely on the hosts’ machinery for translation of viral transcripts. However, for most viruses infecting humans, codon usage preferences (CUPrefs) do not match those of the host. Human papillomaviruses (HPVs) are a showcase to tackle this paradox: they present a large genotypic diversity and a broad range of phenotypic presentations, from asymptomatic infections to productive lesions and cancer. By applying phylogenetic inference and dimensionality reduction methods, we demonstrate first that genes in HPVs are poorly adapted to the average human CUPrefs, the only exception being capsid genes in viruses causing productive lesions. Phylogenetic relationships between HPVs explained only a small proportion of CUPrefs variation. Instead, the most important explanatory factor for viral CUPrefs was infection phenotype, as orthologous genes in viruses with similar clinical presentation displayed similar CUPrefs. Moreover, viral genes with similar spatiotemporal expression patterns also showed similar CUPrefs. Our results suggest that CUPrefs in HPVs reflect either variations in the mutation bias or differential selection pressures depending on the clinical presentation and expression timing. We propose that poor viral CUPrefs may be central to a trade-off between strong viral gene expression and the potential for eliciting protective immune response.

Downregulating viral gene expression: codon usage bias manipulation for the generation of novel influenza A virus vaccines

Vaccination represents the best option to protect humans against influenza virus. However, improving the effectiveness of current vaccines could better stifle the health burden caused by viral infection. Protein synthesis from individual genes can be downregulated by synthetically deoptimizing a gene's codon usage. With more rapid and affordable nucleotide synthesis, generating viruses that contain genes with deoptimized codons is now feasible. Attenuated, vaccine-candidate viruses can thus be engineered with hitherto uncharacterized properties. With eight gene segments, influenza A viruses with variably recoded genomes can produce a spectrum of attenuation that is contingent on the gene segment targeted and the number of codon changes. This review summarizes different targets and approaches to deoptimize influenza A virus codons for novel vaccine generation.

System analysis of synonymous codon usage biases in archaeal virus genomes

Recent studies of geothermally heated aquatic ecosystems have found widely divergent viruses with unusual morphotypes. Archaeal viruses isolated from these hot habitats usually have double-stranded DNA genomes, linear or circular, and can infect members of the Archaea domain. In this study, the synonymous codon usage bias (SCUB) and dinucleotide composition in the available complete archaeal virus genome sequences have been investigated. It was found that there is a significant variation in SCUB among different Archaeal virus species, which is mainly determined by the base composition. The outcome of correspondence analysis (COA) and Spearman׳s rank correlation analysis shows that codon usage of selected archaeal virus genes depends mainly on GC richness of genome, and the gene׳s function, albeit with smaller effects, also contributes to codon usage in this virus. Furthermore, this investigation reveals that aromaticity of each protein is also critical in affecting SCUB of these viral genes although it was less important than that of the mutational bias. Especially, mutational pressure may influence SCUB in SIRV1, SIRV2, ARV1, AFV1, and PhiCh1 viruses, whereas translational selection could play a leading role in HRPV1׳s SCUB. These conclusions not only can offer an insight into the codon usage biases of archaeal virus and subsequently the possible relationship between archaeal viruses and their host, but also may help in understanding the evolution of archaeal viruses and their gene classification, and more helpful to explore the origin of life and the evolution of biology.

•

The SCUB of archaeal virus genes depends mainly on GC richness of genome.

•

The mutational pressure is the main factor that influences SCUB.

•

The aromaticity of each protein is also critical in affecting SCUB.

•

The translational selection could play a leading role in HRPV1׳s SCUB.

•

The mode is helpful to explore the origin of life and the evolution of biology.

Hepatitis A virus adaptation to cellular shutoff is driven by dynamic adjustments of codon usage and results in the selection of populations with altered capsids

Hepatitis A virus (HAV) has a highly biased and deoptimized codon usage compared to the host cell and fails to inhibit host protein synthesis. It has been proposed that an optimal combination of abundant and rare codons controls the translation speed required for the correct capsid folding. The artificial shutoff host protein synthesis results in the selection of variants containing mutations in the HAV capsid coding region critical for folding, stability, and function. Here, we show that these capsid mutations resulted in changes in their antigenicity; in a reduced stability to high temperature, low pH, and biliary salts; and in an increased efficacy of cell entry. In conclusion, the adaptation to cellular shutoff resulted in the selection of large-plaque-producing virus populations.

IMPORTANCE

HAV has a naturally deoptimized codon usage with respect to that of its cell host and is unable to shut down the cellular translation. This fact contributes to the low replication rate of the virus, in addition to other factors such as the highly inefficient internal ribosome entry site (IRES), and explains the outstanding physical stability of this pathogen in the environment mediated by a folding-dependent highly cohesive capsid. Adaptation to artificially induced cellular transcription shutoff resulted in a redeoptimization of its capsid codon usage, instead of an optimization. These genomic changes are related to an overall change of capsid folding, which in turn induces changes in the cell entry process. Remarkably, the adaptation to cellular shutoff allowed the virus to significantly increase its RNA uncoating efficiency, resulting in the selection of large-plaque-producing populations. However, these populations produced much-debilitated virions.

Role of genomic and proteomic tools in the study of host-virus interactions and virus evolution

Viruses have short replication cycles and produce genomic variants within a host, a process that seems to adapt to their specific host and also enable them to infect new hosts. The recent emergence of viral genomic variants from the circulating pool within the host population and re-emergence of the old ones are posing serious threat to agriculture, animal husbandry and humanity as a whole. This review assesses the potential role of genomic and proteomic tools that can monitor not only the course of infection and pathogenesis, but also predict the pandemic or zoonotic epidemic potential of a virus in a previously exposed or immunologically naive biological population.

Sensitive detection of viral transcripts in human tumor transcriptomes

In excess of % of human cancer incidents have a viral cofactor. Epidemiological studies of idiopathic human cancers indicate that additional tumor viruses remain to be discovered. Recent advances in sequencing technology have enabled systematic screenings of human tumor transcriptomes for viral transcripts. However, technical problems such as low abundances of viral transcripts in large volumes of sequencing data, viral sequence divergence, and homology between viral and human factors significantly confound identification of tumor viruses. We have developed a novel computational approach for detecting viral transcripts in human cancers that takes the aforementioned confounding factors into account and is applicable to a wide variety of viruses and tumors. We apply the approach to conducting the first systematic search for viruses in neuroblastoma, the most common cancer in infancy. The diverse clinical progression of this disease as well as related epidemiological and virological findings are highly suggestive of a pathogenic cofactor. However, a viral etiology of neuroblastoma is currently contested. We mapped transcriptomes of neuroblastoma as well as positive and negative controls to the human and all known viral genomes in order to detect both known and unknown viruses. Analysis of controls, comparisons with related methods, and statistical estimates demonstrate the high sensitivity of our approach. Detailed investigation of putative viral transcripts within neuroblastoma samples did not provide evidence for the existence of any known human viruses. Likewise, de-novo assembly and analysis of chimeric transcripts did not result in expression signatures associated with novel human pathogens. While confounding factors such as sample dilution or viral clearance in progressed tumors may mask viral cofactors in the data, in principle, this is rendered less likely by the high sensitivity of our approach and the number of biological replicates analyzed. Therefore, our results suggest that frequent viral cofactors of metastatic neuroblastoma are unlikely.

Many human cancers are caused by infections with tumor viruses and identification of these pathogens is considered a critical contribution to cancer prevention. Deep sequencing enables us to systematically investigate viral nucleotide signatures in order to either verify or exclude the existence of viruses in idiopathic human cancers. We have developed Virana, a novel computational approach for identifying tumor viruses in human cancers that is applicable to a wide variety of tumors and viruses. Virana firstly addresses several important biological confounding factors that may hinder successful detection of these pathogens. We applied our approach in the first systematic search for cancer-causing viruses in metastatic neuroblastoma, the most common form of cancer in infancy. Although the heterogeneous clinical progression of this disease as well as epidemiological and virological findings are suggestive of a pathogenic cofactor, the viral etiology of neuroblastoma is currently contested. We conducted an analysis of experimental controls, comparisons with related approaches, as well as statistical analyses in order to validate our method. In spite of the high sensitivity of our approach, analyses of neuroblastoma transcriptomes did not provide evidence for the existence of any known or unknown human viruses. Our results therefore suggest that frequent viral cofactors of metastatic neuroblastoma are unlikely.

Translation Elongation Rate Measurement of Epstein-Barr Virus Strain GD1

BACKGROUND

Epstein-Barr Virus (EBV) has a great co relationship with human malignancies such as gastric carcinoma. Synonymous codon investigations in viruses could help designing vaccine, to generate immunity. Codon Adaptation Index (CAI) has measured translation elongation rate, among the highly expressed genes. The aim of this study was: usage of "CAI" to measure translation efficiency to know how fast EBV-GD1 could produce its proteins.

METHODS

The complete genomic sequences of human herpes virus 4 strain GD1 have retrieved from <http://www.ncbi.nlm.nih.gov/sites/gquery> (GenBank accession no. AY961628) to extract all protein-coding genes. The sequences have analyzed with DAMBE software.

RESULTS

The results have shown that CAI values for the EBV-GD1 genes were 0.76356 ± 0.02957. The highest and lowest CAI values were 0.82233 and 0.68321 respectively. The results have shown that highly expressed genes mostly had more codon usage bias than low expressed genes.

CONCLUSION

The results provide and introduce not only a system, but also the principles in order to understand the pathogenesis and evolution of EBV-GD1, to open a window, in order to make a better product or vaccine to challenge with the virus.

High codon adaptation in citrus tristeza virus to its citrus host

Background

Citrus tristeza virus (CTV), a member of the genus Closterovirus within the family Closteroviridae, is the causal agent of citrus tristeza disease. Previous studies revealed that the negative selection, RNA recombination and gene flow were the most important forces that drove CTV evolution. However, the CTV codon usage was not studied and thus its role in CTV evolution remains unknown.

Results

A detailed comparative analysis of CTV codon usage pattern was done in this study. Results of the study show that although in general CTV does not have a high degree of codon usage bias, the codon usage of CTV has a high level of resemblance to its host codon usage. In addition, our data indicate that the codon usage resemblance is only observed for the woody plant-infecting closteroviruses but not the closteroviruses infecting the herbaceous host plants, suggesting the existence of different virus-host interactions between the herbaceous plant-infecting and woody plant-infecting closteroviruses.

Conclusion

Based on the results, we suggest that in addition to RNA recombination, negative selection and gene flow, host plant codon usage selection can also affect CTV evolution.

Level of gene expression is a major determinant of protein evolution in the viral order Mononegavirales

Although the rate at which proteins change is a key parameter in molecular evolution, its determinants are poorly understood in viruses. A variety of factors, including gene length, codon usage bias, protein abundance, protein function, and gene expression level, have been shown to affect the rate of protein evolution in a diverse array of organisms. However, the role of these factors in viral evolution has yet to be addressed. The polar 3'-5' stepwise attenuation of transcription in the Mononegavirales, a group of single-strand negative-sense RNA viruses, provides a unique system to explore the determinants of protein evolution in viruses. We analyzed the relative importance of a variety of factors in shaping patterns of sequence variation in full-length genomes from 13 Mononegavirales species. Our analysis suggests that the level of gene expression, and by extension the relative genomic position of each gene, is a key determinant of the protein evolution in these viruses. This appears to be the consequence of selection for translational robustness, but not for translational accuracy, in highly expressed genes. The small genome size and number of proteins encoded by these viruses allowed us to identify other protein-specific factors that may also play a role in virus evolution, such as host-virus interactions and functional constraints. Finally, we explored the evolutionary pressures acting on noncoding regions in Mononegavirales genomes and observed that, despite being less constrained than coding regions, their evolutionary rates are also associated with genomic position.

Characterization of Synonymous Codon Usage in the Newly Identified Duck Plague Virus UL16 Gene

A comparative analysis of the codon usage bias in the newly identified UL16 gene(GenBank accession no.EU195095) of DPV and the UL16 gene of 22 reference herpesviruses was performed. In this study, the synonymous codon usage bias of UL16 gene in the 23 herpesviruses have been analyzed and the results showed obvious differences by the CAI, RSCU, ENC and GC_3s. The results revealed that the synonymous codons with A and T at the third codon positon have widely usage in the codon of UL16 gene of DPV. The ENC-GC_3s plot revealed that the genetic heterogeneity in UL16 gene of herpesviruses was constrained by G+C content at the third codon position. The phylogenetic analysis suggested that DPV was evolutionarily closer to herpesviruses which further clustered into Alphaherpesvirinae. Furthermore the ORF of DPV UL16 gene has sequential rare codons. There were 21 codons showing distinct usage differences between DPV with Escherichia coli, 19 codons showing distinct usage differences between DPV with yeast, and 20 between DPV and Human. Therefore the Escherichia coli, Yeast and Human expression system were suitable for the expression of DPV UL16 gene if some codons could be optimized.

The characteristics of the synonymous codon usage in hepatitis B virus and the effects of host on the virus in codon usage pattern

BACKGROUND

Hepatitis B virus (HBV) infection is one of the main human health problem and causes a large-scale of patients chronic infection worldwide.. As the replication of HBV depends on its host cell system, codon usage pattern for the viral gene might be susceptible to two main selections, namely mutation pressure and translation selection. In this case, a deeper investigation between HBV evolution and host adaptive response might assist control this disease.

RESULT

Relative synonymous codon usage (RSCU) values for the whole HBV coding sequence were studied by Principal component analysis (PCA). The characteristics of the synonymous codon usage patterns, nucleotide contents and the comparison between ENC values of the whole HBV coding sequence indicated that the interaction between virus mutation pressure and host translation selection exists in the processes of HBV evolution. The synonymous codon usage pattern of HBV is a mixture of coincidence and antagonism to that of host cell. But the difference of genetic characteristic of HBV failed to be observed to its different epidemic areas or subtypes, suggesting that geographic factor is limited to influence the evolution of this virus, while genetic characteristic based on HBV genotypes could be divided into three groups, namely (i) genotyps A and E, (ii) genotype B, (iii) genotypes C, D and G.

CONCLUSION

Codon usage patterns from PCA for identification of evolutionary trends in HBV provide an alternative approach to understand the evolution of HBV. Further more, a combined selection of mutation pressure with translation selection on codon usage might shed a light on understanding the evolutionary trends of HBV genotypes.

Bioinformatics analysis of rabbit haemorrhagic disease virus genome

Background

Rabbit haemorrhagic disease virus (RHDV), as the pathogeny of Rabbit haemorrhagic disease, can cause a highly infectious and often fatal disease only affecting wild and domestic rabbits. Recent researches revealed that it, as one number of the Caliciviridae, has some specialties in its genome, its reproduction and so on.

Results

In this report, we firstly analyzed its genome and two open reading frameworks (ORFs) from this aspect of codon usage bias. Our researches indicated that mutation pressure rather than natural is the most important determinant in RHDV with high codon bias, and the codon usage bias is nearly contrary between ORF1 and ORF2, which is maybe one of factors regulating the expression of VP60 (encoding by ORF1) and VP10 (encoding by ORF2). Furthermore, negative selective constraints on the RHDV whole genome implied that VP10 played an important role in RHDV lifecycle.

Conclusions

We conjectured that VP10 might be beneficial for the replication, release or both of virus by inducing infected cell apoptosis initiate by RHDV. According to the results of the principal component analysis for ORF2 of RSCU, we firstly separated 30 RHDV into two genotypes, and the ENC values indicated ORF1 and ORF2 were independent among the evolution of RHDV.

Codon usage roles in human papillomavirus

Human papillomavirus (HPV) genomes, similar to other virus genomes, frequently have a G + C content significantly different from their host species. The HPV genomes show a strong codon usage bias to 18 codons, with 14 showing T at the third position amongst degenerately encoded amino acids. The codon usage pattern in HPV genome plays an important role, which regulates low or non-translational expression of the viral capsid genes and results in very weak protein expression of oncogenes in a wide range of mammalian cells. Codon modification has been proved to be a powerful technology to overcome the translational blockage and weak expression of both HPV capsid genes and oncogenes in different expression systems. Furthermore, keratinocytes are the host cells of HPV infection; the codon usage in HPV capsid genes matches available aminoacyl-tRNAs in differentiated keratinocytes to modulate their protein expression. HPV DNA vaccines with codon optimization have been shown to have higher immunogenicity and induce both strong cellular and humoral responses in animal models, which may be a promising form of therapeutic HPV vaccines.

Comparative analysis of ovine adenovirus 287 and human adenovirus 2 and 5 based on their codon usage

Ovine adenovirus 287 (OAdV287) emerges as one of the most promising gene vectors resulting from its unique biological characteristics. To obtain a more detailed knowledge about the codon usage of OAdV287, a comparative study based on the codon usage of OAdV287 and the prototypes of human adenovirus serotypes 2 and 5 (HAdV2/5) was carried out. Some commonly used indices measuring the codon usage patterns, including effective number of codons, relative synonymous codon usage, and statistical methods, were adopted. Overall, OAdV287 had a more biased and conservative codon usage pattern than that of HAdV2/5. Both mutation pressure and natural selection played important roles in shaping the codon usage patterns of these three adenoviruses. All the preference codons of OAdV287 had A/U ends and were totally different from those of sheep and humans; however, the preference codons of HAdV2/5 mostly had G/C ends and were mostly coincident with those of sheep and humans. The codon usage analysis in this study supplies some clues for further comprehending the unique biological characteristics of OAdV287 as gene vectors.

Analysis of codon usage in type 1 and the new genotypes of duck hepatitis virus

In this study, an abundant (A + U)% and low codon bias were revealed in duck hepatitis virus type 1 (DHV-1) and the new serotype strains isolated from Taiwan, South Korea and Mainland China (DHV-N). 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 these samples. By comparative analysis of the codon usage patterns of 40 ORFs of DHV, we found that all of DHV-1 strains grouped in genotype C; the DHV-N strains isolated in South Korea and China clustered into genotypes B; and the DHV-N strains isolated from Taiwan clustered into genotypes A. The findings revealed that more than one subtype of DHV-1 circulated in East Asia. Furthermore, the results of phylogenetic analyses based on RSCU values and Clustal W method indicated obvious phylogenetic congruities. This suggested that better genome consistency of DHV may exist in nature and phylogenetic analyses based on RSCU values maybe a good method in classifying genotypes of the virus. Our work might give some clues to the features and some evolutionary information of DHV.

The characteristics of the synonymous codon usage in enterovirus 71 virus and the effects of host on the virus in codon usage pattern

To give a new perspective on the evolutionary characteristics shaping the genetic diversity of enterovirus 71 (EV71) and the effects of natural selection from its host on the codon usage pattern of the virus, the relative synonymous codon usage (RSCU) values, codon usage bias (CUB) values, effective number of codons (ENCs) values and nucleotide contents were calculated to implement a comparative analysis to evaluate the dynamics of the virus evolution. The characteristics of the synonymous codon usage patterns and nucleotide contents of EV71 and the comparison between ENC values for the whole coding sequence of EV71 and that of coding sequences for viral proteins of EV71 all indicate that the interaction between mutation pressure from virus and natural selection from host exists in the processes of evolution of EV71. The synonymous codon usage pattern of EV71 is a mixture of coincidence and antagonism to that of host cell. In addition, the genetic diversity of EV71 strains and the preferential selection of some synonymous codons in EV71 strains based on the different epidemic areas were observed, suggesting that geographic and social factors may play roles in influencing the evolution of this virus.

Codon usage bias and the evolution of influenza A viruses. Codon Usage Biases of Influenza Virus

BACKGROUND

The influenza A virus is an important infectious cause of morbidity and mortality in humans and was responsible for 3 pandemics in the 20th century. As the replication of the influenza virus is based on its host's machinery, codon usage of its viral genes might be subject to host selection pressures, especially after interspecies transmission. A better understanding of viral evolution and host adaptive responses might help control this disease.

RESULTS

Relative Synonymous Codon Usage (RSCU) values of the genes from segment 1 to segment 6 of avian and human influenza viruses, including pandemic H1N1, were studied via Correspondence Analysis (CA). The codon usage patterns of seasonal human influenza viruses were distinct among their subtypes and different from those of avian viruses. Newly isolated viruses could be added to the CA results, creating a tool to investigate the host origin and evolution of viral genes. It was found that the 1918 pandemic H1N1 virus contained genes with mammalian-like viral codon usage patterns, indicating that the introduction of this virus to humans was not through in toto transfer of an avian influenza virus.Many human viral genes had directional changes in codon usage over time of viral isolation, indicating the effect of host selection pressures. These changes reduced the overall GC content and the usage of G at the third codon position in the viral genome. Limited evidence of translational selection pressure was found in a few viral genes.

CONCLUSIONS

Codon usage patterns from CA allowed identification of host origin and evolutionary trends in influenza viruses, providing an alternative method and a tool to understand the evolution of influenza viruses. Human influenza viruses are subject to selection pressure on codon usage which might assist in understanding the characteristics of newly emerging viruses.

RCDI/eRCDI: a web-server to estimate codon usage deoptimization

Background

The Relative Codon Deoptimization Index (RCDI) was developed by Mueller et al. (2006) as measure of codon deoptimization by comparing how similar is the codon usage of a gene and the codon usage of a reference genome.

Findings

RCDI/eRCDI is a web application server that calculates the Relative Codon Deoptimization Index and a new expected value for the RCDI (eRCDI). The RCDI is used to estimate the similarity of the codon frequencies of a specific gene in comparison to a given reference genome. The eRCDI is determined by generating random sequences with similar G+C and amino acid composition to the input sequences and may be used as an indicator of the significance of the RCDI values. RCDI/eRCDI is freely available at http://genomes.urv.cat/CAIcal/RCDI.

Conclusions

This web server will be a useful tool for genome analysis, to understand host-virus phylogenetic relationships or to infer the potential host range of a virus and its replication strategy, as well as in experimental virology to ease the step of gene design for heterologous protein expression.

Fine-tuning translation kinetics selection as the driving force of codon usage bias in the hepatitis A virus capsid

Hepatitis A virus (HAV), the prototype of genus Hepatovirus, has several unique biological characteristics that distinguish it from other members of the Picornaviridae family. Among these, the need for an intact eIF4G factor for the initiation of translation results in an inability to shut down host protein synthesis by a mechanism similar to that of other picornaviruses. Consequently, HAV must inefficiently compete for the cellular translational machinery and this may explain its poor growth in cell culture. In this context of virus/cell competition, HAV has strategically adopted a naturally highly deoptimized codon usage with respect to that of its cellular host. With the aim to optimize its codon usage the virus was adapted to propagate in cells with impaired protein synthesis, in order to make tRNA pools more available for the virus. A significant loss of fitness was the immediate response to the adaptation process that was, however, later on recovered and more associated to a re-deoptimization rather than to an optimization of the codon usage specifically in the capsid coding region. These results exclude translation selection and instead suggest fine-tuning translation kinetics selection as the underlying mechanism of the codon usage bias in this specific genome region. Additionally, the results provide clear evidence of the Red Queen dynamics of evolution since the virus has very much evolved to re-adapt its codon usage to the environmental cellular changing conditions in order to recover the original fitness.

Each organism has a specific codon usage signature. Translational selection i.e., selection for the codon adaptation to the tRNA pools, is one of the driving forces of codon bias. In the virus world, this implies an adjustment of the virus codon usage to that of the host cell. Hepatitis A virus appears as an exception to the rule, with a highly deoptimized codon usage, suggesting that translational selection is not the underlying mechanism of its codon bias. However, since the virus lacks a mechanism of cellular protein synthesis inhibition, the deoptimized codon usage may be envisaged as a hawk (cell) and dove (hepatitis A virus) competition strategy for tRNAs and translational selection as well. To confirm this possibility, we artificially induced cell protein synthesis shut-off, thus increasing the tRNA pool availability for the virus, and we took advantage of the quasispecies dynamics to elucidate changes in its codon usage. Virus adaptation to the drug results in a re-deoptimization of codon usage in the capsid region, suggesting a requirement of a slow translation rate, i.e., a translation kinetic selection, instead of a translational selection associated with an optimization of the codon usage. Translation kinetics control is based on the right combination of codons (common and rare) that allows a regulated ribosome traffic rate ensuring the proper protein folding. Capsid folding is critical for a virus transmitted through the fecal-oral route with long extracorporeal periods.

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.

Analysis of synonymous codon usage in the UL24 gene of duck enteritis virus

The analysis on codon usage bias of UL24 gene of duck enteritis virus (DEV) may improve our understanding of the evolution and pathogenesis of DEV and provide a basis for understanding the relevant mechanism for biased usage of synonymous codons and for selecting appropriate expression systems to improve the expression of target genes. The codon usage bias of UL24 genes of DEV and 27 reference herpesviruses were analyzed. The results showed that codon of UL24 gene of DEV was strong bias toward the synonymous codons with A and T at the third codon position. A high level of diversity in codon usage bias existed, and the effective number of codons used in a gene plot revealed that the genetic heterogeneity in UL24 gene of herpesviruses was constrained by the G + C content. The phylogentic analysis suggested that DEV was evolutionarily closer to Alphaherpesvirinae and that there was no significant deviation in codon usage in different virus strains. There were 20 codons showing distinct usage differences between DEV and Escherichia coli, 23 between DEV and Homo sapiens, but only 16 codons between DEV and yeast. Therefore the yeast expression system may be more suitable for the expression of DEV genes.

Convergent host-parasite codon usage between honeybee and bee associated viral genomes

By correlating the codon usage in four insects (the honeybee, red flour beetle, mosquito and fruit fly) with six honeybee host specific viruses, we found that the codon usage patterns of the bee viruses were strongly related to that of the honeybee and only weakly related to the red flour beetle. The insects shared varying degrees of codon usage similarity which roughly follow the known phylogenetic relatedness. All of the codon usage similarity can be described by relatedness-by-descent except for the high codon usage similarity between the honeybee and honeybee associated viruses. This evidence for the convergent evolution of the honeybee viruses toward the codon usage of the honeybee suggests that small host specific viral genomes have the freedom to quickly optimize codon usage to successfully parasitize their preferred host. The codon usage co-evolution of the six host specific honeybee viruses towards the codon usage of the honeybee described in this paper is the first evidence for codon usage correlation between an insect host and a single stranded RNA virus.

Horizontal gene transfer and nucleotide compositional anomaly in large DNA viruses

Background

DNA viruses have a wide range of genome sizes (5 kb up to 1.2 Mb, compared to 0.16 Mb to 1.5 Mb for obligate parasitic bacteria) that do not correlate with their virulence or the taxonomic distribution of their hosts. The reasons for such large variation are unclear. According to the traditional view of viruses as gifted "gene pickpockets", large viral genome sizes could originate from numerous gene acquisitions from their hosts. We investigated this hypothesis by studying 67 large DNA viruses with genome sizes larger than 150 kb, including the recently characterized giant mimivirus. Given that horizontally transferred DNA often have anomalous nucleotide compositions differing from the rest of the genome, we conducted a detailed analysis of the inter- and intra-genome compositional properties of these viruses. We then interpreted their compositional heterogeneity in terms of possible causes, including strand asymmetry, gene function/expression, and horizontal transfer.

Results

We first show that the global nucleotide composition and nucleotide word usage of viral genomes are species-specific and distinct from those of their hosts. Next, we identified compositionally anomalous (cA) genes in viral genomes, using a method based on Bayesian inference. The proportion of cA genes is highly variable across viruses and does not exhibit a significant correlation with genome size. The vast majority of the cA genes were of unknown function, lacking homologs in the databases. For genes with known homologs, we found a substantial enrichment of cA genes in specific functional classes for some of the viruses. No significant association was found between cA genes and compositional strand asymmetry. A possible exogenous origin for a small fraction of the cA genes could be confirmed by phylogenetic reconstruction.

Conclusion

At odds with the traditional dogma, our results argue against frequent genetic transfers to large DNA viruses from their modern hosts. The large genome sizes of these viruses are not simply explained by an increased propensity to acquire foreign genes. This study also confirms that the anomalous nucleotide compositions of the cA genes is sometimes linked to particular biological functions or expression patterns, possibly leading to an overestimation of recent horizontal gene transfers.

Mutation pressure shapes codon usage in the GC-Rich genome of foot-and-mouth disease virus

Foot-and-mouth disease (FMD) is economically the most important viral-induced livestock disease worldwide. In this study, we report the results of a survey of codon usage bias of FMD virus (FMDV) representing all seven serotypes (A, O, C, Asia 1, SAT 1, SAT 2, and SAT 3). Correspondence analysis, a commonly used multivariate statistical approach, was carried out to analyze synonymous codon usage bias. The analysis showed that the overall extent of codon usage bias in FMDV is low. Furthermore, the good correlation between the frequency of G + C at the synonymous third position of sense codons (GC3S) content at silent sites of each sequence and codon usage bias suggested that mutation pressure rather than natural (translational) selection is the most important determinant of the codon bias observed. In addition, other factors, such as the lengths of open reading frame (ORF) and the hydrophobicity of genes also influence the codon usage variation among the genomes of FMDV in a minor way. The result of phylogenetic analyses based on the relative synonymous codon usage (RSCU) values indicated a few obvious phylogenetic incongruities, which suggest that more FMDV genome diversity may exist in nature than is currently indicated. Our work might give some clues to the features of FMDV genome and some evolutionary information of this virus.

Generalized substitution of isoencoding codons shortens the duration of papillomavirus L1 protein expression in transiently gene-transfected keratinocytes due to cell differentiation

Recently we reported that gene codon composition determines differentiation-dependent expression of the PV L1 genes in mouse primary keratinocytes (KCs) in vitro and in vivo (Zhao et al. 2005, Mol. Cell Biol. 25:8643–8655). Here, we investigated whether generalized substitution of isoencoding codons affects the duration of expression of PV L1 genes in mouse and human KCs in day 1 culture transiently transfected with native (Nat) and codon modified (Mod) L1 genes. Following transient transfection, KC continuously transcribed both Nat and Mod PV L1 genes for at least 12 days, with the levels of L1 mRNAs from the Mod L1 genes significantly higher than those from the Nat L1 genes. However, continuous L1 protein expression at day 9 post-transfection was observed for both mouse and human KCs transfected with the Nat L1 genes only. Further, aa-tRNAs prepared from D8 KC cultures enhanced translation of two PV Nat L1 DNAs in RRL lysate and PV Nat L1 mRNAs in D0 cell-free lysate, whereas aa-tRNAs from D0 KCs enhanced translation of PV Mod L1 mRNAs in D8 cell-free lysate. It appears that aa-tRNAs in less-differentiated and differentiated KCs differentially match the PV Nat and Mod L1 mRNAs to regulate their translations in vitro.

A DNA vaccine encoding a codon-optimized human papillomavirus type 16 E6 gene enhances CTL response and anti-tumor activity

The HPV oncoproteins E6 and E7 are consistently expressed in HPV-associated cancer cells and are responsible for their malignant transformation. Therefore, HPV E6 and E7 are ideal target antigens for developing vaccines and immunotherapeutic strategies against HPV-associated neoplasms. Recently, it has been demonstrated that codon optimization of the HPV-16 E7 gene resulted in highly efficient translation of E7 and increased the immunogenicity of E7-specific DNA vaccines. Since vaccines targeting E6 also represent an important strategy for controlling HPV-associated lesions, we developed a codon-optimized HPV-16 E6 DNA vaccine (pNGVL4a-E6/opt) and characterized the E6-specific CD8+ T cell immune responses as well as the protective and therapeutic anti-tumor effects in vaccinated C57BL/6 mice. Our data indicated that transfection of human embryonic kidney cells (293 cells) with pNGVL4a-E6/opt resulted in highly efficient translation of E6. In addition, vaccination with pNGVL4a-E6/opt significantly enhanced E6-specific CD8+ T cell immune responses in C57BL/6 mice. Mice vaccinated with pNGVL4a-E6/opt are able to generate potent protective and therapeutic antitumor effects against challenge with E6-expressing tumor cell line, TC-1. Thus, DNA vaccines encoding a codon-optimized HPV-16 E6 may be a promising strategy for improving the potency of prophylactic and therapeutic HPV vaccines with potential clinical implications.

Gene codon composition determines differentiation-dependent expression of a viral capsid gene in keratinocytes in vitro and in vivo

By establishing mouse primary keratinocytes (KCs) in culture, we were able, for the first time, to express papillomavirus major capsid (L1) proteins by transient transfection of authentic or codon-modified L1 gene expression plasmids. We demonstrate in vitro and in vivo that gene codon composition is in part responsible for differentiation-dependent expression of L1 protein in KCs. L1 mRNA was present in similar amounts in differentiated and undifferentiated KCs transfected with authentic or codon-modified L1 genes and had a similar half-life, demonstrating that L1 protein production is posttranscriptionally regulated. We demonstrate further that KCs substantially change their tRNA profiles upon differentiation. Aminoacyl-tRNAs from differentiated KCs but not undifferentiated KCs enhanced the translation of authentic L1 mRNA, suggesting that differentiation-associated change to tRNA profiles enhances L1 expression in differentiated KCs. Thus, our data reveal a novel mechanism for regulation of gene expression utilized by a virus to direct viral capsid protein expression to the site of virion assembly in mature KCs. Analysis of two structural proteins of KCs, involucrin and keratin 14, suggests that translation of their mRNAs is also regulated, in association with KC differentiation in vitro, by a similar mechanism.

Epstein-Barr virus and cancer

EBV was the first human virus to be directly implicated in carcinogenesis. It infects >90% of the world's population. Although most humans coexist with the virus without serious sequelae, a small proportion will develop tumors. Normal host populations can have vastly different susceptibility to EBV-related tumors as demonstrated by geographical and immunological variations in the prevalence of these cancers. EBV has been implicated in the pathogenesis of Burkitt's lymphoma, Hodgkin's disease, non-Hodgkin's lymphoma, nasopharyngeal carcinoma, and lymphomas, as well as leiomyosarcomas arising in immunocompromised individuals. The presence of this virus has also been associated with epithelial malignancies arising in the gastric region and the breast, although some of this work remains in dispute. EBV uses its viral proteins, the actions of which mimic several growth factors, transcription factors, and antiapoptotic factors, to usurp control of the cellular pathways that regulate diverse homeostatic cellular functions. Recent advances in antiviral therapeutics, application of monoclonal antibodies, and generation of EBV-specific CTLs are beginning to show promise in the treatment of EBV-related disorders.

A synthetic E7 gene of human papillomavirus type 16 that yields enhanced expression of the protein in mammalian cells and is useful for DNA immunization studies

A synthetic E7 gene of human papillomavirus (HPV) type 16 was generated that consists entirely of preferred human codons. Expression analysis of the synthetic E7 gene in human and animal cells showed levels of E7 protein 20- to 100-fold higher than those obtained with wild-type E7. Enhanced expression of E7 protein resulted from highly efficient translation, as well as increased stability of the E7 mRNA due to its codon optimization. Higher levels of E7 protein in cells transfected with synthetic E7 correlated with significant loss of cell viability in various human cell lines. In contrast, lower E7 protein expression driven by the wild-type gene resulted in a slight induction of cell proliferation. Furthermore, mice inoculated with plasmids expressing the synthetic E7 gene produced significantly higher levels of E7 antibodies than littermates injected with wild-type E7, suggesting that synthetic E7 may be useful for DNA immunization studies and the development of genetic vaccines against HPV-16. In view of these results, we hypothesize that HPVs may have retained a pattern of G + C content and codon usage distinct from that of their host cells in response to selective pressure. Thus, the nonhuman codon bias may have been conserved by HPVs to prevent compromising viability of the host cells by excessive viral early protein expression, as well as to evade the immune system.

The extent of codon usage bias in human RNA viruses and its evolutionary origin

Revealing the determinants of codon usage bias is central to the understanding of factors governing viral evolution. Herein, we report the results of a survey of codon usage bias in a wide range of genetically and ecologically diverse human RNA viruses. This analysis showed that the overall extent of codon usage bias in RNA viruses is low and that there is little variation in bias between genes. Furthermore, the strong correlation between base and dinucleotide composition and codon usage bias suggested that mutation pressure rather than natural (translational) selection is the most important determinant of the codon bias observed. However, we also detected correlations between codon usage bias and some characteristics of viral genome structure and ecology, with increased bias in segmented and aerosol-transmitted viruses and decreased bias in vector-borne viruses. This suggests that translational selection may also have some influence in shaping codon usage bias.

Red queen dynamics of protein translation

We explore adaptive theories for the diversity of translational binding based on the genetic code viewed as a primitive mechanism of resistance. Modifying the set of codons bound by tRNA anticodon molecules or changing the specificity of binding, reduces the replication rate of translational parasites such as viruses. Increased translational efficiency of the parasite requires a high degree of specificity of host tRNAs for the parasite codons. This suggests that the genetic code might serve as the first line of defense against infection. We construct a red queen theory for translational diversity: a theory in which host-translational strategies- as defined by the degree of redundancy (a single anticodon binding many codons for a single amino acid) or degeneracy (many anticodons binding many codons for a single amino acid)-are constantly shifting through time to evade parasitism but where neither parasite nor host gain a systematic advantage.

Double-stranded RNA as a not-self alarm signal: to evade, most viruses purine-load their RNAs, but some (HTLV-1, Epstein-Barr) pyrimidine-load

For double-stranded RNA (dsRNA) to signal the presence of foreign (non-self) nucleic acid, self-RNA-self-RNA interactions should be minimized. Indeed, self-RNAs appear to have been fine-tuned over evolutionary time by the introduction of purines in clusters in the loop regions of stem-loop structures. This adaptation should militate against the "kissing" interactions which initiate formation of dsRNA. Our analyses of virus base compositions suggest that, to avoid triggering the host cell's dsRNA surveillance mechanism, most viruses purine-load their RNAs to resemble host RNAs ("stealth" strategy). However, some GC-rich latent viruses (HTLV-1, EBV) pyrimidine-load their RNAs. It is suggested that when virus production begins, these RNAs suddenly increase in concentration and impair host mRNA function by virtue of an excess of complementary "kissing" interactions ("surprise" strategy). Remarkably, the only mRNA expressed in the most fundamental form of EBV latency (the "EBNA-1 program") is purine-loaded. This apparent stealth strategy is reinforced by a simple sequence repeat which prefers purine-rich codons. During latent infection the EBNA-1 protein may evade recognition by cytotoxic T-cells, not by virtue of containing a simple sequence amino acid repeat as has been proposed, but by virtue of the encoding mRNA being purine-loaded to prevent interactions with host RNAs of either genic or non-genic origin.

Influence of genomic G+C content on average amino-acid composition of proteins from 59 bacterial species

The amino-acid composition of 23,490 proteins from 59 bacterial species was analyzed as a function of genomic G+C content. Observed amino-acid frequencies were compared with those expected from a neutral model assuming the absence of selection on average protein composition. Integral membrane proteins and non-integral membrane proteins were analyzed separately. The average deviation from this neutral model shows that there is a selective pressure increasing content in charged amino acids for non-integral membrane proteins, and content in hydrophobic amino acids for integral membrane proteins. Amino-acid frequencies were greatly influenced by genomic G+C content, but the influence was found to be often weaker than predicted. This may be evidence for a selective pressure, maintaining most amino-acid frequencies close to an optimal value. Concordance between the genetic code and protein composition is discussed in the light of this observation.

Repetitive sequence in the Epstein-Barr virus EBNA-3C gene is related to a family of minisatellite arrays in the human genome

A unique feature of the Epstein-Barr virus (EBV) genome is its high content of repetitive sequences. We identified a new human minisatellite element, tentatively designated MEB-1, that is similar to the 10 x "15bp" tandem repeat within the EBV nuclear antigen-3C (EBNA-3C) coding region. Southern blot analysis showed that the human genome has multiple copies of MEB-1-related repeats and that some of them are highly polymorphic. Both MEB-1 and the 10 x "15bp" repeat contain an octamer consensus GC[A/T]GG[A/T]GG, resembling the prokaryotic recombination signal chi. This octamer was also found in another EBV repeat sequence IR3 and the cellular GGA family of repeats that are related to IR3. Since the octamer motif is generally considered to have a role in the generation of a group of minisatellite DNA, these findings suggest that the four viral and cellular repeats have been generated through similar mechanisms promoted by the motif.

Prediction of function in DNA sequence analysis

Recognition of function of newly sequenced DNA fragments is an important area of computational molecular biology. Here we present an extensive review of methods for prediction of functional sites, tRNA, and protein-coding genes and discuss possible further directions of research in this area.

Statistical studies of biomolecular sequences: score-based methods

The massive accumulation of DNA and protein sequence data poses challenges and opportunities in terms of interpretation and analysis. This presentation reviews the method of score-based sequence analysis with the objectives of discerning distinctive segments in single sequences and identifying significant common segments in sequence comparisons. A number of new results are described here for both the theory and its applications. These include distributional theory involving several high scoring segments in single sequences, distribution formulas for general scoring regimes in multiple sequence comparisons, bounds for periodic scoring assignments, sensitivity analysis of genome composition and refinements on predicting exons and genes in DNA sequences.

The unusual nucleotide content of the HIV RNA genome results in a biased amino acid composition of HIV proteins

Extremely high frequencies of the A nucleotide are found in the RNA genomes of the lentivirus group of retroviruses. It is presently unknown what molecular force is responsible for this A-pressure. In this manuscript, we demonstrate a correlation between this 'A-pressure' and the amino acid-usage of the lentivirus family. We compared the amino acid composition of the Gag and Pol proteins of the human immunodeficiency viruses type 1 and 2 (HIV-1 and HIV-2) with that of the second group of human retroviruses; the human T-cell leukemia viruses type I and II (HTLV-I and HTLV-II). Differences in total amino acid content correlate with the preference for A-rich codons in the HIV genome. A pair-wise comparison of homologous amino acid positions in the Pol proteins indicates that both conservative and non-conservative changes can be accounted for by this A-bias. The putative molecular mechanism underlying this A-pressure and the evolutionary consequences are discussed.

Molecular evolution of herpesviruses: genomic and protein sequence comparisons

Phylogenetic reconstruction of herpesvirus evolution is generally founded on amino acid sequence comparisons of specific proteins. These are relevant to the evolution of the specific gene (or set of genes), but the resulting phylogeny may vary depending on the particular sequence chosen for analysis (or comparison). In the first part of this report, we compare 13 herpesvirus genomes by using a new multidimensional methodology based on distance measures and partial orderings of dinucleotide relative abundances. The sequences were analyzed with respect to (i) genomic compositional extremes; (ii) total distances within and between genomes; (iii) partial orderings among genomes relative to a set of sequence standards; (iv) concordance correlations of genome distances; and (v) consistency with the alpha-, beta-, gammaherpesvirus classification. Distance assessments within individual herpesvirus genomes show each to be quite homogeneous relative to the comparisons between genomes. The gammaherpesviruses, Epstein-Barr virus (EBV), herpesvirus saimiri, and bovine herpesvirus 4 are both diverse and separate from other herpesvirus classes, whereas alpha- and betaherpesviruses overlap. The analysis revealed that the most central genome (closest to a consensus herpesvirus genome and most individual herpesvirus sequences of different classes) is that of human herpesvirus 6, suggesting that this genome is closest to a progenitor herpesvirus. The shorter DNA distances among alphaherpesviruses supports the hypothesis that the alpha class is of relatively recent ancestry. In our collection, equine herpesvirus 1 (EHV1) stands out as the most central alphaherpesvirus, suggesting it may approximate an ancestral alphaherpesvirus. Among all herpesviruses, the EBV genome is closest to human sequences. In the DNA partial orderings, the chicken sequence collection is invariably as close as or closer to all herpesvirus sequences than the human sequence collection is, which may imply that the chicken (or other avian species) is a more natural or more ancient host of herpesviruses. In the second part of this report, evolutionary relationships among the 13 herpesvirus genomes are evaluated on the basis of recent methods of amino acid alignment applied to four essential protein sequences. In this analysis, the alignment of the two betaherpesviruses (human cytomegalovirus versus human herpesvirus 6) showed lower scores compared with alignments within alphaherpesviruses (i.e., among EHV1, herpes simplex virus type 1, varicella-zoster virus, pseudorabies virus type 1 and Marek's disease virus) and within gammaherpesviruses (EBV versus herpesvirus saimiri).(ABSTRACT TRUNCATED AT 400 WORDS)