Mathematical Biology

The evolution and maintenance of cooperation is a fundamental problem in evolutionary biology. Because cooperative behaviors impose a cost, Cooperators are vulnerable to exploitation by Defectors that do not pay the cost to cooperate but still benefit from the cooperation of others. The bacteriophage $\Phi_6$ exhibits cooperative and defective phenotypes in infection: during replication, phages produce essential proteins in the host cell cytoplasm. Coinfection between multiple phages is possible. A given phage cannot guarantee exclusive access to its own proteins, so Cooperators contribute to the common pool of proteins while Defectors contribute less and instead appropriate proteins from Cooperators. Previous experimental work found that $\Phi_6$ was trapped in a prisoner's dilemma, predicting that the cooperative phenotype should disappear. Here we propose that environmental feedback, or interplay between phage and host densities, can maintain cooperation in $\Phi_6$ populations by modulating the rate of co-infection and shifting the advantages of cooperation vs. defection. We build and analyze an ODE model and find that for a wide range of parameter values, environmental feedback allows Cooperation to survive.

Evolutionary game theory is a discipline devoted to studying populations of individuals that are subject to evolutionary pressures, and whose success generally depends on the composition of the population. In biological contexts, individuals could be molecules, simple organisms or animals, and evolutionary pressures often take the form of natural selection and mutations. In socioeconomic contexts, individuals could be humans, firms or other institutions, and evolutionary pressures often derive from competition for scarce resources and experimentation.

In this talk I will give a very basic introduction to agent-based evolutionary game theory, a bottom-up approach to modelling and analyzing these systems. The defining feature of this modelling approach is that the individual units of the system and their interactions are explicitly and individually represented in the model. The models thus defined can be usefully formalized as stochastic processes, whose dynamics can be explored using computer simulation and approximated using various mathematical theories.

For over a century, scientists have studied striking spatiotemporal patterns during the continual tooth replacement of reptiles. Aside from the compelling aesthetics of this phenomenon, it is thought that understanding the underlying mechanisms may provide the insight required to trigger adult tooth replacement in humans. Theoretical frameworks have long been proposed to understand the rules behind the observed spatiotemporal order, but have only been analyzed mathematically more recently. Starting from Edmund's observations in crocodiles and proposed theory of replacement waves, we show how a simple model consisting of a row of non-interacting phase oscillators predicts several experimental observations. Next, inspired by the hypothesis put forth by Osborn, we consider a variation of the phase model with ODEs that account for mutual inhibition between tooth sites, and use continuation methods to thoroughly search parameter space for experimentally validated solutions. We then extend the model to a PDE that explicitly accounts for the diffusion of inhibitory signals between teeth, yielding some novel solution types. Using continuation methods once again, we delineate parameter regimes with solutions that closely resemble experimental observations in leopard geckos.

The questions of how healthy colonic crypts maintain their size under the rapid cell turnover in intestinal epithelium, and how homeostasis is disrupted by driver mutations, are central to understanding colorectal tumorigenesis. We propose a three-type stochastic branching process, which accounts for stem, transit-amplifying (TA) and fully differentiated (FD) cells, to model the dynamics of cell populations residing in colonic crypts. Our model is simple in its formulation, allowing us to estimate all but one of the model parameters from the literature. Fitting the single remaining parameter, we find that model results agree well with data from healthy human colonic crypts, capturing the considerable variance in population sizes observed experimentally. Importantly, our model predicts a steady-state population in healthy colonic crypts for relevant parameter values. We show that APC and KRAS mutations, the most significant early alterations leading to colorectal cancer, result in increased steady-state populations in mutated crypts, in agreement with experimental results. Finally, our model predicts a simple condition for unbounded growth of cells in a crypt, corresponding to colorectal malignancy. This is predicted to occur when the division rate of TA cells exceeds their differentiation rate, with implications for therapeutic cancer prevention strategies.

While mathematical models have classically been used in the study of physics and engineering, recently, they have become important tools in other fields such as biology, ecology, and sociology. In this talk I will discuss the use of partial differential equations and dynamical systems to shed light onto social and ecological phenomena. In the first part of this talk, we will focus on an Ecological application. For an efficient wildlife management plan, it is important that we understand (1) why animals move as they do and (2) what movement strategies are robust. I will discuss how reaction-advection-diffusion models can help us shed light into these two issues. The second part of the talk will focus on social applications. I will present a few models in the study of gentrification, urban crime, and protesting activity and discuss how theoretical and numerical analysis have provided intuition into these different social phenomena. Moreover, I will also point out the many benefits of utilizing a mathematical framework when data is not available.