Atlantic Mathematical Biology Workshop Schedule

All talks are in Singer Hall Room 361. Enter the building through the main doors (uphill side) and go up the stairs.

Monday, June 7

8:45 Welcome and Opening
9:00 Plenary Lecture: Gerda de Vries
10:00 Coffee, Tea and Biscuits
10:30 Jeff Picka Random Trajectories in Biological Systems
10:50 Richard Karsten Determining an appropriate Model for Mites on Apple Trees
11:00 Yu Jin Traveling waves for a population model
11:30 Holger Teismann Reproductive value revisited, with an eye on pest control
11:50 Lunch, provided
2:00 Olga Vasilyeva Population dynamics in rivers: steady states and competition
2:20 Plenary Lecture: Andrea Locke
3:20 More Coffee, Tea and Biscuits
3:40 David Drolet Effects of climate warming on metapopulation dynamics of the intertidal amphipod Corophium volutator in the upper Bay of Fundy
4:00 Shengqiang Liu Competition Exclusion in Chemostat Models with Delays
4:20 Humberto Muñoz Novel Polypeptides for the Hydrogel-Based Encapsulation of Mouse ES Cells
4:40 Majid Jaberi-Douraki Dynamics of a honeybee model on regulation of work distribution
5:00 TBA
5:20 Adjourn

Tuesday, June 8

9:00 Plenary Lecture: Abba Gumel
10:00 Coffee, Tea and Biscuits
10:20 Joseph Apaloo Evolutionary Game Theory: ESS, Convergence Stability, and NIS
10:40 Lisa Kanary Model of the Invasive Tunicate Population Spread in Hillsborough Bay, PEI
11:00 Fang Yu
11:20 Huaichun Wang An improved /w/ statistics for covarion and heterotachy tests
11:40 Yijun Lou Modeling malaria control by integrated management strategies
12:00 Xi Hu
12:20 Closing

Abstracts

Joseph Apaloo

Evolutionary Game Theory: ESS, Convergence Stability, and NIS

Ever since Maynard Smith and Price's (1973) pioneering work, evolutionary game theory has advanced from matrix to continuous games, from single to multiple species, from scalar to vector-valued strategies, and from static analyses to adaptive dynamics. Essentially three related stability concepts underlie the theory for predicting the outcomes of natural selection: ESS, convergence stability and NIS. Here, we use the fitness generating function concept (G -function) and adaptive landscapes to illustrate these stability concepts. Explicit consideration of population dynamics and strategy dynamics permits a clear appreciation for each stability concept in terms of ecology (population dynamics) and evolution (strategy dynamics). We conclude by tabulating how most of the evolutionary stability acronyms, definitions and terminologies reduce to one or several aspects of ESS, convergence stability and NIS


Gerda de Vries

Quantitative Analysis of Single Particle Tracking Experiments: Applying Ecological Methods in Cell Biology

A commonly used experimental technique to study the movement of biomolecules in the cell membrane is Single Particle Tracking (SPT). SPT involves tagging biomolecules with a fluorescent label and observing and recording their trajectories over time. A diffusion coefficient then can be extracted from the data from mean square displacement calculations. Although the diffusion coefficient provides an overall measure of mobility, it does not provide insight into the underlying heterogeneity of the membrane environment.

Since the method of data collection from individual biomolecules is analogous to that from individual animals, we propose to use methods from ecology to provide spatial insight. Ecologists regularly quantify animal movement using the concepts of correlated random walk, net squared displacement, and first-passage time. In this talk, we will demonstrate the applicability of these methods in the context of cell biology. In particular, we will show how we can distinguish biomolecule trajectories undergoing a correlated random walk from those that do not, and how we can identify the presence of transient confinement zones in molecular diffusion.


Abba Gumel

Mathematics of the Transmission Dynamics of Chlamydia Trachomatis

Chlamydia trachomatis is a sexually-transmitted disease that continues to inflict major public health and socio-economic burden around the world. The talk will address the problem of the design and analysis of appropriate models for the transmission dynamics of Chlamydia in a heterosexual population. The effect of risk of acquiring or transmitting infection will be assessed.


Jeff Picka

Random Trajectories in Biological Systems

To model the changes in a complex biological system over time, a dynamical system is often used to model the interactions among many individual cells, plants, or animals. To validate any model of this kind, it is necessary to incorporate into the model any uncertainties about the exact mechanism of interactions and about environmental conditions. The conditions required to achieve different degrees of validation will be discussed.


Lisa Kanary1, Andrea Locke2, James Watmough1, Joël Chassé2 and Daniel Bourque2

