Matching Items (49)
Description
Using a simple $SI$ infection model, I uncover the
overall dynamics of the system and how they depend on the incidence
function. I consider both an epidemic and endemic perspective of the
model, but in both cases, three classes of incidence
functions are identified.
In the epidemic form,
power incidences, where the infective portion $I^p$ has $p\in(0,1)$,
cause unconditional host extinction,
homogeneous incidences have host extinction for certain parameter constellations and
host survival for others, and upper density-dependent incidences
never cause host extinction. The case of non-extinction in upper
density-dependent
incidences extends to the case where a latent period is included.
Using data from experiments with rhanavirus and salamanders,
maximum likelihood estimates are applied to the data.
With these estimates,
I generate the corrected Akaike information criteria, which
reward a low likelihood and punish the use of more parameters.
This generates the Akaike weight, which is used to fit
parameters to the data, and determine which incidence functions
fit the data the best.
From an endemic perspective, I observe
that power incidences cause initial condition dependent host extinction for
some parameter constellations and global stability for others,
homogeneous incidences have host extinction for certain parameter constellations and
host survival for others, and upper density-dependent incidences
never cause host extinction.
The dynamics when the incidence function is homogeneous are deeply explored.
I expand the endemic considerations in the homogeneous case
by adding a predator into the model.
Using persistence theory, I show the conditions for the persistence of each of the
predator, prey, and parasite species. Potential dynamics of the system include parasite mediated
persistence of the predator, survival of the ecosystem at high initial predator levels and
ecosystem collapse at low initial predator levels, persistence of all three species, and much more.
overall dynamics of the system and how they depend on the incidence
function. I consider both an epidemic and endemic perspective of the
model, but in both cases, three classes of incidence
functions are identified.
In the epidemic form,
power incidences, where the infective portion $I^p$ has $p\in(0,1)$,
cause unconditional host extinction,
homogeneous incidences have host extinction for certain parameter constellations and
host survival for others, and upper density-dependent incidences
never cause host extinction. The case of non-extinction in upper
density-dependent
incidences extends to the case where a latent period is included.
Using data from experiments with rhanavirus and salamanders,
maximum likelihood estimates are applied to the data.
With these estimates,
I generate the corrected Akaike information criteria, which
reward a low likelihood and punish the use of more parameters.
This generates the Akaike weight, which is used to fit
parameters to the data, and determine which incidence functions
fit the data the best.
From an endemic perspective, I observe
that power incidences cause initial condition dependent host extinction for
some parameter constellations and global stability for others,
homogeneous incidences have host extinction for certain parameter constellations and
host survival for others, and upper density-dependent incidences
never cause host extinction.
The dynamics when the incidence function is homogeneous are deeply explored.
I expand the endemic considerations in the homogeneous case
by adding a predator into the model.
Using persistence theory, I show the conditions for the persistence of each of the
predator, prey, and parasite species. Potential dynamics of the system include parasite mediated
persistence of the predator, survival of the ecosystem at high initial predator levels and
ecosystem collapse at low initial predator levels, persistence of all three species, and much more.
ContributorsFarrell, Alexander E. (Author) / Thieme, Horst R (Thesis advisor) / Smith, Hal (Committee member) / Kuang, Yang (Committee member) / Tang, Wenbo (Committee member) / Collins, James (Committee member) / Arizona State University (Publisher)
Created2017
Description
The increased number of novel pathogens that potentially threaten the human population has motivated the development of mathematical and computational modeling approaches for forecasting epidemic impact and understanding key environmental characteristics that influence the spread of diseases. Yet, in the case that substantial uncertainty surrounds the transmission process during a rapidly developing infectious disease outbreak, complex mechanistic models may be too difficult to be calibrated quick enough for policy makers to make informed decisions. Simple phenomenological models that rely on a small number of parameters can provide an initial platform for assessing the epidemic trajectory, estimating the reproduction number and quantifying the disease burden from the early epidemic phase.
Chapter 1 provides background information and motivation for infectious disease forecasting and outlines the rest of the thesis.
In chapter 2, logistic patch models are used to assess and forecast the 2013-2015 West Africa Zaire ebolavirus epidemic. In particular, this chapter is concerned with comparing and contrasting the effects that spatial heterogeneity has on the forecasting performance of the cumulative infected case counts reported during the epidemic.
In chapter 3, two simple phenomenological models inspired from population biology are used to assess the Research and Policy for Infectious Disease Dynamics (RAPIDD) Ebola Challenge; a simulated epidemic that generated 4 infectious disease scenarios. Because of the nature of the synthetically generated data, model predictions are compared to exact epidemiological quantities used in the simulation.
In chapter 4, these models are applied to the 1904 Plague epidemic that occurred in Bombay. This chapter provides evidence that these simple models may be applicable to infectious diseases no matter the disease transmission mechanism.
Chapter 5, uses the patch models from chapter 2 to explore how migration in the 1904 Plague epidemic changes the final epidemic size.
The final chapter is an interdisciplinary project concerning within-host dynamics of cereal yellow dwarf virus-RPV, a plant pathogen from a virus group that infects over 150 grass species. Motivated by environmental nutrient enrichment due to anthropological activities, mathematical models are employed to investigate the relevance of resource competition to pathogen and host dynamics.
Chapter 1 provides background information and motivation for infectious disease forecasting and outlines the rest of the thesis.
In chapter 2, logistic patch models are used to assess and forecast the 2013-2015 West Africa Zaire ebolavirus epidemic. In particular, this chapter is concerned with comparing and contrasting the effects that spatial heterogeneity has on the forecasting performance of the cumulative infected case counts reported during the epidemic.
In chapter 3, two simple phenomenological models inspired from population biology are used to assess the Research and Policy for Infectious Disease Dynamics (RAPIDD) Ebola Challenge; a simulated epidemic that generated 4 infectious disease scenarios. Because of the nature of the synthetically generated data, model predictions are compared to exact epidemiological quantities used in the simulation.
In chapter 4, these models are applied to the 1904 Plague epidemic that occurred in Bombay. This chapter provides evidence that these simple models may be applicable to infectious diseases no matter the disease transmission mechanism.
Chapter 5, uses the patch models from chapter 2 to explore how migration in the 1904 Plague epidemic changes the final epidemic size.
The final chapter is an interdisciplinary project concerning within-host dynamics of cereal yellow dwarf virus-RPV, a plant pathogen from a virus group that infects over 150 grass species. Motivated by environmental nutrient enrichment due to anthropological activities, mathematical models are employed to investigate the relevance of resource competition to pathogen and host dynamics.
ContributorsPell, Bruce (Author) / Kuang, Yang (Thesis advisor) / Chowell-Puente, Gerardo (Committee member) / Nagy, John (Committee member) / Kostelich, Eric (Committee member) / Gardner, Carl (Committee member) / Arizona State University (Publisher)
Created2016
Description
Glioblastoma multiforme (GBMs) is the most prevalent brain tumor type and causes approximately 40% of all non-metastic primary tumors in adult patients [1]. GBMs are malignant, grade-4 brain tumors, the most aggressive classication as established by the World Health Organization and are marked by their low survival rate; the median survival time is only twelve months from initial diagnosis: Patients who live more than three years are considered long-term survivors [2]. GBMs are highly invasive and their diffusive growth pattern makes it impossible to remove the tumors by surgery alone [3]. The purpose of this paper is to use individual patient data to parameterize a model of GBMs that allows for data on tumor growth and development to be captured on a clinically relevant time scale. Such an endeavor is the rst step to a clinically applicable predictions of GBMs. Previous research has yielded models that adequately represent the development of GBMs, but they have not attempted to follow specic patient cases through the entire tumor process. Using the model utilized by Kostelich et al. [4], I will attempt to redress this deciency. In doing so, I will improve upon a family of models that can be used to approximate the time of development and/or structure evolution in GBMs. The eventual goal is to incorporate Magnetic Resonance Imaging (MRI) data into a parameterized model of GBMs in such a way that it can be used clinically to predict tumor growth and behavior. Furthermore, I hope to come to a denitive conclusion as to the accuracy of the Koteslich et al. model throughout the development of GBMs tumors.
ContributorsManning, Miles (Author) / Kostelich, Eric (Thesis director) / Kuang, Yang (Committee member) / Preul, Mark (Committee member) / Barrett, The Honors College (Contributor) / College of Liberal Arts and Sciences (Contributor)
Created2012-12
Description
Over time, tumor treatment resistance inadvertently develops when androgen de-privation therapy (ADT) is applied to metastasized prostate cancer (PCa). To combat tumor resistance, while reducing the harsh side effects of hormone therapy, the clinician may opt to cyclically alternates the patient’s treatment on and off. This method,known as intermittent ADT, is an alternative to continuous ADT that improves the patient’s quality of life while testosterone levels recover between cycles. In this paper,we explore the response of intermittent ADT to metastasized prostate cancer by employing a previously clinical data validated mathematical model to new clinical data from patients undergoing Abiraterone therapy. This cell quota model, a system of ordinary differential equations constructed using Droop’s nutrient limiting theory, assumes the tumor comprises of castration-sensitive (CS) and castration-resistant (CR)cancer sub-populations. The two sub-populations rely on varying levels of intracellular androgen for growth, death and transformation. Due to the complexity of the model,we carry out sensitivity analyses to study the effect of certain parameters on their outputs, and to increase the identifiability of each patient’s unique parameter set. The model’s forecasting results show consistent accuracy for patients with sufficient data,which means the model could give useful information in practice, especially to decide whether an additional round of treatment would be effective.
ContributorsBennett, Justin Klark (Author) / Kuang, Yang (Thesis director) / Kostelich, Eric (Committee member) / Phan, Tin (Committee member) / School of Mathematical and Statistical Sciences (Contributor, Contributor) / Barrett, The Honors College (Contributor)
Created2021-05
Description
Five immunocompetent C57BL/6-cBrd/cBrd/Cr (albino C57BL/6) mice were injected with GL261-luc2 cells, a cell line sharing characteristics of human glioblastoma multiforme (GBM). The mice were imaged using magnetic resonance (MR) at five separate time points to characterize growth and development of the tumor. After 25 days, the final tumor volumes of the mice varied from 12 mm3 to 62 mm3, even though mice were inoculated from the same tumor cell line under carefully controlled conditions. We generated hypotheses to explore large variances in final tumor size and tested them with our simple reaction-diffusion model in both a 3-dimensional (3D) finite difference method and a 2-dimensional (2D) level set method. The parameters obtained from a best-fit procedure, designed to yield simulated tumors as close as possible to the observed ones, vary by an order of magnitude between the three mice analyzed in detail. These differences may reflect morphological and biological variability in tumor growth, as well as errors in the mathematical model, perhaps from an oversimplification of the tumor dynamics or nonidentifiability of parameters. Our results generate parameters that match other experimental in vitro and in vivo measurements. Additionally, we calculate wave speed, which matches with other rat and human measurements.
ContributorsRutter, Erica (Author) / Stepien, Tracy (Author) / Anderies, Barrett (Author) / Plasencia, Jonathan (Author) / Woolf, Eric C. (Author) / Scheck, Adrienne C. (Author) / Turner, Gregory H. (Author) / Liu, Qingwei (Author) / Frakes, David (Author) / Kodibagkar, Vikram (Author) / Kuang, Yang (Author) / Preul, Mark C. (Author) / Kostelich, Eric (Author) / College of Liberal Arts and Sciences (Contributor)
Created2017-05-31
Description
Background:
Data assimilation refers to methods for updating the state vector (initial condition) of a complex spatiotemporal model (such as a numerical weather model) by combining new observations with one or more prior forecasts. We consider the potential feasibility of this approach for making short-term (60-day) forecasts of the growth and spread of a malignant brain cancer (glioblastoma multiforme) in individual patient cases, where the observations are synthetic magnetic resonance images of a hypothetical tumor.
Results:
We apply a modern state estimation algorithm (the Local Ensemble Transform Kalman Filter), previously developed for numerical weather prediction, to two different mathematical models of glioblastoma, taking into account likely errors in model parameters and measurement uncertainties in magnetic resonance imaging. The filter can accurately shadow the growth of a representative synthetic tumor for 360 days (six 60-day forecast/update cycles) in the presence of a moderate degree of systematic model error and measurement noise.
Conclusions:
The mathematical methodology described here may prove useful for other modeling efforts in biology and oncology. An accurate forecast system for glioblastoma may prove useful in clinical settings for treatment planning and patient counseling.
ContributorsKostelich, Eric (Author) / Kuang, Yang (Author) / McDaniel, Joshua (Author) / Moore, Nina Z. (Author) / Martirosyan, Nikolay L. (Author) / Preul, Mark C. (Author) / College of Liberal Arts and Sciences (Contributor)
Created2011-12-21
Description
A description of numerical and analytical work pertaining to models that describe the growth and progression of glioblastoma multiforme (GBM), an aggressive form of primary brain cancer. Two reaction-diffusion models are used: the Fisher-Kolmogorov-Petrovsky-Piskunov equation and a 2-population model that divides the tumor into actively proliferating and quiescent (or necrotic) cells. The numerical portion of this work (chapter 2) focuses on simulating GBM expansion in patients undergoing treatment for recurrence of tumor following initial surgery. The models are simulated on 3-dimensional brain geometries derived from magnetic resonance imaging (MRI) scans provided by the Barrow Neurological Institute. The study consists of 17 clinical time intervals across 10 patients that have been followed in detail, each of whom shows significant progression of tumor over a period of 1 to 3 months on sequential follow up scans. A Taguchi sampling design is implemented to estimate the variability of the predicted tumors to using 144 different choices of model parameters. In 9 cases, model parameters can be identified such that the simulated tumor contains at least 40 percent of the volume of the observed tumor. In the analytical portion of the paper (chapters 3 and 4), a positively invariant region for our 2-population model is identified. Then, a rigorous derivation of the critical patch size associated with the model is performed. The critical patch (KISS) size is the minimum habitat size needed for a population to survive in a region. Habitats larger than the critical patch size allow a population to persist, while smaller habitats lead to extinction. The critical patch size of the 2-population model is consistent with that of the Fisher-Kolmogorov-Petrovsky-Piskunov equation, one of the first reaction-diffusion models proposed for GBM. The critical patch size may indicate that GBM tumors have a minimum size depending on the location in the brain. A theoretical relationship between the size of a GBM tumor at steady-state and its maximum cell density is also derived, which has potential applications for patient-specific parameter estimation based on magnetic resonance imaging data.
ContributorsHarris, Duane C. (Author) / Kuang, Yang (Thesis advisor) / Kostelich, Eric J. (Thesis advisor) / Preul, Mark C. (Committee member) / Crook, Sharon (Committee member) / Gardner, Carl (Committee member) / Arizona State University (Publisher)
Created2023
Description
Over the past 20 years, the fields of synthetic biology and synthetic biosystems engineering have grown into mature disciplines, leading to significant breakthroughs in cancer research, diagnostics, cell-based medicines, biochemical production, etc. Application of mathematical modelling to biological and biochemical systems have not only given great insight into how these systems function, but also have lent enough predictive power to aid in the forward-engineering of synthetic constructs. However, progress has been impeded by several modes of context-dependence unique to biological and biochemical systems that are not seen in traditional engineering disciplines, resulting in the need for lengthy design-build-test cycles before functional prototypes are generated.In this work, two of these universal modes of context dependence – resource competition and growth feedback –their effects on synthetic gene circuits and potential control mechanisms, are studied and characterized. Results demonstrate that a novel competitive control architecture can be utilized to mitigate the effects of winner-take-all resource competition (a form of context dependence where distinct gene modules influence each other by competing over a shared pool of transcriptional/translational resources) in synthetic gene circuits and restore circuits to their intended function. Application of the fluctuation-dissipation theorem and rigorous stochastic simulations demonstrate that realistic resource constraints present in cells at the transcriptional and translational levels influence noise in gene circuits in a nonmonotonic fashion, either increasing or decreasing noise depending on the transcriptional/translational capacity.
Growth feedback on the other hand links circuit function to cellular growth rate via increased protein dilution rate during exponential growth phase. This in turn can result in the collapse of bistable gene circuits as the accelerated dilution rate forces switches in a high stable state to fall to a low stable state. Mathematical modelling and experimental data demonstrate that application of repressive links can insulate sensitive parts of gene circuits against growth-fluctuations and can in turn increase the robustness of multistable circuits in growth contexts.
The results presented in this work aid in the accumulation of understanding of biological and biochemical context dependence, and corresponding control strategies and design principles engineers can utilize to mitigate these effects.
ContributorsStone, Austin (Author) / Tian, Xiao-jun (Thesis advisor) / Wang, Xiao (Committee member) / Smith, Barbara (Committee member) / Kuang, Yang (Committee member) / Cheng, Albert (Committee member) / Arizona State University (Publisher)
Created2023
Description
There is a need in the ecology literature to have a discussion about the fundamental theories from which population dynamics arises. Ad hoc model development is not uncommon in the field often as a result of a need to publish rapidly and frequently. Ecologists and statisticians like Robert J. Steidl and Kenneth P Burnham have called for a more deliberative approach they call "hard thinking". For example, the phenomena of population growth can be captured by almost any sigmoid function. The question of which sigmoid function best explains a data set cannot be answered meaningfully by statistical regression since that can only speak to the validity of the shape. There is a need to revisit enzyme kinetics and ecological stoichiometry to properly justify basal model selection in ecology. This dissertation derives several common population growth models from a generalized equation. The mechanistic validity of these models in different contexts is explored through a kinetic lens. The behavioral kinetic framework is then put to the test by examining a set of biologically plausible growth models against the 1968-1995 elk population count data for northern Yellowstone. Using only this count data, the novel Monod-Holling growth model was able to accurately predict minimum viable population and life expectancy despite both being exogenous to the model and data set. Lastly, the elk/wolf data from Yellowstone was used to compare the validity of the Rosenzweig-MacArthur and Arditi-Ginzburg models. They both were derived from a more general model which included both predator and prey mediated steps. The Arditi-Ginzburg model was able to fit the training data better, but only the Rosenzweig-MacArthur model matched the validation data. Accounting for animal sexual behavior allowed for the creation of the Monod-Holling model which is just as simple as the logistic differential equation but provides greater insights for conservation purposes. Explicitly acknowledging the ethology of wolf predation helps explain the differences in predictive performances by the best fit Rosenzweig-MacArthur and Arditi-Ginzburg models. The behavioral kinetic framework has proven to be a useful tool, and it has the ability to provide even further insights going forward.
ContributorsPringle, Jack Andrew McCracken (Author) / Anderies, John M (Thesis advisor) / Kuang, Yang (Committee member) / Milner, Fabio (Committee member) / Arizona State University (Publisher)
Created2022
Description
Ecology has been an actively studied topic recently, along with the rapid development of human microbiota-based technology. Scientists have made remarkable progress using bioinformatics tools to identify species and analyze composition. However, a thorough understanding of interspecies interactions of microbial ecosystems is still lacking, which has been a significant obstacle in the further development of related technologies. In this work, a genetic circuit design principle with synthetic biology approaches is developed to form two-strain microbial consortia with different inter-strain interactions. The microbial systems are well-defined and inducible. Co-culture experiment results show that our microbial consortia behave consistently with previous ecological knowledge and thus serves as excellent model systems to simulate ecosystems with similar interactions. Colony patterns also emerge when co-culturing multiple species on solid media. With the engineered microbial consortia, image-processing based methods were developed to quantify the shape of co-culture colonies and distinguish microbial consortia with different interactions. Factors that affect the population ratios were identified through induction and variations in the inoculation process. Further time-lapse experiments revealed the basic rules of colony growth, composition variation, patterning, and how spatial factors impact the co-culture colony.
ContributorsChen, Xingwen (Author) / Wang, Xiao (Thesis advisor) / Kuang, Yang (Committee member) / Tian, Xiaojun (Committee member) / Brafman, David (Committee member) / Plaisier, Christopher (Committee member) / Arizona State University (Publisher)
Created2022