On a Reproductive Stage-Structured Model for the Population Dynamics of the Malaria Vector

A reproductive stage-structured deterministic differential equation model for the population dynamics of the human malaria vector is derived and analysed. The model captures the gonotrophic and behavioural life characteristics of the female Anopheles sp. mosquito and takes into consideration the fact that for the purposes of reproduction, the female Anopheles sp. mosquito must visit and bite humans (or animals) to harvest necessary proteins from blood that it needs for the development of its eggs. Focusing on mosquitoes that feed exclusively on humans, our results indicate the existence of a threshold parameter, the vectorial reproduction number, whose size increases with increasing number of gonotrophic cycles, and is also affected by the female mosquito’s birth rate, its attraction and visitation rate to human residences, and its contact rate with humans. A stability analysis of the model indicates that the mosquito can establish itself in the environment if and only if the value of the vectorial reproduction number exceeds unity and that mosquito eradication is possible if the vectorial reproduction number is less than unity, since, then, the trivial steady state which always exist is unique and is globally and asymptotically stable. When a persistent vector population steady state exists, it is locally and asymptotically stable for a range of reproduction numbers, but can also be driven to instability via a Hopf bifurcation as the reproduction number increases further away from unity. The model derivation identifies and characterizes control parameters relating to activities such as human-mosquito contact and the mosquito’s survival chances between blood meals and egg laying. Our results show that the total mosquito population size increases with increasing number of gonotrophic cycles. Therefore understanding the fundamental aspects of the mosquito’s behaviour provides a pathway for the study of human-mosquito contact and mosquito population control. Control of the mosquito population densities would ultimately lead to malaria control.     

Discrete three-stage population model: persistence and global stability results

A general three-stage discrete-time population model is studied. The inherent net reproductive number for this model is derived. Global stability of the origin is established provided that the inherent net reproductive number is less than one. If it is larger than one the existence of a unique positive fixed point is proved and the persistence of the system is established. Finally, for certain parameter ranges global stability of the positive fixed point is proved.     

Discrete three-stage population model: persistence and global stability results

A general three-stage discrete-time population model is studied. The inherent net reproductive number for this model is derived. Global stability of the origin is established provided that the inherent net reproductive number is less than one. If it is larger than one the existence of a unique positive fixed point is proved and the persistence of the system is established. Finally, for certain parameter ranges global stability of the positive fixed point is proved.     

Modelling the Transmission Dynamics and Control of the Novel 2009 Swine Influenza (H1N1) Pandemic

The paper presents a deterministic compartmental model for the transmission dynamics of swine influenza (H1N1) pandemic in a population in the presence of an imperfect vaccine and use of drug therapy for confirmed cases. Rigorous analysis of the model, which stratifies the infected population in terms of their risk of developing severe illness, reveals that it exhibits a vaccine-induced backward bifurcation when the associated reproduction number is less than unity. The epidemiological consequence of this result is that the effective control of H1N1, when the reproduction number is less than unity, in the population would then be dependent on the initial sizes of the subpopulations of the model. For the case where the vaccine is perfect, it is shown that having the reproduction number less than unity is necessary and sufficient for effective control of H1N1 in the population (in such a case, the associated disease-free equilibrium is globally asymptotically stable). The model has a unique endemic equilibrium when the reproduction number exceeds unity. Numerical simulations of the model, using data relevant to the province of Manitoba, Canada, show that it reasonably mimics the observed H1N1 pandemic data for Manitoba during the first (Spring) wave of the pandemic. Further, it is shown that the timely implementation of a mass vaccination program together with the size of the Manitoban population that have preexisting infection-acquired immunity (from the first wave) are crucial to the magnitude of the expected burden of disease associated with the second wave of the H1N1 pandemic. With an estimated vaccine efficacy of approximately 80%, it is projected that at least 60% of Manitobans need to be vaccinated in order for the effective control or elimination of the H1N1 pandemic in the province to be feasible. Finally, it is shown that the burden of the second wave of H1N1 is expected to be at least three times that of the first wave, and that the second wave would last until the end of January or early February, 2010.     

