Skip to main content


Selection dramatically reduces effective population size in HIV-1 infection

Article metrics

  • 3571 Accesses

  • 15 Citations



In HIV-1 evolution, a 100–100,000 fold discrepancy between census size and effective population size (N e ) has been noted. Although it is well known that selection can reduce N e , high in vivo mutation and recombination rates complicate attempts to quantify the effects of selection on HIV-1 effective size.


We use the inbreeding coefficient and the variance in allele frequency at a linked neutral locus to estimate the reduction in N e due to selection in the presence of mutation and recombination. With biologically realistic mutation rates, the reduction in N e due to selection is determined by the strength of selection, i.e., the stronger the selection, the greater the reduction. However, the dependence of N e on selection can break down if recombination rates are very high (e.g., r ≥ 0.1). With biologically likely recombination rates, our model suggests that recurrent selective sweeps similar to those observed in vivo can reduce within-host HIV-1 effective population sizes by a factor of 300 or more.


Although other factors, such as unequal viral reproduction rates and limited migration between tissue compartments contribute to reductions in N e , our model suggests that recurrent selection plays a significant role in reducing HIV-1 effective population sizes in vivo.


The effective population size, N e , is defined as the size of an idealized population that has the same population genetics properties (generally those properties that measure the magnitude of random genetic drift) as the actual population. Most studies have estimated the within-host N e for HIV-1 during chronic infection to be ~103 [15], though one study estimated N e to be between 105 and 5 × 105 [6]. Even the highest of these estimates is about two orders of magnitude lower than the number of productively infected cells, estimated to be on the order of 107 to 108 cells [7]. Explanations for low N e values include unequal viral reproduction rates [25, 8], structured populations [812], and recurrent selection [25, 8]. The possibility that recurring selection may be reducing viral diversity is unsettling because most of the computational models used to estimate N e assume neutral evolution.

During a selective sweep of a favorable allele, any neutral alleles linked to the selected allele will rise in frequency and become overrepresented in the population. This process, called "hitchhiking", can reduce neutral diversity more than random genetic drift and therefore reduce N e [13]. Although selection has been acknowledged as a possible explanation for the low within-host effective population size during chronic HIV-1 infection [3, 12], high mutation [14, 15] and recombination rates [1620] complicate attempts to study the effects of selection on HIV-1 in vivo. To address these issues, we extended a classic "inbreeding coefficient" method [2123] to derive recurrence equations that account for the combined effects of selection, mutation, and recombination. We then used these equations to quantify the effects of selection on effective size using parameters relevant to HIV-1 evolution in vivo.

Results and Discussion

Overview of the genetic model

Our model follows the basic Wright-Fisher assumptions of a single haploid population of constant size with no subdivision or migration, non-overlapping generations, and random sampling of offspring each generation. We calculated N e in terms of the inbreeding effective size, which is based on the change of the average inbreeding coefficient (F) at a neutral locus (L) that is linked to a locus (S) that is under selection. The inbreeding coefficient is defined as the probability that two individuals are identical by descent (which means they are identical and have a common ancestor). Therefore, for the neutral locus L, two individuals are identical by descent if they are derived from a common ancestor and are identical at locus L, regardless of the status of locus S. Our approach to estimating N e was to determine changes in the inbreeding coefficient at the neutral locus in the presence and absence of selection and recombination. The effective population size was defined as the size of the neutral population that gave changes in the inbreeding coefficient that were equal to those observed in the presence of selection and recombination.

As shown in Figure 1A, in the absence of recombination, an offspring can be derived from a parent in the previous generation with either allele a or A at locus S. An offspring with allele a can be derived by two pathways: from a parent with allele a (without mutation) or a parent with allele A (with an A to a mutation). An offspring with allele A can be derived by two similar pathways. Therefore, F t (the value of F at time t) will be the sum of the probability that two offspring are derived from a certain combination of parents (both with allele A, both with allele a, and one with allele A and the other with allele a) times the probability that the offspring are identical by descent at locus L (see Appendix).

Figure 1

Illustration of the genetic model. A) In the absence of recombination, an offspring with allele a can be derived from a parent with allele a (without mutation) or a parent with allele A (with an A to a mutation); an offspring with allele A can be derived from a parent with allele A (without mutation) or a parent with allele a (with an a to A mutation). B) In the presence of recombination, an offspring with allele a at locus S can be derived from parent(s) in the previous generation by four pathways: 1) Locus S from a parent with allele a without mutation or recombination, (or with recombination between another parent with allele a). 2) Locus S from a parent with allele A following an A to a mutation but no recombination (or with recombination between another parent with allele A). 3) Locus S from a parent with allele a without mutation, but with recombination between another parent with allele A. 4) Locus S from a parent with allele A following an A to a mutation and recombination between another parent with allele a. An offspring with allele A can be derived from parent(s) in four pathways similar to those described above. For the purpose of illustration, only 8 genomes were presented in generation t-1 and t.

In the presence of recombination, loci L and S can be derived from different parents (Figure 1B). An offspring with allele a or A at locus S can be derived from one or more parents in the previous generation by the four pathways illustrated in Figure 1B. As above, F t will be the sum of the probability that the two offspring are derived from a certain combination of pathways (both having locus S from parents with allele A, both having locus S from parents with allele a, one having locus S from a parent with allele a and the other having locus S from a parent with allele A) times the probability that the offspring derived from these pathways are identical by descent at locus L (see APPENDIX).

Effect of selection on effective population size

We used the ratio N/N e to summarize the reduction in N e due to selection from the start of selection at t = 0 until t = t nearlyfixed [the time when the frequency of the advantageous allele reaches (N-1)/N]. This last approximation is helpful because fixation time is asymptotic, with the advantageous allele never reaching 100% in a deterministic model.

In the absence of mutation, the reduction in N e due to selection was most strongly affected by the initial frequency of the advantageous allele, A0 (Figure 2A). In the presence of mutation, the reduction in N e due to selection was most sensitive to the selective advantage, s, of the advantageous allele (Figure 2B). Indeed, for a homogeneous population of N = 107, the N/N e ratio increased 6 – 9 fold with each 10-fold increase in the selection coefficient in the presence of mutation. However, recombination can break the hitchhiking effect of selection on N e (Figure 2C). For example, when r ≤ 10-3, for locus L with U = μ, selection with s = 0.1 reduced N e by ~20 fold. In contrast, when r ≥ 0.1, selection with s = 0.1 had little effect on N e .

Figure 2

Reduction in N e following a selective replacement with and without mutations in loci L and S. A) The effect of different initial frequencies of the advantageous allele. B) The effect of different selection coefficients. C) The effect of different recombination rates. D) The effect of different mutation rates at locus L. E) The effect of different initial inbreeding coefficients. F) The effect of different census population sizes. Panels A and C-F all assume s = 0.1. Solid lines indicate that the N/N e ratios are based on the inbreeding coefficient F t ; dashed lines indicate that the N/N e ratios are based on the variance effective population sizes estimated from our simulations. In the presence of mutation, the dashed lines indicate the N/N e ratios based on the upper and lower estimates of variance effective populations size. Black lines indicate cases with mutation; grey lines indicate cases without mutation. Unless otherwise specified, the following parameters were used: in the absence of mutation, μ = 0, v = 0, U = 0, N = 107, s = 0.1, A0 = 10-6, r = 0, and F0 = FAA,0 = Faa,0 = FAa,0 = 0.1; in the presence of mutation, μ = 2.5 × 10-5, v = μ/3, U = μ, N = 107, s = 0.1, A0 = 0, r = 0, and F0 = Faa,0 = 1, FAA,0 = FAa,0 = 0.

Effective population sizes calculated from the inbreeding coefficient (inbreeding N e ) are usually the same as those calculated from the variance in the allele frequency (variance N e ), though exceptions do occur [2427]. To validate our results, we estimated the effect of selection on N e by calculating the variance in the frequency of the linked neutral allele from simulations using the same genetic model. Values for the inbreeding N e obtained from the calculations above were generally consistent with the estimates of the variance N e derived from these simulations (Figure 2A to 2C). We noted that there was an approximate 3-fold difference in the N e values between the two methods when s = 0.01 (Figure 2B). This is likely due to the fact that the inbreeding N e was estimated using a strict deterministic model; while the variance N e was estimated from simulations of s = 0.01, where genetic drift plays a bigger role.

A very high mutation rate at the neutral locus L (e.g., U = 1000μ) also diminished the reduction in N e due to selection (Figure 2D). In the absence of mutation, the effect of selection was insensitive to changes in the initial homogeneity at locus L (Figure 2E). In the presence of mutation, selection with an initially heterogeneous population at locus L caused greater reductions in N e than selection with an initially homogeneous population. For F0 less than 0.1, however, further increases in the initial heterogeneity (i.e., making F0 even lower) did not lead to further reductions in N e through selection. Interestingly, reductions in N e /N due to selection were insensitive to changes in the census population size, N (Figure 2F).

Effect of recurrent selection on effective population size

For a homogeneous population under recurrent selection, the inbreeding coefficient of the neutral allele decreased until it reached a quasi-steady state, where it fluctuated in a regular "sawtooth" fashion (Figure 3A). The effect of recurrent selection on N e was sensitive to selection strength. For example, for a homogeneous population of N = 107 and U = μ, the decline of F t over time under recurrent selection with s = 0.01 overlapped the neutral curve when N = 9,973,000, while the decline of F t under recurrent selection with s = 0.1 overlapped the neutral curve when N = 28,220 (Figure 3A). In other words, recurrent selection had little effect on N e when s = 0.01, while recurrent selection reduced N e by over 300-fold when s = 0.1 (Figure 3A and 3B). This reduction in N e by recurrent selection could be diminished by high recombination rates (Figure 3B). Although recombination had little impact on the reduction in N e due to selection under a model with s = 0.1 and r ≤ 10-3, with r ≥ 0.1, recombination completely broke the hitchhiking effects of selection on N e with s = 0.1.

Figure 3

Reduction in N e due to recurrent selection. A) The changes of F t over time under selection (left panels) and under neutrality (right panel), in the absence of recombination. B) The effects of different selection coefficients (left panel, r = 0) and different recombination rates (right panel, s = 0.1). The starting parameters were: F0 = Faa,0 = 1, FAA,0 = FAa,0 = 0, A0 = 0, U = μ = 2.5 × 10-5, v = μ/3.


