Matching Items (21)
152291-Thumbnail Image.png
Description
Rabies disease remains enzootic among raccoons, skunks, foxes and bats in the United States. It is of primary concern for public-health agencies to control spatial spread of rabies in wildlife and its potential spillover infection of domestic animals and humans. Rabies is invariably fatal in wildlife if untreated, with a

Rabies disease remains enzootic among raccoons, skunks, foxes and bats in the United States. It is of primary concern for public-health agencies to control spatial spread of rabies in wildlife and its potential spillover infection of domestic animals and humans. Rabies is invariably fatal in wildlife if untreated, with a non-negligible incubation period. Understanding how this latency affects spatial spread of rabies in wildlife is the concern of chapter 2 and 3. Chapter 1 deals with the background of mathematical models for rabies and lists main objectives. In chapter 2, a reaction-diffusion susceptible-exposed-infected (SEI) model and a delayed diffusive susceptible-infected (SI) model are constructed to describe the same epidemic process -- rabies spread in foxes. For the delayed diffusive model a non-local infection term with delay is resulted from modeling the dispersal during incubation stage. Comparison is made regarding minimum traveling wave speeds of the two models, which are verified using numerical experiments. In chapter 3, starting with two Kermack and McKendrick's models where infectivity, death rate and diffusion rate of infected individuals can depend on the age of infection, the asymptotic speed of spread $c^\ast$ for the cumulated force of infection can be analyzed. For the special case of fixed incubation period, the asymptotic speed of spread is governed by the same integral equation for both models. Although explicit solutions for $c^\ast$ are difficult to obtain, assuming that diffusion coefficient of incubating animals is small, $c^\ast$ can be estimated in terms of model parameter values. Chapter 4 considers the implementation of realistic landscape in simulation of rabies spread in skunks and bats in northeast Texas. The Finite Element Method (FEM) is adopted because the irregular shapes of realistic landscape naturally lead to unstructured grids in the spatial domain. This implementation leads to a more accurate description of skunk rabies cases distributions.
ContributorsLiu, Hao (Author) / Kuang, Yang (Thesis advisor) / Jackiewicz, Zdzislaw (Committee member) / Lanchier, Nicolas (Committee member) / Smith, Hal (Committee member) / Thieme, Horst (Committee member) / Arizona State University (Publisher)
Created2013
153500-Thumbnail Image.png
Description
Parkinson's disease (PD) is a neurodegenerative disorder that produces a characteristic set of neuromotor deficits that sometimes includes reduced amplitude and velocity of movement. Several studies have shown that people with PD improved their motor performance when presented with external cues. Other work has demonstrated that high velocity

Parkinson's disease (PD) is a neurodegenerative disorder that produces a characteristic set of neuromotor deficits that sometimes includes reduced amplitude and velocity of movement. Several studies have shown that people with PD improved their motor performance when presented with external cues. Other work has demonstrated that high velocity and large amplitude exercises can increase the amplitude and velocity of movement in simple carryover tasks in the upper and lower extremities. Although the cause for these effects is not known, improvements due to cueing suggest that part of the neuromotor deficit in PD is in the integration of sensory feedback to produce motor commands. Previous studies have documented some somatosensory deficits, but only limited information is available regarding the nature and magnitude of sensorimotor deficits in the shoulder of people with PD. The goals of this research were to characterize the sensorimotor impairment in the shoulder joint of people with PD and to investigate the use of visual feedback and large amplitude/high velocity exercises to target PD-related motor deficits. Two systems were designed and developed to use visual feedback to assess the ability of participants to accurately adjust limb placement or limb movement velocity and to encourage improvements in performance of these tasks. Each system was tested on participants with PD, age-matched control subjects and young control subjects to characterize and compare limb placement and velocity control capabilities. Results demonstrated that participants with PD were less accurate at placing their limbs than age-matched or young control subjects, but that their performance improved over the course of the test session such that by the end, the participants with PD performed as well as controls. For the limb velocity feedback task, participants with PD and age-matched control subjects were less accurate than young control subjects, but at the end of the session, participants with PD and age-matched control subjects were as accurate as the young control subjects. This study demonstrates that people with PD were able to improve their movement patterns based on visual feedback of performance and suggests that this feedback paradigm may be useful in exercise programs for people with PD.
ContributorsSmith, Catherine (Author) / Abbas, James J (Thesis advisor) / Ingalls, Todd (Thesis advisor) / Krishnamurthi, Narayanan (Committee member) / Buneo, Christopher (Committee member) / Rikakis, Thanassis (Committee member) / Arizona State University (Publisher)
Created2015
153468-Thumbnail Image.png
Description
The phycologist, M. R. Droop, studied vitamin B12 limitation in the flagellate Monochrysis lutheri and concluded that its specific growth rate depended on the concentration of the vitamin within the cell; i.e. the cell quota of the vitamin B12. The Droop model provides a mathematical expression to link growth rate

