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

Subtle differences in selective pressures applied on the envelope gene of HIV-1 in pregnant versus non-pregnant women

Pregnancy is associated with modulations of maternal immunity that contribute to foeto-maternal tolerance. To understand whether and how these alterations impact antiviral immunity, a detailed cross-sectional analysis of selective pressures exerted on HIV-1 envelope amino-acid sequences was performed in a group of pregnant (n = 32) and non-pregnant (n = 44) HIV-infected women in absence of treatment with antiretroviral therapy (ART). Independent of HIV-1 subtype, p-distance, dN and dS were all strongly correlated with one another but were not significantly different in pregnant as compared to non-pregnant patients. Differential levels of selective pressure applied on different Env subdomains displayed similar yet non-identical patterns between the two groups, with pressure applied on C1 being significantly lower in constant regions C1 and C2 than in V1, V2, V3 and C3. To draw a general picture of the selection applied on the envelope and compensate for inter-individual variations, we performed a binomial test on selection frequency data pooled from pregnant and non-pregnant women. This analysis uncovered 42 positions, present in both groups, exhibiting statistically-significant frequency of selection that invariably mapped to the surface of the Env protein, with the great majority located within epitopes recognized by Env-specific antibodies or sites associated with the development of cross-reactive neutralizing activity. The median frequency of occurrence of positive selection per site was significantly lower in pregnant versus non-pregnant women. Furthermore, examination of the distribution of positively selected sites using a hypergeometric test revealed that only 2 positions (D137 and S142) significantly differed between the 2 groups. Taken together, these result indicate that pregnancy is associated with subtle yet distinctive changes in selective pressures exerted on the HIV-1 Env protein that are compatible with transient modulations of maternal immunity.

Identification of Owl Monkey CD4 Receptors Broadly Compatible with Early-Stage HIV-1 Isolates

Most HIV-1 variants isolated from early-stage human infections do not use nonhuman primate versions of the CD4 receptor for cellular entry, or they do so poorly. We and others have previously shown that CD4 has experienced strong natural selection over the course of primate speciation, but it is unclear whether this selection has influenced the functional characteristics of CD4 as an HIV-1 receptor. Surprisingly, we find that selection on CD4 has been most intense in the New World monkeys, animals that have never been found to harbor lentiviruses related to HIV-1. Based on this, we sampled CD4 genetic diversity within populations of individuals from seven different species, including five species of New World monkeys. We found that some, but not all, CD4 alleles found in Spix's owl monkeys (Aotus vociferans) encode functional receptors for early-stage human HIV-1 isolates representing all of the major group M clades (A, B, C, and D). However, only some isolates of HIV-1 subtype C can use the CD4 receptor encoded by permissive Spix's owl monkey alleles. We characterized the prevalence of functional CD4 alleles in a colony of captive Spix's owl monkeys and found that 88% of surveyed individuals are homozygous for permissive CD4 alleles, which encode an asparagine at position 39 of the receptor. We found that the CD4 receptors encoded by two other species of owl monkeys (Aotus azarae and Aotus nancymaae) also serve as functional entry receptors for early-stage isolates of HIV-1.

IMPORTANCE

Nonhuman primates, particularly macaques, are used for preclinical evaluation of HIV-1 vaccine candidates. However, a significant limitation of the macaque model is the fact that most circulating HIV-1 variants cannot use the macaque CD4 receptor to enter cells and have to be adapted to these species. This is particularly true for viral variants from early stages of infection, which represent the most relevant vaccine targets. In this study, we found that some individuals from captive owl monkey populations harbor CD4 alleles that are compatible with a broad collection of HIV-1 isolates, including those isolated from early in infection in highly affected populations and representing diverse subtypes.

Pattern of the evolution of HIV-1 enν gene in Côte d׳Ivoire

Cête d׳Ivoire continues to have the highest HIV-1 prevalence rate in West Africa, although the infection number is in constant decline. The external envelope protein of the viruses is a likely site of selection, and responsible for receptor binding and entry into host cells, and therefore constitutes an ideal region with which to investigate the evolutionary processes acting on HIV-1. In this study, we analyse 189 envelope glycoprotein V3 loop region sequences of viruse isolates from 1995 to 2009, from HIV-1 untreated patients living in Cête d׳Ivoire, to decipher the temporal relationship between disease diversity, divergence and selection. Our analyses show that the nonsynonymous and synonymous ratio (dN/dS) was lower than 1 for viral populations analysed within 15 years, which showed the sequences did not undergo adequate immune pressure. The phylogenetic tree of the sequences analysed demonstrated distinctly long internal branches and short external branches, suggesting that only a small number of viruses infected the new host cell at each transmission. In addition to identifying sites under purifying selection, we also identified neutral sites that can cause false positive inference of selection. These sites presented form a resource for future studies of selection pressures acting on HIV-1 enν gene in Cête d׳Ivoire and other West African countries.

Molecular characterization of HIV-1 subtype C gp-120 regions potentially involved in virus adaptive mechanisms

The role of variable regions of HIV-1 gp120 in immune escape of HIV has been investigated. However, there is scant information on how conserved gp120 regions contribute to virus escaping. Here we have studied how molecular sequence characteristics of conserved C3, C4 and V3 regions of clade C HIV-1 gp120 that are involved in HIV entry and are target of the immune response, are modulated during the disease course. We found an increase of “shifting” putative N-glycosylation sites (PNGSs) in the α2 helix (in C3) and in C4 and an increase of sites under positive selection pressure in the α2 helix during the chronic stage of disease. These sites are close to CD4 and to co-receptor binding sites. We also found a negative correlation between electric charges of C3 and V4 during the late stage of disease counteracted by a positive correlation of electric charges of α2 helix and V5 during the same stage. These data allow us to hypothesize possible mechanisms of virus escape involving constant and variable regions of gp120. In particular, new mutations, including new PNGSs occurring near the CD4 and CCR5 binding sites could potentially affect receptor binding affinity and shield the virus from the immune response.

Computational prediction of broadly neutralizing HIV-1 antibody epitopes from neutralization activity data

Broadly neutralizing monoclonal antibodies effective against the majority of circulating isolates of HIV-1 have been isolated from a small number of infected individuals. Definition of the conformational epitopes on the HIV spike to which these antibodies bind is of great value in defining targets for vaccine and drug design. Drawing on techniques from compressed sensing and information theory, we developed a computational methodology to predict key residues constituting the conformational epitopes on the viral spike from cross-clade neutralization activity data. Our approach does not require the availability of structural information for either the antibody or antigen. Predictions of the conformational epitopes of ten broadly neutralizing HIV-1 antibodies are shown to be in good agreement with new and existing experimental data. Our findings suggest that our approach offers a means to accelerate epitope identification for diverse pathogenic antigens.

Low complexity regions (LCRs) contribute to the hypervariability of the HIV-1 gp120 protein

Low complexity regions (LCRs) are sequences of nucleic acids or proteins defined by a compositional bias. Their occurrence has been confirmed in sequences of the three cellular lineages (Bacteria, Archaea and Eucarya), and has also been reported in viral genomes. We present here the results of a detailed computer analysis of the LCRs present in the HIV-1 glycoprotein 120 (gp120) encoded by the viral gene env. The analysis was performed using a sample of 3637 Env polyprotein sequences derived from 4117 completely sequenced and translated HIV-1 genomes available in public databases as of December 2012. We have identified 1229 LCRs located in four different regions of the gp120 protein that correspond to four of the five regions that have been identified as hypervariable (V1, V2, V4 and V5). The remaining 29 LCRs are found in the signal peptide and in the conserved regions C2, C3, C4 and C5. No LCR has been identified in the hypervariable region V3. The LCRs detected in the V1, V2, V4, and V5 hypervariable regions exhibit a high Asn content in their amino acid composition, which very likely correspond to glycosylation sites, which may contribute to the retroviral ability to avoid the immune system. In sharp contrast with what is observed in gp120 proteins lacking LCRs, the glycosylation sites present in LCRs tend to be clustered towards the center of the region forming well-defined islands. The results presented here suggest that LCRs represent a hitherto undescribed source of genomic variability in lentivirus, and that these repeats may represent an important source of antigenic variation in HIV-1 populations. The results reported here may exemplify the evolutionary processes that may have increased the size of primitive cellular RNA genomes and the role of LCRs as a source of raw material during the processes of evolutionary acquisition of new functions.