We examined the combined effects of selection, mutation, and recombination on the effective population size of a neutral locus that is linked to a locus under selection. Consistent with other studies [2123], we found that selection can increase the inbreeding coefficient and reduce the inbreeding effective population size. Without mutation, this reduction is primarily determined by the initial frequency of the advantageous allele, i.e., the lower the initial frequency, the greater the effect. With mutation, this reduction is mostly determined by the strength of selection, i.e., the stronger the selection, the greater the effect. With moderate recombination rates (e.g., r ≤ 10-3), recurrent selection can substantially lower N e , though the hitchhiking effect disappears if the recombination rates are very high (e.g., r ≥ 0.1).

The effective population size of HIV-1 during chronic infection has been shown to be 100- to 100,000-fold lower than within-host census size. Indeed, CTL responses are a driving force of HIV-1 evolution and these responses continuously select for escape mutants during chronic infection [2830]. In a comprehensive study of viral evolution and CTL responses during the first four years of HIV-1 infection in one subject, Liu et al. [30, 31] found that of the 25 epitopes detected in this subject, 17 were largely replaced by mutants over time. The selection coefficients for the CTL escape mutant(s) of a single epitope ranged from 0.2 to 0.4 during acute infection and from 0.0024 to 0.15 during chronic infection, with an average of 0.03 [30]. Humoral and escape-specific CTL responses impose additional selective pressures not quantified in Liu et al. [30, 31]. With low to moderate recombination rates, our model shows that recurrent selection with s = 0.1 reduces viral effective population size by approximately 300-fold. Therefore, during HIV-1 infection, selection alone is likely to reduce the viral effective population size to an N e of ~105. This result is close to the estimate of N e ~5 × 105 that Rouzine and Coffin [6] obtained from a model that accounts for selection. The small discrepancy may be due to their use of a lower mutation rate (10-5 vs. 2.5 × 10-5 in our study) and possible biased sampling of sites with higher underlying mutation rates in their study [5].

With high recombination rates, our model predicts that selection has little effect on N e . Observations of 3 to 13 cross-over events per virion in vitro [1720] suggest an intrinsic recombination rate of 10-4 to 10-3 per adjacent site per generation. However, this range is not relevant to our model since these estimates were obtained using heterozygous virions, which may not be abundant in vivo. While Jung and colleagues [32] have demonstrated that cells in the spleen are infected with multiple viruses (a pre-requisite for the formation of heterozygous virions), they did not determine how often heterozygous virions are formed. More relevant is data in which SCID-HU mice were infected with a 50:50 mixture of two marked strains [19]. Two-to-three weeks after infection, an average of ~0.01% of infected cells carried a phenotypic marker of recombination (present on half of all recombinants). Conservatively assuming a single generation of recombination, we estimate from equation (11, Appendix) that the probability of recombination between their two markers (which were 408 bp apart) was r = ~ p Aa /(p A p a ) = 0.0001/(0.5 × 0.5) = ~4 × 10-4 per virion per generation – a value too low to break the hitchhiking effects of selection in our model. However, we recognize these are approximate values obtained from a somewhat artificial system. HIV-1 evolution studies could benefit from additional studies of marked viruses in animal models and clever retrospective analyses of in vivo data from humans to determine evolutionarily relevant recombination rates.


Genetic model

We assume a Wright-Fisher model with a neutral locus L that is linked to a locus under selection, locus S. The selected locus has two alleles, an advantageous allele, A, with a fitness w = 1 + s, and a disadvantageous allele, a, with a fitness of 1. Allele A mutates to a at rate μ and allele a mutates to A at rate ν, while neutral mutations at locus L occur at rate U. A description of all the characters, parameters, and variables used in this study is listed in Table 1. For the purposes of calculation, we assume the following parameters are known: the initial frequency of allele A (A0); the initial frequency of allele a (a0); and the initial inbreeding coefficient at locus L among all individuals (F0), among individuals with allele A (FAA,0), among individuals with allele a (Faa,0), and between individuals with allele A and those with allele a (FAa,0).

Table 1 Description of characters, parameters and variables.

Parameters for HIV-1

The average mutation rate of HIV-1 has been estimated to be 2.5 × 10-5 per nucleotide per generation [14], although one recent study estimated a higher mutation rate of ~8.5 × 10-5 per site per generation [15]. Assuming that any nucleotide substitution at a defined nucleotide site shifts locus S from the advantageous to the disadvantageous state, we defined μ = 2.5 × 10-5 per generation. Assuming that only a particular nucleotide substitution at this site increases fitness, we set ν = μ/3. Since the census sizes of productively HIV-1 infected cells in vivo exceeds 107 [7, 33], most of the comparisons in this study were with N = 107. Since the accumulation of advantageous alleles in populations is more stochastic as N decreases, we only examined populations with N ≥ 106.

Effect of selection on effective population size without mutation

Under selection, the inbreeding coefficient of the linked neutral locus will increase faster than expected by random genetic drift until the selected advantageous allele is fixed (A t = 100%). Because we are using a deterministic model, fixation time is asymptotic. To quantify the effect of selection, we determined the average time for an advantageous allele to approach fixation, t nearlyfixed , and the value of F at t nearlyfixed . t nearlyfixed can be calculated from t = log ( A t a 0 a t A 0 ) / log ( w ) MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xH8viVGI8Gi=hEeeu0xXdbba9frFj0xb9qqpG0dXdb9aspeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGaemiDaqNaeyypa0JagiiBaWMaei4Ba8Maei4zaCMaeiikaGscfa4aaSaaaeaacqWGbbqqdaWgaaqaaiabdsha0bqabaGaemyyae2aaSbaaeaacqaIWaamaeqaaaqaaiabdggaHnaaBaaabaGaemiDaqhabeaacqWGbbqqdaWgaaqaaiabicdaWaqabaaaaOGaeiykaKIaei4la8IagiiBaWMaei4Ba8Maei4zaCMaeiikaGIaem4DaCNaeiykaKcaaa@46DF@ , where t is the time just before the favored allele A at locus S becomes fixed; i.e., when A t = N 1 N MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xH8viVGI8Gi=hEeeu0xXdbba9frFj0xb9qqpG0dXdb9aspeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGaemyqae0aaSbaaSqaaiabdsha0bqabaGccqGH9aqpjuaGdaWcaaqaaiabd6eaojabgkHiTiabigdaXaqaaiabd6eaobaaaaa@3452@ and a t = 1 N MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xH8viVGI8Gi=hEeeu0xXdbba9frFj0xb9qqpG0dXdb9aspeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaqcfaOaemyyae2aaSbaaeaacqWG0baDaeqaaiabg2da9maalaaabaGaeGymaedabaGaemOta4eaaaaa@326B@ . F was calculated with μ = 0, v = 0, and U = 0. The corresponding N e is defined here as the population size under neutrality that will increase F from F0 to F t n e a r l y f i x e d MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xH8viVGI8Gi=hEeeu0xXdbba9frFj0xb9qqpG0dXdb9aspeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGaemOray0aaSbaaSqaaiabdsha0naaBaaameaacqWGUbGBcqWGLbqzcqWGHbqycqWGYbGCcqWGSbaBcqWG5bqEcqWGMbGzcqWGPbqAcqWG4baEcqWGLbqzcqWGKbazaeqaaaWcbeaaaaa@3DD8@ between t = 0 and t = t nearlyfixed . We determined the corresponding N e under the following conditions: N = 106 to 109; s = 0.01 to 10; A0 = 10-7 to 10-3; and F0 = FAA,0 = Faa,0 = FAa,0 = 10-4 to 0.8 (if F0 = 1, F will not change over time without mutation, regardless of selection).

Effect of selection on effective population size with mutation and recombination

The frequency of the A allele cannot be maintained at 100% with the occurrence of the back mutation from A to a at locus S. Therefore t nearlyfixed was set to the time that A t and a t reached equilibrium, i.e., when A t = At+1. The corresponding N e , the population size under neutrality that will increase F from F0 to F t n e a r l y f i x e d MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xH8viVGI8Gi=hEeeu0xXdbba9frFj0xb9qqpG0dXdb9aspeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGaemOray0aaSbaaSqaaiabdsha0naaBaaameaacqWGUbGBcqWGLbqzcqWGHbqycqWGYbGCcqWGSbaBcqWG5bqEcqWGMbGzcqWGPbqAcqWG4baEcqWGLbqzcqWGKbazaeqaaaWcbeaaaaa@3DD8@ between t = 0 and t = t nearlyfixed , was determined using numerical iteration [Appendix equation (2)]. We determined the corresponding N e under the following conditions: N = 106 to 109; s = 0.01 to 10; A0 = 0 to 10-3; F0 = Faa,0 = 10-4 to 1; FAA,0 = FAa,0 = 0, if A0 = 0 and FAA,0 = FAa,0 = F0, if A0 > 0; μ = 2.5 × 10-5, v = μ/3, U = μ to 1000μ, and r = 0 to 1. With these high advantageous mutation rates and large population sizes (Nv >> 1), individuals with allele a had mutations to allele A in almost every generation, preventing advantageous allele A from being lost from the population due to genetic drift.

Effect of recurrent selection on effective population size

With the fixation of the advantageous allele A, the inbreeding coefficient of locus L will undergo a nearly neutral change unless new alleles linked to locus L become advantageous. To estimate the effect of recurrent selection on N e , we assumed that all loci under selection are linked to locus L in the absence of recombination. We also assumed that each selected locus was under sequential selection, i.e., when the frequency of an advantageous allele reached 99.9% at generation t, we assumed that another locus started to undergo selection (calculated by setting A t = 0, Faa, t= F t , FAA, t= 0, and FAa, t= 0). For simplicity, we assumed that all of the selected loci have the same mutation rate and selection coefficient. We calculated F under recurrent selection under the following conditions: N = 107, A0 = 0, FAA,0 = FAa,0 = 0, F0 = Faa,0 = 1, μ = 2.5 × 10-5, v = μ/3, and U = μ; s = 0.01 to 10; and r = 0 to 1.

Estimate of the effect of selection on variance effective population size by simulation

