Cell-to-cell infection by HIV contributes over half of virus infection

Abstract
eLife digest
Introduction
Results
Discussion
Materials and methods
References
Article and author information
Metrics

Abstract

Cell-to-cell 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 cell-free infection and cell-to-cell infection, the quantitative information that underlies cell-to-cell 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 cell-to-cell and cell-free human immunodeficiency virus type 1 (HIV-1) infections through experimental-mathematical investigation. Our analyses demonstrated that the cell-to-cell 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 cell-free infection would provide only a limited impact on HIV-1 spread.

https://doi.org/10.7554/eLife.08150.001

eLife digest

Viruses such as HIV-1 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 ‘cell-free’ and ‘cell-to-cell’ infection. In cell-free infection, the virus is released into the fluid that surrounds cells and moves from there into the next cell. In cell-to-cell infection the virus instead moves directly between cells across regions where the two cells make contact.

Previous research has suggested that cell-to-cell infection is important for the spread of HIV-1. However, it is not known how much the virus relies on this process, as it is technically challenging to perform experiments that prevent cell-free infection without also stopping cell-to-cell infection.

Iwami, Takeuchi et al. have overcome this problem by combining experiments on laboratory-grown cells with a mathematical model that describes how the different infection methods affect the spread of HIV-1. This revealed that the viruses spread using cell-to-cell infection about 60% of the time, which agrees with results previously found by another group of researchers. Iwami, Takeuchi et al. also found that cell-to-cell 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 laboratory-grown 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 cell-free method of infection would not greatly affect how well HIV-1 is able to infect new cells. It may instead be more effective to develop HIV treatments that prevent cell-to-cell infection by the virus.

https://doi.org/10.7554/eLife.08150.002

Introduction

