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.

The automation and evaluation of nested clade phylogeographic analysis

Nested clade phylogeographic analysis (NCPA) is a popular method for reconstructing the demographic history of spatially distributed populations from genetic data. Although some parts of the analysis are automated, there is no unique and widely followed algorithm for doing this in its entirety, beginning with the data, and ending with the inferences drawn from the data. This article describes a method that automates NCPA, thereby providing a framework for replicating analyses in an objective way. To do so, a number of decisions need to be made so that the automated implementation is representative of previous analyses. We review how the NCPA procedure has evolved since its inception and conclude that there is scope for some variability in the manual application of NCPA. We apply the automated software to three published datasets previously analyzed manually and replicate many details of the manual analyses, suggesting that the current algorithm is representative of how a typical user will perform NCPA. We simulate a large number of replicate datasets for geographically distributed, but entirely random-mating, populations. These are then analyzed using the automated NCPA algorithm. Results indicate that NCPA tends to give a high frequency of false positives. In our simulations we observe that 14% of the clades give a conclusive inference that a demographic event has occurred, and that 75% of the datasets have at least one clade that gives such an inference. This is mainly due to the generation of multiple statistics per clade, of which only one is required to be significant to apply the inference key. We survey the inferences that have been made in recent publications and show that the most commonly inferred processes (restricted gene flow with isolation by distance and contiguous range expansion) are those that are commonly inferred in our simulations. However, published datasets typically yield a richer set of inferences with NCPA than obtained in our random-mating simulations, and further testing of NCPA with models of structured populations is necessary to examine its accuracy.

Implications of climate change for expanding cities worldwide

This paper analyses the trends of the changing environmental effects within growing megacities as their diameters exceed 50–100 km and their populations rise beyond 30 million people. The authors consider how these effects are influenced by climate change, to which urban areas themselves contribute, caused by their increasing greenhouse gas emissions associated with rapidly expanding energy use. Other environmental and social factors are assessed, quantitatively and qualitatively, using detailed modelling of urban mesoscale meteorology, which shows how these factors can lead to large conurbations becoming more vulnerable to climatic and environmental hazards. The paper discusses the likely changes in meteorological and hydrological hazards in urban areas, both as the climate changes and the sizes of urban areas grow. Examples are given of how these risks are being reduced through innovations in warning and response systems, planning and infrastructure design, which should include refuges against extreme natural disasters. Policies are shown to be more effective when they are integrated and based on substantial community involvement. Some conclusions are drawn regarding how policies for the natural and artificial environment and for reducing many kinds of climate and hazard risk are related to future designs and planning of infrastructure and open spaces.

Towards global drought early warning capability: expanding international cooperation for the development of a framework for global drought monitoring and forecasting

Drought is a global problem that has far-reaching impacts and especially 47 on vulnerable populations in developing regions. This paper highlights the need for a Global Drought Early Warning System (GDEWS), the elements that constitute its underlying framework (GDEWF) and the recent progress made towards its development. Many countries lack drought monitoring systems, as well as the capacity to respond via appropriate political, institutional and technological frameworks, and these have inhibited the development of integrated drought management plans or early warning systems. The GDEWS will provide a source of drought tools and products via the GDEWF for countries and regions to develop tailored drought early warning systems for their own users. A key goal of a GDEWS is to maximize the lead time for early warning, allowing drought managers and disaster coordinators more time to put mitigation measures in place to reduce the vulnerability to drought. To address this, the GDEWF will take both a top-down approach to provide global real-time drought monitoring and seasonal forecasting, and a bottom-up approach that builds upon existing national and regional systems to provide continental to global coverage. A number of challenges must be overcome, however, before a GDEWS can become a reality, including the lack of in-situ measurement networks and modest seasonal forecast skill in many regions, and the lack of infrastructure to translate data into useable information. A set of international partners, through a series of recent workshops and evolving collaborations, has made progress towards meeting these challenges and developing a global system.

Landscape structure mediates the effects of a stressor on field vole populations

