An evolutionary role for HIV latency in enhancing viral transmission.

An evolutionary role for HIV latency in enhancing viral transmission.

Products Related to ZikaWestNileDengueMalariaT.BChikungunya, HIV, SARS

Product# 1039 HIV-1 YU2 (M Tropic)Envelope Glycoprotein gp41

Product# 1081 HIV-1 gp120 (ADA)

Product# 1011 HIV-1 gp120 (subtype C)

Product# 1031 HIV-1 gp120 (YU2)

 

SUMMARY

HIV latency is the chief obstacle to eradicating HIV but is widely believed to be an evolutionary accident providing no lentiviral fitness advantage. However, findings of latency being “hardwired” into HIV’s gene-regulatory circuitry appear inconsistent with latency being an evolutionary accident, given HIV’s rapid mutation rate. Here, we propose that latency is an evolutionary “bet-hedging” strategy whose frequency has been optimized to maximize lentiviral transmission by reducing viral extinction during mucosal infections. The model quantitatively fits the available patient data, matches observations of high-frequency latency establishment in cell culture and primates, and generates two counterintuitive but testable predictions. The first prediction is that conventional CD8-depletion experiments in SIV-infected macaques increase latent cells more than viremia. The second prediction is that strains engineered to have higher replicative fitness—via reduced latency—will exhibit lower infectivity in animal-model mucosal inoculations. Therapeutically, the theory predicts treatment approaches that may substantially enhance “activate-and-kill” HIV-cure strategies.

Graphical Abstract

 

INTRODUCTION

HIV actively replicates in CD4+ T lymphocytes but can also enter a long-lived quiescent state termed proviral latency in memory CD4+ T cells (; ). The population of latently infected cells is relatively small in patients (~1 in 106 CD4+ T cells) and does not generate significant viral RNA (). However, latently infected cells provide a critical viral reservoir, which enables lentiviral persistence even during prolonged antiretroviral therapy (ART). Further, if patients interrupt ART, persisting latent viruses reactivate, driving HIV to pre-treatment viral loads within weeks (). Consequently, latency is the chief barrier to a curative HIV therapy.

While latency enables HIV to avoid extinction during ART, the benefit of latency prior to the ART era—during the centuries of natural lentiviral infections—remains unclear. In fact, latency appears to have been deleterious prior to ART since latently infected cells produce no virus and decrease patient viral loads. Given latency’s reduction of lentiviral replicative fitness, the prevailing hypothesis is that latency is an evolutionary accident—an epiphenomenon that only results when lentiviruses infectCD4+ T cells that are transitioning from activated to quiescent memory states (; ; ). Latency is therefore viewed to be an infrequent bystander effect that only occurs after a viral-driven adaptive immune response initiates and CD4+ T lymphocytes begin to form memory subsets. Yet, a recent study in Rhesus macaques indicates that latency reaches high levels within the first 3 days of infection (), which is prior to the generation of an SIV-specific adaptive immune response ().

If latency were a non-beneficial viral trait or epiphenomenon, one would expect it to have been lost due to natural selection or genetic drift, given lentiviruses’ rapid evolutionary rates. Yet, a companion study ( [this issue of Cell]) demonstrates that the ability to establish latency is “hardwired” into HIV’s gene-regulatory circuitry. This study matches recent data showing that ~50% of cell-culture infections—in which adaptive immune responses are absent—result in lentiviral latency (; ). Further, HIV’s auto-regulatory Tat circuit appears optimized to amplify stochastic fluctuations in viral gene expression, producing fluctuations that are sufficient to induce a probabilistic switch to latency (; ; ). In general, stochastic expression noise is thought to be selected against and thus filtered out of regulatory circuits when not beneficial (; ). The persistence of a hardwired latency circuit suggests an unknown selective advantage, which outweighs latency’s putative fitness cost of reducing long-term viral loads.

One possible selective benefit is that—by providing a long-lived viral reservoir—latency could enhance lentiviral survival during unfavorable environmental conditions. Similar “bet-hedging” hypotheses () have been proposed for bacteriophage- γ lysogeny () and bacterial persistence (). However, lentiviral latency would only provide a bet-hedging advantage if there were risks of viral extinction due to environmental fluctuations. In reality, lentiviruses appear in little danger of population crashes, as they evade immune clearance and maintain high viral loads of ~105 particles/ml of blood plasma for years (and lentiviruses clearly did not evolve under pressure from antiretroviral drugs). Further, lentiviruses only infect a small percentage (~1%–2%) of available target cells, making target-cell fluctuations unimportant during chronic infection. Nevertheless, viral loads remain low during one phase of the lentiviral lifecycle: initial mucosal infection.

The probability of successful mucosal infection is low, with <1% of unprotected sex acts between HIV-discordant couples resulting in self-propagating systemic HIV infections (; ; ). When successful infections do occur, they expand from single founder sequences (; ), indicating that only one variant in the transmitted quasispecies avoids extinction. Further, animal models of HIV capture a consistent ~6 day delay from experimental mucosal inoculation to self-propagating infection (; ), which implies that the first days of lentiviral infection provide conditions unsuitable for viral growth.

The unfavorable conditions of early lentiviral infections typically occur in the mucosa, where >90% of HIV infections initiate (). HIV’s evolutionary precursor in non-human primates (SIV) also spreads through mucosal transmission—via sexual activity or fighting with subsequent communal wound licking (). Mucosal challenge experiments in primates with large inoculations provide direct evidence that the mucosa are initially unfavorable to lentiviral growth: large inoculations of ~109 infectious units (by TCID50) initially burn out within ~5 days (). Quantitatively, each initially infected cell lives for ~1 day (), so the number of actively infected cells after 5 days scales with (R0muc)5 — wherein R0muc is the basic reproductive ratio during early mucosal infection. Since actively infected cells crash within ~5 days (), (R0muc)5 approaches 0, implying that R0muc < < 1 during initial mucosal infection.

Here, we quantitatively test the hypothesis that latency provides a bet-hedging advantage that increases the probability of successful lentiviral transmission despite reducing viral loads during systemic infection (Figure 1A). The key point is that increasing the probability of latency (plat) increases the probability that each initially infected cell survives initial mucosal infection. Yet, increasing plat also decreases viral loads in systemically infected hosts, which reduces the inoculum transmitted to new hosts. With a higher per-cell survival rate but fewer initially infected cells, the question is whether latency’s fitness benefits outweigh its costs—which would establish latency as an evolutionarily beneficial trait that is maintained by natural selection.

HIV Latency as a Bet-Hedging Strategy for Maximizing Viral Transmission

(A) Schematic of the lentiviral transmission process. Lentiviral transmission is illustrated as a two-compartment process, beginning with viral inoculation in the mucosa and progressing—in some cases—to systemic infection in the lymphoid tissue, where >98% of CD4+ T cells reside (). The parameter plat reflects the probability that an HIV-infected cell enters latency. An HIV strain incapable of entering latency (plat = 0) would generate increased viral loads during systemic infection, transferring more virions to new hosts. However, the latency-incapable virions would rapidly destroy the small CD4+ T cell population initially present in the mucosa of the new host—reducing the probability of systemic infection (upper). In contrast, an HIV strain capable of entering latency (plat > 0) would generate lower viral loads during systemic infection, transferring fewer virions to new hosts. Yet, the relatively few transferred virions would not destroy all mucosal target cells. By entering long-lived latency in some mucosal cells, the latency-capable strain would increase its probability of surviving initial infection to establish systemic infection (lower).

