Codon and amino acid usage in retroviral genomes is consistent with virus-specific nucleotide pressure

Dengue virus preferentially uses human and mosquito non-optimal codons

Codon optimality refers to the effect that codon composition has on messenger RNA (mRNA) stability and translation level and implies that synonymous codons are not silent from a regulatory point of view. Here, we investigated the adaptation of virus genomes to the host optimality code using mosquito-borne dengue virus (DENV) as a model. We demonstrated that codon optimality exists in mosquito cells and showed that DENV preferentially uses nonoptimal (destabilizing) codons and avoids codons that are defined as optimal (stabilizing) in either human or mosquito cells. Human genes enriched in the codons preferentially and frequently used by DENV are upregulated during infection, and so is the tRNA decoding the nonoptimal and DENV preferentially used codon for arginine. We found that adaptation during single-host passaging in human or mosquito cells results in the selection of synonymous mutations towards DENV's preferred nonoptimal codons that increase virus fitness. Finally, our analyses revealed that hundreds of viruses preferentially use nonoptimal codons, with those infecting a single host displaying an even stronger bias, suggesting that host-pathogen interaction shapes virus-synonymous codon choice.

Interferon-Stimulated Genes that Target Retrovirus Translation

The innate immune system, particularly the interferon (IFN) system, constitutes the initial line of defense against viral infections. IFN signaling induces the expression of interferon-stimulated genes (ISGs), and their products frequently restrict viral infection. Retroviruses like the human immunodeficiency viruses and the human T-lymphotropic viruses cause severe human diseases and are targeted by ISG-encoded proteins. Here, we discuss ISGs that inhibit the translation of retroviral mRNAs and thereby retrovirus propagation. The Schlafen proteins degrade cellular tRNAs and rRNAs needed for translation. Zinc Finger Antiviral Protein and RNA-activated protein kinase inhibit translation initiation factors, and Shiftless suppresses translation recoding essential for the expression of retroviral enzymes. We outline common mechanisms that underlie the antiviral activity of multifunctional ISGs and discuss potential antiretroviral therapeutic approaches based on the mode of action of these ISGs.

Different Patterns of Codon Usage and Amino Acid Composition across Primate Lentiviruses

A common feature of the mammalian Lentiviruses (family Retroviridae) is an RNA genome that contains an extremely high frequency of adenine (31.7-38.2%) while being extremely poor in cytosine (13.9-21.2%). Such a biased nucleotide composition has implications for codon usage, causing a striking difference between the frequency of synonymous codons in Lentiviruses and that in their hosts. To test whether primate Lentiviruses present differences in codon and amino acid composition, we assembled a dataset of genome sequences that includes SIV species infecting Old-World monkeys and African apes, HIV-2, and the four groups of HIV-1. Using principal component analysis, we found that HIV-1 shows a significant enrichment in adenine plus thymine in the third synonymous codon position and in adenine and guanine in the first and second nonsynonymous codon positions. Similarly, we observed an enrichment in adenine and in guanine in nonsynonymous first and second codon positions, which affects the amino acid composition of the proteins Gag, Pol, Vif, Vpr, Tat, Rev, Env, and Nef. This result suggests an effect of natural selection in shaping codon usage. Under the hypothesis that the use of synonyms in HIV-1 could reflect adaptation to that of genes expressed in specific cell types, we found a highly significant correlation between codon usage in HIV-1 and monocytes, which was remarkably higher than that with B and T lymphocytes. This finding is in line with the notion that monocytes represent an HIV-1 reservoir in infected patients, and it could help understand how this reservoir is established and maintained.

An investigation of codon usage pattern analysis in pancreatitis associated genes

BACKGROUND

Pancreatitis is an inflammatory disorder resulting from the autoactivation of trypsinogen in the pancreas. The genetic basis of the disease is an old phenomenon, and evidence is accumulating for the involvement of synonymous/non-synonymous codon variants in disease initiation and progression.

RESULTS

The present study envisaged a panel of 26 genes involved in pancreatitis for their codon choices, compositional analysis, relative dinucleotide frequency, nucleotide disproportion, protein physical properties, gene expression, codon bias, and interrelated of all these factors. In this set of genes, gene length was positively correlated with nucleotide skews and codon usage bias. Codon usage of any gene is dependent upon its AT and GC component; however, AGG, CGT, and CGA encoding for Arg, TCG for Ser, GTC for Val, and CCA for Pro were independent of nucleotide compositions. In addition, Codon GTC showed a correlation with protein properties, isoelectric point, instability index, and frequency of basic amino acids. We also investigated the effect of various evolutionary forces in shaping the codon usage choices of genes.

CONCLUSIONS

This study will enable us to gain insight into the molecular signatures associated with the disease that might help identify more potential genes contributing to enhanced risk for pancreatitis. All the genes associated with pancreatitis are generally associated with physiological function, and mutations causing loss of function, over or under expression leads to an ailment. Therefore, the present study attempts to envisage the molecular signature in a group of genes that lead to pancreatitis in case of malfunction.

Compositional biases in RNA viruses: Causes, consequences and applications

If each of the four nucleotides were represented equally in the genomes of viruses and the hosts they infect, each base would occur at a frequency of 25%. However, this is not observed in nature. Similarly, the order of nucleotides is not random (e.g., in the human genome, guanine follows cytosine at a frequency of ~0.0125, or a quarter the number of times predicted by random representation). Codon usage and codon order are also nonrandom. Furthermore, nucleotide and codon biases vary between species. Such biases have various drivers, including cellular proteins that recognize specific patterns in nucleic acids, that once triggered, induce mutations or invoke intrinsic or innate immune responses. In this review we examine the types of compositional biases identified in viral genomes and current understanding of the evolutionary mechanisms underpinning these trends. Finally, we consider the potential for large scale synonymous recoding strategies to engineer RNA virus vaccines, including those with pandemic potential, such as influenza A virus and Severe Acute Respiratory Syndrome Coronavirus Virus 2. This article is categorized under: RNA in Disease and Development > RNA in Disease RNA Evolution and Genomics > Computational Analyses of RNA RNA Interactions with Proteins and Other Molecules > Protein-RNA Recognition.

Nucleotide Composition and Codon Usage Across Viruses and Their Respective Hosts