1 University of New Brunswick, P.O. Box 4400, Fredericton, NB, E3B 5A3;
2 Gulf Fisheries Centre, Fisheries and Oceans Canada, P.O. Box 5030, Moncton, NB, E1C 9B6

Model of the Invasive Tunicate Population Spread in Hillsborough Bay, PEI

Invasive tunicates have been a problem for mussel farmers in Prince Edward Island since 1997. In this presentation we develop a model to ascertain the conditions necessary for an infestation to spread within Hillsborough Bay from Charlettetown Harbour to Nine Mile Creek. Nine Mile Creek where the majority of mussel seed production and distribution is based.

A system of partial differential equations models the tunicate population spread, taking into consideration parameters based on larval mortality, reproductive and settling rates. Our main result is a study of the conditions necessary for the persistence and spread of tunicates based on the spacing of intermediate settlement substrate. Advection (ocean current) data provided by the Department of Fisheries and Oceans and conditions derived from previous analysis are incorporated into a numerical simulation of the PDE system to demonstrate how the varying amounts of intermediate substrate affect the persistence of a tunicate invasion.


Majid Jaberi-Douraki

Dynamics of a honeybee model on regulation of work distribution

In this talk, we study age-related activities of honeybees (Apis mellifera L.) first by illustrating the life cycle of honeybee which exhibits a combination of individual traits and social cooperation, and then by constructing an age structured model given by a system of difference equations as follows

Hn+1=F(Hn)=(A(Hn)-D)Hn+B       n 1
where the initial conditions H1 is positive, B represents the broods, Hn denotes the population of bees divided into (k+1) classes with respect to their ages, and A(Hn) is a lower triangular matrix representing progression or maturation to the next stage of their life involving hive bees and foraging bees, D indicates the matrix of death rate with subdiagonal entries mi and others zero. In fact the distribution of most activities and behavioural/physiological maturation are regulated through primer and releaser pheromones from the queen, worker bees, and brood; we include such dynamical organization, as important characteristics for the honey bee, by introducing numerous parameters in the model. Using extensive numerical simulations, we conjecture that this system has a unique positive steady state which is globally asymptotically stable. We are able to prove analytically the global stability (basin of attraction) of feasible solution of such dynamical system with biological relevant special cases.


David Drolet and Myriam A Barbeau

Effects of climate warming on metapopulation dynamics of the intertidal amphipod Corophium volutator in the upper Bay of Fundy

We are developing population models to investigate the impact of rising temperature on the dominant invertebrate on mudflats of the Bay of Fundy, the amphipod Corophium volutator. We are using stage-based matrix models using transition probabilities obtained from the literature. The transition probabilities are temperature-dependant and affected by a variety of factors known to influence C. volutator. Single population models will be tested against time series of densities of the different stages at 10 different sites that are currently being sampled, with parameters representing local conditions. Sensitivity analysis will be performed to assess the relative importance of the parameters at the different sites. The models will then be used to evaluate population responses to predicted temperature changes in the area. Finally, the models for the different sites will be integrated by allowing dispersal between mudflats (based on a study that is presently conducted) and the effect of increasing temperature on metapopulation dynamics will be examined.


Humberto Muñoz

Department of Mathematics

Eduardo Martinez-Ceballos

Department of Biological Sciences

Southern University and A&M College

Novel Polypeptides for the Hydrogel-Based Encapsulation of Mouse ES Cells

This work proposes the creation of new polypeptides for the hydrogel-based encapsulation of mouse embryonic stem cells (ES cells). The design of these novel polypeptides will be based on relevant domains of proteins with known ES differentiating activities. The polypeptide-containing hydrogels will be characterized based on their physicochemical properties and their effect on the directed differentiation of ES cells along the three specific cell lineages: ectoderm, mesoderm, and endoderm. Encapsulation conditions will be optimized by applying mathematical algorithms based on biological responses and the physicochemical properties of the hydrogels. The data analysis involves nonlinear optimization and regression analysis methods.


Richard Karsten, Holger Teismann, Rebecca Hammond, Jeff Franklin, Mike Hardman

Determining an appropriate Model for Mites on Apple Trees

