The Galilean turn in population ecology

The standard mathematical models in population ecology assume that a population's growth rate is a function of its environment. In this paper we investigate an alternative proposal according to which the rate of change of the growth rate is a function of the environment and of environmental change. We focus on the philosophical issues involved in such a fundamental shift in theoretical assumptions, as well as on the explanations the two theories offer for some of the key data such as cyclic populations. We also discuss the relationship between this move in population ecology and a similar move from first-order to second-order differential equations championed by Galileo and Newton in celestial mechanics.

How predation and landscape fragmentation affect vole population dynamics

Background: Microtine species in Fennoscandia display a distinct north-south gradient from regular cycles to stable populations. The gradient has often been attributed to changes in the interactions between microtines and their predators. Although the spatial structure of the environment is known to influence predator-prey dynamics of a wide range of species, it has scarcely been considered in relation to the Fennoscandian gradient. Furthermore, the length of microtine breeding season also displays a north-south gradient. However, little consideration has been given to its role in shaping or generating population cycles. Because these factors covary along the gradient it is difficult to distinguish their effects experimentally in the field. The distinction is here attempted using realistic agent-based modelling. Methodology/Principal Findings: By using a spatially explicit computer simulation model based on behavioural and ecological data from the field vole (Microtus agrestis), we generated a number of repeated time series of vole densities whose mean population size and amplitude were measured. Subsequently, these time series were subjected to statistical autoregressive modelling, to investigate the effects on vole population dynamics of making predators more specialised, of altering the breeding season, and increasing the level of habitat fragmentation. We found that fragmentation as well as the presence of specialist predators are necessary for the occurrence of population cycles. Habitat fragmentation and predator assembly jointly determined cycle length and amplitude. Length of vole breeding season had little impact on the oscillations. Significance: There is good agreement between our results and the experimental work from Fennoscandia, but our results allow distinction of causation that is hard to unravel in field experiments. We hope our results will help understand the reasons for cycle gradients observed in other areas. Our results clearly demonstrate the importance of landscape fragmentation for population cycling and we recommend that the degree of fragmentation be more fully considered in future analyses of vole dynamics.

The stochastic nature of predator-prey cycles

We study by numerical simulations the time correlation function of a stochastic lattice model describing the dynamics of coexistence of two interacting biological species that present time cycles in the number of species individuals. Its asymptotic behavior is shown to decrease in time as a sinusoidal exponential function from which we extract the dominant eigenvalue of the evolution operator related to the stochastic dynamics showing that it is complex with the imaginary part being the frequency of the population cycles. The transition from the oscillatory to the nonoscillatory behavior occurs when the asymptotic behavior of the time correlation function becomes a pure exponential, that is, when the real part of the complex eigenvalue equals a real eigenvalue. We also show that the amplitude of the undamped oscillations increases with the square root of the area of the habitat as ordinary random fluctuations. (C) 2009 Elsevier B.V. All rights reserved.

How population loss through habitat boundaries determines the dynamics of a predator-prey system

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Time correlation function in systems with two coexisting biological species

We study a stochastic lattice model describing the dynamics of coexistence of two interacting biological species. The model comprehends the local processes of birth, death, and diffusion of individuals of each species and is grounded on interaction of the predator-prey type. The species coexistence can be of two types: With self-sustained coupled time oscillations of population densities and without oscillations. We perform numerical simulations of the model on a square lattice and analyze the temporal behavior of each species by computing the time correlation functions as well as the spectral densities. This analysis provides an appropriate characterization of the different types of coexistence. It is also used to examine linked population cycles in nature and in experiment.

La toxoplasmose chez les Inuits : investigation de l'écologie de toxoplasma gondii dans l'Arctique canadien

