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

Sequence variability, gene structure, and expression of full-length human endogenous retrovirus H

Recently, we identified and classified 926 human endogenous retrovirus H (HERV-H)-like proviruses in the human genome. In this paper, we used the information to, in silico, reconstruct a putative ancestral HERV-H. A calculated consensus sequence was nearly open in all genes. A few manual adjustments resulted in a putative 9-kb HERV-H provirus with open reading frames (ORFs) in gag, pro, pol, and env. Long terminal repeats (LTRs) differed by 1.1%, indicating proximity to an integration event. The gag ORF was extended upstream of the normal myristylation start site. There was a long leader (including a "pre-gag" ORF) region positioned like the N terminus of murine leukemia virus (MLV) "glyco-Gag," potentially encoding a proline- and serine-rich domain remotely similar to MLV pp12. Another ORF, starting inside the 5' LTR, had no obvious similarity to known protein domains. Unlike other hitherto described gammaretroviruses, the reconstructed Gag had two zinc finger motifs. Alternative splicing of sequences related to the HERV-H consensus was confirmed using dbEST data. env transcripts were most prevalent in colon tumors, but also in normal testis. We found no evidence for full length env transcripts in the dbEST. HERV-H had a markedly skewed nucleotide composition, disfavoring guanine and favoring cytidine. We conclude that the HERV-H consensus shared a gene arrangement common to gammaretroviruses with gag separated by stop codon from pro-pol in the same reading frame, while env resides in another reading frame. There was also alternative splicing. HERV-H consensus yielded new insights in gammaretroviral evolution and will be useful as a model in studies on expression and function.

Evolution and distribution of class II-related endogenous retroviruses

Endogenous retroviruses (ERVs) are widespread in vertebrate genomes and have been loosely grouped into "classes" on the basis of their phylogenetic relatedness to the established genera of exogenous retroviruses. Four of these genera-the lentiviruses, alpharetroviruses, betaretroviruses, and deltaretroviruses-form a well-supported clade in retroviral phylogenies, and ERVs that group with these genera have been termed class II ERVs. We used PCR amplification and sequencing of retroviral fragments from more than 130 vertebrate taxa to investigate the evolution of the class II retroviruses in detail. We confirm that class II retroviruses are largely confined to mammalian and avian hosts and provide evidence for a major novel group of avian retroviruses, and we identify additional members of both the alpha- and the betaretrovirus genera. Phylogenetic analyses demonstrated that the avian and mammalian viruses form distinct monophyletic groups, implying that interclass transmission has occurred only rarely during the evolution of the class II retroviruses. In contrast to previous reports, the lentiviruses clustered as sister taxa to several endogenous retroviruses derived from rodents and insectivores. This topology was further supported by the shared loss of both the class II PR-Pol frameshift site and the class II retrovirus G-patch domain.

New insights in HTLV-I phylogeny by sequencing and analyzing the entire envelope gene

The HTLV-I envelope plays a major role in the process of target cell infection. It is implied in the recognition of the viral receptor(s), penetration of the viral genetic material, and induction of host immunity to the virus. It is thus important to study the genetic variability of the viral env gene as well as its variation in terms of evolution. In a new approach to these features, we sequenced the entire env gene of 65 HTLV-I isolates originating from Gabon, French Guiana, West Indies, and Iran, such isolates representing all major HTLVI phylums but the Australo-Melanesian one. The sequences obtained and all PTLV-I (HTLV-I and STLV-I) env sequences available in the literature were analyzed. Phylogenetic studies using different algorithms (minimum evolution, neighbor joining, maximum parsimony, and maximum likelihood) gave the same clear-cut results. Newly sequenced HTLV-I isolates described in this report allocated in three well-defined subtypes: Cosmopolitan, Central African, and a new distinct one that we termed "Maroni" subtype (present in the Maroni Basin, French Guiana, and West Indies). Clearly, the most divergent PTLV-I strains present in Asia- Australo-Melanesia as well as African and Asian STLV-I derived from the same node in the phylogenetic tree as isolates of the Central African subtype. In addition, we showed that within each PTLV-I subtype, groups of isolates may be characterized by nonrandom and systematically associated mutations.

Composition bias and genome polarity of RNA viruses

