Magnetic resonance imaging (MRI) data of metastatic brain cancer patients at the Barrow Neurological Institute sparked interest in the radiology department due to the possibility that tumor size distributions might mimic a power law or an exponential distribution. In order to consider the question regarding the growth trends of metastatic brain tumors, this thesis analyzes the volume measurements of the tumor sizes from the BNI data and attempts to explain such size distributions through mathematical models. More specifically, a basic stochastic cellular automaton model is used and has three-dimensional results that show similar size distributions of those of the BNI data. Results of the models are investigated using the likelihood ratio test suggesting that, when the tumor volumes are measured based on assuming tumor sphericity, the tumor size distributions significantly mimic the power law over an exponential distribution.

The most advanced social insects, the eusocial insects, form often large societies in which there is reproductive division of labor, queens and workers, have overlapping generations, and cooperative brood care where daughter workers remain in the nest with their queen mother and care for their siblings. The eusocial insects are composed of representative species of bees and wasps, and all species of ants and termites. Much is known about their organizational structure, but remains to be discovered.

The success of social insects is dependent upon cooperative behavior and adaptive strategies shaped by natural selection that respond to internal or external conditions. The objective of my research was to investigate specific mechanisms that have helped shaped the structure of division of labor observed in social insect colonies, including age polyethism and nutrition, and phenomena known to increase colony survival such as egg cannibalism. I developed various Ordinary Differential Equation (ODE) models in which I applied dynamical, bifurcation, and sensitivity analysis to carefully study and visualize biological outcomes in social organisms to answer questions regarding the conditions under which a colony can survive. First, I investigated how the population and evolutionary dynamics of egg cannibalism and division of labor can promote colony survival. I then introduced a model of social conflict behavior to study the inclusion of different response functions that explore the benefits of cannibalistic behavior and how it contributes to age polyethism, the change in behavior of workers as they age, and its biological relevance. Finally, I introduced a model to investigate the importance of pollen nutritional status in a honeybee colony, how it affects population growth and influences division of labor within the worker caste. My results first reveal that both cannibalism and division of labor are adaptive strategies that increase the size of the worker population, and therefore, the persistence of the colony. I show the importance of food collection, consumption, and processing rates to promote good colony nutrition leading to the coexistence of brood and adult workers. Lastly, I show how taking into account seasonality for pollen collection improves the prediction of long term consequences.

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.

Rewired biological pathways and/or rewired microRNA (miRNA)-mRNA interactions might also influence the activity of biological pathways. Here, rewired biological pathways is defined as differential (rewiring) effect of genes on the topology of biological pathways between controls and cases. Similarly, rewired miRNA-mRNA interactions are defined as the differential (rewiring) effects of miRNAs on the topology of biological pathways between controls and cases. In the dissertation, it is discussed that how rewired biological pathways (Chapter 1) and/or rewired miRNA-mRNA interactions (Chapter 2) aberrantly influence the activity of biological pathways and their association with disease.

This dissertation proposes two PageRank-based analytical methods, Pathways of Topological Rank Analysis (PoTRA) and miR2Pathway, discussed in Chapter 1 and Chapter 2, respectively. PoTRA focuses on detecting pathways with an altered number of hub genes in corresponding pathways between two phenotypes. The basis for PoTRA is that the loss of connectivity is a common topological trait of cancer networks, as well as the prior knowledge that a normal biological network is a scale-free network whose degree distribution follows a power law where a small number of nodes are hubs and a large number of nodes are non-hubs. However, from normal to cancer, the process of the network losing connectivity might be the process of disrupting the scale-free structure of the network, namely, the number of hub genes might be altered in cancer compared to that in normal samples. Hence, it is hypothesized that if the number of hub genes is different in a pathway between normal and cancer, this pathway might be involved in cancer. MiR2Pathway focuses on quantifying the differential effects of miRNAs on the activity of a biological pathway when miRNA-mRNA connections are altered from normal to disease and rank disease risk of rewired miRNA-mediated biological pathways. This dissertation explores how rewired gene-gene interactions and rewired miRNA-mRNA interactions lead to aberrant activity of biological pathways, and rank pathways for their disease risk. The two methods proposed here can be used to complement existing genomics analysis methods to facilitate the study of biological mechanisms behind disease at the systems-level.

