Matching Items (106)
Filtering by
- Language: English
Description
The role of climate change, as measured in terms of changes in the climatology of geophysical variables (such as temperature and rainfall), on the global distribution and burden of vector-borne diseases (VBDs) remains a subject of considerable debate. This dissertation attempts to contribute to this debate via the use of mathematical (compartmental) modeling and statistical data analysis. In particular, the objective is to find suitable values and/or ranges of the climate variables considered (typically temperature and rainfall) for maximum vector abundance and consequently, maximum transmission intensity of the disease(s) they cause.
Motivated by the fact that understanding the dynamics of disease vector is crucial to understanding the transmission and control of the VBDs they cause, a novel weather-driven deterministic model for the population biology of the mosquito is formulated and rigorously analyzed. Numerical simulations, using relevant weather and entomological data for Anopheles mosquito (the vector for malaria), show that maximum mosquito abundance occurs when temperature and rainfall values lie in the range [20-25]C and [105-115] mm, respectively.
The Anopheles mosquito ecology model is extended to incorporate human dynamics. The resulting weather-driven malaria transmission model, which includes many of the key aspects of malaria (such as disease transmission by asymptomatically-infectious humans, and enhanced malaria immunity due to repeated exposure), was rigorously analyzed. The model which also incorporates the effect of diurnal temperature range (DTR) on malaria transmission dynamics shows that increasing DTR shifts the peak temperature value for malaria transmission from 29C (when DTR is 0C) to about 25C (when DTR is 15C).
Finally, the malaria model is adapted and used to study the transmission dynamics of chikungunya, dengue and Zika, three diseases co-circulating in the Americas caused by the same vector (Aedes aegypti). The resulting model, which is fitted using data from Mexico, is used to assess a few hypotheses (such as those associated with the possible impact the newly-released dengue vaccine will have on Zika) and the impact of variability in climate variables on the dynamics of the three diseases. Suitable temperature and rainfall ranges for the maximum transmission intensity of the three diseases are obtained.
Motivated by the fact that understanding the dynamics of disease vector is crucial to understanding the transmission and control of the VBDs they cause, a novel weather-driven deterministic model for the population biology of the mosquito is formulated and rigorously analyzed. Numerical simulations, using relevant weather and entomological data for Anopheles mosquito (the vector for malaria), show that maximum mosquito abundance occurs when temperature and rainfall values lie in the range [20-25]C and [105-115] mm, respectively.
The Anopheles mosquito ecology model is extended to incorporate human dynamics. The resulting weather-driven malaria transmission model, which includes many of the key aspects of malaria (such as disease transmission by asymptomatically-infectious humans, and enhanced malaria immunity due to repeated exposure), was rigorously analyzed. The model which also incorporates the effect of diurnal temperature range (DTR) on malaria transmission dynamics shows that increasing DTR shifts the peak temperature value for malaria transmission from 29C (when DTR is 0C) to about 25C (when DTR is 15C).
Finally, the malaria model is adapted and used to study the transmission dynamics of chikungunya, dengue and Zika, three diseases co-circulating in the Americas caused by the same vector (Aedes aegypti). The resulting model, which is fitted using data from Mexico, is used to assess a few hypotheses (such as those associated with the possible impact the newly-released dengue vaccine will have on Zika) and the impact of variability in climate variables on the dynamics of the three diseases. Suitable temperature and rainfall ranges for the maximum transmission intensity of the three diseases are obtained.
ContributorsOkuneye, Kamaldeen O (Author) / Gumel, Abba B (Thesis advisor) / Kuang, Yang (Committee member) / Smith, Hal (Committee member) / Thieme, Horst (Committee member) / Nagy, John (Committee member) / Arizona State University (Publisher)
Created2018
Description
Synthetic manipulation of chromatin dynamics has applications for medicine, agriculture, and biotechnology. However, progress in this area requires the identification of design rules for engineering chromatin systems. In this thesis, I discuss research that has elucidated the intrinsic properties of histone binding proteins (HBP), and apply this knowledge to engineer novel chromatin binding effectors. Results from the experiments described herein demonstrate that the histone binding domain from chromobox protein homolog 8 (CBX8) is portable and can be customized to alter its endogenous function. First, I developed an assay to identify engineered fusion proteins that bind histone post translational modifications (PTMs) in vitro and regulate genes near the same histone PTMs in living cells. This assay will be useful for assaying the function of synthetic histone PTM-binding actuators and probes. Next, I investigated the activity of a novel, dual histone PTM binding domain regulator called Pc2TF. I characterized Pc2TF in vitro and in cells and show it has enhanced binding and transcriptional activation compared to a single binding domain fusion called Polycomb Transcription Factor (PcTF). These results indicate that valency can be used to tune the activity of synthetic histone-binding transcriptional regulators. Then, I report the delivery of PcTF fused to a cell penetrating peptide (CPP) TAT, called CP-PcTF. I treated 2D U-2 OS bone cancer cells with CP-PcTF, followed by RNA sequencing to identify genes regulated by CP-PcTF. I also showed that 3D spheroids treated with CP-PcTF show delayed growth. This preliminary work demonstrated that an epigenetic effector fused to a CPP can enable entry and regulation of genes in U-2 OS cells through DNA independent interactions. Finally, I described and validated a new screening method that combines the versatility of in vitro transcription and translation (IVTT) expressed protein coupled with the histone tail microarrays. Using Pc2TF as an example, I demonstrated that this assay is capable of determining binding and specificity of a synthetic HBP. I conclude by outlining future work toward engineering HBPs using techniques such as directed evolution and rational design. In conclusion, this work outlines a foundation to engineer and deliver synthetic chromatin effectors.
ContributorsTekel, Stefan (Author) / Haynes, Karmella (Thesis advisor) / Mills, Jeremy (Committee member) / Caplan, Michael (Committee member) / Brafman, David (Committee member) / Arizona State University (Publisher)
Created2019
Description
Rabies is an infectious viral disease. It is usually fatal if a victim reaches the rabid stage, which starts after the appearance of disease symptoms. The disease virus attacks the central nervous system, and then it migrates from peripheral nerves to the spinal cord and brain. At the time when the rabies virus reaches the brain, the incubation period is over and the symptoms of clinical disease appear on the victim. From the brain, the virus travels via nerves to the salivary glands and saliva.
A mathematical model is developed for the spread of rabies in a spatially distributed fox population to model the spread of the rabies epizootic through middle Europe that occurred in the second half of the 20th century. The model considers both territorial and wandering rabid foxes and includes a latent period for the infection. Since the model assumes these two kinds of rabid foxes, it is a system of both partial differential and integral equations (with integration
over space and, occasionally, also over time). To study the spreading speeds of the rabies epidemic, the model is reduced to a scalar Volterra-Hammerstein integral equation, and space-time Laplace transform of the integral equation is used to derive implicit formulas for the spreading speed. The spreading speeds are discussed and implicit formulas are given for latent periods of fixed length, exponentially distributed length, Gamma distributed length, and log-normally distributed length. A number of analytic and numerical results are shown pertaining to the spreading speeds.
Further, a numerical algorithm is described for the simulation
of the spread of rabies in a spatially distributed fox population on a bounded domain with Dirichlet boundary conditions. I propose the following methods for the numerical approximation of solutions. The partial differential and integral equations are discretized in the space variable by central differences of second order and by
the composite trapezoidal rule. Next, the ordinary or delay differential equations that are obtained this way are discretized in time by explicit
continuous Runge-Kutta methods of fourth order for ordinary and delay differential systems. My particular interest
is in how the partition of rabid foxes into
territorial and diffusing rabid foxes influences
the spreading speed, a question that can be answered by purely analytic means only for small basic reproduction numbers. I will restrict the numerical analysis
to latent periods of fixed length and to exponentially
distributed latent periods.
The results of the numerical calculations
are compared for latent periods
of fixed and exponentially distributed length
and for various proportions of territorial
and wandering rabid foxes.
The speeds of spread observed in the
simulations are compared
to spreading speeds obtained by numerically solving the analytic formulas
and to observed speeds of epizootic frontlines
in the European rabies outbreak 1940 to 1980.
A mathematical model is developed for the spread of rabies in a spatially distributed fox population to model the spread of the rabies epizootic through middle Europe that occurred in the second half of the 20th century. The model considers both territorial and wandering rabid foxes and includes a latent period for the infection. Since the model assumes these two kinds of rabid foxes, it is a system of both partial differential and integral equations (with integration
over space and, occasionally, also over time). To study the spreading speeds of the rabies epidemic, the model is reduced to a scalar Volterra-Hammerstein integral equation, and space-time Laplace transform of the integral equation is used to derive implicit formulas for the spreading speed. The spreading speeds are discussed and implicit formulas are given for latent periods of fixed length, exponentially distributed length, Gamma distributed length, and log-normally distributed length. A number of analytic and numerical results are shown pertaining to the spreading speeds.
Further, a numerical algorithm is described for the simulation
of the spread of rabies in a spatially distributed fox population on a bounded domain with Dirichlet boundary conditions. I propose the following methods for the numerical approximation of solutions. The partial differential and integral equations are discretized in the space variable by central differences of second order and by
the composite trapezoidal rule. Next, the ordinary or delay differential equations that are obtained this way are discretized in time by explicit
continuous Runge-Kutta methods of fourth order for ordinary and delay differential systems. My particular interest
is in how the partition of rabid foxes into
territorial and diffusing rabid foxes influences
the spreading speed, a question that can be answered by purely analytic means only for small basic reproduction numbers. I will restrict the numerical analysis
to latent periods of fixed length and to exponentially
distributed latent periods.
The results of the numerical calculations
are compared for latent periods
of fixed and exponentially distributed length
and for various proportions of territorial
and wandering rabid foxes.
The speeds of spread observed in the
simulations are compared
to spreading speeds obtained by numerically solving the analytic formulas
and to observed speeds of epizootic frontlines
in the European rabies outbreak 1940 to 1980.
ContributorsAlanazi, Khalaf Matar (Author) / Thieme, Horst R. (Thesis advisor) / Jackiewicz, Zdzislaw (Committee member) / Baer, Steven (Committee member) / Gardner, Carl (Committee member) / Kuang, Yang (Committee member) / Smith, Hal (Committee member) / Arizona State University (Publisher)
Created2018
Description
Neurodegenerative diseases such as Alzheimer’s disease, Parkinson’s disease, or amyotrophic lateral sclerosis are defined by the loss of several types of neurons and glial cells within the central nervous system (CNS). Combatting these diseases requires a robust population of relevant cell types that can be employed in cell therapies, drug screening, or patient specific disease modeling. Human induced pluripotent stem cells (hiPSC)-derived neural progenitor cells (hNPCs) have the ability to self-renew indefinitely and differentiate into the various neuronal and glial cell types of the CNS. In order to realize the potential of hNPCs, it is necessary to develop a xeno-free scalable platform for effective expansion and differentiation. Previous work in the Brafman lab led to the engineering of a chemically defined substrate—vitronectin derived peptide (VDP), which allows for the long-term expansion and differentiation of hNPCs. In this work, we use this substrate as the basis for a microcarrier (MC)-based suspension culture system. Several independently derived hNPC lines were cultured on MCs for multiple passages as well as efficiently differentiated to neurons. Finally, this MC-based system was used in conjunction with a low shear rotating wall vessel (RWV) bioreactor for the integrated, large-scale expansion and neuronal differentiation of hNPCs. Finally, VDP was shown to support the differentiation of hNPCs into functional astrocytes. Overall, this fully defined and scalable biomanufacturing system will facilitate the generation of hNPCs and their derivatives in quantities necessary for basic and translational applications.
ContributorsMorgan, Daylin (Author) / Brafman, David (Thesis advisor) / Stabenfeldt, Sarah (Committee member) / Wang, Xiao (Committee member) / Arizona State University (Publisher)
Created2018
Description
Calcium imaging is a well-established, non-invasive or minimally technique designed to study the electrical signaling neurons. Calcium regulates the release of gliotransmitters in astrocytes. Analyzing astrocytic calcium transients can provide significant insights into mechanisms such as neuroplasticity and neural signal modulation.
In the past decade, numerous methods have been developed to analyze in-vivo calcium imaging data that involves complex techniques such as overlapping signals segregation and motion artifact correction. The hypothesis used to detect calcium signal is the spatiotemporal sparsity of calcium signal, and these methods are unable to identify the passive cells that are not actively firing during the time frame in the video. Statistics regarding the percentage of cells in each frame of view can be critical for the analysis of calcium imaging data for human induced pluripotent stem cells derived neurons and astrocytes.
The objective of this research is to develop a simple and efficient semi-automated pipeline for analysis of in-vitro calcium imaging data. The region of interest (ROI) based image segmentation is used to extract the data regarding intensity fluctuation caused by calcium concentration changes in each cell. It is achieved by using two approaches: basic image segmentation approach and a machine learning approach. The intensity data is evaluated using a custom-made MATLAB that generates statistical information and graphical representation of the number of spiking cells in each field of view, the number of spikes per cell and spike height.
In the past decade, numerous methods have been developed to analyze in-vivo calcium imaging data that involves complex techniques such as overlapping signals segregation and motion artifact correction. The hypothesis used to detect calcium signal is the spatiotemporal sparsity of calcium signal, and these methods are unable to identify the passive cells that are not actively firing during the time frame in the video. Statistics regarding the percentage of cells in each frame of view can be critical for the analysis of calcium imaging data for human induced pluripotent stem cells derived neurons and astrocytes.
The objective of this research is to develop a simple and efficient semi-automated pipeline for analysis of in-vitro calcium imaging data. The region of interest (ROI) based image segmentation is used to extract the data regarding intensity fluctuation caused by calcium concentration changes in each cell. It is achieved by using two approaches: basic image segmentation approach and a machine learning approach. The intensity data is evaluated using a custom-made MATLAB that generates statistical information and graphical representation of the number of spiking cells in each field of view, the number of spikes per cell and spike height.
ContributorsBhandarkar, Siddhi Umesh (Author) / Brafman, David (Thesis advisor) / Stabenfeldt, Sarah (Committee member) / Tian, Xiaojun (Committee member) / Arizona State University (Publisher)
Created2019
Description
Ideas from coding theory are employed to theoretically demonstrate the engineering of mutation-tolerant genes, genes that can sustain up to some arbitrarily chosen number of mutations and still express the originally intended protein. Attention is restricted to tolerating substitution mutations. Future advances in genomic engineering will make possible the ability to synthesize entire genomes from scratch. This presents an opportunity to embed desirable capabilities like mutation-tolerance, which will be useful in preventing cell deaths in organisms intended for research or industrial applications in highly mutagenic environments. In the extreme case, mutation-tolerant genes (mutols) can make organisms resistant to retroviral infections.
An algebraic representation of the nucleotide bases is developed. This algebraic representation makes it possible to convert nucleotide sequences into algebraic sequences, apply mathematical ideas and convert results back into nucleotide terms. Using the algebra developed, a mapping is found from the naturally-occurring codons to an alternative set of codons which makes genes constructed from them mutation-tolerant, provided no more than one substitution mutation occurs per codon. The ideas discussed naturally extend to finding codons that can tolerate t arbitrarily chosen number of mutations per codon. Finally, random substitution events are simulated in both a wild-type green fluorescent protein (GFP) gene and its mutol variant and the amino acid sequence expressed from each post-mutation is compared with the amino acid sequence pre-mutation.
This work assumes the existence of synthetic protein-assembling entities that function like tRNAs but can read k nucleotides at a time, with k greater than or equal to 5. The realization of this assumption is presented as a challenge to the research community.
An algebraic representation of the nucleotide bases is developed. This algebraic representation makes it possible to convert nucleotide sequences into algebraic sequences, apply mathematical ideas and convert results back into nucleotide terms. Using the algebra developed, a mapping is found from the naturally-occurring codons to an alternative set of codons which makes genes constructed from them mutation-tolerant, provided no more than one substitution mutation occurs per codon. The ideas discussed naturally extend to finding codons that can tolerate t arbitrarily chosen number of mutations per codon. Finally, random substitution events are simulated in both a wild-type green fluorescent protein (GFP) gene and its mutol variant and the amino acid sequence expressed from each post-mutation is compared with the amino acid sequence pre-mutation.
This work assumes the existence of synthetic protein-assembling entities that function like tRNAs but can read k nucleotides at a time, with k greater than or equal to 5. The realization of this assumption is presented as a challenge to the research community.
ContributorsAmpofo, Prince Kwame (Author) / Tian, Xiaojun (Thesis advisor) / Kiani, Samira (Committee member) / Kuang, Yang (Committee member) / Arizona State University (Publisher)
Created2019
Description
Cancer is a major health problem in the world today and is expected to become an even larger one in the future. Although cancer therapy has improved for many cancers in the last several decades, there is much room for further improvement. Mathematical modeling has the advantage of being able to test many theoretical therapies without having to perform clinical trials and experiments. Mathematical oncology will continue to be an important tool in the future regarding cancer therapies and management.
This dissertation is structured as a growing tumor. Chapters 2 and 3 consider spheroid models. These models are adept at describing 'early-time' tumors, before the tumor needs to co-opt blood vessels to continue sustained growth. I consider two partial differential equation (PDE) models for spheroid growth of glioblastoma. I compare these models to in vitro experimental data for glioblastoma tumor cell lines and other proposed models. Further, I investigate the conditions under which traveling wave solutions exist and confirm numerically.
As a tumor grows, it can no longer be approximated by a spheroid, and it becomes necessary to use in vivo data and more sophisticated modeling to model the growth and diffusion. In Chapter 4, I explore experimental data and computational models for describing growth and diffusion of glioblastoma in murine brains. I discuss not only how the data was obtained, but how the 3D brain geometry is created from Magnetic Resonance (MR) images. A 3D finite-difference code is used to model tumor growth using a basic reaction-diffusion equation. I formulate and test hypotheses as to why there are large differences between the final tumor sizes between the mice.
Once a tumor has reached a detectable size, it is diagnosed, and treatment begins. Chapter 5 considers modeling the treatment of prostate cancer. I consider a joint model with hormonal therapy as well as immunotherapy. I consider a timing study to determine whether changing the vaccine timing has any effect on the outcome of the patient. In addition, I perform basic analysis on the six-dimensional ordinary differential equation (ODE). I also consider the limiting case, and perform a full global analysis.
This dissertation is structured as a growing tumor. Chapters 2 and 3 consider spheroid models. These models are adept at describing 'early-time' tumors, before the tumor needs to co-opt blood vessels to continue sustained growth. I consider two partial differential equation (PDE) models for spheroid growth of glioblastoma. I compare these models to in vitro experimental data for glioblastoma tumor cell lines and other proposed models. Further, I investigate the conditions under which traveling wave solutions exist and confirm numerically.
As a tumor grows, it can no longer be approximated by a spheroid, and it becomes necessary to use in vivo data and more sophisticated modeling to model the growth and diffusion. In Chapter 4, I explore experimental data and computational models for describing growth and diffusion of glioblastoma in murine brains. I discuss not only how the data was obtained, but how the 3D brain geometry is created from Magnetic Resonance (MR) images. A 3D finite-difference code is used to model tumor growth using a basic reaction-diffusion equation. I formulate and test hypotheses as to why there are large differences between the final tumor sizes between the mice.
Once a tumor has reached a detectable size, it is diagnosed, and treatment begins. Chapter 5 considers modeling the treatment of prostate cancer. I consider a joint model with hormonal therapy as well as immunotherapy. I consider a timing study to determine whether changing the vaccine timing has any effect on the outcome of the patient. In addition, I perform basic analysis on the six-dimensional ordinary differential equation (ODE). I also consider the limiting case, and perform a full global analysis.
ContributorsRutter, Erica Marie (Author) / Kuang, Yang (Thesis advisor) / Kostelich, Eric J (Thesis advisor) / Frakes, David (Committee member) / Gardner, Carl (Committee member) / Jackiewicz, Zdzislaw (Committee member) / Arizona State University (Publisher)
Created2016
Description
Several debilitating neurological disorders, such as Alzheimer's disease, stroke, and spinal cord injury, are characterized by the damage or loss of neuronal cell types in the central nervous system (CNS). Human neural progenitor cells (hNPCs) derived from human pluripotent stem cells (hPSCs) can proliferate extensively and differentiate into the various neuronal subtypes and supporting cells that comprise the CNS. As such, hNPCs have tremendous potential for disease modeling, drug screening, and regenerative medicine applications. However, the use hNPCs for the study and treatment of neurological diseases requires the development of defined, robust, and scalable methods for their expansion and neuronal differentiation. To that end a rational design process was used to develop a vitronectin-derived peptide (VDP)-based substrate to support the growth and neuronal differentiation of hNPCs in conventional two-dimensional (2-D) culture and large-scale microcarrier (MC)-based suspension culture. Compared to hNPCs cultured on ECMP-based substrates, hNPCs grown on VDP-coated surfaces displayed similar morphologies, growth rates, and high expression levels of hNPC multipotency markers. Furthermore, VDP surfaces supported the directed differentiation of hNPCs to neurons at similar levels to cells differentiated on ECMP substrates. Here it has been demonstrated that VDP is a robust growth and differentiation matrix, as demonstrated by its ability to support the expansions and neuronal differentiation of hNPCs derived from three hESC (H9, HUES9, and HSF4) and one hiPSC (RiPSC) cell lines. Finally, it has been shown that VDP allows for the expansion or neuronal differentiation of hNPCs to quantities (>1010) necessary for drug screening or regenerative medicine purposes. In the future, the use of VDP as a defined culture substrate will significantly advance the clinical application of hNPCs and their derivatives as it will enable the large-scale expansion and neuronal differentiation of hNPCs in quantities necessary for disease modeling, drug screening, and regenerative medicine applications.
ContributorsVarun, Divya (Author) / Brafman, David (Thesis advisor) / Nikkhah, Mehdi (Committee member) / Stabenfeldt, Sarah (Committee member) / Arizona State University (Publisher)
Created2016
Description
In recent decades, marine ecologists have conducted extensive field work and experiments to understand the interactions between bacteria and bacteriophage (phage) in marine environments. This dissertation provides a detailed rigorous framework for gaining deeper insight into these interactions. Specific features of the dissertation include the design of a new deterministic Lotka-Volterra model with n + 1 bacteria, n
+ 1 phage, with explicit nutrient, where the jth phage strain infects the first j bacterial strains, a perfectly nested infection network (NIN). This system is subject to trade-off conditions on the life-history traits of both bacteria and phage given in an earlier study Jover et al. (2013). Sufficient conditions are provided to show that a bacteria-phage community of arbitrary size with NIN can arise through the succession of permanent subcommunities, by the successive addition of one new population. Using uniform persistence theory, this entire community is shown to be permanent (uniformly persistent), meaning that all populations ultimately survive.
It is shown that a modified version of the original NIN Lotka-Volterra model with implicit nutrient considered by Jover et al. (2013) is permanent. A new one-to-one infection network (OIN) is also considered where each bacterium is infected by only one phage, and that phage infects only that bacterium. This model does not use the trade-offs on phage infection range, and bacterium resistance to phage. The OIN model is shown to be permanent, and using Lyapunov function theory, coupled with LaSalle’s Invariance Principle, the unique coexistence equilibrium associated with the NIN is globally asymptotically stable provided that the inter- and intra-specific bacterial competition coefficients are equal across all bacteria.
Finally, the OIN model is extended to a “Kill the Winner” (KtW) Lotka-Volterra model
of marine communities consisting of bacteria, phage, and zooplankton. The zooplankton
acts as a super bacteriophage, which infects all bacteria. This model is shown to be permanent.
+ 1 phage, with explicit nutrient, where the jth phage strain infects the first j bacterial strains, a perfectly nested infection network (NIN). This system is subject to trade-off conditions on the life-history traits of both bacteria and phage given in an earlier study Jover et al. (2013). Sufficient conditions are provided to show that a bacteria-phage community of arbitrary size with NIN can arise through the succession of permanent subcommunities, by the successive addition of one new population. Using uniform persistence theory, this entire community is shown to be permanent (uniformly persistent), meaning that all populations ultimately survive.
It is shown that a modified version of the original NIN Lotka-Volterra model with implicit nutrient considered by Jover et al. (2013) is permanent. A new one-to-one infection network (OIN) is also considered where each bacterium is infected by only one phage, and that phage infects only that bacterium. This model does not use the trade-offs on phage infection range, and bacterium resistance to phage. The OIN model is shown to be permanent, and using Lyapunov function theory, coupled with LaSalle’s Invariance Principle, the unique coexistence equilibrium associated with the NIN is globally asymptotically stable provided that the inter- and intra-specific bacterial competition coefficients are equal across all bacteria.
Finally, the OIN model is extended to a “Kill the Winner” (KtW) Lotka-Volterra model
of marine communities consisting of bacteria, phage, and zooplankton. The zooplankton
acts as a super bacteriophage, which infects all bacteria. This model is shown to be permanent.
ContributorsKorytowski, Daniel (Author) / Smith, Hal (Thesis advisor) / Gumel, Abba (Committee member) / Kuang, Yang (Committee member) / Gardner, Carl (Committee member) / Thieme, Horst (Committee member) / Arizona State University (Publisher)
Created2016
Description
There has been important progress in understanding ecological dynamics through the development of the theory of ecological stoichiometry. This fast growing theory provides new constraints and mechanisms that can be formulated into mathematical models. Stoichiometric models incorporate the effects of both food quantity and food quality into a single framework that produce rich dynamics. While the effects of nutrient deficiency on consumer growth are well understood, recent discoveries in ecological stoichiometry suggest that consumer dynamics are not only affected by insufficient food nutrient content (low phosphorus (P): carbon (C) ratio) but also by excess food nutrient content (high P:C). This phenomenon, known as the stoichiometric knife edge, in which animal growth is reduced not only by food with low P content but also by food with high P content, needs to be incorporated into mathematical models. Here we present Lotka-Volterra type models to investigate the growth response of Daphnia to algae of varying P:C ratios. Using a nonsmooth system of two ordinary differential equations (ODEs), we formulate the first model to incorporate the phenomenon of the stoichiometric knife edge. We then extend this stoichiometric model by mechanistically deriving and tracking free P in the environment. This resulting full knife edge model is a nonsmooth system of three ODEs. Bifurcation analysis and numerical simulations of the full model, that explicitly tracks phosphorus, leads to quantitatively different predictions than previous models that neglect to track free nutrients. The full model shows that the grazer population is sensitive to excess nutrient concentrations as a dynamical free nutrient pool induces extreme grazer population density changes. These modeling efforts provide insight on the effects of excess nutrient content on grazer dynamics and deepen our understanding of the effects of stoichiometry on the mechanisms governing population dynamics and the interactions between trophic levels.
ContributorsPeace, Angela (Author) / Kuang, Yang (Thesis advisor) / Elser, James J (Committee member) / Baer, Steven (Committee member) / Tang, Wenbo (Committee member) / Kang, Yun (Committee member) / Arizona State University (Publisher)
Created2014