A mathematical model for within-host Toxoplasma gondii invasion dynamics

Dynamics of toxoplasmosis in the cat's population with an exposed stage and a time delay

We propose a new mathematical model to investigate the effect of the introduction of an exposed stage for the cats who become infected with the T. gondii parasite, but that are not still able to produce oocysts in the environment. The model considers a time delay in order to represent the duration of the exposed stage. Besides the cat population the model also includes the oocysts related to the T. gondii in the environment. The model includes the cats since they are the only definitive host and the oocysts, since they are relevant to the dynamics of toxoplasmosis. The model considers lifelong immunity for the recovered cats and vaccinated cats. In addition, the model considers that cats can get infected through an effective contact with the oocysts in the environment. We find conditions such that the toxoplasmosis disease becomes extinct. We analyze the consequences of considering the exposed stage and the time delay on the stability of the equilibrium points. We numerically solve the constructed model and corroborated the theoretical results.

Emerging disease dynamics in a model coupling within-host and between-host systems

Epidemiological models and immunological models have been studied largely independently. However, the two processes (between- and within-host interactions) occur jointly and models that couple the two processes may generate new biological insights. Particularly, the threshold conditions for disease control may be dramatically different when compared with those generated from the epidemiological or immunological models separately. An example is considered in this paper for an environmentally driven infectious disease such as Toxoplasma gondii. The model explicitly couples the within-host and between-host dynamics. The within-host sub-system is linked to a contaminated environment E via an additional term g(E) to account for the increase in the parasite load V within a host due to the continuous ingestion of parasites from the contaminated environment. The parasite load V can also affect the rate of environmental contamination, which directly contributes to the infection rate of hosts for the between-host sub-system. When the two sub-systems are considered in isolation, the dynamics are standard and simple. That is, either the infection-free equilibrium is stable or a unique positive equilibrium is stable depending on the relevant reproduction number being less or greater than 1. However, when the two sub-systems are explicitly coupled, the full system exhibits more complex dynamics including backward bifurcations; that is, multiple positive equilibria exist with one of which being stable even if the reproduction number is less than 1. The biological implications of such bifurcations are illustrated using an example concerning the spread and control of toxoplasmosis.

Evidence for finely-regulated asynchronous growth of Toxoplasma gondii cysts based on data-driven model selection

Toxoplasma gondii establishes a chronic infection by forming cysts preferentially in the brain. This chronic infection is one of the most common parasitic infections in humans and can be reactivated to develop life-threatening toxoplasmic encephalitis in immunocompromised patients. Host-pathogen interactions during the chronic infection include growth of the cysts and their removal by both natural rupture and elimination by the immune system. Analyzing these interactions is important for understanding the pathogenesis of this common infection. We developed a differential equation framework of cyst growth and employed Akaike Information Criteria (AIC) to determine the growth and removal functions that best describe the distribution of cyst sizes measured from the brains of chronically infected mice. The AIC strongly support models in which T. gondii cysts grow at a constant rate such that the per capita growth rate of the parasite is inversely proportional to the number of parasites within a cyst, suggesting finely-regulated asynchronous replication of the parasites. Our analyses were also able to reject the models where cyst removal rate increases linearly or quadratically in association with increase in cyst size. The modeling and analysis framework may provide a useful tool for understanding the pathogenesis of infections with other cyst producing parasites.

A large portion of people worldwide are chronically infected with T. gondii. Chronic infection with this parasite is characterized by formation of tissue cysts. Bradyzoites slowly replicate within cysts during the chronic stage of infection leading to a corresponding increase in cyst size. Cysts occasionally rupture and release bradyzoites that invade nearby host cells and convert into tachyzoites which can quickly proliferate. Tissue cysts can also be targeted by immune T cells and phagocytes for removal. We developed a differential equation model to investigate the cumulative effects of unknown growth and removal functions on the cyst-size distribution. We then used the AIC to select models that best fit experimental cyst size distribution data obtained from the brains of chronically infected mice. The results suggest that the within-cyst growth of bradyzoites is finely-regulated asynchronous such that the per capita growth rate is inversely proportional to the number of bradyzoites. While it may be surprising that the bradyzoites replicate in a way analogous to a factory producing a product, there may be factors such as nutrient availability, resource allocation, immune response, and other stress factors that may limit replication in cysts.