Diabetes is a disease characterized by reduced insulin action and secretion, leading to elevated blood glucose. In the 1990s, studies showed that intravenous injection of fatty acids led to a sharp negative response in insulin action that subsided hours after the injection. The molecule associated with diminished insulin signalling response was a byproduct of fatty acids, diacylglycerol. This dissertation is focused on the formulation of a model built around the known mechanisms of glucose and fatty acid storage and metabolism within myocytes, as well as downstream effects of diacylglycerol on insulin action. Data from euglycemic-hyperinsulinemic clamp with fatty acid infusion studies are used to validate the qualitative behavior of the model and estimate parameters. The model closely matches clinical data and suggests a new metric to determine quantitative measurements of insulin action downregulation. Analysis and numerical simulation of the long term, piecewise smooth system of ordinary differential equations demonstrates a discontinuous bifurcation implicating nutrient excess as a driver of muscular insulin resistance.

Predicting resistant prostate cancer is critical for lowering medical costs and improving the quality of life of advanced prostate cancer patients. I formulate, compare, and analyze two mathematical models that aim to forecast future levels of prostate-specific antigen (PSA). I accomplish these tasks by employing clinical data of locally advanced prostate cancer patients undergoing androgen deprivation therapy (ADT). I demonstrate that the inverse problem of parameter estimation might be too complicated and simply relying on data fitting can give incorrect conclusions, since there is a large error in parameter values estimated and parameters might be unidentifiable. I provide confidence intervals to give estimate forecasts using data assimilation via an ensemble Kalman Filter. Using the ensemble Kalman Filter, I perform dual estimation of parameters and state variables to test the prediction accuracy of the models. Finally, I present a novel model with time delay and a delay-dependent parameter. I provide a geometric stability result to study the behavior of this model and show that the inclusion of time delay may improve the accuracy of predictions. Also, I demonstrate with clinical data that the inclusion of the delay-dependent parameter facilitates the identification and estimation of parameters.

Glioblastoma multiforme (GBM) is a malignant, aggressive and infiltrative cancer of the central nervous system with a median survival of 14.6 months with standard care. Diagnosis of GBM is made using medical imaging such as magnetic resonance imaging (MRI) or computed tomography (CT). Treatment is informed by medical images and includes chemotherapy, radiation therapy, and surgical removal if the tumor is surgically accessible. Treatment seldom results in a significant increase in longevity, partly due to the lack of precise information regarding tumor size and location. This lack of information arises from the physical limitations of MR and CT imaging coupled with the diffusive nature of glioblastoma tumors. GBM tumor cells can migrate far beyond the visible boundaries of the tumor and will result in a recurring tumor if not killed or removed. Since medical images are the only readily available information about the tumor, we aim to improve mathematical models of tumor growth to better estimate the missing information. Particularly, we investigate the effect of random variation in tumor cell behavior (anisotropy) using stochastic parameterizations of an established proliferation-diffusion model of tumor growth. To evaluate the performance of our mathematical model, we use MR images from an animal model consisting of Murine GL261 tumors implanted in immunocompetent mice, which provides consistency in tumor initiation and location, immune response, genetic variation, and treatment. Compared to non-stochastic simulations, stochastic simulations showed improved volume accuracy when proliferation variability was high, but diffusion variability was found to only marginally affect tumor volume estimates. Neither proliferation nor diffusion variability significantly affected the spatial distribution accuracy of the simulations. While certain cases of stochastic parameterizations improved volume accuracy, they failed to significantly improve simulation accuracy overall. Both the non-stochastic and stochastic simulations failed to achieve over 75% spatial distribution accuracy, suggesting that the underlying structure of the model fails to capture one or more biological processes that affect tumor growth. Two biological features that are candidates for further investigation are angiogenesis and anisotropy resulting from differences between white and gray matter. Time-dependent proliferation and diffusion terms could be introduced to model angiogenesis, and diffusion weighed imaging (DTI) could be used to differentiate between white and gray matter, which might allow for improved estimates brain anisotropy.

Prostate cancer is the second most common kind of cancer in men. Fortunately, it has a 99% survival rate. To achieve such a survival rate, a variety of aggressive therapies are used to treat prostate cancers that are caught early. Androgen deprivation therapy (ADT) is a therapy that is given in cycles to patients. This study attempted to analyze what factors in a group of 79 patients caused them to stick with or discontinue the treatment. This was done using naïve Bayes classification, a machine-learning algorithm. The usage of this algorithm identified high testosterone as an indicator of a patient persevering with the treatment, but failed to produce statistically significant high rates of prediction.

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.

In this research we consider stochastic models of Glioblastoma Multiforme brain tumors. We first look at a model by K. Swanson et al., which describes the dynamics as random diffusion plus deterministic logistic growth. We introduce a stochastic component in the logistic growth in the form of a random growth rate defined by a Poisson process. We show that this stochastic logistic growth model leads to a more accurate evaluation of the tumor growth compared its deterministic counterpart. We also discuss future plans to incorporate individual patient geometry, extend the model to three dimensions and to incorporate effects of different treatments into our model, in collaboration with a local hospital.