Genetic analysis of complex demographic scenarios: Spatially expanding populations of the cane toad, Bufo marinus

Inferring the spatial expansion dynamics of invading species from molecular data is notoriously difficult due to the complexity of the processes involved. For these demographic scenarios, genetic data obtained from highly variable markers may be profitably combined with specific sampling schemes and information from other sources using a Bayesian approach. The geographic range of the introduced toad Bufo marinus is still expanding in eastern and northern Australia, in each case from isolates established around 1960. A large amount of demographic and historical information is available on both expansion areas. In each area, samples were collected along a transect representing populations of different ages and genotyped at 10 microsatellite loci. Five demographic models of expansion, differing in the dispersal pattern for migrants and founders and in the number of founders, were considered. Because the demographic history is complex, we used an approximate Bayesian method, based on a rejection-regression algorithm. to formally test the relative likelihoods of the five models of expansion and to infer demographic parameters. A stepwise migration-foundation model with founder events was statistically better supported than other four models in both expansion areas. Posterior distributions supported different dynamics of expansion in the studied areas. Populations in the eastern expansion area have a lower stable effective population size and have been founded by a smaller number of individuals than those in the northern expansion area. Once demographically stabilized, populations exchange a substantial number of effective migrants per generation in both expansion areas, and such exchanges are larger in northern than in eastern Australia. The effective number of migrants appears to be considerably lower than that of founders in both expansion areas. We found our inferences to be relatively robust to various assumptions on marker. demographic, and historical features. The method presented here is the only robust, model-based method available so far, which allows inferring complex population dynamics over a short time scale. It also provides the basis for investigating the interplay between population dynamics, drift, and selection in invasive species.

Resource availability and the persistence of seed-eating bird populations in agricultural landscapes – a mechanistic modelling approach

1. Reductions in resource availability, associated with land-use change and agricultural intensification in the UK and Europe, have been linked with the widespread decline of many farmland bird species over recent decades. However, the underlying ecological processes which link resource availability and population trends are poorly understood. 2. We construct a spatial depletion model to investigate the relationship between the population persistence of granivorous birds within the agricultural landscape and the temporal dynamics of stubble field availability, an important source of winter food for many of those species. 3. The model is capable of accurately predicting the distribution of a given number of finches and buntings amongst patches of different stubble types in an agricultural landscape over the course of a winter and assessing the relative value of different landscapes in terms of resource availability. 4. Sensitivity analyses showed that the model is relatively robust to estimates of energetic requirements, search efficiency and handling time but that daily seed survival estimates have a strong influence on model fit. Understanding resource dynamics in agricultural landscapes is highlighted as a key area for further research. 5. There was a positive relationship between the predicted number of bird days supported by a landscape over-winter and the breeding population trend for yellowhammer Emberiza citrinella, a species for which survival has been identified as the primary driver of population dynamics, but not for linnet Carduelis cannabina, a species for which productivity has been identified as the primary driver of population dynamics. 6. Synthesis and applications. We believe this model can be used to guide the effective delivery of over-winter food resources under agri-environment schemes and to assess the impacts on granivorous birds of changing resource availability associated with novel changes in land use. This could be very important in the future as farming adapts to an increasingly dynamic trading environment, in which demands for increased agricultural production must be reconciled with objectives for environmental protection, including biodiversity conservation.

In vitro examination of the effect of Orlistat on the ability of the faecal microbiota to utilize dietary lipids

