Adaptive population divergence and directional gene flow across steep elevational gradients in a climate‐sensitive mammal

The ecological effects of climate change have been shown in most major taxonomic groups; however, the evolutionary consequences are less well‐documented. Adaptation to new climatic conditions offers a potential long‐term mechanism for species to maintain viability in rapidly changing environments, but mammalian examples remain scarce. The American pika (Ochotona princeps) has been impacted by recent climate‐associated extirpations and range‐wide reductions in population sizes, establishing it as a sentinel mammalian species for climate change. To investigate evidence for local adaptation and reconstruct patterns of genomic diversity and gene flow across rapidly changing environments, we used a space‐for‐time design and restriction site‐associated DNA sequencing to genotype American pikas along two steep elevational gradients at 30,966 SNPs and employed independent outlier detection methods that scanned for genotype‐environment associations. We identified 338 outlier SNPs detected by two separate analyses and/or replicated in both transects, several of which were annotated to genes involved in metabolic function and oxygen transport. Additionally, we found evidence of directional gene flow primarily downslope from high‐elevation populations, along with reduced gene flow at outlier loci. If this trend continues, elevational range contractions in American pikas will likely be from local extirpation rather than upward movement of low‐elevation individuals; this, in turn, could limit the potential for adaptation within this landscape. These findings are of particular relevance for future conservation and management of American pikas and other elevationally restricted, thermally sensitive species.     

Sap‐feeding insects on forest trees along latitudinal gradients in northern Europe: a climate‐driven patterns

Knowledge of the latitudinal patterns in biotic interactions, and especially in herbivory, is crucial for understanding the mechanisms that govern ecosystem functioning and for predicting their responses to climate change. We used sap‐feeding insects as a model group to test the hypotheses that the strength of plant–herbivore interactions in boreal forests decreases with latitude and that this latitudinal pattern is driven primarily by midsummer temperatures. We used a replicated sampling design and quantitatively collected and identified all sap‐feeding insects from four species of forest trees along five latitudinal gradients (750–1300 km in length, ten sites in each gradient) in northern Europe (59 to 70°N and 10 to 60°E) during 2008–2011. Similar decreases in diversity of sap‐feeding insects with latitude were observed in all gradients during all study years. The sap‐feeder load (i.e. insect biomass per unit of foliar biomass) decreased with latitude in typical summers, but increased in an exceptionally hot summer and was independent of latitude during a warm summer. Analysis of combined data from all sites and years revealed dome‐shaped relationships between the loads of sap‐feeders and midsummer temperatures, peaking at 17 °C in Picea abies, at 19.5 °C in Pinus sylvestris and Betula pubescens and at 22 °C in B. pendula. From these relationships, we predict that the losses of forest trees to sap‐feeders will increase by 0–45% of the current level in southern boreal forests and by 65–210% in subarctic forests with a 1 °C increase in summer temperatures. The observed relationships between temperatures and the loads of sap‐feeders differ between the coniferous and deciduous tree species. We conclude that climate warming will not only increase plant losses to sap‐feeding insects, especially in subarctic forests, but can also alter plant‐plant interactions, thereby affecting both the productivity and the structure of future forest ecosystems.     

American pika in a low‐elevation lava landscape: expanding the known distribution of a temperature‐sensitive species

In 2010, the American pika (Ochotona princeps fenisex) was denied federal protection based on limited evidence of persistence in low‐elevation environments. Studies in nonalpine areas have been limited to relatively few environments, and it is unclear whether patterns observed elsewhere (e.g., Bodie, CA) represent other nonalpine habitats. This study was designed to establish pika presence in a new location, determine distribution within the surveyed area, and evaluate influences of elevation, vegetation, lava complexity, and distance to habitat edge on pika site occupancy. In 2011 and 2012, we conducted surveys for American pika on four distinct subalpine lava flows of Newberry National Volcanic Monument, Oregon, USA. Field surveys were conducted at predetermined locations within lava flows via silent observation and active searching for pika sign. Site habitat characteristics were included as predictors of occupancy in multinomial regression models. Above and belowground temperatures were recorded at a subsample of pika detection sites. Pika were detected in 26% (2011) and 19% (2012) of survey plots. Seventy‐four pika were detected outside survey plot boundaries. Lava complexity was the strongest predictor of pika occurrence, where pika were up to seven times more likely to occur in the most complicated lava formations. Pika were two times more likely to occur with increasing elevation, although they were found at all elevations in the study area. This study expands the known distribution of the species and provides additional evidence for persistence in nonalpine habitats. Results partially support the predictive occupancy model developed for pika at Craters of the Moon National Monument, another lava environment. Characteristics of the lava environment clearly influence pika site occupancy, but habitat variables reported as important in other studies were inconclusive here. Further work is needed to gain a better understanding of the species’ current distribution and ability to persist under future climate conditions.     

