Improving the estimation of the death rate of infected cells from time course data during the acute phase of virus infections: application to acute HIV-1 infection in a humanized mouse model
- Hiroki Ikeda^{1},
- Rob J de Boer^{2},
- Kei Sato^{3},
- Satoru Morita^{4},
- Naoko Misawa^{3},
- Yoshio Koyanagi^{3},
- Kazuyuki Aihara^{5, 6} and
- Shingo Iwami^{1, 7}Email author
https://doi.org/10.1186/1742-4682-11-22
© Ikeda et al.; licensee BioMed Central Ltd. 2014
Received: 13 February 2014
Accepted: 13 May 2014
Published: 21 May 2014
Abstract
Background
Mathematical modeling of virus dynamics has provided quantitative insights into viral infections such as influenza, the simian immunodeficiency virus/human immunodeficiency virus, hepatitis B, and hepatitis C. Through modeling, we can estimate the half-life of infected cells, the exponential growth rate, and the basic reproduction number (R_{0}). To calculate R_{0} from virus load data, the death rate of productively infected cells is required. This can be readily estimated from treatment data collected during the chronic phase, but is difficult to determine from acute infection data. Here, we propose two new models that can reliably estimate the average life span of infected cells from acute-phase data, and apply both methods to experimental data from humanized mice infected with HIV-1.
Methods
Both new models, called as the reduced quasi-steady state (RQS) model and the piece-wise regression (PWR) model, are derived by simplification of a standard model for the acute-phase dynamics of target cells, viruses and infected cells. By having only a limited number of parameters, both models allow us to reliably estimate the death rate of productively infected cells. Simulated datasets with plausible parameter values are generated with the standard model to compare the performance of the new models with that of the major previous model (i.e., the simple exponential model). Finally, we fit models to time course data from HIV-1 infected humanized mice to estimate the several important parameters characterizing their acute infection.
Results and conclusions
The new models provided much better estimates than the previous model because they more precisely capture the de novo infection process. Both models describe the acute phase of HIV-1 infected humanized mice reasonably well, and we estimated an average death rate of infected cells of 0.61 and 0.61, an average exponential growth rate of 0.69 and 0.76, and an average basic reproduction number of 2.30 and 2.38 in the RQS model and the PWR model, respectively. These estimates are fairly close to those obtained in humans.
Keywords
Background
In most viral infections, the initial exponential growth phase is followed by a second exponential phase known as contraction. In hosts exposed to viruses such as influenza and coronaviruses (the causative agents of severe acute respiratory syndrome), the viral load continuously declines during the contraction phase [1, 2]. In contrast, in chronic viral infections, such as human immunodeficiency virus (HIV) and hepatitis C virus (HCV) infections, contraction slows down such that the viral load approaches a steady state, called the virological set point [3, 4]. In both infection types, the expansion and contraction of the viral load have been modeled as single exponential functions, with parameters determined by linear regression of the log transformed data [3, 5–7]. This simple approach is reasonable as long as the conditions, e.g., the availability of target cells or the immune response, hardly change within each phase. Using this approach, the initial growth rate, death rate of the infected cells, and the basic reproduction number (i.e., R_{0}) have been estimated, which has improved our understanding of particular virus infections, and has guided medical treatment [7–9]. For example, once the basic reproduction number is estimated, the critical inhibition, 1–1/R_{0}, induced by vaccines, or by antiviral drugs, to prevent primary virus infection can be calculated [6].
Knowledge of the death rate of infected cells is crucial for properly understanding viral dynamics because the average life time is required for calculating the basic reproduction number. In chronic viral infections such as HIV and HCV, the death rate is estimated from large perturbations of the set point viral load data instigated by potent anti-viral therapy [4, 10–13]. Shortly after effective treatment, the decay rate of viral load approaches the death rate of productively infected cells. However, estimating the death rate of infected cells during the acute phase remains a challenging task, and, as a consequence, it is difficult to accurately estimate the basic reproduction number from viral load data during the early stages of viral infection. In addition to calculation of the basic reproduction number, the death rate per se is useful for evaluating the efficacy of vaccine-induced cellular immune responses during the acute phase of virus infection [14–16]. Therefore, an improved method for estimating the death rate of infected cells during this phase is urgently required.
In this study, we first generated simulated datasets with biologically plausible parameter values, using a population dynamics model of virus population. The time evolution of target cell densities and viral load were modeled during the acute phase. The datasets describing acute infection were subsequently analyzed by two novel mathematical models to evaluate whether the new models could accurately estimate the known parameters. Our proposed models properly described the artificial datasets and delivered better estimates of the parameters and well calculated indices than conventional models (i.e., simple exponential models). Our methods proved especially effective for calculating the death rate of infected cells. We then applied our models to time course data from a human hepatopoietic stem cell-transplanted humanized mouse model infected with HIV type-1 (HIV-1) [17–20], to quantify the infection dynamics during the acute phase. To our knowledge, this is the first report quantifying the dynamics of acute HIV-1 infection in humanized mice. Finally, we discuss how our approach may be combined with animal experiments. Like previous simple exponential models [3, 5–7], our approach is quite general and can be used in several infection models.
Methods
Mathematical models describing the acute phase of virus infection
(I) Reduced quasi-steady state (RQS) model
where the parameters β, δ, p and c represent the conventional rate constants for viral infection of target cells, the death rate of infected cells, the virus production rate in an infected cell, and the clearance rate of virus particles, respectively. The initial expansion of viral load in this model is well approximated by V(t) ≈ V(0) exp (g_{0}t) [3, 5–7], with an exponential growth rate, g_{0} (the Malthusian parameter), given by the positive root of the characteristic equation g_{0}^{2} + (c + δ)g_{0} + cδ - pβT(0) = 0, i.e., ${\mathit{g}}_{0}=\left\{-\left(\mathit{c}+\mathit{\delta}\right)+\sqrt{{\left(\mathit{c}-\mathit{\delta}\right)}^{2}+4\mathit{p}\mathit{\beta}\mathit{T}\left(0\right)}\right\}/2$.
where r = pβ/c is the viral replication rate per target cell, and δ is the death rate of infected cells. Eqs. (1) and (4) together form our first model, that we here call the “reduced quasi-steady state” (RQS) model. The RQS model lumps the 8 parameters of the reduced standard model of Eqs. (1–3) into five parameters, i.e., β, r, δ, T(0), and V(0). Because there is no production of target cells the infection will ultimately be cleared. The five parameters together define several “observables”. First, the basic reproduction number is R_{0} = rT(0)/δ. Second, the initial exponential growth rate is g_{0} ≈ rT(0) - δ ( when δ and g_{0} ≪ c, one can directly calculate the same g_{0} from the characteristic equation). Third, the final level of target cells is given by the epidemiological “final size equation” [24, 25], as the solution of f = exp [ - R_{0}(1 - f)], where f = T_{ min }/T(0) is the fraction of surviving target cells.
(II) Piece wise targets (PWT) model
(III) Piece wise regression (PWR) model
As a control method, we additionally consider a classical method that has been widely adopted in earlier literatures [3, 5–7]. Previously acute infection data have been quantified using piece wise linear regression of the log transformed viral loads before and after the peak in the viral load [3, 5–7]. The ascending and descending slopes of the log viral load roughly correspond to the exponential growth rate, g_{0}, and the death rate, δ, of infected cells, respectively [3, 5–7]. Since g_{0} = rT(0) - δ and R_{0} = rT(0)/δ, these two slopes suffice to estimate the basic reproduction number, R_{0} ≈ 1 + g_{0}/δ, and knowing the initial target cell density, T(0), the viral replication rate per target cell can be estimated from the R_{0}, i.e., r = δR_{0}/T(0) [3, 5–7]. Here we call this classical model as the “piece wise regression” model (PWR). It has been realized before that the down-slope will only reflect the death rate of infected cells, δ, if there is hardly any residual infection of target cells during the contraction phase, i.e., if target cells are markedly depleted [3, 5–7].
Artificial data generated from a population dynamics model of virus infection
To estimate the accuracy of the parameters estimated by our two novel models, we created simulated time course data of target cell densities and viral load during the acute phase of viral infection (lasting approximately 21 days [14–20, 26–28]) assuming biologically plausible parameter values. The artificial datasets were generated with the reduced standard population dynamics model of viral infection, i.e., Eqs. (1–3), in which the target cell dynamics are coupled to the viral dynamics by the infection term.
We added stochastic variation to the “data” generated by this model by adding “observational” noise and/or by varying the parameter values. The log transformed data were perturbed by adding a normally distributed noise variable with zero mean and standard deviation σ (see Results). The datasets describing acute infection were subsequently analyzed by the two novel RQS and PWT models, and by the previous PWR model.
HIV-1 infection in humanized mice
The dynamics of HIV-1 infection during acute infection were quantified in a human hepatopoietic stem cell-transplanted humanized mouse model (NOG-hCD34 mice) [17–20]. Five humanized mice were infected with the CCR5-tropic HIV-1 (strain AD8) [29], and 100 μl of peripheral blood (PB) was routinely collected under anesthesia through the retro-orbital venous plexus at 0, 3, 7, 14, and 21 days post-infection, as previously described [17–20]. The amount of viral RNA in 50 μl of plasma was quantified by RT-PCR (Bio Medical Laboratories, Inc). To estimate target cell densities, the number of memory CD4^{+} T cells was measured by hematometry and flow cytometry, as previously described [17–20]. Briefly, the number of human leukocytes in 10 μl of peripheral blood (PB) was measured in a Celltac α MEK-6450 hematology analyzer (Nihon Kohden, Co.), and the percentage of memory CD4^{+} T cells in human CD45^{+} leukocytes (i.e., CD45^{+} CD3^{+} CD4^{+} CD45RA^{-} cells) was quantified in a FACSCanto II (BD Biosciences) flow cytometer. In the flow cytometry analyses, APC-conjugated anti-CD4 antibody (RPA-4; Biolegend), APC-Cy7-conjugated anti-CD3 antibody (HIT3a; Biolegend), and PE-conjugated anti-CD45 antibody (HI30; Biolegend) were used.
All protocols involving human subjects were reviewed and approved by the Kyoto University institutional review board. Informed written consent from the human subjects was obtained in this study.
Results
Coverage probability of the mathematical models
In the Methods section we formulate two novel mathematical models describing the target cell densities and the viral load during acute infection. We created artificial data with target cell densities and virus loads during acute infection using the reduced standard model for viral infection (i.e., Eqs. (1–3)). The data was generated for one ml of PB with “typical” values of the parameters for HIV-1, i.e., an infection rate β = 10^{-8}per cell per day, a virus production rate p = 4000 particles per day, a death rate of infected cells δ = 0.7 per day, and a clearance rate of c = 23 per day [7–11]. We study whether our simplified models can describe the in silico data, and whether their (lumped) parameters are identifiable. The major biological observables of this model are the initial viral growth rate, ${\mathit{g}}_{0}=\left\{-\left(\mathit{c}+\mathit{\delta}\right)+\sqrt{{\left(\mathit{c}-\mathit{\delta}\right)}^{2}+4\mathit{p}\mathit{\beta}\mathit{T}\left(0\right)}\right\}/2$, the viral replication rate per target cell, r = pβ/c, the death rate of infected cells, δ, and the basic reproduction number, R_{0} = rT(0)/δ.
The coverage probability of the PWR model, the RQS model, and the PWT model
Method | r(viral replication rate) | g_{0}(initial growth rate) | δ(death rate) | R_{0}(basic reproduction number) |
---|---|---|---|---|
Linear regression with the PWR model | 0.25 | 0.87 | 0.12 | 0.25 |
Nonlinear fitting with the RQS model | 0.86 | 0.83 | 0.72 | 0.78 |
Nonlinear fitting with the PWT model | 0.31 | 0.58 | 0.19 | 0.33 |
Sensitivity of parameter estimation to target cell depletion
At the beginning of the infection all models are identical as they all predict exponential growth of the virus load. The models differ around the peak because the previous PWR model assumes an exponential contraction after the peak, whereas the new models allow the peak be formed by the loss of target cells (i.e., by the βT(t)V(t) term in the RQS model, and by the exponential loss of target cells in the PWT model). In both models this loss of target cells continues during the contraction phase. However, if target cells were depleted rapidly such that there would be hardly any infection of target cells during the contraction phase, this difference among the models would vanish, and the contraction phase of the new models would also be dominated by the death rate of infected cells.
As discussed previously [3, 5–7], the PWR model (orange symbols) fails to correctly estimate the death rate of infected cells in the presence of continued de novo infections, i.e., when many target cells survive; see Figure 3A (in which the dashed line denotes the correct death rate at δ = 0.7).
Conversely, the two novel models (RQS: blue symbols and PWT: green symbols) accurately estimate the death rate even if many target cells survive. Note that it is not surprising that the RQS model provides better parameter estimates than the PWT model, because the artificial data was generated by the reduced standard model (Eqs. (1–3)), from which the RQS model was derived by a reasonable quasi-steady state assumption. The parameters are accurately estimated by all methods when target cells are severely depleted (i.e., when f approaches 0), reconfirming that the previous PWR model can accurately quantify the infected cell death when the target cells are severely depleted. This, for example, occurs in CXCR4-tropic simian–human immunodeficiency virus (SHIV) infections (which deplete naïve and memory CD4 T cells during the acute phase) [14, 16, 26–28]. In cases where target cells are not as strongly depleted, the new models do much better than the PWR model in estimating the death rate (Figure 3A), the viral replication rate (Figure 3B) and the basic reproduction number (Figure 3C). All models perform similarly on estimating the exponential growth rate (Figure 3D), because they all assume (implicitly or explicitly) that the number of target cells remains constant during the earliest phases of viral infections.
Application to acute HIV-1 infection in a humanized mouse model
Estimated parameter values and quantities derived from humanized mouse experiments by the novel RQS model
ID | T(0) | V(0) | g _{0} | r | β | δ | R _{0} |
---|---|---|---|---|---|---|---|
cells/ml | RNA copies/ml | day^{-1} | (cell/ml)^{-1}・day^{-1}(×10^{-6}) | (virion /ml)^{-1}・day^{-1}(×10^{-6}) | day^{-1} | --- | |
1 | 274139 | 2461 | 0.467 | 4.46 | 0.10 | 0.76 | 1.62 |
2 | 489096 | 652.1 | 1.065 | 3.37 | 0.03 | 0.58 | 2.83 |
3 | 93238.7 | 717.1 | 0.856 | 12.4 | 1.24 | 0.30 | 3.88 |
4 | 15847.5 | 1233 | 0.429 | 74.1 | 0.57 | 0.75 | 1.58 |
5 | 130580 | 1506 | 0.652 | 10.1 | 0.06 | 0.67 | 1.97 |
Mean | 200580 | 1313 | 0.694 | 20.9 | 0.40 | 0.61 | 2.38 |
Estimated parameter values and quantities derived from humanized mouse experiments by the novel PWT model
ID | T(0) | V(0) | g _{0} | r | t ^{*} | ∆ | δ | R _{0} |
---|---|---|---|---|---|---|---|---|
cells/ml | RNA copies/ml | day^{-1} | (cell/ml)^{-1}・day^{-1}(×10^{-6}) | day | day^{-1} | day^{-1} | --- | |
1 | 282809 | 2325 | 0.48 | 3.91 | 8.7 | 0.07 | 0.62 | 1.77 |
2 | 308636 | 589.7 | 1.06 | 5.86 | 5.8 | 0.13 | 0.75 | 2.42 |
3 | 79882.2 | 693.9 | 0.87 | 15.7 | 3.6 | 0.15 | 0.38 | 3.26 |
4 | 19465.0 | 1013 | 0.50 | 53.9 | 5.6 | 0.07 | 0.55 | 1.90 |
5 | 143159 | 774.8 | 0.87 | 11.4 | 5.1 | 0.08 | 0.76 | 2.14 |
Mean | 166790 | 1079 | 0.76 | 18.2 | 5.8 | 0.10 | 0.61 | 2.30 |
Discussion
We here propose two novel models to quantify the most important parameters characterizing acute viral infections. Both models are major improvements over the previous simple exponential model [3, 5–7] because that model has difficulties estimating the death rate of infected cells when target cells are not depleted after the viral load peak (see Figure 3). The novel models use the observed target cell densities when estimating the parameter values, and by using simulated data we have demonstrated that the new models typically outperform the previous model. Applying the new models to data obtained in humanized mice estimates that the rates at which the virus expands, and at which infected cells die, resemble those measured in humans.
The efficacy of vaccines eliciting cytotoxic immune responses [14–16] could be quantified by our new approach by comparing the estimated death rate of infected cells between normal and vaccinated individuals during acute infection. Additionally, by comparing the death rates estimated during acute infection with those estimated by perturbations with antiviral drugs during the chronic phase [4, 10–13], we might be able to compare the impact of cellular immunity during both phases of an infection. Both examples illustrate that the new models can markedly improve our understanding of viral infection dynamics in vivo.
Our novel approach overcomes the difficulty the previous PWR model had with estimating the death rate of infected cells in situations where target are not severely depleted. Additionally, the choice between the two different models that we propose here can be made on the final densities of the target cells. If the number of target cells reaches a nadir during acute infection, the RQS model seems most appropriate. If the target cells continue to decrease the PWT model should be better.
Thus, whenever one has sufficient time course data from an acute infection, the new models should allow one to estimate the death rate of infected cells, and hence the R_{0}, with reasonable accuracy. Indeed, in animal experiments using rhesus macaques [14–16, 26–28], ferrets [30], and mice [17–20, 31–34], both target cell densities and viral loads have been measured during the acute phase of viral infection. For example, the number of uninfected (and infected) cells in the lungs or respiratory tracts of macaques, ferrets and mice that were experimentally infected with influenza could be measured [30, 31, 34]. Using the CXCR4-tropic SHIV/macaque model, both target cell densities (naïve and memory CD4 T cells) and viral loads from PB have been measured [14, 16, 26–28]. The target cells of simian immunodeficiency virus (SIV), or CCR5-tropic SHIV infection (memory CD4 T cells expressing CCR5), have been measured from gastrointestinal mucosa samples [15]. Thus there are several infection models that can be analyzed with our new models.
Conclusion
In this paper, we developed novel mathematical approaches to estimating parameters from acute viral infection data. We demonstrated that the new models outperform the previous model using simulated data. We quantified the dynamics of acute-phase HIV-1 infections by measuring their time course data in a humanized mouse model. Interestingly we find that the rates at which the virus expands, and at which infected cells die, are similar to those in humans.
Declarations
Acknowledgments
This research is partly supported by the Kyushu University Fund (to H.I.); Grants-in-Aid for Young Scientists B23790500 (to K.S.) and B25800092 (to S.I.) from the Japan Society for the Promotion of Science (JSPS); the Uehara Memorial Foundation (to K.S.); the Shimizu Foundation for Immunological Research Grant (to K.S.); Takeda Science Foundation (to K.S.); Sumitomo Foundation Research Grant (to K.S.); Senshin Medical Research Foundation (to K.S.); JST PRESTO program (to S.I.); the Kyushu University Interdisciplinary Programs in Education and Projects in Research Development (to S.I.); 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 K.S., K.A., and S.I.).
Authors’ Affiliations
References
- Carrat F, Vergu E, Ferguson NM, Lemaitre M, Cauchemez S, Leach S, Valleron AJ: Time lines of infection and disease in human influenza: a review of volunteer challenge studies. Am J Epidemio. 2008, 167 (7): 775-785. 10.1093/aje/kwm375.View ArticleGoogle Scholar
- Peiris JS, Chu CM, Cheng VC, Chan KS, Hung IF, Poon LL, Law KI, Tang BS, Hon TY, Chan CS, Chan KH, Ng JS, Zheng BJ, Ng WL, Lai RW, Guan Y, Yuen KY, HKU/UCH SARS Study Group: Clinical progression and viral load in a community outbreak of coronavirus-associated SARS pneumonia: a prospective study. Lancet. 2003, 361 (9371): 1767-1772. 10.1016/S0140-6736(03)13412-5.View ArticlePubMedGoogle Scholar
- Little SJ, McLean AR, Spina CA, Richman DD, Havlir DV: Viral dynamics of acute HIV-1 infection. J Exp Med. 1999, 190 (6): 841-850. 10.1084/jem.190.6.841.PubMed CentralView ArticlePubMedGoogle Scholar
- Neumann AU, Lam NP, Dahari H, Gretch DR, Wiley TE, Layden TJ, Perelson AS: Hepatitis C viral dynamics in vivo and the antiviral efficacy of interferon-alpha therapy. Science. 1998, 282 (5386): 103-107.View ArticlePubMedGoogle Scholar
- Nowak MA, Lloyd AL, Vasquez GM, Wiltrout TA, Wahl LM, Bischofberger N, Williams J, Kinter A, Fauci AS, Hirsch VM, Lifson JD: Viral dynamics of primary viremia and antiretroviral therapy in simian immunodeficiency virus infection. J Virol. 1997, 71 (10): 7518-7525.PubMed CentralPubMedGoogle Scholar
- Ribeiro RM, Qin L, Chavez LL, Li D, Self SG, Perelson AS: Estimation of the initial viral growth rate and basic reproductive number during acute HIV-1 infection. J Virol. 2010, 84 (12): 6096-6102. 10.1128/JVI.00127-10.PubMed CentralView ArticlePubMedGoogle Scholar
- Nowak MA, May RM: Virus Dynamics. 2000, New York, USA: Oxford University PressGoogle Scholar
- Perelson AS: Modelling viral and immune system dynamics. Nat Rev Immunol. 2002, 2 (1): 28-36. 10.1038/nri700.View ArticlePubMedGoogle Scholar
- Simon V, Ho DD: HIV-1 dynamics in vivo: implications for therapy. Nat Rev Microbiol. 2003, 1 (3): 181-190. 10.1038/nrmicro772.View ArticlePubMedGoogle Scholar
- Perelson AS, Essunger P, Cao Y, Vesanen M, Hurley A, Saksela K, Markowitz M, Ho DD: Decay characteristics of HIV-1-infected compartments during combination therapy. Nature. 1997, 387 (6629): 188-191. 10.1038/387188a0.View ArticlePubMedGoogle Scholar
- Markowitz M, Louie M, Hurley A, Sun E, Di Mascio M, Perelson AS, Ho DD: A novel antiviral intervention results in more accurate assessment of human immunodeficiency virus type 1 replication dynamics and T-cell decay in vivo. J Virol. 2003, 77 (8): 5037-5038. 10.1128/JVI.77.8.5037-5038.2003.PubMed CentralView ArticlePubMedGoogle Scholar
- Dahari H, Guedj J, Perelson AS, Layden TJ: Hepatitis C viral kinetics in the era of direct acting antiviral agents and IL28B. Curr Hepa Rep. 2011, 10 (3): 214-227. 10.1007/s11901-011-0101-7.View ArticleGoogle Scholar
- Iwami S, Koizumi Y, Ikeda H, Kakizoe Y: Quantification of viral infection dynamics in animal experiments. Front Microbiol. 2013, 4: 264-PubMed CentralPubMedGoogle Scholar
- Davenport MP, Zhang L, Shiver JW, Casmiro DR, Ribeiro RM, Perelson AS: Influence of peak viral load on the extent of CD4+ T-cell depletion in simian HIV infection. J Acquir Immune Defic Syndr. 2006, 41 (3): 259-265. 10.1097/01.qai.0000199232.31340.d3.View ArticlePubMedGoogle Scholar
- Wilson DP, Mattapallil JJ, Lay MD, Zhang L, Roederer M, Davenport MP: Estimating the infectivity of CCR5-tropic simian immunodeficiency virus SIV(mac251) in the gut. J Virol. 2007, 81 (15): 8025-8029. 10.1128/JVI.01771-06.PubMed CentralView ArticlePubMedGoogle Scholar
- Petravic J, Ribeiro RM, Casimiro DR, Mattapallil JJ, Roederer M, Shiver JW, Davenport MP: Estimating the impact of vaccination on acute simian-human immunodeficiency virus/simian immunodeficiency virus infections. J Virol. 2008, 82 (23): 11589-11598. 10.1128/JVI.01596-08.PubMed CentralView ArticlePubMedGoogle Scholar
- Nie C, Sato K, Misawa N, Kitayama H, Fujino H, Hiramatsu H, Heike T, Nakahata T, Tanaka Y, Ito M, Koyanagi Y: Selective infection of CD4+ effector memory T lymphocytes leads to preferential depletion of memory T lymphocytes in R5 HIV-1-infected humanized NOD/SCID/IL-2Rgammanull mice. Virology. 2009, 394 (1): 64-72. 10.1016/j.virol.2009.08.011.View ArticlePubMedGoogle Scholar
- Sato K, Izumi T, Misawa N, Kobayashi T, Yamashita Y, Ohmichi M, Ito M, Takaori-Kondo A, Koyanagi Y: Remarkable lethal G-to-A mutations in vif-proficient HIV-1 provirus by individual APOBEC3 proteins in humanized mice. J Virol. 2010, 84 (18): 9546-9556. 10.1128/JVI.00823-10.PubMed CentralView ArticlePubMedGoogle Scholar
- Sato K, Misawa N, Fukuhara M, Iwami S, An DS, Ito M, Koyanagi Y: Vpu augments the initial burst phase of HIV-1 propagation and downregulates BST2 and CD4 in humanized mice. J Virol. 2012, 86 (9): 5000-5013. 10.1128/JVI.07062-11.PubMed CentralView ArticlePubMedGoogle Scholar
- Sato K, Misawa N, Iwami S, Satou Y, Matsuoka M, Ishizaka Y, Ito M, Aihara K, An DS, Koyanagi Y: HIV-1 Vpr accelerates viral replication during acute infection by exploitation of proliferating CD4+ T cells in vivo. PLoS Pathog. 2013, 9 (12): e1003812-10.1371/journal.ppat.1003812.PubMed CentralView ArticlePubMedGoogle Scholar
- Iwasa Y, Michor F, Nowak M: Some basic properties of immune selection. J Theor Biol. 2004, 229 (2): 179-188. 10.1016/j.jtbi.2004.03.013.View ArticlePubMedGoogle Scholar
- Thieme HR: Mathematics in Population Biology. 2003, Princeton, USA: Princeton University PressGoogle Scholar
- Komarova NL, Levy DN, Wodarz D: Effect of synaptic transmission on viral fitness in HIV infection. PLoS One. 2012, 7 (11): e48361-10.1371/journal.pone.0048361.PubMed CentralView ArticlePubMedGoogle Scholar
- Anderson RM: The Kermack-McKendrick epidemic threshold theorem. Bull Math Biol. 1991, 53 (1–2): 3-32.View ArticlePubMedGoogle Scholar
- Diekmann O, Heesterbeek JAP: Mathematical Epidemiology of Infectious Diseases: Model Building, Analysis and Interpretation. 2000, New York, USA: John Wiley & SonsGoogle Scholar
- Igarashi T, Brown CR, Endo Y, Buckler-White A, Plishka R, Bischofberger N, Hirsch V, Martin MA: Macrophage are the principal reservoir and sustain high virus loads in rhesus macaques after the depletion of CD4+ T cells by a highly pathogenic simian immunodeficiency virus/HIV type 1 chimera (SHIV): Implications for HIV-1 infections of humans. Proc Natl Acad Sci U S A. 2001, 98 (2): 658-663. 10.1073/pnas.98.2.658.PubMed CentralView ArticlePubMedGoogle Scholar
- Igarashi T, Endo Y, Nishimura Y, Buckler C, Sadjadpour R, Donau OK, Dumaurier MJ, Plishka RJ, Buckler-White A, Martin MA: Early control of highly pathogenic simian immunodeficiency virus/human immunodeficiency virus chimeric virus infections in rhesus monkeys usually results in long-lasting asymptomatic clinical outcomes. J Virol. 2003, 77 (20): 10829-10840. 10.1128/JVI.77.20.10829-10840.2003.PubMed CentralView ArticlePubMedGoogle Scholar
- Kozyrev IL, Ibuki K, Shimada T, Kuwata T, Takemura T, Hayami M, Miura T: Characterization of less pathogenic infectious molecular clones derived from acute-pathogenic SHIV-89.6p stock virus. Virology. 2001, 282 (1): 6-13. 10.1006/viro.2000.0839.View ArticlePubMedGoogle Scholar
- Theodore TS, Englund G, Buckler-White A, Buckler CE, Martin MA, Peden KW: Construction and characterization of a stable full-length macrophage-tropic HIV type 1 molecular clone that directs the production of high titers of progeny virions. AIDS Res Hum Retroviruses. 1996, 12 (3): 191-194. 10.1089/aid.1996.12.191.View ArticlePubMedGoogle Scholar
- Imai M, Watanabe T, Hatta M, Das SC, Ozawa M, Shinya K, Zhong G, Hanson A, Katsura H, Watanabe S, Li C, Kawakami E, Yamada S, Kiso M, Suzuki Y, Maher EA, Neumann G, Kawaoka Y: Experimental adaptation of an influenza H5 HA confers respiratory droplet transmission to a reassortant H5 HA/H1N1 virus in ferrets. Nature. 2012, 486 (7403): 420-428.PubMed CentralPubMedGoogle Scholar
- Itoh Y, Shinya K, Kiso M, Watanabe T, Sakoda Y, Hatta M, Muramoto Y, Tamura D, Sakai-Tagawa Y, Noda T, Sakabe S, Imai M, Hatta Y, Watanabe S, Li C, Yamada S, Fujii K, Murakami S, Imai H, Kakugawa S, Ito M, Takano R, Iwatsuki-Horimoto K, Shimojima M, Horimoto T, Goto H, Takahashi K, Makino A, Ishigaki H, Nakayama M: In vitro and in vivo characterization of new swine-origin H1N1 influenza viruses. Nature. 2009, 460 (7258): 1021-1025.PubMed CentralPubMedGoogle Scholar
- Mercer DF, Schiller DE, Elliott JF, Douglas DN, Hao C, Rinfret A, Addison WR, Fischer KP, Churchill TA, Lakey JR, Tyrrell DL, Kneteman NM: Hepatitis C virus replication in mice with chimeric human livers. Nat Med. 2001, 7 (8): 927-933. 10.1038/90968.View ArticlePubMedGoogle Scholar
- Miao H, Hollenbaugh JA, Zand MS, Holden-Wiltse J, Mosmann TR, Perelson AS, Wu H, Topham DJ: Quantifying the early immune response and adaptive immune response kinetics in mice infected with influenza A virus. J Virol. 2010, 84 (13): 6687-6698. 10.1128/JVI.00266-10.PubMed CentralView ArticlePubMedGoogle Scholar
- Morita M, Kuba K, Ichikawa A, Nakayama M, Katahira J, Iwamoto R, Watanebe T, Sakabe S, Daidoji T, Nakamura S, Kadowaki A, Ohto T, Nakanishi H, Taguchi R, Nakaya T, Murakami M, Yoneda Y, Arai H, Kawaoka Y, Penninger JM, Arita M, Imai Y: The lipid mediator protectin D1 inhibits influenza virus replication and improves severe influenza. Cell. 2013, 153 (1): 112-125. 10.1016/j.cell.2013.02.027.View ArticlePubMedGoogle Scholar
Copyright
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.