(B) Schematic of the two-compartment model of lentiviral transmission. The two major processes controlling the probability of lentiviral transmission (ptransmission) are: (1) the inoculum of infected cells (I0) and (2) the probability that an infected cell in the inoculum survives initial infection to establish systemic infection (pestab). (Right to left) HIV enters a host mucosal site, but due to the small number of permissive target cells in the early mucosa (prior to day 6), R0 < 1. To successfully establish systemic infection, the virus must avoid extinction until R0 > 1. Critically, the likelihood of an actively infected cell or a free viral particle surviving until day 6 to initiate systemic infection is negligible since virus-producing cells die within 40 hr of infection and viral progeny are cleared from the system ~100-fold more rapidly. In contrast, latently infected cells are long-lived and can reactivate once R0 > 1 to initiate systemic viral expansion. Therefore, despite reducing long-term viral loads, latency may increase ptransmission by increasing viral survival during initial infection. This would make latency evolutionarily beneficial at the population scale.

See also Figure S1

RESULTS AND DISCUSSION

Mathematical Models of Lentiviral Transmission and Rationale for Models

Three classes of mathematical models are developed to quantify the net impact of latency on lentiviral transmission (Figure S1). Each class of models generalizes the well-parameterized basic model of viral dynamics () to include both proviral latency and the conditions of early mucosal infection (i.e., R0muc < 1) during which latency may be critical (Experimental Procedures).

The first class of models tracks initial lentiviral infection in the mucosa alone (Extended Experimental Procedures, Section A). Given the small numbers of infected cells during initial mucosal infection, the established model of mucosal infection is stochastic (). We analyze this experimentally parameterized stochastic model—and a deterministic approximation to this model—to quantify how the probability of viral extinction in the mucosa depends on the probability of latency (plat).

The second class of models extends the single-compartment model into a two-compartment model (Figure 1B) that tracks both initial infection in the mucosa and systemic infection in the lymphoid tissue (Extended Experimental Procedures, Section B). Importantly, the initial and systemic infection model compartments only differ in a single experimentally measured parameter: R0 (Figure 1B and Table S1). Collectively, the models predict an optimal value of plat(poptlat=0.5) that matches latency frequencies measured in cell culture (; ) and is consistent with latency levels measured in mucosal primate infections (). However, the large value of plat does not match the low frequencies of latency observed in chronically infected patients (; ).

The third class of models incorporates a canonical immune response () into the two-compartment model (Extended Experimental Procedures, Section C)—since a key difference between cell-culture models and chronic infection is the presence of an adaptive immune response. Each immune parameter added is either tied to a distinct patient-measured value or has been measured previously in the literature (Table S2). With no added free parameters, the immune model fits all available patient data and predicts the same robust poptlat value.

Latency’s Net Evolutionary Impact Is the Product of Its Impact on Both Initial Infection and Systemic Infection

To calculate the optimal plat value, the two-compartment models track latency’s net evolutionary impact across both mucosal and systemic infections. While the nonlinear models are complex, we decouple latency’s net impact on viral transmission into a product of two factors: (1) the average initial inoculum of infected cells per mucosal inoculation (I0), and (2) the probability that an initially infected cell establishes systemic infection (pestab) (Figure 1A). This product can be derived analytically when the number of infected cells is Poisson distributed and when each infected cell lineage is statistically independent. Under these two assumptions, the probability of lentiviral transmission per-mucosal inoculation (ptransmission) reduces to:

ptransmission = 1 − epestabI0 ≈ pestabI0
[1]

 

The equality in Equation [1] is a direct calculation of the Poisson probability that at least one infected cell in the inoculum I0 establishes systemic infection. Critically, ptransmission < 10−2 since < 1% of lentiviral infections result in self-propagating infections (; ). Given the equality, ptransmission < 10−2 immediately implies that pestab I0 < ~10−2.

Having used the equality to establish that pestab I0 < ~10−2, we can discard the quadratic and higher-order terms in the Taylor Series expansion of epestabI0 with negligible impact. This leads to the subsequent approximation (i.e., linearization) in Equation [1]: ptransmission  pestab I0.

Given Equation [1], the overall goal of determining whether latency’s benefits outweigh its costs reduces to quantifying latency’s impact on pestab and I0.

Latency Increases the Probability that an Initially Infected Cell Survives Mucosal Infection and Establishes Systemic Infection

To quantify latency’s impact on pestab, we begin by tracking lentiviral survival during mucosal infection alone. As noted above, the first 5 days of mucosal infection are characterized by a lack of detectable actively infected cells (; ), indicating that R0 in the mucosa (R0muc) is initially < < 1 (Extended Experimental Procedures, Section D). R0muc < < 1 is also consistent with the infrequency of successful mucosal transmissions (ptransmission < 0.01) and the ~6-day delay before systemic infection when lentiviral infections do establish ().

Both deterministic differential equations models (Figure 2A) and stochastic Monte-Carlo models (Figures S2A and S2B) capture the fitness advantage of latency in the mucosa. Model simulations are performed with R0 < 1 and an inoculated dose of virus that results in a few dozen initially infected cells, matching animal mucosal experiments (; ; ). The quantitative models show that—in the absence of latency—all virions and infected cells are driven extinct in the first 5 days of mucosal infection (Figures 2A, inset, and S2A). In contrast, low levels of latency enable viral survival (Figures 2A and S2B). To test the robustness of these predictions across all R0 < 1 and I0 < 100, a continuous-time branching-process model was developed (). The branching-process model (Extended Experimental Procedures, Section A) directly computes the viral extinction probability as a function of time, providing an efficient alternative to averaging thousands of Monte-Carlo simulations for each R0 and I0. In the absence of latency, the viral extinction probability approaches 1 by day 5 of mucosal infection, except in the small slice when R0 ≈ 1 (Figures S2C and S2D)—which does not match the levels of R0 inferred from animal mucosal challenge experiments ().

An Evolutionary Optimum for Latency

(A) Numerical solutions to Equation [6] showing the dynamics of latently infected cells in early mucosal infection (R0muc = 0.25). As plat increases, the number of surviving latently infected cells increases. (Inset) The dynamics of actively infected cells in early mucosal infection showing that as plat increases, actively infected cells reach extinction more rapidly.

(B) In systemic infection, (R0LT = 10), increases in plat decrease the virus load (and, therefore, the viral dose transmitted to the next host). Dynamics in (A and B) are calculated numerically from Equation [6], using the parameters in Table S1 (r = 0).