Orlistat is an anti-obesity treatment with which several gastrointestinal (GI) side-effects are commonly associated in the initial stages of therapy. There is no physiological explanation as to why two-thirds of those who take the drug experience one or more side-effects. It has been hypothesized that the GI microbiota may protect from or contribute to these GI disturbances. Using in vitro batch culture and human gut model systems, studies were conducted to determine whether increased availability of dietary lipids and/or orlistat affect the composition and/or activity of the faecal microbiota. Results from 24-h batch culture fermentation experiments demonstrated no effect of orlistat in the presence or absence of a dietary lipid (olive oil) on the composition of bacterial communities [as determined by fluorescence in situ hybridization (FISH) and denaturing gradient gel electrophoresis (DGGE) analyses], but did show there was great variability in the lipolytic activities of the microbiotas of individuals, as determined by gas chromatography analysis of long-chain fatty acids in samples. Subsequent studies focused on the effect of orlistat in the presence and absence of lipid in in vitro human gut model systems. Systems were run for 14 days with gut model medium (GMM) only (to steady state, SS), then fed at 12-h intervals with 50 mg orlistat, 2 g olive oil or a mixture of both for 14 days. FISH and DGGE were used to monitor changes in bacterial populations. Bacteria were cultivated from the GMM only (control) systems at SS. All strains isolated were screened for lipolytic activity using tributyrin agar. FISH and DGGE demonstrated that none of the compounds (singly or in combination) added to the systems had any notable effect on microbial population dynamics for any of the donors, although Subdoligranulum populations appeared to be inhibited by orlistat in the presence or absence of lipid. Orlistat had little or no effect on the metabolism of indigenous and added lipids in the fermentation systems, but there was great variability in the way the faecal microbiotas of the donors were able to degrade added lipids. Variability in lipid degradation could be correlated with the number and activity of isolated lipolytic bacteria. The mechanism by which orlistat and the GI microbiota cause side-effects in individuals is unknown, but several hypotheses have been proposed to account for their manifestation. The demonstration of great variability in the lipolytic activity of microbiotas to degrade lipids led to a large-scale cultivation-based study of lipolytic/lipase-positive bacteria present in the human faecal microbiota. Of 4,000 colonies isolated from 15 donors using five different agars, 378 strains were identified that had lipase activity. Molecular identification of strains isolated from five donors demonstrated that lipase activity is more prevalent in the human GI microbiota than previously thought, with members of the phyla Firmicutes, Bacteroidetes and Actinobacteria identified. Molecular identification and characterization of the substrate specificities of the strains will be carried out as part of ongoing work.

Sensitivity analysis of a cyanobacterial growth and movement model under two different flow regimes

Bloom-forming and toxin-producing cyanobacteria remain a persistent nuisance across the world. Modelling cyanobacterial behaviour in freshwaters is an important tool for understanding their population dynamics and predicting the location and timing of the bloom events in lakes, reservoirs and rivers. A new deterministic–mathematical model was developed, which simulates the growth and movement of cyanobacterial blooms in river systems. The model focuses on the mathematical description of the bloom formation, vertical migration and lateral transport of colonies within river environments by taking into account the major factors that affect the cyanobacterial bloom formation in rivers including light, nutrients and temperature. A parameter sensitivity analysis using a one-at-a-time approach was carried out. There were two objectives of the sensitivity analysis presented in this paper: to identify the key parameters controlling the growth and movement patterns of cyanobacteria and to provide a means for model validation. The result of the analysis suggested that maximum growth rate and day length period were the most significant parameters in determining the population growth and colony depth, respectively.

The effects of landscape modifications on the long-term persistence of animal populations

Background: The effects of landscape modifications on the long-term persistence of wild animal populations is of crucial importance to wildlife managers and conservation biologists, but obtaining experimental evidence using real landscapes is usually impossible. To circumvent this problem we used individual-based models (IBMs) of interacting animals in experimental modifications of a real Danish landscape. The models incorporate as much as possible of the behaviour and ecology of four species with contrasting life-history characteristics: skylark (Alauda arvensis), vole (Microtus agrestis), a ground beetle (Bembidion lampros) and a linyphiid spider (Erigone atra). This allows us to quantify the population implications of experimental modifications of landscape configuration and composition. Methodology/Principal Findings: Starting with a real agricultural landscape, we progressively reduced landscape complexity by (i) homogenizing habitat patch shapes, (ii) randomizing the locations of the patches, and (iii) randomizing the size of the patches. The first two steps increased landscape fragmentation. We assessed the effects of these manipulations on the long-term persistence of animal populations by measuring equilibrium population sizes and time to recovery after disturbance. Patch rearrangement and the presence of corridors had a large effect on the population dynamics of species whose local success depends on the surrounding terrain. Landscape modifications that reduced population sizes increased recovery times in the short-dispersing species, making small populations vulnerable to increasing disturbance. The species that were most strongly affected by large disturbances fluctuated little in population sizes in years when no perturbations took place. Significance: Traditional approaches to the management and conservation of populations use either classical methods of population analysis, which fail to adequately account for the spatial configurations of landscapes, or landscape ecology, which accounts for landscape structure but has difficulty predicting the dynamics of populations living in them. Here we show how realistic and replicable individual-based models can bridge the gap between non-spatial population theory and non-dynamic landscape ecology. A major strength of the approach is its ability to identify population vulnerabilities not detected by standard population viability analyses.

