Incomplete inhibition of HIV infection results in more HIV infected lymph node cells by reducing cell death
Abstract
HIV has been reported to be cytotoxic in vitro and in lymph node infection models. Using a computational approach, we found that partial inhibition of transmissions of multiple virions per cell could lead to increased numbers of live infected cells. If the number of viral DNA copies remains above one after inhibition, then eliminating the surplus viral copies reduces cell death. Using a cell line, we observed increased numbers of live infected cells when infection was partially inhibited with the antiretroviral efavirenz or neutralizing antibody. We then used efavirenz at concentrations reported in lymph nodes to inhibit lymph node infection by partially resistant HIV mutants. We observed more live infected lymph node cells, but with fewer HIV DNA copies per cell, relative to no drug. Hence, counterintuitively, limited attenuation of HIV transmission per cell may increase live infected cell numbers in environments where the force of infection is high.
https://doi.org/10.7554/eLife.30134.001eLife digest
The HIVvirus infects cells of the immune system. Once inside, it hijacks the cellular molecular machineries to make more copies of itself, which are then transmitted to new host cells. HIV eventually kills most cells it infects, either in the steps leading to the infection of the cell, or after the cell is already producing virus. HIV can spread between cells in two ways, known as cellfree or celltocell. In the first, individual viruses are released from infected cells and move randomly through the body in the hope of finding new cells to infect. In the second, infected cells interact directly with uninfected cells. The second method is often much more successful at infecting new cells since they are exposed to multiple virus particles.
HIV infections can be controlled by using combinations of antiretroviral drugs, such as efavirenz, to prevent the virus from making more of itself. With a high enough dose, the drugs can in theory completely stop HIV infections, unless the virus becomes resistant to treatment. However, some patients continue to use these drugs even after the virus they are infected with develops resistance. It is not clear what effect taking ineffective, or partially effective, drugs has on how HIV progresses.
Using efavirenz, Jackson, Hunter et al. partially limited the spread of HIV between human cells grown in the laboratory. The experiments mirrored the situation where a partially resistant HIV strain spreads through the body. The results show that the success of cellfree infection is reduced as drug dose increases. Yet paradoxically, in celltocell infection, the presence of drug caused more cells to become infected. This can be explained by the fact that, in celltocell spread, each cell is exposed to multiple copies of the virus. The drug dose reduced the number of viral copies per cell without stopping the virus from infecting completely. The reduced number of viral copies per cell made it more likely that infected cells would survive the infection long enough to produce virus particles themselves.
Viruses that can kill cells, such as HIV, must balance the need to make more of themselves against the speed that they kill their host cell to maximize the number of infected cells. If transmission between cells is too effective and too many virus particles are delivered to the new cell, the virus may not manage to infect new hosts before killing the old ones. These findings highlight this delicate balance. They also indicate a potential issue in using drugs to treat partially resistant virus strains. Without care, these treatments could increase the number of infected cells in the body, potentially worsening the effects of living with HIV.
https://doi.org/10.7554/eLife.30134.002Introduction
HIV infection is known to result in extensive T cell depletion in lymph node environments (Sanchez et al., 2015), where infection is most robust (Brenchley et al., 2004; Doitsh et al., 2010; Doitsh et al., 2014; Finkel et al., 1995; Galloway et al., 2015; Mattapallil et al., 2005). Depletion of HIV infectable target cells, in addition to onset of immune control, is thought to account for the decreased replication ratio of HIV from an initial peak in early infection (Bonhoeffer et al., 1997; Nowak and May, 2000; Perelson, 2002; Phillips, 1996; QuiñonesMateu and Arts, 2006; Ribeiro et al., 2010; Wodarz and Levy, 2007). This is consistent with observations that individuals are most infectious in the initial, acute stage of infection, where the target cell population is relatively intact and produces high viral loads (Hollingsworth et al., 2008; Wawer et al., 2005).
Tcell death occurs by several mechanisms, which are either directly or indirectly mediated by HIV infection. Accumulation of incompletely reverse transcribed HIV transcripts is sensed by interferonγ–inducible protein 16 (Monroe et al., 2014) and leads to pyroptotic death of incompletely infected cells by initiating a cellular defence program involving the activation of caspase 1 (Doitsh et al., 2010; Doitsh et al., 2014; Galloway et al., 2015). HIV proteins Tat and Env have also been shown to lead to cell death of infected cells through CD95mediated apoptosis following Tcell activation (Banda et al., 1992; Westendorp et al., 1995a; Westendorp et al., 1995b1995). Using SIV infection, it has been shown that damage to lymph nodes due to chronic immune activation leads to an environment less conducive to Tcell survival (Zeng et al., 2012). Finally, double strand breaks in the host DNA caused by integration of the reverse transcribed virus results in cell death by the DNAPKmediated activation of the p53 response (Cooper et al., 2013).
The lymph node environment is conducive to HIV infection due to: (1) presence of infectable cells (Deleage et al., 2016; Embretson et al., 1993; TennerRacz et al., 1998); (2) proximity of cells to each other and lack of flow which should enable celltocell HIV spread (Baxter et al., 2014; Dale et al., 2011; Groot et al., 2008; Groppelli et al., 2015; Gummuluru et al., 2002; Hübner et al., 2009; Jolly et al., 2004; Jolly et al., 2011; Münch et al., 2007; Sherer et al., 2007; Sourisseau et al., 2007; Sowinski et al., 2008); (3) decreased penetration of antiretroviral therapy (ART) (Fletcher et al., 2014a). Multiple infections per cell have been reported in celltocell spread of HIV (Baxter et al., 2014; Boullé et al., 2016; Dang et al., 2004; Del Portillo et al., 2011; Dixit and Perelson, 2004; Duncan et al., 2013; Law et al., 2016; Reh et al., 2015; Russell et al., 2013; Sigal et al., 2011; Zhong et al., 2013). In this mode of HIV transmission, an interaction between the infected donor cell and the uninfected target results in directed transmission of large numbers of virions (Baxter et al., 2014; Groppelli et al., 2015; Hübner et al., 2009; Sowinski et al., 2008). This is in contrast to cellfree infection, where freefloating virus finds target cells through diffusion. Both modes occur simultaneously when infected donor cells are cocultured with targets. However, the celltocell route is thought to be the main cause of multiple infections per cell (Hübner et al., 2009). In the lymph nodes, several studies showed multiple infections (Gratton et al., 2000; Jung et al., 2002; Law et al., 2016) while another study did not (Josefsson et al., 2013). One explanation for the divergent results is that different cell subsets are infected to different degrees. For example, T cells were shown not to be multiply infected in the peripheral blood compartment (Josefsson et al., 2011). However, more recent work investigating markers associated with HIV latency in the face of ART found that the average number of HIV DNA copies per cell is greater than one in 3 out of 12 individuals. This occurred in the face of ART in the CD3positive, CD32a high CD4 Tcell subset (Descours et al., 2017). In the absence of suppressive ART, it would be expected that the number of HIV DNA copies per cell would be higher.
Multiple viral integration attempts per cell may increase the probability of death. One consequence of HIVmediated death may be that attenuation of infection may increase viral replication by increasing the number of live targets. Indeed, it has been suggested that more attenuated HIV strains result in more successful infections in terms of the ability of the virus to replicate in the infected individual (Ariën et al., 2005; Nowak and May, 2000; Payne et al., 2014; QuiñonesMateu and Arts, 2006; Wodarz and Levy, 2007).
Here, we experimentally examined the effect of attenuating celltocell spread by using HIV inhibitors. We observed that partially inhibiting infection with drug or antibody resulted in an increase in the number of live infected cells in both a cell line and in lymph node cells. This is, to our knowledge, the first experimental demonstration at the cellular level that attenuation of HIV infection can result in an increase in live infected cells under specific infection conditions.
Results
We introduce a model of infection where each donor to target transmission leads to an infection probability r and death probability q per infection attempt. In our experimental system, one infection attempt is measured as one HIV DNA copy, whether integrated or unintegrated. The probability of successful infection of a target cell given n infection attempts is 1(1 r)^{n} (Sigal et al., 2011). We define L_{n} as the probability of a cell to survive infection in the face of n infection attempts. Assuming infection attempts act independently, L_{n}=(1q)^{n}. The probability of a cell to be infected and not die after it has been exposed to n infection attempts is therefore:
This model makes several simplifying assumptions: (1) all infection attempts have equal probabilities to infect targets. (2) The probability for a cell to die from each transmission is equal between transmissions. (3) Infection attempts act independently, and productive infection and death are independent events. In this model, r and q capture the probabilities for a cell to be infected or die postreverse transcription. For example, mutations which reduce viral fitness by decreasing the probability of HIV to integrate would reduce r, while mutations which reduce the probability of successful reverse transcription would reduce n.
If the number of infection attempts n is Poisson distributed with mean λ, the probability for a cell to be infected is 1e^{rλ} and the probability of a cell to live is L_{n} = e^{qλ} (see Supplementary file 1 for parameters and definitions). As derived in Appendix 1, the probability that a cell is productively infected will be:
Since antiretroviral drugs lead to a reduction in the number of infection attempts by, for example, decreasing the probability of reverse transcription in the case of reverse transcriptase inhibitors, we introduced a drug strength value d, where d = 1 in the absence of drug and d > 1 in the presence of drug. In the presence of drug, λ is decreased to λ/d. The drug therefore tunes λ, and if the antiretroviral regimen is fully suppressive, λ/d is expected to be below what is required for ongoing replication. The probability of a cell to be infected and live given drug strength d is therefore:
Analysis of the probability of a cell to survive and be infected as a function of r and q shows that at each drug strength d/λ, P_{λ} increases as the probability of infection r increases (Figure 1). Hence, the value of r strongly influences the amplitude of P_{λ}. How P_{λ} behaves when drug strength d/λ increases depends on the parameter values of r and q. A subset of parameter values results in a peak in the number of infected cells at intermediate d/λ, decreasing as drug strength increases further (Figure 1). We refer to such a peak in infected numbers as an infection optimum. As q increases, the cost of multiple infections per cell increases, and the infection optimum shifts to higher d/λ values. A fall from the infection optimum at decreasing d/λ is driven by increasing cell death as a result of increasing infection attempts per cell. This slope is therefore shallower, and peaks broader, at low q values (Figure 1).
Our model assumes that cellular infection and death due to an HIV infection attempt are independent processes. This is based on observations that support a role for cell death as a cellular defence mechanism which may occur before productive infection, such as programmed cell death triggered by HIV integration induced DNA damage (Cooper et al., 2013). An alternative model is that HIVmediated cell death depends on productive infection. This would be consistent with cell death due to, for example, expression of viral proteins (Westendorp et al., 1995b1995). Since the concentration of viral proteins may also scale with the number of infections per cell, we derived the mathematical model for such a process in the supplementary mathematical analysis. The models are equivalent, showing that independence of cell death and infection is not a necessary condition for an infection optimum to occur in the presence of inhibitor.
Given that an infection optimum is dependent on parameter values, we next examined whether these parameter values occur experimentally in HIV infection. We therefore first tested for an infection optimum in the RevCEM cell line engineered to express GFP upon HIV Rev protein expression (Wu et al., 2007). We subcloned the cell line to maximize the frequency of GFPpositive cells upon infection (Boullé et al., 2016). We needed to detect the number of infection attempts per cell λ. To estimate this, we used PCR to detect the number of reverse transcribed copies of viral DNA in the cell by splitting each individual infected cell over multiple wells. We then detected the number of wells with HIV DNA by PCR amplification of the reverse transcriptase gene. Hence, the number of positive wells indicated the minimum number of viral DNA copies per cell, since more than one copy can be contained within the same well (Josefsson et al., 2011; Josefsson et al., 2013). We first measured the number of viral DNA copies in ACH2 cells, reported to contain a single inactive HIV integration per genome (Chun et al., 1997; O'Doherty et al., 2002). We sorted a total of 166 ACH2 cells at one cell per well into lysis buffer and subdivided singlecell lysates into four wells (Figure 2—figure supplement 1A). About one quarter of cells showed a PCR product of the expected size. Cells with more than one HIV copy per cell were very rare and may reflect either errors in cell sorting or dividing cells (Figure 2—figure supplement 1B). Similar frequencies were obtained when the ACH2 cell line was subcloned or split over 10 wells (Figure 2—figure supplement 1C). Given that each ACH2 cell contains one HIV DNA copy, the frequency of detection indicated our detection efficiency per HIV DNA copy.
To investigate the effect of multiple infection attempts per cell, we used coculture infection, where infected (donor) cells are coincubated with uninfected (target) cells and lead to celltocell spread. We used approximately 2% infected donor cells as our input, and detected the number of HIV DNA copies per cell by flow cytometric sorting of individual GFPpositive cells followed by splitting each cell lysate over 10 wells. Wells were then amplified by PCR and visualized on an agarose gel (Figure 2A). We assayed 60 cells and obtained a wide distribution of viral DNA copies per cell, which ranged from 0 to 9 copies (Figure 2B). The range of HIV DNA copies per cell fit a Poisson distribution with two means better than either a single mean Gaussian or Poisson distribution. However, the fit of the two mean Poisson distribution did not show two obvious peaks, and instead seemed to fit the data better due to the addition of fit parameters (Figure 2—figure supplement 2). Hence we cannot conclude that the distribution is bimodal. We also detected the HIV copy number in 30 GFPpositive cells infected by cellfree HIV. HIV in cellfree form was obtained by filtering supernatant from HIV producing cells to exclude cells or cell fragments, then infecting target cells with the filtered virus. Infection with this virus is defined here as cellfree infection. In this case, we detected either zero or one HIV copy per cell (Figure 2B inset). The frequency of single HIV DNA copies was 0.23, identical to the measured result in the ACH2 cell line. We computationally corrected the detected number of DNA copies in coculture infection for the sensitivity of our PCR reaction as determined by the ACH2 results (Materials and methods). On average we obtained 15 ± 3 copies per cell after correction.
To tune λ, we added the HIV reverse transcriptase inhibitor efavirenz (EFV) to infections. To calculate d, we used cellfree infection (Figure 2C, see Figure 2—figure supplement 3 for logarithmic yaxis plot), which as verified above, results in single HIV copies per cell. For cellfree infection, we approximate d = 1/Tx, where Tx is defined as the number of infected cells with drug divided by the number of infected cells without drug with single infection attempts (see Materials and methods and [(Sigal et al., 2011]). This is equivalent to 1ε in a commonly used model describing the effect of inhibitors on infection. In this model, ε is drug effectiveness, with the 50% inhibitory drug concentration (IC_{50}) and the Hill coefficient for drug action as parameter values (Canini and Perelson, 2014; Shen et al., 2008). We fit the observed response of infection to EFV using this approach to estimate d across a range of EFV concentrations. Fit of the model to the cellfree data using wild type, EFVsensitive HIV showed a monotonic decrease with IC_{50} = 2.9 nM and Hill coefficient of 2.1 (Figure 2C, black line).
We next dialed in EFV to tune λ/d in coculture infection. To obtain the number of infected target cells, and specifically exclude donor cells or donortarget cell fusions, target cells were marked by the expression of mCherry. Donor cells were stained with the vital stain Cell Trace Far Red (CTFR). The concentration of live infected cells was determined after 2 days in coculture with infected donors. Live infected cells were identified based on the absence of cell death indicator dye DAPI fluorescence, and presence of GFP. The input of infected donor cells was excluded from the count of infected cells based on the absence of mCherry fluorescence. Donortarget cell fusions were excluded by excluding CTFRpositive cells (see Figure 2—figure supplement 4 for gating strategy).
While the percent of infected cells was reduced with drug, the concentration of live infected cells increased (Figure 2—figure supplement 4). We observed a peak in the number of live infected target cells at 4 nM EFV (Figure 2D). We then fit the number of live infected cells using Equation (3), where P_{λ} was multiplied by the input number of target cells per ml (10^{6} cells/ml) to obtain the predicted number of live infected cells per ml of culture. This was done to constrain r in the model, which strongly determines the amplitude of P_{λ/d} as described above. Equation (3) best fit the behaviour of infection when r = 0.22 and q = 0.17, resulting in a peak at 4.8 nM EFV (Figure 2D, black line). Hence an infection optimum is present in the cell line infection system.
In order to determine whether the fitted r and q values were within a reasonable range, we measured these values experimentally. To measure r, we infected with cellfree HIV to avoid the broad distribution of HIV copy numbers observed in celltocell spread, and determined the fraction of live infected cells P_{λ} (Figure 2—figure supplement 5A). We then determined the mean number of HIV copies per cell λ for the same set of experiments corrected by the efficiency of detection (Figure 2—figure supplement 5B). The parameter r was calculated as ln(1P_{λ})/λ (Supplementary file 2). To measure q, we blocked cell division using serum starvation to measure differences in cell concentration due to cell death only, and not due to proliferation of uninfected cells (Figure 2—figure supplement 5C). We then infected with cellfree HIV and measured L_{λ}, defined as the fraction of live cells remaining upon infection with λ HIV DNA copies relative to infection blocked with EFV (see below). To specifically detect the decrease in live cells as a result of events downstream of reverse transcription, we compared infected cells to cells exposed to the same virus concentration but treated with 40 nM EFV, a drug concentration where infection by cellfree virus is negligible (Figure 2—figure supplement 3). q was then calculated as ln(L_{λ})/λ, where L_{λ} was the probability of a cell to live given transmission with λ copies (Supplementary file 2). Measured r and q values were 0.28 ± 0.08 and 0.15 ± 0.07 (mean ±standard deviation), respectively. The solution to Equation (3) using these values showed similar behavior to the solution with the fitted values for wildtype HIV infection, indicating that the fitted values gave a reasonable approximation of the behavior of the system (Figure 2D, dashed green line).
In order to investigate the dynamics of cell depletion due to celltocell HIV spread and its modulation by the addition of an inhibitor, we performed timelapse microscopy over a two day infection window. While infection parameters were different due to the constraints of visualizing cells (Materials and methods), the general trend from the data was deterioration in the number of live cells in the timelapse culture starting at 1 day postinfection when no drug was added. The deterioration in live cell numbers was averted by the addition of EFV (Figure 2—figure supplement 6).
We next investigated whether an infection optimum occurs with EFVresistant HIV. To derive the resistant mutant, we cultured wildtype HIV in our reporter cell line in the presence of EFV. We obtained the L100I partially resistant mutant. We then replaced the reverse transcriptase of the wildtype molecular clone with the mutant reverse transcriptase gene (Materials and methods). We derived d_{mut} for the L100I mutant using cellfree mutant infection (Figure 3A, see Figure 3—figure supplement 1 for logarithmic yaxis plot). The L100I mutant was found to have an IC_{50} = 29 nM EFV and a Hill coefficient of 2.0 (Figure 3A, black line).
We next performed coculture infection (see Figure 3—figure supplement 2 for gating strategy). Similarly to wildtype HIV coculture infection, there was a peak in the number of live infected target cells for the L100I mutant infection. However, the peak in live infected cells was shifted to 40 nM EFV (Figure 3B). Fits were obtained to Equation (3) using d_{mut} values and λ measured for wildtype infection. The fits recapitulated the experimental results when r = 0.29 and q = 0.13, with a fitted peak at 45 nM EFV (Figure 3B, black line). The solution to Equation (3) using the measured values for r and q showed a similar pattern to that obtained with the fitted values (Figure 3B, dashed green line). We note that both wild type and mutant coculture infection has data points above the fit line at the highest drug concentrations. This may be a limitation of our model at drug values much higher than observed at the infection optimum. In this range of drug values, our model predicts a more pronounced decline in the number of infected cells than is observed experimentally.
In order to examine whether a peak in live infected targets can be obtained with an unrelated inhibitor, we used the HIV neutralizing antibody b12. This antibody is effective against celltocell spread of HIV (Baxter et al., 2014; Reh et al., 2015). We obtained a peak in live infected cells at 5 ug/ml b12 (Figure 4). The b12 concentration that resulted in a peak number of live infected cells was the same for wildtype virus and the L100I mutant, showing that L100I mutant fitness gain was EFV specific. In contrast, cellfree infection in the face of b12 showed a sharp and monotonic drop in live infected cells for both wild type and mutant virus (Figure 4—figure supplement 1).
While the RevCEM cell line is a useful tool to illustrate the principles governing the formation of an infection optimum, the sensitivity of such an optimum to parameter values would make its presence in primary HIV target cells speculative. We therefore investigated whether a fitness optimum occurs in primary human lymph node cells, the anatomical site which would be most likely to have a high force of infection. We derived human lymph nodes from HIVnegative individuals from indicated lung resections (Supplementary file 3), cellularized the lymph node tissue using mechanical separation, and infected the resulting lymph node cells with HIV. A fraction of the cells was infected by cellfree virus and used as infected donor cells. We added these to uninfected target cells from the same lymph node to test coculture infection, and detected the number of live infected cells 4 days postinfection with the L100I EFVresistant mutant in the face of EFV. We detected the number of live infected cells by the exclusion of dead cells with the fixable death detection dye eFluor660 followed by single cell staining for HIV Gag using antip24 antibody (Figure 5A).
In each of the lymph nodes tested, we observed a peak in live infected cells at intermediate EFV concentrations. Lymph node cells from participant 205 showed a peak of live infected cells at 100 nM EFV (Figure 5A, first row). The infection optimum in the lymph node cells of study participant 257 was visible as a plateau between 50 and 200 nM EFV. In the presence of EFV, there was a decrease in the fraction of dead cells that was offset by a similar increase in the fraction of live infected cells for lymph nodes from all participants. There were more overall detectable cells with EFV, resulting in differences in the absolute concentrations of live infected cells being larger than the differences in the fractions of live infected cells between EFV and nondrugtreated cells (right two columns in Figure 5A, with absolute number of live infected cells shown in parentheses in the flow cytometry plots). This is most likely due to cells which died early becoming fragments and so being excluded from the total population in the absence of EFV. Peaks in the number of live infected cells in the face of drug may be specific to lymph node derived cells. Celltocell infection of peripheral blood mononuclear cells (PBMC) with wildtype HIV showed a slight peak at a very low EFV concentration in cells from one blood donor, which was not repeated in cells from two other donors (Figure 5—figure supplement 1).
We used a lymph node from study participant 251, where we obtained more cells, to examine the number of HIV DNA copies per cell. Cells from this lymph node showed an infection optimum at 50 nM EFV (Figure 5A, third row). To detect the effect of EFV on integrations per cell, we sorted single cells based on p24positive signal, decrosslinked to remove the fixative (Materials and methods), then divided each cell lysate into four wells. Using fewer wells saved reagents without changing sensitivity, as demonstrated in the ACH2 cell line (Figure 2—figure supplement 1C). We detected HIV DNA copies by PCR 2 days postinfection. We observed multiple DNA copies in EFVuntreated lymph node cells. The number of copies decreased with EFV (Figure 5B). We corrected for sensitivity of detection as quantified in ACH2 cells (Materials and methods). The corrected numbers were 21 HIV DNA copies with no drug, and five copies in the presence of EFV at the infection optimum (Figure 5C, see Figure 5—figure supplement 2 for histograms of raw HIV DNA copy numbers per cell). Hence, the decrease in the number of copies still results in sufficient copies to infect the cell.
Since L100I does not often occur in the absence of other drug resistance mutations according to the Stanford HIV Drug Resistance Database (Rhee et al., 2003), we repeated the experiment with the K103N mutant, a frequently observed mutation in virologic failure with a higher level of resistance to EFV relative to the L100I mutant. We used cellfree infection to obtain drug inhibition per virion at each level of EFV, which we denote d_{103} (Figure 6A, see Figure 6—figure supplement 1 for logarithmic yaxis plot). The fits showed a monotonic decrease with IC_{50} = 26.0 nM and Hill coefficient of 1.5 (Figure 6A, black line). We then proceeded to use the K103N mutant in coculture infection, using cells from two different lymph nodes in different experiments (see Figure 6—figure supplement 2 for results of individual experiments). We observed an infection optimum with EFV in lymph node cells. The peak in the number of live infected cells in the presence of drug was between 80 and 160 nM EFV (Figure 6B). We fit the experimental data with Equation (3) using d_{103} values and the number of DNA copies in the absence of drug measured for L100I infection. We did not calculate the predicted number of infected cells for P_{λ/d} values since the lymph node is a complex environment containing different cell subsets (Sallusto et al., 1999) and the number of infectable target cells at the start of infection is difficult to determine. Hence, we normalized both the experimental number of live infected cells and the P_{λ/d} values from Equation (3) to the maximum value in each case. The fits recapitulated the experimental results when r = 0.91 and q = 0.15, with a fitted peak at 90 nM EFV (Figure 6B, black line). The q value matched the measured result in the cell line, while the r value was much higher. However, the fitted r value in this case is not expected to be accurate since we were unable to constrain it with the number of infected cells relative to the starting number of target cells.
To examine if the observed peak in live cells may be due to EFV alone, we measured cell viability in lymph node cells from one of the study participants used in the above experiment as a function of EFV without infection. No clear dependence on EFV in the absence of infection was detected (Figure 6—figure supplement 3).
Discussion
The optimal virulence concept in ecology proposes that virulence needs to be balanced against host survival for optimal pathogen spread (Bonhoeffer et al., 1996; Bonhoeffer and Nowak, 1994; Gandon et al., 2001; Jensen et al., 2006). At the cellular level, this implies that the number of successfully infected cells may increase when infection virulence is reduced. The current study is, to our knowledge, the first to address this question experimentally at the level of individual cells infected with HIV.
Using a model where cells are infected and die in a probabilistic way, we found that there were two possible outcomes of partially inhibiting infection. In the case where cells were infected by single infection attempts, inhibition always led to a decline in the number of live infected cells, since inhibition reduced the number of infections per cell from one to zero. In contrast, in the case of multiple infection attempts per cell, the possibility existed that inhibition reduced the number of integrating HIV DNA copies, without extinguishing infection of the cell completely. If each HIV DNA copy increases the probability of cell death, reducing the number of HIV DNA copies without eliminating infection should lead to an increased probability of infected cell survival. This would consequently lead to an increase in the number of live infected cells.
We investigated the outcome of partial inhibition of infection in both a cell line and primary lymph node cells. In both systems, we observed that there was a peak in live infected cell number at intermediate inhibitor concentrations. This correlated to a decreased number of viral DNA copies per cell. Further increasing inhibitor concentration led to a decline in live infected cell numbers, and infecting with EFV resistant mutants shifted the peak in live infected number to higher EFV concentrations. Our model as described by Equation (3) reproduced the essential behaviour of the experimental results. Construction of the model assumed independence of productive infection and cell death. However, as shown in the Appendix 1, an equivalent model can be constructed assuming a dependence of cell death on infection. Neither model accurately captures infection dynamics at highdrug concentrations, away from the infection optimum. In this range, where the number of infection attempts per cell is much lower than 1, infection declined more slowly with drug than predicted. The model can be further refined using a distribution for the number of DNA copies per cell. Moreover, the probability of death per HIV DNA copy we denote q may be dependent on how many infection attempts preceded the current infection attempt, and the model can be improved by measuring this dependence.
Physiologically, an infection optimum in the face of an antiretroviral drug may be important in HIV infection of lymph node cells and may be less pronounced in cells from peripheral blood. We used EFV in our study since it is a common component of first line antiretroviral therapy, with frequent drug resistance mutations. However, the infection optimum we describe should occur with other classes of antiretroviral drugs, since all drugs should decrease the multiplicity of infection between cells. In terms of modeling, a future therapy component such as the integrase inhibitor dolutegravir would exert its effect on r and not λ in our model. However, the effect is symmetrical since e^{(λ/d)r} = e^{λ(r/d)}. The more complex outcome of partial inhibition of infection should also be considered in other infections where multiple pathogens infect one cell and host cell death is a possible outcome (Mahamed et al., 2017).
These observations reinforce previous results showing that successful completion of reverse transcription leads to cellular cytotoxicity. In addition to HIV cytotoxicity caused by viral integrations through the mechanism of double strand breaks (Cooper et al., 2013), other mechanisms of HIVinduced death are also present, including IFI16dependent innate immune system sensing of abortive reverse transcripts following nonproductive infection of resting T cells (Doitsh et al., 2014; Monroe et al., 2014). The experiments presented here reflect the effect of partial inhibition on productive infection of HIV target cells, which mostly consist of activated Tcell subsets, not resting T cells. More complex models would be needed to decipher the effect of partial inhibition of HIV infection on resting Tcell numbers and the outcome of this in terms of available Tcell targets in future infection cycles.
This study may have implications for the establishment of viral reservoirs in the context of poorly controlled infections, infections with some degree of drug resistance, or infections where some replication may take place in the face of ART, since infected cell survival is a prerequisite for longterm persistence. The clinical implications of an infection optimum in the presence of EFV with EFVsensitive HIV strains are likely to be negligible, since the drug concentrations at which the infection optimum occurs are extremely low. However, for EFVresistant HIV, the infection optimum shifts to the range of EFV concentrations observed in lymph nodes (~100 nM) (Fletcher et al., 2014b), and can be expected to shift to even higher EFV concentrations with more resistant mutants. As EFV has a longer halflife than the other antiretroviral drugs coformulated with it, it may be the only agent present in partially adherent individuals for substantial periods of time (Taylor et al., 2007). Therefore, partial inhibition of HIV infection with EFV may provide a surprising advantage to EFV resistant mutants, and may allow individuals failing therapy to better transmit drug resistant strains.
Materials and methods
Ethical statement
Request a detailed protocolLymph nodes were obtained from the field of surgery of participants undergoing surgery for diagnostic purposes and/or complications of inflammatory lung disease. Informed consent was obtained from each participant, and the study protocol approved by the University of KwaZuluNatal Institutional Review Board (approval BE024/09). Blood for PBMC was obtained from healthy blood donors under the same study protocol.
Inhibitors, viruses and cell lines
Request a detailed protocolThe following reagents were obtained through the AIDS Research and Reference Reagent Program, National Institute of Allergy and Infectious Diseases, National Institutes of Health: the antiretroviral EFV; RevCEM cells from Y. Wu and J. Marsh; HIV molecular clone pNL43 from M. Martin; ACH2 cells from T. Folks. Cellfree viruses were produced by transfection of HEK293 cells with pNL43 using TransITLT1 (Mirus, Madison, WI ) or Fugene HD (Roche, RischRotkreuz, Switzerland) transfection reagents. Virus containing supernatant was harvested after 2 days of incubation and filtered through a 0.45 µm filter (Corning, New York, NY). b12 antibody was produced from transfecting HEK293 cells with a b12 expression plasmid (expressed under a CMV promoter on a pHAGE6 lentiviral plasmid backbone, gift from A. Balazs), followed by harvesting of cell supernatant and purification at the California Institute of Technology protein expression core. The number of virus genomes in viral stocks was determined using the RealTime HIV1 viral load test (Abbott Diagnostics, Santa Clara, CA). For r and q measurement, 0.45 µm filtered cellfree supernatants from infected RevCEM cells were used, to include any secreted factors which may modulate celldeath. The L100I and K103N mutants were evolved by serial passages of wildtype NL43 in RevCEM cells in the presence of 20 nM EFV. After 16 days of selection, the reverse transcriptase gene was cloned from the proviral DNA and the mutant reverse transcriptase gene was inserted into the NL43 molecular clone. RevCEM clones E7 and G2 used in this study were generated as previously described (Boullé et al., 2016). Briefly, the E7 clone was generated by subcloning RevCEM cells at singlecell density. Surviving clones were subdivided into replicate plates. One of the plates was screened for the fraction of GFP expressing cells upon HIV infection using microscopy, and the clone with the highest fraction of GFPpositive cells was selected. To generate the G2 clone, E7 cells were stably infected with the mCherry gene under the EF1α promoter on a pHAGE2based lentiviral vector (gift from A. Balazs), subcloned, and screened for >99% mCherrypositive cells. All cell lines not authenticated, and mycoplasma negative. Cell culture and experiments were performed in complete RPMI 1640 medium supplemented with LGlutamine, sodium pyruvate, HEPES, nonessential amino acids (Lonza, Basel, Switzerland), and 10% heatinactivated FBS (GE Healthcare BioSciences, Pittsburgh, PA).
Primary cells
Request a detailed protocolLymph node cells were obtained by mechanical separation of lymph nodes and frozen at 5 × 10^{6} cells/ml in a solution of 90% FBS and 10% DMSO with 2.5 μg/ml Amphotericin B (Lonza). Cells were stored in liquid nitrogen until use, then thawed and resuspended at 10^{6} cells/ml in complete RPMI 1640 medium supplemented with LGlutamine, sodium pyruvate, HEPES, nonessential amino acids (Lonza), 10% heatinactivated FBS (Hyclone), and IL2 at 5 ng/ml (PeproTech). Phytohemagglutinin at 10 µg/ml (SigmaAldrich, St Louis, MO) was added to activate cells. PBMCs were isolated by density gradient centrifugation using Histopaque 1077 (SigmaAldrich) and cultured at 10^{6} cells/ml in complete RPMI 1640 medium supplemented with LGlutamine, sodium pyruvate, HEPES, nonessential amino acids (Lonza), 10% heatinactivated FBS (GE Healthcare BioSciences, Pittsburgh, PA), and IL2 at 5 ng/ml (PeproTech, Rocky Hill, NJ). Phytohemagglutinin at 10 µg/ml (SigmaAldrich) was added to activate cells. For both primary cell types, donor cells for coculture infection were cultured for one day then infected by cellfree virus, while target cells were cultured for three days and infected with either cellfree HIV or infected donor cells.
Subcloning of ACH2 cells
Request a detailed protocolCells from the parental ACH2 cell line were diluted to 10 cells/ml in conditioned medium, with conditioned medium generated by culturing ACH2 cells to 10^{6} cells/ml, then filtering through a 0.22 µm filter (Corning). 25 μl of the diluted cell suspension was then distributed to each well of a Greiner μClear 384well plate (mean of 0.5 cells per well). Clones were cultured for 3 weeks, where each week an additional 25 μl of conditioned medium was added to each well. Clones were detected in 5% of wells and two clones, designated D6 and C3, were randomly chosen and further expanded.
Infection
Request a detailed protocolFor a cellfree infection of RevCEM clones, PBMC and lymph node cells, 10^{6} cells/ml were infected with 2 × 10^{8} NL43 viral copies/ml (~20ng p24 equivalent) for 2 days. For coculture infection, infected cells from the cellfree infection were used as the donors and cocultured with 10^{6} cells/ml target cells. For RevCEM clones, 2% infected donor cells were added to uninfected targets and cocultured for 2 days in tissue culture experiments, and 20% infected donor cells were added to uninfected targets and cocultured for 2 days for timelapse experiments. For lymph node cells and PBMCs, a ratio of 1:4 donor to targets cells was used. Infection was over 2 days in PBMC infection and for 4 days for infection of lymph node cells.
Staining and flow cytometry
Request a detailed protocolTo determine the number of live infected cells in reporter cell line experiments, E7 RevCEM reporter cells were infected as above used as donor cells. Prior to coincubation with target cells, donor cells were stained with CellTrace Far Red (CTFR, Thermo Fisher Scientific, Waltham, MA) at 1 µM and washed according to manufacturer’s instructions. The G2 mCherrypositive reporter cells were used as infection targets, and cocultured with 2% infected donor cells for 2 days. The coculture infection was pulsed with 100 ng/ml DAPI (SigmaAldrich) immediately before flow cytometry and the number of live infected targets cells was determined by the number of DAPI negative, CTFR negative and mCherry and GFP double positive cells on a FACSAria Fusion machine (BD Biosciences, Sparks, MD) using the 355, 488 and 633 nm laser lines. For cellfree infections where fewer fluorescence channels were used, a pulse of 300 nM of the farred live cell impermeable dye DRAQ7 (Biolegend, San Diego, CA) immediately before flow cytometry was substituted for DAPI, and live infected cells detected as the number of DRAQ7negative, GFPpositive cells on a FACSCaliber machine using 488 and 633 nm laser lines. Lymph node cells were resuspended in 1 ml of phosphate buffered saline (PBS) and stained at a 1:1000 dilution of the eFlour660 dye (Thermo Fisher Scientific) according to the manufacturer’s instructions. Cells were then fixed and permeabilized using the BD Cytofix/Cytoperm Fixation/Permeabilization kit (BD Biosciences) according to the manufacturer’s instructions. Cells were then stained with antip24 FITC conjugated antibody (KC57, Beckman Coulter, Brea, CA). Live infected lymph node cells were detected as the number of eFluor660negative, p24positive cells. Cells were acquired with a FACSAriaIII or FACSCaliber machine (BD Biosciences) using 488 and 633 nm laser lines. Results were analysed with FlowJo 10.0.8 software. For singlecell sorting to detect the number of HIV DNA copies per cell, cells were singlecell sorted using 85 micron nozzle in a FACSAriaIII machine. GFPpositive, DRAQ7negative RevCEM clones were sorted 1 day postinfection into 96 well plates (Biorad, Hercules, CA) containing 30 µl lysis buffer (2.5 µl 0.1M Dithiothreitol, 5 µl 5% NP40 and 22.5 µl molecular biology grade water [Kurimoto et al., 2007]). For experiments to determine the number of HIV DNA copies to measure r and q, the DRAQ7negative subset was sorted. Fixed, p24positive, eFluor660negative lymph node cells were singlecell sorted two days postinfection into 96well plates containing 5 µl of PKD buffer (Qiagen, Hilden, Germany) with 1:16 proteinase K solution (Qiagen) (Thomsen et al., 2016). Sorted plates were snap frozen and kept at −80°C until ready for PCR. For analysis by flow cytometry, a minimum of 50,000 cells were collected per data point.
Timelapse microscopy and image analysis
Request a detailed protocolFor imaging infection by timelapse microscopy, cell density was reduced to 5 × 10^{4} cells/ml and cells were attached to ployllysine (SigmaAldrich) coated optical sixwell plates (MatTek, Ashland, MA). Infections with and without EFV were imaged in tandem using a Metamorphcontrolled Nikon TiE motorized microscope with a Yokogawa spinning disk with a 20x, 0.75 NA phase objective in a biosafety level three facility. Excitation sources were 488 (GFP) and 561 (mCherry) laser lines and emission was detected through a Semrock Brightline quad band 440–40/521–21/60734/70045 nm filter. Images were captured using an 888 EMCCD camera (Andor, Belfast, UK). Temperature (37°C), humidity and CO_{2} (5%) were controlled using an environmental chamber (OKO Labs, Naples, Italy). Fields of view were captured every 20 min. To facilitate automated image analysis of timelapse experiment data, mCherry expressing G2 clone cells were used as targets and E7 clone cells used as infected donors. The number of live cells was measured as the number of cells expressing mCherry since intracellular mCherry protein is soluble and hence lost upon cell death when cellular membrane integrity is compromised. The number of live infected cells was measured as the number of cells expressing both mCherry and GFP. Three independent experiments were performed. Movies were analyzed using custom code developed with the Matlab R2014a Image Analysis Toolbox. Images in the mCherry channel were thresholded and the imfindcircle function used to detect round objects within the cell radius range. Cell centers were found. GFP signal underwent the same binary thresholding. The number of mCherrypositive 16 pixel^{2} squares around the cell centers was used as the as the number of total target cells at each timepoint, and the number of squares double positive for fluorescence in the GFP channel was used as the number of infected target cells.
Determination of HIV DNA copy number in individual cells
Request a detailed protocol96well plates of cells previously sorted at 1 cell per well were thawed at room temperature and spun down. Fixed cells were decrosslinked by incubating in a thermocycler at 56°C for 1 hr. The lysate from each well was split equally over 10 wells (2.5 µl each well after correction for evaporation) for E7 RevCEM or four wells (6.8 µl each well after correction for evaporation) for lymph nodes, containing 50 µl of Phusion hot start II DNA polymerase (New England Biolabs, Ipswich, MA) PCR reaction mix (10 µl 5X Phusion HF buffer, 1 µl dNTPs, 2.5 µl of the forward primer, 2.5 µl of the reverse primer, 0.5 µl Phusion hot start II DNA polymerase, 2.5 µl of DMSO and molecular biology grade water to 50 µl reaction volume). Two rounds of PCR were performed. The first round reaction amplified a 700 bp region of the reverse transcriptase gene using the forward primer 5’ CCTACACCTGTCAACATAATTGGAAG 3’ and reverse primer 5’ GAATGGAGGTTCTTTCTGATG 3’. Cycling program was 98°C for 30 s, then 34 cycles of 98°C for 10 s, 63°C for 30 s and 72°C for 15 s with a final extension of 72°C for 5 min. 1 µl of the first round product was then transferred into a PCR mix as above, with nested second round primers (forward 5’ TAAAAGCATTAGTAGAAATTTGTACAGA 3’, reverse 5’ GGTAAATCCCCACCTCAACAGATG 3’). The second round PCR amplified a 550 bp product which was then visualized on a 1% agarose gel. PCR reactions were found to work best if sorted plates were thawed no more than once, and plates which underwent repeated freezethaw cycles showed poor amplification.
Correction of raw number of detected DNA copies for detection sensitivity
Request a detailed protocolA stochastic simulation in Matlab was used to generate a distribution for the number of positive wells per cell for each mean number of DNA copies per cell λ. The probability for a DNA copy to be present within a given well and be detected was set as σ/w, where σ was the detection sensitivity calculated as the number of ACH2 with detectable integrations divided by the total number of ACH2 cells assayed (38/166, σ = 0.23), and w was the number of wells. A random number m representing DNA copies per cell from a Poisson distribution with a mean λ was drawn, and a vector R of m random numbers from a uniform distribution was generated. If there existed an element R_{i} of the vector with a value between 0 and σ/w, the first well was occupied. If an element existed with a value between σ/w+γ and 2(σ/w), where γ <<1, the second well was occupied, and if between (σ/w+γ)(n1) and n(σ/w), the nth well was occupied. The sum of wells occupied at least once was determined, and the process repeated j times for each λ, where j was the number of cells in the experimental data. A least squares fit was performed to select λ which best fit the experimental results across well frequencies, and mean and standard deviation for λ was derived by repeating the simulation 10 times.
Fit of the EFV response for single infections using IC_{50} and Hill coefficient
Request a detailed protocolTo obtain d, we normalized Equation (2) by the fraction of infected cells in the absence of drug (Sigal et al., 2011) to obtain Tx = (infected targets with EFV)/ (infected targets no EFV)= ((1 (1 r)^{λ/d}) (1q)^{λ/d})/ ((1 (1 r)^{λ}) (1q)^{λ}). We approximate the result at small r, q to Tx = (1 e^{rλ/d}) e^{qλ/d}/ (1 e^{rλ}) e^{qλ}=eq^{qλ(11/d)} ((1 e^{rλ/d}) / (1 e^{rλ})). Expanding the exponentials we obtain Tx = (1 + qλ(1–1/d)) ((rλ/d)/rλ) = (1 + qλ(1–1/d))(1/d). We note that at λ <1, qλ(1–1/d)<<1, and hence Tx ≅ 1/d. Tx was measured from the experiments to obtain d values at the EFV concentrations used for cellfree infection, where λ < 1. To obtain a fit of d as a function of the concentration of drug that gives halfmaximal inhibition (IC_{50}) and Hill coefficient (h) for EFV, we used the relation for the fraction cells remaining infected in the face of drug (Canini and Perelson, 2014), whose definition is equivalent to Tx at λ < 1:
Measurement of r
Request a detailed protocolCellfree supernatant used in infection was derived as follows: 10^{6} cells/ml were infected with 2 × 10^{8} NL43 viral copies/ml (~20ng p24 equivalent) for 2 days. Thereafter 0.2% of the infected cells from the cellfree infection were added to 10^{6} cells/ml target cells. The infected supernatant from the coculture 2 days postinfection was filtered using a 0.45 µm filter (Corning) and added to cells at a 1:8 dilution, where the dilution was calibrated to result in nonsaturating infection in terms of GFP expression. A fraction of the cells were sorted into lysis buffer at one cell per well 1 day post infection, split over four wells, and PCR performed as described above to determine HIV copy number per cell. The remaining cells from the same infection were used to determine frequency of DRAQ7negative, GFPpositive cells 2 days postinfection using flow cytometry.
Measurement of q
Request a detailed protocolCellfree supernatant used in infection was derived as for r, except that 1 day before harvesting of the viral supernatant from infected cells, infected cells were washed twice with PBS and serumfree growth medium added. At the same time, the target cells for the infection were washed twice with PBS and serumfree growth medium added. Cells were split into two wells, and EFV to a final concentration of 40 nM was added to one of the wells. Supernatant was harvested and filtered as described for r, and added to cells at a 1:2 dilution. A fraction of the cells were sorted into lysis buffer at one cell per well 1 day postinfection, split over four wells, and PCR performed as described above to determine HIV copy number per cell. The remaining cells from the same infection were used to determine the frequency of live and dead cells two days postinfection. The concentration of live cells was measured using the TC20^{TM} automated cell counter (BioRad) with trypan blue staining (Lonza).
Appendix 1
Supplementary Mathematical Analysis for Jackson et al.
I Derivation of infection model assuming independence of infection and cell death
We consider a model of infection where is the probability of one HIV DNA copy to successfully infect the cell and is the probability of the cell to die as a result of the infection attempt. The probability of being infected when virions enter the cell is:
If the number of infection attempts, $n$, is Poisson distributed with mean $\lambda $, the fraction of cells that are productively infected is:
We introduce a drug strength value $d$, where $d=1$ in the absence of drug and $d>1$ in the presence of drug. In the presence of drug, $\lambda $ is decreased to $\lambda /d$. The drug therefore tunes $\lambda $. The probability of a cell to be infected and live given drug strength $d$ is therefore:
Experimentally, we measure $P}_{\lambda$, the number of infected cells, and $L}_{\lambda$, the number of live cells remaining two days postinfection for the cell line.
Assuming a Poisson distribution of the number of infection attempts, the probability of being productively infected is
Similarly, the probability of a cell to live is
II Derivation of infection model assuming dependence of infection and cell death
Cellular infection and death may be the result of independent processes. For example, cell death may be a programmed response to DNA damage induced by the infection attempt and therefore may occur regardless of whether or not the cell is successfully infected. Infection and death may also be dependent. One example of dependence is the cytotoxicity mediated by expressed viral proteins after integration into the host genome.
In this case, the probability of productive infection and survival of the cell given $n$ attempts is
Here, ${r}^{\prime}$ and ${q}^{\prime}$ are the probabilities for infection and death respectively. Their relationship to the experimentally measured $r$ and $q$ are discussed in the next section. ${\left[{r}^{\prime}(1{q}^{\prime})\right]}^{m}$ is the probability that $m$ infection attempts are successful and none of the attempts triggers cell death, while ${(1{r}^{\prime})}^{nm}$ accounts for $nm$ unsuccessful infections. To obtain the total probability of productive infection without cell death given $n$ attempts, we sum over all cases with ${m}^{\prime}\ge 1$ sucessful infections.
Realizing that:
we can simplify the above to:
Averaging this over a Poisson distribution with mean$\lambda $, we find:
III Equivalence of the two models
The two models predict same dependence on $\lambda $ (the difference of two exponentials) but the parameters $r,q$ and ${r}^{\prime},{q}^{\prime}$ play slightly different roles. The models are the same if:
Hence the two models are a reparametrization of the same process. The experimental measurement of $r$ and $q$, as described in the main manuscript, consists of estimating $r$ as the fraction of infected cells when the fraction of dead cells is low due to a low number of measured DNA copies per cell. The death rate $q$ is derived by measuring the number of live cells remaining when the number of HIV DNA copies per cell is known. If the probability of death is dependent on infection, then death is only of infected cells and is a product of the probability to be infected and die (${r}^{\prime}{q}^{\prime}$). Likewise, the underlying probability of infection ${r}^{\prime}$ will account for the fraction dead cells ($q$) and any cells which are infected but not dead ($r(1q)$).
References

Crosslinking CD4 by human immunodeficiency virus gp120 primes T cells for activationinduced apoptosisJournal of Experimental Medicine 176:1099–1106.https://doi.org/10.1084/jem.176.4.1099

Macrophage infection via selective capture of HIV1infected CD4+ T cellsCell Host & Microbe 16:711–721.https://doi.org/10.1016/j.chom.2014.10.010

The curse of the pharaoh: the evolution of virulence in pathogens with long living propagulesProceedings of the Royal Society B: Biological Sciences 263:715–721.https://doi.org/10.1098/rspb.1996.0107

Mutation and the evolution of virulenceProceedings of the Royal Society B: Biological Sciences 258:133–140.https://doi.org/10.1098/rspb.1994.0153

CD4+ T cell depletion during all stages of HIV disease occurs predominantly in the gastrointestinal tractThe Journal of Experimental Medicine 200:749–759.https://doi.org/10.1084/jem.20040874

Viral kinetic modeling: state of the artJournal of Pharmacokinetics and Pharmacodynamics 41:431–443.https://doi.org/10.1007/s1092801493633

Multiploid inheritance of HIV1 during celltocell infectionJournal of Virology 85:7169–7176.https://doi.org/10.1128/JVI.0023111

Defining HIV and SIV reservoirs in lymphoid tissuesPathogens and Immunity 1:68.https://doi.org/10.20411/pai.v1i1.100

Multiplicity of human immunodeficiency virus infections in lymphoid tissueJournal of Virology 78:8942–8945.https://doi.org/10.1128/JVI.78.16.89428945.2004

HIV1 transmission, by stage of infectionThe Journal of Infectious Diseases 198:687–693.https://doi.org/10.1086/590501

HIV1 cell to cell transfer across an Envinduced, actindependent synapseThe Journal of Experimental Medicine 199:283–293.https://doi.org/10.1084/jem.20030648

BookVirus dynamics: mathematical principles of immunology and virology: mathematical principles of immunology and virologyUK: Oxford University Press.

A sensitive, quantitative assay for human immunodeficiency virus type 1 integrationJournal of Virology 76:10942–10950.https://doi.org/10.1128/JVI.76.21.1094210950.2002

Modelling viral and immune system dynamicsNature Reviews Immunology 2:28–36.https://doi.org/10.1038/nri700

Virus fitness: concept, quantification, and application to HIV population dynamicsCurrent Topics in Microbiology and Immunology 299:83–140.

Human immunodeficiency virus reverse transcriptase and protease sequence databaseNucleic Acids Research 31:298–303.https://doi.org/10.1093/nar/gkg100

Retroviruses can establish filopodial bridges for efficient celltocell transmissionNature Cell Biology 9:310–315.https://doi.org/10.1038/ncb1544

Inefficient human immunodeficiency virus replication in mobile lymphocytesJournal of Virology 81:1000–1012.https://doi.org/10.1128/JVI.0162906

Rates of HIV1 transmission per coital act, by stage of HIV1 infection, in Rakai, UgandaThe Journal of Infectious Diseases 191:1403–1409.https://doi.org/10.1086/429411

HIV1 Tat potentiates TNFinduced NFkappa B activation and cytotoxicity by altering the cellular redox stateThe EMBO Journal 14:546–554.

Human immunodeficiency virus evolution towards reduced replicative fitness in vivo and the development of AIDSProceedings of the Royal Society B: Biological Sciences 274:2481–2491.https://doi.org/10.1098/rspb.2007.0413

Revdependent indicator T cell lineCurrent HIV Research 5:394–402.https://doi.org/10.2174/157016207781024018

Lymphoid tissue structure and HIV1 infection: life or death for T cellsTrends in Immunology 33:306–314.https://doi.org/10.1016/j.it.2012.04.002
Decision letter

Wenying ShouReviewing Editor; Fred Hutchinson Cancer Research Center, United States
In the interests of transparency, eLife includes the editorial decision letter and accompanying author responses. A lightly edited version of the letter sent to the authors after peer review is shown, indicating the most substantive concerns; minor comments are not usually included.
Thank you for submitting your article "Incomplete inhibition of HIV infection results in more HIV infected lymph node cells by reducing cell death" for consideration by eLife. Your article has been evaluated by Wenhui Li (Senior Editor) and three reviewers, one of whom is a member of our Board of Reviewing Editors. The reviewers have opted to remain anonymous.
The reviewers have discussed the reviews with one another and the Reviewing Editor has drafted this decision to help you prepare a revised submission.
Summary:
Jackson et al. mathematically predicted and experimentally demonstrated that in lymph node cells, intermediate antiviral (antiHIV) drug concentrations can allow infected cells to live longer by reducing multiplicity of infection. Thus, intermediate drug concentration can presumably facilitate the persistence of HIV. The result is interesting and can have clinical implications.
Reviewer #1's concern is quite substantial. Other reviewers agreed to this concern. "I strongly believe that infection probability r and death probability q should be measured experimentally. This is to eliminate fitting parameters and to ensure quality of the model."
Comments from Reviewers 2 and 3 are attached.
Reviewer #2:
The studies presented by Sigal et al. examine the interrelationship between cell death and the number of DNA copies of HIV that are detected per individual infected cell. Previous studies from a number of groups, including Sigal's, have indicated that celltocell infection leads to infected cells with a higher number of integrated proviruses. The single cell PCR data provided by Sigal et al. support a hypothesis that the celltocell infection not only leads to higher copy number of proviruses, but also generally leads to more rapid cell death. Using drugs or antibodies at low levels that do not fully inhibit infection, they found that low levels of inhibition paradoxically can give rise to a larger number of surviving infected cells. They also appear to find lower rates of cell death under conditions of partial inhibition. The studies are generally well conceived, and the results generally support the conclusions. The overall quality of the data could be improved by showing additional controls and increasing the numbers of cells analyzed to enhance the confidence in the copy number estimates. On the cell death front, it would be helpful for their hypothesis to show a time course with the cell death depicted over time. Lastly, an important detail not mentioned, is whether the use of lymph node cells is critical or not for the phenotypes that they describe. The study may have implications for the establishment of viral reservoirs in the context of poorly controlled infection or infections with some degree of drug resistance.
1) The authors use the ACH2 cell line to test the ability of their assay to detect single copy integrations. Why do they only spread each cell over 4 wells in control studies as opposed to the 10 wells used in their experimental studies? The rationale for using different dilution schemes should be explained. A comparison of the efficiency of PCR at the different dilutions would also be informative.
The ACH2 cell line yields fewer than 1 copy per cell, which may be a limitation of the PCR assay, but also could be reflective of heterogeneity in a cell line that is assumed to be clonal and a uniform karyotype. It is not described how recently the ACH2 cell line that they are using has been cloned. Repeating the ACH2 studies with a subclone would be beneficial.
Figure 2 would be more informative if it showed the results of both cellfree virus infection versus coculture infection. The numbers of cells analyzed in Figure 2C (n=34) does not appear sufficient to provide a robust sense of the distribution of the DNAs in the infected cells. Is this a bimodal distribution? Histogram comparison of cellfree should also be used for comparison.
2) Figure 4.
What happens with partial antibody inhibition of cellfree infection? These studies should be performed with infection from cellfree virus. To further test their hypothesis partial inhibition of cellfree infection should only decrease infected cells.
3) Figure 5.
As controls for cell death measurements, the viability of the uninfected control lymph node cells, treated and untreated, should also be illustrated. In the literature lymph node cells (mostly tonsilar) are very prone to cell death, and it is important to understand to what extent they are measuring virusinduced cell death. In addition, a time course of the cell death observed in the infections may would also be helpful to evaluate their hypothesis that the increase of infected cells in the partially treated cells is due to increased infected cell survival – i.e. decreased death.
4) Are the phenotypes described in the primary lymph node cells also observed for peripheral blood lymphocytes? Or are there differences between the peripheral blood versus the lymph node T cells.
Reviewer #3:
This is an interesting and worthwhile paper that examines the effects of incomplete inhibition of HIV infection with reverse transcription inhibitors. The experimental results presented are in accordance with model predictions. However, the presentation of the model should be improved as detailed below. Also, I found the labeling and numbering of the supplemental figures confusing. I have no substantial concerns and this the paper should be published after minor revision.
1) The terms and concepts in the paper are not clearly defined. In the first paragraph of the Results, you mention each donor to target transmission. It is not clear what a transmission is. Does it refer to viral entry of either a free virus or a virus (or genome) by celltocell transmission? It needs to be defined. Second define what you mean by infection. Does a cell have to produce virus to be considered infected? Does the virus have to integrate or is it sufficient to simply reverse transcribe? Is a latently infected cell an infected cell in your model? Also, when you say drug therapy increases the number of live infected cells do you mean live productively infected cells, live HIV DNA^{+} cells, etc.
2) The assumption that all transmissions have equal probabilities to infect target cells seems to ignore the possibility that some virions carry defective genomes while others do not. This does not seem realistic. In your experimental system is the ratio of HIV RNA to TCID50 close to one so you can ignore defective particles?
3) You assume productive infection and death are independent events. While the events may be independent the probability of death is certainly not independent of whether a cell is productively infected or not. Further the probability of death is time dependent. The probability a cell dies one hour after viral entry (infection?) is clearly quite different than the probability it dies days after infection. You may want to define q as the probability a cell has died by time t after infection and the same for Pλλ where t is 2 days or 4 days for your various experiments.
4) Results, second paragraph. You claim antiretroviral drugs reduce the number of infecting virions. This implies that by infecting you must mean the virus reverse transcribes. In standard viral dynamic models the effects of ART are to reduce the infection probability, i.e., r in your model not λ. Thus it is important to clarify your definitions.
5) You introduce the drug effect as a constant d in the second paragraph of the Results. Later in the paper you make reference to IC_{50} and Hill coefficients for the drug. These need to be tied together in an explicit manner. In viral dynamic modeling the effectiveness of a drug, epsilon, (eps) is introduced where eps=1 is a 100% effective drug, e.g. stops all reverse transcription, and where eps=0 means the drug has no effect. Then a pharmacodynamic model is used to relate eps with the drug concentration, C, e.g. eps= C^n/(IC50^n + C^n). Thus rather than λ/d you could model the drug effect by (1eps)λ and use the above formula to link the effect to drug concentration. See for example, Canini and Perelson, 2014. Equation 3 in your manuscript accomplished the same, but the above seems more straightforward.
6) Results, fourth paragraph. When you measure the number of reverse transcribed copies of viral DNA I assume you are measuring both integrated and unintegrated DNA – please state this explicitly as later you seem to imply you are measuring integrated DNA, e.g. Results, sixth paragraph.
7) The current work focuses only on reverse transcriptase inhibitors. Model therapy uses combination therapy including integrase inhibitors and protease inhibitors. Expanding the discussion of the clinical implications of this work to include combination therapy would be worthwhile.
[Editors' note: further revisions were requested prior to acceptance, as described below.]
Thank you for submitting your article "Incomplete inhibition of HIV infection results in more HIV infected lymph node cells by reducing cell death" for consideration by eLife. Your article has been favorably evaluated by Wenhui Li (Senior Editor) and three reviewers, one of whom, Wenying Shou (Reviewer #1), is a member of our Board of Reviewing Editors. The following individual involved in review of your submission has agreed to reveal his identity: Alan Perelson (Reviewer #3).
The reviewers have discussed the reviews with one another and the Reviewing Editor has drafted this decision to help you prepare a revised submission.
We want you to address the following comments. Note that we are not asking you to do further experiments. For the sake of recordkeeping, we are enclosing all comments.
Reviewer #1:
I am satisfied with their measurements.
In general, I suggest breaking up long sentences (especially the one in Abstract, also the end of the third paragraph in the Introduction) for readability. Sentences that last more than three lines are generally more difficult to read than shorter sentences.
Reviewer #2:
The authors have very rigorously addressed the critiques by performing new experiments and modifying the text. The response to the reviews is particularly comprehensive and thorough. The resulting new data are supportive of the original model and hypothesis. The authors should be commended.
Reviewer #3:
In this revised version on the paper the authors have clarified their definitions by incorporating a supplemental table and additional text describing their model. However, this clarification has raised the following issues that need to be resolved:
1) As defined in the Results, λ is the number of infection attempts, where one attempt is measured as one HIV DNA copy. Thus, λ is an integer. The formula, Equation 1, makes sense when λ is an integer. However, what is determined experimentally is not the integer number of DNA copies but rather the mean number of copies per cell, which is not necessarily an integer. To make sense Equation 1 should be reformulated as a conditional probability of a cell being infected and alive given x attempts, call this P(I  x). What I believe the authors want is the probability of a cell being infected and alive, which is then P(I) =Sum P(I x) (p(x)), where p(x) is the probability of having x infection attempts, and the sum is over all nonnegative values of x. One could then assume, for example, that the number of attempts is Poisson distributed with mean λ. Then the resulting formula for P(I) would involve the mean value of the number of attempts. A better choice would be the empirical distribution that the authors measured. Note that what is actually measured is the number of infected cells that are alive at a certain time, i.e. P(I), and the mean number of attempts, what I called λ above.
2) According to the model q should be independent of the number of attempts. Can this be tested by using different drug concentrations, which would vary the mean number of attempts? If you already have experimental data on this, then include them. Otherwise, just discuss it.
3) In the original manuscript where Equation 1 was derived the authors assumed q was independent of r. In the second paragraph of the Results the authors say q is the probability of a cell to die post reverse transcription, which seems more sensible to me. Under the original assumption of q and r being independent, a cell in which the attempt was unsuccessful would have the same probability of dying as a cell in which reverse transcription succeeded. With q now defined as in the second paragraph of the Results, the fundamental formula given by Equation 1 may not be correct and deserves a more thorough discussion. I think the fundamental process being described is that at each attempt a cell is either infected with probability r, or not infected with probability 1r. If it is infected, i.e. had a successful reverse transcription, then with probability q it remains alive and with probability 1q it dies. Assume there are x attempts. Let z be the number of cells infected after x attempts. Then z is binomially distributed, i.e. z= Bin (x,r). Further, if we are interested in the number of cells infected and alive after x attempts, then P(Ix) = Sum_z=1 to x (x choose r) r^z (1r)^(xz)(1q)^z, where the factor (1q)^z is the probability of a cell surviving after each of the z infections. Simplifying, P(Ix) = [r(1q) + (1r)]^x, which is not the same as Equation 1. The number of attempts, x, is again random and as above one can convert this conditional probability into P(I) by assuming a Poisson (or some other distribution such as the empirically measured one) for the number of attempts p(x).
4) In Supplementary file 1 the authors say the way they measure the probability a cell is alive after λ attempts is by computing the concentration of cells alive with λ attempts divided by the concentration of cells alive with no transmissions. However, since the concentration of cells alive with no transmissions can approach zero as time goes on it is clear that this fraction need not be less than 1. Also, as noted above the authors need to use conditional probabilities and derive formulas that involve the mean number of attempts. I would suggest they compute the ratio of the probability a cell is infected to the probability a cell in not infected and alive, which is what they measure, and see if they can derive a formula for q.
https://doi.org/10.7554/eLife.30134.039Author response
Reviewer #1's concern is quite substantial. Other reviewers agreed to this concern. "I strongly believe that infection probability r and death probability q should be measured experimentally. This is to eliminate fitting parameters and to ensure quality of the model."
In order to address this important point, we designed two sets of experiments, one set for measuring r and the other for q. The definition for r and q is now sharpened as suggested by reviewer 3 in Supplementary file 1 (as described in the response to reviewer 3 comments).
Our general approach was to use cellfree infection to avoid the high multiplicity of infection and broad distribution of HIV copy number observed in celltocell spread. We then measured the fraction of live infected cells (P_{λ}) or cell viability (L_{λ}) postinfection. For each set of experiments we determined the number of infection attempts λ, which was measured as the number of total HIV DNA copies per cell, both integrated and unintegrated. This enabled us to calculate p and q (see below). For q, we compared infected cells to cells infected but treated with 40nM efavirenz, sufficient to reduce reverse transcription to negligible levels in cellfree infection (see Figure 2D). Hence we measured the effect on cell viability of infection from the reverse transcription step onward, as this is what is compared in all our experiments.
Two independent experiments were performed in the cell line per parameter, with each experiment including cell sorting followed by PCR to obtain λ. The specific approach for each parameter was as follows:
Measurement of r (probability of infection per transmission)
1) r was measured using cellfree infection followed by flow cytometry 2 days postinfection to determine frequency of live infected cells P_{λ}, as measured by the fraction of live GFP positive cells normalized by the maximum fraction of cells which can become GFP positive upon infection in our system (70%).
2) λ, the number of viral transmissions per cell, was defined as the average number of detected HIV DNA copies per cell both integrated and unintegrated (see response to reviewer 3 below), as detected by splitting a single cell over multiple wells and performing PCR for HIV DNA.
3) Given that live cells were selected, r = 1(1P_{λ})^{1/λ}
Measurement of q (probability of cell death per transmission)
1) The fraction of live cells remaining postreverse transcription, defined as L_{λ} =[live cells with infection]/[live cells when reverse transcription blocked], was measured with cellfree infection by determining the concentration of live cells using cell counting with trypan blue exclusion. To prevent cell numbers changing due to cell division, proliferation was blocked by serum starvation.
2) λ was defined and measured as with r.
3) q = 1L_{λ}^{1/λ}.
Results of experiments are shown in the new Figure 2—figure supplement 4.
The results are also summarized in Supplementary file 2.
We used the measured values to calculate P_{λ} using Equation 2. The results are presented in Figure 2E, green line.
The measured q value is higher than our fit (0.14 versus 0.08) and the r value lower (0.24 versus 0.40). However, the pattern of the peak is not strongly affected. What this demonstrates is that the measured r and q are in the correct range in the cell line system to explain the observed peak.
This result is now presented in the main text in the Results section:
“In order to determine whether the fitted r and q values were within a reasonable range, we measured these values experimentally. […] The solution to Equation 2 using these values was very similar to the solution with the fitted values, indicating that the fitted values gave a reasonable approximation of the behavior of the system (Figure 2E, green dotted line).”
The methods to measure r and q are presented in the Materials and methods section:
“Measurement of r
Cellfree supernatant used in infection was derived as follows: 10^{6} cells/ml were infected with 2x10^{8} NL43 viral copies/ml (~20ng p24 equivalent) for 2 days. […] The remaining cells from the same infection were used to determine frequency of DRAQ7 negative, GFP positive cells 2 days postinfection using flow cytometry.”
“Measurement of q
Cellfree supernatant used in infection was derived as for r, except that one day before harvesting of the viral supernatant from infected cells, infected cells were washed twice with PBS and serumfree growth medium added. […] The concentration of live cells was measured using the TC20^{TM} automated cell counter (BioRad) with trypan blue staining (Lonza).”
To measure r and q in lymph node cells is not practical. This is because the infection is of a subset of lymph node cells (~10%). Therefore, measurement of λ where most cells are uninfected would require months of running gels. Furthermore, it is very unlikely that we could measure q using the approach above, since a given frequency of cell death in the infectable cell subpopulation would be diluted to a much lower frequency over the whole cell population, and unlikely to be detected above the high background lymph node cell death.
Comments from Reviewers 2 and 3 are attached.
Reviewer #2:
[…] The overall quality of the data could be improved by showing additional controls and increasing the numbers of cells analyzed to enhance the confidence in the copy number estimates.
Additional controls were added as detailed in the replies to the reviewer’s specific comments. The number of cells analyzed increased as follows:
For ACH2 cells, in total an additional 95 cells were analyzed, bringing to total analyzed cells to 166. This is now shown in Figure 2—figure supplement 1.
For the cell line infected with cellfree virus, an additional 12 cells were analyzed, bringing the total analyzed cells to 30. For coculture infection, an additional 26 cells were added, bringing the total analyzed cells to 60. These results are presented in Figure 2C with cellfree data now included as an inset.
With the new results, the mean number of corrected HIV DNA copies per cell was 17.
On the cell death front, it would be helpful for their hypothesis to show a time course with the cell death depicted over time.
We now include timelapse microscopy data showing the depletion of target cells in the absence and presence of efavirenz. The experiments show the total number of live and live infected cells in the optical dish with no drug, and in the presence of efavirenz. These experiments were done before we finalized experimental conditions, so the input of wild type virus infected donor cells (20%) and the corresponding efavirenz concentrations (30nM) were higher. However, the data clearly shows that extensive death due to infection occurs in the absence of drug, and is prevented by efavirenz. The data is presented as Figure 2—figure supplement 5.
This is now described in the Results section:
“In order to investigate the dynamics of cell depletion due to celltocell HIV spread and its modulation by addition of inhibitor, we performed timelapse microscopy over a two day infection window. […] The deterioration in live cell numbers was averted by the addition of EFV (Figure 2—figure supplement 5).”
The method used is described in the Materials and methods:
“Timelapse microscopy and image analysis
For imaging infection by timelapse microscopy, cell density was reduced to 5x10^{4} cells/ml and cells were attached to ployllysine (ΣAldrich) coated optical six well plates (MatTek). […] The number of mCherry positive 16 pixel^{2} squares around the cell centers was used as the as the number of total target cells at each timepoint, and the number of squares double positive for fluorescence in the GFP channel was used as the number of infected target cells.”
Lastly, an important detail not mentioned, is whether the use of lymph node cells is critical or not for the phenotypes that they describe.
In this study we concentrated on lymph nodes, which were reported to be sites of lower drug penetration as well as efficient HIV replication. However, we did try PBMCs, and observed no obvious fitness optimum, as reviewer 2 suggests may be a possibility since PBMCs are less prone to cell death relative to lymph node cells. The first set of 3 independent experiments using blood from one healthy donor showed a small peak at a very low EFV concentration. We then focused on the lower EFV concentration range but could not reproduce the peak in other donors. These results are presented as Figure 5—figure supplement 1.
It is described in the Results section:
“Peaks in the number of live infected cells in the face of drug may be specific to lymph node derived cells. Celltocell infection of peripheral blood mononuclear cells (PBMC) with wild type HIV showed a slight peak at a very low EFV concentration in cells from one blood donor, which was not repeated in cells from two other donors (Figure 5—figure supplement 1).”
We add this observation to the Discussion:
“Physiologically, an infection optimum in the face of an antiretroviral drug may be important in lymph node cell HIV infection and may be less pronounced in infectable cells from peripheral blood.”
The study may have implications for the establishment of viral reservoirs in the context of poorly controlled infection or infections with some degree of drug resistance.
We thank the reviewer for this concise summary and add it to the Discussion:
“This study may have implications for the establishment of viral reservoirs in the context of poorly controlled infections, infections with some degree of drug resistance, or infections where some replication may take place in the face of ART, since infected cell survival is a prerequisite for long term persistence.”
1) The authors use the ACH2 cell line to test the ability of their assay to detect single copy integrations. Why do they only spread each cell over 4 wells in control studies as opposed to the 10 wells used in their experimental studies? The rationale for using different dilution schemes should be explained. A comparison of the efficiency of PCR at the different dilutions would also be informative.
There is no reason for the difference except that we started with 10 well spreads in cell line coculture infections and realized that we can use less reagents to obtain the same data. We now state this in the Results section:
“Using fewer wells saved reagents without changing sensitivity, as demonstrated in the ACH2 cell line (Figure 2—figure supplement 1C).”
We compared efficiency of HIV copy detection in ACH2 cells between 4 and 10 well spreads and did not find a substantial difference, as the frequency in the 10well split was 0.27 while the parental was 0.24. The data is presented in Figure 2—figure supplement 1,
We describe this in the Results section:
“Similar results were obtained when the ACH2 cell line was subcloned or split over 10 wells (Figure 2—figure supplement 1C).”
The ACH2 cell line yields fewer than 1 copy per cell, which may be a limitation of the PCR assay, but also could be reflective of heterogeneity in a cell line that is assumed to be clonal and a uniform karyotype. It is not described how recently the ACH2 cell line that they are using has been cloned. Repeating the ACH2 studies with a subclone would be beneficial.
We subcloned the ACH2 cells as described in the Materials and methods:
“Subcloning of ACH2 cells
Cells from the parental ACH2 cell line were diluted to 10 cells/ml in conditioned medium, with conditioned medium generated by culturing ACH2 cells to 10^{6} cells/ml, then filtering through a 0.22µm filter (Corning). […] Clones were detected in 5% of wells and two clones, designated D6 and C3, were randomly chosen and further expanded.”
We checked two subclones, with 56 cells measured in the first subclone and 24 in the second. Frequency of detection in each was 0.21, which we judged was not substantial relative to the parental line. We describe the results in Figure 2—figure supplement 1C above.
Figure 2 would be more informative if it showed the results of both cellfree virus infection versus coculture infection. The numbers of cells analyzed in Figure 2C (n=34) does not appear sufficient to provide a robust sense of the distribution of the DNAs in the infected cells. Is this a bimodal distribution? Histogram comparison of cellfree should also be used for comparison.
We now include the cellfree distribution as an inset in Figure 2C. Frequency of HIV DNA copies was 0.23, as described in the response to general comment 3. We describe this in the Results section:
“We also detected the HIV copy number in 30 GFP positive cells infected by cellfree HIV. In this case, we detected either zero or one HIV copy per cell (Figure 2C inset). The frequency of single HIV DNA copies was 0.23, identical to the measured result in the ACH2 cell line.”
We checked the fit of the distribution resulting from coculture infection, which we expanded to 60 analyzed cells. While a distribution with two means fit the data better, we could not rule out that this was due to added fit parameters since the fitted distribution itself did not show peaks. Most likely, this means that the peak at 5 raw copies and no cells at 6 copies is a chance event that gives an impression of bimodality. The results are presented in Figure 2—figure supplement 2.
These results are described in the Results section:
“We assayed 60 cells and obtained a wide distribution of viral DNA copies per cell, which ranged from 0 to 9 copies (Figure 2C). […] The frequency of single HIV DNA copies was 0.23, identical to the measured result in the ACH2 cell line.”
2) Figure 4.
What happens with partial antibody inhibition of cellfree infection? These studies should be performed with infection from cellfree virus. To further test their hypothesis partial inhibition of cellfree infection should only decrease infected cells.
We now include the results for antibody inhibition of cellfree infection as Figure 4—figure supplement 1. As can be seen, antibody only decreases infected cells in cellfree infection.
This is described in the Results section:
“In contrast, cellfree infection in the face of b12 showed a sharp and monotonic drop in live infected cells for both wild type and mutant virus (Figure 4—figure supplement 1).”
3) Figure 5.
As controls for cell death measurements, the viability of the uninfected control lymph node cells, treated and untreated, should also be illustrated. In the literature lymph node cells (mostly tonsilar) are very prone to cell death, and it is important to understand to what extent they are measuring virusinduced cell death.
We measured cell death in uninfected lymph node cells, untreated and treated with efavirenz. We did not find that efavirenz alone results in a peak of live cells (Figure 6—figure supplement 3).
This is discussed in the Results section:
“To examine if the observed peak in live cells may be due to EFV alone, we measured cell viability in lymph node cells from one of the study participants used in the above experiment as a function of EFV without infection. No dependence on EFV in the absence of infection was detected (Figure 6—figure supplement 3).”
In addition, a time course of the cell death observed in the infections may would also be helpful to evaluate their hypothesis that the increase of infected cells in the partially treated cells is due to increased infected cell survival – i.e. decreased death.
As discussed in the response to the reviewer’s general comment 3, we now include timelapse measurements of infection treated and untreated with the efavirenz which supports our hypothesis that the increase of infected cells in the partially treated cells is due to increased infected cell survival.
4) Are the phenotypes described in the primary lymph node cells also observed for peripheral blood lymphocytes? Or are there differences between the peripheral blood versus the lymph node T cells.
As discussed in the response to the reviewer’s general comment 4, the work we did perform on PBMCs showed little increase in infection with partial inhibition, consistent with the reviewer’s intuition that this may be a lymph node specific effect.
Reviewer #3:
This is an interesting and worthwhile paper that examines the effects of incomplete inhibition of HIV infection with reverse transcription inhibitors. The experimental results presented are in accordance with model predictions. However, the presentation of the model should be improved as detailed below. Also, I found the labeling and numbering of the supplemental figures confusing. I have no substantial concerns and this the paper should be published after minor revision.
We thank the reviewer for the support. We have clarified the model parameters and how they are measured as detailed in the responses to the Reviewer’s specific comments below to address this very important point. We have ensured that all supplemental figures are labelled in eLife format and sequentially referenced in the text.
1) The terms and concepts in the paper are not clearly defined. In the first paragraph of the Results, you mention each donor to target transmission. It is not clear what a transmission is. Does it refer to viral entry of either a free virus or a virus (or genome) by celltocell transmission? It needs to be defined. Second define what you mean by infection. Does a cell have to produce virus to be considered infected? Does the virus have to integrate or is it sufficient to simply reverse transcribe? Is a latently infected cell an infected cell in your model? Also, when you say drug therapy increases the number of live infected cells do you mean live productively infected cells, live HIV DNA^{+} cells, etc.
The reviewer has identified an important gap in our presentation of the results which we have now corrected. We now present a table of definitions for each parameter in our system and how it is defined and measured. This now appears in the manuscript as Supplementary file 1.
We make extensive use of the sharpened definitions, including in the Results section:
“We introduce a model of infection where each donor to target transmission leads to an infection probability r and death probability q per infection attempt, where each model parameter is defined in Supplementary file […] We define L_{λ} as the probability of a cell to survive infection given λ infection attempts. Assuming infection attempts act independently, L_{λ}=(1q)^{λ}.”
2) The assumption that all transmissions have equal probabilities to infect target cells seems to ignore the possibility that some virions carry defective genomes while others do not. This does not seem realistic. In your experimental system is the ratio of HIV RNA to TCID50 close to one so you can ignore defective particles?
Given that transmission is defined as the mean number of HIV DNA copies per cell and that this is what we measure, any virion which does not successfully reverse transcribe will not be detected. Therefore, the measurement of the number of transmissions per cell λ excludes virions defective in virus cycle stages leading up to an including reverse transcription. The fraction of defective viruses postreverse transcription is captured in the term r, as HIV DNA copies with mutations which make them defective postreverse transcription would not be able to produce viral proteins, our measure of successful infection. This is now described in the Results section:
“r and q capture the probabilities for a cell to be infected or die postreverse transcription. For example, mutations which reduce viral fitness by decreasing the probability of HIV to integrate would reduce r, while mutations which reduce the probability of successful reverse transcription would reduce λ.”
3) You assume productive infection and death are independent events. While the events may be independent the probability of death is certainly not independent of whether a cell is productively infected or not. Further the probability of death is time dependent. The probability a cell dies one hour after viral entry (infection?) is clearly quite different than the probability it dies days after infection. You may want to define q as the probability a cell has died by time t after infection and the same for Pλ where t is 2 days or 4 days for your various experiments.
We thank the reviewer for bringing up the important point that cell death is a timedependent process and the timepoint at which it is quantified needs to be defined. We therefore corrected our definitions to state that the fraction of live infected cells is measured at 2 days postinfection for the cell lines and 4 days postinfection for lymph node cells, as the reviewer suggested. This now appears in Supplementary file 1 as discussed above.
4) Results, second paragraph. You claim antiretroviral drugs reduce the number of infecting virions. This implies that by infecting you must mean the virus reverse transcribes. In standard viral dynamic models the effects of ART are to reduce the infection probability, i.e., r in your model not λ. Thus it is important to clarify your definitions.
The reviewer is correct that you would expect the effect of the reverse transcriptase inhibitor to be included in r. However, our measurable transmission λ is a reverse transcribed copy of HIV DNA. Hence in our model r accounts for the probability of a reverse transcribed viral genome successfully integrating into the host genome and producing viral proteins, which are independent of the action of efavirenz. The way we set up our model, the action of the reverse transcriptase inhibitor is on λ, not r. We address this as described in the response to point 2 of the reviewer’s comments.
5) You introduce the drug effect as a constant d in the second paragraph of the Results. Later in the paper you make reference to IC50 and Hill coefficients for the drug. These need to be tied together in an explicit manner. In viral dynamic modeling the effectiveness of a drug, epsilon, (eps) is introduced where eps=1 is a 100% effective drug, e.g. stops all reverse transcription, and where eps=0 means the drug has no effect. Then a pharmacodynamic model is used to relate eps with the drug concentration, C, e.g. eps= C^n/(IC50^n + C^n). Thus rather than λ/d you could model the drug effect by (1eps)λ and use the above formula to link the effect to drug concentration. See for example, Canini and Perelson, 2014. Equation 3 in your manuscript accomplished the same, but the above seems more straightforward.
The reviewer’s notation is more elegant than what we use since it directly introduces the effect of drug per virus into the equation. The reason we would argue to keep our current notation is that d is on the same scale as λ, which we think makes a ratio such as d/λ more intuitive. Thus d=λ would reduce the infection to single HIV DNA copies per cell. We have stated the conversion between d and ε in the Materials and methods section:
“we used the relation for the fraction cells remaining infected in the face of drug (Canini and Perelson, 2014), whose definition is equivalent to Tx at λ<1:
1d=1=1EFVhEFVh+IC50h (3).”
We also discuss the conversion in the Results section:
“This is equivalent to 1ε in a commonly used model describing the effect of inhibitors on infection, where ε is drug effectiveness, with the 50% inhibitory drug concentration (IC_{50}) and the Hill coefficient for drug action as parameter values (Canini and Perelson, 2014; Shen et al., 2008).”
6) Results, fourth paragraph. When you measure the number of reverse transcribed copies of viral DNA I assume you are measuring both integrated and unintegrated DNA – please state this explicitly as later you seem to imply you are measuring integrated DNA, e.g. Results, sixth paragraph.
We thank the reviewer for pointing out this omission and we now explicitly state that λ measures the number of HIV DNA copies, integrated and unintegrated. This is now reflected in Supplementary file 1 as described above, as well as in the Results section:
“In our experimental system, one infection attempt is measured as one HIV DNA copy, whether integrated or unintegrated.”
7) The current work focuses only on reverse transcriptase inhibitors. Model therapy uses combination therapy including integrase inhibitors and protease inhibitors. Expanding the discussion of the clinical implications of this work to include combination therapy would be worthwhile.
We agree with the reviewer that our analysis can be applied beyond reverse transcriptase inhibitors. We used efavirenz in our study due it being a common principal component of current first line antiretroviral therapy, with common drug resistance mutations. However, the infection optimum we describe should occur with other classes of antiretroviral drugs, since they all should decrease the multiplicity of infection between cells. We now include this point in the Discussion:
“We used EFV in our study since it is a common component of first line antiretroviral therapy, with frequent drug resistance mutations. […]The more complex outcome of partial inhibition of infection should also be considered in other infections where multiple pathogens infect one cell and host cell death is a possible outcome (Mahamed et al., 2017).”
[Editors' note: further revisions were requested prior to acceptance, as described below.]
We want you to address the following comments. Note that we are not asking you to do further experiments. For the sake of recordkeeping, we are enclosing all comments.
Reviewer #1:
I am satisfied with their measurements.
In general, I suggest breaking up long sentences (especially the one in Abstract, also the end of the third paragraph in the Introduction) for readability. Sentences that last more than three lines are generally more difficult to read than shorter sentences.
We have now broken up long sentences as follows.
Abstract:
“Using a computational approach, we found that partial inhibition of transmissions of multiple virions per cell could lead to increased numbers of live infected cells. If the number of viral DNA copies remains above one after inhibition, then eliminating the surplus viral copies reduces cell death.”
Introduction:
“Multiple infections per cell have been reported in celltocell spread of HIV (Baxter et al., 2014; Boull et al., 2016; Dang et al., 2004; Del Portillo et al., 2011; Dixit Perelson, 2004; Duncan et al., 2013; Law et al., 2016; Reh et al., 2015; Russell et al., 2013; Sigal et al., 2011; Zhong et al., 2013). In this mode of HIV transmission, an interaction between the infected donor cell and the uninfected target results in directed transmission of large numbers of virions (Baxter et al., 2014; Groppelli et al., 2015; Hubner et al., 2009; Sowinski et al., 2008).”
Introduction:
“Both modes occur simultaneously when infected donor cells are cocultured with targets. However, the celltocell route is thought to be the main cause of multiple infections per cell (Hubner et al., 2009).”
Introduction:
“However, more recent work investigating markers associated with HIV latency in the face of ART found that the average number of HIV DNA copies per cell is greater than one in 3 out of 12 individuals. This occurred in the face of ART in the CD3 positive, CD32a high CD4 T cell subset (Descours et al., 2017).”
Introduction:
“We observed that partially inhibiting infection with drug or antibody resulted in an increase in the number of live infected cells in both a cell line and in lymph node cells. This is, to our knowledge, the first experimental demonstration at the cellular level that attenuation of HIV infection can result in an increase in live infected cells under specific infection conditions.”
Results:
“To estimate this, we used PCR to detect the number of reverse transcribed copies of viral DNA in the cell by splitting each individual infected cell over multiple wells. We then detected the number of wells with HIV DNA by PCR amplification of the reverse transcriptase gene.”
Results:
“This is equivalent to 1ε in a commonly used model describing the effect of inhibitors on infection. In this model, is drug effectiveness, with the 50 percent inhibitory drug concentration (IC_{50}) and the Hill coefficient for drug action as parameter values (Canini Perelson, 2014; Shen et al., 2008).”
Discussion:
“If each HIV DNA copy increases the probability of cell death, reducing the number of HIV DNA copies without eliminating infection should lead to an increased probability of infected cell survival. This would consequently lead to an increase in the number of live infected cells.”
Reviewer #3:
In this revised version on the paper the authors have clarified their definitions by incorporating a supplemental table and additional text describing their model. However, this clarification has raised the following issues that need to be resolved:
1) As defined in the Results, λ is the number of infection attempts, where one attempt is measured as one HIV DNA copy. Thus, λ is an integer. The formula, Equation 1, makes sense when λ is an integer. However, what is determined experimentally is not the integer number of DNA copies but rather the mean number of copies per cell, which is not necessarily an integer. To make sense Equation 1 should be reformulated as a conditional probability of a cell being infected and alive given x attempts, call this P(I  x). What I believe the authors want is the probability of a cell being infected and alive, which is then P(I) =Sum P(I x) (p(x)), where p(x) is the probability of having x infection attempts, and the sum is over all nonnegative values of x. One could then assume, for example, that the number of attempts is Poisson distributed with mean λ. Then the resulting formula for P(I) would involve the mean value of the number of attempts. A better choice would be the empirical distribution that the authors measured. Note that what is actually measured is the number of infected cells that are alive at a certain time, i.e. P(I), and the mean number of attempts, what I called λ above.
We thank the reviewer for this comment. The reviewer is correct that our previous Equation 1 applies to infection attempts per cell, which are integer values. We have now reformulated the model to use the mean number of infection attempts, as described below and in the new Appendix 1: Supplementary Mathematical Analysis section. We used the Poisson distribution, the first suggestion of the reviewer, as we believe we need much more extensive experimental data to determine the exact empirical distribution.
We now consider the probability of a cell being infected and live given n HIV DNA copies, where n is an integer (Appendix 1Equation 1).
If P(n) is Poisson with mean λ, the fraction of cells that are productively infected is as in Appendix 1Equation 2.
We introduce a drug strength value d, where d = 1 in the absence of drug and d > 1 in the presence of drug. In the presence of drug, λ is decreased to λ/d. The drug therefore tunes λ. The probability of a cell to be infected and live given drug strength d is therefore as in Appendix 1Equation 3.
Appdenxi 1Equation 3 is now used to derive Figure 1, which is qualitatively similar to Figure 1 from the previous submission using Equation 2. The range of q values was modified to show the full range of infection behavior with the new model (Figure 1)
Experimentally, we measure P_{λ}, the number of infected cells, and L_{λ}, the number of live cells remaining two days post infection for the cell line.
Assuming a Poisson distribution of the number of infection attempts, the probability of being productively infected is as in Appendix 1Equation 4.
Similarly, the probability of a cell to live is as in Appendix 1Equation 5.
The derivation of the model above is included in the new Supplementary Mathematical Analysis section. The new definition of P_{λ} and L_{λ} is now updated in Supplementary file 1.
The experimentally determined values for r and q are now recalculated to be r = 0.28 and q = 0.15. This is now updated in Supplementary file 2.
Comparing this model to the model described by the previous model (Equation 2 in the previous submission), we note that the models give very similar results at λ > 1, but that our previous model underestimates the rate of decline in infection at λ < 1, as shown in Author response image 1 for L100I HIV mutant coculture infection of the cell line, where the effect was most pronounced.
To investigate whether the experimental results were above the fit line due to technical factors, we refined our experimental approach to exclude donor cell – target cell fusions, which may give a baseline in the infection. In the previous submission, we judged the fraction of fusions to be negligible and selected on target cells based on mCherry fluorescence. Here, we used the vital stain Cell Trace Far Red (CTFR) to label donor cells. Target cells still expressed mCherry. We therefore selected target cells based on mCherry, but excluded any target cells which were CTFR and mCherry double positive. Due to the extra fluorescence channel, it was necessary to change from our previous death detection dye with farred emission (DRAQ7) to the blue spectrum emission death detection dye (DAPI). For this, we used a different flow cytometer (FACSAria Fusion) which had the 355nm laser for DAPI excitation. The gating strategy is now illustrated in Figure 2— Figure supplement 4 for wild type HIV infection.
It is described in the Results section: “To obtain the number of infected target cells, and specifically exclude donor cells or donortarget cell fusions, target cells were marked by the expression of mCherry. […] Donortarget cell fusions were excluded by excluding CTFR positive cells (see Figure 2— figure supplement 4 for gating strategy).” The results shown in Figure 2— figure supplement 4 are described in the Results section:
“While the percent of infected cells was reduced with drug, the concentration of live infected cells increased (Figure 2—figure supplement 4).”
We also describe the method in the Materials and methods section:
“To determine the number of live infected cells in reporter cell line experiments, E7 RevCEM reporter cells were infected as above used as donor cells. […] The coculture infection was pulsed with 100ng/ml DAPI (Σ) immediately before flow cytometry and the number of live infected targets cells was determined by the number of DAPI negative, CTFR negative and mCherry and GFP double positive cells on a FACSAria Fusion machine (BD Biosciences) using the 355, 488 and 633nm laser lines.”
We used the same approach for infection with the L100I HIV mutant, and the gating strategy is illustrated in the new Figure 3—figure supplement 2.
We describe the new figure in the Results:
“We next performed coculture infection (see Figure 3—figure supplement 2 for gating strategy).”
We also increased the number of data points for wild type and mutant infections in the cell line when we repeated the experiments. The number of data points increased from 5 to 7 for wild type HIV infection, and 6 to 8 for L100I mutant infection of the cell line (3 independent experiments each). For mutant infection, we concentrated on drug concentrations closer to the peak since that was the region with the most pronounced decline in infection. Excluding donortarget cell fusions using new gating strategy resulted in an incremental improvement to the data, as shown in Author response image 2 for wild type mutant HIV infection in the cell line:
These new results showed that the lower than expected decline in infection at high drug concentrations was not a technical artifact due to a baseline of infected donor cells which were fused to targets. This deviation from the model at high drug concentrations away from the peak infection with drug is now noted in the Results section:
“We note that both wild type and mutant coculture infection has data points above the fit line at the highest drug concentrations. This may be a limitation of our model at drug values much higher than observed at the infection optimum. In this range of drug values, our model predicts a more pronounced decline in the number of infected cells than is observed experimentally.”
We also include the observation in the Discussion section:
“Neither model accurately captures infection dynamics at high drug concentrations, away from the infection optimum. In this range, where the number of infection attempts per cell is much lower than 1, infection declined more slowly with drug than predicted.”
We feel this observation may be potentially important to understanding the HIV reservoir. However, the existence of a fitness optimum with drug is a separate issue, and our current model successfully describes how the fitness optimum is formed.
An important conclusion reached when we compared the fits before and after fusion exclusion, as well as the fit of the previous iteration of L100I mutant infection to the current experiments, is that small differences in the experimental data resulted in large differences in the value of r. Hence, r did not seem to be well constrained in the fits. As shown in Figure 1, rdetermines the amplitude of P_{λ}. This is now emphasized in the Results section:
“Hence, the value of r strongly influences the amplitude of P_{λ}.”
To better constrain r as a result of these new observations, we switched from normalizing the number of infected cells and P_{λ} to the maximum value in each case to using the number of live infected cells per ml. For this, the P_{λ} predicted by the model was multiplied by the number of target cells at the start of infection, which was 10^{6} cells/ml. The best fit for the model as well as the result of the calculation with the experimentally determined r and q is shown in updated Figures 2 and 3 as Figure 2D (wild type infection) and Figure 3B (L100I mutant infection).
The great improvement of this approach is that r values are now constrained: r = 0.22 and q = 0.17 for wild type HIV cell line infection and r = 0.29 and q = 0.13 for L100I mutant cell line infection. This now compares well with the experimentally determined r = 0.28 and q = 0.15. The new analysis is described for wild type infection in the Results section:
“While the percent of infected cells was reduced with drug, the concentration of live infected cells increased (Figure 2—figure supplement 4). […] Equation 3 best fit the behaviour of infection when r=0.22 and q=0.17, resulting in a peak at 4.8nM EFV (Figure 2D, black line).”
It is described for cell line infection with mutant L100I virus in the Results section:
“The fits recapitulated the experimental results when r = 0.29 and q = 0.13, with a fitted peak at 45nM EFV (Figure 3B, black line). The solution to Equation 3 using the measured values for r and q showed a similar pattern to that obtained with the fitted values (Figure 3B, dashed green line).”
We applied the new model to the coculture lymph node infection results. However, we could not use the number of live infected cells to constrain r in this case, since we do not know the number of infectable targets. This is because the lymph node is a complex environment containing, besides CD4^{+} T cells of various kinds and potentially infectable APCs, also CD8^{+} T cells and B cells. This is now discussed in the Results section:
“We did not calculate the predicted number of infected cells for P_{λ}/d values since the lymph node is a complex environment containing different cell subsets (Sallusto et al., 1999) and the number of infectable target cells at the start of infection is difficult to determine. Hence, we normalized both the experimental number of live infected cells and the P_{λ}/d values from Equation 3 to the maximum value in each case.”
The results of the new fit are presented in Figure 6B.
For this fit, r = 0.91 and q = 0.15. While q matches the measured cell line result well, the r is high. We believe this because it is poorly constrained in the fit, as discussed in the Results section:
“The fits recapitulated the experimental results when r = 0.91 and q = 0.15, with a fitted peak at 90nM EFV (Figure 6B, black line). […] However, the fitted r value in this case is not expected to be accurate since we were unable to constrain it with the number of infected cells relative to the starting number of target cells.”
2) According to the model q should be independent of the number of attempts. Can this be tested by using different drug concentrations, which would vary the mean number of attempts? If you already have experimental data on this, then include them. Otherwise, just discuss it.
We thank the reviewer for the insightful comments and agree that the probability to die per infection attempt may depend on how many previous infection attempts occurred in the same cell. We discuss this point in the Discussion section:
“Moreover, the probability of death per HIV DNA copy we denote q may be dependent on how many infection attempts preceded the current infection attempt, and the model can be improved by measuring this dependence.”
3) In the original manuscript where Equation 1 was derived the authors assumed q was independent of r. In the second paragraph of the Results the authors say q is the probability of a cell to die post reverse transcription, which seems more sensible to me. Under the original assumption of q and r being independent, a cell in which the attempt was unsuccessful would have the same probability of dying as a cell in which reverse transcription succeeded. With q now defined as in the second paragraph of the Results, the fundamental formula given by Equation 1 may not be correct and deserves a more thorough discussion. I think the fundamental process being described is that at each attempt a cell is either infected with probability r, or not infected with probability 1r. If it is infected, i.e. had a successful reverse transcription, then with probability q it remains alive and with probability 1q it dies. Assume there are x attempts. Let z be the number of cells infected after x attempts. Then z is binomially distributed, i.e. z= Bin (x,r). Further, if we are interested in the number of cells infected and alive after x attempts, then P(Ix) = Sum_z=1 to x (x choose r) r^z (1r)^(xz)(1q)^z, where the factor (1q)^z is the probability of a cell surviving after each of the z infections. Simplifying, P(Ix) = [r(1q) + (1r)]^x, which is not the same as Equation 1. The number of attempts, x, is again random and as above one can convert this conditional probability into P(I) by assuming a Poisson (or some other distribution such as the empirically measured one) for the number of attempts p(x).
We thank the reviewer for the comment which essentially questions whether the probability of death needs to be independent of the probability to infect for an infection optimum with drug to exist. We have now constructed a second model where death is dependent on infection.
In the case where cell death depends on productive infection, the probability of productive infection and survival of the cell given n attempts is shown in Appendix 1Equation 6.
Here, r^{0}and q^{0}are the probabilities for infection and death respectively. Their relationship to the experimentally measured r and q are discussed in the next section. r^{0m}(1 − q^{0})^{m} is the probability that m infection attempts are successful and none of the attempts triggers cell death, while (1 − r^{0})^{n}^{−m }accounts for n − m unsuccessful infections.
Realizing that (Appendix 1Equation 7) we can simplify to Appendix 1Equation 8.
Averaging this over a Poisson distribution with mean λ, we find Appendix 1Equation 9.
The two models are the same if: (Appendix 1Equation 10).
In this case the two models are a reparametrization of the same process. The experimental measurement of r and q, as described in the main manuscript, consists of estimating r as the fraction of infected cells when the fraction of dead cells is low due to a low number of measured DNA copies per cell. q is derived by measuring the number of live cells remaining when the number of HIV DNA copies per cell is known. If the probability of death is dependent on infection, then death is only of infected cells and is a product of the probability to be infected and die (r^{0}q^{0}). Likewise, the underlying probability of infection r^{0}will account for the fraction dead cells (q) and any cells which are infected but not dead (r(1 − q)). We therefore conclude that independence of cell death and successful infection is not a required condition for a fitness optimum. This is now stated in the Results section:
“Our model assumes that cellular infection and death due to an HIV infection attempt are independent processes. […] The models are equivalent, showing that independence of cell death and infection is not a necessary condition for an infection optimum to occur in the presence of inhibitor.”
We also discuss the result in the Discussion section:
“Construction of the model assumed independence of productive infection and cell death. However, as shown in the Supplementary Mathematical Analysis section, an equivalent model can be constructed assuming a dependence of cell death on infection.”
The derivation of the above model is included in the new Supplementary Mathematical Analysis section.
4) In Supplementary file 1 the authors say the way they measure the probability a cell is alive after λ attempts is by computing the concentration of cells alive with λ attempts divided by the concentration of cells alive with no transmissions. However, since the concentration of cells alive with no transmissions can approach zero as time goes on it is clear that this fraction need not be less than 1. Also, as noted above the authors need to use conditional probabilities and derive formulas that involve the mean number of attempts. I would suggest they compute the ratio of the probability a cell is infected to the probability a cell in not infected and alive, which is what they measure, and see if they can derive a formula for q.
If we understand the reviewer’s comment correctly, the reviewer’s concern is that the concentration of cells alive with no transmissions approaches zero because both are measured simultaneously in the same culture. However, we set up the experiment differently, where cells equally exposed to all infection conditions except HIV DNA copies were compared for cell death. This was accomplished by treating one of two separate infected cell cultures with EFV at a concentration which resulted in almost complete inhibition of cellfree infection, the mode of infection used. We have now clarified this in the caption of Figure 2—figure supplement 5. Experimental measurement of r and q:
“The concentration of live cells two days post cellfree infection was measured using trypan blue exclusion and divided by the concentration of live reporter cells in a separate cell culture infected with the identical amount of virus, but with infection was inhibited with 40nM EFV.”
We have also included a clarification in Supplementary file 1, which states how the fraction of live cells remaining L_{λ} is measured.
https://doi.org/10.7554/eLife.30134.040Article and author information
Author details
Funding
Human Frontier Science Program (CDA 00050/2013)
 Alex Sigal