The change in the average inbreeding coefficient is one of several criteria used to estimate effective population size [2427]. To validate our results using a different measure of effective population size, we estimated the effect of selection on N e by calculating the variance in the frequency of the linked neutral allele from simulations using the genetic model described above. The parameters used in these simulations were the same as those used for the calculation for the inbreeding coefficient described above. When simulating selection in the absence of mutation, the simulations were performed under the following conditions: N = 107; s = 0.01 to 10; A0 = 10-7 to 10-3; F0 = FAA,0 = Faa,0 = FAa,0 = 0.1; μ = v = U = 0; and r = 0. When simulating selection in the presence of mutation, the simulations were performed with the following conditions: N = 107; s = 0.01 to 10; A0 = 0; F0 = Faa,0 = 1, FAA,0 = FAa,0 = 0; μ = 2.5 × 10-5, v = μ/3, and U = μ; s = 0.01 to 10; and r = 0 to 1. Since the deterministic model assumes an infinite population size, we only examined a large population size of 107. For each condition, 100,000 simulations were repeated. We calculated the variance of the allele frequency at the linked neutral locus L at the corresponding t nearlyfixed . Under neutrality in the absence of mutation, the allele frequency variance can be calculated by p ( 1 p ) [ 1 ( 1 1 N ) t ] MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xH8viVGI8Gi=hEeeu0xXdbba9frFj0xb9qqpG0dXdb9aspeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGaemiCaaNaeiikaGIaeGymaeJaeyOeI0IaemiCaaNaeiykaKIaei4waSLaeGymaeJaeyOeI0IaeiikaGIaeGymaeJaeyOeI0scfa4aaSaaaeaacqaIXaqmaeaacqWGobGtaaGccqGGPaqkdaahaaWcbeqaaiabdsha0baakiabc2faDbaa@3E87@ [34]. Therefore, the population size under neutrality (N e ) that has the same variance in allele frequency as the population under selection can be determined using numerical iteration. In the presence of mutation, we used simulation to determine the range of the population size under neutrality. These were used to determine the range of allele frequency variances that matched the frequency variance under selection at the corresponding t nearlyfixed .


Recurrence equation for Fin the absence of selection

In the absence of mutation or selection, the inbreeding coefficient is

F t = 1 N + ( 1 1 N ) F t 1 = 1 ( 1 1 N ) t ( 1 F 0 ) MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xI8qiVKYPFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGaemOray0aaSbaaSqaaiabdsha0bqabaGccqGH9aqpjuaGdaWcaaqaaiabigdaXaqaaiabd6eaobaakiabgUcaRiabcIcaOiabigdaXKqbakabgkHiTmaalaaabaGaeGymaedabaGaemOta4eaaiabcMcaPiabdAeagnaaBaaabaGaemiDaqNaeyOeI0IaeGymaedabeaacqGH9aqpcqaIXaqmcqGHsislcqGGOaakcqaIXaqmcqGHsisldaWcaaqaaiabigdaXaqaaiabd6eaobaacqGGPaqkdaahaaqabeaacqWG0baDaaGaeiikaGIaeGymaeJaeyOeI0IaemOray0aaSbaaeaacqaIWaamaeqaaiabcMcaPaaa@4E27@

where t is time in generations and N is the population size [26]. 1 N MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xH8viVGI8Gi=hEeeu0xXdbba9frFj0xb9qqpG0dXdb9aspeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaqcfa4aaSaaaeaacqaIXaqmaeaacqWGobGtaaaaaa@2E88@ gives the probability that two offspring are derived from same parent in which case the probability of them being identical by descendent is 1. ( 1 1 N ) MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xH8viVGI8Gi=hEeeu0xXdbba9frFj0xb9qqpG0dXdb9aspeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGaeiikaGIaeGymaeJaeyOeI0scfa4aaSaaaeaacqaIXaqmaeaacqWGobGtaaGccqGGPaqkaaa@3221@ is the probability that two offspring are derived from different parents in which case the probability of them being identical by descendent is Ft-1. In the presence of mutation, F t = [ 1 N + ( 1 1 N ) F t 1 ] ( 1 U ) 2 MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xH8viVGI8Gi=hEeeu0xXdbba9frFj0xb9qqpG0dXdb9aspeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGaemOray0aaSbaaSqaaiabdsha0bqabaGccqGH9aqpcqGGBbWwjuaGdaWcaaqaaiabigdaXaqaaiabd6eaobaakiabgUcaRiabcIcaOiabigdaXiabgkHiTKqbaoaalaaabaGaeGymaedabaGaemOta4eaaOGaeiykaKIaemOray0aaSbaaSqaaiabdsha0jabgkHiTiabigdaXaqabaGccqGGDbqxcqGGOaakcqaIXaqmcqGHsislcqWGvbqvcqGGPaqkdaahaaWcbeqaaiabikdaYaaaaaa@467C@ [35]. To obtain F t in terms of F0, let α = 1 N × ( 1 U ) 2 MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xH8viVGI8Gi=hEeeu0xXdbba9frFj0xb9qqpG0dXdb9aspeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGaeqySdeMaeyypa0tcfa4aaSaaaeaacqaIXaqmaeaacqWGobGtaaGccqGHxdaTcqGGOaakcqaIXaqmcqGHsislcqWGvbqvcqGGPaqkdaahaaWcbeqaaiabikdaYaaaaaa@392F@ , and β = ( 1 1 N ) × ( 1 U ) 2 MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xH8viVGI8Gi=hEeeu0xXdbba9frFj0xb9qqpG0dXdb9aspeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGaeqOSdiMaeyypa0JaeiikaGIaeGymaeJaeyOeI0scfa4aaSaaaeaacqaIXaqmaeaacqWGobGtaaGaeiykaKIccqGHxdaTcqGGOaakcqaIXaqmcqGHsislcqWGvbqvcqGGPaqkdaahaaWcbeqaaiabikdaYaaaaaa@3CC0@ . This gives

F1 = α + βF0

F2 = α + βF1 = α + β × (α + βF0) = α + αβ + β2F0

F3 = α + βF2 = α + β × (α + αβ + β2F0) = α + αβ + αβ2 + β3F0

F t = α + αβ + αβ2 + αβ3 + ... + αβt-1 + βtF0.

The formula, 1 + x + ... + xn-1 = (1-xn)/(1-x), gives the following:

F t = α (1 - βt)/(1 - β) + βtF0, (2)

As t approaches infinity, F converges to the equilibrium F ^ = α ( 1 β ) 1 1 + 2 N U MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xH8viVGI8Gi=hEeeu0xXdbba9frFj0xb9qqpG0dXdb9aspeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGafmOrayKbaKaacqGH9aqpjuaGdaWcaaqaaiabeg7aHbqaaiabcIcaOiabigdaXiabgkHiTiabek7aIjabcMcaPaaakiabgIKi7MqbaoaalaaabaGaeGymaedabaGaeGymaeJaey4kaSIaeGOmaiJaemOta4Kaemyvaufaaaaa@3DD2@ , as shown previously by Kimura and Crow [35].

Recurrence equations for Fin the presence of a selected locus without recombination

In the presence of selection, the F value at time t is the sum of the probability of two offspring being derived from parents having alleles AA, aa, or Aa at locus S multiplied by the probability that the offspring will be identical by descent at locus L. In other words:

F t = { p A , t 2 [ 1 N A t 1 + ( 1 1 N A t 1 ) F A A , t 1 ] + p a , t 2 [ 1 N a t 1 + ( 1 1 N a t 1 ) F a a , t 1 ] + 2 p A , t p a , t F A a , t 1 } ( 1 U ) 2 . MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xI8qiVKYPFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGaemOray0aaSbaaSqaaiabdsha0bqabaGccqGH9aqpdaGadaqaaiabdchaWnaaDaaaleaacqWGbbqqcqGGSaalcqWG0baDaeaacqaIYaGmaaGccqGGBbWwjuaGdaWcaaqaaiabigdaXaqaaiabd6eaojabdgeabnaaBaaabaGaemiDaqNaeyOeI0IaeGymaedabeaaaaGccqGHRaWkcqGGOaakcqaIXaqmcqGHsisljuaGdaWcaaqaaiabigdaXaqaaiabd6eaojabdgeabnaaBaaabaGaemiDaqNaeyOeI0IaeGymaedabeaaaaGccqGGPaqkcqWGgbGrdaWgaaWcbaGaemyqaeKaemyqaeKaeiilaWIaemiDaqNaeyOeI0IaeGymaedabeaakiabc2faDjabgUcaRiabdchaWnaaDaaaleaacqWGHbqycqGGSaalcqWG0baDaeaacqaIYaGmaaGccqGGBbWwjuaGdaWcaaqaaiabigdaXaqaaiabd6eaojabdggaHnaaBaaabaGaemiDaqNaeyOeI0IaeGymaedabeaaaaGccqGHRaWkcqGGOaakcqaIXaqmcqGHsisljuaGdaWcaaqaaiabigdaXaqaaiabd6eaojabdggaHnaaBaaabaGaemiDaqNaeyOeI0IaeGymaedabeaaaaGccqGGPaqkcqWGgbGrdaWgaaWcbaGaemyyaeMaemyyaeMaeiilaWIaemiDaqNaeyOeI0IaeGymaedabeaakiabc2faDjabgUcaRiabikdaYiabdchaWnaaBaaaleaacqWGbbqqcqGGSaalcqWG0baDaeqaaOGaemiCaa3aaSbaaSqaaiabdggaHjabcYcaSiabdsha0bqabaGccqWGgbGrdaWgaaWcbaGaemyqaeKaemyyaeMaeiilaWIaemiDaqNaeyOeI0IaeGymaedabeaaaOGaay5Eaiaaw2haaiabcIcaOiabigdaXiabgkHiTiabdwfavjabcMcaPmaaCaaaleqabaGaeGOmaidaaOGaeiOla4caaa@9498@