Landscape effects on gene flow for a climate‐sensitive montane species,the American pika

Climate change is arguably the greatest challenge to conservation of our time. Most vulnerability assessments rely on past and current species distributions to predict future persistence but ignore species' abilities to disperse through landscapes, which may be particularly important in fragmented habitats and crucial for long‐term persistence in changing environments. Landscape genetic approaches explore the interactions between landscape features and gene flow and can clarify how organisms move among suitable habitats, but have suffered from methodological uncertainties. We used a landscape genetic approach to determine how landscape and climate‐related features influence gene flow for American pikas (Ochotona princeps) in Crater Lake National Park. Pikas are heat intolerant and restricted to cool microclimates; thus, range contractions have been predicted as climate changes. We evaluated the correlation between landscape variables and genetic distance using partial Mantel tests in a causal modelling framework, and used spatially explicit simulations to evaluate methods of model optimization including a novel approach based on relative support and reciprocal causal modelling. We found that gene flow was primarily restricted by topographic relief, water and west‐facing aspects, suggesting that physical restrictions related to small body size and mode of locomotion, as well as exposure to relatively high temperatures, limit pika dispersal in this alpine habitat. Our model optimization successfully identified landscape features influencing resistance in the simulated data for this landscape, but underestimated the magnitude of resistance. This is the first landscape genetic study to address the fundamental question of what limits dispersal and gene flow in the American pika.     

Genomic and phenotypic differentiation of Arabidopsis thaliana along altitudinal gradients in the North Italian Alps

Altitudinal gradients in mountain regions are short‐range clines of different environmental parameters such as temperature or radiation. We investigated genomic and phenotypic signatures of adaptation to such gradients in five Arabidopsis thaliana populations from the North Italian Alps that originated from 580 to 2350 m altitude by resequencing pools of 19–29 individuals from each population. The sample includes two pairs of low‐ and high‐altitude populations from two different valleys. High‐altitude populations showed a lower nucleotide diversity and negative Tajima's D values and were more closely related to each other than to low‐altitude populations from the same valley. Despite their close geographic proximity, demographic analysis revealed that low‐ and high‐altitude populations split between 260 000 and 15 000 years before present. Single nucleotide polymorphisms whose allele frequencies were highly differentiated between low‐ and high‐altitude populations identified genomic regions of up to 50 kb length where patterns of genetic diversity are consistent with signatures of local selective sweeps. These regions harbour multiple genes involved in stress response. Variation among populations in two putative adaptive phenotypic traits, frost tolerance and response to light/UV stress was not correlated with altitude. Taken together, the spatial distribution of genetic diversity reflects a potentially adaptive differentiation between low‐ and high‐altitude populations, whereas the phenotypic differentiation in the two traits investigated does not. It may resemble an interaction between adaptation to the local microhabitat and demographic history influenced by historical glaciation cycles, recent seed dispersal and genetic drift in local populations.     

Habitat continuity and stepping‐stone oceanographic distances explain population genetic connectivity of the brown alga Cystoseira amentacea

Effective predictive and management approaches for species occurring in a metapopulation structure require good understanding of interpopulation connectivity. In this study, we ask whether population genetic structure of marine species with fragmented distributions can be predicted by stepping‐stone oceanographic transport and habitat continuity, using as model an ecosystem‐structuring brown alga, Cystoseira amentacea var. stricta. To answer this question, we analysed the genetic structure and estimated the connectivity of populations along discontinuous rocky habitat patches in southern Italy, using microsatellite markers at multiple scales. In addition, we modelled the effect of rocky habitat continuity and ocean circulation on gene flow by simulating Lagrangian particle dispersal based on ocean surface currents allowing multigenerational stepping‐stone dynamics. Populations were highly differentiated, at scales from few metres up to thousands of kilometres. The best possible model fit to explain the genetic results combined current direction, rocky habitat extension and distance along the coast among rocky sites. We conclude that a combination of variable suitable habitat and oceanographic transport is a useful predictor of genetic structure. This relationship provides insight into the mechanisms of dispersal and the role of life‐history traits. Our results highlight the importance of spatially explicit modelling of stepping‐stone dynamics and oceanographic directional transport coupled with habitat suitability, to better describe and predict marine population structure and differentiation. This study also suggests the appropriate spatial scales for the conservation, restoration and management of species that are increasingly affected by habitat modifications.     