The genetic material of the three domains of life (Bacteria, Archaea, and Eukaryota) is always double-stranded DNA, and their GC content (molar content of guanine plus cytosine) varies between ≈ 13% and ≈ 75%. Nucleotide composition is the simplest way of characterizing genomes. Despite this simplicity, it has several implications. Indeed, it is the main factor that determines, among other features, dinucleotide frequencies, repeated short DNA sequences, and codon and amino acid usage. Which forces drive this strong variation is still a matter of controversy. For rather obvious reasons, most of the studies concerning this huge variation and its consequences, have been done in free-living organisms. However, no recent comprehensive study of all known viruses has been done (that is, concerning all available sequences). Viruses, by far the most abundant biological entities on Earth, are the causative agents of many diseases. An overview of these entities is important also because their genetic material is not always double-stranded DNA: indeed, certain viruses have as genetic material single-stranded DNA, double-stranded RNA, single-stranded RNA, and/or retro-transcribing. Therefore, one may wonder if what we have learned about the evolution of GC content and its implications in prokaryotes and eukaryotes also applies to viruses. In this contribution, we attempt to describe compositional properties of ∼ 10,000 viral species: base composition (globally and according to Baltimore classification), correlations among non-coding regions and the three codon positions, and the relationship of the nucleotide frequencies and codon usage of viruses with the same feature of their hosts. This allowed us to determine how the base composition of phages strongly correlate with the value of their respective hosts, while eukaryotic viruses do not (with fungi and protists as exceptions). Finally, we discuss some of these results concerning codon usage: reinforcing previous results, we found that phages and hosts exhibit moderate to high correlations, while for eukaryotes and their viruses the correlations are weak or do not exist.

Mouse APOBEC3 Restriction of Retroviruses

Apolipoprotein B mRNA editing enzyme, catalytic peptide 3 (APOBEC3) proteins are critical host proteins that counteract and prevent the replication of retroviruses. Unlike the genome of humans and other species, the mouse genome encodes a single Apobec3 gene, which has undergone positive selection, as reflected by the allelic variants found in different inbred mouse strains. This positive selection was likely due to infection by various mouse retroviruses, which have persisted in their hosts for millions of years. While mouse retroviruses are inhibited by APOBEC3, they nonetheless still remain infectious, likely due to the actions of different viral proteins that counteract this host factor. The study of viruses in their natural hosts provides important insight into their co-evolution.

Modulation of HIV-1 Gag/Gag-Pol frameshifting by tRNA abundance

A hallmark of translation in human immunodeficiency virus type 1 (HIV-1) is a –1 programmed ribosome frameshifting event that produces the Gag-Pol fusion polyprotein. The constant Gag to Gag-Pol ratio is essential for the virion structure and infectivity. Here we show that the frameshifting efficiency is modulated by Leu-tRNA^Leu that reads the UUA codon at the mRNA slippery site. This tRNA^Leu isoacceptor is particularly rare in human cell lines derived from T-lymphocytes, the cells that are targeted by HIV-1. When UUA decoding is delayed, the frameshifting follows an alternative route, which maintains the Gag to Gag-Pol ratio constant. A second potential slippery site downstream of the first one is normally inefficient but can also support –1-frameshifting when altered by a compensatory resistance mutation in response to current antiviral drug therapy. Together these different regimes allow the virus to maintain a constant –1-frameshifting efficiency to ensure successful virus propagation.

KHNYN is essential for the zinc finger antiviral protein (ZAP) to restrict HIV-1 containing clustered CpG dinucleotides

CpG dinucleotides are suppressed in most vertebrate RNA viruses, including HIV-1, and introducing CpGs into RNA virus genomes inhibits their replication. The zinc finger antiviral protein (ZAP) binds regions of viral RNA containing CpGs and targets them for degradation. ZAP does not have enzymatic activity and recruits other cellular proteins to inhibit viral replication. We found that KHNYN, a protein with no previously known function, interacts with ZAP. KHNYN overexpression selectively inhibits HIV-1 containing clustered CpG dinucleotides and this requires ZAP and its cofactor TRIM25. KHNYN requires both its KH-like domain and NYN endonuclease domain for antiviral activity. Crucially, depletion of KHNYN eliminated the deleterious effect of CpG dinucleotides on HIV-1 RNA abundance and infectious virus production and also enhanced the production of murine leukemia virus. Overall, we have identified KHNYN as a novel cofactor for ZAP to target CpG-containing retroviral RNA for degradation.

Analysis of codon usage patterns and influencing factors in Nipah virus

•

Nipah virus (NiV) genes are AT-rich.

•

Codon usage bias of NiV genes is low.

•

Patterns of codon usage bias differ across the genomes of NiV.

•

Both mutation pressure and natural selection influenced codon usage bias of NiV genes.

Codon usage bias (CUB) is the unequal usage of synonymous codons of an amino acid in which some codons are used more often than others and is widely used in understanding molecular biology, genetics, and functional regulation of gene expression. Nipah virus (NiV) is an emerging zoonotic paramyxovirus that causes fatal disease in both humans and animals. NiV was first identified during an outbreak of a disease in Malaysia in 1998 and then occurred periodically since 2001 in India, Bangladesh, and the Philippines. We used bioinformatics tools to analyze the codon usage patterns in a genome-wide manner among 11 genomes of NiV as no work was reported yet. The compositional properties revealed that the overall GC and AT contents were 41.96 and 58.04%, respectively i.e. Nipah virus genes were AT-rich. Correlation analysis between overall nucleotide composition and its 3^rd codon position suggested that both mutation pressure and natural selection might influence the CUB across Nipah genomes. Neutrality plot revealed natural selection might have played a major role while mutation pressure had a minor role in shaping the codon usage bias in NiV genomes.

Increasing the CpG dinucleotide abundance in the HIV-1 genomic RNA inhibits viral replication

BACKGROUND

The human immunodeficiency virus type 1 (HIV-1) structural protein Gag is necessary and sufficient to form viral particles. In addition to encoding the amino acid sequence for Gag, the underlying RNA sequence could encode cis-acting elements or nucleotide biases that are necessary for viral replication. Furthermore, RNA sequences that inhibit viral replication could be suppressed in gag. However, the functional relevance of RNA elements and nucleotide biases that promote or repress HIV-1 replication remain poorly understood.

RESULTS

To characterize if the RNA sequence in gag controls HIV-1 replication, the matrix (MA) region was codon modified, allowing the RNA sequence to be altered without affecting the protein sequence. Codon modification of nucleotides (nt) 22-261 or 22-378 in gag inhibited viral replication by decreasing genomic RNA (gRNA) abundance, gRNA stability, Gag expression, virion production and infectivity. Comparing the effect of these point mutations to deletions of the same region revealed that the mutations inhibited infectious virus production while the deletions did not. This demonstrated that codon modification introduced inhibitory sequences. There is a much lower than expected frequency of CpG dinucleotides in HIV-1 and codon modification introduced a substantial increase in CpG abundance. To determine if they are necessary for inhibition of HIV-1 replication, codons introducing CpG dinucleotides were mutated back to the wild type codon, which restored efficient Gag expression and infectious virion production. To determine if they are sufficient to inhibit viral replication, CpG dinucleotides were inserted into gag in the absence of other changes. The increased CpG dinucleotide content decreased HIV-1 infectivity and viral replication.