Dynamics of envelope evolution in clade C SHIV-infected pig-tailed macaques during disease progression analyzed by ultra-deep pyrosequencing

Understanding the evolution of the human immunodeficiency virus type 1 (HIV-1) envelope during disease progression can provide tremendous insights for vaccine development, and simian-human immunodeficiency virus (SHIV) infection of non-human primate provides an ideal platform for such studies. A newly developed clade C SHIV, SHIV-1157ipd3N4, which was able to infect rhesus macaques, closely resembled primary HIV-1 in transmission and pathogenesis, was used to infect several pig-tailed macaques. One of the infected animals subsequently progressed to AIDS, whereas one remained a non-progressor. The viral envelope evolution in the infected animals during disease progression was analyzed by a bioinformatics approach using ultra-deep pyrosequencing. Our results showed substantial envelope variations emerging in the progressor animal after the onset of AIDS. These envelope variations impacted the length of the variable loops and charges of different envelope regions. Additionally, multiple mutations were located at the CD4 and CCR5 binding sites, potentially affecting receptor binding affinity, viral fitness and they might be selected at late stages of disease. More importantly, these envelope mutations are not random since they had repeatedly been observed in a rhesus macaque and a human infant infected by either SHIV or HIV-1, respectively, carrying the parental envelope of the infectious molecular clone SHIV-1157ipd3N4. Moreover, similar mutations were also observed from other studies on different clades of envelopes regardless of the host species. These recurring mutations in different envelopes suggest that there may be a common evolutionary pattern and selection pathway for the HIV-1 envelope during disease progression.

Genes in the terminal regions of orthopoxvirus genomes experience adaptive molecular evolution

BACKGROUND

Orthopoxviruses are dsDNA viruses with large genomes, some encoding over 200 genes. Genes essential for viral replication are located in the center of the linear genome and genes encoding host response modifiers and other host interacting proteins are located in the terminal regions. The central portion of the genome is highly conserved, both in gene content and sequence, while the terminal regions are more diverse. In this study, we investigated the role of adaptive molecular evolution in poxvirus genes and the selective pressures that act on the different regions of the genome. The relative fixation rates of synonymous and non-synonymous mutations (the d(N)/d(S) ratio) are an indicator of the mechanism of evolution of sequences, and can be used to identify purifying, neutral, or diversifying selection acting on a gene. Like highly conserved residues, amino acids under diversifying selection may be functionally important. Many genes experiencing diversifying selection are involved in host-pathogen interactions, such as antigen-antibody interactions, or the "host-pathogen arms race."

RESULTS

We analyzed 175 gene families from orthopoxviruses for evidence of diversifying selection. 79 genes were identified as experiencing diversifying selection, 25 with high confidence. Many of these genes are located in the terminal regions of the genome and function to modify the host response to infection or are virion-associated, indicating a greater role for diversifying selection in host-interacting genes. Of the 79 genes, 20 are of unknown function, and implicating diversifying selection as an important mechanism in their evolution may help characterize their function or identify important functional residues.

CONCLUSIONS

We conclude that diversifying selection is an important mechanism of orthopoxvirus evolution. Diversifying selection in poxviruses may be the result of interaction with host defense mechanisms.

Change of positive selection pressure on HIV-1 envelope gene inferred by early and recent samples

HIV-1 infection has been on the rise in Japan recently, and the main transmission route has changed from blood transmission in the 1980s to homo- and/or hetero-sexual transmission in the 2000s. The lack of early viral samples with clinical information made it difficult to investigate the possible virological changes over time. In this study, we sequenced 142 full-length env genes collected from 16 Japanese subjects infected with HIV-1 in the 1980s and in the 2000s. We examined the diversity change in sequences and potential adaptive evolution of the virus to the host population. We used a codon-based likelihood method under the branch-site and clade models to detect positive selection operating on the virus. The clade model was extended to account for different positive selection pressures in different viral populations. The result showed that the selection pressure was weaker in the 2000s than in the 1980s, indicating that it might have become easier for the HIV to infect a new host and to develop into AIDS now than 20 years ago and that the HIV may be becoming more virulent in the Japanese population. The study provides useful information on the surveillance of HIV infection and highlights the utility of the extended clade models in analysis of virus populations which may be under different selection pressures.

An introduction to a novel population genetic approach for HIV characterization

The rapid evolution of the HIV genome is influenced in part by host selection pressure, which may cause parallel evolution among strains under shared selection pressures. To understand the mechanisms behind HIV-host immune escape across host populations, researchers have compared signatures of positive selection pressure on HIV codons across HIV subtypes and across phylogenetic groups of isolates within major subtypes, all relying on a criterion of phylogenetic separation. The HIV codon sites that retain diversity, evolve convergently among sets of hosts (cohorts) and diverge between cohorts may be phylogenetically undiagnostic (reveal little information about the relationship of the strains) and thus undetectable on a tree. We propose a new approach to characterizing genetic divergence among isolates using existing population genetic methods to better understand HIV response to host selection pressures. The approach combines population genetic statistical methods with codon analysis to identify putative amino acid sites evolving convergently. To illustrate the approach, we compared the C2-V3-C3 region of the envelope protein of HIV-1 clade B isolates between Haiti and USA hosts. This region showed no phylogenetic separation between host populations. Still, we identified codon sites in the C2-V3-C3 HIV-1 region that may have evolved differently between the two host populations. The sites are localized in human leukocyte antigen (HLA) class I binding epitopes, N-glycosylation motifs or both and are limited to the C2 and C3 regions. Our method provides a potential means to reveal candidate sites actively involved in HIV-1 immune escape that would otherwise be missed if a requisite for phylogenetic distinctiveness was made a priori. This strategy may prove to be a helpful way to characterize HIV genetic variation among hosts with suspected selection pressure differences, like progressors versus non-progressors.

Dynamic features of the selective pressure on the human immunodeficiency virus type 1 (HIV-1) gp120 CD4-binding site in a group of long term non progressor (LTNP) subjects

The characteristics of intra-host human immunodeficiency virus type 1 (HIV-1) env evolution were evaluated in untreated HIV-1-infected subjects with different patterns of disease progression, including 2 normal progressor [NP], and 5 Long term non-progressor [LTNP] patients. High-resolution phylogenetic analysis of the C2-C5 env gene sequences of the replicating HIV-1 was performed in sequential samples collected over a 3–5 year period; overall, 301 HIV-1 genomic RNA sequences were amplified from plasma samples, cloned, sequenced and analyzed. Firstly, the evolutionary rate was calculated separately in the 3 codon positions. In all LTNPs, the 3^rd codon mutation rate was equal or even lower than that observed at the 1^st and 2^nd positions (p = 0.016), thus suggesting strong ongoing positive selection. A Bayesian approach and a maximum-likelihood (ML) method were used to estimate the rate of virus evolution within each subject and to detect positively selected sites respectively. A great number of N-linked glycosylation sites under positive selection were identified in both NP and LTNP subjects. Viral sequences from 4 of the 5 LTNPs showed extensive positive selective pressure on the CD4-binding site (CD4bs). In addition, localized pressure in the area of the IgG-b12 epitope, a broad neutralizing human monoclonal antibody targeting the CD4bs, was documented in one LTNP subject, using a graphic colour grade 3-dimensional visualization. Overall, the data shown here documenting high selective pressure on the HIV-1 CD4bs of a group of LTNP subjects offers important insights for planning novel strategies for the immune control of HIV-1 infection.

