Abstract
Celltocell viral infection, in which viruses spread through contact of infected cell with surrounding uninfected cells, has been considered as a critical mode of virus infection. However, since it is technically difficult to experimentally discriminate the two modes of viral infection, namely cellfree infection and celltocell infection, the quantitative information that underlies celltocell infection has yet to be elucidated, and its impact on virus spread remains unclear. To address this fundamental question in virology, we quantitatively analyzed the dynamics of celltocell and cellfree human immunodeficiency virus type 1 (HIV1) infections through experimentalmathematical investigation. Our analyses demonstrated that the celltocell infection mode accounts for approximately 60% of viral infection, and this infection mode shortens the generation time of viruses by 0.9 times and increases the viral fitness by 3.9 times. Our results suggest that even a complete block of the cellfree infection would provide only a limited impact on HIV1 spread.
https://doi.org/10.7554/eLife.08150.001eLife digest
Viruses such as HIV1 replicate by invading and hijacking cells, forcing the cells to make new copies of the virus. These copies then leave the cell and continue the infection by invading and hijacking new cells. There are two ways that viruses may move between cells, which are known as ‘cellfree’ and ‘celltocell’ infection. In cellfree infection, the virus is released into the fluid that surrounds cells and moves from there into the next cell. In celltocell infection the virus instead moves directly between cells across regions where the two cells make contact.
Previous research has suggested that celltocell infection is important for the spread of HIV1. However, it is not known how much the virus relies on this process, as it is technically challenging to perform experiments that prevent cellfree infection without also stopping celltocell infection.
Iwami, Takeuchi et al. have overcome this problem by combining experiments on laboratorygrown cells with a mathematical model that describes how the different infection methods affect the spread of HIV1. This revealed that the viruses spread using celltocell infection about 60% of the time, which agrees with results previously found by another group of researchers. Iwami, Takeuchi et al. also found that celltocell infection increases how quickly viruses can infect new cells and replicate inside them, and improves the fitness of the viruses.
The environment around cells in humans and other animals is different to that found around laboratorygrown cells, and so more research will be needed to check whether this difference affects which method of infection the virus uses. If the virus does spread in a similar way in the body, then blocking the cellfree method of infection would not greatly affect how well HIV1 is able to infect new cells. It may instead be more effective to develop HIV treatments that prevent celltocell infection by the virus.
https://doi.org/10.7554/eLife.08150.002Introduction
In in vitro cell cultures and in infected individuals, viruses may display two types of replication strategies: cellfree infection and celltocell infection (Sattentau, 2008; Martin and Sattentau, 2009; TalbertSlagle et al., 2014). Both transmission means require the assembly of infectious virus particles (Monel et al., 2012), which are released in the extracellular medium for cellfree transmission, or concentrated in the confined space of celltocell contacts between an infected cell and bystander target cells in the case of celltocell transmission. It has been shown that most enveloped virus species, including human immunodeficiency virus type 1 (HIV1), a causative agent of AIDS, spread via celltocell infection, and it is considered that the replication efficacy of celltocell infection is much higher than that of cellfree infection (Sattentau, 2008; Martin and Sattentau, 2009; TalbertSlagle et al., 2014). However, it is technically impossible to let viruses execute only celltocell infection. In addition, since these two infection processes occur in a synergistic (i.e., nonlinear) manner, the additive (i.e., linear) idea that ‘total infection’ minus ‘cellfree infection’ is equal to ‘celltocell infection’ does not hold true universally. Hence, it was difficult to estimate and compare the efficacies of cellfree and celltocell infection, and different reports provided different estimates (Dimitrov et al., 1993; Carr et al., 1999; Chen et al., 2007; Sourisseau et al., 2007; Zhong et al., 2013). Thus, the quantitative information that underlies celltocell infection has yet to be elucidated and its impact on virus spread remains unclear.
In this study, through coupled experimental and mathematical investigation, we demonstrate that the efficacy of celltocell HIV1 infection is 1.4fold higher than that of cellfree infection (i.e., celltocell infection accounts for approximately 60% of total infection). We also show that the celltocell infection shortens the generation time of viruses by 0.9 times, and increases the viral fitness by 3.9 times. These findings strongly suggest that the celltocell infection plays a critical role in efficient and rapid spread of viral infection. Furthermore, we discuss the role of the celltocell infection in HIV1 infected individuals, based on in silico simulation with our estimated parameters.
Results
Adaptation of a mathematical model to explicitly consider cellfree and celltocell infection
A static cell culture system (i.e., a conventional cell culture system) allows viruses to perform both cellfree and celltocell infection. On the other hand, Sourisseau et al. have reported that the celltocell infection can be prevented by mildly shaking the cell culture infected with viruses (Sourisseau et al., 2007). Consistent with the previous report (Sourisseau et al., 2007), we verified that shaking did not induce nonspecific consequences on HIV1 infection (Figure 2—figure supplement 1). To quantitatively estimate the efficacy of the cellfree infection and that of the celltocell infection respectively, we adopted this experimental method (see ‘Materials and methods’). Static cultures of Jurkat cells, an HIV1susceptible human CD4^{+} Tcell line, allow HIV1 to propagate both by the cellfree and celltocell infection, while under shaking conditions, Jurkat cells allows HIV1 to replicate only by the cellfree infection (Figure 1A).
Previous mathematical models, which have been widely used for data analyses, essentially describe only the cellfree infection (Nowak and May, 2000; Perelson, 2002; Iwami et al., 2012a, 2012b) or implicitly both infection modes (Komarova and Wodarz, 2013; Komarova et al., 2013a, 2013b). Here we used the following revised model including both infection modes explicitly:
where T(t) and I(t) are the numbers of uninfected and infected cells per ml of a culture, respectively, and V(t) is the viral load measured by the amount of HIV1 p24 per ml of culture supernatant. The target cells (we used Jurkat cells) grow at a rate g with the carrying capacity of T_{max} (the maximum number of cells in the cell culture flask). The parameters β, δ, p and c represent the cellfree infection rate, the death rate of infected cells, the virus production rate, and the clearance rate of virions, respectively. Note that c, g, and δ include the removal of virus, and of the uninfected and infected cells, due to the experimental samplings. In our earlier works (Iwami et al., 2012a, 2012b; Fukuhara et al., 2013; Kakizoe et al., 2015), we have shown that the approximating punctual removal as a continuous exponential decay has minimal impact on the model parameters and provides an appropriate fit to the experimental data. In addition, we introduce the parameter ω, describing the infection rate via celltocell contacts (Sourisseau et al., 2007; Sattentau, 2008; Sigal et al., 2011). In the shaking cell culture system, we fixed ω = 0 because the shaking inhibits the formation of celltocell contacts completely (Sourisseau et al., 2007). In previous reports, Komarova et al. used a quasiequilibrium approximation for the number of free virus, and incorporated the dynamics of V(t) into that of I(t) in Komarova and Wodarz (2013), Komarova et al. (2013a), and Komarova et al. (2013b). However, in cell culture system, the clearance of virions usually is not much larger than the death rate of infected cells, like in vivo (see below). This fact does not validate the quasiequilibrium approximation, and it may affect the quantification of the dynamics of the celltocell and cellfree infection. We introduced the above full model, relying on a carefully designed experiment, to accurately extract the quantitative information that underlies HIV1 infection. Furthermore, our experimental datasets include all timeseries of the number of uninfected, infected cell, and virions. Thus, our coupled experimental and mathematical investigations with a sufficient datasets allowed us to estimate all parameters in Equations 1–3, and to compute the basic reproduction number, generation time, and Malthus coefficient (see below).
Data fitting to quantify the cellfree and celltocell contribution to HIV spread
Correctly estimated parameter sets with possible variation are required to reproduce model prediction for pure celltocell infection in silico. However, point estimation of the model parameter set by a conventional ordinary least square method does not capture possible variations of kinetic parameters and model prediction. To assess the variability of kinetic parameters and model prediction, we perform Bayesian estimation for the whole dataset using Markov Chain Monte Carlo (MCMC) sampling (see ‘Materials and methods’ and Supplementary file 1), and simultaneously fit Equations 1–3 with ω > 0 and ω = 0 to the concentration of p24negative and positive Jurkat cells and the amount of p24 viral protein in the static and shaking cell cultures, respectively. Here we note that g and T_{max} were separately estimated and fixed to be 0.47 ± 0.10 for the static culture and 0.54 ± 0.09 for the shaking culture per day, and (1.51 ± 0.02) × 10^{6} and (1.22 ± 0.02) × 10^{6} cells per flask of medium from the cell growth experiments, respectively (see ‘Materials and methods’, Figure 2—figure supplement 2 and Supplementary file 2). In addition, we used c value of 2.3 per day, which is estimated from daily harvesting of viruses (i.e., the amount of p24 have to be reduced by around 90% per day by the daily mediumreplacement procedure).
The remaining four common parameters β, ω, δ and p, along with the six initial values for T(0), I(0) and V(0) in the static and the shaking cell cultures, were determined by fitting the model to the data. Experimental measurements, which were below the detection limit, were excluded in the fitting. The estimated parameters of the model and derived quantities are given in Table 1, and the estimated initial values are summarized in Supplementary file 3. The typical behavior of the model using these bestfit parameter estimates is shown together with the data in Figure 2, which reveals that Equations 1–3 describe these in vitro data very well. The shadowed regions correspond to 95% posterior predictive intervals, the dashed lines give the bestfit solution (mean) for Equations 1–3, and the dots show the experimental datasets. This suggests that the parameters that were estimated are representative for the various processes underlying the HIV1 kinetics including the celltocell and cellfree infection.
Our model (i.e., Equations 1–3) applied to timecourse experimental data under static and shaking conditions (i.e., Figure 2A and Figure 2B, respectively) allowed to extract the kinetic parameters in the model (see Table 1), in particular the rate constant for the cellfree infection (β) and the rate constant for the celltocell infection (ω). However, from the estimated values of β and ω, we could not directly compare the efficiency of the two infection modes, because of the different units of measure of these parameters (p24/day for β, and cells/day for ω). To quantify each infection mode and overcome the above difficulty, we derived the basic reproduction number R_{0} (Perelson and Nelson, 1999; Nowak and May, 2000; Iwami et al., 2012b), an index reflecting the average number of newly infected cells produced from any one infected cell (see mathematical appendix in ‘Materials and methods’). Note that secondly infected cells are produced from both the cellfree and celltocell infection. Interestingly, in spite of nonlinear interaction between the two modes of virus transmission, our derivation of R_{0} revealed that the secondly infected cells were the sum of the basic reproduction number through the cellfree infection R_{cf} = βpT_{max}/δc and the basic reproduction number through the celltocell infection R_{cc} = ωT_{max}/δ, (i.e., R_{0} = R_{cf} + R_{cc}) (see Figure 1B). Using all accepted MCMC parameter estimates from the timecourse experimental datasets, we calculated that on average the mean of the total basic reproductive number is R_{0} = 5.83 ± 0.94 (average ± standard deviation), and the mean number of secondly infected cells through the cellfree infection and the celltocell infection are R_{cf} = 2.44 ± 0.23 and R_{cc} = 3.39 ± 0.91, respectively (see Table 1). The distributions of calculated R_{0}, R_{cf}, and R_{cc}, are shown in Figure 3A–C, respectively. These estimates indicate that the contribution of the celltocell infection is almost 60% on average (i.e., R_{cc}/(R_{cc} + R_{cf}) = 0.57 ± 0.07: Table 1) and this mode of infection is predominant during the HIV1 spread in Jurkat cells. In Figure 3D, the distributions of calculated ratio are shown. Interestingly, this estimation is consistent with that by Komarova and Wodarz (2013), Komarova et al. (2013a), and Komarova et al. (2013b), although they did not take into account the difference of the death rate in the shaking and static conditions.
Advantage of celltocell infection
We also derived the viral generation time, defined as the time it takes for a population of virions to infect cells and reproduce (Perelson and Nelson, 1999), from Equations 1–3 in the static and shaking cell cultures (see mathematical appendix in ‘Materials and methods’). In the presence and absence of the celltocell infection (i.e., for the static and shaking cell cultures, respectively), the mean generation time is calculated as 1/δ + R_{cf}/cR_{0} = 2.22 ± 0.32 days and 1/δ + 1/c = 2.47 ± 0.32 days, respectively (see Table 2). Thus, celltocell infection shortens the generation time by on average 0.90 times, and enables HIV1 to efficiently infect target cells (Sato et al., 1992; Carr et al., 1999). Furthermore, we calculated the Malthus coefficient, defined as the fitness of virus (Nowak and May, 2000; Nowak, 2006) (or the speed of virus infection) (see mathematical appendix in ‘Materials and methods’). In the presence and absence of the celltocell infection, the Malthus coefficient is calculated as 1.86 ± 0.37 and 0.49 ± 0.05 per day, respectively (see Table 2). Thus, celltocell infection increases the HIV1 fitness by 3.80fold (corresponding to 944fold higher viral load 5 days after the infection) and plays an important role in the rapid spread of HIV1. Thus, the efficient viral spread via the celltocell infection is relevant, especially at the beginning of virus infection.
Virtual experiments of celltocell infection in silico
While the shaking culture prevents the celltocell infection, it is technically difficult to completely block the cellfree infection. Here, using our estimated kinetic parameters (Table 1 and Supplementary file 3), we carried out a ‘virtual experiment’ eliminating the contribution of the cellfree infection using all accepted MCMC estimated parameter values, allowing to estimate only the celltocell infection, in silico (see Figure 4). Our simulated mean values (represented by solid lines) of the celltocell infection of HIV1 are consistently located between the time course of experimental data under the static conditions (closed circles, including both the cellfree and cell to cell infections) and those under the shaking conditions (open circles, reflecting only the cellfree infection). The shadowed regions correspond to 95% posterior predictive intervals. In terms of the dynamics of infected cells and virus production, the simulated values corresponding to celltocell virus propagation, are closer to experimental data from the coupled cellfree and celltocell infection, than to data from the cellfree infection only. This shows that the cellfree infection, which contributes approximately 40% to the whole HIV1 infection process, plays a limited role on the virus spread. In other words, even if we could completely block the cellfree infection, the celltocell infection would still effectively spread viruses (Sigal et al., 2011). We address this point in ‘Discussion’.
Discussion
Through experimentalmathematical investigation, here we quantitatively elucidated the dynamics of the celltocell and cellfree HIV1 infection modes (Figure 2 and Table 1). We derived the basic reproduction number, R_{0}, and divided it into the numbers of secondly infected cells through the cellfree infection, R_{cf}, and the celltocell infection, R_{cc}, respectively (Figure 1B and mathematical appendix in ‘Materials and methods’). Based on our calculated values of these three indexes, we found that about 60% of the viral infection is attributed to the celltocell infection in the in vitro cell culture system (Table 1), which is consistent with previous estimation by Komarova and Wodarz (2013), Komarova et al. (2013a), and Komarova et al. (2013b). In addition, we revealed that the celltocell infection effectively promotes the virus infection by reducing the generation time (×0.9 times), and by increasing the Malthus coefficient (×3.80 times) (Table 2).
When we consider the significance of the celltocell infection in patients infected with HIV1, it should be noted that the environment of immune cells including CD4^{+} Tcells in vivo is radically different from the conditions of in vitro cell cultures. For instance, lymphocytes are closely packed in lymphoid tissues such as lymph nodes, and thereby, the frequency for the infected cell to contact with adjacent uninfected cells in vivo would be much higher than that in in vitro cell cultures. In addition, Murooka et al. have directly demonstrated that HIV1infected cells converge to lymph nodes and can be vehicles for viral dissemination in vivo (Murooka et al., 2012). Moreover, certain studies have suggested that celltocell viral spread is resistant to antiviral immunity such as neutralizing antibodies and cytotoxic T lymphocytes (Martin and Sattentau, 2009). Therefore, these notions strongly suggest that the contribution of the celltocell infection for viral propagation in vivo may be much higher than that estimated from the in vitro cell culture system.
As another significance of celltocell viral spread, Sigal et al. have suggested that the celltocell infection permits viral replication even under the antiretroviral therapy (Sigal et al., 2011). This is attributed to the fact that the multiplicity of infection per cell is tremendously higher than that reached by an infectious viral particle. However, in the previous report (Sigal et al., 2011), the contribution of the celltocell infection remained unclear. To further understand the role of the celltocell infection, we quantified the contributions of the celltocell and cellfree infection modes (Table 1). Interestingly, we found that the celltocell infection mode is predominant during the infection. Furthermore, our virtual experiments showed that a complete block of the cellfree infection, which is highly susceptible to current antiviral drugs, provides only a limited impact on the whole HIV1 infection (Figure 3). Taken together, our findings further support that the celltocell infection can be a barrier to prevent the cure of HIV1 infection, which is discussed in Sigal et al. (2011). However, it should be noted that some papers have shown that celltocell spread cannot overcome the action of most antiHIV1 drugs (Titanji et al., 2013; Agosto et al., 2014). To fully elucidate this issue, further investigations will be needed.
In addition to HIV1, other viruses such as herpes simplex virus, measles virus, and human hepatitis C virus drive their dissemination via celltocell infection (Sattentau, 2008; TalbertSlagle et al., 2014). Although the impact of celltocell viral spread is a topic of broad interest in virology, it was difficult to explore this issue by conventional virological experiments, because an infected cell is simultaneously capable of achieving celltocell infection along with producing infectious viral particles. By applying mathematical modeling to the experimental data, here we estimated the sole dynamics of cellfree infection in the cell culture system. The synergistic strategy of experiments with mathematical modeling is a powerful approach to quantitatively elucidate the dynamics of virus infection in a way that is inaccessible through conventional experimental approaches.
Materials and methods
Cell culture and HIV1 infection
Request a detailed protocolJurkat cell line (Watanabe et al., 2012) was cultured in the culture medium: RPMI 1640 (Sigma, St. Louis, MO) containing 2% fetal calf serum and antibiotics. The virus solution was prepared as previously described (Sato et al., 2010, 2013, 2014; Iwami et al., 2012a). Briefly, 30 μg of pNL43 plasmid (Adachi et al., 1986) (GenBank accession no. M19921.2) was transfected into 293T cells by the calciumphosphate method. At 48 hr posttransfection, the culture supernatant was harvested, centrifuged, and then filtered through a 0.45μmporesize filter to produce virus solution. The infectivity of virus solution was titrated as previously described (Iwami et al., 2012a). Briefly, the virus solution obtained was serially diluted and then inoculated onto phytohemagglutininstimulated human peripheral blood mononuclear cells in a 96well plate in triplicate. At 14 days postinfection, the endpoint was determined by using an HIV1 p24 antigen enzymelinked immunosorbent assay (ELISA) kit (ZetptoMetrix, Buffalo, NY) according to the manufacture's procedure, and virus infectivity was calculated as the 50% tissue culture infectious doses (TCID_{50}) according to the ReedMuench method.
For HIV1 infection, 3 × 10^{5} of Jurkat cells were infected with HIV1 (multiplicity of infection 0.1) at 37°C for 2 hr. The infected cells were washed three times with the culture medium, and then suspended with 3 ml of culture medium and seeded into a 25cm^{2} flask (Nunc, Rochester, NY). For the static infection, the infected cell culture was kept in a 37°C/5% CO_{2} incubator as usual. For the shaking infection, the infected cell culture was handled as previously described (Sourisseau et al., 2007). Briefly, the cell culture was kept on a Petit rocker Model2230 (Wakenyaku, Japan) placed in 37°C/5% CO_{2} incubator, and was gently shaken at 40 movements per min. The amount of virus particles in the culture supernatant and the number of infected cells were routinely measured as follows: a portion (300 μl) of the infected cell culture was routinely harvested, and the amount of released virions in the culture supernatant was quantified by using an HIV1 p24 antigen ELISA kit (ZetptoMetrix) according to the manufacture's procedure. The cell number was counted by using a Scepter handled automated cell counter (Millipore, Germany) according to the manufacture's protocol. The percentage of infected cells was measured by flow cytometry. Flow cytometry was performed with a FACSCalibur (BD Biosciences, San Jose, CA) as previously described (Sato et al., 2010; Sato et al., 2011, 2013, 2014; Iwami et al., 2012a), and the obtained data were analyzed with CellQuest software (BD Biosciences). For flow cytometry analysis, a fluorescein isothiocyanatelabeled antiHIV1 p24 antibody (KC57; Beckman Coulter, Pasadena, CA) was used. The representative dot plots are shown in Figure 2—figure supplement 3. The data is available upon request. The remaining cell culture was centrifuged and then resuspended with 3 ml of fresh culture medium. It should be noted that the procedure for HIV1 infection was performed at time t = −2 day in the figures. Because there is no viral protein production in the first day after infection, each in vitro experimental quantity was measured daily from t = 0 day (i.e., 2 days after HIV1 inoculation). The detection threshold of each value are the followings: cell number (cell counting), 3000 cells/ml; % p24positive cells (flow cytometry), 0.3%; and p24 antigen in culture supernatant (p24 antigen ELISA), 80 pg/ml.
Parameter estimation
Request a detailed protocolA statistical model adopted in the Bayesian inference assumes measurement error to follow normal distribution with mean zero and unknown variance (error variance). A distribution of error variance is also inferred with the Gamma distribution as its prior distribution. Posterior predictive parameter distribution as an output of MCMC computation represents parameter variability. Distributions of model parameters and initial values were inferred directly by MCMC computations. On the other hand, distributions of the basic reproduction numbers and the other quantities were calculated from the inferred parameter sets (Figure 3 for graphical representation). A set of computations for Equations 1–3 with estimated parameter sets gives a distribution of outputs (virus load and cell density) as model prediction. To investigate variation of model prediction, global sensitivity analyses were performed. The range of possible variation is drawn in Figure 2 as 95% confidence interval. Technical details of MCMC computations are summarized in Supplementary file 1.
Quantification of Jurkat cell growth
Request a detailed protocolWe here estimate the growth kinetics of Jurkat cells, which have been commonly used for HIV1 studies, under the normal (i.e., mockinfected) condition with the following mathematical model:
where the variable T(t) is the number of Jurkat cells at time t and the parameters g and T_{max} are the growth rate of the cells (i.e., Log2/g is the doubling time) and the carrying capacity of the cell culture flask, respectively. Nonlinear leastsquares regression (FindMinimum package of Mathematica9.0) was performed to fit Equation 4 to the timecourse numbers of Jurkat cells in the normal condition. The fitted parameter values are listed in Supplementary file 2 and the model behavior using these bestfit parameter estimates is presented together with the data in Figure 2—figure supplement 2.
Mathematical appendix
Request a detailed protocolThe linearized equation of Equations 1–3 at the virusfree steady state, (T_{max}, 0, 0), is given as follows:
Let b(t) be the number of newly produced infected cells in the linear phase:
Applying the variation of constants formula to Equations 5, 6, we have
Inserting Equation 9 into Equation 8 to exchange the order of integrals, we have
where
From Equation 7 and Equations 9, 10, we arrive at the following renewal equation:
where h(t) is given by
and the kernel Ψ(x) is given by
In the above expression, ϕ_{j}(x) denotes the probability density function given by
and, ∗ denotes the convolution of functions. From the general theory of the basic reproduction number (Inaba, 2012), R_{0} for the reproduction of infected cells is given by
where R_{cf} and R_{cc} denote the reproduction numbers for infected cells mediated by the cellfree and celltocell infection, respectively.
Next we consider the reproduction process of viruses. Let ρ(t):= pI(t) be the number of newly produced viruses at time t. From Equations 8, 9, we obtain
where
Inserting Equation 11 into Equation 12, we again arrive at the following renewal equation:
where q(t) is given by
Note that the reproduction kernel Ψ(x) for the virus reproduction is the same as the kernel for the cell reproduction. Thus the probability density function of the virus reproduction is given by
Then the generation time for the virus reproduction, denoted by G, is calculated as follows:
where G_{cf} : = 1/δ + 1/c and G_{cc} : = 1/δ are the generation times for virus reproduction mediated by the cellfree and celltocell infection, respectively.
The Malthusian coefficient for the virus reproduction must be given as the dominant real root of the EulerLotka equation as
where $\widehat{{\varphi}_{j}}$ denotes the Laplace transformation of a function ϕ_{j}. That is,
Therefore the EulerLotka equation can be calculated explicitly as follows:
which is reduced to a quadratic equation,
If R_{0} > 1, Equation 13 has a unique positive root, which is no other than the Malthusian coefficient for the virus reproduction, so it is calculated as,
References

