Advertisement

  • Loading metrics

Open Access

Peer-reviewed

Research Article

Adaptation to Human Populations Is Revealed by Within-Host Polymorphisms in HIV-1 and Hepatitis C Virus

  • Art F. Y Poon ,

    To whom correspondence should be addressed. E-mail: afpoon@ucsd.edu

    Affiliation Department of Pathology, University of California San Diego, La Jolla, California, United States of America

    ⨯
  • Sergei L. Kosakovsky Pond,

    Affiliation Department of Pathology, University of California San Diego, La Jolla, California, United States of America

    ⨯
  • Phil Bennett,

    Affiliation Science Park, University of Warwick, Coventry, United Kingdom

    ⨯
  • Douglas D Richman,

    Affiliations Department of Pathology, University of California San Diego, La Jolla, California, United States of America , School of Medicine, University of California San Diego, La Jolla, California, United States of America , Veterans Affairs San Diego Healthcare System, San Diego, California, United States of America

    ⨯
  • Andrew J. Leigh Brown,

    Affiliation Institute of Evolutionary Biology, University of Edinburgh, Edinburgh, Scotland, United Kingdom

    ⨯
  • Simon D. W Frost

    Affiliation Department of Pathology, University of California San Diego, La Jolla, California, United States of America

    ⨯