CONCLUSIONS

The HIV-1 RNA sequence contains low abundance of CpG dinucleotides. Increasing the abundance of CpG dinucleotides inhibits multiple steps of the viral life cycle, providing a functional explanation for why CpG dinucleotides are suppressed in HIV-1.

HIV-1 tolerates changes in A-count in a small segment of the pol gene

Background

The HIV-1 RNA genome has a biased nucleotide composition with a surplus of As. Several hypotheses have been put forward to explain this striking phenomenon, but the A-count of the HIV-1 genome has thus far not been systematically manipulated. The reason for this reservation is the likelihood that known and unknown sequence motifs will be affected by such a massive mutational approach, thus resulting in replication-impaired virus mutants. We present the first attempt to increase and decrease the A-count in a relatively small polymerase (pol) gene segment of HIV-1 RNA.

Results

To minimize the mutational impact, a new mutational approach was developed that is inspired by natural sequence variation as present in HIV-1 isolates. This phylogeny-instructed mutagenesis allowed us to create replication-competent HIV-1 mutants with a significantly increased or decreased local A-count. The local A-count of the wild-type (wt) virus (40.2%) was further increased to 46.9% or reduced to 31.7 and 26.3%. These HIV-1 variants replicate efficiently in vitro, despite the fact that the pol changes cause a quite profound move in HIV–SIV sequence space.

Conclusions

Extrapolating these results to the complete 9 kb RNA genome, we may cautiously suggest that the A-rich signature does not have to be maintained. This survey also provided clues that silent codon changes, in particular from G-to-A, determine the subtype-specific sequence signatures.

Differential Codon Usage Pattern of HIV-1 tat Gene in Those with Slower and Faster Rates of Disease Progression

Codon usage has been identified as one of the most important factors that influence gene expression. The frequencies with which the different codons are used vary significantly between different organisms and also between the genes within the same organism. HIV has a remarkable nucleotide composition with an above average percentage of "A" nucleotides resulting in a codon usage pattern different from that of the human host. In this study, we compared the codon usage pattern of HIV-1 genes among different groups of HIV disease progressors to understand the influence of differential codon usage pattern on the pathogenic manifestation in the host.

Biased Nucleotide Composition and Differential Codon Usage Pattern in HIV-1 and HIV-2

HIV-1 and HIV-2 are closely related retroviruses with differences in pathogenicity and geographic distribution. HIV-2 infection is associated with slower disease progression and transmission, longer latency period, low or undetectable plasmatic viral loads, and reduced likelihood of progression to AIDS, compared to HIV-1. In this investigation, we analyzed HIV-2 genes and genomes and compared them with that of HIV-1 belonging to various subtypes. Comparative analysis of the effective number of codons (ENC) for each of the nine genes of the two viruses revealed that the tat gene of HIV-2 had a higher ENC value compared to HIV-1 tat, reflecting lower levels of expression of HIV-2 tat. Lower levels of tat protein particularly during the early stages of infection could result in a lower viral load, lower viral set point, and delayed progression of disease in HIV-2-infected individuals compared to HIV-1-infected subjects. Furthermore, the GC3 composition of the regulatory genes of HIV-2 was ≥50%, suggesting a firm effort by these viruses to adapt themselves to evolutionary survival. We hypothesize that differential codon usage could be one of the possible factors that could contribute to the diminished pathogenicity of HIV-2 in the host as compared to HIV-1.

The thermodynamics of Pr55Gag-RNA interaction regulate the assembly of HIV

The interactions that occur during HIV Pr55^Gag oligomerization and genomic RNA packaging are essential elements that facilitate HIV assembly. However, mechanistic details of these interactions are not clearly defined. Here, we overcome previous limitations in producing large quantities of full-length recombinant Pr55^Gag that is required for isothermal titration calorimetry (ITC) studies, and we have revealed the thermodynamic properties of HIV assembly for the first time. Thermodynamic analysis showed that the binding between RNA and HIV Pr55^Gag is an energetically favourable reaction (ΔG<0) that is further enhanced by the oligomerization of Pr55^Gag. The change in enthalpy (ΔH) widens sequentially from: (1) Pr55^Gag-Psi RNA binding during HIV genome selection; to (2) Pr55^Gag-Guanosine Uridine (GU)-containing RNA binding in cytoplasm/plasma membrane; and then to (3) Pr55^Gag-Adenosine(A)-containing RNA binding in immature HIV. These data imply the stepwise increments of heat being released during HIV biogenesis may help to facilitate the process of viral assembly. By mimicking the interactions between A-containing RNA and oligomeric Pr55^Gag in immature HIV, it was noted that a p6 domain truncated Pr50^{Gag Δp6} is less efficient than full-length Pr55^Gag in this thermodynamic process. These data suggest a potential unknown role of p6 in Pr55^Gag-Pr55^Gag oligomerization and/or Pr55^Gag-RNA interaction during HIV assembly. Our data provide direct evidence on how nucleic acid sequences and the oligomeric state of Pr55^Gag regulate HIV assembly.

Formation of any virus particle will require energy, yet the precise biophysical properties that drive the formation of HIV particles remain undefined. Isothermal titration calorimetry (ITC) is a biophysical technique that is the gold standard to reveal parameters governing biochemical and biophysical reaction. However, ITC requires large amount of proteins for analysis. As large quantities of full-length recombinant HIV Pr55^Gag proteins have not been available in the past 30 years due to technical limitation, a comprehensive thermodynamic analysis of full-length HIV Pr55^Gag has not been possible. Here, we have generated sufficient amount of full-length recombinant HIV Pr55^Gag protein for isothermal titration calorimetry analysis. Our analyses have shown that the major interactions amongst HIV proteins and RNA sequences during viral assembly are energetically favourable reactions. In other words, HIV Pr55^Gag proteins and viral RNA have evolved to overcome the energy barrier for virus formation by utilising energy obtained from protein-RNA interactions in order to facilitate the viral assembly process. Furthermore, HIV also use the oligomeric states of HIV Pr55^Gag proteins and the RNA sequences as means to regulate the viral assembly process.

Impact of the biased nucleotide composition of viral RNA genomes on RNA structure and codon usage