An SIRVS epidemic model with pulse vaccination strategy

The aim of this paper is to analyze an SIRVS epidemic model in which pulse vaccination strategy (PVS) is included. We are interested in finding the basic reproductive number of the model which determine whether or not the disease dies out. The global attractivity of the disease-free periodic solution (DFPS for short) is obtained when the basic reproductive number is less than unity. The disease is permanent when the basic reproductive number is greater than unity, i.e., the epidemic will turn out to endemic. Our results indicate that the disease will go to extinction when the vaccination rate reaches some critical value.     

Global Stability of Equilibria in a Two-Sex HPV Vaccination Model

Human papillomavirus (HPV) is the primary cause of cervical carcinoma and its precursor lesions, and is associated with a variety of other cancers and diseases. A prophylactic quadrivalent vaccine against oncogenic HPV 16/18 and warts-causing genital HPV 6/11 types is currently available in several countries. Licensure of a bivalent vaccine against oncogenic HPV 16/18 is expected in the near future. This paper presents a two-sex, deterministic model for assessing the potential impact of a prophylactic HPV vaccine with several properties. The model is based on the susceptible-infective-removed (SIR) compartmental structure. Important epidemiological thresholds such as the basic and effective reproduction numbers and a measure of vaccine impact are derived. We find that if the effective reproduction number is greater than unity, there is a locally unstable infection-free equilibrium and a unique, globally asymptotically stable endemic equilibrium. If the effective reproduction number is less than unity, the infection-free equilibrium is globally asymptotically stable, and HPV will be eliminated.     

Extinction times and moment closure in the stochastic logistic process

We investigate the statistics of extinction times for an isolated population, with an initially modest number M of individuals, whose dynamics are controlled by a stochastic logistic process (SLP). The coefficient of variation in the extinction time V is found to have a maximum value when the death and birth rates are close in value. For large habitat size K we find that Vmax is of order K1/4 / M1/2, which is much larger than unity so long as M is small compared to K1/2. We also present a study of the SLP using the moment closure approximation (MCA), and discuss the successes and failures of this method. Regarding the former, the MCA yields a steady-state distribution for the population when the death rate is low. Although not correct for the SLP model, the first three moments of this distribution coincide with those calculated exactly for an adjusted SLP in which extinction is forbidden. These exact calculations also pinpoint the breakdown of the MCA as the death rate is increased.     

An age-structured within-host HIV-1 infection model with virus-to-cell and cell-to-cell transmissions

In this paper, a within-host HIV-1 infection model with virus-to-cell and direct cell-to-cell transmission and explicit age-since-infection structure for infected cells is investigated. It is shown that the model demonstrates a global threshold dynamics, fully described by the basic reproduction number. By analysing the corresponding characteristic equations, the local stability of an infection-free steady state and a chronic-infection steady state of the model is established. By using the persistence theory in infinite dimensional system, the uniform persistence of the system is established when the basic reproduction number is greater than unity. By means of suitable Lyapunov functionals and LaSalle's invariance principle, it is shown that if the basic reproduction number is less than unity, the infection-free steady state is globally asymptotically stable; if the basic reproduction number is greater than unity, the chronic-infection steady state is globally asymptotically stable. Numerical simulations are carried out to illustrate the feasibility of the theoretical results.     

Equilibria and dynamics of selection at a diallelic autosomal locus in the Nagylaki-Crow continuous model of a monoecious population

Let birth rates and death rates be constant, birth rates positive, fertilities additive, and each birth rate not larger than twice any other birth rate. Global convergence to equilibria is proved for the model in the title. There is at most one polymorphic equilibrium or there are a continuum of equilibria. The phase portraits are given. If there is a polymorphic equilibrium, then the largest negatively invariant set in the state space is a continuous curve connecting the two fixation equilibria. This curve coincides with the Hardy-Weinberg manifold exactly when the death rate is additive. Disregarding extinction, the polymorphic equilibria are the same for the continuous model as for the corresponding discrete model exactly when the death rate is additive.     