Spatial patterns of office employment in the New York region

This study analyzes the regional spatial dynamics of the New York region for a period of roughly twenty years and places the effects of the 9/11 terrorist attacks in the context of longer-term regional dynamics. The analysis reveals that office-using industries are still heavily concentrated in Manhattan despite ongoing decentralization in many of these industries over the last twenty years. Financial services tend to be highly concentrated in Manhattan whereas administrative and support services are the least concentrated of the six major office-using industry groups. Although office employment has been by and large stagnant in Manhattan for at least two decades, growth of output per worker has outpaced the CMSA as well as the national average. This productivity differential is mainly attributable to competitive advantages of office-using industries in Manhattan and not to differences in industry composition. Finally, the zip-code level analysis of the Manhattan core area yielded further evidence of the existence of significant spillover effects at the small-scale level.

Sources of variation in the ampicillin-resistant Escherichia coli concentration in the feces of organic broiler chickens

Currently, there are limited published data for the population dynamics of antimicrobial-resistant commensal bacteria. This study was designed to evaluate both the proportions of the Escherichia coli populations that are resistant to ampicillin at the level of the individual chicken on commercial broiler farms and the feasibility of obtaining repeated measures of fecal E. coli concentrations. Short-term temporal variation in the concentration of fecal E. coli was investigated, and a preliminary assessment was made of potential factors involved in the shedding of high numbers of ampicillin-resistant E. coli by growing birds in the absence of the use of antimicrobial drugs. Multilevel linear regression modeling revealed that the largest component of random variation in log-transformed fecal E. coli concentrations was seen between sampling occasions for individual birds. The incorporation of fixed effects into the model demonstrated that the older, heavier birds in the study were significantly more likely (P = 0.0003) to shed higher numbers of ampicillin-resistant E. coli. This association between increasing weight and high shedding was not seen for the total fecal E. coli population (P = 0.71). This implies that, in the absence of the administration of antimicrobial drugs, the proportion of fecal E. coli that was resistant to ampicillin increased as the birds grew. This study has shown that it is possible to collect quantitative microbiological data on broiler farms and that such data could make valuable contributions to risk assessments concerning the transfer of resistant bacteria between animal and human populations.

Representing the acquisition and use of energy by individuals in agent-based models of animal populations

Summary 1. Agent-based models (ABMs) are widely used to predict how populations respond to changing environments. As the availability of food varies in space and time, individuals should have their own energy budgets, but there is no consensus as to how these should be modelled. Here, we use knowledge of physiological ecology to identify major issues confronting the modeller and to make recommendations about how energy budgets for use in ABMs should be constructed. 2. Our proposal is that modelled animals forage as necessary to supply their energy needs for maintenance, growth and reproduction. If there is sufficient energy intake, an animal allocates the energy obtained in the order: maintenance, growth, reproduction, energy storage, until its energy stores reach an optimal level. If there is a shortfall, the priorities for maintenance and growth/reproduction remain the same until reserves fall to a critical threshold below which all are allocated to maintenance. Rates of ingestion and allocation depend on body mass and temperature. We make suggestions for how each of these processes should be modelled mathematically. 3. Mortality rates vary with body mass and temperature according to known relationships, and these can be used to obtain estimates of background mortality rate. 4. If parameter values cannot be obtained directly, then values may provisionally be obtained by parameter borrowing, pattern-oriented modelling, artificial evolution or from allometric equations. 5. The development of ABMs incorporating individual energy budgets is essential for realistic modelling of populations affected by food availability. Such ABMs are already being used to guide conservation planning of nature reserves and shell fisheries, to assess environmental impacts of building proposals including wind farms and highways and to assess the effects on nontarget organisms of chemicals for the control of agricultural pests. Keywords: bioenergetics; energy budget; individual-based models; population dynamics.

Fertilisers and insect herbivores: a meta-analysis

Despite the importance of a thorough understanding of the effect of synthetic fertiliser on insect population dynamics, existing literature is conflicting and an area of intense debate. Here, a categorical random-effects meta-analysis and a vote count meta-analysis are employed to examine the effects of nitrogen(N), phosphorus (P), potassium (K) and NPK fertiliser on insect population dynamics. In agreement with the general consensus, insects were found to respond positively, overall, to fertilisers. Sucking insects showed a much stronger response to fertilisers than chewing insects. The environment in which a study is conducted can have a marked effect on insect responses to fertiliser, with natural environments showing the potential for buffering effects of nitrogen fertilisers in particular. As well as highlighting the potential shortfall in the amount of research investigating particularly the effects of potassium and phosphorus, this study provides an invaluable flag post in the ongoing research investigating fertiliser effects on ecosystems.