We are interested in the influence of nucleotide composition on the fundamental characteristics of the virus RNA genome. Most RNA viruses have genomes with a distinct nucleotide composition, e.g. ranging from minimally 12.9 % to maximally 40.3 % (C- and U-count, respectively, in coronavirus HKU). We present a global analysis of diverse virus types, including plus-strand, minus-strand and double-strand RNA viruses, for the impact of this nucleotide preference on the predicted structure of the RNA genome that is packaged in virion particles and on the codon usage in the viral open reading frames. Several virus-specific features will be described, but also some general conclusions were drawn. Without exception, the virus-specific nucleotide bias was enriched in the unpaired, single-stranded regions of the RNA genome, thus creating an even more striking virus-specific signature. We present a simple mechanism that is based on elementary aspects of RNA structure folding to explain this general trend. In general, the nucleotide bias was the major determinant of the virus-specific codon usages, thus limiting a role for codon selection and translational control. We will discuss molecular and evolutionary scenarios that may be responsible for the diverse nucleotide biases of RNA viruses.

Analysis of the synonymous codon usage bias in recently emerged enterovirus D68 strains

Understanding the codon usage pattern of a pathogen and relationship between pathogen and host's codon usage patterns has fundamental and applied interests. Enterovirus D68 (EV-D68) is an emerging pathogen with a potentially high public health significance. In the present study, the synonymous codon usage bias of 27 recently emerged, and historical EV-D68 strains was analyzed. In contrast to previously studied enteroviruses (enterovirus 71 and poliovirus), EV-D68 and human host have a high discrepancy between favored codons. Analysis of viral synonymous codon usage bias metrics, viral nucleotide/dinucleotide compositional parameters, and viral protein properties showed that mutational pressure is more involved in shaping the synonymous codon usage bias of EV-D68 than translation selection. Computation of codon adaptation indices allowed to estimate expression potential of the EV-D68 genome in several commonly used laboratory animals. This approach requires experimental validation and may provide an auxiliary tool for the rational selection of laboratory animals to model emerging viral diseases. Enterovirus D68 genome compositional and codon usage data can be useful for further pathogenesis, animal model, and vaccine design studies.

Nucleotide composition of the Zika virus RNA genome and its codon usage

BACKGROUND

RNA viruses have genomes with a distinct nucleotide composition and codon usage. We present the global characteristics of the RNA genome of Zika virus (ZIKV), an emerging pathogen within the Flavivirus genus. ZIKV was first isolated in 1947 in Uganda, caused a widespread epidemic in South and Central America and the Caribbean in 2015 and has recently been associated with microcephaly in newborns.

METHODS

The nearly 11 kb positive-stranded RNA genome of ZIKV was analyzed for its nucleotide composition, also in the context of the folded RNA molecule. Nucleotide trends were investigated along the genome length by skew analyses and we analyzed the codons used for translation of the ZIKV proteins.

RESULTS

ZIKV RNA has a biased nucleotide composition in being purine-rich and pyrimidine-poor. This preference for purines is a general characteristic of the mosquito-borne and tick-borne flaviviruses. The virus-specific nucleotide bias is further enriched in the unpaired, single-stranded regions of the structured ZIKV RNA genome, thus further imposing this ZIKV-specific signature. The codons used for translation of the ZIKV proteins is also unusual, but we show that it is the underlying bias in nucleotide composition of the viral RNA that largely dictates these codon preferences.

CONCLUSIONS

The ZIKV RNA genome has a biased nucleotide composition that dictates the codon usage of this flavivirus. We discuss the evolutionary scenarios and molecular mechanisms that may be responsible for these distinctive ZIKV RNA genome features.

Innate immune system activation by viral RNA: How to predict it?

The immune system is able to identify foreign pathogens via different pathways. In the case of viral infection, recognition of the viral RNA is a crucial step, and many efforts have been made to understand which features of viral RNA are detected by the immune system. The biased viral RNA composition, measured as host-virus nucleotidic divergence, or CpG enrichment, has been proposed as salient signal. Peculiar structural features of these RNA could also be related to the immune system activation. Here, we gather multiple datasets and proceed to a meta-analysis to uncover the best predictors of immune system activation by viral RNA. "A" nucleotide content and Minimum Folding Energy are good predictors, and are more easily generalized than more complex indicators suggested previously. As RNA composition and structure are highly correlated, we suggest further experiments on synthetic sequences to identify the viral RNA sensing mechanisms by immune system receptors.

Inferring the hosts of coronavirus using dual statistical models based on nucleotide composition

Many coronaviruses are capable of interspecies transmission. Some of them have caused worldwide panic as emerging human pathogens in recent years, e.g., severe acute respiratory syndrome coronavirus (SARS-CoV) and Middle East respiratory syndrome coronavirus (MERS-CoV). In order to assess their threat to humans, we explored to infer the potential hosts of coronaviruses using a dual-model approach based on nineteen parameters computed from spike genes of coronaviruses. Both the support vector machine (SVM) model and the Mahalanobis distance (MD) discriminant model achieved high accuracies in leave-one-out cross-validation of training data consisting of 730 representative coronaviruses (99.86% and 98.08% respectively). Predictions on 47 additional coronaviruses precisely conformed to conclusions or speculations by other researchers. Our approach is implemented as a web server that can be accessed at http://bioinfo.ihb.ac.cn/seq2hosts.

Homology-independent metrics for comparative genomics

A mainstream procedure to analyze the wealth of genomic data available nowadays is the detection of homologous regions shared across genomes, followed by the extraction of biological information from the patterns of conservation and variation observed in such regions. Although of pivotal importance, comparative genomic procedures that rely on homology inference are obviously not applicable if no homologous regions are detectable. This fact excludes a considerable portion of “genomic dark matter” with no significant similarity — and, consequently, no inferred homology to any other known sequence — from several downstream comparative genomic methods. In this review we compile several sequence metrics that do not rely on homology inference and can be used to compare nucleotide sequences and extract biologically meaningful information from them. These metrics comprise several compositional parameters calculated from sequence data alone, such as GC content, dinucleotide odds ratio, and several codon bias metrics. They also share other interesting properties, such as pervasiveness (patterns persist on smaller scales) and phylogenetic signal. We also cite examples where these homology-independent metrics have been successfully applied to support several bioinformatics challenges, such as taxonomic classification of biological sequences without homology inference. They where also used to detect higher-order patterns of interactions in biological systems, ranging from detecting coevolutionary trends between the genomes of viruses and their hosts to characterization of gene pools of entire microbial communities. We argue that, if correctly understood and applied, homology-independent metrics can add important layers of biological information in comparative genomic studies without prior homology inference.

On the nucleotide composition and structure of retroviral RNA genomes