(C) Schematic flowchart of the derivation of the (optimal) latency probability poptlat that maximizes ptransmission. Red text indicates key assumptions made at each step of the derivation. For example, R0muc < < 1 implies that the vast majority of latently infected cells during initial infection are produced in the first generation, leading to the approximation LR0>1initplatI0. The results of the analytic derivation quantify the tradeoff of latency: increasing plat linearly increases pestab but decreases I0 by the factor (1-plat). Since this tradeoff is almost equally balanced, the optimal latency probability, poptlat, approximately equals 0.5.

(D) Normalized probability of host-to-host transmission (ptransmission) as a function of plat. Results shown are obtained either analytically, from Equation [5] (magenta line), or numerically using the plateau levels of actively infected cells (I) and latently infected cells (L) simulated in A and B (magenta dots). As in C, the probability of transmission is maximized when plat ~0.5.

(E) Normalized probability of host-to-host transmission when systemic infections emerge from non-latent routes (e.g., dendritic cells) with probability fnonlatent > 0 (Equations [S12 and S13]). The maximum probability of transmission occurs at slightly lower plat values, but poptlat is still large.

See also Figure S2.

For completeness, the surviving number of mucosally infected cells was directly computed using a Wright-Fisher model (; Extended Experimental Procedures, Section A). The Wright-Fisher simulations demonstrate that the surviving number of mucosally infected cells increases approximately linearly with plat for each I0 (Figures S2E–S2G). This linear dependence can also be derived analytically. Given that R0muc < < 1 during initial mucosal infection, the majority of latently infected cells are produced in the first generation of infection (Extended Experimental Procedures, Section A). Since these cells are unlikely to reactivate during the short duration of initial infection, the number of latently infected cells that survive mucosal infection is ≈ platI0, the latent fraction of the inoculum. Thus, both simulations and analytics indicate that increasing plat approximately linearly increases the number of infected cells that survive initial mucosal infection.

Given that latency appears to increase viral survival in the early mucosa, we next tested whether latency increases the probability of systemic infection, which mainly occurs in the lymphoid tissue where >98% of CD4+ T cells reside (). To do so, the Wright-Fisher model was extended into a two-compartment model that directly captures the two typical stages of lentiviral infection: early mucosal infection and systemic (lymphoid) infection (Extended Experimental Procedures, Section B). Only a single parameter value is assumed to differ between the early mucosal and systemic infection compartments. While R0muc is parameterized to be <1, R0 during systemic infection in the lymphoid tissue (R0LT) is set to 10 to match its value in chronically infected patients ().

The two-compartment model fits the available human and animal data of early infections, showing that: (1) only a small fraction of mucosal infections result in systemic infections (), (2) successful systemic infections emerge after ~5–7 days (), and (3) systemic infections initiate from single “founder” infected cells (; ). More importantly, the two-compartment model directly shows that latency increases the probability (pestab) of systemic infection—with pestab maximized when plat > 0.6 (Figure S2H; Extended Experimental Procedures, Section E).

Latency Decreases the Inoculum in a New Host

While increasing plat increases the probability of systemic lymphoid infection for any given inoculum of initially infected cells (I0), the probability of lentiviral infection also depends on I0 itself. Critically, I0 is proportional to the viral load of the transmitting patient (Extended Experimental Procedures, Equation S4). Thus, we can quantify latency’s impact on I0 by measuring latency’s impact on viral loads in systemically infected patients.

To track latency’s effect on systemic viral loads, we simulated the deterministic model in the lymphoid compartment alone (i.e., R0 = 10). Initial mucosal infection was not tracked in these simulations because of the data showing that systemic infections emerge from single “founder” viruses independent of the inoculum (; ). These data indicate that mucosal dynamics affect the probability of systemic infection, but not the level once established. Thus, we assumed the existence of a single founder infected cell and solved Equation [6] numerically. Assuming successful systemic establishment, the systemic infection model shows that increasing plat decreases long-term viral loads (Figure 2B). Consequently, increasing the frequency of latency (plat) decreases infection inocula (I0) at the population scale.

The Evolutionarily Optimal Probability of Latency Is ~0.5

Given Equation [1], if latency’s benefit to pestab exceeds its cost to I0, then latency increases the probability of lentiviral transmission (ptransmission). Mathematically, this net evolutionary benefit of latency can only occur if the (evolutionarily optimal) value of plat that maximizes ptransmission is greater than 0. Here, we test whether the maximizing value of plat is greater than 0, directly quantifying latency’s net evolutionary benefit.

We first derive pestab as a function of plat. After initial mucosal infection, only latently infected cells persist, with the number of surviving latently infected cells defined to be LR0>1init. As noted above, due to R0muc < < 1, the majority of mucosal latent infections emerge in the first generation of infection, making LR0>1initplatI0 (Figures 2A, S2F, and S2G). At least one of these surviving infected cells must be reactivated (with probability preact) to establish systemic infection. Thus, the per-inoculum probability of establishing systemic infection is:

pestab=(LR0>1initI0)preactplatpreact
[2]

 

Equation [2] emerges from the result that only latently infected cells survive initial infection in the mucosa (Figures 2A and S2A–S2E). To demonstrate robustness, below we introduce a “leakage” probability (fnonlatent) that reflects the fraction of systemic infections that are established by non-latent cells—including Langerhans dendritic cells, actively infected cells, and free virions.

We next solve for I0 as a function of plat. As noted above, the average infectious dose (i.e., I0) that can be transmitted to a new individual is directly proportional to the time integral of the viral load— ∫ V(t)dt, Equation [S4]—over the duration of systemic infection (). Analytically solving this time integral yields (Extended Experimental Procedures, Section B):

I0const(plat)[(1plat)RLT01]
[3]

 

The constant term in Equation [3] only implies constant in plat—it may depend on other parameters. Further, Equation [3] is solved under the assumption that latently infected cells rarely reactivate prior to cell death (i.e., r < < dL in Table S1). This conservative assumption reduces the optimal level of latency by presuming that latently infected cells generally die before contributing to viral loads. Given this maximal fitness cost, latency reduces the reproductive ratio during systemic infection, R0LT, by the factor (1 − plat).

By combining Equations [13], ptransmission emerges as a function of plat (Figure 2C):

ptransmissionpestabI0const(plat)preactplat[(1plat)RLT01]
[4]

 

Equation [4] shows that, for each value of R0LT, the probability of viral transmission has an optimum at a specific plat. To analytically derive this optimum, we make the simplifying assumption that preact is constant in plat. This makes ptransmissionplat·[(1plat)RLT01]. Differentiating the simplified transmission probability with respect to plat yields the following optimal probability of latency, denoted poptlat:

poptlat1(1/RLT0)2
[5]

 

Strikingly, for a typical value of R0LT ~10 (), poptlat0.5 is the probability of latency that maximizes lentiviral transmission (Figure 2C).

In agreement with these analytic derivations, numerical solutions also show that ptransmission has an optimum at plat ≈ 0.5 (Figure 2D). The numerical simulations are generated by directly calculating ∫ V(t)dt in model runs, rather than approximating it via Equation [3]. Sensitivity analyses show that this optimum at plat ≈ 0:5 exists across the entire observed range of R0LT values (Figure 2D).