A dynamic causal model of the coupling between pulse stimulation and neural activity

We present a dynamic causal model that can explain context-dependent changes in neural responses, in the rat barrel cortex, to an electrical whisker stimulation at different frequencies. Neural responses were measured in terms of local field potentials. These were converted into current source density (CSD) data, and the time series of the CSD sink was extracted to provide a time series response train. The model structure consists of three layers (approximating the responses from the brain stem to the thalamus and then the barrel cortex), and the latter two layers contain nonlinearly coupled modules of linear second-order dynamic systems. The interaction of these modules forms a nonlinear regulatory system that determines the temporal structure of the neural response amplitude for the thalamic and cortical layers. The model is based on the measured population dynamics of neurons rather than the dynamics of a single neuron and was evaluated against CSD data from experiments with varying stimulation frequency (1–40 Hz), random pulse trains, and awake and anesthetized animals. The model parameters obtained by optimization for different physiological conditions (anesthetized or awake) were significantly different. Following Friston, Mechelli, Turner, and Price (2000), this work is part of a formal mathematical system currently being developed (Zheng et al., 2005) that links stimulation to the blood oxygen level dependent (BOLD) functional magnetic resonance imaging (fMRI) signal through neural activity and hemodynamic variables. The importance of the model described here is that it can be used to invert the hemodynamic measurements of changes in blood flow to estimate the underlying neural activity.

Within-host evolution of Staphylococcus aureus during asymptomatic carriage

Background Staphylococcus aureus is a major cause of healthcare associated mortality, but like many important bacterial pathogens, it is a common constituent of the normal human body flora. Around a third of healthy adults are carriers. Recent evidence suggests that evolution of S. aureus during nasal carriage may be associated with progression to invasive disease. However, a more detailed understanding of within-host evolution under natural conditions is required to appreciate the evolutionary and mechanistic reasons why commensal bacteria such as S. aureus cause disease. Therefore we examined in detail the evolutionary dynamics of normal, asymptomatic carriage. Sequencing a total of 131 genomes across 13 singly colonized hosts using the Illumina platform, we investigated diversity, selection, population dynamics and transmission during the short-term evolution of S. aureus. Principal Findings We characterized the processes by which the raw material for evolution is generated: micro-mutation (point mutation and small insertions/deletions), macro-mutation (large insertions/deletions) and the loss or acquisition of mobile elements (plasmids and bacteriophages). Through an analysis of synonymous, non-synonymous and intergenic mutations we discovered a fitness landscape dominated by purifying selection, with rare examples of adaptive change in genes encoding surface-anchored proteins and an enterotoxin. We found evidence for dramatic, hundred-fold fluctuations in the size of the within-host population over time, which we related to the cycle of colonization and clearance. Using a newly-developed population genetics approach to detect recent transmission among hosts, we revealed evidence for recent transmission between some of our subjects, including a husband and wife both carrying populations of methicillin-resistant S. aureus (MRSA). Significance This investigation begins to paint a picture of the within-host evolution of an important bacterial pathogen during its prevailing natural state, asymptomatic carriage. These results also have wider significance as a benchmark for future systematic studies of evolution during invasive S. aureus disease.

British container breeding mosquitoes: the impact of urbanisation and climate change on community composition and phenology

The proliferation of artificial container habitats in urban areas has benefitted urban adaptable mosquito species globally. In areas where mosquitoes transmit viruses and parasites, it can promote vector population productivity and fuel mosquito-borne disease outbreaks. In Britain, storage of water in garden water butts is increasing, potentially expanding mosquito larval habitats and influencing population dynamics and mosquito-human contact. Here we show that the community composition, abundance and phenology of mosquitoes breeding in experimental water butt containers were influenced by urbanisation. Mosquitoes in urban containers were less species-rich but present in significantly higher densities (100.4±21.3) per container than those in rural containers (77.7±15.1). Urban containers were dominated by Culex pipiens (a potential vector of West Nile Virus [WNV]) and appear to be increasingly exploited by Anopheles plumbeus (a human-biting potential WNV and malaria vector). Culex phenology was influenced by urban land use type, with peaks in larval abundances occurring earlier in urban than rural containers. Among other factors, this was associated with an urban heat island effect which raised urban air and water temperatures by 0.9°C and 1.2°C respectively. Further increases in domestic water storage, particularly in urban areas, in combination with climate changes will likely alter mosquito population dynamics in the UK.