Retroviral RNA genomes display a rich variety in their nucleotide composition. For instance, the single-stranded RNA genome of human T cell leukemia virus (HTLV-1) is C-rich and G-poor and that of the human immunodeficiency virus (HIV-1) is A-rich and C-poor. Animal retroviruses add further variation to this unexplained, but many times remarkable virus-specific property. We previously described that the nucleotide bias is even more extreme in the unpaired regions of the structured HIV-1 RNA genome, which has been probed by SHAPE technology. We now document that the same trend is apparent for the MFold-predicted RNA structure of HIV-1 RNA and subsequently investigated the predicted structures of the RNA genomes of other retroviruses. We conclude that all virus-specific signatures are enhanced for the unpaired nucleotides in the RNA genome. Consequently, the differences in nucleotide count between the diverse human and animal retroviruses are further exposed in the single stranded genome regions. We used a skew analysis to visualize these striking differences in nucleotide usage. Evolutionary events responsible for these nucleotide signatures will be discussed.

Large-scale nucleotide optimization of simian immunodeficiency virus reduces its capacity to stimulate type I interferon in vitro

Lentiviral RNA genomes present a strong bias in their nucleotide composition with extremely high frequencies of A nucleotide in human immunodeficiency virus type 1 (HIV-1) and simian immunodeficiency virus (SIV). Based on the observation that human optimization of RNA virus gene fragments may abolish their ability to stimulate the type I interferon (IFN-I) response, we identified the most biased sequences along the SIV genome and showed that they are the most potent IFN-I stimulators. With the aim of designing an attenuated SIV genome based on a reduced capacity to activate the IFN-I response, we synthesized artificial SIV genomes whose biased sequences were optimized toward macaque average nucleotide composition without altering their regulatory elements or amino acid sequences. A synthetic SIV optimized with 169 synonymous mutations in gag and pol genes showed a 100-fold decrease in replicative capacity. Interestingly, a synthetic SIV optimized with 70 synonymous mutations in pol had a normal replicative capacity. Its ability to stimulate IFN-I was reduced when infected cells were cocultured with reporter cells. IFN regulatory factor 3 (IRF3) transcription factor was required for IFN-I stimulation, implicating cytosolic sensors in the detection of SIV-biased RNA in infected cells. No reversion of introduced mutations was observed for either of the optimized viruses after 10 serial passages. In conclusion, we have designed large-scale nucleotide-modified SIVs that may display attenuated pathogenic potential.

IMPORTANCE

In this study, we synthesized artificial SIV genomes in which the most hyperbiased sequences were optimized to bring them closer to the nucleotide composition of the macaque SIV host. Interestingly, we generated a stable synthetic SIV optimized with 70 synonymous mutations in pol gene, which had a normal replicative capacity but a reduced ability to stimulate type I IFN. This demonstrates the possibility to rationally change viral nucleotide composition to design replicative and genetically stable lentiviruses with attenuated pathogenic potentials.

Human Retrovirus Codon Usage from tRNA Point of View: Therapeutic Insights

The purpose of this study was to investigate the balance between transfer ribonucleic acid (tRNA) supply and demand in retrovirus-infected cells, seeking the best targets for antiretroviral therapy based on the hypothetical tRNA Inhibition Therapy (TRIT). Codon usage and tRNA gene data were retrieved from public databases. Based on logistic principles, a therapeutic score (T-score) was calculated for all sense codons, in each retrovirus-host system. Codons that are critical for viral protein translation, but not as critical for the host, have the highest T-score values. Theoretically, inactivating the cognate tRNA species should imply a severe reduction of the elongation rate during viral mRNA translation. We developed a method to predict tRNA species critical for retroviral protein synthesis. Four of the best TRIT targets in HIV-1 and HIV-2 encode Large Hydrophobic Residues (LHR), which have a central role in protein folding. One of them, codon CUA, is also a TRIT target in both HTLV-1 and HTLV-2. Therefore, a drug designed for inactivating or reducing the cytoplasmatic concentration of tRNA species with anticodon TAG could attenuate significantly both HIV and HTLV protein synthesis rates. Inversely, replacing codons ending in UA by synonymous codons should increase the expression, which is relevant for DNA vaccine design.

Causes and implications of codon usage bias in RNA viruses

Choice of synonymous codons depends on nucleotide/dinucleotide composition of the genome (termed mutational pressure) and relative abundance of tRNAs in a cell (translational pressure). Mutational pressure is commonly simplified to genomic GC content; however mononucleotide and dinucleotide frequencies in different genomes or mRNAs may vary significantly, especially in RNA viruses. A series of in silico shuffling algorithms were developed to account for these features and analyze the relative impact of mutational pressure components on codon usage bias in RNA viruses. Total GC content was a poor descriptor of viral genome composition and causes of codon usage bias. Genomic nucleotide content was the single most important factor of synonymous codon usage. Moreover, the choice between compatible amino acids (e.g., leucine and isoleucine) was strongly affected by genomic nucleotide composition. Dinucleotide composition at codon positions 2-3 had additional effect on codon usage. Together with mononucleotide composition bias, it could explain almost the entire codon usage bias in RNA viruses. On the other hand, strong dinucleotide content bias at codon position 3-1 found in some viruses had very little effect on codon usage. A hypothetical innate immunity sensor for CpG in RNA could partially explain the codon usage bias, but due to dependence of virus translation upon biased host translation machinery, experimental studies are required to further explore the source of dinucleotide bias in RNA viruses.

The A-nucleotide preference of HIV-1 in the context of its structured RNA genome

A bipartition of HIV-1 RNA genome sequences into single- and double-stranded nucleotides is possible based on the secondary structure model of a complete 9 kb genome. Subsequent analysis revealed that the well-known lentiviral property of A-accumulation is profoundly present in single-stranded domains, yet absent in double-stranded domains. Mutational rate analysis by means of an unrestricted model of nucleotide substitution suggests the presence of an evolutionary equilibrium to preserve this biased nucleotide distribution.

The biased nucleotide composition of the HIV genome: a constant factor in a highly variable virus

Viruses often deviate from their hosts in the nucleotide composition of their genomes. The RNA genome of the lentivirus family of retroviruses, including human immunodeficiency virus (HIV), contains e.g. an above average percentage of adenine (A) nucleotides, while being extremely poor in cytosine (C). Such a deviant base composition has implications for the amino acids that are encoded by the open reading frames (ORFs), both in the requirement of specific tRNA species and in the preference for amino acids encoded by e.g. A-rich codons. Nucleotide composition does obviously affect the secondary and tertiary structure of the RNA genome and its biological functions, but it does also influence phylogenetic analysis of viral genome sequences, and possibly the activity of the integrated DNA provirus. Over time, the nucleotide composition of the HIV-1 genome is exceptionally conserved, varying by less than 1% per base position per isolate within either group M, N, or O during 1983–2009. This extreme stability of the nucleotide composition may possibly be achieved by negative selection, perhaps conserving semi-stable RNA secondary structure as reverse transcription would be significantly affected for a less A-rich genome where secondary structures are expected to be more stable and thus more difficult to unfold.