European Research Council (Stg. 260686)
 Richard A Neher
DELTAS Africa Initiative (Graduate Fellowship (DEL15–006))
 Isabella Markham Ferreira
National Research Foundation (Graduate Fellowship)
 Laurelle Jackson
 Jessica Hunter
National Institutes of Health (R21MH104220)
 Alex Sigal
Poliomyelitis Research Foundation (Graduate Fellowship)
 Isabella Markham Ferreira
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
This work was supported by National Institutes of Health Grant R21MH104220. AS was supported by a Human Frontiers Science Program Career Development Award CDA 00050/2013. RAN is supported by the European Research Council through grant Stg. 260686. LJ and JH are supported by a fellowship from the South African National Research Foundation. IMF is supported through a SubSaharan African Network for TB/HIV Research Excellence (SANTHE, a DELTAS Africa Initiative (grant #DEL15–006)) fellowship, and a Poliomyelitis Research Foundation fellowship 17/59. Open access publication of this article has been made possible through support from the Victor Daitz Information Gateway, an initiative of the Victor Daitz Foundation and the University of KwaZuluNatal.
Ethics
Human subjects: Lymph nodes were obtained from the field of surgery of participants undergoing surgery for diagnostic purposes and/or complications of inflammatory lung disease. Informed consent was obtained from each participant, and the study protocol approved by the University of KwaZuluNatal Institutional Review Board (approval BE024/09).
Reviewing Editor
 Wenying Shou, Fred Hutchinson Cancer Research Center, United States
Publication history
 Received: July 18, 2017
 Accepted: March 8, 2018
 Accepted Manuscript published: March 20, 2018 (version 1)
 Version of Record published: April 10, 2018 (version 2)
Copyright
© 2018, Jackson et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics

 1,603
 Page views

 209
 Downloads

 6
 Citations
Article citation count generated by polling the highest count across the following sources: PubMed Central, Crossref, Scopus.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading

 Computational and Systems Biology
Improving muscle function has great potential to improve the quality of life. To identify novel regulators of skeletal muscle metabolism and function, we performed a proteomic analysis of gastrocnemius muscle from 73 genetically distinct inbred mouse strains, and integrated the data with previously acquired genomics and >300 molecular/phenotypic traits via quantitative trait loci mapping and correlation network analysis. These data identified thousands of associations between protein abundance and phenotypes and can be accessed online (https://muscle.coffeeprot.com/) to identify regulators of muscle function. We used this resource to prioritize targets for a functional genomic screen in human bioengineered skeletal muscle. This identified several negative regulators of muscle function including UFC1, an E2 ligase for protein UFMylation. We show UFMylation is upregulated in a mouse model of amyotrophic lateral sclerosis, a disease that involves muscle atrophy. Furthermore, in vivo knockdown of UFMylation increased contraction force, implicating its role as a negative regulator of skeletal muscle function.

 Computational and Systems Biology
Intracellular states probed by gene expression profiles and metabolic activities are intrinsically noisy, causing phenotypic variations among cellular lineages. Understanding the adaptive and evolutionary roles of such variations requires clarifying their linkage to population growth rates. Extending a cell lineage statistics framework, here we show that a population’s growth rate can be expanded by the cumulants of a fitness landscape that characterize how fitness distributes in a population. The expansion enables quantifying the contribution of each cumulant, such as variance and skewness, to population growth. We introduce a function that contains all the essential information of cell lineage statistics, including mean lineage fitness and selection strength. We reveal a relation between fitness heterogeneity and population growth rate response to perturbation. We apply the framework to experimental cell lineage data from bacteria to mammalian cells, revealing distinct levels of growth rate gain from fitness heterogeneity across environments and organisms. Furthermore, third or higher order cumulants’ contributions are negligible under constant growth conditions but could be significant in regrowing processes from growtharrested conditions. We identify cellular populations in which selection leads to an increase of fitness variance among lineages in retrospective statistics compared to chronological statistics. The framework assumes no particular growth models or environmental conditions, and is thus applicable to various biological phenomena for which phenotypic heterogeneity and cellular proliferation are important.