Abstract

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 ≥ 3,000) and HCV-infected patients (n ≥ 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.

Author Summary

The rapid accumulation of genetic variation in human viruses, such as human immunodeficiency virus type 1 (HIV-1) and hepatitis C virus (HCV), enables these pathogens to elude the immune system and forestalls the development of effective vaccines. This variation may be shaped by selection due to host-specific immune responses, such that the total virus population mirrors the immune diversity of the host population. However, the often-neglected viral genetic variation within hosts may also play an important role in shaping variation in the total virus population. We carry out an innovative analysis of bulk HIV-1 and HCV sequences isolated from over 4,000 human patients, exploiting “mixtures” (i.e., ambiguous nucleotides) as a site-specific marker of within-host genetic variation. We find that nonsynonymous mixtures are disproportionately abundant at codon positions under strong host-specific immune selection. Because existing models of virus evolution provide no explanation for this outcome, we have formulated a new model supplemented with stochastic simulations to demonstrate that the rapid transmission of viruses through diverse selective environments creates a positive correlation between nonsynonymous variation within and among hosts.

Citation: Poon AFY, Kosakovsky Pond SL, Bennett P, Richman DD, Leigh Brown AJ, Frost SDW (2007) Adaptation to Human Populations Is Revealed by Within-Host Polymorphisms in HIV-1 and Hepatitis C Virus. PLoS Pathog 3(3): e45. https://doi.org/10.1371/journal.ppat.0030045

Editor: Richard A. Koup, National Institute of Allergy and Infectious Diseases, United States of America

Received: October 11, 2006; Accepted: February 11, 2007; Published: March 30, 2007

This is an open-access article distributed under the terms of the Creative Commons Public Domain declaration which stipulates that, once placed in the public domain, this work may be freely reproduced, distributed, transmitted, modified, built upon, or otherwise used by anyone for any lawful purpose.

Funding: This work was supported by the US National Institutes of Health (grant numbers AI43638, AI47745, AI57167), the Universitywide AIDS Research Program (grant number IS02-SD-701), and UCSD Centers for AIDS Research (CFAR)/National Institute of Allergy and Infectious Diseases (NIAID) Developmental Awards to SDWF and SLKP (grant number AI36214).

Competing interests: The authors have declared that no competing interests exist.

Abbreviations: CTL, cytotoxic T-lymphocyte; HCV, hepatitis C virus; HIV-1, human immunodeficiency virus type 1; HLA, human leukocyte antigen; PR, protease; RT, reverse transcriptase

Introduction

The cellular immune response mediated by CD8+ cytotoxic T-lymphocytes (CTLs) performs a critical role in the immune control of human viruses such as human immunodeficiency virus (HIV-1) [1] and hepatitis C virus (HCV) [2]. Consequently, the major histocompatibility (MHC) loci that encode the human leukocyte antigen (HLA) class I molecules, which recognize and bind CTL epitopes in viral proteins, are among the most highly polymorphic genes in the human population [3]. Nevertheless, the CTL response often fails to control the infection completely because of mutations that occur within HLA-restricted CTL epitopes, enabling the virus to escape binding and recognition [4]. Because epitopes are often located in functionally conserved regions of the viral genome, escape mutations may become costly to maintain in the absence of a selective HLA allele [5,6]. Thus, when an escape variant is transmitted between HLA-mismatched individuals, competitive growth frequently selects for reversion of the mutation to wild-type, as demonstrated experimentally in simian immunodeficiency virus–infected rhesus macaques [7] and in a comparative study of HIV-1–infected human patients [8].

Consequently, host-specific selection for escape or reversion may play an important role in shaping genetic variation in the circulating virus population [1,2,5,9,10]. For instance, population-based analyses of HIV-1 [9] and HCV [11] sequences have found several significant associations between divergent sites within CTL epitopes and the selective HLA alleles in the host population, suggesting that the frequency of escape polymorphisms in the circulating virus population are directly shaped by the immune diversity of the host population. Furthermore, the viral load of HIV-1–infected individuals has been found to be positively correlated with the frequency of the corresponding HLA supertypes in the host population, implying that the total virus population is adapting to the most frequent HLA supertypes [12]. If escape variants are readily transmitted between hosts, then a host with a common HLA supertype is more likely to encounter a virus that has already escaped its immune response [13], conferring a selective advantage to rare HLA supertypes. However, the virus genotype that becomes transmitted to the next host does not necessarily represent the ultimate outcome of adaptation to the previous host. Escape variants that have been transmitted into a host lacking a selective HLA allele can persist over long periods of time before reversion, or fail to revert at all over the duration of the study [8,14]. A delay or absence of reversion may be attributable to weak selection, when the fitness of the escape variant is either intrinsically high, or it has acquired compensatory mutations.

To evaluate the role of CTL-mediated selection in shaping the genetic variation of human viruses, we have carried out a large-scale analysis of HIV-1 and HCV protein-coding sequences isolated from human hosts. Previous analyses of clonal HIV-1 subtype B envelope [5,15] and protease (PR) [16] sequences have shown that across codon positions, genetic variation within hosts is positively correlated with variation among hosts. These correlations suggest that the genetic variation at both levels of the virus population is being shaped by a common set of biological constraints. However, the use of clonal sequences to characterize within-host variation restricted these analyses to small samples of hosts (n ≤ 12). In addition, quantifying the influence of selection on genetic variation within and among hosts is potentially confounded by variation in mutation rates among codon positions. Because mutation is the ultimate source of all genetic variation, site-specific variation at either level will be roughly proportional to the local mutation rate, which can yield a positive correlation in the absence of selection [17]. Indeed, this effect constitutes the basis for several tests of non-neutral evolution in genetic sequences [18–20].

To address the problem of limited sample size, we exploit “sequencing mixtures” as a site-specific marker of genetic variation within hosts. A sequencing mixture occurs when multiple distinct peaks occur above the same position in a sequencing electropherogram [21]; by convention, mixtures are encoded in sequences by ambiguous nucleotide characters (International Union of Pure and Applied Chemistry symbols “M”, “R”, “W”, “S”, “Y”, and “K”). Because mixtures can indicate the presence of a nucleotide polymorphism in the population, population-based (or “bulk”) sequencing is employed to detect minority variants that occur at frequencies above 10%–25% [21–23]. Although population-based sequencing may fail to detect mixtures below this threshold, transient polymorphisms under selection are more likely to be sampled at intermediate frequencies. This application of mixtures is particularly relevant to viruses with extremely high mutation rates such as HIV-1 and HCV, for which population-based sequences are exceedingly abundant. In this study, we use mixtures to quantify the effect of selection on within-host variation in population-based sequences of RT-PCR–amplified viral RNA from blood plasma isolated from over 4,000 HIV-1– or HCV-infected patients.

To remove the confounding effect of variation in mutation rates, we normalized the nonsynonymous variation per codon position by the synonymous variation, for either level of the virus population. Thus, we calculated the site-specific difference between the frequencies of nonsynonymous (mN) and synonymous mixtures (mS), and estimated the analogous difference between the rates of nonsynonymous (dN) and synonymous substitution (dS). Our estimates of mN and dN were both scaled by the expected number of nonsynonymous sites at each codon position; likewise, estimates mS and dS were scaled by the expected number of synonymous sites in the codon. The difference in substitution rates (dN − dS) is a conventional summary statistic for diversifying selection among hosts, i.e., host-specific selection causing nonsynonymous variation to accumulate among individual virus populations. We propose that the difference in mixture frequencies (mN − mS) can be employed as a summary statistic characterizing selection within each host. For instance, mN − mS > 0 can represent transient nonsynonymous polymorphisms undergoing directional selection (which drives the fixation of a specific variant within the host). Using these quantities, we will show that the distribution of mixtures in our samples of HIV-1 and HCV sequences cannot be explained by variation in mutation rates alone, and that host-specific selection is an important force shaping variation at both levels of the total virus population.

Because existing models of virus evolution seldom account for genetic variation both within and among hosts (but see [24,25]), we formulate a novel yet simple model that invokes both host-specific selection and rapid transmission between hosts to explain the observed patterns of genetic variation within and among hosts infected by HIV-1 or HCV. Bolstered by stochastic simulations, our model specifies the conditions that yield this outcome, and quantitatively predicts the joint effect of selection and transmission on the genetic composition of the circulating virus population. We find that when host-specific selection for escape and reversion is unbalanced and the transmission rate is high, then the frequency of escape variants becomes considerably skewed from expectations derived from the immune diversity of the host population. Failing to account for this effect may lead to erroneous conclusions on the overall importance of CTL-mediated selection in directing the evolution of the total virus population, or the relative contribution of specific CTL epitopes. Furthermore, the design of an effective vaccine to human viruses such as HIV-1 or HCV is highly contingent upon our ability to anticipate the response of an infection to CTL-mediated selection.

Results

Sequencing Mixtures Reveal CTL Selection

We screened for sequencing mixtures in population-based sequences of HIV-1 PR (n = 3,458) and reverse transcriptase (RT, n = 1,997) isolated from 3,004 and 1,989 treatment-naive individuals, respectively, and HCV sequences of envelope protein E1 (n = 2,691) and the hyper-variable region HVR1 of envelope protein E2 (n = 346). Although many sequences had at least one mixture (55% HIV-1, 63% HCV), there were relatively few mixtures per sequence on average (0.015 mixtures per codon position in HIV-1, 0.011 in HCV), suggesting that only a small number of codon positions had mixtures at detectable (20%–80%) frequencies in a given host (Figure S1). We found substantial variation among codon positions in mixture frequencies (Figure S2), which was greater for nonsynonymous (coefficient of variation = 1.98 HIV-1, 1.28 HCV) than synonymous mixtures (0.95 HIV-1, 1.06 HCV). There was no significant correlation between nonsynonymous and synonymous mixture frequencies per codon position in either HIV-1 (RT, Pearson's ρ = 0.04, p-value = 0.52; PR, ρ = 0.13, p-value = 0.21) or HCV gene sequences (E1, ρ = 0.01, p-value = 0.75; E2, ρ = −0.13, p-value = 0.18), indicating that the variation in mixture frequencies among codon positions was not simply due to local mutation rates.

The difference between nonsynonymous and synonymous mixture frequencies (mN − mS) was highly correlated with the difference between nonsynonymous and synonymous substitution rates (dN − dS) per codon position for both HIV-1 and HCV gene sequences (Figure 1A). This positive correlation between dN − dS and mN − mS remained significant for both E1 and E2 gene sequences even when different genotypes of HCV were analyzed separately. Overall, the quantity dN − dS assumed a negative value when averaged across the gene sequence, implying that nonsynonymous variation at the majority of codon positions was largely neutral or deleterious throughout the host population. Nevertheless, we detected significant diversifying selection (dN − dS > 0) at nine codon positions in HIV-1 PR (12, 13, 19, 35, 37, 63, 64, 77, and 93) and eight positions in RT (35, 39, 102, 122, 135, 200, 211, and 245) after correcting for the false-discovery rate [26] (α = 0.05); likewise, significant diversifying selection was attributed to several codon positions in HCV E1 and E2 (HVR1) sequences, which varied by genotype.

thumbnail
Download:
Figure 1. Genetic Variation within Hosts Is Shaped by Host-Specific Selection for CTL Escape

(A) The difference in nonsynonymous and synonymous mixture frequencies within hosts (mN − mS) is positively correlated with diversifying selection among hosts (dN − dS) per codon position. Each point corresponds to a unique codon position in the respective gene sequence. Dashed lines indicate the mean value for each quantity, which is consistently negative in dN − dS, implying purifying selection overall. Solid lines indicate a linear fit to the data. HCV genotypes are plotted separately as shown in the figure legends. A single outlier caused by a rare substitution lies outside the plot region for HIV-1 RT, but does not influence the significance of this correlation (Pearson's ρ = 0.619, p-value < 3 × 10−16).

(B) Selection for CTL escape elevates the frequency of nonsynonymous mixtures (solid circles) relative to synonymous mixtures (open triangles) at anchor residues within known A2-supertype–restricted epitopes in HIV-1 PR and RT and HCV E1 (predicted). Asterisks indicate anchor residues associated with disproportionately high frequencies of nonsynonymous mixtures.

https://doi.org/10.1371/journal.ppat.0030045.g001

For specific CTL epitopes in HIV-1 PR, RT, and HCV E1 sequences, we observed disproportionately higher frequencies of nonsynonymous mixtures at the anchor residues (Figure 1B) critical for MHC binding. In contrast, the profile of synonymous mixture frequencies within these epitopes lacked any distinct peaks in association with anchor residues. Overall, the median difference between the frequencies of nonsynonymous and synonymous mixtures was significantly greater at known HLA-B–restricted epitopes (median mN − mS = −0.2% mixtures per sequence per site) than in the remainder of the HIV-1 RT sequence (−0.5%; Wilcoxon rank-sum test, p-value = 0.007). We also found that mN − mS was greater at the anchor residues of HLA-B–restricted epitopes (median = −0.2%) than in an equivalent random sample of codon positions from HIV-1 RT on average (median = −0.4%), but this difference was only marginally significant (p-value = 0.11). In contrast, the median was not significantly greater at the known HLA-A–restricted epitopes within RT (Wilcoxon rank-sum test, p-value = 0.22), consistent with previous studies suggesting that HLA-B alleles assume a dominant role in the CTL control of HIV-1 [9,27]. In HIV-1 PR, the median excess in nonsynonymous mixtures was considerably greater within the single known HLA-B–restricted epitope (median = 0.7%) than in the rest of the gene sequence (median = −0.4%), but this difference was only marginally significant due to the small sample of codon positions (Wilcoxon rank-sum test, p-value = 0.1). Again, there was no significant difference in median values between HLA-A–restricted epitopes and the remainder of the PR sequence (Wilcoxon rank-sum test, p-value = 0.55).

Similarly, in the HCV E1 sequences, we found that the median excess of nonsynonymous mixtures was significantly greater within the two known HLA-B–restricted epitopes (median = 0.9%) than in an equivalent random sample of codon positions (median = −0.2%; Wilcoxon rank-sum test, p-value = 0.023). However, the median value for known HLA-A–restricted epitopes in HCV E1 was significantly less (median = −0.5%) than that in the remaining codon positions (median = −0.1%; Wilcoxon rank-sum test, p-value = 0.003). There were only two known CTL epitopes in the HCV E2 HVR1 sequence, both classified as HLA-A–restricted. We found no significant association between the quantity mN − mS and codon positions located within these epitopes (Wilcoxon rank-sum test, p-value = 0.87). In sum, nonsynonymous mixtures tend to accumulate disproportionately at codon positions under CTL selection, preferentially within HLA-B–restricted epitopes.

Simulation Results

A surplus of nonsynonymous mixtures within CTL epitopes represents transient polymorphisms that are eventually driven to fixation in the host by selection for escape or reversion [28]. This implies that the probability of sampling nonsynonymous sequencing mixture should decline with the intensity of host-specific selection at that codon position. As a result, host-specific selection would produce negative correlation between mN − mS and dN − dS across codon positions in the range dN − dS > 0, contrary to what we have observed in HIV-1 and HCV gene sequences. This paradox can be reconciled by incorporating the early transmission of unfavorable variants into a model of virus evolution (Figure 2). When selection and transmission act on similar time scales, the composition of the circulating virus population (i.e., the source of new infections) will not necessarily match the diversity of HLA genotypes in the host population. Suppose that an escape variant is transmitted from a host with a rare HLA genotype to a new host with a common HLA genotype. If the escape variant cannot outcompete the wild-type virus in the absence of a CTL response, then selection will favor reversion [7,8]. But the selective advantage of the wild-type virus may be so narrow that a substantial probability remains of transmitting the original escape variant [8,14]. Under such conditions, the severe bottleneck upon transmission could fix either the wild-type or escape variant in the new individual population (Figure 2). Because the next host will likely have the common HLA genotype, this transmission event can recreate the selective conditions requiring a transient nonsynonymous polymorphism to occur.

thumbnail
Download:
Figure 2. Effect of Transmission Rate on the Frequency of Mixtures

This schematic depicts the transmission chain of a virus population, where each host is represented by an enclosed graph that represents the evolving frequency of a CTL escape variant over time. The hosts either possess an HLA allele which favors the escape variant (HLA+, orange-shaded boxes) or the wild-type virus (HLA−, white-shaded boxes). A severe transmission bottleneck causes the population in the next host to be initially fixed for either the wild-type or escape variant (filled circle). If selection for escape or reversion is sufficiently strong (upper schematic in blue), then the favored virus genotype will tend to become fixed within the host before transmission occurs (open circle). Under such conditions, transient polymorphisms will only occur whenever the virus is transmitted between hosts of opposite type. On the other hand, if transmission and selection occur on similar time scales (lower schematic in red), then the host type does not necessarily predict which virus genotype becomes transmitted, causing transient polymorphisms to become more abundant (starred boxes).

https://doi.org/10.1371/journal.ppat.0030045.g002

To investigate this hypothesis, we implemented a simulation of allele frequency evolution within individual virus populations with ongoing transmission through a succession of hosts. Each individual virus population was represented by a single locus containing either an escape variant (at frequency p) or the wild-type allele. We assumed that transmission of the virus to a new host involved a severe bottleneck, such that the next population was initially fixed for either the escape variant (with probability p) or wild-type allele. Viral fitness in a given host was determined by a single MHC locus, at which an allele restricting the wild-type virus (HLA+) was present at a frequency q in the host population. We observed that the mean frequency of within-host polymorphisms fpoly:0.2 ≤ p ≤ 0.8 converged over time to an equilibrium value, which declined with stronger host-specific selection if the transmission rate was low (Figure 3A). On the other hand, if the transmission rate was high, then fpoly increased with stronger selection and thereby became positively correlated with genetic variation among hosts.

thumbnail
Download:
Figure 3. Factors Influencing Within-Host Polymorphisms and the Global Frequency of Escape Variants

(A) A contour plot depicting the mean effect of selection and transmission rate on displacing the frequency of detectable polymorphisms from the neutral expectation (Δfpoly, refer to the color key), as estimated from simulations. (The expectation E[fpoly] is jointly determined by the forward and back mutation rates, μ and ν, and population size, N.) The x-axis corresponds to the log-transformed transmission rate, log10k. The y-axis represents the mean log-transformed selection coefficient, E(log10s) = qlog10(sesc) + (1 − q)log10(srev).

(B) A 10-fold disparity in selection intensities sesc = 0.02, srev = 0.002) causes π˄ to substantially exceed q with increasing transmission rate, k. Each set of points represents mean estimates of π˄ from simulations (with virus population size N = 5,000 and μ = ν = 10−4). Dashed lines indicate predicted values from the deterministic model, which performs poorly when k is too high (i.e., when transmissions occur rapidly, allele frequencies are almost always near zero or one where stochastic variation is greatest [31]). The typical range of q is indicated by the shaded plot region.

https://doi.org/10.1371/journal.ppat.0030045.g003

By sustaining high levels of polymorphism within hosts, a joint increase in selection and transmission rate may also cause the frequency of the escape mutation in the circulating virus population (π = E(p)) to depart substantially from the expected value at equilibrium in the absence of polymorphism ( π˄ = q, i.e., individual virus populations fix alleles matching host HLA genotypes). In our simulations, if selection favoring escape in HLA+ hosts was sufficiently stronger than selection for reversion in HLA− hosts, then π˄ became substantially greater than q at equilibrium (Figure 3B). On the other hand, if selection favoring reversion in HLA− hosts was greater, then the equilibrium value of π˄ was deflected in the opposite direction, below q (not shown). This departure of π˄ from q became more pronounced with increasing transmission rates. Unequal mutation rates between the virus alleles could also contribute to this effect (Figure S3). An escape allele may therefore predominate the circulating virus population even when the selective HLA allele in the host population is rare. In other words, an individual possessing a rare HLA allele may nevertheless stand a high chance of becoming infected by a matched escape variant if selection for reversion is weak and the transmission rate is high.

Deterministic Model of Viral Evolution

This process sustaining high levels of nonsynonymous polymorphism at codon positions under host-specific selection is related to the maintenance of genetic variation in a subdivided population by local adaptation [29,30] and can be illustrated with a simple deterministic model. We use the following differential equation [31]: to describe the mean rate of change in p within a given host, where s is the selection coefficient, and μ and ν are the forward and back mutation rates, respectively. Initial conditions for Equation 1 were defined to reflect the severe bottleneck imposed by transmission of the virus (i.e., p(0) = 0 or p(0) = 1). Assuming that transmission occurs after a constant time interval (τ), the expected value of π after n transmissions is obtained from the recurrence equation: where pHLA+ and pHLA− are approximate solutions of Equation 1 for evolution of p in HLA+ and HLA− hosts, respectively (Protocol S1). Equation 2 has an equilibrium solution: which reduces to π˄ = q when μ = ν and selection for escape and reversion is symmetric between host types (sesc = srev). As τ approaches ∞, π˄ also converges towards q because the evolution of the escape allele within hosts is resolved before transmission (i.e., and ). Conversely, as τ approaches zero, π˄ converges towards a quantity determined by the ratio of ν and μ (Protocol S2). The behavior of π˄ at these limits implies the existence of an intermediate waiting time to transmission (τmax), which maximizes the departure of π˄ from q. An approximation of τmax indicates that it is on the order of max(sesc,srev)−1 when selection is stronger than mutation (Protocol S3). Thus, our model confirms that the greatest departure of π˄ from the expectation q occurs when the mean transmission rate corresponds to the overall intensity of selection.

We found a strong correspondence between this model and simulations (Pearson's ρ = 0.92, p-value < 10−15; Figure S4) with all incongruous cases being caused by stochastic effects due to effective population sizes within hosts of N = 102 or below. The effective population size for HIV-1 is estimated to be on the order of 103 and greater, while the total census population size is typically several orders of magnitude larger [32–34], and the census size for HCV is approximately 10-fold greater still. Hence, this model is a reasonably accurate representation of evolution within realistic HIV-1 and HCV populations.

Discussion

In this study, we have described a novel pattern in the genetic variation of two human viruses, and formulated a simple population genetic model, supplemented with stochastic simulations, to explain it. However, because of the limited availability of population-based sequences that have not been stripped of sequencing mixtures, we were required to restrict our analysis to the RT and PR coding region of HIV-1, in which mixtures provide useful information on the evolution of resistance [21]. Although we focused our investigation on subtype B sequences isolated from treatment-naive individuals, we had no direct control over the sequencing and base-calling conditions of this data set. On the other hand, we obtained unprocessed sequencing electropherogram data of the HCV E1 envelope coding region, such that we could uniformly apply our own methods across all sequences. We were also unable to control for the circumstances under which sequences were isolated from either HIV-1– or HCV-infected patients, e.g., days since infection or seroconversion, regionality of patient populations. Even so, these sampling issues would not bias inferences based on site-by-site comparisons of sequence variation (e.g., relative mixture frequencies). We were able to recover an exceptionally clear and consistent signal of a link between within-host and among-host genetic variation among codon positions in HIV-1 and HCV sequences. This pattern represents strong evidence for CTL-mediated selection in both viruses, specifically targeting with HLA-B–restricted epitopes.

The rapid accumulation of genetic variation in HIV-1 and HCV enables these viruses to elude the immune system and forestalls the development of effective vaccines. Identifying the factors that shape genetic diversity in these human viruses remains a formidable challenge. Because these viruses possess exceptionally high mutation rates, extensive genetic variation accumulates within hosts that may be shaped by ongoing host-specific adaptation. However, the development of models of virus evolution within hosts has been largely independent of “dynamical” models of the transmission and spread of viruses across host cells and individuals [25]. As a result, few models of virus evolution integrate the evolution within hosts with viral dynamics at the level of the host population, which could otherwise reveal emergent properties of evolution within hosts. For example, there is an extensive literature characterizing selection in HIV-1 [10,35–47] by comparing inferred rates of nonsynonymous and synonymous substitutions, but these studies employ methods that do not explicitly distinguish between within- and among-host variation (but see [19,48]).

However, empirical evidence indicates that aspects of the host population can influence patterns of evolution within hosts, and vice versa. For instance, Ross and Rodrigo [10] found evidence that the magnitude and persistence of site-specific diversifying selection within patients was correlated with the rate of progression to acquired immune deficiency syndrome (AIDS), which may influence long-term epidemiological dynamics in the host population. Moore et al. [9] found significant associations between divergent codon positions within CTL epitopes in HIV-1 RT and HLA allelic variation in the host population, which implied that CTL-mediated selection within hosts was influencing the evolution of the total virus population. More recently, Kosakovsky Pond et al. [48] developed a customized phylogenetic analysis to detect significant turnover in codon positions under diversifying selection in HIV-1 PR and RT sequences among human populations with distinct HLA frequencies. They also found that many nonsynonymous substitutions that were mapped to terminal branches of the tree (i.e., occurring within hosts) were absent from internal branches, suggesting that adaptations within individual virus populations were not always maintained at the level of the total virus population [48].

These observations motivate the theoretical development of models of viral evolution that capture the interaction between the within-host and among-host levels of genetic variation. Recently, Grenfell et al. [24] sought to unify the characteristic shape of phylogenetic trees for different virus pathogens with the evolutionary processes within hosts. For instance, phylogenetic trees derived from HIV-1 or HCV sequences sampled from the host population tend to be more “balanced”, reflecting the epidemiological spread of the virus [24]. In contrast, trees derived from influenza A virus hemagglutinin sequences are less balanced, containing a persistent “backbone” that continually spawns short-lived lineages [49]. They proposed that this variation in tree shape, which indirectly manifests the genetic variation among hosts, was driven by the rate at which variants with a selective advantage in the previous host were being transmitted to the subsequent host. Our model complements this previous work by directly evaluating the influence of within-host evolution on the accumulation of nonsynonymous substitutions that differentiate individual virus populations, and the reciprocal effect of this divergence among hosts on variation within hosts. As a result, we can obtain quantitative predictions on how selection within hosts and the transmission rate will influence the frequency of escape variants in the total virus population. The model also predicts that variation in the mean surplus of nonsynonymous mixtures (quantified by the summary statistic mN − mS) per gene indicates divergent intensities of host-specific selection. Similarly, the characteristic transmission rates and overall intensity of selection of different viruses (e.g., HIV-1, HCV, influenza A virus) may revealed by a divergence in the mean surplus of nonsynonymous mixtures per virus. We did not attempt to infer differences between genes or viruses from the absolute frequencies of mixtures in the current data set due to the potential variation in sequencing protocols (as discussed above). Nevertheless, our model should motivate investigators in viral evolution to provide access to raw sequencing data, including annotation of variables that could influence the detection of polymorphisms (e.g., lab sequencing protocol, automated sequencer type and manufacturer).

Based on the distribution of relative mixture frequencies (i.e., site-by-site comparisons within genes), our model indicates that the genetic variation of HIV-1 and HCV is being shaped by the ongoing transmission of unfavorable variants, skewing the frequency of an escape variant in the total virus population towards the direction that host-specific selection is strongest. This unexplored imprint of within-host evolution, manifested as a site-specific surplus of nonsynonymous mixtures within CTL epitopes, can strongly influence the overall composition of the circulating virus population, in addition to founder effects. Because we observed this phenomenon in both HIV-1 and HCV, it may be a common feature of viruses that exhibit both prolific genetic variation within hosts and substantial rates of transmission.

Materials and Methods

HIV-1 and HCV sequence data.

We obtained treatment-naive HIV-1 subtype B sequences from the HIV Drug Resistance Database at Stanford University (Stanford HIVDB) [50]. At the time of analysis, there were 3,458 PR and 1,997 RT sequences meeting our criteria, representing 3,004 and 1,989 patients, respectively. By restricting the data set to treatment-naive individuals, we sought to minimize the confounding effects of selection for drug-resistant variants. Further screening for antiviral resistance was carried out by aligning each sequence to its closest subtype reference sequence (obtained from the Los Alamos National Laboratory [LANL] HIV sequence database; [51]) and scoring for resistance according to the Stanford HIVDB mutation scores using customized scripts in HyPhy [52]. Assuming worst-case resolution of ambiguous nucleotides (i.e., maximized scores), 149 RT and 58 PR sequences with at least low-level resistance (score ≥ 15) were discarded from the data sets. All 297 nucleotide sites from PR sequences were included in our analyses. RT sequences were truncated to nucleotide sites 1 to 741 to exclude poorly sampled tail regions from the analyses.

In addition, we obtained 2,691 chromatogram traces generated from ABI 310 and Beckman CEQ 8000 automated sequencers, covering the core E1 region of HCV. For the majority of traces, each corresponded to a unique isolate from a patient for the initial diagnosis and genotyping of an HCV infection. All trace files were converted to standard chromatogram format and processed with the base-calling program Phred [53]. Potential sequencing mixtures were identified by screening the uncalled peak output using a custom Python script. An uncalled peak was classified as representing a minority variant if: (1) it was located within ±1 trace points of a called peak; and the area under the uncalled peak was (2) at least 20% of the called peak area; (3) at least 10% the mean area of the last ten called peaks; and (4) at least two times greater than the mean area of the last five uncalled peaks. All sequences were truncated to the E1 coding region spanning the nucleotide sites 1 to 399. We also obtained 346 published population-based RT-PCR sequences from Genbank (see Accession Numbers) spanning the hyper-variable region HVR1 of HCV envelope protein E2 [54–57].

Site-specific estimation of substitution rates.

Sequences were aligned using ClustalW [58] and manually adjusted with Se-Al version 2.0 [59] (alignments available upon request). We used neighbor-joining [60] with Tamura-Nei [61] distance to reconstruct the phylogeny from each sequence alignment. Pairwise distances from each phylogeny indicated that the sequences were highly divergent (Figure S5). To estimate the number of nonsynonymous and synonymous substitutions with branch corrections at each codon position, we employed the single-likelihood ancestor counting method [62] as implemented in HyPhy [52,63] using the default settings. Ambiguious nucleotides were resolved to the consensus codon at that position in order to remove any possible influence of mixture frequencies on estimates of substitution rates. We tested for significant positive selection (dN > dS) by applying a continuous extension of the binomial distribution to model the probability that a given proportion of substitutions are nonsynonymous, given the proportion of sites that are nonsynonymous at the codon position [63].

Association with CTL epitopes.

For analyzing associations between nonsynonymous mixture frequencies and epitopes within HIV-1 PR and RT, we applied the CTL epitope definitions from the LANL HIV immunology database [64]. Similarly, we applied the CTL epitope definitions from the LANL HCV immunology database for analyzing associations within HCV E1 and E2 (HVR1) [65].

Simulations of virus evolution.

We implemented a simulation of virus evolution in a host population using an iterative Moran process [66]. Both virus and host populations were each modeled by a single two-allele locus, representing the immune escape and HLA genotypes, respectively. Instantaneous rates for the unit increase and decrease of escape allele frequency within a host were: where j is the number of wild-type alleles in an ideal population of constant size N, and λ1 and λ2 are the wild-type and escape virus growth rates. If the host was HLA−, we set λ1 = 1 and λ2 such that the selection coefficient for reversion srev = (λ1 − λ2) . Otherwise, we set λ1 < λ2 so that sesc = (λ2 − λ1) > 0. After an exponentially distributed waiting time (τ) with rate k, a randomly selected individual from K = 103 hosts was replaced. This new host was HLA+ with probability q (and HLA− otherwise), and infected by wild-type virus with probability jτ/N, where jτ is obtained from the iterative application of j+ and j−, and the total number of events occurring in the time interval τ was determined by random draws from an exponential distribution with the rate (j+ + j−). Otherwise, it was infected by the escape mutant virus. This new infection was therefore initially fixed for either the wild-type or escape virus genotype, assuming a severe bottleneck upon transmission between hosts.

Simulations were run for 200 × K transmissions, which was sufficient for π to converge to an equilibrium for all parameter values evaluated. We recorded the frequency of the escape allele in the individual virus population (p = 1 − j/N), from which we calculated the mean frequency among hosts (π = E(p)). Given the empirical detection threshold of minority variants from population-based sequencing, an individual virus population was considered to be detectably polymorphic if 0.2 < p <0.8. Unique parameter values were assigned to 100 replicate simulations by Latin hypercube sampling from their respective ranges: q = (0, 0.5); μ,ν = (10−5, 10−3); sesc,srev = (0.002, 0.2); N = (102, 105); and k = (0.00137, 0.0137), such that transmissions occur after approximately 0.2 to 2 years (where τ is in units of days).

To compare the simulation results to our deterministic model, we used the numerical integration function in Mathematica 5.1 (Wolfram Research, http://www.wolfram.com) to calculate the expectation of Equation 3 assuming that the waiting time τ was exponentially distributed with rate parameter k.

Supporting Information

Figure S1. Histograms for the Frequency of Nonsynonymous and Synonymous Mixtures per Sequence

The range of frequencies for HCV E2 (HVR1) has been truncated at ten mixtures per sequence for clarity, although a small number of sequences contain as many as 18 mixtures. In HIV-1 PR and RT and HCV E2 (HVR1), there is an excess of mixture-free sequences, possibly due to an under-reporting bias of mixtures which are often interpreted as sequencing errors. HCV E1 sequences were obtained directly from unprocessed trace files and were not subject to this bias. The level of dispersion in the observed frequency distributions was evaluated by fitting Poisson and negative binomial models using a generalized linear models procedure. Goodness-of-fit, quantified by Akaike's information criterion, was improved by the negative binomial model in all cases, and estimates of the dispersion parameter confirmed overdispersion of mixture frequencies in HIV-1 PR and RT and HCV E2 (HVR1).

https://doi.org/10.1371/journal.ppat.0030045.sg001

(13 KB PDF)

Figure S2. Mixture Frequency Distributions

The histograms depict frequency distributions for nonsynonymous (above) and synonymous (below) mixtures per codon position. Note that the histograms for HCV E1 and E2 (HVR1) are on different scales. There is conspicuously greater variation among codon positions in nonsynonymous mixture frequencies, more notably in HIV-1 sequences. Codon positions associated with peaks in the frequency of nonsynonymous mixtures are indicated above each distribution by the alignment consensus amino acid and residue number.

https://doi.org/10.1371/journal.ppat.0030045.sg002

(169 KB PDF)

Figure S3. Effect of Mutation Rate Asymmetry on the Frequency of Escape Mutations

A contour plot depicting the difference π˄ − q as a function of the disparity in selection coefficients and mutation rates between HLA+ and HLA− hosts (r = log10(sesc) − log10(srev) + log10μ − log10ν) and the transmission rate (log10k). When the net effect of mutation and selection is equivalent between HLA+ and HLA− hosts (r = 0), then π˄ converges to q and is independent of variation in transmission rate. In contrast, when there is a net imbalance in mutation and selection (r ≠ 0), there is a departure of π˄ from q; this departure becomes greater with increasing transmission rates. Viral population size has no apparent effect on the difference π˄ − q. Each open circle corresponds to a replicate simulation with unique parameter values set by Latin hypercube sampling.

https://doi.org/10.1371/journal.ppat.0030045.sg003

(30 KB PDF)

Figure S4. Comparison of Stochastic Simulation and Model Predictions

Scatterplot illustrating correspondence between predicted value of π from the deterministic model (x-axis) and values obtained from simulations at equilibrium ( π˄, y-axis). A solid line is drawn at x = y to indicate an exact match between model and simulation frequencies. Dashed lines above and below the x = y axis enclose variation in frequencies within a ±10% interval. Disparity between the model and simulations is caused by a lack of stochastic factors in the model. Replacing the unidirectional mutation approximations of the model (pHLA+ and pHLA−) by the exact formula has no visible effect on the correspondence between the model and simulations.

https://doi.org/10.1371/journal.ppat.0030045.sg004

(9 KB PDF)

Figure S5. Frequency Distributions of Pair-Wise Distances of HIV-1 and HCV Sequences

There is a substantial amount of divergence among the vast majority of HIV-1 sequences in the reconstructed phylogenetic trees, with only <0.1% of pairwise distances below 0.01. This is consistent with the low number of HIV-1 PR and RT sequences that were re-sampled from the same patient. HCV E1 and E2 (HVR1) sequences were highly divergent on average. A small proportion of pairwise distances between HCV E1 sequences (1.1%), particularly in subtype 4d, were below 0.05. Similarly, about 3% of pairwise distances between HCV E2 (HVR1) sequences were below a threshold of 0.25. Hence, a minority of HCV sequences may have represented multiple isolates from patients, but were too few overall to influence the outcome of our analyses.

https://doi.org/10.1371/journal.ppat.0030045.sg005

(14 KB PDF)

Protocol S1. Approximation of Allele Frequency Evolution

https://doi.org/10.1371/journal.ppat.0030045.sd001

(73 KB PDF)

Protocol S2. Limit Behavior of Deterministic Model

https://doi.org/10.1371/journal.ppat.0030045.sd002

(41 KB PDF)

Protocol S3. Approximation of Optimal Waiting Time to Transmission

https://doi.org/10.1371/journal.ppat.0030045.sd003

(60 KB PDF)

Accession Numbers

GenBank (http://www.ncbi.nlm.nih.gov/Genbank/index.html) accession numbers for the HCV E1 envelope protein-coding sequences used in our study are AY766700–AY768365. GenBank accession numbers for the E2 envelope protein (HVR1) sequences used in our study are the following: AY390002, AY390005, AY390008, AY390010, AY390013, AY390016, AY390019, AY390022, AY390024, AY390027, AY390030, AY390032, AY742960–AY743049, AY309923–AY309954, AY314963–AY314969, AY390002–AY390035, AY564735–AY564784, and AY935999–AY936132.

Acknowledgments

We also wish to thank Robert Shafer and the Stanford HIV Drug Resistance Database team for the development and maintenance of an excellent resource for HIV resistance analyses.

Author Contributions

DDR, AJLB, and SDWF conceived and designed the experiments. AFYP performed the experiments and wrote the paper. AFYP and SLKP analyzed the data. AFYP, SLKP, and PB contributed reagents/materials/analysis tools.

References

  1. 1. Allen TM, Altfeld M, Geer SC, Kalife ET, Moore C, et al. (2005) Selective escape from CD8+ T-cell responses represents a major driving force of human immunodeficiency virus type 1 (HIV-1) sequence diversity and reveals constraints on HIV-1 evolution. J Virol 79: 13239–13249.
  2. 2. Bowen DG, Walker CM (2005) Adaptive immune responses in acute and chronic hepatitis C virus infection. Nature 436: 946–952.
  3. 3. Klein J (1986) Natural history of the major histocompatibility complex. New York: Wiley. 775 p.
    • 4. Pircher H, Moskkophidis D, Rohrer U, Bürki K, Hengartner H, et al. (1990) Viral escape by selection of cytotoxic T cell-resistant virus variants in vivo. Nature 346: 629–633.
    • 5. Yusim K, Kesmir C, Gaschen B, Addo MM, Altfeld M, et al. (2002) Clustering patterns of cytotoxic T-lymphocyte epitopes in human immunodeficiency virus type 1 (HIV-1) proteins reveal imprints of immune evasion on HIV-1 global variation. J Virol 76: 8757–8768.
    • 6. Fernandez CS, Stratov I, De Rose R, Walsh K, Dale CJ, et al. (2005) Rapid viral escape at an immunodominant simian-human immunodeficiency virus cytotoxic T-lymphocyte epitope exacts a dramatic fitness cost. J Virol 79: 5721–5731.
    • 7. Friedrich TC, Dodds EJ, Yant LJ, Vojnov L, Rudersdorf R, et al. (2004) Reversion of CTL escape-variant immunodeficiency viruses in vivo. Nature Med 10: 275–281.
    • 8. Leslie AJ, Pfafferott KJ, Chetty P, Draenert R, Addo MM, et al. (2004) HIV evolution: CTL escape mutation and reversion after transmission. Nature Med 10: 282–289.
    • 9. Moore CB, John M, James IR, Christiansen FT, Witt CS, et al. (2002) Evidence of HIV-1 adaptation to HLA-restricted immune responses at a population level. Science 296: 1439–1443.
    • 10. Ross HA, Rodrigo AG (2002) Immune-mediated positive selection drives human immunodeficiency virus type 1 molecular variation and predicts disease duration. J Virol 76: 11715–11720.
    • 11. Gaudieri S, Rauch A, Park LP, Freitas E, Herrmann S, et al. (2006) Evidence of viral adaptation to HLA class I-restricted immune pressure in chronic hepatitis C virus infection. J Virol 80: 11094–11104.
    • 12. Trachtenberg E, Korber B, Sollars C, Kepler TB, Hraber PT, et al. (2003) Advantage of a rare HLA supertype in HIV disease progression. Nature Med 9: 928–935.
    • 13. Goulder P, Brander C, Tang YH, Tremblay C, Colbert RA, et al. (2001) Evolution and transmission of stable CTL escape mutations in HIV infection. Nature 412: 334–338.
    • 14. Barouch DH, Powers J, Truitt DM, Kishko MG, Arthur JC, et al. (2005) Dynamic immune responses maintain cytotoxic T lymphocyte epitope mutations in transmitted simian immunodeficiency virus variants. Nature Immunol 6: 247–252.
    • 15. Shankarappa R, Margolick JB, Gange SJ, Rodrigo AG, Upchurch D, et al. (1999) Consistent viral evolutionary changes associated with the progression of human immunodeficiency virus type 1 infection. J Virol 73: 10489–502.
    • 16. Rouzine IM, Coffin JM (1999) Search for the mechanism of genetic variation in the pro gene of human immunodeficiency virus. J Virol 73: 8167–8178.
    • 17. Kimura M (1968) Evolutionary rate at the molecular level. Nature 217: 624–626.
    • 18. Hudson RR, Kreitman M, Aguadé M (1987) A test of neutral molecular evolution based on nucleotide data. Genetics 116: 153–159.
    • 19. Williamson S (2003) Adaptation in the env gene of HIV-1 and evolutionary theories of disease progression. Mol Biol Evol 20: 1318–1325.
    • 20. Smith NGC, Eyre-Walker A (2002) Adaptive protein evolution in Drosophila. Nature 415: 1022–1024.
    • 21. Larder BA, Kohli A, Kellam P, Kemp SD, Kronick M, et al. (1993) Quantitative detection of HIV-1 drug resistance mutations by automated DNA sequencing. Nature 365: 671–673.
    • 22. Leitner T, Halapi E, Scarlatti G, Rossi P, Albert J, et al. (1993) Analysis of heterogeneous viral populations by direct DNA-sequencing. Biotechniques 15: 120–127.
    • 23. Schuurman R, Demeter L, Reichelderfer P, Tijnagel J, de Groot T, et al. (1999) Worldwide evaluation of DNA sequencing approaches for identification of drug resistance mutations in the human immunodeficiency virus type 1 reverse transcriptase. J Clin Microbiol 37: 2291–2296.
    • 24. Grenfell BT, Pybus OG, Gog JR, Wood JLN, Daly JM, et al. (2004) Unifying the epidemiological and evolutionary dynamics of pathogens. Science 303: 327–332.
    • 25. Kelly JK, Williamson S, Orive ME, Smith MS, Holt RD (2003) Linking dynamical and population genetic models of persistent viral infection. Am Nat 162: 14–28.
    • 26. Benjamini Y, Hochberg Y (1995) Controlling the false discovery rate: A practical and powerful approach to multiple testing. J R Stat Soc B 57: 289–300.
    • 27. Kiepiela P, Leslie AJ, Honeyborne I, Ramduth D, Thobakgale C, et al. (2004) Dominant influence of HLA-B in mediating the potential co-evolution of HIV and HLA. Nature 432: 769–774.
    • 28. John M, Mallal S (2005) CTL responses to HIV and SIV: Wrestling with smoke. Nature Immunol 6: 232–234.
    • 29. Moran PAP (1962) The statistical processes of evolutionary theory. Oxford: Clarendon Press.
      • 30. Felsenstein J (1976) The theoretical population genetics of variable selection and migration. Ann Rev Genet 10: 253–280.
      • 31. Crow JF, Kimura M (1970) An introduction to population genetics theory. New York: Harper and Row.
        • 32. Leigh Brown AJ (1997) Analysis of HIV-1 env gene sequences reveals evidence for a low effective number in the viral population. Proc Natl Acad Sci USA 94: 1862–1865.
        • 33. Nijhuis M, Boucher CAB, Schipper P, Leitner T, Schuurman R, et al. (1998) Stochastic processes strongly influence HIV-1 evolution during suboptimal protease-inhibitor therapy. Proc Natl Acad Sci USA 95: 14441–14446.
        • 34. Achaz G, Palmer S, Kearney M, Maldarelli F, Mellors JW, et al. (2004) A robust measure of HIV-1 population turnover within chronically infected individuals. Mol Biol Evol 21: 1902–1912.
        • 35. Bonhoeffer S, Holmes EC, Nowak MA (1995) Causes of HIV diversity. Nature 376: 125.
        • 36. Liu SL, Schacker T, Musey L, Shriner D, McElrath MJ, et al. (1997) Divergent patterns of progression to AIDS after infection from the same source: Human immunodeficiency virus type 1 evolution and antiviral responses. J Virol 71: 4284–4295.
        • 37. Yamaguchi Y, Gojobori T (1997) Evolutionary mechanisms and population dynamics of the third variable region of HIV within single hosts. Proc Natl Acad Sci USA 94: 1264–1269.
        • 38. Yamaguchi-Kabata Y, Gojobori T (2000) Reevaluation of amino acid variability of the human immunodeficiency virus type 1 gp120 envelope glycoprotein and prediction of new discontinuous epitopes. J Virol 74: 4335–4350.
        • 39. Chen L, Perlina A, Lee CJ (2004) Positive selection detection in 40,000 human immunodeficiency virus (HIV) type 1 sequences automatically identifies drug resistance and positive fitness mutations in HIV protease and reverse transcriptase. J Virol 78: 3722–3732.
        • 40. Simmonds P, Balfe P, Ludlam C, Bishop J, Leigh Brown A (1990) Analysis of sequence diversity in hypervariable regions of the external glycoprotein of human immunodeficiency virus type 1. J Virol 64: 5840–5850.
        • 41. Seibert SA, Howell CY, Hughes MK, Hughes AL (1995) Natural selection on the gag, pol, and env genes of human immunodeficiency virus 1 (HIV-1). Mol Biol Evol 12: 803–813.
        • 42. Nielsen R, Yang Z (1998) Likelihood methods for detecting positively selected amino acid sites and applications to the HIV-1 envelope gene. Genetics 148: 929–936.
        • 43. Crandall KA, Kelsey CR, Imamichi H, Lane HC, Salzman NP (1999) Parallel evolution of drug resistance in HIV: Failure of nonsynonymous/synonymous substitution rate ratio to detect selection. Mol Biol Evol 16: 372–382.
        • 44. Cichutek K, Merget H, Norley S, Linde R, Kreuz W, et al. (1992) Development of a quasispecies of human immunodeficiency virus type 1 in vivo. Proc Natl Acad Sci USA 89: 7365–7369.
        • 45. Wolfs TFW, Zwart W, Bakker M, Goudsmit J (1992) HIV-1 genomic RNA diversification following sexual and parenteral virus transmission. Virology 189: 103–110.
        • 46. McDonald RA, Mayers DL, Chung RC, Wagner KF, Ratto-Kim S, et al. (1997) Evolution of human immunodefiency virus type 1 env sequence variation in patients with diverse rates of disease progression and T-cell function. J Virol 71: 1871–1879.
        • 47. Zanotto PMA, Kallas EG, de Souza RF, Holmes EC (1999) Genealogical evidence for positive selection in the nef gene of HIV-1. Genetics 153: 1077–1089.
        • 48. Kosakovsky Pond SL, Frost SDW, Grossman Z, Gravenor MB, Richman DD, et al. (2006) Adaptation to different human populations by HIV-1 revealed by codon-based analyses. PLoS Comput Biol 2: e62.
        • 49. Ferguson NM, Galvani AP, Bush RM (2003) Ecological and immunological determinants of influenza evolution. Nature 422: 428–433.
        • 50. Rhee SY, Gonzales MJ, Kantor R, Betts BJ, Jaideep R, et al. (2003) Human immunodeficiency virus reverse transcriptase and protease sequence database. Nucl Acids Res 31: 298–303.
        • 51. Leitner T, Korber B, Daniels M, Calef C, Foley B (2005) HIV-1 subtype and circulating recombinant form (CRF) reference sequences, 2005. In: Leitner T, Foley B, Hahn B, Marx P, McCutchan F, et al., editors. HIV Sequence Compendium 2005. Los Alamos (New Mexico): Los Alamos National Laboratory.
          • 52. Kosakovsky Pond SL, Frost SDW, Muse SV (2005) HyPhy: Hypothesis testing using phylogenies. Bioinformatics 21: 676–679.
          • 53. Ewing B, Hillier L, Wendl MC, Green P (1998) Base-calling of automated sequencer traces using Phred. I. Accuracy assessment. Genome Res 8: 175–185.
          • 54. Alfonso V, Flichman DM, Sookoian S, Mbayed VA, Campos RH (2004) Evolutionary study of HVR1 of E2 in chronic hepatitis C virus infection. J Gen Virol 85: 39–46.
          • 55. Ben Othman S, Bouzgarrou N, Achour A, Bourlet T, Pozzetto B, et al. (2004) High prevalence and incidence of hepatitis C virus infections among dialysis patients in the East-Centre of Tunisia. Pathol Biol 52: 323–327.
          • 56. Nicot F, Legrand-Abravanel F, Sandres-Saune K, Boulestin A, Dubois M, et al. (2005) Heterogeneity of hepatitis C virus genotype 4 strains circulating in south-western France. J Gen Virol 86: 107–114.
          • 57. Pasquier C, Njouom R, Ayouba A, Dubois M, Sartre MT, et al. (2005) Distribution and heterogeneity of hepatitis C genotypes in hepatitis patients in Cameroon. J Med Virol 77: 390–398.
          • 58. Thompson JD, Higgins DG, Gibson TJ (1994) CLUSTAL W: Improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. Nucl Acids Res 22: 4673–4680.
          • 59. Rambaut A (1996) Se-Al: Sequence alignment editor. Available: http://evolve.zoo.ox.ac.uk. Accessed 27 February 2007.
            • 60. Saitou N, Nei M (1987) The neighbor-joining method: A new method for reconstructing phylogenetic trees. Mol Biol Evol 4: 406–425.
            • 61. Tamura K, Nei M (1993) Estimation of the number of nucleotide substitutions in the control region of mitochondrial DNA in humans and chimpanzees. Mol Biol Evol 10: 512–526.
            • 62. Suzuki Y, Gojobori T (1999) A method for detecting positive selection at single amino acid sites. Mol Biol Evol 16: 1315–1328.
            • 63. Kosakovsky Pond SL, Frost SDW (2005) Not so different after all: A comparison of methods for detecting amino acid sites under selection. Mol Biol Evol 22: 1208–1222.
            • 64. Korber BTM, Brander C, Haynes BF, Koup R, Moore JP, et al. (2005) HIV Molecular Immunology 2005. Los Alamos (New Mexico): Los Alamos National Laboratory.
              • 65. Yusim K, Richardson R, Tao N, Szinger J, Funkhouser R, et al. (2005) The Los Alamos hepatitis C immunology database. Appl Bioinformatics 4: 217–225.
              • 66. Moran PAP (1958) Random processes in genetics. Math Proc Camb Phil Soc 54: 60–71.