Here, At-1 and at-1 are the frequencies of the advantageous and disadvantageous alleles at locus S at generation t-1. p A , t = w A t 1 w A t 1 + a t 1 MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xH8viVGI8Gi=hEeeu0xXdbba9frFj0xb9qqpG0dXdb9aspeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGaemiCaa3aaSbaaSqaaiabdgeabjabcYcaSiabdsha0bqabaGccqGH9aqpjuaGdaWcaaqaaiabdEha3jabdgeabnaaBaaabaGaemiDaqNaeyOeI0IaeGymaedabeaaaeaacqWG3bWDcqWGbbqqdaWgaaqaaiabdsha0jabgkHiTiabigdaXaqabaGaey4kaSIaemyyae2aaSbaaeaacqWG0baDcqGHsislcqaIXaqmaeqaaaaaaaa@43F2@ and p a , t = a t 1 w A t 1 + a t 1 MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xH8viVGI8Gi=hEeeu0xXdbba9frFj0xb9qqpG0dXdb9aspeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGaemiCaa3aaSbaaSqaaiabdggaHjabcYcaSiabdsha0bqabaGccqGH9aqpjuaGdaWcaaqaaiabdggaHnaaBaaabaGaemiDaqNaeyOeI0IaeGymaedabeaaaeaacqWG3bWDcqWGbbqqdaWgaaqaaiabdsha0jabgkHiTiabigdaXaqabaGaey4kaSIaemyyae2aaSbaaeaacqWG0baDcqGHsislcqaIXaqmaeqaaaaaaaa@42FB@ give the probabilities that an offspring at generation t is derived from a parent at generation t-1 with allele A or a, respectively. F AA , F Aa , and F aa give the probabilities that parents with the indicated alleles will be identical by descent at locus L. Given that both parents have allele A or a at locus S, the 1 N A t 1 MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xH8viVGI8Gi=hEeeu0xXdbba9frFj0xb9qqpG0dXdb9aspeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaqcfa4aaSaaaeaacqaIXaqmaeaacqWGobGtcqWGbbqqdaWgaaqaaiabdsha0jabgkHiTiabigdaXaqabaaaaaaa@3302@ and 1 N a t 1 MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xH8viVGI8Gi=hEeeu0xXdbba9frFj0xb9qqpG0dXdb9aspeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaqcfa4aaSaaaeaacqaIXaqmaeaacqWGobGtcqWGHbqydaWgaaqaaiabdsha0jabgkHiTiabigdaXaqabaaaaaaa@3342@ terms respectively give the probabilities that two offspring have the same parent (in which case the probability of being identical by descent at locus L, in the absence of mutation is 1). The 1 1 N A t 1 MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xH8viVGI8Gi=hEeeu0xXdbba9frFj0xb9qqpG0dXdb9aspeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaqcfaOaeGymaeJaeyOeI0YaaSaaaeaacqaIXaqmaeaacqWGobGtcqWGbbqqdaWgaaqaaiabdsha0jabgkHiTiabigdaXaqabaaaaaaa@34DF@ and 1 1 N a t 1 MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xH8viVGI8Gi=hEeeu0xXdbba9frFj0xb9qqpG0dXdb9aspeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaqcfaOaeGymaeJaeyOeI0YaaSaaaeaacqaIXaqmaeaacqWGobGtcqWGHbqydaWgaaqaaiabdsha0jabgkHiTiabigdaXaqabaaaaaaa@351F@ terms give the probability that the two offspring came from different parents (in which case the probabilities of identity by descent at locus L, in the absence of mutation, are FAA, t-1 and Faa, t-1 respectively). The term (1-U)2 accounts for the fact that two individuals cannot be identical by descent if there is a mutation at the neutral locus L.

If the parameters w, μ, v, U, A0, a0, F0, FAA,0, Faa,0, and FAa,0 are known, F1 can be calculated using equation (3). In addition, FAA,1, Faa,1, and FAa,1 can be calculated using the following equations:

F A A , t = { ( p A A , t ) 2 [ 1 N A t 1 + ( 1 1 N A t 1 ) F A A , t 1 ] + ( p a A , t ) 2 [ 1 N a t 1 + ( 1 1 N a t 1 ) F a a , t 1 ] + 2 p A A , t p a A , t F A a , t 1 } ( 1 U ) 2 , MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xI8qiVKYPFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaqbaeaabiqaaaqaaiabdAeagnaaBaaaleaacqWGbbqqcqWGbbqqcqGGSaalcqWG0baDaeqaaOGaeyypa0Jaei4EaSNaeiikaGIaemiCaa3aaSbaaSqaaiabdgeabjabgkziUkabdgeabjabcYcaSiabdsha0bqabaGccqGGPaqkdaahaaWcbeqaaiabikdaYaaakiabcUfaBLqbaoaalaaabaGaeGymaedabaGaemOta4Kaemyqae0aaSbaaeaacqWG0baDcqGHsislcqaIXaqmaeqaaaaakiabgUcaRiabcIcaOiabigdaXiabgkHiTKqbaoaalaaabaGaeGymaedabaGaemOta4Kaemyqae0aaSbaaeaacqWG0baDcqGHsislcqaIXaqmaeqaaaaakiabcMcaPiabdAeagnaaBaaaleaacqWGbbqqcqWGbbqqcqGGSaalcqWG0baDcqGHsislcqaIXaqmaeqaaOGaeiyxa0Laey4kaSIaeiikaGIaemiCaa3aaSbaaSqaaiabdggaHjabgkziUkabdgeabjabcYcaSiabdsha0bqabaGccqGGPaqkdaahaaWcbeqaaiabikdaYaaakiabcUfaBLqbaoaalaaabaGaeGymaedabaGaemOta4Kaemyyae2aaSbaaeaacqWG0baDcqGHsislcqaIXaqmaeqaaaaakiabgUcaRiabcIcaOiabigdaXiabgkHiTKqbaoaalaaabaGaeGymaedabaGaemOta4Kaemyyae2aaSbaaeaacqWG0baDcqGHsislcqaIXaqmaeqaaaaakiabcMcaPiabdAeagnaaBaaaleaacqWGHbqycqWGHbqycqGGSaalcqWG0baDcqGHsislcqaIXaqmaeqaaOGaeiyxa0fabaGaey4kaSIaeGOmaiJaemiCaa3aaSbaaSqaaiabdgeabjabgkziUkabdgeabjabcYcaSiabdsha0bqabaGccqWGWbaCdaWgaaWcbaGaemyyaeMaeyOKH4QaemyqaeKaeiilaWIaemiDaqhabeaakiabdAeagnaaBaaaleaacqWGbbqqcqWGHbqycqGGSaalcqWG0baDcqGHsislcqaIXaqmaeqaaOGaeiyFa0NaeiikaGIaeGymaeJaeyOeI0IaemyvauLaeiykaKYaaWbaaSqabeaacqaIYaGmaaGccqGGSaalaaaaaa@A815@
F a a , t = { ( p a a , t ) 2 [ 1 N a t 1 + ( 1 1 N a t 1 ) F a a , t 1 ] + ( p A a , t ) 2 [ 1 N A t 1 + ( 1 1 N A t 1 ) F A A , t 1 ] + 2 p A a , t p a a , t F A a , t 1 } ( 1 U ) 2 MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xI8qiVKYPFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaqbaeaabiqaaaqaaiabdAeagnaaBaaaleaacqWGHbqycqWGHbqycqGGSaalcqWG0baDaeqaaOGaeyypa0Jaei4EaSNaeiikaGIaemiCaa3aaSbaaSqaaiabdggaHjabgkziUkabdggaHjabcYcaSiabdsha0bqabaGccqGGPaqkdaahaaWcbeqaaiabikdaYaaakiabcUfaBLqbaoaalaaabaGaeGymaedabaGaemOta4Kaemyyae2aaSbaaeaacqWG0baDcqGHsislcqaIXaqmaeqaaaaakiabgUcaRiabcIcaOiabigdaXiabgkHiTKqbaoaalaaabaGaeGymaedabaGaemOta4Kaemyyae2aaSbaaeaacqWG0baDcqGHsislcqaIXaqmaeqaaaaakiabcMcaPiabdAeagnaaBaaaleaacqWGHbqycqWGHbqycqGGSaalcqWG0baDcqGHsislcqaIXaqmaeqaaOGaeiyxa0Laey4kaSIaeiikaGIaemiCaa3aaSbaaSqaaiabdgeabjabgkziUkabdggaHjabcYcaSiabdsha0bqabaGccqGGPaqkdaahaaWcbeqaaiabikdaYaaakiabcUfaBLqbaoaalaaabaGaeGymaedabaGaemOta4Kaemyqae0aaSbaaeaacqWG0baDcqGHsislcqaIXaqmaeqaaaaakiabgUcaRiabcIcaOiabigdaXiabgkHiTKqbaoaalaaabaGaeGymaedabaGaemOta4Kaemyqae0aaSbaaeaacqWG0baDcqGHsislcqaIXaqmaeqaaaaakiabcMcaPiabdAeagnaaBaaaleaacqWGbbqqcqWGbbqqcqGGSaalcqWG0baDcqGHsislcqaIXaqmaeqaaOGaeiyxa0fabaGaey4kaSIaeGOmaiJaemiCaa3aaSbaaSqaaiabdgeabjabgkziUkabdggaHjabcYcaSiabdsha0bqabaGccqWGWbaCdaWgaaWcbaGaemyyaeMaeyOKH4QaemyyaeMaeiilaWIaemiDaqhabeaakiabdAeagnaaBaaaleaacqWGbbqqcqWGHbqycqGGSaalcqWG0baDcqGHsislcqaIXaqmaeqaaOGaeiyFa0NaeiikaGIaeGymaeJaeyOeI0IaemyvauLaeiykaKYaaWbaaSqabeaacqaIYaGmaaaaaaaa@A8AB@


FAa,t= (pAA,tpaa,tFAa,t-1 + paa,tpaA,tFaa,t-1 + pAa,tpAA,tFAA,t-1 + pAa,tpaA,tFAa,t-1)(1 - U)2

Where pxy, tis the probability that a sampled offspring is descended from a parent with allele x given that the offspring has allele y at locus S. The reasoning behind equations (4) – (6) is similar to that for equation (3). For each offspring, pxy, tcan be calculated as the probability of the parent having allele x at locus S multiplied by the probability that x mutates to y (or fails to mutate, if x = y), divided by the probability that the offspring is y. In other words,