This review will discuss all aspects of the lentiviral genome composition, both of the RNA and of its derived double-stranded DNA genome, with a focus on HIV-1, the nucleotide composition over time, the effects of artificially humanized codons as well as contributions of immune system pressure on HIV nucleotide bias.

Global co-existence of two evolutionary lineages of parvovirus B19 1a, different in genome-wide synonymous positions

Parvovirus B19 (B19V) can cause infection in humans. To date, three genotypes of B19V, with subtypes, are known, of which genotype 1a is the most prevalent genotype in the Western world. We sequenced the genome of B19V strains of 65 asymptomatic, recently infected Dutch blood donors, to investigate the spatio-temporal distribution of B19V strains, in the years 2003-2009. The sequences were compared to B19V sequences from Dutch patients with fifth disease, and to global B19V sequences as available from GenBank. All Dutch B19V strains belonged to genotype 1a. Phylogenetic analysis of the strains from Dutch blood donors showed that two groups of genotype 1a co-exist. A clear-cut division into the two groups was also found among the B19V strains from Dutch patients, and among the B19V sequences in GenBank. The two groups of genotype 1a co-exist around the world and do not appear to differ in their ability to cause disease. Strikingly, the two groups of B19V predominantly differ in synonymous mutations, distributed throughout the entire genome of B19V. We propose to call the two groups of B19V genotype 1a respectively subtype 1a1 and 1a2.

The biased nucleotide composition of HIV-1 triggers type I interferon response and correlates with subtype D increased pathogenicity

The genome of human immunodeficiency virus (HIV) has an average nucleotide composition strongly biased as compared to the human genome. The consequence of such nucleotide composition on HIV pathogenicity has not been investigated yet. To address this question, we analyzed the role of nucleotide bias of HIV-derived nucleic acids in stimulating type-I interferon response in vitro. We found that the biased nucleotide composition of HIV is detected in human cells as compared to humanized sequences, and triggers a strong innate immune response, suggesting the existence of cellular immune mechanisms able to discriminate RNA sequences according to their nucleotide composition or to detect specific secondary structures or linear motifs within biased RNA sequences. We then extended our analysis to the entire genome scale by testing more than 1300 HIV-1 complete genomes to look for an association between nucleotide composition of HIV-1 group M subtypes and their pathogenicity. We found that subtype D, which has an increased pathogenicity compared to the other subtypes, has the most divergent nucleotide composition relative to the human genome. These data support the hypothesis that the biased nucleotide composition of HIV-1 may be related to its pathogenicity.

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.

HIV-1 modulates the tRNA pool to improve translation efficiency

Despite its poorly adapted codon usage, HIV-1 replicates and is expressed extremely well in human host cells. HIV-1 has recently been shown to package non-lysyl transfer RNAs (tRNAs) in addition to the tRNA(Lys) needed for priming reverse transcription and integration of the HIV-1 genome. By comparing the codon usage of HIV-1 genes with that of its human host, we found that tRNAs decoding codons that are highly used by HIV-1 but avoided by its host are overrepresented in HIV-1 virions. In particular, tRNAs decoding A-ending codons, required for the expression of HIV's A-rich genome, are highly enriched. Because the affinity of Gag-Pol for all tRNAs is nonspecific, HIV packaging is most likely passive and reflects the tRNA pool at the time of viral particle formation. Codon usage of HIV-1 early genes is similar to that of highly expressed host genes, but codon usage of HIV-1 late genes was better adapted to the selectively enriched tRNA pool, suggesting that alterations in the tRNA pool are induced late in viral infection. If HIV-1 genes are adapting to an altered tRNA pool, codon adaptation of HIV-1 may be better than previously thought.

Local sequence targeting in the AID/APOBEC family differentially impacts retroviral restriction and antibody diversification

Nucleic acid cytidine deaminases of the activation-induced deaminase (AID)/APOBEC family are critical players in active and innate immune responses, playing roles as target-directed, purposeful mutators. AID specifically deaminates the host immunoglobulin (Ig) locus to evolve antibody specificity, whereas its close relative, APOBEC3G (A3G), lethally mutates the genomes of retroviral pathogens such as HIV. Understanding the basis for the target-specific action of these enzymes is essential, as mistargeting poses significant risks, potentially promoting oncogenesis (AID) or fostering drug resistance (A3G). AID prefers to deaminate cytosine in WRC (W = A/T, R = A/G) motifs, whereas A3G favors deamination of CCC motifs. This specificity is largely dictated by a single, divergent protein loop in the enzyme family that recognizes the DNA sequence. Through grafting of this substrate-recognition loop, we have created enzyme variants of A3G and AID with altered local targeting to directly evaluate the role of sequence specificity on immune function. We find that grafted loops placed in the A3G scaffold all produced efficient restriction of HIV but that foreign loops in the AID scaffold compromised hypermutation and class switch recombination. Local targeting, therefore, appears alterable for innate defense against retroviruses by A3G but important for adaptive antibody maturation catalyzed by AID. Notably, AID targeting within the Ig locus is proportionally correlated to its in vitro ability to target WRC sequences rather than non-WRC sequences. Although other mechanisms may also contribute, our results suggest that local sequence targeting by AID/APOBEC3 enzymes represents an elegant example of co-evolution of enzyme specificity with its target DNA sequence.

Use of nucleotide composition analysis to infer hosts for three novel picorna-like viruses

Nearly complete genome sequences of three novel RNA viruses were acquired from the stool of an Afghan child. Phylogenetic analysis indicated that these viruses belong to the picorna-like virus superfamily. Because of their unique genomic organization and deep phylogenetic roots, we propose these viruses, provisionally named calhevirus, tetnovirus-1, and tetnovirus-2, as prototypes of new viral families. A newly developed nucleotide composition analysis (NCA) method was used to compare mononucleotide and dinucleotide frequencies for RNA viruses infecting mammals, plants, or insects. Using a large training data set of 284 representative picornavirus-like genomic sequences with defined host origins, NCA correctly identified the kingdom or phylum of the viral host for >95% of picorna-like viruses. NCA predicted an insect host origin for the 3 novel picorna-like viruses. Their presence in human stool therefore likely reflects ingestion of insect-contaminated food. As metagenomic analyses of different environments and organisms continue to yield highly divergent viral genomes NCA provides a rapid and robust method to identify their likely cellular hosts.