Toxoplasma gondii, un protozoaire très répandu dans le monde, peut infecter de nombreuses espèces homéothermes incluant les mammifères et les oiseaux qui développent alors une toxoplasmose. L’impact de la toxoplasmose en termes de santé publique est majeur, particulièrement chez les personnes immunodéprimées et les foetus. Les niveaux d’infection humaine dans certaines régions de l’Arctique Canadien sont parmi les plus élevés au monde et ce, malgré l’absence de félidés qui sont les seuls hôtes capables d’excréter T. gondii. Plusieurs études ont suggéré la consommation de viande crue de mammifères marins et notamment de phoques comme source d’infection des Inuits. Notre travail de recherche visait à comprendre les mécanismes de dispersion de T. gondii dans les écosystèmes aquatiques menant à la contamination du milieu marin de l’Arctique par des oocystes, et à évaluer l’importance de cette voie de dispersion dans l’infection des phoques et conséquemment dans celle des Inuits. Notre hypothèse était que les oocystes de T. gondii, excrétés durant l’hiver par des félidés dans le Subarctique et transportés par les rivières pendant la fonte printanière, contaminaient les estuaires de l’Arctique Canadien. Dans un premier temps, une étude transversale de séroprévalence chez les phoques de l’Arctique Canadien a montré que ces populations étaient infectées par T. gondii et pouvaient ainsi a priori constituer une source d’infection pour les Inuit. Des variations spatio-temporelles de la séroprévalence étaient observées suggérant un lien potentiel avec des variations dans la contamination environnementale par les oocystes. Un schéma conceptuel explicitant les mécanismes de transport et de devenir des oocystes de T. gondii, du phénomène de la fonte de la neige jusqu’à l’exposition des organismes marins, a été proposé dans le chapitre suivant. Des interactions entre les différents mécanismes identifiés, qui agissent sur des échelles spatio-temporelles variées, devraient favoriser l’apparition de concentrations relativement élevées aux estuaires permettant ainsi l’exposition et potentiellement l’infection de phoques. Pour évaluer la contamination environnementale par les oocystes excrétés par la population de lynx du bassin versant de l’Arctique Canadien (les seuls félidés majoritairement distribués dans ce vaste territoire), nous avons mené une étude sérologique de type transversale dans cette population. Cette étude a permis de montrer que des lynx étaient infectés par T. gondii et a également suggéré que la dynamique des cycles de populations lynx-lièvres pouvait être un processus important dans la transmission de T. gondii. Finalement, la modélisation du transport hydrique des oocystes a indiqué que les concentrations hypothétiques d’oocystes dans l’eau de la fonte pourraient être suffisantes pour permettre l’exposition au niveau des estuaires de bivalves filtreurs, qui sont des proies pour les phoques et donc potentiellement des sources infectieuses pour ces derniers. Dans des écosystèmes nordiques en pleine mutation, la compréhension des mécanismes de transmission d’agents pathogènes d’origine hydrique comme T. gondii est plus que nécessaire, notamment dans le but de protéger les populations fragilisées de ces régions.

Stochastic lattice gas model describing the dynamics of the SIRS epidemic process

We study a stochastic process describing the onset of spreading dynamics of an epidemic in a population composed of individuals of three classes: susceptible (S), infected (I), and recovered (R). The stochastic process is defined by local rules and involves the following cyclic process: S -> I -> R -> S (SIRS). The open process S -> I -> R (SIR) is studied as a particular case of the SIRS process. The epidemic process is analyzed at different levels of description: by a stochastic lattice gas model and by a birth and death process. By means of Monte Carlo simulations and dynamical mean-field approximations we show that the SIRS stochastic lattice gas model exhibit a line of critical points separating the two phases: an absorbing phase where the lattice is completely full of S individuals and an active phase where S, I and R individuals coexist, which may or may not present population cycles. The critical line, that corresponds to the onset of epidemic spreading, is shown to belong in the directed percolation universality class. By considering the birth and death process we analyze the role of noise in stabilizing the oscillations. (C) 2009 Elsevier B.V. All rights reserved.

(Table 2) Carbon and nitrogen isotopic ratios in blood sampled from arctic foxes (Vulpes lagopus) on Bylot Island between 2003-2008