The stability of ecosystems: a brief overview of the paradox of enrichment

In theory, enrichment of resource in a predator-prey model leads to destabilization of the system, thereby collapsing the trophic interaction, a phenomenon referred to as "the paradox of enrichment". After it was first proposed by Rosenzweig (1971), a number of subsequent studies were carried out on this dilemma over many decades. In this article, we review these theoretical and experimental works and give a brief overview of the proposed solutions to the paradox. The mechanisms that have been discussed are modifications of simple predator-prey models in the presence of prey that is inedible, invulnerable, unpalatable and toxic. Another class of mechanisms includes an incorporation of a ratio-dependent functional form, inducible defence of prey and density-dependent mortality of the predator. Moreover, we find a third set of explanations based on complex population dynamics including chaos in space and time. We conclude that, although any one of the various mechanisms proposed so far might potentially prevent destabilization of the predator-prey dynamics following enrichment, in nature different mechanisms may combine to cause stability, even when a system is enriched. The exact mechanisms, which may differ among systems, need to be disentangled through extensive field studies and laboratory experiments coupled with realistic theoretical models.

Agriculture and the life histories of Mauritius kestrels

Habitat modification for agriculture is one of the greatest current threats to global biodiversity. Studies show large-scale population declines and short-term demographic impacts, but knowledge of the long-term effects of agriculture on individuals remains poor. This thesis examines the short- and long-term impact of agriculture on a reintroduced population of the Mauritius kestrel Falco punctatus, a tropical forest-dwelling raptor endemic to the island of Mauritius, that also utilises agricultural habitats. This population is a particularly appropriate model system, because complete life history data exists for individuals over a 22-year period, alongside detailed habitat and climate data. Agriculture has a short-term detrimental effect on Mauritius kestrel breeding success by exacerbating the seasonal decline in fledgling production. This is partly driven by the habitat-specific composition of the prey community that kestrels exploit to feed their chicks. The fledglings from agricultural territories tend to recruit in agricultural territories. This is largely due to poor natal dispersal and fine-scale spatial autocorrelation in the habitat matrix. Breeders do not respond to agriculture in the breeding territory by dispersing, unless the pair bond is broken. Therefore, individuals originating in agricultural territories tend to recruit, and remain in, agricultural territories throughout their lives. In addition to this, females from agricultural natal territories have shorter lifespans, schedule their peak reproductive output earlier in life, and exhibit more rapid senescence than non-agricultural females. The combination of this long-term effect and the adult experience of agriculture imposed by life history and environmental constraints, leads to a lower mean lifetime reproductive rate compared to females originating in non-agricultural habitats. These results demonstrate that agriculture experienced in early life has a lifelong effect on individuals. The effects can persist in time and space, with potentially delayed effects on population dynamics. These findings are important for understanding species’ responses to agricultural expansion.

Sublethal neonicotinoid insecticide exposure reduces solitary bee reproductive success

1. Pollinating insects provide crucial and economically important ecosystem services to crops and wild plants, but pollinators, particularly bees, are globally declining as a result of various driving factors, including the prevalent use of pesticides for crop protection. Sublethal pesticide exposure negatively impacts numerous pollinator lifehistory traits, but its influence on reproductive success remains largely unknown. Such information is pivotal, however, to our understanding of the long-term effects on population dynamics. 2 We investigated the influence of field-realistic trace residues of the routinely used neonicotinoid insecticides thiamethoxam and clothianidin in nectar substitutes on the entire life-time fitness performance of the red mason bee Osmia bicornis. 3 We show that chronic, dietary neonicotinoid exposure has severe detrimental effects on solitary bee reproductive output. Neonicotinoids did not affect adult bee mortality; however, monitoring of fully controlled experimental populations revealed that sublethal exposure resulted in almost 50% reduced total offspring production and a significantly male-biased offspring sex ratio. 4 Our data add to the accumulating evidence indicating that sublethal neonicotinoid effects on non-Apis pollinators are expressed most strongly in a rather complex, fitness-related context. Consequently, to fully mitigate long-term impacts on pollinator population dynamics, present pesticide risk assessments need to be expanded to include whole life-cycle fitness estimates, as demonstrated in the present study using O. bicornis as a model.