Bayesian estimation of recent migration rates after a spatial expansion

Approximate Bayesian computation (ABC) is a highly flexible technique that allows the estimation of parameters under demographic models that are too complex to be handled by full-likelihood methods. We assess the utility of this method to estimate the parameters of range expansion in a two-dimensional stepping-stone model, using samples from either a single deme or multiple demes. A minor modification to the ABC procedure is introduced, which leads to an improvement in the accuracy of estimation. The method is then used to estimate the expansion time and migration rates for five natural common vole populations in Switzerland typed for a sex-linked marker and a nuclear marker. Estimates based on both markers suggest that expansion occurred < 10,000 years ago, after the most recent glaciation, and that migration rates are strongly male biased.

Extreme host plant conservatism during at least 20 million years of host plant pursuit by oak gallwasps

Diversification of insect herbivores is often associated with coevolution between plant toxins and insect countermeasures, resulting in a specificity that restricts host plant shifts. Gall inducers, however, bypass plant toxins and the factors influencing host plant associations in these specialized herbivores remain unclear. We reconstructed the evolution of host plant associations in Western Palaearctic oak gallwasps (Cynipidae: Cynipini), a species-rich lineage of specialist herbivores on oak (Quercus). (1) Bayesian analyses of sequence data for three genes revealed extreme host plant conservatism, with inferred shifts between major oak lineages (sections Cerris and Quercus) closely matching the minimum required to explain observed diversity. It thus appears that the coevolutionary demands of gall induction constrain host plant shifts, both in cases of mutualism (e.g., fig wasps, yucca moths) and parasitism (oak gallwasps). (2) Shifts between oak sections occurred independently in sexual and asexual generations of the gallwasp lifecycle, implying that these can evolve independently. (3) Western Palaearctic gallwasps associated with sections Cerris and Quercus diverged at least 20 million years ago (mya), prior to the arrival of oaks in the Western Palaearctic from Asia 5-7 mya. This implies an Asian origin for Western Palaearctic gallwasps, with independent westwards range expansion by multiple lineages.

Statistical evaluation of alternative models of human evolution

An appropriate model of recent human evolution is not only important to understand our own history, but it is necessary to disentangle the effects of demography and selection on genome diversity. Although most genetic data support the view that our species originated recently in Africa, it is still unclear if it completely replaced former members of the Homo genus, or if some interbreeding occurred during its range expansion. Several scenarios of modern human evolution have been proposed on the basis of molecular and paleontological data, but their likelihood has never been statistically assessed. Using DNA data from 50 nuclear loci sequenced in African, Asian and Native American samples, we show here by extensive simulations that a simple African replacement model with exponential growth has a higher probability (78%) as compared with alternative multiregional evolution or assimilation scenarios. A Bayesian analysis of the data under this best supported model points to an origin of our species approximate to 141 thousand years ago (Kya), an exit out-of-Africa approximate to 51 Kya, and a recent colonization of the Americas approximate to 10.5 Kya. We also find that the African replacement model explains not only the shallow ancestry of mtDNA or Y-chromosomes but also the occurrence of deep lineages at some autosomal loci, which has been formerly interpreted as a sign of interbreeding with Homo erectus.

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.

The population genetic structure of clonal organisms generated by exponentially bounded and fat-tailed dispersal

Long distance dispersal (LDD) plays an important role in many population processes like colonization, range expansion, and epidemics. LDD of small particles like fungal spores is often a result of turbulent wind dispersal and is best described by functions with power-law behavior in the tails ("fat tailed"). The influence of fat-tailed LDD on population genetic structure is reported in this article. In computer simulations, the population structure generated by power-law dispersal with exponents in the range of -2 to -1, in distinct contrast to that generated by exponential dispersal, has a fractal structure. As the power-law exponent becomes smaller, the distribution of individual genotypes becomes more self-similar at different scales. Common statistics like G(ST) are not well suited to summarizing differences between the population genetic structures. Instead, fractal and self-similarity statistics demonstrated differences in structure arising from fat-tailed and exponential dispersal. When dispersal is fat tailed, a log-log plot of the Simpson index against distance between subpopulations has an approximately constant gradient over a large range of spatial scales. The fractal dimension D-2 is linearly inversely related to the power-law exponent, with a slope of similar to -2. In a large simulation arena, fat-tailed LDD allows colonization of the entire space by all genotypes whereas exponentially bounded dispersal eventually confines all descendants of a single clonal lineage to a relatively small area.

Host and habitat filtering in seedling root-associated fungal communities: taxonomic and functional diversity are altered in ‘novel’ soils