Inter-individual variation in diet within generalist animal populations is thought to be a widespread phenomenon but its potential causes are poorly known. Inter-individual variation can be amplified by the availability and use of allochthonous resources, i.e., resources coming from spatially distinct ecosystems. Using a wild population of arctic fox as a study model, we tested hypotheses that could explain variation in both population and individual isotopic niches, used here as proxy for the trophic niche. The arctic fox is an opportunistic forager, dwelling in terrestrial and marine environments characterized by strong spatial (arctic-nesting birds) and temporal (cyclic lemmings) fluctuations in resource abundance. First, we tested the hypothesis that generalist foraging habits, in association with temporal variation in prey accessibility, should induce temporal changes in isotopic niche width and diet. Second, we investigated whether within-population variation in the isotopic niche could be explained by individual characteristics (sex and breeding status) and environmental factors (spatiotemporal variation in prey availability). We addressed these questions using isotopic analysis and Bayesian mixing models in conjunction with linear mixed-effects models. We found that: i) arctic fox populations can simultaneously undergo short-term (i.e., within a few months) reduction in both isotopic niche width and inter-individual variability in isotopic ratios, ii) individual isotopic ratios were higher and more representative of a marine-based diet for non-breeding than breeding foxes early in spring, and iii) lemming population cycles did not appear to directly influence the diet of individual foxes after taking their breeding status into account. However, lemming abundance was correlated to proportion of breeding foxes, and could thus indirectly affect the diet at the population scale.

(Table 1) Nesting characteristics of greater snow geese (Chen caerulescens atlanticus) on Bylot Island between 1995-2005

1. Habitat heterogeneity and predator behaviour can strongly affect predator-prey interactions but these factors are rarely considered simultaneously, especially when systems encompass multiple predators and prey. 2. In the Arctic, greater snow geese Anser caerulescens atlanticus L. nest in two structurally different habitats: wetlands that form intricate networks of water channels, and mesic tundra where such obstacles are absent. In this heterogeneous environment, goose eggs are exposed to two types of predators: the arctic fox Vulpes lagopus L. and a diversity of avian predators. We hypothesized that, contrary to birds, the hunting ability of foxes would be impaired by the structurally complex wetland habitat, resulting in a lower predation risk for goose eggs. 3. In addition, lemmings, the main prey of foxes, show strong population cycles. We thus further examined how their fluctuations influenced the interaction between habitat heterogeneity and fox predation on goose eggs. 4. An experimental approach with artificial nests suggested that foxes were faster than avian predators to find unattended goose nests in mesic tundra whereas the reverse was true in wetlands. Foxes spent 3-5 times more time between consecutive attacks on real goose nests in wetlands than in mesic tundra. Their attacks on goose nests were also half as successful in wetlands than in mesic tundra whereas no difference was found for avian predators. 5. Nesting success in wetlands (65%) was higher than in mesic tundra (56%) but the difference between habitats increased during lemming crashes (15%) compared to other phases of the cycle (5%). Nests located at the edge of wetland patches were also less successful than central ones, suggesting a gradient in accessibility of goose nests in wetlands for foxes. 6. Our study shows that the structural complexity of wetlands decreases predation risk from foxes but not avian predators in arctic-nesting birds. Our results also demonstrate that cyclic lemming populations indirectly alter the spatial distribution of productive nests due to a complex interaction between habitat structure, prey-switching and foraging success of foxes.

Monarchs across the Pacific: the Columbus hypothesis revisited

The 'Columbus hypothesis' suggests that the annual north-south return migration of Danaus plexippus in North America is a very recently evolved behaviour, less than 200 years old. This hypothesis rests, in part, on an analysis of the 19th century spread of the monarch across the Pacific that assumes a continuous east to west movement and is based predominantly on one publication. We review all the contemporary literature and present new analysis of the data. The movement of the monarch across the Pacific in the second half of the 19th century is best explained by a model which involves no more than three spot introductions, directly or indirectly aided by human movement, followed by natural spread of the monarch across island groups. Contemporary records refer to 'boom' and 'bust' population cycles on newly settled islands, which may have led to high rates of monarch movement. We see no evidence in the records to suggest an east to west sweep by monarch populations as suggested by the Columbus hypothesis. (C) 2004 The Linnean Society of London.

