Spatiotemporal Dynamics of Virus Infection Spreading in Tissues

Virus spreading in tissues is determined by virus transport, virus multiplication in host cells and the virus-induced immune response. Cytotoxic T cells remove infected cells with a rate determined by the infection level. The intensity of the immune response has a bell-shaped dependence on the concentration of virus, i.e., it increases at low and decays at high infection levels. A combination of these effects and a time delay in the immune response determine the development of virus infection in tissues like spleen or lymph nodes. The mathematical model described in this work consists of reaction-diffusion equations with a delay. It shows that the different regimes of infection spreading like the establishment of a low level infection, a high level infection or a transition between both are determined by the initial virus load and by the intensity of the immune response. The dynamics of the model solutions include simple and composed waves, and periodic and aperiodic oscillations. The results of analytical and numerical studies of the model provide a systematic basis for a quantitative understanding and interpretation of the determinants of the infection process in target organs and tissues from the image-derived data as well as of the spatiotemporal mechanisms of viral disease pathogenesis, and have direct implications for a biopsy-based medical testing of the chronic infection processes caused by viruses, e.g. HIV, HCV and HBV.

[1]  E. Wherry,et al.  Molecular and cellular insights into T cell exhaustion , 2015, Nature Reviews Immunology.

[2]  R. Zinkernagel,et al.  Regulation of the Immune Response by Antigen , 2001, Science.

[3]  A. Ducrot,et al.  Asymptotic behaviour of travelling waves for the delayed Fisher–KPP equation , 2014 .

[4]  Gui-Quan Sun,et al.  Pattern formation of an epidemic model with diffusion , 2012, Nonlinear Dynamics.

[5]  Todd M. Allen,et al.  Innate immune reconstitution with suppression of HIV-1. , 2016, JCI insight.

[6]  Peter Turchin,et al.  Complex Population Dynamics , 2003 .

[7]  A. Bertoletti,et al.  PD‐1/PD‐L1 pathway and T‐cell exhaustion in chronic hepatitis virus infection , 2010, Journal of viral hepatitis.

[8]  Vitaly Volpert,et al.  Elliptic Partial Differential Equations: Volume 2: Reaction-Diffusion Equations , 2014 .

[9]  Anna Medyukhina,et al.  Image‐based systems biology of infection , 2015, Cytometry. Part A : the journal of the International Society for Analytical Cytology.

[10]  Xiao-Qiang Zhao,et al.  Bistable Traveling Waves for Monotone Semiflows with Applications , 2011, 1102.4556.

[11]  R. D. Driver Introduction to Delay Differential Equations , 1977 .

[12]  T. M. Khlebodarova,et al.  Mathematical model of the Tat-Rev regulation of HIV-1 replication in an activated cell predicts the existence of oscillatory dynamics in the synthesis of viral components , 2014, BMC Genomics.

[13]  R. Zinkernagel,et al.  On the role of viruses in the evolution of immune responses. , 1985, British medical bulletin.

[14]  S. Petrovskii,et al.  Delay driven spatiotemporal chaos in single species population dynamics models. , 2016, Theoretical population biology.

[15]  G. Webb,et al.  A mathematical model of cell-to-cell spread of HIV-1 that includes a time delay , 2003, Journal of mathematical biology.

[16]  Zhen Jin,et al.  Influence of isolation degree of spatial patterns on persistence of populations , 2016 .

[17]  A. Bertoletti,et al.  Immune Response in Hepatitis B Virus Infection. , 2015, Cold Spring Harbor perspectives in medicine.

[18]  M. Macal,et al.  Innate and Adaptive Immune Regulation During Chronic Viral Infections. , 2015, Annual review of virology.

[19]  G. Fogel,et al.  HIV DNA Is Frequently Present within Pathologic Tissues Evaluated at Autopsy from Combined Antiretroviral Therapy-Treated Patients with Undetectable Viral Loads , 2016, Journal of Virology.

[20]  R. Germain,et al.  The human condition: an immunological perspective , 2011, Nature Immunology.

[21]  Rolf M. Zinkernagel,et al.  Virus persistence in acutely infected immunocompetent mice by exhaustion of antiviral cytotoxic effector T cells , 1993, Nature.

[22]  R. Germain,et al.  Spatiotemporal basis of innate and adaptive immunity in secondary lymphoid tissue. , 2014, Annual review of cell and developmental biology.

[23]  R. Zinkernagel,et al.  On Immunity Against Infections and Vaccines: Credo 2004 , 2004, Scandinavian journal of immunology.