Large Optimal Latency Probability Is Robust to Changes in Model Assumptions

The main prediction of a large poptlat value remains valid even if one removes key mathematical assumptions. In particular, the two-compartment Wright-Fisher model (Extended Experimental Procedures, Section B) inverts the assumption that preact is constant in plat, allowing preact to strongly decrease in plat. Even in this extreme scenario—in which latency has a substantial fitness cost beyond its reduction of viral loads during systemic infection— poptlat>1/3 (Figure S2I). Similarly, the large poptlat value remains valid when one relaxes the assumption that only latently infected cells seed systemic infections. To show this, we analytically re-calculated poptlat when a fraction (fnonlatent) of successful infections are established via non-latent routes (Extended Experimental Procedures, Section E). Even if 80% of lentiviral transmissions are established via non-latent routes, poptlat=0.1. More generally, as long as fnonlatent is less than 100%, latency remains evolutionarily beneficial (Figures 2E and S2J).

Strikingly, relaxing other model assumptions increases the large poptlat value. For example, relaxing the assumption that latently infected cells die prior to reactivation (i.e., r < < dL) reduces the cost of latency during systemic infection and therefore increases the optimal latency probability. In fact, if r  dL, poptlat=1 (Extended Experimental Procedures, Section E). Further, if lentiviral transmissibility saturates at high viral loads ()—so that latency’s decrease of steady-state viral loads does not decrease I0—then poptlat would again equal 1, due to the absence of a cost to latency (Extended Experimental Procedures, Section E).

Simplified Two-Compartment Model Fits the High Frequencies of Latency Measured in Experimental Models

The predicted value of poptlat~0.5 matches the latency frequencies of 50% () or higher () measured in cell culture. poptlat~0.5 is also consistent with a recent in vivo study in Rhesus macaques, in which a large reservoir of latently infected cells is documented on day 3 of mucosal infection (). However, poptlat~0.5 is inconsistent with the low latency frequencies measured in chronically infected patients. Only 1 in 106–107 patient CD4+ T cells appear to be latently infected (; ). This has led to estimates of plat ~10−5 − 10−4 (; ). While more recent studies indicate that the latency frequency in patient cells is ~60-fold higher (), this still leaves plat < < 0.5 during chronic infection. Below, we show that the dichotomy between latency’s high frequency in early infection and cell culture and latency’s low frequency in chronic infection can be explained by the onset of the adaptive immune response.

Mathematical Models Incorporating the Immune Response Are Required to Explain the Divergent Latency Frequencies between Experimental Models and Patients

Unlike early mucosal infections or cell-culture infections, chronic lentiviral infections contain an HIV-specific adaptive immune response (). Previous work has shown that this adaptive immune response must be incorporated into the basic model of viral dynamics (; ) to fit the 2–3 log drop in viral loads between the viral peak during acute infection and the viral set point established during chronic infection (). We hypothesized that incorporating a canonical adaptive immune response (; ) would also be necessary to observe the reduced level of latently infected cells documented during chronic infection.

A substantial body of literature suggests that the model assumptions that plat and r are constant must be relaxed to account for the adaptive immune response. In particular, the activation levels of CD4+ T cells appear to increase during chronic infection in vivo, as is measured by the expression levels of three activation markers () and the increased turnover rates of CD4+ T cells (). While the exact mechanism is unknown, one potential driver of CD4+ T cell activation is the body’s homeostatic response to the depletion of CD4+ T cells during acute infection (). Another potential mechanism is CD8+ T cells’ secreting activating cytokines such as TNF-α (). Whatever the mechanism, cellular activation factors sharply decrease plat and sharply activate HIV transcription (; ; ), for example, by accumulating transcription factors (e.g., NF-κB) that activate the HIV LTR promoter. Further, in the companion study (), mathematical modeling shows that cellular activation levels bias HIV circuit output (i.e., plat and r), even though latency is hardwired into the circuit.

Since an adaptive immune response is associated with an increase in CD4+ T cell activation levels () that reduces plat and increases r (; ; ), we hypothesized that the adaptive-immune response could be responsible for the reduced plat levels in chronically infected patients (Figure 3A). This hypothesis was quantitatively tested by allowing plat and r to vary as functions of the effector CD8+ T cell concentration, E[t] (Extended Experimental Procedures, Section C). Before the initiation of the adaptive-immune response (i.e., before chronic infection), the model naturally generates high latency probabilities of ~0.5 and low reactivation rates, as in the simplified models above. However, after the viremia peak, cellular activation () and cell death () become substantial, increasing r(E[t]) to high levels and decreasing plat(E [t]) to low levels (Figure 3B). As a result, the immune model mechanistically explains the divergent latency frequencies measured between experimental models (cell culture and non-human primates) and chronically infected patients (Figure 3B).

Incorporating the Immune Response Explains the Divergent HIV-Latency Frequencies between Experimental Models and Patients

(A) Extended model of systemic HIV infection, which includes CD8+ T cells (E) that kill actively infected cells (or suppress viral replication) and activate latently infected cells (Equations [S9] and [S10]).

(B) The latency probability (plat) and reactivation rate (r) change dramatically around the time of the viremia peak due to the immune response (e.g., due to bystander cytokine activation by immune cells, Equation [S10]). Inclusion of immune cells into the model is capable of interpreting the low incidence of latently infected cells in chronically infected patients.

Models Incorporating the Immune Response Fit Available Patient Data while Retaining the Robust Optimal Latency Prediction

While the immune-response model interprets the low levels of plat measured during chronic infection, validation against all available patient data is a critical test of the model. Thus, wetested whether the model could recapitulate extant patient data on: (1) viral loads before ART (), (2) effector T cell concentrations before ART (), (3) latently infected cells before ART (), and (4) latently infected cells after ART (). Strikingly, the extended immune-response model is able to fit these four data plateaus (Figure 4A), using established parameter estimates (Table S2). In particular, the immune- response model reproduces the depressed latent reservoir of ~106 cells measured in chronically infected patients. Further, the model captures the ~1 log drop in the latent reservoir under ART (Figure 4A), because ART leads to antigen depletion. This causes the immune-cell population to contract and the reactivation rate r(t) to decrease to its low background level. To be sure that these fits were not artifacts due to model complexity, we also tested simplified immune response models (Extended Experimental Procedures, Section E). While these simplified models fit the four steady-state plateaus, they cannot reproduce the pre-steady-state kinetics measured in patients (Figure S3). In contrast, the full immune model fits both steady-state and pre-steady-state kinetics (Figure 4A, inset), including the viral decay kinetics measured in patients who undergo ART ().

The Extended Immune-Response Model Fits the Available In Vivo Data and Does Not Change the Optimal Latency Probability for Resting Cells, poptlat(0)