Spatio-temporal landscape heterogeneity has rarely been considered in population-level impact assessments. Here we test whether landscape heterogeneity is important by examining the case of a pesticide applied seasonally to orchards which may affect non-target vole populations, using a validated ecologically realistic and spatially explicit agent-based model. Voles thrive in unmanaged grasslands and untreated orchards but are particularly exposed to applied pesticide treatments during dispersal between optimal habitats. We therefore hypothesised that vole populations do better (1) in landscapes containing more grassland and (2) where areas of grassland are closer to orchards, but (3) do worse if larger areas of orchards are treated with pesticide. To test these hyposeses we made appropriate manipulations to a model landscape occupied by field voles. Pesticide application reduced model population sizes in all three experiments, but populations subsequently wholly or partly recovered. Population depressions were, as predicted, lower in landscapes containing more unmanaged grassland, in landscapes with reduced distance between grassland and orchards, and in landscapes with fewer treated orchards. Population recovery followed a similar pattern except for an unexpected improvement in recovery when the area of treated orchards was increased. Outside the period of pesticide application, orchards increase landscape connectivity and facilitate vole dispersal and so speed population recovery. Overall our results show that accurate prediction of population impact cannot be achieved without taking account of landscape structure. The specifics of landscape structure and habitat connectivity are likely always important in mediating the effects of stressors.

Accumulated metal speciation in earthworm populations with multi-generational exposure to metalliferous soils: cell fractionation and high energy synchrotron analyses

Predicting metal bioaccumulation and toxicity in soil organisms is complicated by site-specific biotic and abiotic parameters. In this study we exploited tissue fractionation and digestion techniques, combined with X-ray absorption spectroscopy (XAS), to investigate the whole-body and subcellular distributions, ligand affinities, and coordination chemistry of accumulated Pb and Zn in field populations of the epigeic earthworm Lumbricus rubellus inhabiting three contrasting metalliferous and two unpolluted soils. Our main findings were (i) earthworms were resident in soils with concentrations of Pb and Zn ranging from 1200 to 27 000 mg kg(-1) and 200 to 34 000 mg kg(-1), respectively; (ii) Pb and Zn primarily accumulated in the posterior alimentary canal in nonsoluble subcellular fractions of earthworms; (iii) site-specific differences in the tissue and subcellular partitioning profiles of populations were observed, with earthworms from a calcareous site partitioning proportionally more Pb to their anterior body segments and Zn to the chloragosome-rich subcellular fraction than their acidic-soil inhabiting counterparts; (iv) XAS indicated that the interpopulation differences in metal partitioning between organs were not accompanied by qualitative differences in ligand-binding speciation, because crystalline phosphate-containing pyromorphite was a predominant chemical species in the whole-worm tissues of all mine soil residents. Differences in metal (Pb, Zn) partitioning at both organ and cellular levels displayed by field populations with protracted histories of metal exposures may reflect their innate ecophysiological responses to essential edaphic variables, such as Ca2+ status. These observations are highly significant in the challenging exercise of interpreting holistic biomarker data delivered by "omic" technologies.

Comparison of subcellular partitioning, distribution, and internal speciation of Cu between Cu-tolerant and naive populations of Dendrodrilus rubidus Savigny

When considering contaminated site ecology and ecological risk assessment a key question is whether organisms that appear unaffected by accumulation of contaminants are tolerant or resistant to those contaminants. A population of Dendrodrilus rubidus Savigny earthworms from the Coniston Copper Mines, an area of former Cu mining, exhibit increased tolerance and accumulation of Cu relative to a nearby non-Cu exposed population. Distribution of total Cu between different body parts (posterior, anterior, body wall) of the two populations was determined after a 14 day exposure to 250 mg Cu kg(-1) in Cu-amended soil. Cu concentrations were greater in Coniston earthworms but relative proportions of Cu in different body parts were the same between populations. Cu speciation was determined using extended X-ray absorption fine structure spectroscopy (EXAFS). Cu was coordinated to 0 atoms in the exposure soil but to S atoms in the earthworms. There was no difference in this speciation between the different earthworm populations. In another experiment earthworms were exposed to a range of Cu concentrations (200-700 mg Cu kg(-1)). Subcellular partitioning of accumulated Cu was determined. Coniston earthworms accumulated more Cu but relative proportions of Cu in the different fractions (cytosol > granular > tissue fragments, cell membranes, and intact cells) were the same between populations. Results suggest that Coniston D. rubidus are able to survive in the Cu-rich Coniston Copper Mines soil through enlargement of the same Cu storage reservoirs that exist in a nearby non-Cu exposed population.

Identifying management zones in agricultural fields using spatially constrained classification of soil and ancillary data

Site-specific management requires accurate knowledge of the spatial variation in a range of soil properties within fields. This involves considerable sampling effort, which is costly. Ancillary data, such as crop yield, elevation and apparent electrical conductivity (ECa) of the soil, can provide insight into the spatial variation of some soil properties. A multivariate classification with spatial constraint imposed by the variogram was used to classify data from two arable crop fields. The yield data comprised 5 years of crop yield, and the ancillary data 3 years of yield data, elevation and ECa. Information on soil chemical and physical properties was provided by intensive surveys of the soil. Multivariate variograms computed from these data were used to constrain sites spatially within classes to increase their contiguity. The constrained classifications resulted in coherent classes, and those based on the ancillary data were similar to those from the soil properties. The ancillary data seemed to identify areas in the field where the soil is reasonably homogeneous. The results of targeted sampling showed that these classes could be used as a basis for management and to guide future sampling of the soil.