The role of the humoral immune response in the molecular evolution of the envelope C2, V3 and C3 regions in chronically HIV-2 infected patients

Background

This study was designed to investigate, for the first time, the short-term molecular evolution of the HIV-2 C2, V3 and C3 envelope regions and its association with the immune response. Clonal sequences of the env C2V3C3 region were obtained from a cohort of eighteen HIV-2 chronically infected patients followed prospectively during 2–4 years. Genetic diversity, divergence, positive selection and glycosylation in the C2V3C3 region were analysed as a function of the number of CD4+ T cells and the anti-C2V3C3 IgG and IgA antibody reactivity

Results

The mean intra-host nucleotide diversity was 2.1% (SD, 1.1%), increasing along the course of infection in most patients. Diversity at the amino acid level was significantly lower for the V3 region and higher for the C2 region. The average divergence rate was 0.014 substitutions/site/year, which is similar to that reported in chronic HIV-1 infection. The number and position of positively selected sites was highly variable, except for codons 267 and 270 in C2 that were under strong and persistent positive selection in most patients. N-glycosylation sites located in C2 and V3 were conserved in all patients along the course of infection. Intra-host variation of C2V3C3-specific IgG response over time was inversely associated with the variation in nucleotide and amino acid diversity of the C2V3C3 region. Variation of the C2V3C3-specific IgA response was inversely associated with variation in the number of N-glycosylation sites.

Conclusion

The evolutionary dynamics of HIV-2 envelope during chronic aviremic infection is similar to HIV-1 implying that the virus should be actively replicating in cellular compartments. Convergent evolution of N-glycosylation in C2 and V3, and the limited diversification of V3, indicates that there are important functional constraints to the potential diversity of the HIV-2 envelope. C2V3C3-specific IgG antibodies are effective at reducing viral population size limiting the number of virus escape mutants. The C3 region seems to be a target for IgA antibodies and increasing N-linked glycosylation may prevent HIV-2 envelope recognition by these antibodies. Our results provide new insights into the biology of HIV-2 and its relation with the human host and may have important implications for vaccine design.

Retrocyclin RC-101 overcomes cationic mutations on the heptad repeat 2 region of HIV-1 gp41

Retrocyclin RC-101, a theta-defensin with lectin-like properties, potently inhibits infection by many HIV-1 subtypes by binding to the heptad repeat 2 (HR2) region of glycoprotein 41 (gp41) and preventing six-helix bundle formation. In the present study, we used in silico computational exploration to identify residues of HR2 that interacted with RC-101, and then analyzed the HIV-1 sequence database at Los Alamos National Laboratory (New Mexico, USA) for residue variations in the heptad repeat 1 (HR1) and HR2 segments that could plausibly impart in vivo resistance. Docking RC-101 to gp41 peptides in silico confirmed its strong preference for HR2 over HR1, and implicated residues crucial for its ability to bind HR2. We mutagenized these residues in pseudotyped HIV-1 JR.FL reporter viruses, and subjected them to single-round replication assays in the presence of 1.25-10 microg x mL(-1) RC-101. Apart from one mutant that was partially resistant to RC-101, the other pseudotyped viruses with single-site cationic mutations in HR2 manifested absent or impaired infectivity or retained wild-type susceptibility to RC-101. Overall, these data suggest that most mutations capable of rendering HIV-1 resistant to RC-101 will also exert deleterious effects on the ability of HIV-1 to initiate infections - an interesting and novel property for a potential topical microbicide.

The quest for natural selection in the age of comparative genomics

Continued genome sequencing has fueled progress in statistical methods for understanding the action of natural selection at the molecular level. This article reviews various statistical techniques (and their applicability) for detecting adaptation events and the functional divergence of proteins. As large-scale automated studies become more frequent, they provide a useful resource for generating biological null hypotheses for further experimental and statistical testing. Furthermore, they shed light on typical patterns of lineage-specific evolution of organisms, on the functional and structural evolution of protein families and on the interplay between the two. More complex models are being developed to better reflect the underlying biological and chemical processes and to complement simpler statistical models. Linking molecular processes to their statistical signatures in genomes can be demanding, and the proper application of statistical models is discussed.

QUASI analysis of the HIV-1 envelope sequences in the Los Alamos National Laboratory HIV sequence database: pattern and distribution of positive selection sites and their frequencies over years

The envelope (env) protein of human immunodeficiency virus type 1 (HIV-1) plays a crucial role in virus entry and is a central target for HIV vaccine design. Using the QUASI program, we analyzed the conserved regions of all currently available env sequences in the Los Alamos National Laboratory HIV Sequence Database and identified positive selection (PS) sites that are likely to be restricted by host immune responses. We found that PS sites are dispersed across conserved regions of env sequence, and that the C3, C4, and C5 regions were the most targeted. Several regions were identified as being PS free and were mainly distributed in the C1 and C2 regions. When comparing individual QUASI PS site frequencies across clades and geographical regions with the overall frequency of the entire env database, the env sequences from North America showed significantly lower PS site frequency, while those from Asia were significantly higher using Student's t test. The QUASI PS site frequency of env proteins from viruses isolated from different years showed that the PS site frequencies of the env population increased over time. Our study provides an overview of PS sites across the conserved regions of HIV-1 env sequences.

Adaptation to human populations is revealed by within-host polymorphisms in HIV-1 and hepatitis C virus

CD8(+) cytotoxic T-lymphocytes (CTLs) perform a critical role in the immune control of viral infections, including those caused by human immunodeficiency virus type 1 (HIV-1) and hepatitis C virus (HCV). As a result, genetic variation at CTL epitopes is strongly influenced by host-specific selection for either escape from the immune response, or reversion due to the replicative costs of escape mutations in the absence of CTL recognition. Under strong CTL-mediated selection, codon positions within epitopes may immediately "toggle" in response to each host, such that genetic variation in the circulating virus population is shaped by rapid adaptation to immune variation in the host population. However, this hypothesis neglects the substantial genetic variation that accumulates in virus populations within hosts. Here, we evaluate this quantity for a large number of HIV-1- (n > or = 3,000) and HCV-infected patients (n > or = 2,600) by screening bulk RT-PCR sequences for sequencing "mixtures" (i.e., ambiguous nucleotides), which act as site-specific markers of genetic variation within each host. We find that nonsynonymous mixtures are abundant and significantly associated with codon positions under host-specific CTL selection, which should deplete within-host variation by driving the fixation of the favored variant. Using a simple model, we demonstrate that this apparently contradictory outcome can be explained by the transmission of unfavorable variants to new hosts before they are removed by selection, which occurs more frequently when selection and transmission occur on similar time scales. Consequently, the circulating virus population is shaped by the transmission rate and the disparity in selection intensities for escape or reversion as much as it is shaped by the immune diversity of the host population, with potentially serious implications for vaccine design.

Genetic and biologic characterization of HIV type 1 subtype C isolates from south Brazil