Mathematical Analysis of a Two Strain HIV/AIDS Model with Antiretroviral Treatment

A two strain HIV/AIDS model with treatment which allows AIDS patients with sensitive HIV-strain to undergo amelioration is presented as a system of non-linear ordinary differential equations. The disease-free equilibrium is shown to be globally asymptotically stable when the associated epidemic threshold known as the basic reproduction number for the model is less than unity. The centre manifold theory is used to show that the sensitive HIV-strain only and resistant HIV-strain only endemic equilibria are locally asymptotically stable when the associated reproduction numbers are greater than unity. Qualitative analysis of the model including positivity, boundedness and persistence of solutions are presented. The model is numerically analysed to assess the effects of treatment with amelioration on the dynamics of a two strain HIV/AIDS model. Numerical simulations of the model show that the two strains co-exist whenever the reproduction numbers exceed unity. Further, treatment with amelioration may result in an increase in the total number of infective individuals (asymptomatic) but results in a decrease in the number of AIDS patients. Further, analysis of the reproduction numbers show that antiretroviral resistance increases with increase in antiretroviral use.     

Malaria Drug Resistance: The Impact of Human Movement and Spatial Heterogeneity

Human habitat connectivity, movement rates, and spatial heterogeneity have tremendous impact on malaria transmission. In this paper, a deterministic system of differential equations for malaria transmission incorporating human movements and the development of drug resistance malaria in an \(n\) patch system is presented. The disease-free equilibrium of the model is globally asymptotically stable when the associated reproduction number is less than unity. For a two patch case, the boundary equilibria (drug sensitive-only and drug resistance-only boundary equilibria) when there is no movement between the patches are shown to be locally asymptotically stable when they exist; the co-existence equilibrium is locally asymptotically stable whenever the reproduction number for the drug sensitive malaria is greater than the reproduction number for the resistance malaria. Furthermore, numerical simulations of the connected two patch model (when there is movement between the patches) suggest that co-existence or competitive exclusion of the two strains can occur when the respective reproduction numbers of the two strains exceed unity. With slow movement (or low migration) between the patches, the drug sensitive strain dominates the drug resistance strain. However, with fast movement (or high migration) between the patches, the drug resistance strain dominates the drug sensitive strain.     

The general mixing of addicts and needles in a variable-infectivity needle-sharing environment

 In this paper we develop and analyse a model for the spread of HIV/AIDS amongst a population of injecting drug users. The model we discuss focuses on the transmission of HIV through the sharing of contaminated drug injection equipment and in particular we examine the mixing of addicts and needles when the AIDS incubation period is divided into three distinct infectious stages. The impact of this assumption is to greatly increase the complexity of the HIV transmission mechanism. We begin the paper with a brief literature review followed by the derivation of a model which incorporates three classes of infectious addicts and three classes of infectious needles and where a general probability structure is used to represent the interaction of addicts and needles of varying levels of infectivity. We find that if the basic reproductive number is less than or equal to unity then there exists a globally stable disease free equilibrium. The model possesses an endemic equilibrium solution if the basic reproductive number exceeds unity. We then conduct a brief simulation study of our model. We find that the spread of disease is heavily influenced by the way addicts and needles of different levels of infectivity interact. Received: 20 September 2001 / Revised version: 21 December 2001 / Published online: 17 May 2002     

Modelling time to population extinction when individual reproduction is autocorrelated

In nature, individual reproductive success is seldom independent from year to year, due to factors such as reproductive costs and individual heterogeneity. However, population projection models that incorporate temporal autocorrelations in individual reproduction can be difficult to parameterise, particularly when data are sparse. We therefore examine whether such models are necessary to avoid biased estimates of stochastic population growth and extinction risk, by comparing output from a matrix population model that incorporates reproductive autocorrelations to output from a standard age‐structured matrix model that does not. We use a range of parameterisations, including a case study using moose data, treating probabilities of switching reproductive class as either fixed or fluctuating. Expected time to extinction from the two models is found to differ by only small amounts (under 10%) for most parameterisations, indicating that explicitly accounting for individual reproductive autocorrelations is in most cases not necessary to avoid bias in extinction estimates.     