The advantages of going large: genome‐wide SNPs clarify the complex population history and systematics of the threatened western pond turtle

As the field of phylogeography has matured, it has become clear that analyses of one or a few genes may reveal more about the history of those genes than the populations and species that are the targets of study. To alleviate these concerns, the discipline has moved towards larger analyses of more individuals and more genes, although little attention has been paid to the qualitative or quantitative gains that such increases in scale and scope may yield. Here, we increase the number of individuals and markers by an order of magnitude over previously published work to comprehensively assess the phylogeographical history of a well‐studied declining species, the western pond turtle (Emys marmorata). We present a new analysis of 89 independent nuclear SNP markers and one mitochondrial gene sequence scored for rangewide sampling of >900 individuals, and compare these to smaller‐scale, rangewide genetic and morphological analyses. Our enlarged SNP data fundamentally revise our understanding of evolutionary history for this lineage. Our results indicate that the gains from greatly increasing both the number of markers and individuals are substantial and worth the effort, particularly for species of high conservation concern such as the pond turtle, where accurate assessments of population history are a prerequisite for effective management.     

Sampling related individuals within ponds biases estimates of population structure in a pond‐breeding amphibian

Effective conservation and management of pond‐breeding amphibians depends on the accurate estimation of population structure, demographic parameters, and the influence of landscape features on breeding‐site connectivity. Population‐level studies of pond‐breeding amphibians typically sample larval life stages because they are easily captured and can be sampled nondestructively. These studies often identify high levels of relatedness between individuals from the same pond, which can be exacerbated by sampling the larval stage. Yet, the effect of these related individuals on population genetic studies using genomic data is not yet fully understood. Here, we assess the effect of within‐pond relatedness on population and landscape genetic analyses by focusing on the barred tiger salamanders (Ambystoma mavortium) from the Nebraska Sandhills. Utilizing genome‐wide SNPs generated using a double‐digest RADseq approach, we conducted standard population and landscape genetic analyses using datasets with and without siblings. We found that reduced sample sizes influenced parameter estimates more than the inclusion of siblings, but that within‐pond relatedness led to the inference of spurious population structure when analyses depended on allele frequencies. Our landscape genetic analyses also supported different models across datasets depending on the spatial resolution analyzed. We recommend that future studies not only test for relatedness among larval samples but also remove siblings before conducting population or landscape genetic analyses. We also recommend alternative sampling strategies to reduce sampling siblings before sequencing takes place. Biases introduced by unknowingly including siblings can have significant implications for population and landscape genetic analyses, and in turn, for species conservation strategies and outcomes.     

Genetic diversity and recombination in natural populations of the nematode‐trapping fungus Arthrobotrys oligospora from China

Nematophagous fungi can trap and capture nematodes and other small invertebrates. This unique ability has made them ideal organisms from which to develop biological control agents against plant‐ and animal‐parasitic nematodes. However, effective application of biocontrol agents in the field requires a comprehensive understanding about the ecology and population genetics of the nematophagous fungi in natural environments. Here, we genotyped 228 strains of the nematode‐trapping fungus Arthrobotrys oligospora using 12 single nucleotide polymorphic markers located on eight random DNA fragments. The strains were from different ecological niches and geographical regions from China. Our analyses identified that ecological niche separations contributed significantly, whereas geographic separation contributed relatively little to the overall genetic variation in our samples of A. oligospora. Interestingly, populations from stressful environments seemed to be more variable and showed more evidence for recombination than those from benign environments at the same geographic areas. We discussed the implications of our results to the conservation and biocontrol application of A. oligospora in agriculture and forestry.     

Seascape genetics along environmental gradients in the Arabian Peninsula: insights from ddRAD sequencing of anemonefishes

Understanding the processes that shape patterns of genetic structure across space is a central aim of landscape genetics. However, it remains unclear how geographical features and environmental variables shape gene flow, particularly for marine species in large complex seascapes. Here, we evaluated the genomic composition of the two‐band anemonefish Amphiprion bicinctus across its entire geographical range in the Red Sea and Gulf of Aden, as well as its close relative, Amphiprion omanensis endemic to the southern coast of Oman. Both the Red Sea and the Arabian Sea are complex and environmentally heterogeneous marine systems that provide an ideal scenario to address these questions. Our findings confirm the presence of two genetic clusters previously reported for A. bicinctus in the Red Sea. Genetic structure analyses suggest a complex seascape configuration, with evidence of both isolation by distance (IBD) and isolation by environment (IBE). In addition to IBD and IBE, genetic structure among sites was best explained when two barriers to gene flow were also accounted for. One of these coincides with a strong oligotrophic–eutrophic gradient at around 16–20?N in the Red Sea. The other agrees with a historical bathymetric barrier at the straight of Bab al Mandab. Finally, these data support the presence of interspecific hybrids at an intermediate suture zone at Socotra and indicate complex patterns of genomic admixture in the Gulf of Aden with evidence of introgression between species. Our findings highlight the power of recent genomic approaches to resolve subtle patterns of gene flow in marine seascapes.     