Climatic and land use changes have significant consequences for the distribution of tree species, both through natural dispersal processes and following management prescriptions. Responses to these changes will be expressed most strongly in seedlings near current species range boundaries. In northern temperate forest ecosystems, where changes are already being observed, ectomycorrhizal fungi contribute significantly to successful tree establishment. We hypothesised that communities of fungal symbionts might therefore play a role in facilitating, or limiting, host seedling range expansion. To test this hypothesis, ectomycorrhizal communities of interior Douglas-fir and interior lodgepole pine seedlings were analysed in a common greenhouse environment following growth in five soils collected along an ecosystem gradient. Currently, Douglas-fir’s natural distribution encompasses three of the five soils, whereas lodgepole pine’s extends much further north. Host filtering was evident amongst the 29 fungal species encountered: 7 were shared, 9 exclusive to Douglas-fir and 13 exclusive to lodgepole pine. Seedlings of both host species formed symbioses with each soil fungal community, thus Douglas-fir did so even where those soils came from outside its current distribution. However, these latter communities displayed significant taxonomic and functional differences to those found within the host distribution, indicative of habitat filtering. In contrast, lodgepole pine fungal communities displayed high functional similarity across the soil gradient. Taxonomic and/or functional shifts in Douglas-fir fungal communities may prove ecologically significant during the predicted northward migration of this species; especially in combination with changes in climate and management operations, such as seed transfer across geographical regions for forestry purposes.

Can site and landscape-scale environmental attributes buffer bird populations against weather events?

Projected impacts of climate change on the populations and distributions of species pose a challenge for conservationists. In response, a number of adaptation strategies to enable species to persist in a changing climate have been proposed. Management to maximise the quality of habitat at existing sites may reduce the magnitude or frequency of climate-driven population declines. In addition large-scale management of landscapes could potentially improve the resilience of populations by facilitating inter-population movements. A reduction in the obstacles to species’ range expansion, may also allow species to track changing conditions better through shifts to new locations, either regionally or locally. However, despite a strong theoretical base, there is limited empirical evidence to support these management interventions. This makes it difficult for conservationists to decide on the most appropriate strategy for different circumstances. Here extensive data from long-term monitoring of woodland birds at individual sites are used to examine the two-way interactions between habitat and both weather and population count in the previous year. This tests the extent to which site-scale and landscape-scale habitat attributes may buffer populations against variation in winter weather (a key driver of woodland bird population size) and facilitate subsequent population growth. Our results provide some support for the prediction that landscape-scale attributes (patch isolation and area of woodland habitat) may influence the ability of some woodland bird species to withstand weather-mediated population declines. These effects were most apparent among generalist woodland species. There was also evidence that several, primarily specialist, woodland species are more likely to increase following population decline where there is more woodland at both site and landscape scales. These results provide empirical support for the concept that landscape-scale conservation efforts may make the populations of some woodland bird species more resilient to climate change. However in isolation, management is unlikely to provide a universal benefit to all species.

The ancient Britons: groundwater fauna survived extreme climate change over tens of millions of years across NW Europe

Global climate changes during the Cenozoic (65.5–0 Ma) caused major biological range shifts and extinctions. In northern Europe, for example, a pattern of few endemics and the dominance of wide-ranging species is thought to have been determined by the Pleistocene (2.59–0.01 Ma) glaciations. This study, in contrast, reveals an ancient subsurface fauna endemic to Britain and Ireland. Using a Bayesian phylogenetic approach, we found that two species of stygobitic invertebrates (genus Niphargus) have not only survived the entire Pleistocene in refugia but have persisted for at least 19.5 million years. Other Niphargus species form distinct cryptic taxa that diverged from their nearest continental relative between 5.6 and 1.0 Ma. The study also reveals an unusual biogeographical pattern in the Niphargus genus. It originated in north-west Europe approximately 87 Ma and underwent a gradual range expansion. Phylogenetic diversity and species age are highest in north-west Europe, suggesting resilience to extreme climate change and strongly contrasting the patterns seen in surface fauna. However, species diversity is highest in south-east Europe, indicating that once the genus spread to these areas (approximately 25 Ma), geomorphological and climatic conditions enabled much higher diversification. Our study highlights that groundwater ecosystems provide an important contribution to biodiversity and offers insight into the interactions between biological and climatic processes.

"Ever-growing Amman", Jordan: Urban expansion, social polarisation and contemporary urban planning issues