Adaptation of timing of life history traits and population dynamic responses to climate change in spatially structured populations

Changes in the seasonal timing of life history events are documented effects of climate change. We used a general model to study how dispersal and competitive interactions affect eco-evolutionary responses to changes in the temporal distribution of resources over the season. Specifically, we modeled adaptation of the timing of reproduction and population dynamic responses in two competing populations that disperse between two habitats characterized by an early and late resource peak. We investigated three scenarios of environmental change: (1) food peaks advance in both habitats, (2) in the late habitat only and (3) in the early habitat only. At low dispersal rates the evolutionarily stable timing of reproduction closely matched the local resource peak and the environmental change typically caused population decline. Larger dispersal rates rendered less intuitive eco-evolutionary population responses. First, dispersal caused mismatch between evolutionarily stable timing of reproduction and local resource peaks and as a result, reproductive output for subpopulations could increase as well as decrease when resource availability underwent temporal shifts. Second, population responses were contingent on competition between populations. This could accelerate population declines and cause extinctions or even reverse population trends from negative to positive compared to the low dispersal case. When dispersal rate was large and the early resource peak was advanced available niche space was reduced. Hence, even when a population survived the environmental change and obtained positive equilibrium population density, subsequent adaptation of competing populations could drive it to extinction due to convergent evolution and competitive exclusion. These results shed new light on the role of competition and dispersal for the evolution of timing of life history events and provide guidelines for understanding short and long-term population response to climate change.     

Assessing seasonal demographic covariation to understand environmental‐change impacts on a hibernating mammal

Natural populations are exposed to seasonal variation in environmental factors that simultaneously affect several demographic rates (survival, development and reproduction). The resulting covariation in these rates determines population dynamics, but accounting for its numerous biotic and abiotic drivers is a significant challenge. Here, we use a factor‐analytic approach to capture partially unobserved drivers of seasonal population dynamics. We use 40 years of individual‐based demography from yellow‐bellied marmots (Marmota flaviventer) to fit and project population models that account for seasonal demographic covariation using a latent variable. We show that this latent variable, by producing positive covariation among winter demographic rates, depicts a measure of environmental quality. Simultaneously, negative responses of winter survival and reproductive‐status change to declining environmental quality result in a higher risk of population quasi‐extinction, regardless of summer demography where recruitment takes place. We demonstrate how complex environmental processes can be summarized to understand population persistence in seasonal environments.     

Early is better: seasonal egg fitness and timing of reproduction in a zooplankton life-history model

Timing of reproduction influences future prospects of offspring and therefore the reproductive value of parents. Early offspring are often more valuable than later ones when food availability and predation risk fluctuate seasonally. Marine zooplankton have evolved a diversity of life history strategies in response to seasonality. We present a state-dependent life history model for the annual and herbivorous high-latitude copepod Calanoides acutus . Individuals are characterised by four states; developmental stage, structural size, energy reserves and vertical location. There are two habitats, a surface habitat with seasonal predation risk and food availability, and a safer deep habitat with no food and low metabolism (diapause). Optimal life histories (diapause and energy allocation strategies) are found by dynamic programming. Seasonal egg fitness (reproductive value) emerges from the model and peak values are typically before the feeding season. Disentangling the fitness components, we conclude that seasonality in egg fitness is caused both by environmental seasonality in food and predation risk and by time-constraints on development and diapause preparation. Realised egg production, as predicted from population simulations, does not match the seasonal peak in offspring fitness but is delayed relative to peak egg fitness. We term this an 'internal life history mismatch' as constraints and tradeoffs cause sub-optimal birth dates for most eggs whereas mothers maximise their reproductive value by high number of eggs rather than few and optimally timed eggs. The earliest eggs have a disproportionately high contribution to population recruitment, emphasising the importance of early eggs and the need to understand seasonal patterns in offspring fitness.     