The molecular and biological properties of HIV-1 subtype C strains from South Brazil were investigated. We sequenced gag and env fragments of viruses from 22 HIV-1-infected individuals from Porto Alegre City, which has the highest frequency of subtype C in the country. The sequences were then compared with other subtype B, C, and F strains isolated in Brazil and other countries using phylogenetic methods. Amino acid signatures were identified and correlated with phenotypic characteristics. We identified six strains with subtype C (27.3%), eight subtype B (36.4%), one subtype F (4.5%), six C/B recombinants (27.3%), and one B/F recombinant (4.5%). The Brazilian subtype C sequences formed a unique phylogenetic group and presented 6 and 18 specific amino acid signatures in gag and env, respectively. Three distinct patterns of C/B recombinants presented characteristic Brazilian amino acid substitutions. Subtype C viruses were predominantly R5 and non-syncytium-inducing, while C/B recombinants were R5/X4 and syncytium-inducing viruses. These findings suggest that subtype C viruses circulating in Brazil are the result of a unique introduction into the country. Recombination events between subtypes B and C have been occurring frequently for more than 10 years in South Brazil. Biological characterization confirms the hypothesis that subtype C is distinct from the others in the evolution of coreceptor utilization.

Bursts of nonsynonymous substitutions in HIV-1 evolution reveal instances of positive selection at conservative protein sites

The fixation of a new allele can be driven by Darwinian positive selection or can be due to random genetic drift. Identifying instances of positive selection is a difficult task, because its impact is routinely obscured by the action of negative selection. The nature of the genetic code dictates that positive selection in favor of an amino acid replacement should often cause a burst of two or three nucleotide substitutions at a single codon site, because a large fraction of amino acid replacements cannot be achieved after just one nucleotide substitution. Here, we study pairs of successive nonsynonymous substitutions at one codon in the course of evolution of HIV-1 genes within HIV-1 populations inhabiting infected individuals. Such pairs are more numerous and more clumped than expected if different substitutions were independent and than what is observed for pairs of successive synonymous substitutions. Bursts of nonsynonymous substitutions in HIV-1 evolution cannot be explained by mutational biases and must, therefore, be due to positive selection. Both reversals, exact or imprecise, of fixed deleterious mutations and acquisitions of amino acids with new properties are responsible for the bursts. Temporal clumping is strongest at codon sites with a low overall rate of nonsynonymous evolution, implying that a substantial fraction of replacements of conservative amino acids are driven by positive selection. We identified many conservative sites of HIV-1 proteins that occasionally experience positive selection.

The HIV Env variant N283 enhances macrophage tropism and is associated with brain infection and dementia

HIV infects tissue macrophages and brain microglia, which express lower levels of CD4 and CCR5 than CD4+ T cells in peripheral blood. Mechanisms that enhance HIV tropism for macrophages in the CNS and other tissues are not well understood. Here, we identify an HIV envelope glycoprotein (Env) variant in the CD4-binding site of gp120, Asn 283 (N283), that is present at a high frequency in brain tissues from AIDS patients with HIV-associated dementia (HAD). N283 increases gp120 affinity for CD4 by decreasing the gp120-CD4 dissociation rate, enhancing the capacity of HIV Envs to use low levels of CD4 for virus entry and increasing viral replication in macrophages and microglia. Structural modeling suggests that the enhanced ability of Envs with N283 to use low levels of CD4 is due to a hydrogen bond formed with Gln 40 of CD4. N283 is significantly more frequent in brain-derived Envs from HAD patients (41%; n=330) compared with non-HAD patients (8%; n=151; P<0.001). These findings suggest that the macrophage-tropic HIV Env variant N283 is associated with brain infection and dementia in vivo, representing an example of a HIV variant associated with a specific AIDS-related complication.

Evolution of BK virus based on complete genome data

The human polyomavirus BK virus (BKV) is ubiquitous in humans, infecting children asymptomatically. BKV is the only primate polyomavirus that has subtypes (I-IV) distinguishable by immunological reactivity. Nucleotide (nt) variations in a major capsid protein (VP1) gene region (designated the epitope region), probably responsible for antigenic diversity, have been used to classify BKV isolates into subtypes. Here, with all the protein-encoding gene sequences, we attempted to elucidate the evolutionary relationships among 28 BKV isolates belonging to subtypes I, III, and IV (no isolate belonging to subtype II, a minor one, was included). First, using the GTR + Gamma + I model, maximum likelihood trees were reconstructed for individual viral genes as well as for concatenated viral genes. On the resultant trees, the 28 BKV isolates were consistently divided into three clades corresponding to subtypes I, III, and IV, although bootstrap probabilities are not always high. Then we used more sophisticated likelihood models, one of which takes account of codon structure, to elucidate the phylogenetic relationships among BKV subtypes, but the phylogeny of the deep branchings remained ambiguous. Furthermore, the possibility of positive selection in the evolution of BKV was examined using the nonsynonymous/synonymous rate ratio as a measure of selection. An analysis based on entire genes could not detect any strong evidence for positive selection, but that based on the epitope region identified a few sites potentially under positive selection (these sites were among those showing subtype linked polymorphisms).

Minimization of genetic distances by the consensus, ancestral, and center-of-tree (COT) sequences for HIV-1 variants within an infected individual and the design of reagents to test immune reactivity

Eliciting maximal immune responses to highly divergent viruses is a challenge and a focus in AIDS vaccine development. Another challenge is to identify the immune correlates of protective immunity. Recent AIDS vaccine design approaches attempt to use reconstructed centralized viral sequences that minimize genetic differences to circulating viruses. Using these approaches, we derive and analyze consensus (CON), ancestral (ANC), and center-of-tree (COT) sequences to represent intra-individual HIV-1 env variants encoding a range of diversities and phylogenetic structures. Each reconstructed sequence significantly minimized genetic distances to extant sequences throughout the first 5 years of infection of an individual. Interestingly, ANC sequences diverged and were not significantly better than extant sequences in minimizing genetic distances at later stages of infection and disease, likely due to the development of a substantially asymmetric phylogeny. COT or CON sequences derived from autologous virus samplings may be useful for increasing the sensitivity of assessments of immune reactivity against HIV.

Correlation between SIV Tat evolution and AIDS progression in cerebrospinal fluid of morphine-dependent and control macaques infected with SIV and SHIV

Morphine abuse has been associated with higher virus replication and accelerated disease progression in a non-human primate model of AIDS. In our previous report, we have shown that 50% of morphine-addicted macaques progress rapidly and that 2/3 of the rapid progressors exhibit severe neuropathogenesis. In this report, we examined the sequence evolution of the SIV Tat protein, known to participate in AIDS neuropathology, in the cerebrospinal fluid (CSF) of morphine-dependent and control macaques over the first 20 weeks of infection. The CSF SIV Tat evolution was found to be inversely related with disease progression, and the highly neuropathogenic inoculum clone sequence was the prevalent CSF form in rapid progressors. Divergence from the inoculum clone was significantly greater in both morphine-dependent normal progressors and control macaques than in the morphine-dependent rapid progressors. Furthermore, we also found evidence of a trend that morphine alters the type of mutation, resulting in an enhanced ratio of transitions to transversions (Ts:Tv). Rapid disease exacerbates this trend and appears to influence the distribution of nonsynonymous changes in the first exon of SIV tat, with a clear majority of mutations occurring in the C-terminal half of the protein where the known functionally important domains reside. Thus, morphine abuse may change the nature and extent of mutations that drive viral evolution.

A novel simian immunodeficiency virus from black mangabey (Lophocebus aterrimus) in the Democratic Republic of Congo

