Evolutionary mechanisms and population dynamics of the third variable envelope region of HIV within single hosts

Negligible genetic diversity of mycobacterium tuberculosis host immune system protein targets: evidence of limited selective pressure

A common theme in medical microbiology is that the amount of amino acid sequence variation in proteins that are targets of the host immune system greatly exceeds that found in metabolic enzymes or other housekeeping proteins. Twenty-four Mycobacterium tuberculosis genes coding for targets of the host immune system were sequenced in 16 strains representing the breadth of genomic diversity in the species. Of the 24 genes, 19 were invariant and only six polymorphic nucleotide sites were identified in the 5 genes that did have variation. The results document the highly unusual circumstance that prominent M. tuberculosis antigenic proteins have negligible structural variation worldwide. The data are best explained by a combination of three factors: (i) evolutionarily recent global dissemination in humans, (ii) lengthy intracellular quiescence, and (iii) active replication in relatively few fully immunocompetent hosts. The very low level of amino acid diversity in antigenic proteins may be cause for optimism in the difficult fight to control global tuberculosis.

Reevaluation of amino acid variability of the human immunodeficiency virus type 1 gp120 envelope glycoprotein and prediction of new discontinuous epitopes

To elucidate the evolutionary mechanisms of the human immunodeficiency virus type 1 gp120 envelope glycoprotein at the single-site level, the degree of amino acid variation and the numbers of synonymous and nonsynonymous substitutions were examined in 186 nucleotide sequences for gp120 (subtype B). Analyses of amino acid variabilities showed that the level of variability was very different from site to site in both conserved (C1 to C5) and variable (V1 to V5) regions previously assigned. To examine the relative importance of positive and negative selection for each amino acid position, the numbers of synonymous and nonsynonymous substitutions that occurred at each codon position were estimated by taking phylogenetic relationships into account. Among the 414 codon positions examined, we identified 33 positions where nonsynonymous substitutions were significantly predominant. These positions where positive selection may be operating, which we call putative positive selection (PS) sites, were found not only in the variable loops but also in the conserved regions (C1 to C4). In particular, we found seven PS sites at the surface positions of the alpha-helix (positions 335 to 347 in the C3 region) in the opposite face for CD4 binding. Furthermore, two PS sites in the C2 region and four PS sites in the C4 region were detected in the same face of the protein. The PS sites found in the C2, C3, and C4 regions were separated in the amino acid sequence but close together in the three-dimensional structure. This observation suggests the existence of discontinuous epitopes in the protein's surface including this alpha-helix, although the antigenicity of this area has not been reported yet.

Codon-Substitution Models for Heterogeneous Selection Pressure at Amino Acid Sites

Comparison of relative fixation rates of synonymous (silent) and nonsynonymous (amino acid-altering) mutations provides a means for understanding the mechanisms of molecular sequence evolution. The nonsynonymous/synonymous rate ratio (ω = dN/dS) is an important indicator of selective pressure at the protein level, with ω = 1 meaning neutral mutations, ω < 1 purifying selection, and ω > 1 diversifying positive selection. Amino acid sites in a protein are expected to be under different selective pressures and have different underlying ω ratios. We develop models that account for heterogeneous ω ratios among amino acid sites and apply them to phylogenetic analyses of protein-coding DNA sequences. These models are useful for testing for adaptive molecular evolution and identifying amino acid sites under diversifying selection. Ten data sets of genes from nuclear, mitochondrial, and viral genomes are analyzed to estimate the distributions of ω among sites. In all data sets analyzed, the selective pressure indicated by the ω ratio is found to be highly heterogeneous among sites. Previously unsuspected Darwinian selection is detected in several genes in which the average ω ratio across sites is <1, but in which some sites are clearly under diversifying selection with ω > 1. Genes undergoing positive selection include the β-globin gene from vertebrates, mitochondrial protein-coding genes from hominoids, the hemagglutinin (HA) gene from human influenza virus A, and HIV-1 env, vif, and pol genes. Tests for the presence of positively selected sites and their subsequent identification appear quite robust to the specific distributional form assumed for ω and can be achieved using any of several models we implement. However, we encountered difficulties in estimating the precise distribution of ω among sites from real data sets.

Intrapatient variability of HIV type 1 group O ANT70 during a 10-year follow-up