Ecological drivers of a vector borne pathogen: distribution and abundance of Borrelia burgdorferi sensu lato and its vector Ixodes ricinus in Scotland

Vector-borne disease emergence in recent decades has been associated with different environmental drivers including changes in habitat, hosts and climate. Lyme borreliosis is among the most important vector-borne diseases in the Northern hemisphere and is an emerging disease in Scotland. Transmitted by Ixodid tick vectors between large numbers of wild vertebrate host species, Lyme borreliosis is caused by bacteria from the Borrelia burgdorferi sensu lato species group. Ecological studies can inform how environmental factors such as host abundance and community composition, habitat and landscape heterogeneity contribute to spatial and temporal variation in risk from B. burgdorferi s.l. In this thesis a range of approaches were used to investigate the effects of vertebrate host communities and individual host species as drivers of B. burgdorferi s.l. dynamics and its tick vector Ixodes ricinus. Host species differ in reservoir competence for B. burgdorferi s.l. and as hosts for ticks. Deer are incompetent transmission hosts for B. burgdorferi s.l. but are significant hosts of all life-stages of I. ricinus. Rodents and birds are important transmission hosts of B. burgdorferi s.l. and common hosts of immature life-stages of I. ricinus. In this thesis, surveys of woodland sites revealed variable effects of deer density on B. burgdorferi prevalence, from no effect (Chapter 2) to a possible ‘dilution’ effect resulting in lower prevalence at higher deer densities (Chapter 3). An invasive species in Scotland, the grey squirrel (Sciurus carolinensis), was found to host diverse genotypes of B. burgdorferi s.l. and may act as a spill-over host for strains maintained by native host species (Chapter 4). Habitat fragmentation may alter the dynamics of B. burgdorferi s.l. via effects on the host community and host movements. In this thesis, there was lack of persistence of the rodent associated genospecies of B. burgdorferi s.l. within a naturally fragmented landscape (Chapter 3). Rodent host biology, particularly population cycles and dispersal ability are likely to affect pathogen persistence and recolonization in fragmented habitats. Heterogeneity in disease dynamics can occur spatially and temporally due to differences in the host community, habitat and climatic factors. Higher numbers of I. ricinus nymphs, and a higher probability of detecting a nymph infected with B. burgdorferi s.l., were found in areas with warmer climates estimated by growing degree days (Chapter 2). The ground vegetation type associated with the highest number of I. ricinus nymphs varied between studies in this thesis (Chapter 2 & 3) and does not appear to be a reliable predictor across large areas. B. burgdorferi s.l. prevalence and genospecies composition was highly variable for the same sites sampled in subsequent years (Chapter 2). This suggests that dynamic variables such as reservoir host densities and deer should be measured as well as more static habitat and climatic factors to understand the drivers of B. burgdorferi s.l. infection in ticks. Heterogeneity in parasite loads amongst hosts is a common finding which has implications for disease ecology and management. Using a 17-year data set for tick infestations in a wild bird community in Scotland, different effects of age and sex on tick burdens were found among four species of passerine bird (Chapter 5). There were also different rates of decline in tick burdens among bird species in response to a long term decrease in questing tick pressure over the study. Species specific patterns may be driven by differences in behaviour and immunity and highlight the importance of comparative approaches. Combining whole genome sequencing (WGS) and population genetics approaches offers a novel approach to identify ecological drivers of pathogen populations. An initial analysis of WGS from B. burgdorferi s.s. isolates sampled 16 years apart suggests that there is a signal of measurable evolution (Chapter 6). This suggests demographic analyses may be applied to understand ecological and evolutionary processes of these bacteria. This work shows how host communities, habitat and climatic factors can affect the local transmission dynamics of B. burgdorferi s.l. and the potential risk of infection to humans. Spatial and temporal heterogeneity in pathogen dynamics poses challenges for the prediction of risk. New tools such as WGS of the pathogen (Chapter 6) and blood meal analysis techniques will add power to future studies on the ecology and evolution of B. burgdorferi s.l.