(A) Dynamics of cell compartments during systemic infection calculated from Equations [S9] and [S10]. Antiretroviral therapy (ART) initiated during steadystate infection causes a decline of the latent reservoir (L). The saturation of the fall in the latent reservoir is due to the decline in immune cells (E) during ART. (Data points across human patients) Virus load prior to ART () (green triangles); latent cells prior to ART () and after highly active ART () (cyan triangles); effector CD8 T cells () (red triangles). For each data set (triangles), box-and-whisker plots show the upper and lower quartiles of the patient data. (Blowout) Virus load after the onset of ART () (green triangles, error bars show SD).

(B) Normalized transmission rate ptransmission as a function of plat(0) calculated from the dynamics in A and Equation [1]. Two cases are shown for comparison: with immune cells (E, green triangles) and without immune cells (E = N = 0, blue curve). Inclusion of immune cells into the model only weakly affects the prediction of a large optimal latency probability for resting cells, poptlat(0)~0.5. Model parameters in A and B are in Tables S1 and S2 (with R0LT = 15 and plat(0) = 0.5 in A). See also Figure S3.

Critically, the level of the adaptive immune response does not change the prediction of the simplified model (i.e., the model without an immune response) that the initial latency probability plat(0) has a large optimum of ~0.5 (Figures 4B and S3). As a result, the prediction of the high optimal latency probability is directly applicable to natural lentiviral hosts even if they exhibit depressed immune responses. Further, as in the simplified models lacking an immune response, the large poptlat value is preserved even when a large fraction of systemic infections are mediated by non-latent cells (Extended Experimental Procedures, Section E). The optimal latency prediction is also robust to perturbations of epidemiological assumptions, such as the monotonic dependence of lentiviral transmission on viral loads (Extended Experimental Procedures, Section E). Overall, the robustness of poptlat in the immune model matches the robustness of poptlat in the simplified models.

Experimental Depletion of CD8+ T Cells in SIV-Infected Macaques Will Increase the Latent Reservoir ~3 Logs More Than Viremia

The immune model argues that CD8+ T cells depress the latent reservoir during chronic infection—either directly (e.g., through secreted cytokines) or indirectly (e.g., through activation of downstream cell types that secrete factors). Thus, a direct test of the model can be achieved by depleting CD8+ T cells with anti-CD8 antibodies. CD8 depletion should increase the latency probability (plat) toward its original high value of ~0.5 and concomitantly decrease the reactivation rate (r) toward its original low value. In fact, the model quantitatively predicts the outcome of this experiment. Whereas previous CD8 depletion studies have already measured an ~1–3 log increase in the number of actively infected cells following CD8 depletion in SIV-infected Rhesus macaques (; ; ), the model predicts that the latent reservoir will increase by ~5 logs following CD8 depletion (Figure 5A). Thus, the increase in the latent reservoir would be ~3 logs greater than the increase in actively infected cells and viremia (Figure 5B). A corollary prediction is that CD8 depletion during early pre-peak infection (), prior to a high-level adaptive immune response, will only increase the latent reservoir ~2- to 3-fold and will thus be harder to reliably measure (Figure S4). Notably, these experimental tests of the model require viral outgrowth assays () since directly measuring proviral DNA will only report on actively infected cells, which outnumber latently infected cells by orders of magnitude. A viral outgrowth assay post-CD8 depletion would provide quantitative verification of the model and would consequently test the model’s output that latency is a viral bet-hedging strategy tuned by natural selection.

Depletion of CD8+ T Cells in SIV-Infected Macaques Is Predicted to Increase the Latent Reservoir Significantly More Than Viremia

(A) Predicted dynamics in systemic infection for the extended model (Equations [S9] and [S10]). Data points and parameters are as in Figure 4, with the upper and lower quartiles of the patient data (triangles) shown in box-and-whisker plots.

(B) The ratio of virions to latently infected cells will be inverted following CD8+ T cell depletion (post-depletion corresponds to day 125 in A). The dramatic 2-log increase in viremia has been observed, as shown by the data points at 1 week post-depletion in  and . The dashed horizontal line at 10−3 RNA/ml/cell corresponds to a 1:1 ratio of latently and actively infected cell counts. Blue bars correspond to the parameters and compartment sizes in the simulation example in A. The maximal expected errors (vertical bars) are estimated from the whisker box borders in A (the two middle quartiles). Since the dynamic balance between actively infected cells and latently infected cells is modulated by plat and r, the depletion of immune cells affecting plat and r is predicted to change this balance and disproportionately increase the latent reservoir.

See also Figures S4 and S5.

Viral Strains Engineered to Have Higher Replicative Fitness—via Reduced Latency—Will Exhibit Lower Infectivity in Animal-Model Mucosal Inoculations

A more direct experimental test of the model would involve mucosal challenge experiments using recombinant SIV strains engineered to have substantially reduced latency probabilities. Engineering strains with reduced latency efficiencies appears possible since different HIV-1 clades are already known to exhibit different latency frequencies. These clade-specific differences appear to be driven by cis elements within the HIV-1 LTR (; ). The model directly predicts that the reduced-latency recombinants will establish self-propagating systemic infections less frequently than the wild-type strains maintaining high latency frequencies. Further, these reduced latency strains could be quantitatively tested for increased replicative fitnesses via competitive growth assays with wild-type strains. If decreasing latency both increased replicative fitness and decreased successful lentiviral transmission, this would directly show that proviral latency provides a bet-hedging advantage that increases viral transmission despite reducing steady-state viral loads.

Proviral Latency Contrasted with Alternate Mechanisms of Initial Viral Survival

A natural question is whether alternatives to latently infected CD4+ T cells exist that also increase the probability of initial viral survival in the mucosa. One proposed non-latent route is dendritic cell migration from the mucosa to the target-cell rich lymphoid tissue (; ). More specifically, Langerhans dendritic cells present in the mucosa can be infected by HIV and are prone to migration to the lymphoid tissue, where they can support subsequent dissemination of HIV by cis transfer (). Yet, Langerhans cells’ dissemination of HIV may be partially blocked by neutralizing antibodies (). Follicular dendritic cells may provide another route of viral survival; however, these cells do not migrate to the mucosa (). In contrast to dendritic cells, proviral latent cells are neither impacted by neutralizing antibodies (being quiescent) nor blocked by the mucosal barrier, which has been proposed to be a viral bottleneck (). Latency can thus act as a type of “Trojan horse” for the virus. More fundamentally, even if alternative routes of initial viral survival exist, the results of this study (i.e., poptlat>0) remain robust as long as latency seeds some fraction of systemic infections (Figures 2E and S2J).

Suppressing Latent Reactivation in the First Week of Infection Could Substantially Reduce the Latent Reservoir, Enhancing “Kick-and-Kill” Therapy

The model presents a potential therapeutic strategy that exploits the need for latently infected cells to reactivate to both establish systemic infection and dramatically increase the size of the latent reservoir (Figure S5). Thus, if the early reactivation rate were reduced—for example, by suppressing antigen- presenting cell (APC) migration () or HIV transcriptional reactivation ()—systemic infection would be rendered less likely and the latent reservoir size would be substantially decreased (Figure S5). While a caveat of this proposed approach is detection and treatment within the first week of infection, similar early treatments have been achieved; for a review, see . Critically, a substantially smaller latent reservoir of ~102 cells would require the reactivation of far fewer latent cells by imperfect “shock-and-kill” strategies (; ). As a result, suppression of reactivation during the first week of infection followed by shock and kill could substantially enhance the chances of HIV eradication.