I have observed a relationship between GC content in coding sequences of RNA viruses and their genome polarity. Positive-stranded RNA viruses have significantly higher GC contents than negative-stranded RNA viruses. Coding sequences of all negative-stranded RNA viruses are biased toward high A in coding strands (high T in genomes), while two distinct patterns were observed among positive-stranded RNA genomes. This finding suggests that RNA viruses with different genome polarity are under different mutational pressure, which may be a consequence of the difference in the strategies of viral genome expression and replication. The GC content directly affects the viral codon adaptation index using highly expressed human genes as the reference set, which may theoretically predict the efficiency of viral gene expression in human cells.

Intrinsic obstacles to human immunodeficiency virus type 1 coreceptor switching

The natural evolution of human immunodeficiency virus type 1 infection often includes a switch in coreceptor preference late in infection from CCR5 to CXCR4, a change associated with expanded target cell range and worsened clinical prognosis. Why coreceptor switching takes so long is puzzling, since it requires as few as one to two mutations. Here we report three obstacles that impede the CCR5-to-CXCR4 switch. Coreceptor switch variants were selected by target cell replacement in vitro. Most switch variants showed diminished replication compared to their parental R5 isolate. Transitional intermediates were more sensitive to both CCR5 and CXCR4 inhibitors than either the parental R5 virus or the final R5X4 (or rare X4) variant. The small number of mutations in viruses selected for CXCR4 use were distinctly nonrandom, with a dominance of charged amino acid substitutions encoded by G-to-A transitions, changes in N-linked glycosylation sites, and isolate-specific mutation patterns. Diminished replication fitness, less-efficient coreceptor use, and unique mutational pathways may explain the long delay from primary infection until the emergence of CXCR4-using viruses.

Recent insights into HIV-1 Vif

The lentiviruses, including HIV-1 (but excluding equine infectious anemia virus), encode a viral infectivity factor (Vif) protein. Circumstantial evidence suggested that Vif acts to neutralize an inhibitory host defense mechanism, but progress in the field was limited because the identity of the cellular target was unknown. The recent identification of the elusive host cell factor let loose a flood of advances. These findings have revealed a novel innate defense mechanism against retroviruses. In infected cells, the cellular cytidine deaminase APOBEC3G, a relative of the activation-induced deaminase (AID), is encapsidated into assembling virions. The enzyme lies in the virion, waiting to wreak havoc on the viral genome in the next round of virus replication--unless it is first caught by Vif.

Mutational patterns correlate with genome organization in SARS and other coronaviruses

Focused efforts by several international laboratories have resulted in the sequencing of the genome of the causative agent of severe acute respiratory syndrome (SARS), novel coronavirus SARS-CoV, in record time. Using cumulative skew diagrams, I found that mutational patterns in the SARS-CoV genome were strikingly different from other coronaviruses in terms of mutation rates, although they were in general agreement with the model of the coronavirus lifecycle. These findings might be relevant for the development of sequence-based diagnostics and the design of agents to treat SARS.

Single-strand specificity of APOBEC3G accounts for minus-strand deamination of the HIV genome

HIV-1 deleted for the vif accessory gene encapsidates the cellular cytidine deaminase APOBEC3G. Upon infection, the encapsidated APOBEC3G induces G-->A mutations in the viral reverse transcripts. The G-->A mutations result either from C-->U deamination of the minus strand or deamination of both strands followed by repair of the plus strand. We report here that minus-strand deamination occurred over the length of the virus genome, preferentially at CCCA sequences, with a graded frequency in the 5'-->3' direction. APOBEC3G induced previously undetected C-->T mutations in the 5' U3 and the primer-binding site, both of which become transiently single-stranded during reverse transcription. In vitro, APOBEC3G bound and deaminated single-stranded DNA (ssDNA) but not double-stranded DNA (dsDNA) or DNA-RNA hybrids. We propose that the requirement for ssDNA accounts for the minus-strand mutations, the 5'-->3' graded frequency of deamination and the rare C-->T mutations.

Contributions of vesicular stomatitis virus to the understanding of RNA virus evolution

Vesicular stomatitis virus has been a preferred system to study evolution for several decades. New approaches to antiviral treatment, such as lethal mutagenesis, stem from investigations done with VSV. Recent work has shed new light in the way we view neutrality, a fundamental concept to understand the evolutionary history of RNA viruses.