Relating traits to diversification: a simple test

We describe a simple comparative method for determining whether rates of diversification are correlated with continuous traits in species-level phylogenies. This involves comparing traits of species with net speciation rate (number of nodes linking extant species with the root divided by the root to tip evolutionary distance), using a phylogenetically corrected correlation. We use simulations to examine the power of this test. We find that the approach has acceptable power to uncover relationships between speciation and a continuous trait and is robust to background random extinction; however, the power of the approach is reduced when the rate of trait evolution is decreased. The test has low power to relate diversification to traits when extinction rate is correlated with the trait. Clearly, there are inherent limitations in using only data on extant species to infer correlates of extinction; however, this approach is potentially a powerful tool in analyzing correlates of speciation.     

Chancroid transmission dynamics: a mathematical modeling approach

Mathematical models have long been used to better understand disease transmission dynamics and how to effectively control them. Here, a chancroid infection model is presented and analyzed. The disease-free equilibrium is shown to be globally asymptotically stable when the reproduction number is less than unity. High levels of treatment are shown to reduce the reproduction number suggesting that treatment has the potential to control chancroid infections in any given community. This result is also supported by numerical simulations which show a decline in chancroid cases whenever the reproduction number is less than unity.     

Birth seasonality and calf mortality in a large population of Asian elephants

In seasonal environments, many species concentrate their reproduction in the time of year most likely to maximize offspring survival. Asian elephants (Elephas maximus) inhabit regions with seasonal climate, but females can still experience 16‐week reproductive cycles throughout the year. Whether female elephants nevertheless concentrate births on periods with maximum offspring survival prospects remains unknown. We investigated the seasonal timing of births, and effects of birth month on short‐ and long‐term mortality of Asian elephants, using a unique demographic data set of 2350 semicaptive, longitudinally monitored logging elephants from Myanmar experiencing seasonal variation in both workload and environmental conditions. Our results show variation in birth rate across the year, with 41% of births occurring between December and March. This corresponds to the cool, dry period and the beginning of the hot season, and to conceptions occurring during the resting, nonlogging period between February and June. Giving birth during the peak December to March period improves offspring survival, as the odds for survival between age 1 and 5 years are 44% higher for individuals born during the high birth rate period than those conceived during working months. Our results suggest that seasonal conditions, most likely maternal workload and/or climate, limit conception rate and calf survival in this population through effects on maternal stress, estrus cycles, or access to mates. This has implications for improving the birth rate and infant survival in captive populations by limiting workload of females of reproductive age. As working populations are currently unsustainable and supplemented through the capture of wild elephants, it is imperative to the conservation of Asian elephants to understand and alleviate the effects of seasonal conditions on vital rates in the working population in order to reduce the pressure for further capture from the wild.     

Oscillations in a size-structured prey-predator model

This article introduces a predator–prey model with the prey structured by body size, based on reports in the literature that predation rates are prey-size specific. The model is built on the foundation of the one-species physiologically structured models studied earlier. Three types of equilibria are found: extinction, multiple prey-only equilibria and possibly multiple predator–prey coexistence equilibria. The stabilities of the equilibria are investigated. Comparison is made with the underlying ODE Lotka–Volterra model. It turns out that the ODE model can exhibit sustain oscillations if there is an Allee effect in the net reproduction rate, that is the net reproduction rate grows for some range of the prey’s population size. In contrast, it is shown that the structured PDE model can exhibit sustain oscillations even if the net reproductive rate is strictly declining with prey population size. We find that predation, even size-non-specific linear predation can destabilize a stable prey-only equilibrium, if reproduction is size specific and limited to individuals of large enough size. Furthermore, we show that size-specific predation can also destabilize the predator–prey equilibrium in the PDE model. We surmise that size-specific predation allows for temporary prey escape which is responsible for destabilization in the predator–prey dynamics.