p A A , t = w × A t 1 w × A t 1 + a t 1 ( 1 μ ) / A t , p a A , t = a t 1 w × A t 1 + a t 1 v / A t MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xI8qiVKYPFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaqbaeqabeGaaaqaaiabdchaWnaaBaaaleaacqWGbbqqcqGHsgIRcqWGbbqqcqGGSaalcqWG0baDaeqaaOGaeyypa0tcfa4aaSaaaeaacqWG3bWDcqGHxdaTcqWGbbqqdaWgaaqaaiabdsha0jabgkHiTiabigdaXaqabaaabaGaem4DaCNaey41aqRaemyqae0aaSbaaeaacqWG0baDcqGHsislcqaIXaqmaeqaaiabgUcaRiabdggaHnaaBaaabaGaemiDaqNaeyOeI0IaeGymaedabeaaaaGccqGGOaakcqaIXaqmcqGHsislcqaH8oqBcqGGPaqkcqGGVaWlcqWGbbqqdaWgaaWcbaGaemiDaqhabeaakiabcYcaSaqaaiabdchaWnaaBaaaleaacqWGHbqycqGHsgIRcqWGbbqqcqGGSaalcqWG0baDaeqaaOGaeyypa0tcfa4aaSaaaeaacqWGHbqydaWgaaqaaiabdsha0jabgkHiTiabigdaXaqabaaabaGaem4DaCNaey41aqRaemyqae0aaSbaaeaacqWG0baDcqGHsislcqaIXaqmaeqaaiabgUcaRiabdggaHnaaBaaabaGaemiDaqNaeyOeI0IaeGymaedabeaaaaGaemODayNccqGGVaWlcqWGbbqqdaWgaaWcbaGaemiDaqhabeaaaaaaaa@767C@


p a a , t = a t 1 w × A t 1 + a t 1 ( 1 v ) / a t , p A a , t = w × A t 1 w × A t 1 + a t 1 μ / a t MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xI8qiVKYPFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaqbaeqabeGaaaqaaiabdchaWnaaBaaaleaacqWGHbqycqGHsgIRcqWGHbqycqGGSaalcqWG0baDaeqaaOGaeyypa0tcfa4aaSaaaeaacqWGHbqydaWgaaqaaiabdsha0jabgkHiTiabigdaXaqabaaabaGaem4DaCNaey41aqRaemyqae0aaSbaaeaacqWG0baDcqGHsislcqaIXaqmaeqaaiabgUcaRiabdggaHnaaBaaabaGaemiDaqNaeyOeI0IaeGymaedabeaaaaGccqGGOaakcqaIXaqmcqGHsislcqWG2bGDcqGGPaqkcqGGVaWlcqWGHbqydaWgaaWcbaGaemiDaqhabeaakiabcYcaSaqaaiabdchaWnaaBaaaleaacqWGbbqqcqGHsgIRcqWGHbqycqGGSaalcqWG0baDaeqaaOGaeyypa0tcfa4aaSaaaeaacqWG3bWDcqGHxdaTcqWGbbqqdaWgaaqaaiabdsha0jabgkHiTiabigdaXaqabaaabaGaem4DaCNaey41aqRaemyqae0aaSbaaeaacqWG0baDcqGHsislcqaIXaqmaeqaaiabgUcaRiabdggaHnaaBaaabaGaemiDaqNaeyOeI0IaeGymaedabeaaaaGaeqiVd0MccqGGVaWlcqWGHbqydaWgaaWcbaGaemiDaqhabeaaaaaaaa@777C@

where A t and a t are given by

A t = w × A t 1 w × A t 1 + a t 1 ( 1 μ ) + a t 1 w × A t 1 + a t 1 v MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xI8qiVKYPFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGaemyqae0aaSbaaSqaaiabdsha0bqabaGccqGH9aqpjuaGdaWcaaqaaiabdEha3jabgEna0kabdgeabnaaBaaabaGaemiDaqNaeyOeI0IaeGymaedabeaaaeaacqWG3bWDcqGHxdaTcqWGbbqqdaWgaaqaaiabdsha0jabgkHiTiabigdaXaqabaGaey4kaSIaemyyae2aaSbaaeaacqWG0baDcqGHsislcqaIXaqmaeqaaaaakiabcIcaOiabigdaXiabgkHiTiabeY7aTjabcMcaPiabgUcaRKqbaoaalaaabaGaemyyae2aaSbaaeaacqWG0baDcqGHsislcqaIXaqmaeqaaaqaaiabdEha3jabgEna0kabdgeabnaaBaaabaGaemiDaqNaeyOeI0IaeGymaedabeaacqGHRaWkcqWGHbqydaWgaaqaaiabdsha0jabgkHiTiabigdaXaqabaaaaOGaemODayhaaa@60D1@
a t = a t 1 w × A t 1 + a t 1 ( 1 v ) + w A t 1 w × A t 1 + a t 1 μ MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xI8qiVKYPFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGaemyyae2aaSbaaSqaaiabdsha0bqabaGccqGH9aqpjuaGdaWcaaqaaiabdggaHnaaBaaabaGaemiDaqNaeyOeI0IaeGymaedabeaaaeaacqWG3bWDcqGHxdaTcqWGbbqqdaWgaaqaaiabdsha0jabgkHiTiabigdaXaqabaGaey4kaSIaemyyae2aaSbaaeaacqWG0baDcqGHsislcqaIXaqmaeqaaaaakiabcIcaOiabigdaXiabgkHiTiabdAha2jabcMcaPiabgUcaRKqbaoaalaaabaGaem4DaCNaemyqae0aaSbaaeaacqWG0baDcqGHsislcqaIXaqmaeqaaaqaaiabdEha3jabgEna0kabdgeabnaaBaaabaGaemiDaqNaeyOeI0IaeGymaedabeaacqGHRaWkcqWGHbqydaWgaaqaaiabdsha0jabgkHiTiabigdaXaqabaaaaOGaeqiVd0gaaa@5EFA@

When A1, a1, F1, FAA,1, Faa,1, and FAa,1 are available, we can calculate F2 using equation (3), and FAA,2, Faa,2, and FAa,2 using equations (4) to (6). Therefore, F t can be obtained by iteration.

Recurrence equations for Fin the presence of a selected locus and recombination

Assuming that loci L and S recombine with a probability r per generation, we obtain

P A A , t = p A , t ( 1 r ) + p A , t r N p A , t 1 N 1 MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xI8qiVKYPFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGaemiuaa1aaSbaaSqaaiabdgeabjabdgeabjabcYcaSiabdsha0bqabaGccqGH9aqpcqWGWbaCdaWgaaWcbaGaemyqaeKaeiilaWIaemiDaqhabeaakiabcIcaOiabigdaXiabgkHiTiabdkhaYjabcMcaPiabgUcaRiabdchaWnaaBaaaleaacqWGbbqqcqGGSaalcqWG0baDaeqaaOGaemOCaixcfa4aaSaaaeaacqWGobGtcqWGWbaCdaWgaaqaaiabdgeabjabcYcaSiabdsha0bqabaGaeyOeI0IaeGymaedabaGaemOta4KaeyOeI0IaeGymaedaaaaa@4FB8@
P a a , t = p a , t ( 1 r ) + p a , t r N p a , t 1 N 1 MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xI8qiVKYPFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGaemiuaa1aaSbaaSqaaiabdggaHjabdggaHjabcYcaSiabdsha0bqabaGccqGH9aqpcqWGWbaCdaWgaaWcbaGaemyyaeMaeiilaWIaemiDaqhabeaakiabcIcaOiabigdaXiabgkHiTiabdkhaYjabcMcaPiabgUcaRiabdchaWnaaBaaaleaacqWGHbqycqGGSaalcqWG0baDaeqaaOGaemOCaixcfa4aaSaaaeaacqWGobGtcqWGWbaCdaWgaaqaaiabdggaHjabcYcaSiabdsha0bqabaGaeyOeI0IaeGymaedabaGaemOta4KaeyOeI0IaeGymaedaaaaa@50F8@
P A a , t = p a A , t = p A , t r N p a , t N 1 MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xI8qiVKYPFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGaemiuaa1aaSbaaSqaaiabdgeabjabdggaHjabcYcaSiabdsha0bqabaGccqGH9aqpcqWGWbaCdaWgaaWcbaGaemyyaeMaemyqaeKaeiilaWIaemiDaqhabeaakiabg2da9iabdchaWnaaBaaaleaacqWGbbqqcqGGSaalcqWG0baDaeqaaOGaemOCaixcfa4aaSaaaeaacqWGobGtcqWGWbaCdaWgaaqaaiabdggaHjabcYcaSiabdsha0bqabaaabaGaemOta4KaeyOeI0IaeGymaedaaaaa@4ACE@

where pxy, tis the probability that an individual at generation t has a neutral locus L derived from an individual with allele x at locus S, and a selected locus S derived from an individual with allele y at locus S (x and y can be A or a). This probability is the sum of the probability of no recombination and the probability of recombination between individuals with indicated allele at locus S. The -1's in the (Npx, t- 1) and (N - 1) terms above account for the fact that a haploid individual cannot recombine with itself.

Similar to equations (3) – (6), with recombination,