In order to understand primate lentivirus evolution, characterization of additional simian immunodeficiency virus (SIV) strains is essential. Here, an SIV from a black mangabey (Lophocebus aterrimus) originating from the Democratic Republic of Congo was analysed phylogenetically. The monkey had cross-reactive antibodies against human immunodeficiency virus type 1 (HIV-1) and HIV-2. The viral pol region sequence was amplified by nested PCR and sequence analysis confirmed that it was related to known SIV sequences. This is the first report to characterize genetically an SIV from the monkey genus Lophocebus. Phylogenetic analysis of the pol region revealed that this novel SIV, designated SIVbkm, fell into the SIVsyk and SIVgsn virus group, containing viruses isolated from the genus Cercopithecus, and suggests that cross-species transmission has occurred between species of the genera Lophocebus and Cercopithecus.

Enhanced mucosal and systemic immune response with intranasal immunization of mice with HIV peptides entrapped in PLG microparticles in combination with Ulex Europaeus-I lectin as M cell target

The predominant route of HIV infection is through the sexual transmission via M cells. Most of the peptide and protein vaccines show poor transport across the epithelial barrier and are commonly administered by parenteral route. In the present study four HIV peptides from envelope (gp 41-LZ (leucine zipper), gp 41-FD (fusion domain) and gp120-C2) and regulatory (Nef) region in poly lactic-co-glycolide (PLG) micro-particle delivery were evaluated in mice of outbred and with different genetic background to compare immune response versus MHC restriction. Out of the combinational and single routes of immunization attempted, the single route maintained the IgG, IgA and sIgA in sera and washes for longer duration as compared to combinational routes in which the response was declined. The study demonstrated that single intranasal immunization offered significantly higher immune response (p<0.05) over oral and rectal mucosal routes in terms of inducing systemic as well as mucosal response. Also, the specific activity measurement of IgA and IgG in sera and sIgA in washes were correlating to the antibody titers. However, the intramuscular route of immunization generated systemic response only. The entrapment of plant lectin UEA-1 a ligand specific for M cells in micro-particle further enhanced the immune response in all the mucosal routes. The IgG isotypes generated were of IgG1 and IgG2a/2b in sera for all the peptides. The T cell proliferation response study with and without UEA-1 lectin in micro-particles showed significantly high (p<0.05) stimulation index (SI) with intranasal immunization for all the peptides from cells collected from spleen (SP), peyer's patches (PP) and lamina propria (LP) with SI in the order LP cells>PP>or=SP. The cytokine measurement profile of IL-2, IFN-gamma and IL-6 and low levels of IL-4 in the cultural supernatants of SP, PP and LP showed mixed CD4(+) Th1 and Th2 immune response. The p24 assay showed high percent inhibition of HIV-IIIB virus with sera and washes obtained from intranasal route. Thus, overall the study highlighted the combination of UEA-1 lectin with HIV peptides in micro-particles through intranasal immunization generated systemic as well as mucosal immune response.

Detecting selection in noncoding regions of nucleotide sequences

We present a maximum-likelihood method for examining the selection pressure and detecting positive selection in noncoding regions using multiple aligned DNA sequences. The rate of substitution in noncoding regions relative to the rate of synonymous substitution in coding regions is modeled by a parameter zeta. When a site in a noncoding region is evolving neutrally zeta = 1, while zeta > 1 indicates the action of positive selection, and zeta < 1 suggests negative selection. Using a combined model for the evolution of noncoding and coding regions, we develop two likelihood-ratio tests for the detection of selection in noncoding regions. Data analysis of both simulated and real viral data is presented. Using the new method we show that positive selection in viruses is acting primarily in protein-coding regions and is rare or absent in noncoding regions.

Molecular evolution of the hepatitis delta virus antigen gene: recombination or positive selection?

We present the statistical analysis of diversifying selective pressures on the hepatitis D antigen gene (HDAg). Thirty-three distinct HDAg sequences from subtypes I, II, and III were tested for positive selection using maximum likelihood methods based on models of codon substitution that allow variable selective pressures across sites. Such methods have been shown to be sufficiently accurate and successful in detecting positive selection in a variety of viral and nonviral protein-coding genes. About 11% of codon sites in HDAg were estimated to be under diversifying selection. Remarkably, most of the residues predicted to evolve under positive selection were located in the immunogenic domain and the N-terminus region with reported antigenic activity. These sites are potential targets of the host's immune response. Identification of residues mutating to escape immune recognition may help to distinguish the most virulent strains and aid vaccine design. Possible interplay between positive selection and recombination on the gene is discussed but no significant evidence for recombination was found.

Accuracy and power of statistical methods for detecting adaptive evolution in protein coding sequences and for identifying positively selected sites

The parsimony method of Suzuki and Gojobori (1999) and the maximum likelihood method developed from the work of Nielsen and Yang (1998) are two widely used methods for detecting positive selection in homologous protein coding sequences. Both methods consider an excess of nonsynonymous (replacement) substitutions as evidence for positive selection. Previously published simulation studies comparing the performance of the two methods show contradictory results. Here we conduct a more thorough simulation study to cover and extend the parameter space used in previous studies. We also reanalyzed an HLA data set that was previously proposed to cause problems when analyzed using the maximum likelihood method. Our new simulations and a reanalysis of the HLA data demonstrate that the maximum likelihood method has good power and accuracy in detecting positive selection over a wide range of parameter values. Previous studies reporting poor performance of the method appear to be due to numerical problems in the optimization algorithms and did not reflect the true performance of the method. The parsimony method has a very low rate of false positives but very little power for detecting positive selection or identifying positively selected sites.

Evidence for heterogeneous selective pressures in the evolution of the env gene in different human immunodeficiency virus type 1 subtypes

Recent studies have demonstrated the emergence of human immunodeficiency virus type 1 (HIV-1) subtypes with various levels of fitness. Using heterogeneous maximum-likelihood models of adaptive evolution implemented in the PAML software package, with env sequences representing each HIV-1 group M subtype, we examined the various intersubtype selective pressures operating across the env gene. We found heterogeneity of evolutionary mechanisms between the different subtypes with a category of amino acid sites observed that had undergone positive selection for subtypes C, F1, and G, while these sites had undergone purifying selection in all other subtypes. Also, amino acid sites within subtypes A and K that had undergone purifying selection were observed, while these sites had undergone positive selection in all other subtypes. The presence of such sites indicates heterogeneity of selective pressures within HIV-1 group M subtype evolution that may account for the various levels of fitness of the subtypes.

Influence of random genetic drift on human immunodeficiency virus type 1 env evolution during chronic infection

Human immunodeficiency virus type 1 (HIV-1) has high replication and mutation rates that generate large census populations and high levels of genetic variation. We examined the roles of natural selection, population growth, random genetic drift, and recombination in shaping the variation in 1509 C2-V5 env sequences derived from nine men with chronic HIV-1 infection. These sequences were obtained from clinical visits that reflect the first 6-13.7 years of infection. Pairwise comparisons of nonsynonymous and synonymous distances, Tajima's D test, Fu and Li's D* test, and a test of recurrent mutation revealed evidence for episodes of nonneutral evolution in a total of 22 out of 145 blood samples, representing six of the nine individuals. Using three coalescent-based maximum-likelihood estimators, we found viral effective population sizes in all nine individuals to be approximately 10(3). We also show that a previous estimate of the effective population size of approximately 10(5) based on rare haplotype frequencies decreases to approximately 10(3) upon correcting a biased sampling procedure. We conclude that the genetic variation in these data sets can be explained by a predominance of random genetic drift of neutral mutations with brief episodes of natural selection that were frequently masked by recombination.

Linkage of amino acid variation and evolution of human immunodeficiency virus type 1 gp120 envelope glycoprotein (subtype B) with usage of the second receptor