Bryozoan populations reflect nutrient enrichment and productivity gradients in rivers

1. The hypothesis that nutrient enrichment will affect bryozoan abundance was tested using two complementary investigations; a field-based method determining bryozoan abundance in 20 rivers of different nutrient concentrations by deploying statoblast (dormant propagule) traps and an experimental laboratory microcosm study measuring bryozoan growth and mortality. These two methods confirmed independently that increased nutrient concentrations in water promote increases in the biomass of freshwater bryozoans. 2. Statoblasts of the genus Plumatella were recorded in all rivers, regardless of nutrient concentrations, demonstrating that freshwater bryozoans are widespread. Concentrations of Plumatella statoblasts were high in rivers with high nutrient concentrations relative to those with low to moderate nutrient concentrations. Regression analyses indicated that phosphorus concentrations, in particular, significantly influenced statoblast concentrations. 3. Concentrations of Lophopus crystallinus statoblasts were also higher in sites characterised by high nutrient concentrations. Logistic regression analysis revealed that the presence of L. crystallinus statoblasts was significantly associated with decreasing altitude and increasing phosphorus concentrations. This apparently rare species was found in nine rivers (out of 20), seven of which were new sites for L. crystallinus. 4. Growth rates of Fredericella sultana in laboratory microcosms increased with increasing nutrient concentration and high mortality rates were associated with low nutrient concentrations. 5. Our results indicate that bryozoans respond to increasing nutrient concentrations by increased growth, resulting in higher biomasses in enriched waters. We also found that an important component of bryozoan diets can derive from food items lacking chlorophyll a. Finally, bryozoans may be used as independent proxies for inferring trophic conditions, a feature that may be especially valuable in reconstructing historical environments by assessing the abundance of statoblasts in sediment cores.

Inherited resistance to arsenate toxicity in two populations of Lumbricus rubellus

No unequivocal evidence exists of genetically inherited resistance to metals/metalloids in field populations of earthworms. We studied cocoon production in adult Lumbricus rubellus Hoffmeister collected from an abandoned arsenic and copper mine (Devon Great Consols, Devon, UK), and abandoned tungsten mine (Carrock Fell, Cumbria, UK) and an uncontaminated cultured population. The earthworms were kept in uncontaminated soil for nine weeks. From a total of 42 L. rubellus from each site, Devon Great Consols adults produced 301 cocoons, of which 42 were viable; Carrock Fell 60 cocoons, of which 11 were viable; and the reference population 101 cocoons, of which 62 were viable. The hatchlings were collected and stored at 4degreesC at weekly intervals. After 12 weeks, all hatchlings were transferred to clean soil and maintained at 15degreesC for 20 weeks until they showed evidence of a clitellum. In toxicity trials, F1 generation L. rubellus were exposed to 2,000 mg As/kg as sodium arsenate or 300 mg Cu/kg as copper chloride for 28 d. The F1 generation L. rubellus from Devon Great Consols mine demonstrated resistance to arsenate but not copper. All L. rubellus from Devon Great Consols kept in soil treated with sodium arsenate remained in good condition over the 28-d period but lost condition rapidly and suffered high mortality in soil treated with copper chloride. The control population suffered high mortality in soil treated with sodium arsenate and copper chloride. Previous work has shown that field-collected adults demonstrate resistance to both arsenate and Cu toxicity under these conditions. Thus, while arsenate resistance may be demonstrated in F1 generation L. rubellus from one of the contaminated sites, Cu resistance is not. The F1 adults and F2 cocoons did not have significantly higher levels of As than the control population, with no residual As tissue burden, suggesting that resistance to As in these populations may be inherited.

Dynamic scaling of bred vectors in spatially extended chaotic systems

We unfold a profound relationship between the dynamics of finite-size perturbations in spatially extended chaotic systems and the universality class of Kardar-Parisi-Zhang (KPZ). We show how this relationship can be exploited to obtain a complete theoretical description of the bred vectors dynamics. The existence of characteristic length/time scales, the spatial extent of spatial correlations and how to time it, and the role of the breeding amplitude are all analyzed in the light of our theory. Implications to weather forecasting based on ensembles of initial conditions are also discussed.