In in vitro cell cultures and in infected individuals, viruses may display two types of replication strategies: cell-free infection and cell-to-cell infection (Sattentau, 2008; Martin and Sattentau, 2009; Talbert-Slagle 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 cell-free transmission, or concentrated in the confined space of cell-to-cell contacts between an infected cell and bystander target cells in the case of cell-to-cell transmission. It has been shown that most enveloped virus species, including human immunodeficiency virus type 1 (HIV-1), a causative agent of AIDS, spread via cell-to-cell infection, and it is considered that the replication efficacy of cell-to-cell infection is much higher than that of cell-free infection (Sattentau, 2008; Martin and Sattentau, 2009; Talbert-Slagle et al., 2014). However, it is technically impossible to let viruses execute only cell-to-cell 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 ‘cell-free infection’ is equal to ‘cell-to-cell infection’ does not hold true universally. Hence, it was difficult to estimate and compare the efficacies of cell-free and cell-to-cell 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 cell-to-cell 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 cell-to-cell HIV-1 infection is 1.4-fold higher than that of cell-free infection (i.e., cell-to-cell infection accounts for approximately 60% of total infection). We also show that the cell-to-cell 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 cell-to-cell infection plays a critical role in efficient and rapid spread of viral infection. Furthermore, we discuss the role of the cell-to-cell infection in HIV-1 infected individuals, based on in silico simulation with our estimated parameters.

Results

Adaptation of a mathematical model to explicitly consider cell-free and cell-to-cell infection

A static cell culture system (i.e., a conventional cell culture system) allows viruses to perform both cell-free and cell-to-cell infection. On the other hand, Sourisseau et al. have reported that the cell-to-cell 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 HIV-1 infection (Figure 2—figure supplement 1). To quantitatively estimate the efficacy of the cell-free infection and that of the cell-to-cell infection respectively, we adopted this experimental method (see ‘Materials and methods’). Static cultures of Jurkat cells, an HIV-1-susceptible human CD4⁺ T-cell line, allow HIV-1 to propagate both by the cell-free and cell-to-cell infection, while under shaking conditions, Jurkat cells allows HIV-1 to replicate only by the cell-free infection (Figure 1A).

Figure 1

Download asset Open asset

Cell culture systems and the basic reproduction number under cell-to-cell and cell-free infection.

(A) Static and shaking cultures of Jurkat cells. The static and shaking cell cultures allow human immunodeficiency virus type 1 (HIV-1) to perform both cell-free and cell-to-cell infection, and only cell-free infection, respectively. (B) The basic reproduction number, R₀, is defined as the number of the secondly infected cells produced from a typical infected cell during its infectious period. In the presence of the cell-to-cell and cell-free infection, the basic reproduction number consists of two sub-reproduction numbers through the cell-free infection, R_cf, and through the cell-to-cell infection, R_cc, respectively.

https://doi.org/10.7554/eLife.08150.003

Previous mathematical models, which have been widely used for data analyses, essentially describe only the cell-free 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:

\frac{d T (t)}{d t} = g T (t) (1 - \frac{T (t) + I (t)}{T_{m a x}}) - β T (t) V (t) - ω T (t) I (t),

\frac{d I (t)}{d t} = β T (t) V (t) + ω T (t) I (t) - δ I (t),

\frac{d V (t)}{d t} = p I (t) - c V (t),

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 HIV-1 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 cell-free 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 cell-to-cell 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 cell-to-cell contacts completely (Sourisseau et al., 2007). In previous reports, Komarova et al. used a quasi-equilibrium 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 quasi-equilibrium approximation, and it may affect the quantification of the dynamics of the cell-to-cell and cell-free infection. We introduced the above full model, relying on a carefully designed experiment, to accurately extract the quantitative information that underlies HIV-1 infection. Furthermore, our experimental datasets include all time-series 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 cell-free and cell-to-cell contribution to HIV spread

Correctly estimated parameter sets with possible variation are required to reproduce model prediction for pure cell-to-cell 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 p24-negative 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⁶ and (1.22 ± 0.02) × 10⁶ 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 medium-replacement 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 best-fit 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 best-fit 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 HIV-1 kinetics including the cell-to-cell and cell-free infection.

Table 1

Parameters estimated by mathematical-experimental analysis

https://doi.org/10.7554/eLife.08150.004

Parameter name	Symbol	Unit	Exp. 1	Exp. 2	Exp. 3	Ave. ± S.D.
Parameters obtained from simultaneous fit to time-course experimental dataset
Rate constant for cell-free infection	β	10⁻⁶ × (p24 day)⁻¹	5.59* (3.54–8.41)†	3.27 (2.05–5.01)	‡3.70 (2.28–5.77)	4.18 ± 1.41
Rate constant for cell-to-cell infection	ω	10⁻⁶ × (cell day)⁻¹	0.88 (0.45–1.39)	1.25 (0.70–1.97)	1.13 (0.64–1.79)	1.09 ± 0.33
Production rate of total viral protein	p	day⁻¹	0.37 (0.22–0.59)	0.59 (0.34–0.92)	0.54 (0.31–0.86)	0.50 ± 0.16
Death rate of infected cells	δ	day⁻¹	0.45 (0.32–0.64)	0.54 (0.38–0.75)	0.50 (0.36–0.68)	0.50 ± 0.10
Quantities derived from fitted values
Basic reproduction number through cell-free infection	R_cf	–	2.88 (2.34–3.53)	2.27 (1.98–2.66)	2.43 (2.04–2.95)	2.44 ± 0.23
Basic reproduction number through cell-to-cell infection	R_cc	–	2.95 (1.48–4.70)	3.65 (1.77–6.05)	3.39 (1.82–5.38)	3.39 ± 0.91
Basic reproduction number	R₀	–	5.83 (4.20–7.75)	5.92 (3.99–8.46)	5.83 (4.21–7.89)	5.83 ± 0.94
Contribution of cell-to-cell infection	$\frac{R_{c c}}{R_{c f} + R_{c c}}$	–	0.50 (0.34–0.63)	0.60 (0.44–0.72)	0.57 (0.43–0.70)	0.57 ± 0.07

*

Mean value.
†

95% confidence interval.
‡

Average and standard deviation of merged values in experiment 1, 2, and 3.

Figure 2 with 3 supplements see all

Download asset Open asset

Dynamics of HIV-1 infection in Jurkat cells through cell-free and cell-to-cell infection.

Jurkat cells were inoculated with HIV-1 (at multiplicity of infection 0.1) in the static and shaking cell cultures. Panels A and B show the time-course of experimental data for the numbers of uninfected cells (top) and infected cells (middle), and the amount of viral protein p24 (bottom) in the static and shaking cell culture systems, respectively. The shadow regions correspond to 95% posterior predictive intervals, the dashed curves give the best-fit solution (mean) for Equations 1–3 to the time-course dataset. All data in each experiment were fitted simultaneously. In panels A and B, the results of three independent experiments are respectively shown.

https://doi.org/10.7554/eLife.08150.005

Our model (i.e., Equations 1–3) applied to time-course 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 cell-free infection (β) and the rate constant for the cell-to-cell 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₀ (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 cell-free and cell-to-cell infection. Interestingly, in spite of nonlinear interaction between the two modes of virus transmission, our derivation of R₀ revealed that the secondly infected cells were the sum of the basic reproduction number through the cell-free infection R_cf = βpT_max/δc and the basic reproduction number through the cell-to-cell infection R_cc = ωT_max/δ, (i.e., R₀ = R_cf + R_cc) (see Figure 1B). Using all accepted MCMC parameter estimates from the time-course experimental datasets, we calculated that on average the mean of the total basic reproductive number is R₀ = 5.83 ± 0.94 (average ± standard deviation), and the mean number of secondly infected cells through the cell-free infection and the cell-to-cell infection are R_cf = 2.44 ± 0.23 and R_cc = 3.39 ± 0.91, respectively (see Table 1). The distributions of calculated R₀, R_cf, and R_cc, are shown in Figure 3A–C, respectively. These estimates indicate that the contribution of the cell-to-cell 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 HIV-1 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.

Figure 3

Download asset Open asset

Distribution of the basic reproduction numbers, generation time, and Malthus coefficient.

The distribution of the basic reproduction number, R₀, the number of secondary infected cells through the cell-free infection, R_cf, and the cell-to-cell infection, R_cc, calculated from all accepted Markov Chain Monte Carlo (MCMC) parameter estimates are shown in A, B, and C, respectively. The contribution of the cell-to-cell infection (i.e., R_cc/(R_cf + R_cc)) is distributed as in D. For each plot, the last 15,000 MCMC samples among the total 50,000 samples are used. a.u., arbitrary unit.

https://doi.org/10.7554/eLife.08150.009

Advantage of cell-to-cell 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 cell-to-cell infection (i.e., for the static and shaking cell cultures, respectively), the mean generation time is calculated as 1/δ + R_cf/cR₀ = 2.22 ± 0.32 days and 1/δ + 1/c = 2.47 ± 0.32 days, respectively (see Table 2). Thus, cell-to-cell infection shortens the generation time by on average 0.90 times, and enables HIV-1 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 cell-to-cell infection, the Malthus coefficient is calculated as 1.86 ± 0.37 and 0.49 ± 0.05 per day, respectively (see Table 2). Thus, cell-to-cell infection increases the HIV-1 fitness by 3.80-fold (corresponding to 944-fold higher viral load 5 days after the infection) and plays an important role in the rapid spread of HIV-1. Thus, the efficient viral spread via the cell-to-cell infection is relevant, especially at the beginning of virus infection.

Table 2

Generation time and Malthus coefficient of virus infection

https://doi.org/10.7554/eLife.08150.010

Cell culture system	Exp. 1	Exp. 2	Exp. 3	Ave. ± S.D.
Generation time of HIV-1
Static cell culture	2.51* days	2.08 days	2.22 days	(2.22 ± 0.32)‡ days
Static cell culture	(1.78–3.38) days	(1.54–2.78) days	(1.69–2.93) days	–
Shaking cell culture	2.73† days	2.34 days	2.47 days	(2.47 ± 0.32) days
Shaking cell culture	(1.99–3.59) days	(1.77–3.06) days	(1.91–3.18) days	–
Malthus coefficient of HIV-1
Static cell culture	1.61 day⁻¹	2.03 day⁻¹	1.86 day⁻¹	(1.86 ± 0.37) day⁻¹
Static cell culture	(1.10–2.27) day⁻¹	(1.32–3.01) day⁻¹	(1.26–2.72) day⁻¹	–
Shaking cell culture	0.57 day⁻¹	0.46 day⁻¹¹	0.49 day⁻¹	(0.49 ± 0.05) day⁻¹
Shaking cell culture	(0.47–0.67) day⁻¹	(0.38–0.56) day⁻¹	(0.39–0.61) day⁻¹	–

*

Mean value.
†

95% confidence interval.
‡

Average and standard deviation of merged values in experiment 1, 2, and 3.
HIV-1, human immunodeficiency virus type 1.

Virtual experiments of cell-to-cell infection in silico

While the shaking culture prevents the cell-to-cell infection, it is technically difficult to completely block the cell-free infection. Here, using our estimated kinetic parameters (Table 1 and Supplementary file 3), we carried out a ‘virtual experiment’ eliminating the contribution of the cell-free infection using all accepted MCMC estimated parameter values, allowing to estimate only the cell-to-cell infection, in silico (see Figure 4). Our simulated mean values (represented by solid lines) of the cell-to-cell infection of HIV-1 are consistently located between the time course of experimental data under the static conditions (closed circles, including both the cell-free and cell to cell infections) and those under the shaking conditions (open circles, reflecting only the cell-free 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 cell-to-cell virus propagation, are closer to experimental data from the coupled cell-free and cell-to-cell infection, than to data from the cell-free infection only. This shows that the cell-free infection, which contributes approximately 40% to the whole HIV-1 infection process, plays a limited role on the virus spread. In other words, even if we could completely block the cell-free infection, the cell-to-cell infection would still effectively spread viruses (Sigal et al., 2011). We address this point in ‘Discussion’.

Figure 4

Download asset Open asset

Simulating cell-to-cell infection of HIV-1.

Using our estimated parameters, the pure cell-to-cell infection is simulated in silico (solid curves). The simulated values are located between the time course of experimental data under the static conditions (closed circles) and those under the shaking conditions (open circles). The shadowed regions correspond to 95% posterior predictive intervals.

https://doi.org/10.7554/eLife.08150.011

Discussion

Through experimental-mathematical investigation, here we quantitatively elucidated the dynamics of the cell-to-cell and cell-free HIV-1 infection modes (Figure 2 and Table 1). We derived the basic reproduction number, R₀, and divided it into the numbers of secondly infected cells through the cell-free infection, R_cf, and the cell-to-cell 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 cell-to-cell 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 cell-to-cell 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 cell-to-cell infection in patients infected with HIV-1, it should be noted that the environment of immune cells including CD4⁺ T-cells 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 HIV-1-infected cells converge to lymph nodes and can be vehicles for viral dissemination in vivo (Murooka et al., 2012). Moreover, certain studies have suggested that cell-to-cell viral spread is resistant to anti-viral immunity such as neutralizing antibodies and cytotoxic T lymphocytes (Martin and Sattentau, 2009). Therefore, these notions strongly suggest that the contribution of the cell-to-cell infection for viral propagation in vivo may be much higher than that estimated from the in vitro cell culture system.

As another significance of cell-to-cell viral spread, Sigal et al. have suggested that the cell-to-cell infection permits viral replication even under the anti-retroviral 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 cell-to-cell infection remained unclear. To further understand the role of the cell-to-cell infection, we quantified the contributions of the cell-to-cell and cell-free infection modes (Table 1). Interestingly, we found that the cell-to-cell infection mode is predominant during the infection. Furthermore, our virtual experiments showed that a complete block of the cell-free infection, which is highly susceptible to current antiviral drugs, provides only a limited impact on the whole HIV-1 infection (Figure 3). Taken together, our findings further support that the cell-to-cell infection can be a barrier to prevent the cure of HIV-1 infection, which is discussed in Sigal et al. (2011). However, it should be noted that some papers have shown that cell-to-cell spread cannot overcome the action of most anti-HIV-1 drugs (Titanji et al., 2013; Agosto et al., 2014). To fully elucidate this issue, further investigations will be needed.

In addition to HIV-1, other viruses such as herpes simplex virus, measles virus, and human hepatitis C virus drive their dissemination via cell-to-cell infection (Sattentau, 2008; Talbert-Slagle et al., 2014). Although the impact of cell-to-cell 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 cell-to-cell infection along with producing infectious viral particles. By applying mathematical modeling to the experimental data, here we estimated the sole dynamics of cell-free 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 HIV-1 infection

Request a detailed protocol

Jurkat 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 pNL4-3 plasmid (Adachi et al., 1986) (GenBank accession no. M19921.2) was transfected into 293T cells by the calcium-phosphate method. At 48 hr post-transfection, the culture supernatant was harvested, centrifuged, and then filtered through a 0.45-μm-pore-size 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 phytohemagglutinin-stimulated human peripheral blood mononuclear cells in a 96-well plate in triplicate. At 14 days postinfection, the endpoint was determined by using an HIV-1 p24 antigen enzyme-linked 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₅₀) according to the Reed-Muench method.

For HIV-1 infection, 3 × 10⁵ of Jurkat cells were infected with HIV-1 (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 25-cm² flask (Nunc, Rochester, NY). For the static infection, the infected cell culture was kept in a 37°C/5% CO₂ 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 Model-2230 (Wakenyaku, Japan) placed in 37°C/5% CO₂ 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 HIV-1 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 isothiocyanate-labeled anti-HIV-1 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 HIV-1 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 HIV-1 inoculation). The detection threshold of each value are the followings: cell number (cell counting), 3000 cells/ml; % p24-positive cells (flow cytometry), 0.3%; and p24 antigen in culture supernatant (p24 antigen ELISA), 80 pg/ml.

Parameter estimation

Request a detailed protocol

A 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 protocol

We here estimate the growth kinetics of Jurkat cells, which have been commonly used for HIV-1 studies, under the normal (i.e., mock-infected) condition with the following mathematical model:

\frac{d T (t)}{d t} = g T (t) (1 - \frac{T (t)}{T_{\max}}),

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 least-squares regression (FindMinimum package of Mathematica9.0) was performed to fit Equation 4 to the time-course numbers of Jurkat cells in the normal condition. The fitted parameter values are listed in Supplementary file 2 and the model behavior using these best-fit parameter estimates is presented together with the data in Figure 2—figure supplement 2.

Mathematical appendix

Request a detailed protocol

The linearized equation of Equations 1–3 at the virus-free steady state, (T_max, 0, 0), is given as follows:

\frac{d I (t)}{d t} = β T_{\max} V (t) + ω T_{\max} I (t) - δ I (t),

\frac{d V (t)}{d t} = p I (t) - c V (t) .

Let b(t) be the number of newly produced infected cells in the linear phase:

b (t) : = β T_{\max} V (t) + ω T_{\max} I (t) .

Applying the variation of constants formula to Equations 5, 6, we have

V (t) = V (0) e^{- c t} + \int_{0}^{t} e^{- c (t - s)} p I (s) d s,

I (t) = I (0) e^{- δ t} + \int_{0}^{t} e^{- δ (t - z)} b (z) d z .

Inserting Equation 9 into Equation 8 to exchange the order of integrals, we have

V (t) = g (t) + p \int_{0}^{t} \int_{0}^{x} e^{- c (x - θ) - δ θ} d θ b (t - x) d x,

where

g (t) ∶ = V (0) e^{- c t} + \int_{0}^{t} e^{- c (t - s)} p I (0) e^{- δ t} d s .

From Equation 7 and Equations 9, 10, we arrive at the following renewal equation:

b (t) = h (t) + \int_{0}^{t} Ψ (x) b (t - x) d x,

where h(t) is given by

h (t) ∶ = ω T_{m a x} I (0) e^{- δ t} + β T_{m a x} g (t),

and the kernel Ψ(x) is given by

Ψ (x) ∶ = β T_{m a x} p \int_{0}^{x} e^{- δ θ - c (x - θ)} d θ + ω T_{m a x} e^{- δ x},

= \frac{β T_{m a x} p}{δ c} (ϕ_{1} * ϕ_{2}) (x) + \frac{ω T_{m a x}}{δ} ϕ_{1} (x) .

In the above expression, ϕ_j(x) denotes the probability density function given by

ϕ_{1} (x) = δ e^{- δ x}, ϕ_{2} (x) = c e^{- c x},

and, ∗ denotes the convolution of functions. From the general theory of the basic reproduction number (Inaba, 2012), R₀ for the reproduction of infected cells is given by

R_{0} = \int_{0}^{\infty} Ψ (x) d x = \frac{β T_{m a x} p}{δ c} + \frac{ω T_{m a x}}{δ} = R_{c f} + R_{c c},

where R_cf and R_cc denote the reproduction numbers for infected cells mediated by the cell-free and cell-to-cell 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

ρ (t) = p I (0) e^{- δ t} + \int_{0}^{t} e^{- δ (t - z)} (β T_{m a x} p V (z) + ω T_{m a x} ρ (z)) d z,

where

V (t) = V (0) e^{- c t} + \int_{0}^{t} e^{- c (t - s)} ρ (s) d s .

Inserting Equation 11 into Equation 12, we again arrive at the following renewal equation:

ρ (t) = q (t) + \int_{0}^{t} Ψ (x) ρ (t - x) d x,

where q(t) is given by

q (t) ∶ = p I (0) e^{- δ t} + \int_{0}^{t} e^{- δ (t - z)} p β T_{m a x} V (0) e^{- c z} d z .

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

ψ (x) ∶ = \frac{Ψ (x)}{R_{0}} = \frac{R_{c f}}{R_{0}} (ϕ_{1} * ϕ_{2}) (x) + \frac{R_{c c}}{R_{0}} ϕ_{1} (x) .

Then the generation time for the virus reproduction, denoted by G, is calculated as follows:

G ∶ = \int_{0}^{\infty} t ψ (t) d t = \frac{R_{c f}}{R_{0}} G_{c f} + \frac{R_{c c}}{R_{0}} G_{c c} \leq G_{c f},

where G_cf : = 1/δ + 1/c and G_cc : = 1/δ are the generation times for virus reproduction mediated by the cell-free and cell-to-cell infection, respectively.

The Malthusian coefficient for the virus reproduction must be given as the dominant real root of the Euler-Lotka equation as

\int_{0}^{\infty} e^{- λ x} Ψ (x) d x = \frac{β T_{m a x} p}{δ c} \hat{ϕ_{1}} (λ) \hat{ϕ_{2}} (λ) + \frac{ω T_{m a x}}{δ} \hat{ϕ_{1}} (λ) = 1,

where $\hat{ϕ_{j}}$ denotes the Laplace transformation of a function ϕ_j. That is,

\hat{ϕ_{1}} (λ) = \int_{0}^{\infty} e^{- λ x} ϕ_{1} (x) d s = \frac{δ}{δ + λ}, \hat{ϕ_{2}} (λ) = \int_{0}^{\infty} e^{- λ x} ϕ_{2} (x) d s = \frac{c}{c + λ} .

Therefore the Euler-Lotka equation can be calculated explicitly as follows:

\frac{β T_{m a x} p}{δ c} \frac{δ c}{(δ + λ) (c + λ)} + \frac{ω T_{m a x}}{δ} \frac{δ}{δ + λ} = 1,

which is reduced to a quadratic equation,

λ^{2} + δ c (G_{c c} + (1 - R_{c c}) (G_{c f} - G_{c c})) λ + δ c (1 - R_{0}) = 0 .

If R₀ > 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,

λ = \frac{- δ c (G_{c c} + (1 - R_{c c}) (G_{c f} - G_{c c})) + \sqrt{δ^{2} c^{2} {(G_{c c} + (1 - R_{c c}) (G_{c f} - G_{c c}))}^{2} - 4 δ c (1 - R_{0})}}{2} .

References

1. Adachi A
2. Gendelman HE
3. Koenig S
4. Folks T
5. Willey R
6. Rabson A
7. Martin MA
(1986)
Production of acquired immunodeficiency syndrome-associated retrovirus in human and nonhuman cells transfected with an infectious molecular clone

Journal of Virology 59:284–291.
- Google Scholar
1. Agosto LM
2. Zhong P
3. Munro J
4. Mothes W
(2014) Highly active antiretroviral therapies are effective against HIV-1 cell-to-cell transmission
PLOS Pathogens 10:e1003982.

https://doi.org/10.1371/journal.ppat.1003982
- Google Scholar
1. Carr JM
2. Hocking H
3. Li P
4. Burrell CJ
(1999) Rapid and efficient cell-to-cell transmission of human immunodeficiency virus infection from monocyte-derived macrophages to peripheral blood lymphocytes
Virology 265:319–329.

https://doi.org/10.1006/viro.1999.0047
- Google Scholar
1. Chen P
2. Hübner W
3. Spinelli MA
4. Chen BK
(2007) Predominant mode of human immunodeficiency virus transfer between T cells is mediated by sustained Env-dependent neutralization-resistant virological synapses
Journal of Virology 81:12582–12595.

https://doi.org/10.1128/JVI.00381-07
- Google Scholar
1. Dimitrov DS
2. Willey RL
3. Sato H
4. Chang LJ
5. Blumenthal R
6. Martin MA
(1993)
Quantitation of human immunodeficiency virus type 1 infection kinetics

Journal of Virology 67:2182–2190.
- Google Scholar
1. Fukuhara M
2. Iwami S
3. Sato K
4. Nishimura Y
5. Shimizu H
6. Aihara K
7. Koyanagi Y
(2013) Quantification of the dynamics of enterovirus 71 infection by experimental-mathematical investigation
Journal of Virology 87:701–705.

https://doi.org/10.1128/JVI.01453-12
- Google Scholar
1. Inaba H
(2012) On a new perspective of the basic reproduction number in heterogeneous environments
Journal of Mathematical Biology 65:309–348.

https://doi.org/10.1007/s00285-011-0463-z
- Google Scholar
1. Iwami S
2. Holder BP
3. Beauchemin CA
4. Morita S
5. Tada T
6. Sato K
7. Igarashi T
8. Miura T
(2012a) Quantification system for the viral dynamics of a highly pathogenic simian/human immunodeficiency virus based on an in vitro experiment and a mathematical model
Retrovirology 9:18.

https://doi.org/10.1186/1742-4690-9-18
- Google Scholar
1. Iwami S
2. Sato K
3. De Boer RJ
4. Aihara K
5. Miura T
6. Koyanagi Y
(2012b) Identifying viral parameters from in vitro cell cultures
Frontiers in Microbiology 3:319.

https://doi.org/10.3389/fmicb.2012.00319
- Google Scholar
1. Kakizoe Y
2. Morita S
3. Nakaoka S
4. Takeuchi Y
5. Sato K
6. Miura T
7. Beauchemin CA
8. Iwami S
(2015) A conservation law for virus infection kinetics in vitro
Journal of Theoretical Biology 376:39–47.

https://doi.org/10.1016/j.jtbi.2015.03.034
- Google Scholar
(2013a) Relative contribution of free-virus and synaptic transmission to the spread of HIV-1 through target cell populations
Biology Letters 9:20121049.

https://doi.org/10.1098/rsbl.2012.1049
- Google Scholar
(2013b) Synaptic transmission and the susceptibility of HIV infection to anti-viral drugs
Scientific Reports 3:2103.

https://doi.org/10.1038/srep02103
- Google Scholar
1. Komarova NL
2. Wodarz D
(2013) Virus dynamics in the presence of synaptic transmission
Mathematical Biosciences 242:161–171.

https://doi.org/10.1016/j.mbs.2013.01.003
- Google Scholar
1. Martin N
2. Sattentau Q
(2009) Cell-to-cell HIV-1 spread and its implications for immune evasion
Current Opinion in HIV and AIDS 4:143–149.

https://doi.org/10.1097/COH.0b013e328322f94a
- Google Scholar
1. Monel B
2. Beaumont E
3. Vendrame D
4. Schwartz O
5. Brand D
6. Mammano F
(2012) HIV cell-to-cell transmission requires the production of infectious virus particles and does not proceed through env-mediated fusion pores
Journal of Virology 86:3924–3933.

https://doi.org/10.1128/JVI.06478-11
- Google Scholar
1. Murooka TT
2. Deruaz M
3. Marangoni F
4. Vrbanac VD
5. Seung E
6. von Andrian UH
7. Tager AM
8. Luster AD
9. Mempel TR
(2012) HIV-infected T cells are migratory vehicles for viral dissemination
Nature 490:283–287.

https://doi.org/10.1038/nature11398
- Google Scholar
Book
1. Nowak M
2. May RM
(2000)
Virus dynamics: mathematical principles of immunology and virology: mathematical principles of immunology and virology

Oxford University Press.
- Google Scholar
Book
1. Nowak MA
(2006)
Evolutionary dynamics

Harvard University Press.
- Google Scholar
1. Perelson AS
2. Nelson PW
(1999) Mathematical analysis of HIV-1 dynamics in vivo
SIAM Review 41:3–44.

https://doi.org/10.1137/S0036144598335107
- Google Scholar
1. Perelson AS
(2002) Modelling viral and immune system dynamics
Nature Reviews Immunology 2:28–36.

https://doi.org/10.1038/nri700
- Google Scholar
(1992) Cell-to-cell spread of HIV-1 occurs within minutes and may not involve the participation of virus particles
Virology 186:712–724.

https://doi.org/10.1016/0042-6822(92)90038-Q
- Google Scholar
1. Sato K
2. Izumi T
3. Misawa N
4. Kobayashi T
5. Yamashita Y
6. Ohmichi M
7. Ito M
8. Takaori-Kondo A
9. Koyanagi Y
(2010) Remarkable lethal G-to-A mutations in vif-proficient HIV-1 provirus by individual APOBEC3 proteins in humanized mice
Journal of Virology 84:9546–9556.

https://doi.org/10.1128/JVI.00823-10
- Google Scholar
1. Sato K
2. Misawa N
3. Nie C
4. Satou Y
5. Iwakiri D
6. Matsuoka M
7. Takahashi R
8. Kuzushima K
9. Ito M
10. Takada K
11. Koyanagi Y
(2011) A novel animal model of Epstein-Barr virus-associated hemophagocytic lymphohistiocytosis in humanized mice
Blood 117:5663–5673.

https://doi.org/10.1182/blood-2010-09-305979
- Google Scholar
1. Sato K
2. Misawa N
3. Iwami S
4. Satou Y
5. Matsuoka M
6. Ishizaka Y
7. Ito M
8. Aihara K
9. An DS
10. Koyanagi Y
(2013) HIV-1 Vpr accelerates viral replication during acute infection by exploitation of proliferating CD4+ T cells in vivo
PLOS Pathogens 9:e1003812.

https://doi.org/10.1371/journal.ppat.1003812
- Google Scholar
1. Sato K
2. Takeuchi JS
3. Misawa N
4. Izumi T
5. Kobayashi T
6. Kimura Y
7. Iwami S
8. Takaori-Kondo A
9. Hu WS
10. Aihara K
11. Ito M
12. An DS
13. Pathak VK
14. Koyanagi Y
(2014) APOBEC3D and APOBEC3F potently promote HIV-1 diversification and evolution in humanized mouse model
PLOS Pathogens 10:e1004453.

https://doi.org/10.1371/journal.ppat.1004453
- Google Scholar
1. Sattentau Q
(2008) Avoiding the void: cell-to-cell spread of human viruses
Nature Reviews Microbiology 6:815–826.

https://doi.org/10.1038/nrmicro1972
- Google Scholar
1. Sigal A
2. Kim JT
3. Balazs AB
4. Dekel E
5. Mayo A
6. Milo R
7. Baltimore D
(2011) Cell-to-cell spread of HIV permits ongoing replication despite antiretroviral therapy
Nature 477:95–98.

https://doi.org/10.1038/nature10347
- Google Scholar
(2007) Inefficient human immunodeficiency virus replication in mobile lymphocytes
Journal of Virology 81:1000–1012.

https://doi.org/10.1128/JVI.01629-06
- Google Scholar
1. Talbert-Slagle K
2. Atkins KE
3. Yan KK
4. Khurana E
5. Gerstein M
6. Bradley EH
7. Berg D
8. Galvani AP
9. Townsend JP
(2014) Cellular superspreaders: an epidemiological perspective on HIV infection inside the body
PLOS Pathogens 10:e1004092.

https://doi.org/10.1371/journal.ppat.1004092
- Google Scholar
(2013) Protease inhibitors effectively block cell-to-cell spread of HIV-1 between T cells
Retrovirology 10:161.

https://doi.org/10.1186/1742-4690-10-161
- Google Scholar
1. Watanabe T
2. Urano E
3. Miyauchi K
4. Ichikawa R
5. Hamatake M
6. Misawa N
7. Sato K
8. Ebina H
9. Koyanagi Y
10. Komano J
(2012) The hematopoietic cell-specific Rho GTPase inhibitor ARHGDIB/D4GDI limits HIV type 1 replication
AIDS Research and Human Retroviruses 28:913–922.

https://doi.org/10.1089/AID.2011.0180
- Google Scholar
1. Zhong P
2. Agosto LM
3. Ilinskaya A
4. Dorjbal B
5. Truong R
6. Derse D
7. Uchil PD
8. Heidecker G
9. Mothes W
(2013) Cell-to-cell transmission can overcome multiple donor and target cell barriers imposed on cell-free HIV
PLOS ONE 8:e53138.

https://doi.org/10.1371/journal.pone.0053138
- Google Scholar

Article and author information

Author details

Shingo Iwami
1. Mathematical Biology Laboratory, Department of Biology, Faculty of Sciences, Kyushu University, Fukuoka, Japan
2. PRESTO, Japan Science and Technology Agency, Saitama, Japan
3. CREST, Japan Science and Technology Agency, Saitama, Japan
Contribution
SI, Conception and design, Analysis and interpretation of data, Drafting or revising the article, Contributed unpublished essential data or reagents

Contributed equally with
Junko S Takeuchi

For correspondence
siwami@kyushu-u.org

Competing interests
The authors declare that no competing interests exist.
Junko S Takeuchi

Laboratory of Viral Pathogenesis, Institute for Virus Research, Kyoto University, Kyoto, Japan

Contribution
JST, Conception and design, Acquisition of data, Contributed unpublished essential data or reagents

Contributed equally with
Shingo Iwami

Competing interests
The authors declare that no competing interests exist.
Shinji Nakaoka

Graduate School of Medicine, University of Tokyo, Tokyo, Japan

Contribution
SN, Analysis and interpretation of data, Drafting or revising the article

Competing interests
The authors declare that no competing interests exist.
Fabrizio Mammano
1. INSERM-Genetics and Ecology of viruses, Hospital Saint Louis, Paris, France
2. Université Paris Diderot, Sorbonne Paris Cité, Paris, France
Contribution
FM, Conception and design, Analysis and interpretation of data, Drafting or revising the article, Contributed unpublished essential data or reagents

Competing interests
The authors declare that no competing interests exist.

"This ORCID iD identifies the author of this article:" 0000-0002-9193-7696
François Clavel
1. INSERM-Genetics and Ecology of viruses, Hospital Saint Louis, Paris, France
2. Université Paris Diderot, Sorbonne Paris Cité, Paris, France
Contribution
FC, Conception and design, Contributed unpublished essential data or reagents

Competing interests
The authors declare that no competing interests exist.
Hisashi Inaba

Graduate School of Mathematical Sciences, University of Tokyo, Tokyo, Japan

Contribution
HI, Analysis and interpretation of data, Drafting or revising the article

Competing interests
The authors declare that no competing interests exist.
Tomoko Kobayashi

Laboratory for Animal Health, Department of Animal Science, Faculty of Agriculture, Tokyo University of Agriculture, Kanagawa, Japan

Contribution
TK, Acquisition of data, Contributed unpublished essential data or reagents

Competing interests
The authors declare that no competing interests exist.
Naoko Misawa

Laboratory of Viral Pathogenesis, Institute for Virus Research, Kyoto University, Kyoto, Japan

Contribution
NM, Acquisition of data, Contributed unpublished essential data or reagents

Competing interests
The authors declare that no competing interests exist.
Kazuyuki Aihara
1. Institute of Industrial Science, University of Tokyo, Tokyo, Japan
2. Graduate School of Information Science and Technology, University of Tokyo, Tokyo, Japan
Contribution
KA, Conception and design, Analysis and interpretation of data, Drafting or revising the article

Competing interests
The authors declare that no competing interests exist.
Yoshio Koyanagi

Laboratory of Viral Pathogenesis, Institute for Virus Research, Kyoto University, Kyoto, Japan

Contribution
YK, Conception and design, Drafting or revising the article, Contributed unpublished essential data or reagents

Competing interests
The authors declare that no competing interests exist.
Kei Sato
1. CREST, Japan Science and Technology Agency, Saitama, Japan
2. Laboratory of Viral Pathogenesis, Institute for Virus Research, Kyoto University, Kyoto, Japan
Contribution
KS, Conception and design, Acquisition of data, Analysis and interpretation of data, Drafting or revising the article, Contributed unpublished essential data or reagents

For correspondence
ksato@virus.kyoto-u.ac.jp

Competing interests
The authors declare that no competing interests exist.

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); Grants-in-Aid 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 World-Leading 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 (H27-ShinkoJitsuyoka-General-016) (to SI, SN, HI, KA, and KS); Agence Nationale de Recherches sur le Sida et les Hepatites Virales (ANRS) (to FM); a Grant-in-Aid for Scientific Research on Innovative Areas 24115008 from the Ministry of Education, Culture, Sports, Science and Technology of Japan (to YK); JSPS Core-to-Core 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).

Copyright

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.