To clarify the relationship between the amino acid variations of the gp120 of human immunodeficiency virus type 1 (HIV-1) and the chemokine receptors that are used as the second receptor for HIV, we evaluated amino acid site variation of gp120 between the X4 strains (use CXCR4) and the R5 strains (use CCR5) from 21 sequences of subtype B. Our analysis showed that residues 306 and 322 in the V3 loop and residue 440 in the C4 region were associated with usage of the second receptor. The polymorphism at residue 440 is clearly associated with the usage of the second receptor: The amino acid at position 440 was a basic amino acid in the R5 strains, and a nonbasic and smaller amino acid in the X4 strains, while the V3 loop of the X4 strains was more basic than that of the R5 strains. This suggests that residue 440 in the C4 region, which is close to the V3 loop in the three-dimensional structure, is critical in determining which second receptor is used. Analysis of codon frequency suggests that, in almost all cases, the difference at residue 440 between basic amino acids in the R5 strains and nonbasic amino acids in the X4 strains could be due to a single nucleotide change. These findings predict that the evolutionary changes in amino acid residue 440 may be correlated with evolutionary changes in the V3 loop. One possibility is that a change in electric charge at residue 440 compensates for a change in electric charge in the V3 loop. The amino acid polymorphism at position 440 can be useful to predict the cell tropism of a strain of HIV-1 subtype B.

A statistical characterization of consistent patterns of human immunodeficiency virus evolution within infected patients

Within-patient HIV populations evolve rapidly because of a high mutation rate, short generation time, and strong positive selection pressures. Previous studies have identified "consistent patterns" of viral sequence evolution. Just before HIV infection progresses to AIDS, evolution seems to slow markedly, and the genetic diversity of the viral population drops. This evolutionary slowdown could be caused either by a reduction in the average viral replication rate or because selection pressures weaken with the collapse of the immune system. The former hypothesis (which we denote "cellular exhaustion") predicts a simultaneous reduction in both synonymous and nonsynonymous evolution, whereas the latter hypothesis (denoted "immune relaxation") predicts that only nonsynonymous evolution will slow. In this paper, we present a set of statistical procedures for distinguishing between these alternative hypotheses using DNA sequences sampled over the course of infection. The first component is a new method for estimating evolutionary rates that takes advantage of the temporal information in longitudinal DNA sequence samples. Second, we develop a set of probability models for the analysis of evolutionary rates in HIV populations in vivo. Application of these models to both synonymous and nonsynonymous evolution affords a comparison of the cellular-exhaustion and immune-relaxation hypotheses. We apply the procedures to longitudinal data sets in which sequences of the env gene were sampled over the entire course of infection. Our analyses (1) statistically confirm that an evolutionary slowdown occurs late in infection, (2) strongly support the immune-relaxation hypothesis, and (3) indicate that the cessation of nonsynonymous evolution is associated with disease progression.

[Sequence variation of HIV and bioinformatics]

The envelope glycoprotein of human immunodeficiency virus type 1 (HIV-1) interacts with receptors on the target cell and mediates virus entry by fusing the viral and cell membranes. To maintain the viral infectivity, amino acids that interact with receptors are expected to be more conserved than the other sites on the protein surface. In contrast to the functional constraint of amino acids for the receptor binding, some amino acid changes in this protein may produce antigenic variations that enable the virus to escape from recognition of the host immune system. Therefore, both positive selection (higher fitness) and negative selection (lower fitness) against amino acid changes are taking place during evolution of surface proteins of parasites To elucidate the evolutionary mechanisms of the whole HIV-1 gp120 envelope glycoprotein at the single site level, we collected and analyzed all available sequence data for the protein. By analyzing 186 sequences of the HIV-1 gp120 (subtype B), we reevaluated amino acid variability at the single site level, and estimated the numbers of synonymous and nonsynonymous substitutions at each codon position to detect positive and negative selection. We identified 33 amino acid positions which may be under positive selection. Some of these positions may form discontinuous epitopes. We also analyzed amino acid sequences to find amino acid positions responsible for usage of the second receptor. We found that, in addition to the V3 loop, amino acid variation at residue 440 in C4 region is clearly linked with the usage of CXCR 4.

Widespread adaptive evolution in the human immunodeficiency virus type 1 genome

We investigated variable selective pressures among amino acid sites in HIV-1 genes. Selective pressure at the amino acid level was measured by using the nonsynonymous/synonymous substitution rate ratio (omega = dN/dS). To identify amino acid sites under positive selection with omega > 1, we applied maximum likelihood models that allow variable omega ratios among sites to analyze genomic sequences of 26 HIV-1 lineages including subtypes A, B, and C. Likelihood ratio tests detected sites under positive selection in each of the major genes in the genome: env, gag, pol, vif, and vpr. Positive selection was also detected in nef, tat, and vpu, although those genes are very small. The majority of positive selection sites is located in gp160. Positive selection was not detected if omega was estimated as an average across all sites, indicating the lack of power of the averaging approach. Candidate positive selection sites were mapped onto the available protein tertiary structures and immunogenic epitopes. We measured the physiochemical properties of amino acids and found that those at positive selection sites were more diverse than those at variable sites. Furthermore, amino acid residues at exposed positive selection sites were more physiochemically diverse than at buried positive selection sites. Our results demonstrate genomewide diversifying selection acting on the HIV-1.

Comparative study of adaptive molecular evolution in different human immunodeficiency virus groups and subtypes

Molecular adaptation, as characterized by the detection of positive selection, was quantified in a number of genes from different human immunodeficiency virus type 1 (HIV-1) group M subtypes, group O, and an HIV-2 subtype using the codon-based maximum-likelihood method of Yang and coworkers (Z. H. Yang, R. Nielsen, N. Goldman, and A. M. K. Pedersen, Genetics 155:431-449, 2000). The env gene was investigated further since it exhibited the strongest signal for positive selection compared to those of the other two major HIV genes (gag and pol). In order to investigate the pattern of adaptive evolution across env, the location and strength of positive selection in different HIV-1 sequence alignments was compared. The number of sites having a significant probability of being positively selected varied among these different alignment data sets, ranging from 25 in HIV-1 group M subtype A to 40 in HIV-1 group O. Strikingly, there was a significant tendency for positively selected sites to be located at the same position in different HIV-1 alignments, ranging from 10 to 16 shared sites for the group M intersubtype comparisons and from 6 to 8 for the group O to M comparisons, suggesting that all HIV-1 variants are subject to similar selective forces. As the host immune response is believed to be the dominant driving force of adaptive evolution in HIV, this result would suggest that the same sites are contributing to viral persistence in diverse HIV infections. Thus, the positions of the positively selected sites were investigated in reference to the inferred locations of different epitope types (antibody, T helper, and cytotoxic T lymphocytes) and the positions of N and O glycosylation sites. We found a significant tendency for positively selected sites to fall outside T-helper epitopes and for positively selected sites to be strongly associated with N glycosylation sites.

Unifying the epidemiological and evolutionary dynamics of pathogens

A key priority for infectious disease research is to clarify how pathogen genetic variation, modulated by host immunity, transmission bottlenecks, and epidemic dynamics, determines the wide variety of pathogen phylogenies observed at scales that range from individual host to population. We call the melding of immunodynamics, epidemiology, and evolutionary biology required to achieve this synthesis pathogen "phylodynamics." We introduce a phylodynamic framework for the dissection of dynamic forces that determine the diversity of epidemiological and phylogenetic patterns observed in RNA viruses of vertebrates. A central pillar of this model is the Evolutionary Infectivity Profile, which captures the relationship between immune selection and pathogen transmission.