[24]  Zhen Jin,et al.  Effects of time delay and space on herbivore dynamics: linking inducible defenses of plants to herbivore outbreak , 2015, Scientific Reports.

[25]  J. M. Murray,et al.  Turing Patterns from Dynamics of Early HIV Infection , 2012, Bulletin of mathematical biology.

[26]  Wan-Tong Li,et al.  Existence and stability of traveling wave fronts in reaction advection diffusion equations with nonlocal delay , 2007 .

[27]  A. Oxenius,et al.  Antigen amount dictates CD8+ T‐cell exhaustion during chronic viral infection irrespective of the type of antigen presenting cell , 2012, European journal of immunology.

[28]  V. Volpert,et al.  Reaction-diffusion waves in biology. , 2009, Physics of life reviews.

[29]  R. Germain,et al.  The future of immunoimaging — Deeper, bigger, more precise, and definitively more colorful , 2013, European journal of immunology.

[30]  Martin Meier-Schellersheim,et al.  Pathogenesis of HIV infection: what the virus spares is as important as what it destroys , 2006, Nature Medicine.

[31]  Ronald N Germain,et al.  Visualization and dynamic analysis of host-pathogen interactions. , 2014, Current opinion in immunology.

[32]  S. Trofimchuk,et al.  Slowly oscillating wavefronts of the KPP-Fisher delayed equation , 2012, 1206.0484.

[33]  Stephanie Forrest,et al.  A spatial model of the efficiency of T cell search in the influenza-infected lung. , 2016, Journal of theoretical biology.

[34]  Alan S. Perelson,et al.  Viral kinetic modeling: state of the art , 2014, Journal of Pharmacokinetics and Pharmacodynamics.

[35]  Ian Williams,et al.  Hybrid Spreading Mechanisms and T Cell Activation Shape the Dynamics of HIV-1 Infection , 2015, PLoS Comput. Biol..

[36]  C. Mims The Pathogenesis of the Acute Exanthems , 1996 .

[37]  R. Cheynier,et al.  The infiltration kinetics of simian immunodeficiency virus-specific T cells drawn to sites of high antigenic stimulation determines local in vivo viral escape , 2001, Proceedings of the National Academy of Sciences of the United States of America.

[38]  Jianhong Wu,et al.  Nonlocality of Reaction-Diffusion Equations Induced by Delay: Biological Modeling and Nonlinear Dynamics , 2004 .

[39]  Vitaly Volpert,et al.  Elliptic Partial Differential Equations , 2011 .

[40]  R. Dunia,et al.  Mathematical modeling of viral infection dynamics in spherical organs , 2013, Journal of mathematical biology.

[41]  A. Perelson,et al.  Quantification of Hepatitis C Virus Cell-to-Cell Spread Using a Stochastic Modeling Approach , 2015, Journal of Virology.

[42]  A. Bertoletti,et al.  Adaptive immunity in HBV infection. , 2016, Journal of hepatology.

[43]  Vitaly Volpert,et al.  Traveling Wave Solutions of Parabolic Systems , 1994 .

[44]  Li Li,et al.  Patch invasion in a spatial epidemic model , 2015, Appl. Math. Comput..

[45]  D. Richman,et al.  Spatiotemporal dynamics of HIV propagation. , 2002, Journal of theoretical biology.

[46]  Alan S. Perelson,et al.  Mathematical Analysis of HIV-1 Dynamics in Vivo , 1999, SIAM Rev..

[47]  Henri Berestycki,et al.  The non-local Fisher–KPP equation: travelling waves and steady states , 2009 .

[48]  Gui-Quan Sun,et al.  Mathematical modeling of population dynamics with Allee effect , 2016, Nonlinear Dynamics.

[49]  Xiao-Qiang Zhao,et al.  Global Asymptotic Stability of Traveling Waves in Delayed Reaction-Diffusion Equations , 2000, SIAM J. Math. Anal..

[50]  Jianhong Wu,et al.  Existence, Uniqueness and Asymptotic Stability of Traveling Wavefronts in A Non-Local Delayed Diffusion Equation , 2007 .

[51]  Wenzhang Huang,et al.  Periodic Traveling Waves for Diffusion Equations with Time Delayed and Non-local Responding Reaction , 2007 .

[52]  O. Pybus,et al.  HIV Maintains an Evolving and Dispersed Population in Multiple Tissues during Suppressive Combined Antiretroviral Therapy in Individuals with Cancer , 2016, Journal of Virology.