Implications for Alternate Antiviral Therapy Approaches

A further implication of the result that latency is a hardwired, evolutionarily maintained trait is that it may be easier to control HIV by increasing, rather than purging, the latent reservoir (; ; ). Current shock-and-kill therapies are fighting natural selection in attempting to reactivate each of ~105 latent cells. In contrast, discovering a non-toxic compound that switches 90%–95% of actively infected cells to latency would drive HIV’s basic reproductive ratio (R0) below 1, making HIV infection unsustainable. While still a hypothetical avenue, enhancing viral latency may provide a viable alternative if shock-and-kill strategies fail to achieve their goal of complete eradication.

EXPERIMENTAL PROCEDURES

A Simplified Two-Compartment Model to Quantify the Net Impact of Latency on Lentiviral Transmission

All models described in the main text are variations of the well-parameterized basic model of viral dynamics () expanded to include latent infections (, ; , ). Absent an immune response, the deterministic form of the models is captured by the following ordinary differential equations:

UninfectedtargetcellsdTdt=breplenishmentdTTnatural deathkVTinfection

 

Actively infected cellsdIdt=(1plat)kVTactive infectiondIIdeath+rLreactivation
[6]

 

Latently infected cellsdLdt=platkVTlatent infectiondLLdeathrLreactivation

 

VirusdVdt=ndIIproductioncVclearance

 

In the model above, uninfected “target” cells (T) are produced at rate b, decay at rate dT, and can be infected by virus particles (V) at rate k. Upon viral infection, target cells become either latently infected cells (L) with probability plat or become actively infected (virus-producing) cells (I) with probability 1 − plat. Latently infected cells reactivate into actively infected cells at rate r or die at the (slow) rate dL. Actively infected cells produce “burst sizes” of n virions as they die at rate dI. Virions decay at the relatively fast rate c. All parameter values are given in Table S1; Table S2 contains parameters for the model extended to include an adaptive immune response (Extended Experimental Procedures, Section C).

Critically, the infection models can be simplified by re-parameterizing the equations interms of thebasic reproductiveratio:R0=bkn/cdT.This “non-dimensionalization” enables us to capture the disparate dynamics between mucosal infection (Figure 2A) and systemic infection (Figure 2B) by simulating the same model for both infection stages and only varying a single parameter, R0. Further, R0muc is experimentally bounded to be < < 1 from the viral dynamics during initial infection (), and R0LTis similarly measured to be ~10 during systemic infection (). As a result, no assumptions about unknown parameter values are needed to obtain the optimal latency probability (poptlat). More directly, Equation [5] shows that (poptlat) only depends on R0LT (for detailed derivations and tests of the models, see Extended Experimental Procedures).

Highlights

  • Mathematical model proposes evolutionary basis for HIV latency

  • Hardwired latency circuit enhances HIV transmission across target-cell-poor mucosa

  • Predicted optimal latency rate for HIV transmission matches measured levels

  • Model predictions are testable in primates by modulating latency rates or CD8 levels

Supplementary Material

01

Click here to view.(736K, pdf)

ACKNOWLEDGMENTS

Weare grateful to Lani Wu and Stephen Altschuler for input and discussions; to Alan Perelson, Warner Greene, Eric Verdin, Anand Pai, Brandon Razooky, and John Coffin for helpful comments; and to Abhyudai Singh for generating preliminary data. We are also grateful to the graphics department at the Gladstone Institutes for artistic expertise and help with figure schematics. This work was supported by the Alfred P. Sloan Foundation, the Wyss Institute Technology Development Fellowship, the NIH Director’s Pioneer Award Program (DP1 OD017181), as well as NIH awards R21AI109611, F32AI102520, and U19AI096113 as part of the Delaney Collaboratory for AIDS Research and Eradication (CARE).

Footnotes

 

SUPPLEMENTAL INFORMATION

Supplemental Information includes Extended Experimental Procedures, five figures, and two tables and can be found with this article online at http://dx.doi.org/10.1016/j.cell.2015.02.017.

 