Two climate‐sensitive tree‐ring chronologies from Arnhem Land,monsoonal Australia

The ecology of the Australian monsoon tropics is fundamentally shaped by dry conditions between May and October followed by highly variable rainfall over the months of November to April. Due to its crucial ecological importance, a better understanding of past hydroclimate variability in the region is of great interest to land managers and custodians in this region. Short instrumental records also make highly resolved terrestrial palaeoclimate records for northern Australia prior to 1900 CE of considerable scientific importance. Here, we present two new well‐replicated Callitris intratropica ring‐width chronologies from Arnhem Land in northern Australia, one of which extends the tree‐ring record in the region by another 86 years, back to 1761. Both chronologies have clearly defined regional patterns of correlations with temperature, precipitation, potential evapotranspiration and two drought indices (the self‐calibrating Palmer Drought Severity Index (PDSI) and the Standardised Precipitation Evapotranspiration Index (SPEI)) across the lower latitudes of the Northern Territory. Results indicate considerable scope for hydroclimatic reconstructions based on C. intratropica for transitional periods into and out of the wettest time of the year. This suggests that such reconstructions would reflect variability in the duration of the wet period. While precipitation or streamflow reconstructions may be possible for both these transitional periods, drought reconstructions will be best focused on the months of March–May at the end of the wet period. Hydroclimate reconstructions would provide important baseline information for understanding the rate and magnitude of current regional climate change for these ecologically and culturally important transitional periods.     

metapop2: Re‐implementation of software for the analysis and management of subdivided populations using gene and allelic diversity

Management programmes often have to make decisions based on the analysis of the genetic properties and diversity of populations. Expected heterozygosity (or gene diversity) and population structure parameters are often used to make recommendations for conservation, such as avoidance of inbreeding or migration across subpopulations. Allelic diversity, however, can also provide complementary and useful information for conservation programmes, as it is highly sensitive to population bottlenecks, and is more related to long‐term selection response than heterozygosity. Here we present a completely revised and updated re‐implementation of the software metapop for the analysis of diversity in subdivided populations, as well as a tool for the management and dynamic estimation of optimal contributions in conservation programmes. This new update includes computation of allelic diversity for population analysis and management, as well as a simulation mode to forecast the consequences of taking different management strategies over time. Furthermore, the new implementation in C++ includes code optimization and improved memory usage, allowing for fast analysis of large data sets including single nucleotide polymorphism markers, as well as enhanced cross‐software and cross‐platform compatibility.     

Extinction,reintroduction, and restoration of a lizard meta‐population equilibrium in the missouri ozarks

In this issue of Molecular Ecology, Neuwald & Templeton (2013) report on a 22‐year study of natural populations of Collared Lizards (Crotaphytus collaris) that evolved on isolated on rock outcrops (‘glades’) in the Ozark Mountains in eastern Missouri. This ecosystem was originally maintained by frequent fires that kept the forest understory open, but fire‐suppression was adopted as official policy in about 1945, which led to a loss of native biodiversity, including local extinctions of some lizard populations. Policies aimed at restoring biodiversity included controlled burns and re‐introductions of lizards to some glades, which began in 1984. Populations were monitored from 1984–2006, and demographic and genetic data collected from 1 679 lizards were used to documents shifts in meta‐population dynamics over four distinct phases of lizard recovery: 1–an initial translocation of lizards drawn from the same source populations onto three glades that were likely part of one meta‐population; 2–a period of isolation and genetic drift associated with the absence of fires; 3–a period of rapid colonization and population increase following restoration of fire; and 4–stabilization of the meta‐population under regular prescribed burning. This study system thus provides a rare opportunity to characterize the dynamics of a landscape‐scale management strategy on the restoration of the meta‐population of a reintroduced species; long‐term case studies of the extinction, founding, increase, and stabilization of a well‐defined meta‐population, based on both demographic and population genetic data, are rare in the conservation, ecological, and evolutionary literature.     

Modelling lake cyanobacterial blooms: Disentangling the climate‐driven impacts of changing mixed depth and water temperature