F t = { ( p A A , t 2 + p A a , t 2 + 2 p A A , t p A a , t ) [ 1 N A t 1 + ( 1 1 N A t 1 ) F A A , t 1 ] + ( p a A , t 2 + p a a , t 2 + 2 p a A , t p a a , t ) [ 1 N a t 1 + ( 1 1 N a t 1 ) F a a , t 1 ] + 2 ( p A A , t p a A , t + p A A , t p a a , t + p a A , t p A a , t + p a a , t p A a , t ) F A a , t 1 } ( 1 U ) 2 = { p A , t 2 [ 1 N A t 1 + ( 1 1 N A t 1 ) F A A , t 1 ] + p a , t 2 [ 1 N a t 1 + ( 1 1 N a t 1 ) F a a , t 1 ] + 2 p A , t p a , t F A a , t 1 } ( 1 U ) 2 MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xI8qiVKYPFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaqbaeaabiqaaaqaaiabdAeagnaaBaaaleaacqWG0baDaeqaaOGaeyypa0ZaaiWaaeaafaqaaeWabaaabaGaeiikaGIaemiCaa3aa0baaSqaaiabdgeabjabdgeabjabcYcaSiabdsha0bqaaiabikdaYaaakiabgUcaRiabdchaWnaaDaaaleaacqWGbbqqcqWGHbqycqGGSaalcqWG0baDaeaacqaIYaGmaaGccqGHRaWkcqaIYaGmcqWGWbaCdaWgaaWcbaGaemyqaeKaemyqaeKaeiilaWIaemiDaqhabeaakiabdchaWnaaBaaaleaacqWGbbqqcqWGHbqycqGGSaalcqWG0baDaeqaaOGaeiykaKIaei4waSvcfa4aaSaaaeaacqaIXaqmaeaacqWGobGtcqWGbbqqdaWgaaqaaiabdsha0jabgkHiTiabigdaXaqabaaaaOGaey4kaSIaeiikaGIaeGymaeJaeyOeI0scfa4aaSaaaeaacqaIXaqmaeaacqWGobGtcqWGbbqqdaWgaaqaaiabdsha0jabgkHiTiabigdaXaqabaaaaOGaeiykaKIaemOray0aaSbaaSqaaiabdgeabjabdgeabjabcYcaSiabdsha0jabgkHiTiabigdaXaqabaGccqGGDbqxaeaacqGHRaWkcqGGOaakcqWGWbaCdaqhaaWcbaGaemyyaeMaemyqaeKaeiilaWIaemiDaqhabaGaeGOmaidaaOGaey4kaSIaemiCaa3aa0baaSqaaiabdggaHjabdggaHjabcYcaSiabdsha0bqaaiabikdaYaaakiabgUcaRiabikdaYiabdchaWnaaBaaaleaacqWGHbqycqWGbbqqcqGGSaalcqWG0baDaeqaaOGaemiCaa3aaSbaaSqaaiabdggaHjabdggaHjabcYcaSiabdsha0bqabaGccqGGPaqkcqGGBbWwjuaGdaWcaaqaaiabigdaXaqaaiabd6eaojabdggaHnaaBaaabaGaemiDaqNaeyOeI0IaeGymaedabeaaaaGccqGHRaWkcqGGOaakcqaIXaqmcqGHsisljuaGdaWcaaqaaiabigdaXaqaaiabd6eaojabdggaHnaaBaaabaGaemiDaqNaeyOeI0IaeGymaedabeaaaaGccqGGPaqkcqWGgbGrdaWgaaWcbaGaemyyaeMaemyyaeMaeiilaWIaemiDaqNaeyOeI0IaeGymaedabeaakiabc2faDbqaaiabgUcaRiabikdaYiabcIcaOiabdchaWnaaBaaaleaacqWGbbqqcqWGbbqqcqGGSaalcqWG0baDaeqaaOGaemiCaa3aaSbaaSqaaiabdggaHjabdgeabjabcYcaSiabdsha0bqabaGccqGHRaWkcqWGWbaCdaWgaaWcbaGaemyqaeKaemyqaeKaeiilaWIaemiDaqhabeaakiabdchaWnaaBaaaleaacqWGHbqycqWGHbqycqGGSaalcqWG0baDaeqaaOGaey4kaSIaemiCaa3aaSbaaSqaaiabdggaHjabdgeabjabcYcaSiabdsha0bqabaGccqWGWbaCdaWgaaWcbaGaemyqaeKaemyyaeMaeiilaWIaemiDaqhabeaakiabgUcaRiabdchaWnaaBaaaleaacqWGHbqycqWGHbqycqGGSaalcqWG0baDaeqaaOGaemiCaa3aaSbaaSqaaiabdgeabjabdggaHjabcYcaSiabdsha0bqabaGccqGGPaqkcqWGgbGrdaWgaaWcbaGaemyqaeKaemyyaeMaeiilaWIaemiDaqNaeyOeI0IaeGymaedabeaaaaaakiaawUhacaGL9baacqGGOaakcqaIXaqmcqGHsislcqWGvbqvcqGGPaqkdaahaaWcbeqaaiabikdaYaaaaOqaaiabg2da9maacmaabaGaemiCaa3aa0baaSqaaiabdgeabjabcYcaSiabdsha0bqaaiabikdaYaaakiabcUfaBLqbaoaalaaabaGaeGymaedabaGaemOta4Kaemyqae0aaSbaaeaacqWG0baDcqGHsislcqaIXaqmaeqaaaaakiabgUcaRiabcIcaOiabigdaXiabgkHiTKqbaoaalaaabaGaeGymaedabaGaemOta4Kaemyqae0aaSbaaeaacqWG0baDcqGHsislcqaIXaqmaeqaaaaakiabcMcaPiabdAeagnaaBaaaleaacqWGbbqqcqWGbbqqcqGGSaalcqWG0baDcqGHsislcqaIXaqmaeqaaOGaeiyxa0Laey4kaSIaemiCaa3aa0baaSqaaiabdggaHjabcYcaSiabdsha0bqaaiabikdaYaaakiabcUfaBLqbaoaalaaabaGaeGymaedabaGaemOta4Kaemyyae2aaSbaaeaacqWG0baDcqGHsislcqaIXaqmaeqaaaaakiabgUcaRiabcIcaOiabigdaXiabgkHiTKqbaoaalaaabaGaeGymaedabaGaemOta4Kaemyyae2aaSbaaeaacqWG0baDcqGHsislcqaIXaqmaeqaaaaakiabcMcaPiabdAeagnaaBaaaleaacqWGHbqycqWGHbqycqGGSaalcqWG0baDcqGHsislcqaIXaqmaeqaaOGaeiyxa0Laey4kaSIaeGOmaiJaemiCaa3aaSbaaSqaaiabdgeabjabcYcaSiabdsha0bqabaGccqWGWbaCdaWgaaWcbaGaemyyaeMaeiilaWIaemiDaqhabeaakiabdAeagnaaBaaaleaacqWGbbqqcqWGHbqycqGGSaalcqWG0baDcqGHsislcqaIXaqmaeqaaaGccaGL7bGaayzFaaGaeiikaGIaeGymaeJaeyOeI0IaemyvauLaeiykaKYaaWbaaSqabeaacqaIYaGmaaaaaaaa@579B@
F A A , t = { ( p A A ( 1 μ ) + p A a v A t ) 2 [ 1 N A t 1 + ( 1 1 N A t 1 ) F A A t 1 ] +( p a A ( 1 μ ) + p a a v A t ) 2 [ 1 N a t 1 + ( 1 1 N a t 1 ) F a a t 1 ] + 2[ p A A ( 1 μ ) × p a A ( 1 μ ) + p A A ( 1 μ ) × p a a v + p a A ( 1 μ ) × p A a v + p A a v × p a a v ] A t 2 F A a , t 1 } ( 1 U ) 2 MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xI8qiVKYPFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaqbaeaabmqaaaqaaiabdAeagnaaBaaaleaacqWGbbqqcqWGbbqqcqGGSaalcqWG0baDaeqaaOGaeyypa0Jaei4EaSNaeiikaGscfa4aaSaaaeaacqWGWbaCdaWgaaqaaiabdgeabjabdgeabbqabaGaeiikaGIaeGymaeJaeyOeI0IaeqiVd0MaeiykaKIaey4kaSIaemiCaa3aaSbaaeaacqWGbbqqcqWGHbqyaeqaaiabdAha2bqaaiabdgeabnaaBaaabaGaemiDaqhabeaaaaGccqGGPaqkdaahaaWcbeqaaiabikdaYaaakiabcUfaBLqbaoaalaaabaGaeGymaedabaGaemOta4Kaemyqae0aaSbaaeaacqWG0baDcqGHsislcqaIXaqmaeqaaaaakiabgUcaRiabbIcaOiabigdaXiabgkHiTKqbaoaalaaabaGaeGymaedabaGaemOta4Kaemyqae0aaSbaaeaacqWG0baDcqGHsislcqaIXaqmaeqaaaaakiabbMcaPiabdAeagnaaBaaaleaacqWGbbqqcqWGbbqqcqqGSaalcqqGGaaicqWG0baDcqGHsislcqaIXaqmaeqaaOGaeeyxa0fabaGaee4kaSIaeeikaGscfa4aaSaaaeaacqWGWbaCdaWgaaqaaiabdggaHjabdgeabbqabaGaeiikaGIaeGymaeJaeyOeI0IaeqiVd0MaeiykaKIaey4kaSIaemiCaa3aaSbaaeaacqWGHbqycqWGHbqyaeqaaiabdAha2bqaaiabdgeabnaaBaaabaGaemiDaqhabeaaaaGccqqGPaqkdaahaaWcbeqaaiabbkdaYaaakiabbUfaBLqbaoaalaaabaGaeGymaedabaGaemOta4Kaemyyae2aaSbaaeaacqWG0baDcqGHsislcqaIXaqmaeqaaaaakiabgUcaRiabbIcaOiabigdaXiabgkHiTKqbaoaalaaabaGaeGymaedabaGaemOta4Kaemyyae2aaSbaaeaacqWG0baDcqGHsislcqaIXaqmaeqaaaaakiabbMcaPiabdAeagnaaBaaaleaacqWGHbqycqWGHbqycqqGSaalcqqGGaaicqWG0baDcqGHsislcqaIXaqmaeqaaOGaeeyxa0fabaGaee4kaSscfa4aaSaaaeaacqqGYaGmcqqGBbWwcqWGWbaCdaWgaaqaaiabdgeabjabdgeabbqabaGaeiikaGIaeGymaeJaeyOeI0IaeqiVd0MaeiykaKIaey41aqRaemiCaa3aaSbaaeaacqWGHbqycqWGbbqqaeqaaiabcIcaOiabigdaXiabgkHiTiabeY7aTjabcMcaPiabgUcaRiabdchaWnaaBaaabaGaemyqaeKaemyqaeeabeaacqGGOaakcqaIXaqmcqGHsislcqaH8oqBcqGGPaqkcqGHxdaTcqWGWbaCdaWgaaqaaiabdggaHjabdggaHbqabaGaemODayNaey4kaSIaemiCaa3aaSbaaeaacqWGHbqycqWGbbqqaeqaaiabcIcaOiabigdaXiabgkHiTiabeY7aTjabcMcaPiabgEna0kabdchaWnaaBaaabaGaemyqaeKaemyyaegabeaacqWG2bGDcqGHRaWkcqWGWbaCdaWgaaqaaiabdgeabjabdggaHbqabaGaemODayNaey41aqRaemiCaa3aaSbaaeaacqWGHbqycqWGHbqyaeqaaiabdAha2jabc2faDbqaaiabdgeabnaaDaaabaGaemiDaqhabaGaeGOmaidaaaaakiabdAeagnaaBaaaleaacqWGbbqqcqWGHbqycqGGSaalcqWG0baDcqGHsislcqaIXaqmaeqaaOGaeiyFa0NaaGzaVlabcIcaOiabigdaXiabgkHiTiabdwfavjabcMcaPmaaCaaaleqabaGaeGOmaidaaaaaaaa@FAB9@
F a a , t = { ( p A a ( 1 v ) + p A A μ a t ) 2 [ 1 N A t 1 + ( 1 1 N A t 1 ) F A A t 1 ] +( p a a ( 1 v ) + p a A μ a t ) 2 [ 1 N a t 1 + ( 1 1 N a t 1 ) F a a t 1 ] + 2[ p A a ( 1 v ) × p a a ( 1 v ) + p A a ( 1 v ) × p a A μ + p a a ( 1 v ) × p A A μ + p A A μ × p a A μ ] a t 2 F A a , t 1 } ( 1 U ) 2 MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xI8qiVKYPFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaqbaeaabmqaaaqaaiabdAeagnaaBaaaleaacqWGHbqycqWGHbqycqGGSaalcqWG0baDaeqaaOGaeyypa0Jaei4EaSNaeiikaGscfa4aaSaaaeaacqWGWbaCdaWgaaqaaiabdgeabjabdggaHbqabaGaeiikaGIaeGymaeJaeyOeI0IaemODayNaeiykaKIaey4kaSIaemiCaa3aaSbaaeaacqWGbbqqcqWGbbqqaeqaaiabeY7aTbqaaiabdggaHnaaBaaabaGaemiDaqhabeaaaaGccqGGPaqkdaahaaWcbeqaaiabikdaYaaakiabcUfaBLqbaoaalaaabaGaeGymaedabaGaemOta4Kaemyqae0aaSbaaeaacqWG0baDcqGHsislcqaIXaqmaeqaaaaakiabgUcaRiabbIcaOiabigdaXiabgkHiTKqbaoaalaaabaGaeGymaedabaGaemOta4Kaemyqae0aaSbaaeaacqWG0baDcqGHsislcqaIXaqmaeqaaaaakiabbMcaPiabdAeagnaaBaaaleaacqWGbbqqcqWGbbqqcqqGSaalcqqGGaaicqWG0baDcqGHsislcqaIXaqmaeqaaOGaeeyxa0fabaGaee4kaSIaeeikaGscfa4aaSaaaeaacqWGWbaCdaWgaaqaaiabdggaHjabdggaHbqabaGaeiikaGIaeGymaeJaeyOeI0IaemODayNaeiykaKIaey4kaSIaemiCaa3aaSbaaeaacqWGHbqycqWGbbqqaeqaaiabeY7aTbqaaiabdggaHnaaBaaabaGaemiDaqhabeaaaaGccqqGPaqkdaahaaWcbeqaaiabbkdaYaaakiabbUfaBLqbaoaalaaabaGaeGymaedabaGaemOta4Kaemyyae2aaSbaaeaacqWG0baDcqGHsislcqaIXaqmaeqaaaaakiabgUcaRiabbIcaOiabigdaXiabgkHiTKqbaoaalaaabaGaeGymaedabaGaemOta4Kaemyyae2aaSbaaeaacqWG0baDcqGHsislcqaIXaqmaeqaaaaakiabbMcaPiabdAeagnaaBaaaleaacqWGHbqycqWGHbqycqqGSaalcqqGGaaicqWG0baDcqGHsislcqaIXaqmaeqaaOGaeeyxa0fabaGaee4kaSscfa4aaSaaaeaacqqGYaGmcqqGBbWwcqWGWbaCdaWgaaqaaiabdgeabjabdggaHbqabaGaeiikaGIaeGymaeJaeyOeI0IaemODayNaeiykaKIaey41aqRaemiCaa3aaSbaaeaacqWGHbqycqWGHbqyaeqaaiabcIcaOiabigdaXiabgkHiTiabdAha2jabcMcaPiabgUcaRiabdchaWnaaBaaabaGaemyqaeKaemyyaegabeaacqGGOaakcqaIXaqmcqGHsislcqWG2bGDcqGGPaqkcqGHxdaTcqWGWbaCdaWgaaqaaiabdggaHjabdgeabbqabaGaeqiVd0Maey4kaSIaemiCaa3aaSbaaeaacqWGHbqycqWGHbqyaeqaaiabcIcaOiabigdaXiabgkHiTiabdAha2jabcMcaPiabgEna0kabdchaWnaaBaaabaGaemyqaeKaemyqaeeabeaacqaH8oqBcqGHRaWkcqWGWbaCdaWgaaqaaiabdgeabjabdgeabbqabaGaeqiVd0Maey41aqRaemiCaa3aaSbaaeaacqWGHbqycqWGbbqqaeqaaiabeY7aTjabc2faDbqaaiabdggaHnaaDaaabaGaemiDaqhabaGaeGOmaidaaaaakiabdAeagnaaBaaaleaacqWGbbqqcqWGHbqycqGGSaalcqWG0baDcqGHsislcqaIXaqmaeqaaOGaeiyFa0NaaGzaVlabcIcaOiabigdaXiabgkHiTiabdwfavjabcMcaPmaaCaaaleqabaGaeGOmaidaaaaaaaa@FBF9@
F A a , t = { ( p A A ( 1 μ ) + p A a v A t ) × ( p A a ( 1 v ) + p A A μ a t ) [ 1 N A t 1 + ( 1 1 N A t 1 ) F A A t 1 ] +( p a A ( 1 μ ) + p a a v A t ) × ( p a a ( 1 v ) + p a A μ a t )[ 1 N a t 1 + ( 1 1 N a t 1 ) F a a t 1 ] + 2 p A A p a A μ ( 1 μ ) + 2 p A a p a a v ( 1 v ) + ( p A A p a a + p a A p A a ) [ ( 1 μ ) ( 1 v ) + μ v ] A t a t F A a t 1 } ( 1 U ) 2 . MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xI8qiVKYPFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaqbaeaabmqaaaqaaiabdAeagnaaBaaaleaacqWGbbqqcqWGHbqycqGGSaalcqWG0baDaeqaaOGaeyypa0Jaei4EaSNaeiikaGscfa4aaSaaaeaacqWGWbaCdaWgaaqaaiabdgeabjabdgeabbqabaGaeiikaGIaeGymaeJaeyOeI0IaeqiVd0MaeiykaKIaey4kaSIaemiCaa3aaSbaaeaacqWGbbqqcqWGHbqyaeqaaiabdAha2bqaaiabdgeabnaaBaaabaGaemiDaqhabeaaaaGccqGGPaqkcqGHxdaTcqGGOaakjuaGdaWcaaqaaiabdchaWnaaBaaabaGaemyqaeKaemyyaegabeaacqGGOaakcqaIXaqmcqGHsislcqWG2bGDcqGGPaqkcqGHRaWkcqWGWbaCdaWgaaqaaiabdgeabjabdgeabbqabaGaeqiVd0gabaGaemyyae2aaSbaaeaacqWG0baDaeqaaaaakiabcMcaPiabcUfaBLqbaoaalaaabaGaeGymaedabaGaemOta4Kaemyqae0aaSbaaeaacqWG0baDcqGHsislcqaIXaqmaeqaaaaakiabgUcaRiabbIcaOiabigdaXiabgkHiTKqbaoaalaaabaGaeGymaedabaGaemOta4Kaemyqae0aaSbaaeaacqWG0baDcqGHsislcqaIXaqmaeqaaaaakiabbMcaPiabdAeagnaaBaaaleaacqWGbbqqcqWGbbqqcqqGSaalcqqGGaaicqWG0baDcqGHsislcqaIXaqmaeqaaOGaeeyxa0fabaGaee4kaSIaeeikaGscfa4aaSaaaeaacqWGWbaCdaWgaaqaaiabdggaHjabdgeabbqabaGaeiikaGIaeGymaeJaeyOeI0IaeqiVd0MaeiykaKIaey4kaSIaemiCaa3aaSbaaeaacqWGHbqycqWGHbqyaeqaaiabdAha2bqaaiabdgeabnaaBaaabaGaemiDaqhabeaaaaGccqqGPaqkcqGHxdaTcqqGOaakjuaGdaWcaaqaaiabdchaWnaaBaaabaGaemyyaeMaemyyaegabeaacqGGOaakcqaIXaqmcqGHsislcqWG2bGDcqGGPaqkcqGHRaWkcqWGWbaCdaWgaaqaaiabdggaHjabdgeabbqabaGaeqiVd0gabaGaemyyae2aaSbaaeaacqWG0baDaeqaaaaakiabbMcaPiabbUfaBLqbaoaalaaabaGaeGymaedabaGaemOta4Kaemyyae2aaSbaaeaacqWG0baDcqGHsislcqaIXaqmaeqaaaaakiabgUcaRiabbIcaOiabigdaXiabgkHiTKqbaoaalaaabaGaeGymaedabaGaemOta4Kaemyyae2aaSbaaeaacqWG0baDcqGHsislcqaIXaqmaeqaaaaakiabbMcaPiabdAeagnaaBaaaleaacqWGHbqycqWGHbqycqqGSaalcqqGGaaicqWG0baDcqGHsislcqaIXaqmaeqaaOGaeeyxa0fabaGaee4kaSscfa4aaSaaaeaacqaIYaGmcqWGWbaCdaWgaaqaaiabdgeabjabdgeabbqabaGaemiCaa3aaSbaaeaacqWGHbqycqWGbbqqaeqaaiabeY7aTjabcIcaOiabigdaXiabgkHiTiabeY7aTjabcMcaPiabgUcaRiabikdaYiabdchaWnaaBaaabaGaemyqaeKaemyyaegabeaacqWGWbaCdaWgaaqaaiabdggaHjabdggaHbqabaGaemODayNaeiikaGIaeGymaeJaeyOeI0IaemODayNaeiykaKIaey4kaSIaeiikaGIaemiCaa3aaSbaaeaacqWGbbqqcqWGbbqqaeqaaiabdchaWnaaBaaabaGaemyyaeMaemyyaegabeaacqGHRaWkcqWGWbaCdaWgaaqaaiabdggaHjabdgeabbqabaGaemiCaa3aaSbaaeaacqWGbbqqcqWGHbqyaeqaaiabcMcaPiabcUfaBjabcIcaOiabigdaXiabgkHiTiabeY7aTjabcMcaPiabcIcaOiabigdaXiabgkHiTiabdAha2jabcMcaPiabgUcaRiabeY7aTjabdAha2jabc2faDbqaaiabdgeabnaaBaaabaGaemiDaqhabeaacqWGHbqydaWgaaqaaiabdsha0bqabaaaaOGaemOray0aaSbaaSqaaiabdgeabjabdggaHjabbYcaSiabbccaGiabdsha0jabgkHiTiabigdaXaqabaGccqqG9bqFcqGGOaakcqaIXaqmcqGHsislcqWGvbqvcqGGPaqkdaahaaWcbeqaaiabikdaYaaaaaGccqGGUaGlaaa@2345@