Genetic and functional analysis of full-length human immunodeficiency virus type 1 env genes derived from brain and blood of patients with AIDS

The genetic evolution of human immunodeficiency virus type 1 (HIV-1) in the brain is distinct from that in lymphoid tissues, indicating tissue-specific compartmentalization of the virus. Few primary HIV-1 envelope glycoproteins (Envs) from uncultured brain tissues have been biologically well characterized. In this study, we analyzed 37 full-length env genes from uncultured brain biopsy and blood samples from four patients with AIDS. Phylogenetic analysis of intrapatient sequence sets showed distinct clustering of brain relative to blood env sequences. However, no brain-specific signature sequence was identified. Furthermore, there was no significant difference in the number or positions of N-linked glycosylation sites between brain and blood env sequences. The patterns of coreceptor usage were heterogeneous, with no clear distinction between brain and blood env clones. Nine Envs used CCR5 as a coreceptor, one used CXCR4, and two used both CCR5 and CXCR4 in cell-to-cell fusion assays. Eight Envs could also use CCR3, CCR8, GPR15, STRL33, Apj, and/or GPR1, but these coreceptors did not play a major role in virus entry into microglia. Recognition of epitopes by the 2F5, T30, AG10H9, F105, 17b, and C11 monoclonal antibodies varied among env clones, reflecting genetic and conformational heterogeneity. Envs from two patients contained 28 to 32 N-glycosylation sites in gp120, compared to around 25 in lab strains and well-characterized primary isolates. These results suggest that HIV-1 Envs in brain cannot be distinguished from those in blood on the basis of coreceptor usage or the number or positions of N-glycosylation sites, indicating that other properties underlie neurotropism. The study also demonstrates characteristics of primary HIV-1 Envs from uncultured tissues and implies that Env variants that are glycosylated more extensively than lab strains and well-characterized primary isolates should be considered during development of vaccines and neutralizing antibodies.

Parasitic exploitation as an engine of diversity

Parasitic exploitation occurs within and between a wide variety of taxa in a plethora of diverse contexts. Theoretical and empirical analyses indicate that parasitic exploitation can generate substantial genetic and phenotypic polymorphism within species. Under some circumstances, parasitic exploitation may also be an important factor causing reproductive isolation and promoting speciation. Here we review research relevant to the relationship between parasitic exploitation, within species-polymorphism, and speciation in some of the major arenas in which such exploitation has been studied. This includes research on the vertebrate major histocompatibility loci, plant-pathogen interactions, the evolution of sexual reproduction, intragenomic conflict, sexual conflict, kin mimicry and social parasitism, tropical forest diversity and the evolution of language. We conclude by discussing some of the issues raised by comparing the effect of parasitic exploitation on polymorphism and speciation in different contexts.

Effect of recombination on the accuracy of the likelihood method for detecting positive selection at amino acid sites

Maximum-likelihood methods based on models of codon substitution accounting for heterogeneous selective pressures across sites have proved to be powerful in detecting positive selection in protein-coding DNA sequences. Those methods are phylogeny based and do not account for the effects of recombination. When recombination occurs, such as in population data, no unique tree topology can describe the evolutionary history of the whole sequence. This violation of assumptions raises serious concerns about the likelihood method for detecting positive selection. Here we use computer simulation to evaluate the reliability of the likelihood-ratio test (LRT) for positive selection in the presence of recombination. We examine three tests based on different models of variable selective pressures among sites. Sequences are simulated using a coalescent model with recombination and analyzed using codon-based likelihood models ignoring recombination. We find that the LRT is robust to low levels of recombination (with fewer than three recombination events in the history of a sample of 10 sequences). However, at higher levels of recombination, the type I error rate can be as high as 90%, especially when the null model in the LRT is unrealistic, and the test often mistakes recombination as evidence for positive selection. The test that compares the more realistic models M7 (beta) against M8 (beta and omega) is more robust to recombination, where the null model M7 allows the positive selection pressure to vary between 0 and 1 (and so does not account for positive selection), and the alternative model M8 allows an additional discrete class with omega = d(N)/d(S) that could be estimated to be >1 (and thus accounts for positive selection). Identification of sites under positive selection by the empirical Bayes method appears to be less affected than the LRT by recombination.

Turnover of env variable region 1 and 2 genotypes in subjects with late-stage human immunodeficiency virus type 1 infection

The env gene of human immunodeficiency virus type 1 (HIV-1) includes some of the most genetically diverse regions of the viral genome, which are called variable regions 1 through 5 (V1 through V5). We have developed a heteroduplex tracking assay to detect changes in variable regions 1 and 2 of env (V1/V2-HTA). Using sequences from two molecular clones as probes, we have studied the nature of longitudinal virus population changes in a cohort of HIV-1-infected subjects. Viral sequences present in 21 subjects with late-stage HIV-1 infection were initially screened for stability of the virus population by V1/V2-HTA. The virus populations at entry comprised an average of five coexisting V1/V2 genotypic variants (as identified by HTA). Eight of the 21 subjects were examined in detail because of the dynamic behavior of their env variants over an approximately 9-month period. In each of these cases we detected a single discrete transition of V1/V2 genotypes based on monthly sampling. The major V1/V2 genotypes (those present at >10% abundance) from the eight subjects were cloned and sequenced to define the nature of V1/V2 variability associated with a discrete transition. Based on a comparison of V1/V2 genotypic variants present at entry with the newly emerged variants we categorized the newly emerged variants into two groups: variants without length differences and variants with length differences. Variants without length differences had fewer nucleotide substitutions, with the changes biased to either V1 or V2, suggestive of recent evolutionary events. Variants with length differences included ones with larger numbers of changes that were distributed, suggestive of recall of older genotypes. Most length differences were located in domains where the codon motif AVT (V = A, G, C) had become enriched and fixed. Finally, recombination events were detected in two subjects, one of which resulted in the reassortment of V1 and V2 regions. We suggest that turnover in V1/V2 populations was largely driven by selection on either V1 or V2 and that escape was accomplished either through changes focused in the region under selection or by the appearance of a highly divergent variant.

Kaposi's sarcoma-associated herpesvirus cytotoxic T lymphocytes recognize and target Darwinian positively selected autologous K1 epitopes

Kaposi's sarcoma-associated herpesvirus (KSHV) is the infectious cause of Kaposi's sarcoma (KS) and certain lymphoproliferations particularly in the context of human immunodeficiency virus (HIV) type 1-induced immunosuppression. The introduction of effective therapies to treat HIV has led to a decline in the incidence of KS, suggesting that immune responses may play a role in controlling KSHV infection and pathogenesis. Cytotoxic-T-lymphocyte (CTL) activity against KSHV proteins has been demonstrated; however, the identification of KSHV CTL epitopes remains elusive and problematic. Although the herpesvirus genomic layout is generally conserved, KSHV encodes a unique hypervariable protein, K1, with intense biological selection pressure at specific amino acid sites. To investigate whether this variability is partly driven by cellular immunity, we designed K1 peptides that match only the unique viral sequence for every individual studied here (autologous peptides). We identified functional CTL epitopes within K1's most variable areas, and we show that a given individual responds only to autologous peptides and not to peptides from other individuals. Furthermore, these epitopes are highly conserved sequences within KSHV isolates from a specific strain but are not conserved between different strains. We conclude that CTL recognition contributes to K1, and therefore to KSHV, evolution.

Network analysis of human and simian immunodeficiency virus sequence sets reveals massive recombination resulting in shorter pathways