REFERENCES

  • Archin NM, Liberty AL, Kashuba AD, Choudhary SK, Kuruc JD, Crooks AM, Parker DC, Anderson EM, Kearney MF, Strain MC, et al. Administration of vorinostat disrupts HIV-1 latency in patients on antiretroviral therapy. Nature. 2012;487:482–485. [PMC free article] [PubMed] []
  • Arkin A, Ross J, McAdams HH. Stochastic kinetic analysis of developmental pathway bifurcation in phage lambda-infected Escherichia coli cells. Genetics. 1998;149:1633–1648. [PMC free article] [PubMed] []
  • Balaban NQ. Persistence: mechanisms for triggering and enhancing phenotypic variability. Curr. Opin. Genet. Dev. 2011;21:768–775. [PubMed] []
  • Batada NN, Hurst LD. Evolution of chromosome organization driven by selection for reduced gene expression noise. Nat. Genet. 2007;39:945–949. [PubMed] []
  • Burnett JC, Miller-Jensen K, Shah PS, Arkin AP, Schaffer DV. Control of stochastic gene expression by host factors at the HIV promoter. PLoS Pathog. 2009;5:e1000260. [PMC free article] [PubMed] []
  • Calvanese V, Chavez L, Laurent T, Ding S, Verdin E. Dual-color HIV reporters trace a population of latently infected cells and enable their purification. Virology. 2013;446:283–292. [PMC free article] [PubMed] []
  • Chun T-W, Stuyver L, Mizell SB, Ehler LA, Mican JAM, Baseler M, Lloyd AL, Nowak MA, Fauci AS. Presence of an inducible HIV-1 latent reservoir during highly active antiretroviral therapy. Proc. Natl. Acad. Sci. USA. 1997a;94:13193–13197. [PMC free article] [PubMed] []
  • Chun TW, 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. 1997b;387:183–188. [PubMed] []
  • Chun TW, Engel D, Mizell SB, Ehler LA, Fauci AS. Induction of HIV-1 replication in latently infected CD4+ T cells using a combination of cytokines. J. Exp. Med. 1998;188:83–91. [PMC free article] [PubMed] []
  • Coffin J, Swanstrom R. HIV pathogenesis: dynamics and genetics of viral populations and infected cells. Cold Spring Harb. Perspect. Med. 2013;3:a012526. [PMC free article] [PubMed] []
  • Cohen D. Optimizing reproduction in a randomly varying environment. J. Theor. Biol. 1966;12:119–129. [PubMed] []
  • Dahabieh MS, Ooms M, Simon V, Sadowski I. A doubly fluorescent HIV-1 reporter shows that the majority of integrated HIV-1 is latent shortly after infection. J. Virol. 2013;87:4716–4727. [PMC free article] [PubMed] []
  • Dar RD, Hosmane NN, Arkin MR, Siliciano RF, Weinberger LS. Screening for noise in gene expression identifies drug synergies. Science. 2014;344:1392–1396. [PMC free article] [PubMed] []
  • De Boer RJ, Perelson AS. Target cell limited and immune control models of HIV infection: a comparison. J. Theor. Biol. 1998;190:201–214. [PubMed] []
  • Deeks SG. HIV: Shock and kill. Nature. 2012;487:439–440. [PubMed] []
  • Doitsh G, Cavrois M, Lassen KG, Zepeda O, Yang Z, Santiago ML, Hebbeler AM, Greene WC. Abortive HIV infection mediates CD4 T cell depletion and inflammation in human lymphoid tissue. Cell. 2010;143:789–801. [PMC free article] [PubMed] []
  • Eisele E, Siliciano RF. Redefining the viral reservoirs that prevent HIV-1 eradication. Immunity. 2012;37:377–388. [PMC free article] [PubMed] []
  • Finzi D, Hermankova M, Pierson T, Carruth LM, Buck C, Chaisson RE, Quinn TC, Chadwick K, Margolick J, Brookmeyer R, et al. Identification of a reservoir for HIV-1 in patients on highly active antiretroviral therapy. Science. 1997;278:1295–1300. [PubMed] []
  • Finzi D, Blankson J, Siliciano JD, Margolick JB, Chadwick K, Pierson T, Smith K, Lisziewicz J, Lori F, Flexner C, et al. Latent infection of CD4+ T cells provides a mechanism for lifelong persistence of HIV-1, even in patients on effective combination therapy. Nat. Med. 1999;5:512–517. [PubMed] []
  • Fraser HB, Hirsh AE, Giaever G, Kumm J, Eisen MB. Noise minimization in eukaryotic gene expression. PLoS Biol. 2004;2:e137. [PMC free article] [PubMed] []
  • Fraser C, Hollingsworth TD, Chapman R, de Wolf F, Hanage WP. Variation in HIV-1 set-point viral load: epidemiological analysis and an evolutionary hypothesis. Proc. Natl. Acad. Sci. USA. 2007;104:17441–17446. [PMC free article] [PubMed] []
  • Gray RH, Wawer MJ, Brookmeyer R, Sewankambo NK, Serwadda D, Wabwire-Mangen F, Lutalo T, Li X, vanCott T, Quinn TC Rakai Project Team. Probability of HIV-1 transmission per coital act in monogamous, heterosexual, HIV-1-discordant couples in Rakai, Uganda. Lancet. 2001;357:1149–1153. [PubMed] []
  • Grimmett G, Stirzaker D. Probability and Random Processes. Second Edition. Oxford University Press; 1992. []
  • Haaland RE, Hawkins PA, Salazar-Gonzalez J, Johnson A, Tichacek A, Karita E, Manigart O, Mulenga J, Keele BF, Shaw GM, et al. Inflammatory genital infections mitigate a severe genetic bottleneck in heterosexual transmission of subtype A and C HIV-1. PLoS Pathog. 2009;5:e1000274. [PMC free article] [PubMed] []
  • Haase AT. Early events in sexual transmission of HIV and SIV and opportunities for interventions. Annu. Rev. Med. 2011;62:127–139. [PubMed] []
  • Han Y, Wind-Rotolo M, Yang HC, Siliciano JD, Siliciano RF. Experimental approaches to the study of HIV-1 latency. Nat. Rev. Microbiol. 2007;5:95–106. [PubMed] []
  • Hartl DL, Clark AG. Principles of population genetics. Fourth Edition. Sunderland, Mass: Sinauer Associates; 2007. []
  • Ho YC, Shan L, Hosmane NN, Wang J, Laskey SB, Rosenbloom DI, Lai J, Blankson JN, Siliciano JD, Siliciano RF. Replication-competent noninduced proviruses in the latent reservoir increase barrier to HIV-1 cure. Cell. 2013;155:540–551. [PMC free article] [PubMed] []
  • Jeeninga RE, Westerhout EM, van Gerven ML, Berkhout B. HIV-1 latency in actively dividing human T cell lines. Retrovirology. 2008;5:37. [PMC free article] [PubMed] []
  • Jin X, Bauer DE, Tuttleton SE, Lewin S, Gettie A, Blanchard J, Irwin CE, Safrit JT, Mittler J, Weinberger L, et al. Dramatic rise in plasma viremia after CD8(+) T cell depletion in simian immunodeficiency virus-infected macaques. J. Exp. Med. 1999;189:991–998. [PMC free article] [PubMed] []
  • Kahn JO, Walker BD. Acute human immunodeficiency virus type 1 infection. N. Engl. J. Med. 1998;339:33–39. [PubMed] []
  • Kearney M, Maldarelli F, Shao W, Margolick JB, Daar ES, Mellors JW, Rao V, Coffin JM, Palmer S. Human immunodeficiency virus type 1 population genetics and adaptation in newly infected individuals. J. Virol. 2009;83:2715–2727. [PMC free article] [PubMed] []
  • Keele BF, Giorgi EE, Salazar-Gonzalez JF, Decker JM, Pham KT, Salazar MG, Sun C, Grayson T, Wang S, Li H, et al. Identification and characterization of transmitted and early founder virus envelopes in primary HIV-1 infection. Proc. Natl. Acad. Sci. USA. 2008;105:7552–7557. [PMC free article] [PubMed] []
  • Kuroda MJ, Schmitz JE, Charini WA, Nickerson CE, Lifton MA, Lord CI, Forman MA, Letvin NL. Emergence of CTL coincides with clearance of virus during primary simian immunodeficiency virus infection in rhesus monkeys. J. Immunol. 1999;162:5127–5133. [PubMed] []
  • Li Q, Duan L, Estes JD, Ma ZM, Rourke T, Wang Y, Reilly C, Carlis J, Miller CJ, Haase AT. Peak SIV replication in resting memory CD4+ T cells depletes gut lamina propria CD4+ T cells. Nature. 2005;434:1148–1152. [PubMed] []
  • Markowitz M, Louie M, Hurley A, Sun E, Di Mascio M, Perelson AS, Ho DD. A novel antiviral intervention results in more accurate assessment of human immunodeficiency virus type 1 replication dynamics and T-cell decay in vivo. J. Virol. 2003;77:5037–5038. [PMC free article] [PubMed] []
  • Matano T, Shibata R, Siemon C, Connors M, Lane HC, Martin MA. Administration of an anti-CD8 monoclonal antibody interferes with the clearance of chimeric simian/human immunodeficiency virus during primary infections of rhesus macaques. J. Virol. 1998;72:164–169. [PMC free article] [PubMed] []
  • Metzner KJ, Jin X, Lee FV, Gettie A, Bauer DE, Di Mascio M, Perelson AS, Marx PA, Ho DD, Kostrikis LG, Connor RI. Effects of in vivo CD8(+) T cell depletion on virus replication in rhesus macaques immunized with a live, attenuated simian immunodeficiency virus vaccine. J. Exp. Med. 2000;191:1921–1931. [PMC free article] [PubMed] []
  • Miller CJ, Li Q, Abel K, Kim EY, Ma ZM, Wietgrefe S, La Franco-Scheuch L, Compton L, Duan L, Shore MD, et al. Propagation and dissemination of infection after vaginal transmission of simian immunodeficiency virus. J. Virol. 2005;79:9217–9227. [PMC free article] [PubMed] []
  • Mohri H, Bonhoeffer S, Monard S, Perelson AS, Ho DD. Rapid turnover of T lymphocytes in SIV-infected rhesus macaques. Science. 1998;279:1223–1227. [PubMed] []
  • Murphy K. Janeway’s Immunobiology. Eighth Edition. London, New York: Garland Science; 2011. []
  • Nowak M, May R. Virus dynamics: mathematical principles of immunology and virology. Oxford, New York: Oxford University Press; 2000. []
  • Pearson JE, Krapivsky P, Perelson AS. Stochastic theory of early viral infection: continuous versus burst production of virions. PLoS Comput. Biol. 2011;7:e1001058. [PMC free article] [PubMed] []
  • Peressin M, Proust A, Schmidt S, Su B, Lambotin M, Biedma ME, Laumond G, Decoville T, Holl V, Moog C. Efficient transfer of HIV-1 in trans and in cis from Langerhans dendritic cells and macrophages to autologous T lymphocytes. AIDS. 2014;28:667–677. [PubMed] []
  • Pierson T, McArthur J, Siliciano RF. Reservoirs for HIV-1: mechanisms for viral persistence in the presence of antiviral immune responses and antiretroviral therapy. Annu. Rev. Immunol. 2000;18:665–708. [PubMed] []
  • Razooky BS, Pai A, Aull K, Rouzine IM, Weinberger LS. A hardwired HIV latency program. Cell. 2015;160:990–1001. this issue. [PMC free article] [PubMed] []
  • Richman DD, Margolis DM, Delaney M, Greene WC, Hazuda D, Pomerantz RJ. The challenge of finding a cure for HIV infection. Science. 2009;323:1304–1307. [PubMed] []
  • Rong L, Perelson AS. Modeling HIV persistence, the latent reservoir, and viral blips. J. Theor. Biol. 2009a;260:308–331. [PMC free article] [PubMed] []
  • Rong L, Perelson AS. Modeling latently infected cell activation: viral and latent reservoir persistence, and viral blips in HIV-infected patients on potent therapy. PLoS Comput. Biol. 2009b;5:e1000533. [PMC free article] [PubMed] []
  • Santiago ML, Range F, Keele BF, Li Y, Bailes E, Bibollet-Ruche F, Fruteau C, Noë R, Peeters M, Brookfield JF, et al. Simian immunodeficiency virus infection in free-ranging sooty mangabeys (Cercocebus atys atys) from the Taï Forest, Côte d’Ivoire: implications for the origin of epidemic human immunodeficiency virus type 2. J. Virol. 2005;79:12515–12527. [PMC free article] [PubMed] []
  • Schmitz JE, Kuroda MJ, Santra S, Sasseville VG, Simon MA, Lifton MA, Racz P, Tenner-Racz K, Dalesandro M, Scallon BJ, et al. Control of viremia in simian immunodeficiency virus infection by CD8+ lymphocytes. Science. 1999;283:857–860. [PubMed] []
  • Sedaghat AR, Siliciano JD, Brennan TP, Wilke CO, Siliciano RF. Limits on replenishment of the resting CD4+ T cell reservoir for HIV in patients on HAART. PLoS Pathog. 2007;3:e122. [PMC free article] [PubMed] []
  • Sedaghat AR, Siliciano RF, Wilke CO. Low-level HIV-1 replication and the dynamics of the resting CD4+ T cell reservoir for HIV-1 in the setting of HAART. BMC Infect. Dis. 2008;8:2. [PMC free article] [PubMed] []
  • Siliciano RF, Greene WC. HIV latency. Cold Spring Harbor Perspect. Med. 2011;1:a007096. [PMC free article] [PubMed] []
  • Stafford MA, Corey L, Cao Y, Daar ES, Ho DD, Perelson AS. Modeling plasma virus concentration during primary HIV infection. J. Theor. Biol. 2000;203:285–301. [PubMed] []
  • Su B, Xu K, Lederle A, Peressin M, Biedma ME, Laumond G, Schmidt S, Decoville T, Proust A, Lambotin M, et al. Neutralizing antibodies inhibit HIV-1 transfer from primary dendritic cells to autologous CD4 T lymphocytes. Blood. 2012;120:3708–3717. [PubMed] []
  • Turnbull EL, Wong M, Wang S, Wei X, Jones NA, Conrod KE, Aldam D, Turner J, Pellegrino P, Keele BF, et al. Kinetics of expansion of epitope-specific T cell responses during primary HIV-1 infection. J. Immunol. 2009;182:7131–7145. [PubMed] []
  • van der Sluis RM, Pollakis G, van Gerven ML, Berkhout B, Jeeninga RE. Latency profiles of full length HIV-1 molecular clone variants with a subtype specific promoter. Retrovirology. 2011;8:73. [PMC free article] [PubMed] []
  • Wawer MJ, Gray RH, Sewankambo NK, Serwadda D, Li X, Laeyendecker O, Kiwanuka N, Kigozi G, Kiddugavu M, Lutalo T, et al. Rates of HIV-1 transmission per coital act, by stage of HIV-1 infection, in Rakai, Uganda. J. Infect. Dis. 2005;191:1403–1409. [PubMed] []
  • Weinberger AD, Weinberger LS. Stochastic fate selection in HIV-infected patients. Cell. 2013;155:497–499. [PubMed] []
  • Weinberger LS, Burnett JC, Toettcher JE, Arkin AP, Schaffer DV. Stochastic gene expression in a lentiviral positive-feedback loop: HIV-1 Tat fluctuations drive phenotypic diversity. Cell. 2005;122:169–182. [PubMed] []
  • Weinberger LS, Dar RD, Simpson ML. Transient-mediated fate determination in a transcriptional circuit of HIV. Nat. Genet. 2008;40:466–470. [PubMed] []
  • Whitney JB, Hill AL, Sanisetty S, Penaloza-MacMaster P, Liu J, Shetty M, Parenteau L, Cabral C, Shields J, Blackmore S, et al. Rapid seeding of the viral reservoir prior to SIV viraemia in rhesus monkeys. Nature. 2014;512:74–77. [PMC free article] [PubMed] []
  • Wu L, KewalRamani VN. Dendritic-cell interactions with HIV: infection and viral dissemination. Nat. Rev. Immunol. 2006;6:859–868. [PMC free article] [PubMed] []
  • Zhang Z, Schuler T, Zupancic M, Wietgrefe S, Staskus KA, Reimann KA, Reinhart TA, Rogan M, Cavert W, Miller CJ, et al. Sexual transmission and propagation of SIV and HIV in resting and activated CD4+ T cells. Science. 1999;286:1353–1357. [PubMed] []

 

 

 

 

 

Jawahar Raina

Leave a comment

Please note: comments must be approved before they are published.