In Nova Scotia's Annapolis Valley, infestations of mites, primarily the European red mite can cause serious economic losses in apple orchards. Efficient control of these mites, either by spraying pesticides or introducing predators, is an active area of research. In this talk, we discuss our attempts to determine an appropriate mathematical model that can capture the essential characteristics of the population dynamics. Experiments and observations have indicated that the fecundity and survival of individuals in the populations may depend on a long list of factors including age, population density, predation, temperature, day of year, rainfall and health of the host tree. Including all these factors in a model results in a model with a large number of parameters that must be determined from the observations. However, accurate observations of mites in orchards are difficult to obtain. Observed time series of mite populations exhibit strong generational cycles within a larger seasonal cycle. But the number of complete time series are few in number and the observations are limited to a few age groups and few observations during the season. We examine a range of possible models for the mite population including partial differential equations in age and time, delay differential equations, systems of ODEs, and Lesley-Matrix models. Using "goodness of fit" measurements we determine the models with the fewest parameters that best replicate the observations. We also generate artificial time series to determine if improved observations would change our results.


Holger Teismann, Richard Karsten, Angela Vogels

Reproductive value revisited, with an eye on pest control

While Fisher's reproductive value is a well-established notion in mathematical demography, it is not necessarily easily understood. In this talk we will review the basic theory behind the reproductive value, and revisit some of the motivation and explanations that have been given. We will then try to illuminate the connection to optimal control problems, which we consider in the form of pest control problems. This talk is 'survey-ish' in nature; its latter part is inspired by recent work by M. Kuhn et al. (Theor. Pop. Biol. 77, 164-170, 2010).


Shengqiang Liu

Harbin Institute of Technology, China

Competition Exclusion in Chemostat Models with Delays

We consider the global asymptotic behaviors of three exploitative competition models among n-species for a chemostat with general response functions and differential removal rates. Three different classes of delays in describing the conversion process of nutrient to new cells are studied separately: 1. infinite-distributed type; 2. finite-distributed type; 3. discrete type. By carefully constructing the proper Lyapunov functionals, it is shown that competitive exclusions hold for these models with a series of general growth response functions. Previous results in the literatures are significantly improved and extended.

Joint work with Haitao Song, Lin Wang and Xinxin Wang


Olga Vasilyeva

University of Ottawa

Population dynamics in rivers: steady states and competition

We study diffusion-reaction-advection models describing population dynamics of logistically growing aquatic organisms subject to a constant downsteam drift, with reflecting upstream and outflow downstream boundary conditions. In the case of a single population, we determine conditions for existence, uniqueness and stability of non-trivial steady state solutions, and we analyze the dependence of such solutions on advection speed, growth rate and length of the habitat. Such analysis offers a possible explanation of the "drift paradox" in our context. In the case of two competing species, we use a diffusion-advection version of the Lotka-Volterra competition model. Combining numerical and and analytical techniques, we describe the effect of advection on competition outcomes.


Huaichun Wang

Department of Mathematics and Statistics

Dalhousie University

An improved /w/ statistics for covarion and heterotachy tests

The /w/ statistics is a quantity that compares amino acid substitution patterns between two monophyletic groups of protein sequences. It is defined as the difference between the fraction of varied sites in both groups and the fraction of varied sites in each group. The /w/ test has been used to distinguish a covarion process from equal rates and rates variation across sites processes. Using simulation we show that the /w/ test is effective only for small datasets and for datasets that have small substitution rates in the groups. Using site entropy as a measure of variability of the sequence site we generalize the /w/ statistics by assigning those sites that are varied in both groups but have a large entropy difference to those sites that are varied only in one group, thus modifying the fractions of varied sites in each group and in both groups. We show that the /w/’ test has more power to detect heterotachy processes (covarion, bivariate rate shift and branch length mixture) in large and variable data. We also show that the difference in correlation coefficient of the site entropies between the two groups can be used for separating the heterotachy models and the rate-across-sites model.


Yijun Lou

Modeling malaria control by integrated management strategies

Malaria creates serious health and economic problems which call for integrated management strategies to disrupt interactions among the vector mosquitoes, parasite, and humans. In order to evaluate the potential impacts of control methods on malaria transmission, we propose a mathematical model which consists of five ordinary differential equations. The global dynamics of this model is rigorously analyzed. We derive biologically plausible and insightful quantities (reproduction numbers) that completely determine the long-term behavior of solutions. Strategies for malaria control are analyzed to demonstrate the biological application of our theoretical work.


Yu Jin

University of Alberta

Traveling waves for a population model

In this talk, I will give a class of coupled cooperative reaction-diffusion systems, in which one population (or subpopulation) diffuses while the other is sedentary. The shooting method is used to prove the existence of the bistable travelling wave. Global attractivity (with phase shift) and uniqueness (up to translation) of the traveling wave are obtained via the dynamical system approach. The results are applied to some specific examples of reaction-diffusion population models.