Amman the primate capital city of the Hashemite Kingdom of Jordan currently has a population in excess of 2 million, but in 1924 it consisted of little more than a collection of dwellings and some 2000-3000 inhabitants. The present paper sets out to document and explain the phenomenal expansion of "ever-growing Amman". The physical geography of the urban region and the early growth of the city are considered at the outset and this leads directly to consideration of the highly polarised social structuring that characterises contemporary Amman. In doing this, original data derived from the recent Greater Amman Municipality's Geographical Information System are presented. In this respect, the essential modernity of the city is exemplified. The employment and industrial bases of the city and a range of pressing contemporary issues are then considered, including transport and congestion, the provision of urban water under conditions of water stress and privatisation, and urban and regional development planning for the city. The paper concludes by emphasizing the growing regional and international geopolitical salience of the city of Amman at the start of the 21st century. (C) 2008 Elsevier Ltd. All rights reserved.

Acclimation of photosynthesis to elevated CO2 in onion (Allium cepa) grown at a range of temperatures

Onion (Allium cepa) was grown in the field within temperature gradient tunnels (providing about -2.5degreesC to +2.5degreesC from outside temperatures) maintained at either 374 or 532 mumol mol(-1) CO2. Plant leaf area was determined non-destructively at 7 day intervals until the time of bulbing in 12 combinations of temperature and CO2 concentration. Gas exchange was measured in each plot at the time of bulbing, and the carbohydrate content of the leaf (source) and bulb (sink) was determined. Maximum rate of leaf area expansion increased with mean temperature. Leaf area duration and maximum rate of leaf area expansion were not significantly affected by CO2. The light-saturated rates of leaf photosynthesis (A(sat)) were greater in plants grown at normal than at elevated CO2 concentrations at the same measurement CO2 concentration. Acclimation of photosynthesis decreased with an increase in growth temperature, and with an increase in leaf nitrogen content at elevated CO2. The ratio of intercellular to atmospheric CO2 (C-i/C-a ratio) was 7.4% less for plants grown at elevated compared with normal CO2. A(sat) in plants grown at elevated CO2 was less than in plants grown at normal CO2 when compared at the same C-i Hence, acclimation of photosynthesis was due both to stomatal acclimation and to limitations to biochemical CO2 fixation. Carbohydrate content of the onion bulbs was greater at elevated than at normal CO2. In contrast, carbohydrate content was less at elevated compared with normal CO2 in the leaf sections in which CO2 exchange was measured at the same developmental stage. Therefore, acclimation of photosynthesis in fully expanded onion leaves was detected despite the absence of localised carbohydrate accumulation in these field-grown crops.

Genesis and expansion of metazoan transcription factor gene classes

We know little about the genomic events that led to the advent of a multicellular grade of organization in animals, one of the most dramatic transitions in evolution. Metazoan multicellularity is correlated with the evolution of embryogenesis, which presumably was underpinned by a gene regulatory network reliant on the differential activation of signaling pathways and transcription factors. Many transcription factor genes that play critical roles in bilaterian development largely appear to have evolved before the divergence of cnidarian and bilaterian lineages. In contrast, sponges seem to have a more limited suite of transcription factors, suggesting that the developmental regulatory gene repertoire changed markedly during early metazoan evolution. Using whole- genome information from the sponge Amphimedon queenslandica, a range of eumetazoans, and the choanoflagellate Monosiga brevicollis, we investigate the genesis and expansion of homeobox, Sox, T- box, and Fox transcription factor genes. Comparative analyses reveal that novel transcription factor domains ( such as Paired, POU, and T- box) arose very early in metazoan evolution, prior to the separation of extant metazoan phyla but after the divergence of choanoflagellate and metazoan lineages. Phylogenetic analyses indicate that transcription factor classes then gradually expanded at the base of Metazoa before the bilaterian radiation, with each class following a different evolutionary trajectory. Based on the limited number of transcription factors in the Amphimedon genome, we infer that the genome of the metazoan last common ancestor included fewer gene members in each class than are present in extant eumetazoans. Transcription factor orthologues present in sponge, cnidarian, and bilaterian genomes may represent part of the core metazoan regulatory network underlying the origin of animal development and multicellularity.

Acclimation of photosynthesis to elevated CO2 in onion (allium cepa) grown at a range of temperatures