The remarkable frequency of human immunodeficiency virus type 1 genetic recombination

The genetic diversity of human immunodeficiency virus type 1 (HIV-1) results from a combination of point mutations and genetic recombination, and rates of both processes are unusually high. This review focuses on the mechanisms and outcomes of HIV-1 genetic recombination and on the parameters that make recombination so remarkably frequent. Experimental work has demonstrated that the process that leads to recombination--a copy choice mechanism involving the migration of reverse transcriptase between viral RNA templates--occurs several times on average during every round of HIV-1 DNA synthesis. Key biological factors that lead to high recombination rates for all retroviruses are the recombination-prone nature of their reverse transcription machinery and their pseudodiploid RNA genomes. However, HIV-1 genes recombine even more frequently than do those of many other retroviruses. This reflects the way in which HIV-1 selects genomic RNAs for coencapsidation as well as cell-to-cell transmission properties that lead to unusually frequent associations between distinct viral genotypes. HIV-1 faces strong and changeable selective conditions during replication within patients. The mode of HIV-1 persistence as integrated proviruses and strong selection for defective proviruses in vivo provide conditions for archiving alleles, which can be resuscitated years after initial provirus establishment. Recombination can facilitate drug resistance and may allow superinfecting HIV-1 strains to evade preexisting immune responses, thus adding to challenges in vaccine development. These properties converge to provide HIV-1 with the means, motive, and opportunity to recombine its genetic material at an unprecedented high rate and to allow genetic recombination to serve as one of the highest barriers to HIV-1 eradication.

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.

The A-rich RNA sequences of HIV-1 pol are important for the synthesis of viral cDNA

The bias of A-rich codons in HIV-1 pol is thought to be a record of hypermutations in viral genomes that lack biological functions. Bioinformatic analysis predicted that A-rich sequences are generally associated with minimal local RNA structures. Using codon modifications to reduce the amount of A-rich sequences within HIV-1 genomes, we have reduced the flexibility of RNA sequences in pol to analyze the functional significance of these A-rich 'structurally poor' RNA elements in HIV-1 pol. Our data showed that codon modification of HIV-1 sequences led to a suppression of virus infectivity by 5-100-fold, and this defect does not correlate with, viral entry, viral protein expression levels, viral protein profiles or virion packaging of genomic RNA. Codon modification of HIV-1 pol correlated with an enhanced dimer stability of the viral RNA genome, which was associated with a reduction of viral cDNA synthesis both during HIV-1 infection and in a cell free reverse transcription assay. Our data provided direct evidence that the HIV-1 A-rich pol sequence is not merely an evolutionary artifact of enzyme-induced hypermutations, and that HIV-1 has adapted to rely on A-rich RNA sequences to support the synthesis of viral cDNA during reverse transcription, highlighting the utility of using 'structurally poor' RNA domains in regulating biological process.

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.

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.

APOBEC3 inhibits mouse mammary tumour virus replication in vivo

Genomes of all mammals encode apobec3 genes, which are thought to have a function in intrinsic cellular immunity to several viruses including human immunodeficiency virus type 1 (HIV-1). APOBEC3 (A3) proteins are packaged into virions and inhibit retroviral replication in newly infected cells, at least in part by deaminating cytidines on the negative strand DNA intermediates. However, the role of A3 in innate resistance to mouse retroviruses is not understood. Here we show that A3 functions during retroviral infection in vivo and provides partial protection to mice against infection with mouse mammary tumour virus (MMTV). Both mouse A3 and human A3G proteins interacted with the MMTV nucleocapsid in an RNA-dependent fashion and were packaged into virions. In addition, mouse A3-containing and human A3G-containing virions showed a marked decrease in titre. Last, A3(-/-) mice were more susceptible to MMTV infection, because virus spread was more rapid and extensive than in their wild-type littermates.

Potential evolutionary influences on overlapping reading frames in the bovine leukemia virus pXBL region

Bovine leukemia virus contains a pXBL region encoding the 3' parts of four regulatory proteins (Tax, Rex, G4, R3) in overlapping reading frames. Here we report the pXBL polymorphisms of 30 isolates from four countries. Rates of overall and synonymous substitutions were consistently lower, and nucleotide/amino acid composition bias and codon bias higher, in more-overlapped than in less-overlapped regions. Ratios of nonsynonymous/synonymous substitutions were lowest in the tax gene and its subregions. The 5' parts of the four genes showed selection patterns corresponding to their genomic context outside of the pXBL region. Longer G4 variants due to a natural stop codon mutation had additional triple overlap with reduced sequence variability. These data support the concept that a higher level of overlapping in coding regions correlates with greater evolutionary constraint. Tax, the most conserved among the four regulatory proteins, showed purifying selection consistent with its importance in the viral life cycle.

Evolution of CXCR4-using human immunodeficiency virus type 1 SF162 is associated with two unique envelope mutations

CCR5-using human immunodeficiency virus type 1 (HIV-1) isolates typically gain CXCR4 use via multiple mutations in V3 and often V1/V2 regions of envelope, and patterns of mutations are distinct for each isolate. Here, we report that multiple CXCR4-using variants of a parental CCR5-using HIV-1 isolate, SF162, obtained by either target cell selection or CCR5 inhibition have a common mutation pattern characterized by the same two V3 mutations and that these mutations preexisted in some of the SF162 stocks. These results imply that SF162 has a single pathway for acquiring CXCR4 use and that prolonged culture is sufficient to select for R5X4 variants.

Host-related nucleotide composition and codon usage as driving forces in the recent evolution of the Astroviridae

The evolutionary history of the Astroviridae comprises the ancient separation between avian and mammalian astrovirus lineages followed by diversification among mammalian astroviruses. The latter process included several cross-species transmissions. We found that the recent, but not the ancient, evolution of astroviruses was associated with a switch in nucleotide composition and codon usage among non-human mammalian versus human/avian astroviruses. Virus and hosts phylogenies based on codon usage agreed with each other and matched the hosts' evolutionary emergence order. This recent switch in driving forces acting at the synonymous level points to the adaptation of codon usage by viruses to that of their hosts after cross-species transmissions. This is the first demonstration of nucleotide composition and codon usage being active driving forces during the recent evolutionary history of a virus group in the host-parasite system.

Tax & rex: overlapping genes of the Deltaretrovirus group