[53]  Sergei Trofimchuk,et al.  Slowly oscillating wave solutions of a single species reaction–diffusion equation with delay , 2008 .

[54]  C. Beauchemin Probing the effects of the well-mixed assumption on viral infection dynamics. , 2005, Journal of theoretical biology.

[55]  Sotiris A. Prokopiou,et al.  Multiscale Modeling of the Early CD8 T-Cell Immune Response in Lymph Nodes: An Integrative Study , 2014, Comput..

[56]  A. Haase,et al.  Large number of rebounding/founder HIV variants emerge from multifocal infection in lymphatic tissues after treatment interruption , 2015, Proceedings of the National Academy of Sciences.

[57]  Sebastian Bonhoeffer,et al.  Spatial models of virus-immune dynamics. , 2005, Journal of theoretical biology.

[58]  Elsje Pienaar,et al.  A computational tool integrating host immunity with antibiotic dynamics to study tuberculosis treatment. , 2015, Journal of theoretical biology.

[59]  B. Sua,et al.  Mathematical modelling of immune response in tissues , 2009 .

[60]  Hal L. Smith,et al.  Abstract functional-differential equations and reaction-diffusion systems , 1990 .

[61]  Jaap Goudsmit,et al.  Ongoing HIV dissemination during HAART , 1999, Nature Medicine.

[62]  Z. Grossman,et al.  Pathogen-sensing, regulatory T cells, and responsiveness-tuning collectively regulate foreign- and self-antigen mediated T-cell responses. , 2013, Cold Spring Harbor symposia on quantitative biology.

[63]  Massimo Bernaschi,et al.  Modeling lymphocyte homing and encounters in lymph nodes , 2009, BMC Bioinformatics.

[64]  U. Schwarz,et al.  Multiscale modeling of virus replication and spread , 2016, FEBS letters.

[65]  Klaus W. Schaaf Asymptotic behavior and traveling wave solutions for parabolic functional-differential equations , 1987 .

[66]  Z. Grossman,et al.  Dynamic tuning of lymphocytes: physiological basis, mechanisms, and function. , 2015, Annual review of immunology.

[67]  S. Petrovskii,et al.  Long-term transients and complex dynamics of a stage-structured population with time delay and the Allee effect. , 2016, Journal of theoretical biology.

[68]  Hal L. Smith,et al.  An introduction to delay differential equations with applications to the life sciences / Hal Smith , 2010 .

[69]  G. Bocharov Understanding Complex Regulatory Systems: Integrating Molecular Biology and Systems Analysis , 2005, Transfusion Medicine and Hemotherapy.

[70]  J. Rawlings,et al.  Image-Guided Modeling of Virus Growth and Spread , 2008, Bulletin of mathematical biology.

[71]  S. Iwami,et al.  Dynamics of HIV infection in lymphoid tissue network , 2016, Journal of mathematical biology.

[72]  S. Trofimchuk,et al.  Traveling waves for a model of the Belousov–Zhabotinsky reaction , 2011, 1103.0176.

[73]  A. Sigal,et al.  HIV Cell-to-Cell Spread Results in Earlier Onset of Viral Gene Expression by Multiple Infections per Cell , 2016, PLoS pathogens.

[74]  Jianhong Wu Theory and Applications of Partial Functional Differential Equations , 1996 .

[75]  Zhen Jin,et al.  Pattern transitions in spatial epidemics: Mechanisms and emergent properties , 2016, Physics of Life Reviews.

[76]  Sang-Jun Ha,et al.  Visualizing Antigen-Specific and Infected Cells in Situ Predicts Outcomes in Early Viral Infection , 2009, Science.

[77]  Roy M. Anderson,et al.  Underwhelming the Immune Response: Effect of Slow Virus Growth on CD8+-T-Lymphocyte Responses , 2004, Journal of Virology.

[78]  Juan Zhang,et al.  Pattern formation of a spatial predator-prey system , 2012, Appl. Math. Comput..

[79]  G. Kochs,et al.  Temporal and Spatial Resolution of Type I and III Interferon Responses In Vivo , 2010, Journal of Virology.

[80]  F. Fenner,et al.  The pathogenesis of the acute exanthems; an interpretation based on experimental investigations with mousepox; infectious ectromelia of mice. , 1948, Lancet.

[81]  Sergei I. Trofimchuk,et al.  Global continuation of monotone wavefronts , 2012, J. Lond. Math. Soc..