where A t and a t are calculated using equations (7) and (8).


  1. 1.

    Drummond AJ, Nicholls GK, Rodrigo AG, Solomon W: Estimating mutation parameters, population history and genealogy simultaneously from temporally spaced sequence data. Genetics. 2002, 161 (3): 1307-1320.

  2. 2.

    Leigh-Brown AJ: Analysis of HIV-1 env gene sequences reveals evidence for a low effective number in the viral population. Proc Natl Acad Sci USA. 1997, 94 (5): 1862-1865. 10.1073/pnas.94.5.1862.

  3. 3.

    Achaz G, Palmer S, Kearney M, Maldarelli F, Mellors JW, Coffin JM, Wakeley J: A Robust Measure of HIV-1 Population Turnover within Chronically Infected Individuals. Mol Biol Evol. 2004, 21 (10): 1902-1912. 10.1093/molbev/msh196.

  4. 4.

    Seo T-K, Thorne JL, Hasegawa M, Kishino H: Estimation of Effective Population Size of HIV-1 Within a Host: A Pseudomaximum-Likelihood Approach. Genetics. 2002, 160 (4): 1283-1293.

  5. 5.

    Shriner D, Shankarappa R, Jensen MA, Nickle DC, Mittler JE, Margolick JB, Mullins JI: Influence of random genetic drift on human immunodeficiency virus type 1 env evolution during chronic infection. Genetics. 2004, 166 (3): 1155-1164. 10.1534/genetics.166.3.1155.

  6. 6.

    Rouzine IM, Coffin JM: Linkage disequilibrium test implies a large effective population number for HIV in vivo. Proc Natl Acad Sci USA. 1999, 96 (19): 10758-10763. 10.1073/pnas.96.19.10758.

  7. 7.

    Haase AT: Population biology of HIV-1 infection: viral and CD4+ T cell demographics and dynamics in lymphatic tissues. Annu Rev Immunol. 1999, 17: 625-656. 10.1146/annurev.immunol.17.1.625.

  8. 8.

    Kouyos RD, Althaus CL, Bonhoeffer S: Stochastic or deterministic: what is the effective population size of HIV-1?. Trends Microbiol. 2006, 14 (12): 507-511. 10.1016/j.tim.2006.10.001.

  9. 9.

    Haddad DN, Birch C, Middleton T, Dwyer DE, Cunningham AL, Saksena NK: Evidence for late stage compartmentalization of HIV-1 resistance mutations between lymph node and peripheral blood mononuclear cells. AIDS. 2000, 14 (15): 2273-2281. 10.1097/00002030-200010200-00008.

  10. 10.

    Fulcher JA, Hwangbo Y, Zioni R, Nickle D, Lin X, Heath L, Mullins JI, Corey L, Zhu T: Compartmentalization of human immunodeficiency virus type 1 between blood monocytes and CD4+ T cells during infection. J Virol. 2004, 78 (15): 7883-7893. 10.1128/JVI.78.15.7883-7893.2004.

  11. 11.

    Wong JK, Ignacio CC, Torriani F, Havler D, Fitch NJS, Richman DD: In vivo compartmentalization of human immunodeficiency virus: evidence from the examination of pol sequences from autopsy tissues. J Virol. 1997, 71 (3): 2059-2071.

  12. 12.

    Shriner D, Liu Y, Nickle DC, Mullins JI: Evolution of intrahost HIV-1 genetic diversity during chronic infection. Int J Org Evolution. 2006, 60 (6): 1165-1176.

  13. 13.

    Barton NH: Genetic hitchhiking. Philos Trans R Soc Lond B Biol Sci. 2000, 355 (1403): 1553-1562. 10.1098/rstb.2000.0716.

  14. 14.

    Mansky LM: Forward mutation rate of Human Immunodeficiency Virus Type 1 in a T lymphoid cell line. AIDS Res Hum Retroviruses. 1996, 12 (4): 307-314.

  15. 15.

    O'Neil PK, Sun G, Yu H, Ron Y, Dougherty JP, Preston BD: Mutational analysis of HIV-1 long terminal repeats to explore the relative contribution of reverse transcriptase and RNA polymerase II to viral mutagenesis. J Biol Chem. 2002, 277 (41): 38053-38061. 10.1074/jbc.M204774200.

  16. 16.

    Shriner D, Rodrigo AG, Nickle DC, Mullins JI: Pervasive genomic recombination of HIV-1 in vivo. Genetics. 2004, 167 (4): 1573-1583. 10.1534/genetics.103.023382.

  17. 17.

    Zhuang J, Jetzt AE, Sun G, Yu H, Klarmann G, Ron Y, Preston BD, Dougherty JP: Human Immunodeficiency Virus Type 1 Recombination: Rate, Fidelity, and Putative Hot Spots. J Virol. 2002, 76 (22): 11273-11282. 10.1128/JVI.76.22.11273-11282.2002.

  18. 18.

    Jetzt AE, Yu H, Klarmann GJ, Ron Y, Preston BD, Dougherty JP: High rate of recombination throughout the human immunodeficiency virus type 1 genome. J Virol. 2000, 74 (3): 1234-1240. 10.1128/JVI.74.3.1234-1240.2000.

  19. 19.

    Levy DN, Aldrovandi GM, Kutsch O, Shaw GM: Dynamics of HIV-1 recombination in its natural target cells. Proc Natl Acad Sci USA. 2004, 101 (12): 4204-4209. 10.1073/pnas.0306764101.

  20. 20.

    Rhodes T, Wargo H, Hu WS: High rates of human immunodeficiency virus type 1 recombination: near-random segregation of markers one kilobase apart in one round of viral replication. J Virol. 2003, 77 (20): 11193-11200. 10.1128/JVI.77.20.11193-11200.2003.

  21. 21.

    Robertson A: Inbreeding in artificial selection programmes. Genet Res. 1961, 2: 189-194.

  22. 22.

    Woolliams JA, Bijma P: Predicting rates of inbreeding in populations undergoing selection. Genetics. 2000, 154 (4): 1851-1864.

  23. 23.

    Wray NR, Thompson R: Prediction of rates of inbreeding in selected populations. Genet Res. 1990, 55 (1): 41-54.

  24. 24.

    Ewens WJ: On the concept of the effective population size. Theor Popul Biol. 1982, 21: 373-378. 10.1016/0040-5809(82)90024-7.

  25. 25.

    Kimura M, Crow JF: The measurement of effective population numbers. Evolution. 1963, 17: 279-288. 10.2307/2406157.

  26. 26.

    Wright S: Evolution in Mendelian Populations. Genetics. 1931, 16: 97-159.

  27. 27.

    Hartl DL, Clark AG: Principles of Population Genetics. 1989, Sunderland, MA: Sinauer Associates, 2

  28. 28.

    Allen TM, Altfeld M, Geer SC, Kalife ET, Moore C, O'Sullivan KM, Desouza I, Feeney ME, Eldridge RL, Maier EL, et al: 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. 2005, 79 (21): 13239-13249. 10.1128/JVI.79.21.13239-13249.2005.

  29. 29.

    Jones NA, Wei X, Flower DR, Wong M, Michor F, Saag MS, Hahn BH, Nowak MA, Shaw GM, Borrow P: Determinants of human immunodeficiency virus type 1 escape from the primary CD8+ cytotoxic T lymphocyte response. J Exp Med. 2004, 200 (10): 1243-1256. 10.1084/jem.20040511.

  30. 30.

    Liu Y, McNevin J, Cao J, Zhao H, Genowati I, Wong K, McLaughlin S, McSweyn M, Diem K, Stevens C, et al: Selection on the human immunodeficiency virus type 1 proteome following primary infection. J Virol. 2006, 80 (19): 9519-9529. 10.1128/JVI.00575-06.

  31. 31.

    Liu Y, McNevin J, Zhao H, Tebit DM, Troyer RM, McSweyn M, Ghosh AK, Shriner D, Arts EJ, McElrath MJ, et al: Evolution of HIV-1 CTL epitopes: Fitness-Balanced Escape. J Virol. 2007, 81 (22): 12179-12188. 10.1128/JVI.01277-07.

  32. 32.

    Jung A, Maier R, Vartanian J-P, Bocharov G, Jung V, Fischer U, Meese E, Wain-Hobson S, Meyerhans A: Multiply infected spleen cells in HIV patients. Nature. 2002, 418 (6894): 144-10.1038/418144a.

  33. 33.

    Chun T-W, Carruth L, Finzi D, Shen X, DiGiuseppe JA, Taylor H, Hermankova M, Chadwick K, Margolick J, Quinn TC, et al: Quantification of latent tissue reservoirs and total body viral load in HIV-1 infection. Nature. 1997, 387: 183-188. 10.1038/387183a0.

  34. 34.

    Nei M: Molecular Evolutionary Genetics. 1987, New York: Columbia University Press

  35. 35.

    Kimura M, Crow JF: The Number of Alleles That Can Be Maintained in a Finite Population. Genetics. 1964, 49: 725-738.