Production of acquired immunodeficiency syndromeassociated retrovirus in human and nonhuman cells transfected with an infectious molecular cloneJournal of Virology 59:284–291.

Quantitation of human immunodeficiency virus type 1 infection kineticsJournal of Virology 67:2182–2190.

On a new perspective of the basic reproduction number in heterogeneous environmentsJournal of Mathematical Biology 65:309–348.https://doi.org/10.1007/s002850110463z

Identifying viral parameters from in vitro cell culturesFrontiers in Microbiology 3:319.https://doi.org/10.3389/fmicb.2012.00319

A conservation law for virus infection kinetics in vitroJournal of Theoretical Biology 376:39–47.https://doi.org/10.1016/j.jtbi.2015.03.034

Virus dynamics in the presence of synaptic transmissionMathematical Biosciences 242:161–171.https://doi.org/10.1016/j.mbs.2013.01.003

Celltocell HIV1 spread and its implications for immune evasionCurrent Opinion in HIV and AIDS 4:143–149.https://doi.org/10.1097/COH.0b013e328322f94a

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

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

Avoiding the void: celltocell spread of human virusesNature Reviews Microbiology 6:815–826.https://doi.org/10.1038/nrmicro1972

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

The hematopoietic cellspecific Rho GTPase inhibitor ARHGDIB/D4GDI limits HIV type 1 replicationAIDS Research and Human Retroviruses 28:913–922.https://doi.org/10.1089/AID.2011.0180
Decision letter

