Set point viral load in HIV patients ranges over several orders of magnitude and is a key determinant of disease progression in HIV. A number of recent studies have reported high heritability of set point viral load implying that viral genetic factors contribute substantially to the overall variation in viral load. The high heritability is surprising given the diversity of host factors associated with controlling viral infection. Here we develop an analytical model that describes the temporal changes of the distribution of set point viral load as a function of heritability. This model shows that high heritability is the most parsimonious explanation for the observed variance of set point viral load. Our results thus not only reinforce the credibility of previous estimates of heritability but also shed new light onto mechanisms of viral pathogenesis.
The time course of viral load in HIV infected patients follows a characteristic pattern. During primary infection the viral load rapidly grows to very high levels. The peak viremia is attained within the first few weeks of infection. Thereafter the viral load declines rapidly over a period of several months and eventually settles down at a much lower level referred to as the viral set point. Set point viral load (spVL) is a central characteristic of the course of the disease. Firstly, the virus load measurements do fluctuate in patients, the time average of the viral load remains remarkably close to the spVL in most of patients over the time scale of several years [1,2]. Secondly, higher spVL is associated with faster disease progression [3]. The stability of spVL within patients is in strong contrast to the enormous variation in spVL observed between patients. While variation in spVL between patients ranges over 3–4 orders of magnitude [3–6], the time trend over longitudinal viral load measurements typically changes by less then 0. 1 log per year [1,2]. Given that spVL is a key predictor of disease progression, there is considerable interest in identifying the host and viral genetic factors underlying the variation in spVL. A well known example for the influence of naturally occurring variation in human genetic factors on viral load is the Δ32 deletion in the CCR5 gene [7]. Moreover polymorphisms in HLA-B and C alleles have been associated with variance in virus load and genome-wide association studies (GWAS) showed that about 20% of the variance in log spVL can be attributed to specific single nucleotide polymorphisms [8–11]. 20% is likely a lower bound for the overall contribution of host genetic factors, because GWAS generally suffer from the problem that they can only identify common genetic variants with strong effects and do not account for epistatic effects between host genes [12]. Natural variation in the virus can also affect spVL. For example the transmission of a nef-deficient virus through a contaminated blood sample resulted in a low viral load in the recipients [13]. Moreover, several studies have reported a correlation between predicted replicative capacity and viral load [14–16]. As this prediction is based only on the viral genotype a patient carries, this implies that naturally occurring variation in viruses does affect viral load. A number of recent studies attempted to estimate the contribution of the viral genotype to the variation in spVL by quantifying the statistical association of viral load between donors and recipients either directly in donor-recipient pairs or through phylogenetic analysis [17–22]; for reviews see Müller et al. [23] and Fraser et al. [24]. A meta-analysis of previously published donor-recipient studies correcting for various co-factors such as age and sex yielded a heritability of 33% with a 95% confidence interval of 20–46% [24]. The two studies that inferred heritability based on phylogenetic methods provided the most extreme estimates with 5. 7% reported by Hodcroft et al. [22] and 59% reported by Alizon et al. [21]. While the phylogenetic approaches have an advantage over the donor-recipient based approaches in that they can use much larger patient populations, it is currently unclear to what extent the underlying assumptions of the phylogenetic approaches of no selection and high frequency of sampling affect the robustness of these results. The discrepant estimates call for a better quantitative understanding of the underlying factors determining heritability of log spVL in HIV. To this end we develop here a quantitative model that describes the change of the distribution of log spVL in a patient population in relation to heritability over a full transmission cycle. The model extends the approach of Shirreff et al. [25] and is similar in spirit to the integral projection models in ecology that are used to describe the temporal changes of distributions of a continuous phenotypic trait in populations [26–28]. In contrast to many applications in ecology, the application to distributions of log spVL has the advantage that all relevant processes and populations for which data are available, are numerically well approximated by a Gaussian function. This fact enables us to obtain complete analytical understanding of how spVL changes through time based on a model parametrized by available data.
We consider the change of the spVL distribution over one full reproduction cycle on the epidemiological level, i. e. from the current to the next generation of patients. We divide the patient population into “carriers” (HIV infected individuals prior to selection for transmission), “donors” (individuals that have been selected for transmission) and “recipients” (individuals that have just been infected by donors). Furthermore, we divide the reproduction cycle into three steps: (i) selection of donors from the carriers with replacement according to their transmission potential, (ii) transmission from donors to recipients, and (iii) intrahost evolution of the virus from the start of infection to the next transmission. Finally, we explicitly distinguish between factors contributing to set point viral load with regard to being transmissible (i. e. viral genetic factors) versus being non-transmissible (i. e. host genetic factors, environmental factors, or any interaction between host, virus, and the environment). A schematic overview over the effects of these steps on the distribution of log spVL is shown in Fig. 1. In the Supplementary Materials we show how the change of the spVL distribution can be computed for any distribution over a full transmission cycle. If all populations and processes are well approximated by Gaussian functions, then an approximation to the resulting log spVL distributions can be computed analytically (see Methods and Supplementary Materials). Assuming that the population is in equilibrium we obtain for the mean, M̃C, variance, ṼC, and heritability, h2, the following expressions: M ˜ C = μ o + μ i (1+ ν e + ν o V ˜ C − ν e), (1) V ˜ C = ν t + ν i 2 (1+ 1+4 ν e + ν o ν t + ν i) + ν e, (2) h 2 =1− ν e V ˜ C. (3) Here the parameters μo and νo characterize the transmission potential [6], i. e. the overall probability of a patient to transmit the infection as a function of log spVL (Fig. 1 (o) – (i) ). This transmission potential is given by the product of the rate of transmission per contact and the disease duration. As the former increases and the latter decrease with increasing spVL, the transmission potential has a maximum at intermediate levels of spVL [6]. The parameter νe gives the variance of the contribution of host/environmental effects on log spVL. The parameter νt describes the variance due to the bottleneck at transmission from donor to recipient, as a founder strain is selected randomly from the diverse population in the donor (Fig. 1 (i) – (ii) ). The parameters μi and νi describe the mean and variance of the contribution of intrahost evolution to log spVL (Fig. 1 (ii) – (iii) ). Our model assumes that the bottleneck at transmission is neutral with regard to selection on set point viral load. Note, that the assumption is without loss of generality. This is important because there is evidence for selection at transmission [29], although it is unclear whether selection acts on spVL. Any selective effect at transmission, however, can be subsumed into the parameter μi. Hence, the effect of selection is effectively incorporated in our model. The parameter for the mean contribution by the host/environment, μe, does not appear in equations 1 or 2. This is because the equations refer to the phenotypic value of spVL, i. e. the sum of the genetic contributions of the virus and the contributions from the host/environment. Any large environmental/host effect on the mean can always be compensated by correspondingly strong genetic effect of the virus on the mean but with opposite sign. The above results are applicable if, (i) if the population is approximately in equilibrium, and (ii) all populations and processes are numerically well approximated by Gaussian functions. Assumption (i) has been discussed in detail previously [6,25,30,31]. In essence, this assumption is supported by three observations. Firstly, the mean of the spVL distribution coincides with the optimum of the transmission potential (see Fig. 2 and Fraser et al. [6]). Secondly, the rate of change of spVL has decreased over the last 25 years [31]. Thirdly, the rate of evolution is sufficiently rapid such that a spVL that is optimal for transmission could have evolved over the course of the epidemic [25]. These findings suggest that the distribution of set point viral load is indeed approximately in equilibrium, which in turn makes it is plausible to assume that the environmental and genetic factors determining set point viral load are also in equilibrium. Regarding assumption (ii), we note that a Gaussian function describes a distribution or process by a main effect (mean) and some variational noise (variance). Thus in absence of any better knowledge, a Gaussian distribution is a natural starting point to describe any process and simply represents a second order approximation to an unknown distribution. We can assess the validity of describing the distributions of spVL in carriers and the transmission potential graphically using available data. Inspection of Fig. 2A and Figure S1 in Supplementary S1 Text shows that the viral load amongst carriers is indeed numerically well approximated by a Gaussian with mean log spVL, MC ≈ 4. 5, and variance in log spVL, VC ≈ 0. 5. Also the fit of a Gaussian to the transmission potential (see Fig. 2B) is a very good approximation (mean μo ≈ 4. 6 and variance νo ≈ 1. 0), even though the transmission potential as estimated by Fraser et al. [6] is slightly right-skewed. There are no data to inform the shape of the processes of transmission and intrahost evolution. Using a description that has a mean effect with some variation around this mean is natural. Nonetheless, we test the effect of numerical deviations from a Gaussian with the following simulations. Firstly, we use the exact right-skewed transmission potential as given by Fraser et al. [6]. The analytical approximations for the distribution of the population in equilibrium remain excellent when the substantial deviations of the transmission potential from a Gaussian are incorporated (see Figure S2). Secondly, we study the robustness towards deviations from Gaussian functions in the processes describing intrahost evolution and the transmission bottleneck. Even when both processes are strongly skewed, the analytical approximations for mean and variance are excellent (typically less than 2% deviation, see Figure S3 in Supplementary S1 Text). To assess what heritability values are compatible with the observed mean and variance of log SPVL in the carrier population we take a simple approach that is in essence Approximate Bayesian Computing with rejection sampling. To this end we define plausible prior distributions for the parameters of the model. Sampling randomly from the priors we determine the resulting means and variances of log spVL in carriers and reject sets of parameters that lead to means and variances outside a defined permissible range. The set of accepted parameters gives the posterior distribution. For the range of permissible mean log spVL we assume 4 < M̃C < 5, which is compatible but somewhat larger than the observed range in the studies reported by Fraser et al. [6] and Geskus et al. [5] (see Fig. 2A and Supplementary Materials, Section E). For the permissible range of variances of log spVL we assume that 0. 3 < ṼC < 0. 8, which again is compatible but somewhat larger than the values reported by Fraser et al. [6] and Geskus et al. [5] (see Supplementary Materials, Section E). We use uniform priors for all parameters. The parameters μo and νo, which describe mean and variance of the transmission potential, have thus far only been estimated only by a single peer reviewed study (Fraser et al. [6] and Fig. 2B; see also [32]). To account for uncertainty in the estimates of these parameters we use 4 < μo < 5 and 0. 5 < νo < 1. 5. Estimates for remaining parameters cannot be easily derived from the existing literature. To account for uncertainty in these parameters we assume 0 < νe < 1; −1 < μi < 1; 0 < νi < 0. 3 and 0 < νt < 0. 3. Fig. 3 shows the posterior parameter distribution from the rejection sampling. Different colors in the scatter plots indicate different levels of mean heritability at given parameter combinations. The contour lines show the density of posterior distribution. The key result shown in the figure is that the majority of accepted parameter values result in high values of heritability (purple to orange color at contour lines of highest posterior densities). While low values of heritability are also compatible with the observed mean and variance of log spVL, they occur rarely in the posterior sample and are at the edges of the prior distributions (red to blue areas). The center of mass of the posterior sample is in areas with high heritability, higher in fact then what would seem compatible with current estimates of heritability and host genetic factors. There are two factors not included in this analysis: measurement error of spVL and prior knowledge of the host contribution to spVL. Increasing the measurement accuracy of spVL would increase heritability estimates based on both donor-recipient pairs and phylogenetic inference. Incorporating prior knowledge of the host genetic contribution would set an upper bound on the estimates of heritability in our analysis. Thus accounting for these two factors bring the center of mass of the heritability distribution closer to the measured values of heritability. The figure also highlights that generally wider priors would not change the posterior distribution because parameter values at the upper end of priors are never accepted. The change of the mean capacity of the virus to induce spVL through intrahost evolution, μi, is restricted to values smaller than 0. 6 and decreases with increasing variance generated by the host/environment effects, νe. Increasing νe corresponds to decreasing heritability (see eq. 3) and thus high levels of μi require high levels of heritability. The center of mass of the posterior sample suggests that the most parsimonious explanation of the observed mean and variance of log spVL implies both small intrahost evolution and high heritability. One criticism leveled against the transmission potential as quantified in [6] is that it does not appropriately reflect transmissions occuring during the acute or the AIDS phase. In the Supplementary Material, Section F. 3, we show that our quantitative results are robust towards using a corrected transmission potential.
The above analysis shows that the most parsimonious explanation of the observed distribution of spVL in HIV carrier populations requires high heritability of spVL. Although low heritability values are also compatible with the observed distribution of spVL in HIV carrier populations, parameter combinations resulting in these low values have a small probability and occur at the edge of the realistic parameter range. The skepticism with which the estimated heritability values have been met in the field suggests that the general expectation is that heritability of spVL should be low. In contrast, our analysis shows that high heritability values are not only compatible with, but are also the more parsimonious explanation of the observed distribution in spVL in HIV carrier populations. Low heritability only occurs if the processes of intrahost evolution and the transmission bottleneck have a weak effect on spVL, i. e. if the parameters μi, νi and νt are small. An intuition can be obtained by noting that in equilibrium the variance generating and variance eliminating processes balance out. The transmission potential only exerts weak selection on log spVL and therefore only marginally reduces variance. The decrease of variance by selection for transmission has to be compensated by an increase in variance by intrahost evolution and the transmission bottleneck. For too low heritability, the genetic variance generated by intrahost evolution and transmission bottlenecks would overwhelm the reduction of variance due to selection by the transmission potential. While there are to our knowledge no data that allow to estimate the variance generated at transmission, νt, the posterior distributions of μi and νi are broadly compatible with the observed changes of virus load within patients [1,2, 33]. Taken together our analysis suggests that the most parsimonious explanation of the distribution of log spVL is high h2 but low νt, μi and νi. Hence, heritability is high while the processes of intrahost evolution and transmission bottleneck have a small effect on the capacity of the virus to modulate log spVL. High heritability implies a substantial genetic control of the log spVL by the virus. The observation that at the same time the contribution of intrahost evolution to spVL is small raises an interesting question: How can a strongly heritable trait show little intrahost evolution? Given the otherwise ample evidence for rapid intrahost evolution of HIV such as escape from drugs or the immune response, the absences of intrahost evolution of spVL is surprising. Generally a trait is expected to respond to selection, if (i) the trait is heritable, (ii) there is phenotypic variation of the trait in a population, and (iii) the trait is linked to fitness. That spVL is heritable has been reported previously [23,24] and our analysis reinforces the credibility of these findings. That there is phenotypic variation in the control of spVL by the virus is plausible given the large genetic variation of the virus population within an individual. What remains is whether it is conceivable that the capacity of a viral genotype to induce spVL is only weakly linked to fitness. One hypothesis that could reconcile high heritability with little intrahost evolution is that variation in viral load between patients is in part due to virus-induced activation of target cells. Difference in activation rate of target cells has previously been argued to account for a substantial part of the variation in viral load [4]. Furthermore, if target cell activation is at least partially under the control of the virus, then this control may indeed be weakly linked to intrahost fitness. If the target cell activation is systemic (i. e. not locally confined to the inducing virus) then increased target cell activation increases the pool of susceptible cells, but the benefit of increased target cell activation is not confined to the producer virus. As a result selection for virus induced activation rate is expected to be neutral or nearly neutral [34]. Indeed, an explicit model of the evolution of log spVL for a virus induced control of target cell activation can reconcile high heritability with absence of intrahost evolution [35]. Our modeling approach is based on describing how the distribution of a continuous phenotypic trait, here log spVL, changes in a population over a full cycle of reproduction. This approach is closely related to the method of integral projection models, which has been developed and widely applied in ecology and population biology [26–28,36,37]. The approach can in principle describe how arbitrary distributions change over time as a function of processes such as selection and reproduction. Here we are able to obtain a full analytical description of the temporal change of the spVL distribution, because all relevant distributions and processes can be well approximated by Gaussian functions. We also show that our analytical results remain robust even for substantial deviations numerical deviations from Gaussian functions (see Supplementary Materials, Section F). Moreover, the model can be parametrized on the basis of available data. There are ample data for mean and variance of spVL and also most of the parameters can be confined to plausible ranges based on the literature. Our study clearly supports that high heritability is compatible with the observed distribution of log spVL in HIV carriers. High heritability of spVL does not preclude that also the host genotype has a considerable effect on virus load. However, it does lead to the expectation that over the course of infection the capacity to induce higher spVL should increase considerably unless this capacity is only weakly linked to intrahost fitness. This sheds new light onto the mechanisms controlling viral load. There should be identifiable genetic variation in the virus population that is associated with viral load, and moreover, the loci associated with control of viral load should be weakly linked to intrahost fitness. Genome-wide association studies mapping viral genetic polymorphisms to variance in log spVL seem a natural approach to test this prediction. A recent study by Bartha et al. [38] was unable to identify any statistical associations, but was powered only to detect individual non-synonymous mutations with an effect size of >4% on heritability. Larger studies will thus be required to identify whether and which viral polymorphisms are associated with set point viral load.
The spVL in a patient is generally determined by viral genetic factors, host genetic factors, the environment and interactions between these factors. Since only the virus is transmitted from donors to recipients, we subsume all non-transmissible effects such as the host genetic factors, environmental effects and all interactions between host, virus and the environment generically under “environmental effects”, e. The transmissible effects due to the viral genotype are the “genotypic effects”, g. The “phenotype” spVL is then given by g+e. We assume that the distribution of log spVL in the carrier population is given by a normal distribution N (M C, V C), where MC and VC are the mean and variance, respectively. The transmission potential, defined as the overall probability of transmission of an HIV carrier integrated over the entire course of the disease, is assumed to be a function of log spVL which can be well approximated by a normal distribution N (μ o, ν o) (see Fig. 2). Here μo is the log spVL at which the transmission potential is maximal and νo characterizes how strongly the transmission potential selects for transmission at μo. We assume that g and e are independent and normally distributed in the carrier population with N (m C, v C) and N (μ e, ν e), respectively. Here mC and vC are the variables that describe the mean and variance of the distribution of viral genotypes in the carrier population. Note that here the independence of g and e refers to the quantitative contribution of virus and host to spVL. Importantly, this independence does not imply an absence of virus genotype by host genotype interactions, such as an interaction between a particular viral epitope and a host HLA molecule. Genotype by genotype interactions are non-transmissible and thus subsumed in e. The parameters μe and νe describe mean and variance of the distribution of environmental effects, which comprise host effects, interactions and any non-transmissible effect. The distribution of phenotype log spVL in the carrier population is then given by a normal distribution with mean and variance, M C = m C + μ e, and V C = v C + ν e. (4) Selection for transmission acts on log spVL, i. e. on the sum of the genotypic and environmental effects, and is given by the transmission potential. Specifically, the probability of transmission for a given log spVL, ϕ, is given by (see Fig. 2B), S (ϕ) = 1 2 π ν o e - (ϕ - μ o) 2 2 ν o. (5) Applying the above transmission potential to the carrier population, we find that the genotype and phenotype in the donor population are again normally distributed (see Supplementary Materials, Equations B6 and B7). The donor genotype has mean and variance, m D = m C (ν e + ν o) + (μ o - μ e) v C v C + ν e + ν o, and v D = v C (ν e + ν o) ν e + ν o + v C. (6) The donor phenotype has mean and variance (see Supplementary Materials, Equations B8 and B9), M D = M C ν o + μ o V C ν o + V C, and V D = V C ν o ν o + V C. (7) Note, that the mean and variance of the environmental effects (i. e. the host effect) is not given by the differences between the phenotypic and genotypic values, because environment and genotype in the donors are correlated. This is because selection for transmission acts on the sum of environmental and genotypic effects. In other words selection for transmission selects a subset of viral genotypes and host genotypes, and host and viral genotypes are correlated, because selection operates on their combined effect. When the virus is transmitted from the donor to the recipient population, the virus is “harvested” from a non-random distribution of environmental effects (and thus also from a non-random set of hosts). The harvested virus is then redistributed over a random set of new hosts/environments in the recipient population. Thus all environmental effects in the donor population are erased at transmission and the environmental contribution in the recipients is redrawn from N (μ e, ν e). To account for the fact that the virus population experiences a strong bottleneck from recipient to donor, we assume that the viral genotype is not transmitted exactly from donor to recipient but instead is assumed to be randomly drawn out of a distribution of genotypes in the donor patient. Assuming that this distribution is normal with mean mD and variance νt we obtain that both genotype and phenotype in the donor population are normally distributed. The recipient genotype has mean and variance, m R = m D, and v R = v D + ν t. (8) The recipient phenotype has mean and variance, M R = m R + μ e 0, and V R = v R + ν e 0, (9) where μ e 0 and ν e 0 are the mean and variance of the host/environmental effects prior to infection. The environmental effects are redrawn randomly, because they are not inherited from one transmission to the next. Note, that we assume here that the bottleneck at transmission is neutral. This assumption does not imply that there is no selection at the transmission stage, but rather that the bottleneck is neutral with regard to the spVL that the transmitted strains will eventually cause. Any selection at and after transmission on the viral genotypic contribution to log spVL is subsumed in the next step, intrahost evolution. After transmission the virus population in the recipient may change in a directed fashion according to intrahost evolution. Assuming that the overall change of the viral genotype due to intrahost evolution can be approximated by a normal distribution we find that the distribution of genotypes and phenotypes in the next generation of carriers, C′ is again normal. The distribution of the genotypes has a mean and variance m C' = m R + μ i, and v C' = v R + ν i. (10) The parameter μi thus describes any genetic change in the virus that affects log spVL across all patients in the same way. The parameter νi describes genetic changes that affect log spVL in a manner that is specific to the patient, i. e. it describes the effect of changes of log spVL due to genetic interactions between virus and host. As the environmental effects comprise the immune response by the host, the mean and variance in environmental effects may change in coevolution with the virus through μ e i and ν e i, respectively. Thus we obtain for mean and variance of the distribution of phenotypes, M C' = m C' + μ e 0 + μ e i, and V C' = v C' + ν e 0 + ν e i. (11) Note, that any selection for spVL at the transmission bottleneck can now be interpreted as a genotypic change that occurs during intrahost evolution. Thus the overall model is appropriate both for non-selective and selective bottlenecks. Heritability, h2, is defined as fraction of genotypic variance relative to phenotypic variance in the carrier population [39]. Thus we have, h 2 = v C V C = 1 - ν e V C. (12) Heritability can be estimated in a parent-offspring regression [39], where h2 is equal to the regression slope b. Donor-recipient pairs can be seen as parent-offspring pairs, where care must be taken since the donors are not randomly selected from the carrier population but are selected according to the transmission potential. Since, however, we are measuring the heritability of spVL and donors are selected based on spVL, the regression of recipients on selected donors is equal to heritability of spVL in carriers [24,39]. We now have a complete analytical description how mean and variance of log spVL change from the current to the next generation of carriers. The fact that the log spVL that maximizes the transmission potential and the mean of the distribution of log spVL in the carrier populations (see Fig. 2 and Fraser et al. [6]) are both around 4. 5, we can assume that the process is roughly at equilibrium. In equilibrium we have that the mean and variance of the distribution of phenotypes does not change, i. e. MC′ = MC and VC′ = VC. This will be fulfilled if the genetic and environmental contributions are also at equilibrium, implying in particular that μ e = μ e 0 + μ e i and ν e = ν e 0 + ν e i (see Supplementary Materials section C. 1). Using Equation 12 we can express the equilibrium mean and variance of log spVL as a function of νe, the variance of the contribution of the host/environment to log spVL (see Supplementary Materials, Equations C9 and C10), M ˜ C = μ o + μ i (1+ ν e + ν o V ˜ C − ν e), (13) V ˜ C = ν t + ν i 2 (1+ 1+4 ν e + ν o ν t + ν i) + ν e. (14) or as a function the heritability h2 (see Supplementary Materials, Equations C12 and C13), M ˜ C = μ o + μ i (1+ (1− h 2) V ˜ C + ν o h 2 V ˜ C), (15) V ˜ C = ν t + ν i 2 (h 2) 2 (1+ 1+ 4 (h 2) 2 ν o ν t + ν i). (16)
Following an initial peak in viremia, the viral load in HIV infected patients settles down to a set point which remains more or less stable during chronic HIV infection. This set point viral load is one of the key factors determining the rate of disease progression. The extent to which it is determined by the virus versus host genetics is thus central to developing a better understanding of disease progression. Here we develop an analytical model that describes the changes of the distribution of set point viral load in the HIV carrier population over a full cycle of transmission. Applying this model to patient data we find that the most parsimonious explanation for the observed large variation of set point viral load across HIV patients is that set point viral load is highly heritable from donors to recipients. This implies that set point viral load is to a considerable extent under the genetic control of the virus.
High Heritability Is Compatible with the Broad Distribution of Set Point Viral Load in HIV Carriers