Role of noise in population dynamics cycles

Noise is an intrinsic feature of population dynamics and plays a crucial role in oscillations called phase-forgetting quasicycles by converting damped into sustained oscillations. This function of noise becomes evident when considering Langevin equations whose deterministic part yields only damped oscillations. We formulate here a consistent and systematic approach to population dynamics, leading to a Fokker-Planck equation and the associate Langevin equations in accordance with this conceptual framework, founded on stochastic lattice-gas models that describe spatially structured predator-prey systems. Langevin equations in the population densities and predator-prey pair density are derived in two stages. First, a birth-and-death stochastic process in the space of prey and predator numbers and predator-prey pair number is obtained by a contraction method that reduces the degrees of freedom. Second, a van Kampen expansion in the inverse of system size is then performed to get the Fokker-Planck equation. We also study the time correlation function, the asymptotic behavior of which is used to characterize the transition from the cyclic coexistence of species to the ordinary coexistence.

Influence of the tidal and rainfall cycles on the population structure and density of Mesacanthion hirsutum gerlach (nematoda, thoracostomopsidae) on a tropical sandy beach (Tamandaré Bay, Pernambuco, Brazil)

A estrutura populacional (machos, fêmeas e juvenis) e densidade de Mesacanthion hirsutum Gerlach foram estudadas durante as marés (baixa, enchente, alta e vazante) de dois ciclos de maré consecutivos, em quatro meses diferentes do ano (Maio, Julho, Setembro e Novembro). As variações de densidade de Mesacanthion hirsutum mostraram associação com o ciclo de chuvas, com densidades mais baixas durante Julho e Setembro e significativamente maiores em Maio e Novembro. A estrutura populacional constituiu-se em sua maior parte por juvenis indicando uma reprodução continua durante todo o período de estudo. Não foram encontradas diferenças significativas entre os períodos claros e escuros do dia, contudo maiores densidades foram detectadas durante as marés altas e vazantes demonstrando que a espécie pode estar se dispersando através da coluna d'água e/ou migrando dentro do sedimento.

Cumulative live birth rates after one or more complete cycles of IVF : a population-based study of linked cycle data from 178 898 women

© The Author 2016. Published by Oxford University Press on behalf of the European Society of Human Reproduction and Embryology. All rights reserved. For Permissions, please email: journals.permissions@oup.com.

Circulation of antibodies against yellow fever virus in a simian population in the area of Porto Primavera Hydroelectric Plant, São Paulo, Brazil

Yellow fever (YF) is an acute viral infectious disease transmitted by mosquitoes which occurs in two distinct epidemiological cycles: sylvatic and urban. In the sylvatic cycle, the virus is maintained by monkey's infection and transovarian transmission in vectors. Surveillance of non-human primates is required for the detection of viral circulation during epizootics, and for the identification of unaffected or transition areas. An ELISA (enzyme-linked immunosorbent assay) was standardized for estimation of the prevalence of IgG antibodies against yellow fever virus in monkey sera (Alouatta caraya) from the reservoir area of Porto Primavera Hydroelectric Plant, in the state of São Paulo, Brazil. A total of 570 monkey sera samples were tested and none was reactive to antibodies against yellow fever virus. The results corroborate the epidemiology of yellow fever in the area. Even though it is considered a transition area, there were no reports to date of epizootics or yellow fever outbreaks in humans. Also, entomological investigations did not detect the presence of vectors of this arbovirus infection. ELISA proved to be fast, sensitive, an adequate assay, and an instrument for active search in the epidemiological surveillance of yellow fever allowing the implementation of prevention actions, even before the occurrence of epizootics.