Arup K ChakrabortyReviewing Editor; Massachusetts Institute of Technology, United States
eLife posts the editorial decision letter and author response on a selection of the published articles (subject to the approval of the authors). An edited version of the letter sent to the authors after peer review is shown, indicating the substantive concerns or comments; minor concerns are not usually shown. Reviewers have the opportunity to discuss the decision before the letter is sent (see review process). Similarly, the author response typically shows only responses to the major concerns raised by the reviewers.
Thank you for submitting your work entitled “Celltocell infection by HIV contributes over half of virus infection” for peer review at eLife. Your submission has been favorably evaluated by Naama Barkai (Senior Editor), a Reviewing Editor, and two reviewers. One of the two reviewers, Vitaly Ganusov, has agreed to share his identity.
The reviewers have discussed the reviews with one another and the Reviewing Editor has drafted this decision to help you prepare a revised submission.
While others have previously shown that celltocell spread is a rapid and efficient form of infection, the novelty of this work is in its attempt to quantify the contributions of celltocell spread and cellfree infection in the same culture. However, there are three major concerns about the connections between the experiments and assumptions made in the modeling that need to be addressed in order to ascertain that the conclusions are correct. These major concerns are:
1) The parameter, ω, which measures the infection rate by celltocell spread is obtained by fitting the model to data, and it is assumed that ω=0 in the experiments with shaking. β, the cellfree infection rate, is assumed to be the same under both infection conditions. In the data, shaking reduces infection, consistent with decreased celltocell spread. However, it is also consistent with several other possibilities. For example, reduction of both celltocell spread and cellfree infection due to many possible factors including lower production rates or infectivity of virus. If this is true, a similar attenuation of viral dynamics would be observed upon shaking but for a completely different reason, and invalidate the assumptions made to fit the data. To prove that the attenuated viral dynamics are strictly the result of a loss of celltocell spread, a control experiment is needed. It is up to you to choose the appropriate control experiment that establishes this. We think that the following experiment might be appropriate:
Produce cellfree virus by filtering out any cells, then add it to either shaking or static Jurkat cell cultures and measure the percent infected cells under each condition at a few time points before two days (i.e. before the second cycle of infection). If shaking only affects celltocell transmission, the results should be identical for shaken or static cultures in this experiment.
2) The parameter values are taken to be constant in time, but this may not be the case. Consider the parameter g, as an example. Cells that reach a maximum density may go into growth arrest and there is a lag time to get them out of it. This would mean that g is not constant. Is there experimental evidence that the parameters do not change with time. Also, at high levels of infection, which appears to be true in the experiments, is the cell culture still stable? If not, g would not be a constant.
3) It seems that the simulations are carried out without virus removal (without c) for 24 hours, and then an instantaneous removal occurs at 90% (not continuous removal of the virus as 2.3 per day). This setup does not seem congruent with the experiments. The simulations should be carried out in a similar way to the experimental setup to make sure that the parameter values are properly estimated. In addition, is it certain that cells/infected cells are not removed in the experiments, as it seems is assumed in the model? Finally, more quantitative estimates of agreement between the fitted models and experimental data are required. For example, by visual inspection of Figure 2, one could argue that the model does not describe the loss of infected cells and the virus after the peak well (Figure 2A), and in Figure 2B, it is not clear how good the fits are. Error bars need to be provided for the estimated parameters in Table 1.
[Editors' note: further revisions were requested prior to acceptance, as described below.]
Thank you for resubmitting your work entitled “Celltocell infection by HIV contributes over half of virus infection” for further consideration at eLife. Your revised article has been favorably evaluated by Naama Barkai (Senior Editor), a Reviewing Editor, and one of the original reviewers. The manuscript has been improved but there are some remaining issues that need to be addressed before acceptance, as outlined below:
We have read your response and revised manuscript, and believe that you have not addressed the first major and essential point. As discussed in the first report, the assumption you make is that shaking eliminates celltocell spread but does not affect cellfree infection. The premise of the whole paper is based on this assumption. According to the data presented in Figure 2, this may not be the case: Day 0 infection as shown in the middle panels of Figure 2A (static) and Figure 2B (shaking) is strictly the result of cellfree infection, as you indicate in the Methods. In Figure 2A (static), the number of cell infected is between 10^{3} to 10^{4}/ml. In Figure 2B (shaking), the number of infected cells – by what is supposed to be the same cellfree infection input according to the methods – seems to be undetectable at day 0, with the first measurements appearing at a later time. This indicates that shaking, at least as performed by the method used in this paper, decreases cellfree infection. We requested that you perform a control experiment to show that this is not the case, and that shaking does not affect cellfree infection. We had also suggested a simple experiment for you to consider: add cellfree virus to static and shaking cultures and show that the infection rate is unaffected by shaking at the first infection cycle (subsequent cycles would involve celltocell spread due to the coculture, and would not be relevant). Rather than carry out this (or some other) control experiment, you point to a past paper that showed this. But, given the data in Figure 2 (see above), it is unclear that the assumption that shaking does not affect cell free infection is correct in your experiment. Since the rest of the work presented in the paper relies on this assumption, a control experiment demonstrating this assumption to be true is necessary.
https://doi.org/10.7554/eLife.08150.015Author response
While others have previously shown that celltocell spread is a rapid and efficient form of infection, the novelty of this work is in its attempt to quantify the contributions of celltocell spread and cellfree infection in the same culture. However, there are three major concerns about the connections between the experiments and assumptions made in the modeling that need to be addressed in order to ascertain that the conclusions are correct. These major concerns are:
1) The parameter, ω, which measures the infection rate by celltocell spread is obtained by fitting the model to data, and it is assumed that ω=0 in the experiments with shaking. β, the cellfree infection rate, is assumed to be the same under both infection conditions. In the data, shaking reduces infection, consistent with decreased celltocell spread. However, it is also consistent with several other possibilities. For example, reduction of both celltocell spread and cellfree infection due to many possible factors including lower production rates or infectivity of virus. If this is true, a similar attenuation of viral dynamics would be observed upon shaking but for a completely different reason, and invalidate the assumptions made to fit the data. To prove that the attenuated viral dynamics are strictly the result of a loss of celltocell spread, a control experiment is needed. It is up to you to choose the appropriate control experiment that establishes this. We think that the following experiment might be appropriate:
Produce cellfree virus by filtering out any cells, then add it to either shaking or static Jurkat cell cultures and measure the percent infected cells under each condition at a few time points before two days (i.e. before the second cycle of infection). If shaking only affects celltocell transmission, the results should be identical for shaken or static cultures in this experiment.
In the original publication (Sourisseau et al., J Virol, 2007) describing the experimental model that we used here, the authors carefully verified that shaking did not induce nonspecific consequences on HIV infection. Namely, they verified that static and shaking conditions did not induce differences in cell viability, cell growth rate, expression level of receptor/coreceptors and adhesion molecules (Figure 2 in Sourisseau et al., J Virol, 2007). The authors also verified the absence of effect on virus release and on the efficiency of cellfree virus infection (Figure 3 in Sourisseau et al., J Virol, 2007). The previous demonstration of the absence of nonspecific consequences on HIV replication is now indicated in the manuscript (subsection “Adaptation of a mathematical model to explicitly consider cellfree and celltocell infection”).
2) The parameter values are taken to be constant in time, but this may not be the case. Consider the parameter g, as an example. Cells that reach a maximum density may go into growth arrest and there is a lag time to get them out of it. This would mean that g is not constant. Is there experimental evidence that the parameters do not change with time. Also, at high levels of infection, which appears to be true in the experiments, is the cell culture still stable? If not, g would not be a constant.
Thank you for pointing this out, we would like to explain our assumption in detail, especially for the growth rate. In Equation (1), we assumed the growth of target cells is described by g{1−(T(t)+I(t))/T_{max}}T(t). This is called “logistic growth”. In this formulation, g{1−(T(t)+I(t))/T_{max}} represents the average growth rate per cell. That is, if the density of total cells (i.e. T(t)+I(t)) is low, then the growth rate is approximately g, but if the density is high (i.e., around maximum density, T_{max}), then the growth rate becomes 0, which models cell growth arrest. Because the term of 1−(T(t)+I(t))/T_{max} is a function of time, the average growth rate gradually decrease from 1 to 0 as the density increases from low to maximum, which shows a lag time. Thus, the logistic formulation well captures the densitydependent cell growth in cell culture as we previously showed in Fukuhara et al. (2013). In contrast, for other parameters, we simply assumed they are constant. However, at least, in cell culture system, many papers including our own studies (Fukuhara et al., 2013; Iwami et al., Retrovirology, 2012; Iwami et al., Front Microbiol., 2012; Komarova et al., 2013; Kakizoe et al., 2015; Beauchemin et al., 2008) has shown that the simple constant assumption sufficiently reproduced the infection dynamics and quantified parameters reasonably. Therefore, we would like to keep our parameters constant in order to avoid the complexity of timedependent parameters and to allow the use of our novel model for the celltocell and cellfree HIV1 infections.
3) It seems that the simulations are carried out without virus removal (without c) for 24 hours, and then an instantaneous removal occurs at 90% (not continuous removal of the virus as 2.3 per day). This setup does not seem congruent with the experiments. The simulations should be carried out in a similar way to the experimental setup to make sure that the parameter values are properly estimated. In addition, is it certain that cells/infected cells are not removed in the experiments, as it seems is assumed in the model? Finally, more quantitative estimates of agreement between the fitted models and experimental data are required. For example, by visual inspection of Figure 2, one could argue that the model does not describe the loss of infected cells and the virus after the peak well (Figure 2A), and in Figure 2B, it is not clear how good the fits are. Error bars need to be provided for the estimated parameters in Table 1.
In addressing this very relevant and important comment, we have made a number of changes to our manuscript but also to our analysis itself. First, we would like to explain our assumption in Equations (13) for the removal due to the experimental samplings. For each of the daily measurements of the virus concentration, the medium in our experiments was harvested, reducing the viral concentration by 90%. This removal can be captured using Equation (3) by approximating the punctual removal of virus at each sampling time as a continuous, exponential decay of the viral load over the period between samples. In such a case, parameter c corresponds to the sum of the rate of virus loss due to harvesting of the medium plus the rate of loss due to degradation of the extracellular virus (which is negligibly small). In addition, on a daily basis, 10% of the cells in the culture were harvested to measure the number of target cells and infected cells. Similarly, the removals of the target and infected cells were included in the value of g and δ, respectively. We clarified this assumption by adding a few sentences to the subsection “Adaptation of a mathematical model to explicitly consider cellfree and celltocell infection”.
As an example, we have expanded our analysis using Equations (13) to also include an analysis of the experimental data using the model with punctual removal, i.e. we determined bestfit parameters for this punctual model and present graphs of its agreement to both data sets over time (Author response image 1). Generally, the punctual removal model describes the experimental data well, and the estimated values for the model parameters are similar to those found using Equations (13); β=3.49×10^{−6}, ω=1.17×10^{−6}, p=0.53, δ=0.50 in Eqs.(13), and β=6.65×10^{−6}, ω=1.69×10^{−6}, p=0.13, δ=0.43 in the punctual model. In our earlier works (Fukuhara et al., 2013; Iwami et al., Retrovirology, 2012; Iwami et al., Front Microbiol., 2012; Kakizoe et al., 2015) we have also shown that approximating punctual removal as a continuous exponential decay has minimal impact on the model parameters and provides an appropriate fit to the experimental data. And also, unfortunately, we could not define the basic reproduction numbers such as R_{cf}, R_{cc}, and R_{0} for the punctual model because of its discontinuity. Quantifying these values and comparing them are key findings of this paper. Therefore, we used the exponential decay approximation (i.e. Equations (13) here. Nevertheless, we do feel that inclusion of the above discussion about the punctual removal is a valuable addition to our manuscript and improves the completeness of our work (see the second paragraph of the subsection “Adaptation of a mathematical model to explicitly consider cellfree and celltocell infection”).
Furthermore, for more quantitative estimates of agreement between the fitted models and experimental data, we performed Bayesian estimation for the whole dataset using Markov Chain Monte Carlo (MCMC) sampling. The Bayesian estimation enables us to assess the variability of kinetic parameters and model prediction as posterior predictive intervals (see Results, Methods and Supplementary file 1). We thank the reviewers for leading us in this direction.
[Editors' note: further revisions were requested prior to acceptance, as described below.]
We have read your response and revised manuscript, and believe that you have not addressed the first major and essential point. As discussed in the first report, the assumption you make is that shaking eliminates celltocell spread but does not affect cellfree infection. The premise of the whole paper is based on this assumption. According to the data presented in Figure 2, this may not be the case: Day 0 infection as shown in the middle panels of Figure 2A (static) and Figure 2B (shaking) is strictly the result of cellfree infection, as you indicate in the Methods. In Figure 2A (static), the number of cell infected is between 10^{3} to 10^{4}/ml. In Figure 2B (shaking), the number of infected cells – by what is supposed to be the same cellfree infection input according to the methods – seems to be undetectable at day 0, with the first measurements appearing at a later time. This indicates that shaking, at least as performed by the method used in this paper, decreases cellfree infection. We requested that you perform a control experiment to show that this is not the case, and that shaking does not affect cellfree infection. We had also suggested a simple experiment for you to consider: add cellfree virus to static and shaking cultures and show that the infection rate is unaffected by shaking at the first infection cycle (subsequent cycles would involve celltocell spread due to the coculture, and would not be relevant). Rather than carry out this (or some other) control experiment, you point to a past paper that showed this. But, given the data in Figure 2 (see above), it is unclear that the assumption that shaking does not affect cell free infection is correct in your experiment. Since the rest of the work presented in the paper relies on this assumption, a control experiment demonstrating this assumption to be true is necessary.
According to the suggestion raised by the editors (“it is unclear that the assumption that shaking does not affect cell free infection is correct in your experiment”), we carried out an additional experiment and the data is shown as Figure 2–figure supplement 1. As shown, we have clearly verified that the shaking procedure does not affect the efficacy of cellfree infection.
https://doi.org/10.7554/eLife.08150.016Article and author information
Author details
Funding
Japan Science and Technology Agency
 Shingo Iwami
The funder had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
This work was supported in part by JST PRESTO program (to SI); JST CREST program (to SI, HI, KA, and KS); GrantsinAid for Young Scientists B25800092 (to SI) and B25871132 (to SN) from the Japan Society for the Promotion of Science (JSPS); JSPS KAKENHI Grant Number 10192783 and 15KT0107 (to SI), 25400194 (to HI) and 15K07166 (to KS); Inamori Foundation (to SI); the Aihara Innovative Mathematical Modeling Project, JSPS, through the ‘Funding Program for WorldLeading Innovative R & D on Science and Technology (FIRST Program)’, initiated by Council for Science and Technology Policy (to SI, SN, HI, KA, and KS); the Japan Agency for Medical Research and Development, AMED (H27ShinkoJitsuyokaGeneral016) (to SI, SN, HI, KA, and KS); Agence Nationale de Recherches sur le Sida et les Hepatites Virales (ANRS) (to FM); a GrantinAid for Scientific Research on Innovative Areas 24115008 from the Ministry of Education, Culture, Sports, Science and Technology of Japan (to YK); JSPS CoretoCore program, A. Advanced Research Networks (to YK); Research on HIV/AIDS from AMED 15Afk0410013h0001 (to YK); Takeda Science Foundation (to KS); Sumitomo Foundation Research Grant (to KS); Senshin Medical Research Foundation (to KS); Imai Memorial Trust for AIDS Research (to KS); Ichiro Kanehara Foundation (to KS); Kanae Foundation for the Promotion of Medical Science (to KS); Suzuken Memorial Foundation (to KS); Uehara Memorial Foundation (to KS).
Reviewing Editor
 Arup K Chakraborty, Massachusetts Institute of Technology, United States
Publication history
 Received: April 16, 2015
 Accepted: September 4, 2015
 Version of Record published: October 6, 2015 (version 1)
Copyright
© 2015, Iwami 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

 4,627
 Page views

 763
 Downloads

 76
 Citations
Article citation count generated by polling the highest count across the following sources: Crossref, Scopus, PubMed Central.