Download references


We thank James I. Mullins for his guidance and critical comments on an early draft, two anonymous reviewers for constructive criticisms, and Reneé Ireton for editing the final version. The authors were supported by grants from the NIH (R03 AI055394, R01 HL072631, P01 AI57005, R01 AI058894, and RO1 AI047734), the University of Washington Center for AIDS Research (NIH grant P30 AI27757), and a gift from the Frank H. Jernigan Charitable Foundation.

Author information

Correspondence to Yi Liu.

Additional information

Authors' contributions

YL and JM jointly conceived the study. YL derived the equations, wrote the computer code, performed the computational experiments, and drafted the manuscript. JM advised on the study design, participated in the analysis of the mathematical and computational data, and helped draft the manuscript. Both authors have read and approved the paper.

Authors’ original submitted files for images

Below are the links to the authors’ original submitted files for images.

Authors’ original file for figure 1

Authors’ original file for figure 2

Authors’ original file for figure 3

Rights and permissions

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and Permissions

About this article

Cite this article

Liu, Y., Mittler, J.E. Selection dramatically reduces effective population size in HIV-1 infection. BMC Evol Biol 8, 133 (2008) doi:10.1186/1471-2148-8-133

Download citation


  • Recombination Rate
  • Effective Population Size
  • Recurrent Selection
  • Inbreeding Coefficient
  • Initial Frequency