Onion (Allium cepa) was grown in the field within temperature gradient tunnels (providing about -2.5 degrees C to +2.5 degrees C from outside temperatures) maintained at either 374 or 532 mumol mol (-1) CO2. Plant leaf area was determined non-destructively at 7 day intervals until the time of bulbing in 12 combinations of temperature and CO2 concentration. Gas exchange was measured in each plot at the time of bulbing, and the carbohydrate content of the leaf (source) and bulb (sink) was determined. Maximum rate of leaf area expansion increased with mean temperature. Leaf area duration and maximum rate of leaf area expansion were not significantly affected by CO2. The light-saturated rates of leaf photosynthesis (A(sat)) were greater in plants grown at normal than at elevated CO2 concentrations at the same measurement CO2 concentration. Acclimation of photosynthesis decreased with an increase in growth temperature, and with an increase in leaf nitrogen content at elevated CO2. The ratio of intercellular to atmospheric CO2 (C-i/C-a ratio) was 7.4% less for plants grown at elevated compared with normal CO2. A(sat) in plants grown at elevated CO2 was less than in plants grown at normal CO2 when compared at the same C-i Hence, acclimation of photosynthesis was due both to stomatal acclimation and to limitations to biochemical CO2 fixation. Carbohydrate content of the onion bulbs was greater at elevated than at normal CO2. In contrast, carbohydrate content was less at elevated compared with normal CO2 in the leaf sections in which CO2 exchange was measured at the same developmental stage. Therefore, acclimation of photosynthesis in fully expanded onion leaves was detected despite the absence of localised carbohydrate accumulation in these field-grown crops.

A differential thermal expansion approach to crystal structure determination from powder diffraction data

Differential thermal expansion over the range 90-210 K has been applied successfully to determine the crystal structure of chlorothiazide from synchrotron powder diffraction data using direct methods. Key to the success of the approach is the use of a multi-data-set Pawley refinement to extract a set of reflection intensities that is more 'single-crystal-like' than those extracted from a single data set. The improvement in reflection intensity estimates is quantified by comparison with reference single-crystal intensities. (C) 2008 International Union of Crystallography Printed in Singapore - all rights reserved

Structures and negative thermal expansion properties of the one-dimensional cyanides, CuCN, AgCN and AuCN

The behaviour of the lattice parameters of HTCuCN (high-temperature form), AgCN and AuCN have been investigated as a function of temperature over the temperature range 90–490 K. All materials show one-dimensional negative thermal expansion (NTE) along the ––(M––CN)–– chain direction c (ac(HT-CuCN) ¼32.1 10–6 K1, ac(AgCN)¼23.910–6 K1 and ac(AuCN) ¼9.3106 K1 over the temperature range 90–490 K). The origin of this behaviour has been studied using RMC modelling of Bragg and total neutron diffraction data from AgCN and AuCN at 10 and 300 K. These analyses yield details of the local motions within the chains responsible for NTE. The low-temperature form of CuCN, LT-CuCN, has been studied using single-crystal X-ray diffraction. In this form of CuCN, wavelike distortions of the ––(Cu––CN)–– chains occur in the static structure, which are reminiscent of the motions seen in the RMC modelling of AgCN and AuCN, which are responsible for the NTE behaviour.

Relationship between phonons and thermal expansion in Zn(CN)2 and Ni(CN)2 from inelastic neutron scattering and ab initio calculations

Zn(CN)2 and Ni(CN)2 are known for exhibiting anomalous thermal expansion over a wide temperature range. The volume thermal expansion coefficient for the cubic, three dimensionally connected material, Zn(CN)2, is negative (alpha(V) = −51  10(-6) K-1) while for Ni(CN)2, a tetragonal material, the thermal expansion coefficient is negative in the two dimensionally connected sheets (alpha(a) = −7  10(-6) K-1), but the overall thermal expansion coefficient is positive (alpha(V) = 48  10(-6) K-1). We have measured the temperature dependence of phonon spectra in these compounds and analyzed them using ab initio calculations. The spectra of the two compounds show large differences that cannot be explained by simple mass renormalization of the modes involving Zn (65.38 amu) and Ni (58.69 amu) atoms. This reflects the fact that the structure and bonding are quite different in the two compounds. The calculated pressure dependence of the phonon modes and of the thermal expansion coefficient, alpha(V), are used to understand the anomalous behavior in these compounds. Our ab initio calculations indicate that phonon modes of energy approx. 2 meV are major contributors to negative thermal expansion (NTE) in both the compounds. The low-energy modes of approx.8 and 13 meV in Zn(CN)2 also contribute significantly to the NTE in Zn(CN)2 and Ni(CN)2, respectively. The measured temperature dependence of the phonon spectra has been used to estimate the total anharmonicity of both compounds. For Zn(CN)2, the temperature-dependent measurements (total anharmonicity), along with our previously reported pressure dependence of the phonon spectra (quasiharmonic), is used to separate the explicit temperature effect at constant volume (intrinsic anharmonicity).