Bovine leukemia virus and human T-cell leukemia viruses I and II, members of the Deltaretrovirus group, have two regulatory genes, tax and rex, that are coded in overlapping reading frames. We found that sequence variations in the rex gene of each virus result in amino acid differences significantly more often than variations in the tax gene. For all three viruses the highest ratio of non-synonymous to synonymous changes was found in the rex gene. In the overlapping regions of tax and rex, the second codon position of Rex corresponds to the third codon position of Tax. Nucleotide C was present in all genes of the three viruses at the highest frequency and this bias was most pronounced in the rex gene. More specifically we found that the C bias and nucleotide variation is greatest at the second codon position of Rex and the third codon position of Tax in the area of tax/rex overlap. Changes in the second codon position of Rex always resulted in amino acid change whereas changes in the third codon position of Tax resulted in amino acid changes less than a third of the time. Analysis of the amino acid frequencies in both proteins shows that there is a disproportionately large percentage of the amino acids alanine, proline, serine and threonine (the four amino acids whose second codon position is C) in Rex. These findings led us to hypothesize that the Rex protein can withstand more amino acid changes than can the Tax protein suggesting that the Tax protein experiences higher evolutionary constraints and is the more conserved of the two proteins.

Factors influencing synonymous codon and amino acid usage biases in Mimivirus

Synonymous codon and amino acid usage biases have been investigated in 903 Mimivirus protein-coding genes in order to understand the architecture and evolution of Mimivirus genome. As expected for an AT-rich genome, third codon positions of the synonymous codons of Mimivirus carry mostly A or T bases. It was found that codon usage bias in Mimivirus genes is dictated both by mutational pressure and translational selection. Evidences show that four factors such as mean molecular weight (MMW), hydropathy, aromaticity and cysteine content are mostly responsible for the variation of amino acid usage in Mimivirus proteins. Based on our observation, we suggest that genes involved in translation, DNA repair, protein folding, etc., have been laterally transferred to Mimivirus a long ago from living organism and with time these genes acquire the codon usage pattern of other Mimivirus genes under selection pressure.

Evidence for the acquisition of multi-drug resistance in an HIV-1 clinical isolate via human sequence transduction

Insertions in HIV-1 reverse transcriptase's fingers subdomain can enhance chain terminator excision and confer resistance to multiple nucleoside analogs. Inserts that resemble flanking sequences likely arise by local sequence duplication. However, a remarkable variety of non-repeat fingers insertions have been observed. Here, molecular epidemiology, sequence analyses and mechanistic modeling were employed to show that one Japanese isolate's RT fingers insert likely resulted from non-homologous recombination between virus and host sequences and the transductive copying of 37 nucleotides from human chromosome 17. These findings provide evidence that human sequence transduction can, at least rarely, contribute to genetic and phenotypic variation in pandemic HIV.

Synonymous codon usage in adenoviruses: influence of mutation, selection and protein hydropathy

Trends in synonymous codon usage in adenoviruses have been examined through the multivariate statistical analysis on the annotated protein-coding regions of 22 adenoviral species, for which complete genome sequences are available. One of the major determinants of such trends is the G+C content at third codon positions of the genes, the average value of which varied from one viral genome to other depending on the overall mutational bias of the species. G3S and C3S interacted synergistically along the first principal axis of correspondence analysis on the Relative Synonymous Codon Usage of adenoviral genes, but antagonistically along the second principal axis. The intra-genomic variation in codon usage pattern in adenoviruses is generally influenced by asymmetrical mutational bias in two DNA strands. Other major determinants of the trends are the natural selection, putatively operative at the level of translation and quite interestingly, hydropathy of the encoded proteins. The trends in codon usage, though characterized by distinct virus-specific mutational bias, do not exhibit any sign of host-specificity. Significant variations are observed in synonymous codon choice in structural and nonstructural genes of adenoviruses.

Use of endogenous retroviral sequences (ERVs) and structural markers for retroviral phylogenetic inference and taxonomy

BACKGROUND

Endogenous retroviral sequences (ERVs) are integral parts of most eukaryotic genomes and vastly outnumber exogenous retroviruses (XRVs). ERVs with a relatively complete structure were retrieved from the genetic archives of humans and chickens, diametrically opposite representatives of vertebrate retroviruses (over 3300 proviruses), and analyzed, using a bioinformatic program, RetroTector, developed by us. This rich source of proviral information, accumulated in a local database, and a collection of XRV sequences from the literature, allowed the reconstruction of a Pol based phylogenetic tree, more extensive than previously possible. The aim was to find traits useful for classification and evolutionary studies of retroviruses. Some of these traits have been used by others, but they are here tested in a wider context than before.

RESULTS

In the ERV collection we found sequences similar to the XRV-based genera: alpha-, beta-, gamma-, epsilon- and spumaretroviruses. However, the occurrence of intermediates between them indicated an evolutionary continuum and suggested that taxonomic changes eventually will be necessary. No delta or lentivirus representatives were found among ERVs. Classification based on Pol similarity is congruent with a number of structural traits. Acquisition of dUTPase occurred three times in retroviral evolution. Loss of one or two NC zinc fingers appears to have occurred several times during evolution. Nucleotide biases have been described earlier for lenti-, delta- and betaretroviruses and were here confirmed in a larger context.

CONCLUSION

Pol similarities and other structural traits contribute to a better understanding of retroviral phylogeny. "Global" genomic properties useful in phylogenies are i.) translational strategy, ii.) number of Gag NC zinc finger motifs, iii.) presence of Pro N-terminal dUTPase (dUTPasePro), iv.) presence of Pro C-terminal G-patch and v.) presence of a GPY/F motif in the Pol integrase (IN) C-terminal domain. "Local" retroviral genomic properties useful for delineation of lower level taxa are i.) host species range, ii.) nucleotide compositional bias and iii.) LTR lengths.

Natural variation in Vif: differential impact on APOBEC3G/3F and a potential role in HIV-1 diversification

The HIV-1 Vif protein counteracts the antiviral activity exhibited by the host cytidine deaminases APOBEC3G and APOBEC3F. Here, we show that defective vif alleles can readily be found in HIV-1 isolates and infected patients. Single residue changes in the Vif protein sequence are sufficient to cause the loss of Vif-induced APOBEC3 neutralization. Interestingly, not all the detected defects lead to a complete inactivation of Vif function since some mutants retained selective neutralizing activity against APOBEC3F but not APOBEC3G or vice versa. Concordantly, independently hypermutated proviruses with distinguishable patterns of G-to-A substitution attributable to cytidine deamination induced by APOBEC3G, APOBEC3F, or both enzymes were present in individuals carrying proviruses with completely or partly defective Vif variants. Natural variation in Vif function may result in selective and partial neutralization of cytidine deaminases and thereby promote viral sequence diversification within HIV-1 infected individuals.