The phycologist, M. R. Droop, studied vitamin B12 limitation in the flagellate Monochrysis lutheri and concluded that its specific growth rate depended on the concentration of the vitamin within the cell; i.e. the cell quota of the vitamin B12. The Droop model provides a mathematical expression to link growth rate to the intracellular concentration of a limiting nutrient. Although the Droop model has been an important modeling tool in ecology, it has only recently been applied to study cancer biology. Cancer cells live in an ecological setting, interacting and competing with normal and other cancerous cells for nutrients and space, and evolving and adapting to their environment. Here, the Droop equation is used to model three cancers.

First, prostate cancer is modeled, where androgen is considered the limiting nutrient since most tumors depend on androgen for proliferation and survival. The model's accuracy for predicting the biomarker for patients on intermittent androgen deprivation therapy is tested by comparing the simulation results to clinical data as well as to an existing simpler model. The results suggest that a simpler model may be more beneficial for a predictive use, although further research is needed in this field prior to implementing mathematical models as a predictive method in a clinical setting.

Next, two chronic myeloid leukemia models are compared that consider Imatinib treatment, a drug that inhibits the constitutively active tyrosine kinase BCR-ABL. Both models describe the competition of leukemic and normal cells, however the first model also describes intracellular dynamics by considering BCR-ABL as the limiting nutrient. Using clinical data, the differences in estimated parameters between the models and the capacity for each model to predict drug resistance are analyzed.

Last, a simple model is presented that considers ovarian tumor growth and tumor induced angiogenesis, subject to on and off anti-angiogenesis treatment. In this environment, the cell quota represents the intracellular concentration of necessary nutrients provided through blood supply. Mathematical analysis of the model is presented and model simulation results are compared to pre-clinical data. This simple model is able to fit both on- and off-treatment data using the same biologically relevant parameters.
ContributorsEverett, Rebecca Anne (Author) / Kuang, Yang (Thesis advisor) / Nagy, John (Committee member) / Milner, Fabio (Committee member) / Crook, Sharon (Committee member) / Jackiewicz, Zdzislaw (Committee member) / Arizona State University (Publisher)
Created2015
153054-Thumbnail Image.png
Description
During attempted fixation, the eyes are not still but continue to produce so called "fixational eye movements", which include microsaccades, drift, and tremor. Microsaccades are thought to help prevent and restore vision loss during fixation, and to correct fixation errors, but how they contribute to these functions remains a matter

During attempted fixation, the eyes are not still but continue to produce so called "fixational eye movements", which include microsaccades, drift, and tremor. Microsaccades are thought to help prevent and restore vision loss during fixation, and to correct fixation errors, but how they contribute to these functions remains a matter of debate. This dissertation presents the results of four experiments conducted to address current controversies concerning the role of microsaccades in visibility and oculomotor control.

The first two experiments set out to correlate microsaccade production with the visibility of foveal and peripheral targets of varied spatial frequencies, during attempted fixation. The results indicate that microsaccades restore the visibility of both peripheral targets and targets presented entirely within the fovea, as a function of their spatial frequency characteristics.

The last two experiments set out to determine the role of microsaccades and drifts on the correction of gaze-position errors due to blinks in human and non-human primates, and to characterize microsaccades forming square-wave jerks (SWJs) in non-human primates. The results showed that microsaccades, but not drifts, correct gaze-position errors due to blinks, and that SWJ production and dynamic properties are equivalent in human and non-human primates.