相似文献   

Pelagic phytoplankton community change‐points across nutrient gradients and in response to invasive mussels

相似文献   

Signatures of local adaptation along environmental gradients in a range‐expanding damselfly (Ischnura elegans)

Insect distributions are shifting rapidly in response to climate change and are undergoing rapid evolutionary change. We investigate the molecular signatures underlying local adaptation in the range‐expanding damselfly, Ischnura elegans. Using a landscape genomic approach combined with generalized dissimilarity modelling (GDM), we detect selection signatures on loci via allelic frequency change along environmental gradients. We analyse 13,612 single nucleotide polymorphisms (SNPs), derived from restriction site‐associated DNA sequencing (RADseq), in 426 individuals from 25 sites spanning the I. elegans distribution in Sweden, including its expanding northern range edge. Environmental association analysis (EAA) and the magnitude of allele frequency change along the range expansion gradient revealed significant signatures of selection in relation to high maximum summer temperature, high mean annual precipitation and low wind speeds at the range edge. SNP annotations with significant signatures of selection revealed gene functions associated with ongoing range expansion, including heat shock proteins (HSP40 and HSP70), ion transport (V‐ATPase) and visual processes (long‐wavelength‐sensitive opsin), which have implications for thermal stress response, salinity tolerance and mate discrimination, respectively. We also identified environmental thresholds where climate‐mediated selection is likely to be strong, and indicate that I. elegans is rapidly adapting to the climatic environment during its ongoing range expansion. Our findings empirically validate an integrative approach for detecting spatially explicit signatures of local adaptation along environmental gradients.     

Fine-scale genetic population structure in a mobile marine mammal: inshore bottlenose dolphins in Moreton Bay, Australia

Highly mobile marine species in areas with no obvious geographic barriers are expected to show low levels of genetic differentiation. However, small‐scale variation in habitat may lead to resource polymorphisms and drive local differentiation by adaptive divergence. Using nuclear microsatellite genotyping at 20 loci, and mitochondrial control region sequencing, we investigated fine‐scale population structuring of inshore bottlenose dolphins (Tursiops aduncus) inhabiting a range of habitats in and around Moreton Bay, Australia. Bayesian structure analysis identified two genetic clusters within Moreton Bay, with evidence of admixture between them (F_ST = 0.05, P = 0.001). There was only weak isolation by distance but one cluster of dolphins was more likely to be found in shallow southern areas and the other in the deeper waters of the central northern bay. In further analysis removing admixed individuals, southern dolphins appeared genetically restricted with lower levels of variation (AR = 3.252, π = 0.003) and high mean relatedness (r = 0.239) between individuals. In contrast, northern dolphins were more diverse (AR = 4.850, π = 0.009) and were mixing with a group of dolphins outside the bay (microsatellite‐based STRUCTURE analysis), which appears to have historically been distinct from the bay dolphins (mtDNA Φ_ST = 0.272, P < 0.001). This study demonstrates the ability of genetic techniques to expose fine‐scale patterns of population structure and explore their origins and mechanisms. A complex variety of inter‐related factors including local habitat variation, differential resource use, social behaviour and learning, and anthropogenic disturbances are likely to have played a role in driving fine‐scale population structure among bottlenose dolphins in Moreton Bay.     

Morphometric and genetic differentiation among populations of flat‐headed cusimanse (Crossarchus platycephalus) in Nigeria

Geographic barriers can partition genetic diversity among populations and drive evolutionary divergence between populations, promoting the speciation process and affecting conservation goals. We integrated morphological and genomic data to assess the distribution of variation in the flat‐headed cusimanse (Crossarchus platycephalus), a species of least conservation concern, on either side of the River Niger in Nigeria. Ecological disturbances affect the conservation status of many other animals in this region. The two populations were differentiated in the snout and fore limbs, with greater morphological diversity in the western population. We used Restriction site Associated DNA sequencing (RAD‐seq) and identified two genotypic clusters in a STRUCTURE analysis. Individuals from the eastern population are almost entirely assigned to one cluster, whereas genotypes from the western population are a mixture of the two clusters. The population from west of the River Niger also had higher heterozygosity. The morphological and population genetic data are therefore in agreement that the population from west of the River Niger is more diverse than the eastern population, and the eastern population contains a subset of the genetic variation found in the western population. Our results demonstrate that combining morphological and genotypic measures of diversity can provide a congruent picture of the distribution of intraspecific variation. The results also suggest that future work should explore the role of the River Niger as a natural barrier to migration in Nigeria.