HIV-1 ANT70 is the first HIV-1 group O virus isolate obtained from a 25-year-old Cameroonian woman, who seroconverted in March 1987. This individual has remained asymptomatic and clinically healthy (clinical stage WHO 1, CDC II) even though she did not receive any antiretroviral therapy for HIV-1 before 97 months post-seroconversion. CD4+ T cell counts declined steadily to 200/microl at 70 months postseroconversion. The HIV-1 ANT70 nucleotide and amino acid sequence diversity of the V3C3-encoding env fragment within this individual was followed over a 10-year period. RT-PCR, cloning, sequencing, and genetic analyses were performed on eight plasma follow-up samples. Extensive increasing intra- and intersample variation was observed. This is the first long-term (>10 years) follow-up of the genetic variability of an HIV-1 group O-infected individual. As the course of the disease in the HIV-1 ANT70-infected woman was similar in many aspects to that of group M-infected individuals, it remains to be elucidated whether the changes observed in the V3 loop are critical for disease progression.

Host-specific modulation of the selective constraints driving human immunodeficiency virus type 1 env gene evolution

To address the evolution of human immunodeficiency virus type 1 (HIV-1) within a single host, we analyzed the HIV-1 C2-V5 env regions of both cell-free genomic-RNA- and proviral-DNA-derived clones. Sequential samples were collected over a period of 3 years from six untreated subjects (three typical progressors [TPs] and three slow progressors [SPs], all with a comparable length of infection except one. The evolutionary analysis of the C2-V5 env sequences performed on 506 molecular clones (253 RNA- and 253 DNA-derived sequences) highlighted a series of differences between TPs and SPs. In particular, (i) clonal sequences from SPs (DNA and RNA) showed lower nucleotide similarity than those from TPs (P = 0. 0001), (ii) DNA clones from SPs showed higher intra- and intersample nucleotide divergence than those from TPs (P < 0.05), (iii) higher host-selective pressure was generally detectable in SPs (DNA and RNA sequences), and (iv) the increase in the genetic distance of DNA and RNA sequences over time was paralleled by an increase in both synonymous (Ks) and nonsynonymous (Ka) substitutions in TPs but only in nonsynonymous substitutions in SPs. Several individual peculiarities of the HIV-1 evolutionary dynamics emerged when the V3, V4, and V5 env regions of both TPs and SPs were evaluated separately. These peculiarities, probably reflecting host-specific features of selective constraints and their continuous modulation, are documented by the dynamics of Ka/Ks ratios of hypervariable env domains.

Evolution of envelope sequences from the genital tract and peripheral blood of women infected with clade A human immunodeficiency virus type 1

The development of viral diversity during the course of human immunodeficiency virus type 1 (HIV-1) infection may significantly influence viral pathogenesis. The paradigm for HIV-1 evolution is based primarily on studies of male cohorts in which individuals were presumably infected with a single virus variant of subtype B HIV-1. In this study, we evaluated virus evolution based on sequence information of the V1, V2, and V3 portions of HIV-1 clade A envelope genes obtained from peripheral blood and cervical secretions of three women with genetically heterogeneous viral populations near seroconversion. At the first sample following seroconversion, the number of nonsynonymous substitutions per potential nonsynonymous site (dn) significantly exceeded substitutions at potential synonymous sites (ds) in plasma viral sequences from all individuals. Generally, values of dn remained higher than values of ds as sequences from blood or mucosa evolved. Mutations affected each of the three variable regions of the envelope gene differently; insertions and deletions dominated changes in V1, substitutions involving charged amino acids occurred in V2, and sequential replacement of amino acids over time at a small subset of positions distinguished V3. The relationship among envelope nucleotide sequences obtained from peripheral blood mononuclear cells, plasma, and cervical secretions was evaluated for each individual by both phylogenetic and phenetic analyses. In all subjects, sequences from within each tissue compartment were more closely related to each other than to sequences from other tissues (phylogenetic tissue compartmentalization). At time points after seroconversion in two individuals, there was also greater genetic identity among sequences from the same tissue compartment than among sequences from different tissue compartments (phenetic tissue compartmentalization). Over time, temporal phylogenetic and phenetic structure was detectable in mucosal and plasma viral samples from all three women, suggesting a continual process of migration of one or a few infected cells into each compartment followed by localized expansion and evolution of that population.

CC-chemokines enhance the replication of T-tropic strains of HIV-1 in CD4(+) T cells: role of signal transduction

This study demonstrates that several CC-chemokines, including those that inhibit entry and replication of macrophage-tropic strains of HIV, increase the replication of T cell (T)-tropic strains in CD4(+) T cells. Enhancement of T-tropic HIV replication is observed at early stages of replication, requires signaling through inhibitory guanine nucleotide-binding regulatory (Gi) proteins, and is associated with increased cell surface colocalization of CD4 and the T-tropic HIV coreceptor CXCR4. These findings may further our understanding of the factors that influence the replication and spread of T-tropic strains of HIV in vivo and suggest that the use of cell signaling CC-chemokines as therapeutic agents for the purpose of limiting HIV replication in vivo should be approached with caution.

High frequency of the GWG (Pro Trp) envelope variant of HIV-1 in Southeast Brazil

The HIV-1 variant that contains the GWG amino acid sequence in the crown of the principal neutralizing determinant (PND) has been detected in a few patients in Japan, France, and Brazil by direct sequencing. We describe for the first time the use of restriction fragment length polymorphism (RFLP) and limited DNA sequencing of the C2-V3 region of the HIV-1 envelope (env) gene to determine the prevalence of the variant in 75 HIV-1-infected Brazilian patients. Overall prevalence of the GWG sequence as indicated by RFLP was 57% (43 of 75). The prevalence in females (72%) was higher than that in males (32%) and newborns (40%). Two GFG sequences and 1 GLG sequence were also detected. This finding is relevant for the planning of vaccines and for studies of the epidemiology of HIV-1 in Brazil.

Likelihood models for detecting positively selected amino acid sites and applications to the HIV-1 envelope gene

Several codon-based models for the evolution of protein-coding DNA sequences are developed that account for varying selection intensity among amino acid sites. The "neutral model" assumes two categories of sites at which amino acid replacements are either neutral or deleterious. The "positive-selection model" assumes an additional category of positively selected sites at which nonsynonymous substitutions occur at a higher rate than synonymous ones. This model is also used to identify target sites for positive selection. The models are applied to a data set of the V3 region of the HIV-1 envelope gene, sequenced at different years after the infection of one patient. The results provide strong support for variable selection intensity among amino acid sites The neutral model is rejected in favor of the positive-selection model, indicating the operation of positive selection in the region. Positively selected sites are found in both the V3 region and the flanking regions.

Variation in the glycoprotein and VP35 genes of Marburg virus strains

Marburg virus, the prototype of the family Filoviridae, differs genetically, serologically, and morphologically from Ebola viruses. To better define the genetic variation within the species, VP35 and glycoprotein (GP) genes of representative human isolates from four known episodes of Marburg virus hemorrhagic fever were analyzed. The percentage nucleotide differences in the GP gene coding regions of Marburg viruses (0.1-21%) was nearly equal to the percentage amino acid changes (0-23%), while the percentage nucleotide differences in VP35 coding regions (0.3-20.9%) were higher than the percentage amino acid changes (0.9-6.1%), indicating a greater number of nonsynonymous changes occurring in the GP gene. The higher variation in the GP gene and the corresponding protein, especially those changes in the variable middle region of the GP, suggests that the variability may be the result of responses to natural host pressures. Analysis of the GP gene open reading frame shows a nonrandom distribution of nonsynonymous mutations that may indicate positive Darwinian selection is operating within the variable region. A heptad repeat region and an adjoining predicted fusion peptide are found in the C-terminal third of Marburg virus GPs, as has been previously shown for Ebola virus, and are similar to those found in transmembrane glycoproteins of retroviruses, paramyxoviruses, coronaviruses, and influenza viruses. Comparative analyses showed that there are two lineages within the Marburg virus species of filoviruses. The most recent isolate from Kenya (1987) represents a separate genetic lineage within the Marburg virus species (21-23% amino acid difference). However, this lineage likely does not represent a separate Marburg subtype, as the extent of divergence is less than that separating Ebola virus subtypes.

Co-evolution of human immunodeficiency virus and cytotoxic T-lymphocyte responses

After more than a decade of intensive research, the precise role of human immunodeficiency virus (HIV)-specific cytotoxic T lymphocytes (CTL) in determining the course of the infection remains open to argument. It is established that HIV-specific CTL appear early in the infection and are temporally associated with the clearance of culturable virus from the blood; that CTL are generally detectable at very high levels throughout the asymptomatic phase and decline at the time of progression to AIDS; and that CTL-mediated killing is sufficiently fast to prevent production of new virions by HIV-infected cells. However, viral turnover is high throughout the course of the infection, and infected individuals progress inexorably to disease in spite of the CTL response. In order to address the question of whether CTL play an active part in influencing the course of HIV infection, one approach has been to seek evidence for CTL-mediated selection pressure on the virus. Several clear examples of CTL epitope-specific mutations selected to fixation are described. We argue that CTL escape is a common event which occurs at all stages of the infection. Detailed longitudinal studies are required to detect CTL escape and to understand the complexities contributed by factors such as a polyvalent CTL response and the presence of epitope variants which antagonise the CTL response. In conclusion, there is strong evidence of a dynamic process in which CTL impose important selection constraints upon HIV from which the virus attempts to escape; ultimately, at the time of disease progression, the tenuous control of CTL over the virus is lost.