These combined findings suggest that microsaccades, like saccades, serve multiple and non-exclusive functional roles in vision and oculomotor control, as opposed to having a single specialized function.
ContributorsCostela, Francisco M (Author) / Crook, Sharon M (Committee member) / Martinez-Conde, Susana (Committee member) / Macknik, Stephen L. (Committee member) / Baer, Stephen (Committee member) / McCamy, Michael B (Committee member) / Arizona State University (Publisher)
Created2014
149953-Thumbnail Image.png
Description
The theme for this work is the development of fast numerical algorithms for sparse optimization as well as their applications in medical imaging and source localization using sensor array processing. Due to the recently proposed theory of Compressive Sensing (CS), the $\ell_1$ minimization problem attracts more attention for its ability

The theme for this work is the development of fast numerical algorithms for sparse optimization as well as their applications in medical imaging and source localization using sensor array processing. Due to the recently proposed theory of Compressive Sensing (CS), the $\ell_1$ minimization problem attracts more attention for its ability to exploit sparsity. Traditional interior point methods encounter difficulties in computation for solving the CS applications. In the first part of this work, a fast algorithm based on the augmented Lagrangian method for solving the large-scale TV-$\ell_1$ regularized inverse problem is proposed. Specifically, by taking advantage of the separable structure, the original problem can be approximated via the sum of a series of simple functions with closed form solutions. A preconditioner for solving the block Toeplitz with Toeplitz block (BTTB) linear system is proposed to accelerate the computation. An in-depth discussion on the rate of convergence and the optimal parameter selection criteria is given. Numerical experiments are used to test the performance and the robustness of the proposed algorithm to a wide range of parameter values. Applications of the algorithm in magnetic resonance (MR) imaging and a comparison with other existing methods are included. The second part of this work is the application of the TV-$\ell_1$ model in source localization using sensor arrays. The array output is reformulated into a sparse waveform via an over-complete basis and study the $\ell_p$-norm properties in detecting the sparsity. An algorithm is proposed for minimizing a non-convex problem. According to the results of numerical experiments, the proposed algorithm with the aid of the $\ell_p$-norm can resolve closely distributed sources with higher accuracy than other existing methods.
ContributorsShen, Wei (Author) / Mittlemann, Hans D (Thesis advisor) / Renaut, Rosemary A. (Committee member) / Jackiewicz, Zdzislaw (Committee member) / Gelb, Anne (Committee member) / Ringhofer, Christian (Committee member) / Arizona State University (Publisher)
Created2011
150711-Thumbnail Image.png
Description
In vertebrate outer retina, changes in the membrane potential of horizontal cells affect the calcium influx and glutamate release of cone photoreceptors via a negative feedback. This feedback has a number of important physiological consequences. One is called background-induced flicker enhancement (BIFE) in which the onset of dim background enhances

In vertebrate outer retina, changes in the membrane potential of horizontal cells affect the calcium influx and glutamate release of cone photoreceptors via a negative feedback. This feedback has a number of important physiological consequences. One is called background-induced flicker enhancement (BIFE) in which the onset of dim background enhances the center flicker response of horizontal cells. The underlying mechanism for the feedback is still unclear but competing hypotheses have been proposed. One is the GABA hypothesis, which states that the feedback is mediated by gamma-aminobutyric acid (GABA), an inhibitory neurotransmitter released from horizontal cells. Another is the ephaptic hypothesis, which contends that the feedback is non-GABAergic and is achieved through the modulation of electrical potential in the intersynaptic cleft between cones and horizontal cells. In this study, a continuum spine model of the cone-horizontal cell synaptic circuitry is formulated. This model, a partial differential equation system, incorporates both the GABA and ephaptic feedback mechanisms. Simulation results, in comparison with experiments, indicate that the ephaptic mechanism is necessary in order for the model to capture the major spatial and temporal dynamics of the BIFE effect. In addition, simulations indicate that the GABA mechanism may play some minor modulation role.
ContributorsChang, Shaojie (Author) / Baer, Steven M. (Thesis advisor) / Gardner, Carl L (Thesis advisor) / Crook, Sharon M (Committee member) / Kuang, Yang (Committee member) / Ringhofer, Christian (Committee member) / Arizona State University (Publisher)
Created2012
150809-Thumbnail Image.png
Description
Dopamine (DA) is a neurotransmitter involved in attention, goal oriented behavior, movement, reward learning, and short term and working memory. For the past four decades, mathematical and computational modeling approaches have been useful in DA research, and although every modeling approach has limitations, a model is an efficient way to