The intrinsic recombination rate of human immunodeficiency virus (HIV) exceeds the point mutation rate by a factor of 10. As the majority of infected cells in vivo harbour multiple proviruses, the stage is set for rampant recombination. Therefore, it may be presumed that phylogenic relationships and mutation frequencies will probably be affected by recombination. However, the proportion of homoplasies arising from recombination and mutation is not known. By studying the evolution of the hypervariable regions of the simian immunodeficiency virus envelope gene among four macaques, it is shown that homoplasies arise more from recombination than from point mutation. When recombination is accounted for, the minimum number of substitutions in a sequence set may be reduced by as much as 45 %. In fact, the true number of point mutations in a set of HIV sequences tends to the number of discrete substitutions. Hence, lineages are younger than anticipated previously, although not in proportion to the ratio of the intrinsic recombination/point mutation rate. Recombination also inflates codon polymorphisms.

Analysis of the functional relationship between V3 loop and gp120 context with regard to human immunodeficiency virus coreceptor usage using naturally selected sequences and different viral backbones

The human immunodeficiency virus type 1 (HIV-1) gp120 V3 loop plays a predominant role in chemokine receptor usage; however, other linear and nonlinear gp120 domains are involved in this step of the HIV-1 replication cycle. At present, the functional relationship between V3 and these domains with regard to coreceptor usage is unclear. To gain insights into the nature of this relationship in naturally selected viral variants, we developed a recombinant strategy based on two different gp120 backbones derived from CXCR4 (X4)- and CCR5 (R5)-tropic viral strains, respectively. Using this recombinant model system, we evaluated the phenotype patterns conferred to chimeric viruses by exogenous V3 loops from reference molecular clones and samples from infected subjects. In 13 of 17 recombinants (76%), a comparable phenotype was observed independently of the gp120 backbone, whereas in a minority of the recombinant viruses (4/17, 24%) viral infectivity depended on the gp120 context. No case of differential tropism using identical V3 sequence in the two gp120 contexts was observed. Site-directed mutagenesis experiments were performed to evaluate the phenotypic impact of specific V3 motifs. The data indicate that while the interaction of HIV-1 with chemokine receptors is driven by V3 loop and influenced by its evolutionary potential, the gp120 context plays a role in influencing the replication competence of the variants, suggesting that compensatory mutations occurring at sites other than V3 are necessary in some cases.

Human cytomegalovirus glycoprotein N (gpUL73-gN) genomic variants: identification of a novel subgroup, geographical distribution and evidence of positive selective pressure

Human cytomegalvirus (HCMV) ORF UL73 is a polymorphic locus, encoding the viral glycoprotein gpUL73-gN, a component of the gC-II envelope complex. The previously identified gN genomic variants, denoted gN-1, gN-2, gN-3 and gN-4, were further investigated in this work by analysing a large panel of HCMV clinical isolates collected from all over the world (223 samples). Sequencing and phylogenetic analysis confirmed the existence of the four gN genotypes, but also allowed the identification of a novel subgroup belonging to the gN-3 genotype, which was designated gN-3b. The number of non-synonymous (d(N)) and synonymous (d(S)) nucleotide substitutions and their ratio (d(N)/d(S)) were estimated among the gN genotypes to evaluate the possibility of positive selection. Results showed that the four variants evolved by neutral (random) selection, but that the gN-3 and gN-4 genotypes are maintained by positive selective pressure. The 223 HCMV clinical isolates were subdivided according to their geographical origin, and four main regions of gN prevalence were identified: Europe, China, Australia and Northern America. The gN variants were found to be widespread and represented within the regions analysed without any significant difference, and no new genotype was detected. Finally, for clinical and epidemiological purposes, a rapid and low-cost method for genetic grouping of the HCMV clinical isolates was developed based on the RFLP revealed by SacI, ScaI and SalI digestion of the PCR-amplified UL73 sequence. This technique enabled us to distinguish all four gN genomic variants and also their subtypes.

Molecular evolution meets the genomics revolution

Changes in technology in the past decade have had such an impact on the way that molecular evolution research is done that it is difficult now to imagine working in a world without genomics or the Internet. In 1992, GenBank was less than a hundredth of its current size and was updated every three months on a huge spool of tape. Homology searches took 30 minutes and rarely found a hit. Now it is difficult to find sequences with only a few homologs to use as examples for teaching bioinformatics. For molecular evolution researchers, the genomics revolution has showered us with raw data and the information revolution has given us the wherewithal to analyze it. In broad terms, the most significant outcome from these changes has been our newfound ability to examine the evolution of genomes as a whole, enabling us to infer genome-wide evolutionary patterns and to identify subsets of genes whose evolution has been in some way atypical.

Immune-mediated positive selection drives human immunodeficiency virus type 1 molecular variation and predicts disease duration

Using likelihood-based evolutionary methods, we demonstrate that the broad genetic diversity of human immunodeficiency virus type 1 (HIV-1) in an infected individual is a consequence of site-specific positive selection for diversity, a likely consequence of immune recognition. In particular, the extent of positive selection appears to be a good predictor of disease duration. Positively selected sites along HIV-1 partial env sequences are numerous but not distributed uniformly. In a sample of eight patients studied longitudinally, the proportion of sites per sample under positive selection was a statistically significant predictor of disease duration. Among long-term progressors, positive selection persisted at sites over time and appears to be associated with helper T-cell epitopes. In contrast, sites under positive selection shifted from one longitudinal sample to the next in short-term progressors. Our study is consistent with the hypothesis that a broad and persistent immunologic response is associated with a slower rate of disease progression. In contrast, narrow, shifting immune responses characterize short-term progressors.

Selection for human immunodeficiency virus type 1 recombinants in a patient with rapid progression to AIDS

Although human immunodeficiency virus type 1 (HIV-1) recombinants have been found with high frequency, little is known about the forces that select for these viruses or their importance to pathogenesis. Here we document the emergence and dynamics of 11 distinct HIV-1 recombinants in a man who was infected with two subtype B HIV-1 strains and progressed rapidly to AIDS without developing substantial cellular or humoral immune responses. Although numerous frequency oscillations were observed, a single recombinant lineage eventually came to dominate the population. Numerical simulations indicate that the successive recombinant forms displaced each other too rapidly to be explained by any simple model of random genetic drift or sampling variation. All of the recombinants, including several resulting from independent recombination events, possessed the same sequence motif in the V3 loop, suggesting intense selection on this segment of the viral envelope protein. The outgrowth of the predominant V3 loop recombinants was not, however, associated with changes in coreceptor utilization. The final variant was instead notable for having lost 3 of 14 potential glycosylation sites. We also observed high ratios of synonymous-to-nonsynonymous nucleotide changes-suggestive of purifying selection-in all viral populations, with particularly high ratios in newly arising recombinants. Our study, therefore, illustrates the unusual and important patterns of viral adaptation that can occur in a patient with weak immune responses. Although it is hard to tease apart cause and effect in a single patient, the correlation with disease progression in this patient suggests that recombination between divergent viruses, with its ability to create chimeras with increased fitness, can accelerate progression to AIDS.

Simulation study of the reliability and robustness of the statistical methods for detecting positive selection at single amino acid sites

Inferring positive selection at single amino acid sites is of biological and medical importance. Parsimony-based and likelihood-based methods have been developed for this purpose, but the reliabilities of these methods are not well understood. Because the evolutionary models assumed in these methods are only rough approximations to reality, it is desirable that the methods are not very sensitive to violation of the assumptions made. In this study we show by computer simulation that the likelihood-based method is sensitive to violation of the assumptions and produces many false-positive results under certain conditions, whereas the parsimony-based method tends to be conservative. These observations, together with those from previous studies, suggest that the positively selected sites inferred by the parsimony-based method are more reliable than those inferred by the likelihood-based method.