Dopamine (DA) is a neurotransmitter involved in attention, goal oriented behavior, movement, reward learning, and short term and working memory. For the past four decades, mathematical and computational modeling approaches have been useful in DA research, and although every modeling approach has limitations, a model is an efficient way to generate and explore hypotheses. This work develops a model of DA dynamics in a representative, single DA neuron by integrating previous experimental, theoretical and computational research. The model consists of three compartments: the cytosol, the vesicles, and the extracellular space and forms the basis of a new mathematical paradigm for examining the dynamics of DA synthesis, storage, release and reuptake. The model can be driven by action potentials generated by any model of excitable membrane potential or even from experimentally induced depolarization voltage recordings. Here the model is forced by a previously published model of the excitable membrane of a mesencephalic DA neuron in order to study the biochemical processes involved in extracellular DA production. After demonstrating that the model exhibits realistic dynamics resembling those observed experimentally, the model is used to examine the functional changes in presynaptic mechanisms due to application of cocaine. Sensitivity analysis and numerical studies that focus on various possible mechanisms for the inhibition of DAT by cocaine provide insight for the complex interactions involved in DA dynamics. In particular, comparing numerical results for a mixed inhibition mechanism to those for competitive, non-competitive and uncompetitive inhibition mechanisms reveals many behavioral similarities for these different types of inhibition that depend on inhibition parameters and levels of cocaine. Placing experimental results within this context of mixed inhibition provides a possible explanation for the conflicting views of uptake inhibition mechanisms found in experimental neuroscience literature.
ContributorsTello-Bravo, David (Author) / Crook, Sharon M (Thesis advisor) / Greenwood, Priscilla E (Thesis advisor) / Baer, Steven M. (Committee member) / Castaneda, Edward (Committee member) / Castillo-Chavez, Carlos (Committee member) / Arizona State University (Publisher)
Created2012
156080-Thumbnail Image.png
Description
While techniques for reading DNA in some capacity has been possible for decades,

the ability to accurately edit genomes at scale has remained elusive. Novel techniques

have been introduced recently to aid in the writing of DNA sequences. While writing

DNA is more accessible, it still remains expensive, justifying the increased interest in

in

While techniques for reading DNA in some capacity has been possible for decades,

the ability to accurately edit genomes at scale has remained elusive. Novel techniques

have been introduced recently to aid in the writing of DNA sequences. While writing

DNA is more accessible, it still remains expensive, justifying the increased interest in

in silico predictions of cell behavior. In order to accurately predict the behavior of

cells it is necessary to extensively model the cell environment, including gene-to-gene

interactions as completely as possible.

Significant algorithmic advances have been made for identifying these interactions,

but despite these improvements current techniques fail to infer some edges, and

fail to capture some complexities in the network. Much of this limitation is due to

heavily underdetermined problems, whereby tens of thousands of variables are to be

inferred using datasets with the power to resolve only a small fraction of the variables.

Additionally, failure to correctly resolve gene isoforms using short reads contributes

significantly to noise in gene quantification measures.

This dissertation introduces novel mathematical models, machine learning techniques,

and biological techniques to solve the problems described above. Mathematical

models are proposed for simulation of gene network motifs, and raw read simulation.

Machine learning techniques are shown for DNA sequence matching, and DNA

sequence correction.

Results provide novel insights into the low level functionality of gene networks. Also

shown is the ability to use normalization techniques to aggregate data for gene network

inference leading to larger data sets while minimizing increases in inter-experimental

noise. Results also demonstrate that high error rates experienced by third generation

sequencing are significantly different than previous error profiles, and that these errors can be modeled, simulated, and rectified. Finally, techniques are provided for amending this DNA error that preserve the benefits of third generation sequencing.
ContributorsFaucon, Philippe Christophe (Author) / Liu, Huan (Thesis advisor) / Wang, Xiao (Committee member) / Crook, Sharon M (Committee member) / Wang, Yalin (Committee member) / Sarjoughian, Hessam S. (Committee member) / Arizona State University (Publisher)
Created2017
156933-Thumbnail Image.png
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

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.
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
154838-Thumbnail Image.png
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

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.
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