The impact of global climate change on genetic diversity within populations and species

Ecotypes and evolutionary significant units in endangered North African gazelles

Conservation planning of threatened taxa relies upon accurate data on systematics, ecological traits and suitable habitats. The genus Gazella includes taxa with distinct morphologies and ecological traits, but close phylogenetic relationships. The North African Gazella cuvieri and Gazella leptoceros loderi share morphological and physiological characters but the former is darker and found in mountain areas, while the latter is lighter and associated with sand dunes. Here we aim to assess the genetic distinctiveness of these taxa, to characterize their ecological niches and to identify potential occurrence areas, by analysing 327 samples across North-West Africa. Phylogenetic analyses based on mitochondrial (CYTB) and five nuclear gene fragments (KCAS, LAC, SPTBN1, PRKCI and THYR) show that both taxa comprise a single monophyletic group. However, ecological niche-based modelling suggests that populations of these taxa occupy distinct geographic areas and specific environments. Predicted areas of sympatry were restricted, as a consequence of local sharp transitions in climatic traits. The lack of genetic differentiation between these taxa suggests they should be lumped into G. cuvieri, while ecological and morphological differences indicate they correspond to distinct ecotypes. Conservation planning of G. cuvieri should consider the preservation of both mountain and lowland ecotypes to maintain the overall adaptive potential of the species. This integrative approach provides valuable insights in identifying evolutionary units and should be extended to other gazelles.

Population genetic structure and hybridization patterns in the cryptic sister species Chironomus riparius and Chironomus piger across differentially polluted freshwater systems

Chironomids are an integral and functionally important part of many freshwater ecosystems. Yet, to date, there is limited understanding of their microevolutionary processes under chemically polluted natural environments. In this study, we investigated the genetic variation within populations of the ecotoxicological model species Chironomus riparius and its cryptic sister species Chironomus piger at 18 metal-contaminated and reference sites in northwestern Portugal. Microsatellite analysis was conducted on 909 samples to answer if metal contamination affects genetic variation in natural chironomid populations as previously suggested from controlled laboratory experiments. Similarly high levels of genetic diversity and significant but weak genetic substructuring were found across all sites and temporal replicates, with no effects of metal contamination on the genetic variation or species' abundance, although C. piger tended to be less frequent at highly contaminated sites. Our results indicate that high levels of gene flow and population dynamic processes may overlay potential pollutant effects. At least for our study species, we conclude that the "genetic erosion hypothesis", which suggests that chemical pollution will reduce genome-wide genetic variability in affected populations, does not hold under natural conditions. Interestingly, our study provides evidence of successful hybridization between the two sister species under natural conditions.

Climate-induced glacier and snow loss imperils alpine stream insects

Climate warming is causing rapid loss of glaciers and snowpack in mountainous regions worldwide. These changes are predicted to negatively impact the habitats of many range-restricted species, particularly endemic, mountaintop species dependent on the unique thermal and hydrologic conditions found only in glacier-fed and snow melt-driven alpine streams. Although progress has been made, existing understanding of the status, distribution, and ecology of alpine aquatic species, particularly in North America, is lacking, thereby hindering conservation and management programs. Two aquatic insects - the meltwater stonefly (Lednia tumana) and the glacier stonefly (Zapada glacier) - were recently proposed for listing under the U.S. Endangered Species Act due to climate-change-induced habitat loss. Using a large dataset (272 streams, 482 total sites) with high-resolution climate and habitat information, we describe the distribution, status, and key environmental features that limit L. tumana and Z. glacier across the northern Rocky Mountains. Lednia tumana was detected in 113 streams (175 sites) within Glacier National Park (GNP) and surrounding areas. The probability of L. tumana occurrence increased with cold stream temperatures and close proximity to glaciers and permanent snowfields. Similarly, densities of L. tumana declined with increasing distance from stream source. Zapada glacier was only detected in 10 streams (24 sites), six in GNP and four in mountain ranges up to ~600 km southwest. Our results show that both L. tumana and Z. glacier inhabit an extremely narrow distribution, restricted to short sections of cold, alpine streams often below glaciers predicted to disappear over the next two decades. Climate warming-induced glacier and snow loss clearly imperils the persistence of L. tumana and Z. glacier throughout their ranges, highlighting the role of mountaintop aquatic invertebrates as sentinels of climate change in mid-latitude regions.

Bioclimatic transect networks: Powerful observatories of ecological change

Transects that traverse substantial climate gradients are important tools for climate change research and allow questions on the extent to which phenotypic variation associates with climate, the link between climate and species distributions, and variation in sensitivity to climate change among biomes to be addressed. However, the potential limitations of individual transect studies have recently been highlighted. Here, we argue that replicating and networking transects, along with the introduction of experimental treatments, addresses these concerns. Transect networks provide cost‐effective and robust insights into ecological and evolutionary adaptation and improve forecasting of ecosystem change. We draw on the experience and research facilitated by the Australian Transect Network to demonstrate our case, with examples, to clarify how population‐ and community‐level studies can be integrated with observations from multiple transects, manipulative experiments, genomics, and ecological modeling to gain novel insights into how species and systems respond to climate change. This integration can provide a spatiotemporal understanding of past and future climate‐induced changes, which will inform effective management actions for promoting biodiversity resilience.

Evolutionary consequences of historical metal contamination for natural populations of Chironomus riparius (Diptera: Chironomidae)

Populations inhabiting metal-impacted freshwater systems located nearby industrial and urban areas may be under intense selection. The present study aims to address two fundamental microevolutionary aspects of metal contamination in the midge Chironomus riparius (Meigen): Are populations inhabiting historically metal contaminated sites genetically adapted to metals? And, are populations from these sites genetically eroded? To answer these questions, C. riparius populations were sampled from three sites with well-known histories of metal contamination and three nearby-located references. Genetic adaptation to metals was investigated through acute and chronic exposures to cadmium (Cd), after rearing all populations for at least six generations under laboratory clean conditions. Genetic diversity was estimated based on the allelic variation of seven microsatellite markers. Results showed higher acute tolerance to Cd in populations originating from metal contaminated sites compared to their respective references and significant differences in two out of three pairwise comparisons. However, there was a mismatch between acute and chronic tolerance to Cd with results of the partial life-cycle tests suggesting fitness costs under control clean conditions in two metal-adapted populations. Despite no evidences of genetic erosion in populations sampled from metal contaminated sites, our results suggest genetically inherited tolerance to Cd in populations inhabiting historically contaminated sites. These findings lend support to the use of C. riparius as a model organism in evolutionary toxicology and highlight the importance of coupling measures of neutral genetic diversity with assessments of chemical tolerance of populations for a better understanding of contaminant-induced adaptation and evolutionary processes.

Adaptive differentiation coincides with local bioclimatic conditions along an elevational cline in populations of a lichen-forming fungus

Background

Many fungal species occur across a variety of habitats. Particularly lichens, fungi forming symbioses with photosynthetic partners, have evolved remarkable tolerances for environmental extremes. Despite their ecological importance and ubiquity, little is known about the genetic basis of adaption in lichen populations. Here we studied patterns of genome-wide differentiation in the lichen-forming fungus Lasallia pustulata along an altitudinal gradient in the Mediterranean region. We resequenced six populations as pools and identified highly differentiated genomic regions. We then detected gene-environment correlations while controlling for shared population history and pooled sequencing bias, and performed ecophysiological experiments to assess fitness differences of individuals from different environments.

Results

We detected two strongly differentiated genetic clusters linked to Mediterranean and temperate-oceanic climate, and an admixture zone, which coincided with the transition between the two bioclimates. High altitude individuals showed ecophysiological adaptations to wetter and more shaded conditions. Highly differentiated genome regions contained a number of genes associated with stress response, local environmental adaptation, and sexual reproduction.

Conclusions

Taken together our results provide evidence for a complex interplay between demographic history and spatially varying selection acting on a number of key biological processes, suggesting a scenario of ecological speciation.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-017-0929-8) contains supplementary material, which is available to authorized users.

Population genomics reveals high gene flow in grass rockfish (Sebastes rastrelliger)

Variation at the genomic level is important for understanding a species' demographic and phylogenetic history, in addition to identifying candidate regions that may be under selection and important to local adaptation. We used restriction-site associated DNA (RAD) sequencing to sample single nucleotide polymorphisms (SNPs) from the grass rockfish genome (Sebastes rastrelliger) to assess range-wide patterns of population structure and test for signatures of selection (i.e. outlier loci). Two different filtering criteria yielded 6572 and 10,393 SNPs that were genotyped in over 70% of 110 individuals from six sampling areas using RAD loci. Global F_ST across sampling sites is 0.0009-0.002 for all loci. Two methods for the detection of outlier loci failed to identify any outlier loci. Clustering analyses and discriminant analysis of principal components were not able to detect any population structure in our dataset. These results suggest high levels of gene flow for this species across its range and that this level of gene flow makes it difficult to detect outlier loci in the current dataset.

Understanding and monitoring the consequences of human impacts on intraspecific variation

Intraspecific variation is a major component of biodiversity, yet it has received relatively little attention from governmental and nongovernmental organizations, especially with regard to conservation plans and the management of wild species. This omission is ill-advised because phenotypic and genetic variations within and among populations can have dramatic effects on ecological and evolutionary processes, including responses to environmental change, the maintenance of species diversity, and ecological stability and resilience. At the same time, environmental changes associated with many human activities, such as land use and climate change, have dramatic and often negative impacts on intraspecific variation. We argue for the need for local, regional, and global programs to monitor intraspecific genetic variation. We suggest that such monitoring should include two main strategies: (i) intensive monitoring of multiple types of genetic variation in selected species and (ii) broad-brush modeling for representative species for predicting changes in variation as a function of changes in population size and range extent. Overall, we call for collaborative efforts to initiate the urgently needed monitoring of intraspecific variation.

Accounting for adaptive capacity and uncertainty in assessments of species' climate-change vulnerability

Climate-change vulnerability assessments (CCVAs) are valuable tools for assessing species' vulnerability to climatic changes, yet failure to include measures of adaptive capacity and to account for sources of uncertainty may limit their effectiveness. We took a more comprehensive approach that incorporates exposure, sensitivity, and capacity to adapt to climate change. We applied our approach to anadromous steelhead trout (Oncorhynchus mykiss) and nonanadromous bull trout (Salvelinus confluentus), threatened salmonids within the Columbia River Basin (U.S.A.). We quantified exposure on the basis of scenarios of future stream temperature and flow, and we represented sensitivity and capacity to adapt to climate change with metrics of habitat quality, demographic condition, and genetic diversity. Both species were found to be highly vulnerable to climate change at low elevations and in their southernmost habitats. However, vulnerability rankings varied widely depending on the factors (climate, habitat, demographic, and genetic) included in the CCVA and often differed for the 2 species at locations where they were sympatric. Our findings illustrate that CCVA results are highly sensitive to data inputs and that spatial differences can complicate multispecies conservation. Based on our results, we suggest that CCVAs be considered within a broader conceptual and computational framework and be used to refine hypotheses, guide research, and compare plausible scenarios of species' vulnerability to climate change.

Climate change and alpine stream biology: progress, challenges, and opportunities for the future

In alpine regions worldwide, climate change is dramatically altering ecosystems and affecting biodiversity in many ways. For streams, receding alpine glaciers and snowfields, paired with altered precipitation regimes, are driving shifts in hydrology, species distributions, basal resources, and threatening the very existence of some habitats and biota. Alpine streams harbour substantial species and genetic diversity due to significant habitat insularity and environmental heterogeneity. Climate change is expected to affect alpine stream biodiversity across many levels of biological resolution from micro- to macroscopic organisms and genes to communities. Herein, we describe the current state of alpine stream biology from an organism-focused perspective. We begin by reviewing seven standard and emerging approaches that combine to form the current state of the discipline. We follow with a call for increased synthesis across existing approaches to improve understanding of how these imperiled ecosystems are responding to rapid environmental change. We then take a forward-looking viewpoint on how alpine stream biologists can make better use of existing data sets through temporal comparisons, integrate remote sensing and geographic information system (GIS) technologies, and apply genomic tools to refine knowledge of underlying evolutionary processes. We conclude with comments about the future of biodiversity conservation in alpine streams to confront the daunting challenge of mitigating the effects of rapid environmental change in these sentinel ecosystems.

Effects of Pleistocene sea-level fluctuations on mangrove population dynamics: a lesson from Sonneratia alba

Background

A large-scale systematical investigation of the influence of Pleistocene climate oscillation on mangrove population dynamics could enrich our knowledge about the evolutionary history during times of historical climate change, which in turn may provide important information for their conservation.

Results

In this study, phylogeography of a mangrove tree Sonneratia alba was studied by sequencing three chloroplast fragments and seven nuclear genes. A low level of genetic diversity at the population level was detected across its range, especially at the range margins, which was mainly attributed to the steep sea-level drop and associated climate fluctuations during the Pleistocene glacial periods. Extremely small effective population size (Ne) was inferred in populations from both eastern and western Malay Peninsula (44 and 396, respectively), mirroring the fragility of mangrove plants and their paucity of robustness against future climate perturbations and human activity. Two major genetic lineages of high divergence were identified in the two mangrove biodiversity centres: the Indo-Malesia and Australasia regions. The estimated splitting time between these two lineages was 3.153 million year ago (MYA), suggesting a role for pre-Pleistocene events in shaping the major diversity patterns of mangrove species. Within the Indo-Malesia region, a subdivision was implicated between the South China Sea (SCS) and the remaining area with a divergence time of 1.874 MYA, corresponding to glacial vicariance when the emerged Sunda Shelf halted genetic exchange between the western and eastern coasts of the Malay Peninsula during Pleistocene sea-level drops. Notably, genetic admixture was observed in populations at the boundary regions, especially in the two populations near the Malacca Strait, indicating secondary contact between divergent lineages during interglacial periods. These interregional genetic exchanges provided ample opportunity for the re-use of standing genetic variation, which could facilitate mangrove establishment and adaptation in new habitats, especially in the context of global climate changes.

Conclusion

Phylogeogrpahic analysis in this study reveal that Pleistocene sea-level fluctuations had profound influence on population differentiation of the mangrove tree S. alba. Our study highlights the fragility of mangrove plants and offers a guide for the conservation of coastal mangrove communities experiencing ongoing changes in sea-level.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-016-0849-z) contains supplementary material, which is available to authorized users.

The role of genetic diversity and past-history selection pressures in the susceptibility of Chironomus riparius populations to environmental stress

Natural populations experiencing intense selection and genetic drift may exhibit limited potential to adapt to environmental change. The present study addresses the following aspects of the "genetic erosion" hypothesis in the midge Chironomus riparius: does long-term mercury (Hg) contamination affect the Hg tolerance of midge populations inhabiting such impacted areas? If so, is there any fitness cost under changing environmental conditions? And does genetic impoverishment influence the susceptibility of C. riparius to cope with environmental stressful conditions? For this end, we tested the acute and chronic tolerance to Hg and salinity in four C. riparius populations differing in their levels of genetic diversity (assessed through microsatellite markers) and past-histories of Hg exposure. Results showed that the midge population collected from a heavily Hg-contaminated site had higher Hg tolerance compared to the population collected from a closely-located reference site suggesting directional selection for Hg-tolerant traits in its native environment despite no genetic erosion in the field. No increased susceptibility under changing environmental conditions of salinity stress was observed. Moreover, results also showed that populations with higher genetic diversity performed better in the partial life-cycle assays providing evidence on the key role that genetic diversity plays as mediator of populations' susceptibility to environmental stress. Our findings are discussed in terms of the suitability of C. riparius as a model organism in evolutionary toxicology studies as well as the validity of ecotoxicological assessments using genetically eroded laboratory populations.

Genetically informed ecological niche models improve climate change predictions

We examined the hypothesis that ecological niche models (ENMs) more accurately predict species distributions when they incorporate information on population genetic structure, and concomitantly, local adaptation. Local adaptation is common in species that span a range of environmental gradients (e.g., soils and climate). Moreover, common garden studies have demonstrated a covariance between neutral markers and functional traits associated with a species' ability to adapt to environmental change. We therefore predicted that genetically distinct populations would respond differently to climate change, resulting in predicted distributions with little overlap. To test whether genetic information improves our ability to predict a species' niche space, we created genetically informed ecological niche models (gENMs) using Populus fremontii (Salicaceae), a widespread tree species in which prior common garden experiments demonstrate strong evidence for local adaptation. Four major findings emerged: (i) gENMs predicted population occurrences with up to 12-fold greater accuracy than models without genetic information; (ii) tests of niche similarity revealed that three ecotypes, identified on the basis of neutral genetic markers and locally adapted populations, are associated with differences in climate; (iii) our forecasts indicate that ongoing climate change will likely shift these ecotypes further apart in geographic space, resulting in greater niche divergence; (iv) ecotypes that currently exhibit the largest geographic distribution and niche breadth appear to be buffered the most from climate change. As diverse agents of selection shape genetic variability and structure within species, we argue that gENMs will lead to more accurate predictions of species distributions under climate change.

Limited Dispersal and Significant Fine - Scale Genetic Structure in a Tropical Montane Parrot Species

Tropical montane ecosystems are biodiversity hotspots harbouring many endemics that are confined to specific habitat types within narrow altitudinal ranges. While deforestation put these ecosystems under threat, we still lack knowledge about how heterogeneous environments like the montane tropics promote population connectivity and persistence. We investigated the fine-scale genetic structure of the two largest subpopulations of the endangered El Oro parakeet (Pyrrhura orcesi) endemic to the Ecuadorian Andes. Specifically, we assessed the genetic divergence between three sites separated by small geographic distances but characterized by a heterogeneous habitat structure. Although geographical distances between sites are small (3-17 km), we found genetic differentiation between all sites. Even though dispersal capacity is generally high in parrots, our findings indicate that dispersal is limited even on this small geographic scale. Individual genotype assignment revealed similar genetic divergence across a valley (~ 3 km distance) compared to a continuous mountain range (~ 13 km distance). Our findings suggest that geographic barriers promote genetic divergence even on small spatial scales in this endangered endemic species. These results may have important implications for many other threatened and endemic species, particularly given the upslope shift of species predicted from climate change.

Advances using molecular data in insect systematics

The size of molecular datasets has been growing exponentially since the mid 1980s, and new technologies have now dramatically increased the slope of this increase. New datasets include genomes, transcriptomes, and hybrid capture data, producing hundreds or thousands of loci. With these datasets, we are approaching a consensus on the higher level insect phylogeny. Huge datasets can produce new challenges in interpreting branch support, and new opportunities in developing better models and more sophisticated partitioning schemes. Dating analyses are improving as we recognize the importance of careful fossil calibration selection. With thousands of genes now available, coalescent methods have come of age. Barcode libraries continue to expand, and new methods are being developed for incorporating them into phylogenies with tens of thousands of individuals.

Phylogeographic Insights into a Peripheral Refugium: The Importance of Cumulative Effect of Glaciation on the Genetic Structure of Two Endemic Plants

Quaternary glaciations and mostly last glacial maximum have shaped the contemporary distribution of many species in the Alps. However, in the Maritime and Ligurian Alps a more complex picture is suggested by the presence of many Tertiary paleoendemisms and by the divergence time between lineages in one endemic species predating the Late Pleistocene glaciation. The low number of endemic species studied limits the understanding of the processes that took place within this region. We used species distribution models and phylogeographical methods to infer glacial refugia and to reconstruct the phylogeographical pattern of Silene cordifolia All. and Viola argenteria Moraldo & Forneris. The predicted suitable area for last glacial maximum roughly fitted current known distribution. Our results suggest that separation of the major clades predates the last glacial maximum and the following repeated glacial and interglacial periods probably drove differentiations. The complex phylogeographical pattern observed in the study species suggests that both populations and genotypes extinction was minimal during the last glacial maximum, probably due to the low impact of glaciations and to topographic complexity in this area. This study underlines the importance of cumulative effect of previous glacial cycles in shaping the genetic structure of plant species in Maritime and Ligurian Alps, as expected for a Mediterranean mountain region more than for an Alpine region.

Population genetic structure of cotton pink bollworm, Pectinophora gossypiella (Saunders) (Lepidoptera: Gelechiidae) using mitochondrial cytochrome oxidase I (COI) gene sequences from India

Pink bollworm (PBW), Pectinophora gossypiella is one of the most destructive pest's globally inflicting huge economic losses in cotton even during later stages of crop growth. In the present investigation, the population genetic structure, distribution, and genetic diversity of P. gossypiella in cotton growing zones of India using partial mitochondrial DNA cytochrome oxidase-I (COI) gene was addressed. The overall haplotype (Hd), number of nucleotide differences (K), and nucleotide diversity (π) were 0.3028, 0.327, and 0.00047, respectively which suggest that entire population exhibited low level of genetic diversity. Zone-wise clustering of population revealed that central zone recorded low level of Hd (0.2730) as compared to north (0.3619) and south (0.3028) zones. The most common haplotype (H1) reported in all 19 locations could be proposed as ancestral/original haplotype. This haplotype with one mutational step formed star-like phylogeny connected with 11 other haplotypes. The phylogenetic relationship studies revealed that most haplotypes of populations are closely related to each other. Haplotype 5 was exclusively present in Dharwad (South zone) shared with populations of Hanumangarh and Bathinda (North zone). The result indicated that there is no isolation by distance effect among the Indian populations of PBW. The present study reports a low genetic diversity among PBW populations of India and H1, as ancestral haplotype from which other haplotypes have evolved suggests that the migration and dispersal over long distance and invasiveness are major factors.

Multilocus Analyses Reveal Postglacial Demographic Shrinkage of Juniperus morrisonicola (Cupressaceae), a Dominant Alpine Species in Taiwan

Postglacial climate changes alter geographical distributions and diversity of species. Such ongoing changes often force species to migrate along the latitude/altitude. Altitudinal gradients represent assemblage of environmental, especially climatic, variable factors that influence the plant distributions. Global warming that triggered upward migrations has therefore impacted the alpine plants on an island. In this study, we examined the genetic structure of Juniperus morrisonicola, a dominant alpine species in Taiwan, and inferred historical, demographic dynamics based on multilocus analyses. Lower levels of genetic diversity in north indicated that populations at higher latitudes were vulnerable to climate change, possibly related to historical alpine glaciers. Neither organellar DNA nor nuclear genes displayed geographical subdivisions, indicating that populations were likely interconnected before migrating upward to isolated mountain peaks, providing low possibilities of seed/pollen dispersal across mountain ranges. Bayesian skyline plots suggested steady population growth of J. morrisonicola followed by recent demographic contraction. In contrast, most lower-elevation plants experienced recent demographic expansion as a result of global warming. The endemic alpine conifer may have experienced dramatic climate changes over the alternation of glacial and interglacial periods, as indicated by a trend showing decreasing genetic diversity with the altitudinal gradient, plus a fact of upward migration.

Phylogeography of the genus Podococcus (Palmae/Arecaceae) in Central African rain forests: Climate stability predicts unique genetic diversity

The tropical rain forests of Central Africa contain high levels of species diversity. Paleovegetation or biodiversity patterns suggested successive contraction/expansion phases on this rain forest cover during the last glacial maximum (LGM). Consequently, the hypothesis of the existence of refugia e.g. habitat stability that harbored populations during adverse climatic periods has been proposed. Understory species are tightly associated to forest cover and consequently are ideal markers of forest dynamics. Here, we used two central African rain forest understory species of the palm genus, Podococcus, to assess the role of past climate variation on their distribution and genetic diversity. Species distribution modeling in the present and at the LGM was used to estimate areas of climatic stability. Genetic diversity and phylogeography were estimated by sequencing near complete plastomes for over 120 individuals. Areas of climatic stability were mainly located in mountainous areas like the Monts de Cristal and Monts Doudou in Gabon, but also lowland coastal forests in southeast Cameroon and northeast Gabon. Genetic diversity analyses shows a clear North-South structure of genetic diversity within one species. This divide was estimated to have originated some 500,000years ago. We show that, in Central Africa, high and unique genetic diversity is strongly correlated with inferred areas of climatic stability since the LGM. Our results further highlight the importance of coastal lowland rain forests in Central Africa as harboring not only high species diversity but also important high levels of unique genetic diversity. In the context of strong human pressure on coastal land use and destruction, such unique diversity hotspots need to be considered in future conservation planning.

Loss of Genetic Diversity and Increased Subdivision in an Endemic Alpine Stonefly Threatened by Climate Change

Much remains unknown about the genetic status and population connectivity of high-elevation and high-latitude freshwater invertebrates, which often persist near snow and ice masses that are disappearing due to climate change. Here we report on the conservation genetics of the meltwater stonefly Lednia tumana (Ricker) of Montana, USA, a cold-water obligate species. We sequenced 1530 bp of mtDNA from 116 L. tumana individuals representing "historic" (>10 yr old) and 2010 populations. The dominant haplotype was common in both time periods, while the second-most-common haplotype was found only in historic samples, having been lost in the interim. The 2010 populations also showed reduced gene and nucleotide diversity and increased genetic isolation. We found lower genetic diversity in L. tumana compared to two other North American stonefly species, Amphinemura linda (Ricker) and Pteronarcys californica Newport. Our results imply small effective sizes, increased fragmentation, limited gene flow, and loss of genetic variation among contemporary L. tumana populations, which can lead to reduced adaptive capacity and increased extinction risk. This study reinforces concerns that ongoing glacier loss threatens the persistence of L. tumana, and provides baseline data and analysis of how future environmental change could impact populations of similar organisms.

Integrating phylogeography and species distribution models: cryptic distributional responses to past climate change in an endemic rodent from the central Chile hotspot

AIM

Biodiversity losses under the species level may have been severely underestimated in future global climate change scenarios. Therefore, it is important to characterize the diversity units at this level, as well as to understand their ecological responses to climatic forcings. We have chosen an endemic rodent from a highly endangered ecogeographic area as a model to look for distributional responses below the species level: Phyllotis darwini.

LOCATION

The central Chile biodiversity hotspot: This area harbours a high number of endemic species, and it is known to have experienced vegetational displacements between two mountain systems during and after the Last Glacial Maximum.

METHODS

We have characterized cryptic lineages inside P. darwini in a classic phylogeographic approach; those intraspecific lineages were considered as relevant units to construct distribution models at Last Glacial Maximum and at present, as border climatic conditions. Differences in distribution between border conditions for each lineage were interpreted as distributional responses to post-glacial climate change.

RESULTS

The species is composed of two major phylogroups: one of them has a broad distribution mainly across the valley but also in mountain ranges, whereas the other displays a disjunct distribution across both mountain ranges and always above 1500 m. The lineage distribution model under LGM climatic conditions suggests that both lineages were co-distributed in the southern portion of P. darwini's current geographic range, mainly at the valley and at the coast.

MAIN CONCLUSIONS

Present distribution of lineages in P. darwini is the consequence of a cryptic distributional response to climate change after LGM: postglacial northward colonization, with strict altitudinal segregation of both phylogroups.

Contemporary genetic structure and postglacial demographic history of the black scorpionfish, Scorpaena porcus, in the Mediterranean and the Black Seas

Understanding the distribution of genetic diversity in the light of past demographic events linked with climatic shifts will help to forecast evolutionary trajectories of ecosystems within the current context of climate change. In this study, mitochondrial sequences and microsatellite loci were analysed using traditional population genetic approaches together with Bayesian dating and the more recent approximate Bayesian computation scenario testing. The genetic structure and demographic history of a commercial fish, the black scorpionfish, Scorpaena porcus, was investigated throughout the Mediterranean and Black Seas. The results suggest that the species recently underwent population expansions, in both seas, likely concomitant with the warming period following the Last Glacial Maximum, 20 000 years ago. A weak contemporaneous genetic differentiation was identified between the Black Sea and the Mediterranean Sea. However, the genetic diversity was similar for populations of the two seas, suggesting a high number of colonizers entered the Black Sea during the interglacial period and/or the presence of a refugial population in the Black Sea during the glacial period. Finally, within seas, an east/west genetic differentiation in the Adriatic seems to prevail, whereas the Black Sea does not show any structured spatial genetic pattern of its population. Overall, these results suggest that the Black Sea is not that isolated from the Mediterranean, and both seas revealed similar evolutionary patterns related to climate change and changes in sea level.

Crossing regimes of temperature dependence in animal movement

A pressing challenge in ecology is to understand the effects of changing global temperatures on food web structure and dynamics. The stability of these complex ecological networks largely depends on how predator-prey interactions may respond to temperature changes. Because predators and prey rely on their velocities to catch food or avoid being eaten, understanding how temperatures may affect animal movement is central to this quest. Despite our efforts, we still lack a mechanistic understanding of how the effect of temperature on metabolic processes scales up to animal movement and beyond. Here, we merge a biomechanical approach, the Metabolic Theory of Ecology and empirical data to show that animal movement displays multiple regimes of temperature dependence. We also show that crossing these regimes has important consequences for population dynamics and stability, which depend on the parameters controlling predator-prey interactions. We argue that this dependence upon interaction parameters may help explain why experimental work on the temperature dependence of interaction strengths has so far yielded conflicting results. More importantly, these changes in the temperature dependence of animal movement can have consequences that go well beyond ecological interactions and affect, for example, animal communication, mating, sensory detection, and any behavioral modality dependent on the movement of limbs. Finally, by not taking into account the changes in temperature dependence reported here we might not be able to properly forecast the impact of global warming on ecological processes and propose appropriate mitigation action when needed.

Host and Environmental Specificity in Bacterial Communities Associated to Two Highly Invasive Marine Species (Genus Asparagopsis)

As habitats change due to global and local pressures, population resilience, and adaptive processes depend not only on their gene pools but also on their associated bacteria communities. The hologenome can play a determinant role in adaptive evolution of higher organisms that rely on their bacterial associates for vital processes. In this study, we focus on the associated bacteria of the two most invasive seaweeds in southwest Iberia (coastal mainland) and nearby offshore Atlantic islands, Asparagopsis taxiformis and Asparagopsis armata. Bacterial communities were characterized using 16S rRNA barcoding through 454 next generation sequencing and exploratory shotgun metagenomics to provide functional insights and a backbone for future functional studies. The bacterial community composition was clearly different between the two species A. taxiformis and A. armata and between continental and island habitats. The latter was mainly due to higher abundances of Acidimicrobiales, Sphingomonadales, Xanthomonadales, Myxococcales, and Alteromonadales on the continent. Metabolic assignments for these groups contained a higher number of reads in functions related to oxidative stress and resistance to toxic compounds, more precisely heavy metals. These results are in agreement with their usual association with hydrocarbon degradation and heavy-metals detoxification. In contrast, A. taxiformis from islands contained more bacteria related to oligotrophic environments which might putatively play a role in mineralization of dissolved organic matter. The higher number of functional assignments found in the metagenomes of A. taxiformis collected from Cape Verde Islands suggest a higher contribution of bacteria to compensate nutrient limitation in oligotrophic environments. Our results show that Asparagopsis-associated bacterial communities have host-specificity and are modulated by environmental conditions. Whether this environmental effect reflects the host's selective requirements or the locally available bacteria remains to be addressed. However, the known functional capacities of these bacterial communities indicate their potential for eco-physiological functions that could be valuable for the host fitness.

Distribution and population genetic variation of cryptic species of the Alpine mayfly Baetis alpinus (Ephemeroptera: Baetidae) in the Central Alps

BACKGROUND

Many species contain evolutionarily distinct groups that are genetically highly differentiated but morphologically difficult to distinguish (i.e., cryptic species). The presence of cryptic species poses significant challenges for the accurate assessment of biodiversity and, if unrecognized, may lead to erroneous inferences in many fields of biological research and conservation.

RESULTS

We tested for cryptic genetic variation within the broadly distributed alpine mayfly Baetis alpinus across several major European drainages in the central Alps. Bayesian clustering and multivariate analyses of nuclear microsatellite loci, combined with phylogenetic analyses of mitochondrial DNA, were used to assess population genetic structure and diversity. We identified two genetically highly differentiated lineages (A and B) that had no obvious differences in regional distribution patterns, and occurred in local sympatry. Furthermore, the two lineages differed in relative abundance, overall levels of genetic diversity as well as patterns of population structure: lineage A was abundant, widely distributed and had a higher level of genetic variation, whereas lineage B was less abundant, more prevalent in spring-fed tributaries than glacier-fed streams and restricted to high elevations. Subsequent morphological analyses revealed that traits previously acknowledged as intraspecific variation of B. alpinus in fact segregated these two lineages.

CONCLUSIONS

Taken together, our findings indicate that even common and apparently ecologically well-studied species may consist of reproductively isolated units, with distinct evolutionary histories and likely different ecology and evolutionary potential. These findings emphasize the need to investigate hidden diversity even in well-known species to allow for appropriate assessment of biological diversity and conservation measures.

Buffering and Amplifying Interactions among OAW (Ocean Acidification & Warming) and Nutrient Enrichment on Early Life-Stage Fucus vesiculosus L. (Phaeophyceae) and Their Carry Over Effects to Hypoxia Impact

Ocean acidification and warming (OAW) are occurring globally. Additionally, at a more local scale the spreading of hypoxic conditions is promoted by eutrophication and warming. In the semi-enclosed brackish Baltic Sea, occasional upwelling in late summer and autumn may expose even shallow-water communities including the macroalga Fucus vesiculosus to particularly acidified, nutrient-rich and oxygen-poor water bodies. During summer 2014 (July-September) sibling groups of early life-stage F. vesiculosus were exposed to OAW in the presence and absence of enhanced nutrient levels and, subsequently to a single upwelling event in a near-natural scenario which included all environmental fluctuations in the Kiel Fjord, southwestern Baltic Sea, Germany (54°27 ´N, 10°11 ´W). We strove to elucidate the single and combined impacts of these potential stressors, and how stress sensitivity varies among genetically different sibling groups. Enhanced by a circumstantial natural heat wave, warming and acidification increased mortalities and reduced growth in F. vesiculosus germlings. This impact, however, was mitigated by enhanced nutrient conditions. Survival under OAW conditions strongly varied among sibling groups hinting at a substantial adaptive potential of the natural Fucus populations in the Western Baltic. A three-day experimental upwelling caused severe mortality of Fucus germlings, which was substantially more severe in those sibling groups which previously had been exposed to OAW. Our results show that global (OAW), regional (nutrient enrichment) and local pressures (upwelling), both alone and co-occurring may have synergistic and antagonistic effects on survival and/or growth of Fucus germlings. This result emphasizes the need to consider combined stress effects.

New view of population genetics of zooplankton: RAD-seq analysis reveals population structure of the North Atlantic planktonic copepod Centropages typicus

Detection of population genetic structure of zooplankton at medium-to-small spatial scales in the absence of physical barriers has remained challenging and controversial. The large population sizes and high rates of gene flow characteristic of zooplankton have made resolution of geographical differentiation very difficult, especially when using few genetic markers and assuming equilibrium conditions. Next-generation sequencing now allows simultaneous sampling of hundreds to thousands of genetic markers; new analytical approaches allow studies under nonequilibrium conditions and directional migration. Samples of the North Atlantic Ocean planktonic copepod, Centropages typicus, were analysed using restriction site-associated DNA (RAD) sequencing on a PROTON platform. Although prior studies revealed no genetic differentiation of populations across the geographical range of the species, analysis of RAD tags showed significant structure across the North Atlantic Ocean. We also compared the likelihood for models of connectivity among NW Atlantic populations under various directional flow scenarios that replicate oceanographic conditions of the sampled domain. High-density marker sampling with RAD sequencing markedly outperformed other technical and analytical approaches in detection of population genetic structure and characterization of connectivity of this high geneflow zooplankton species.

Drastic underestimation of amphipod biodiversity in the endangered Irano-Anatolian and Caucasus biodiversity hotspots

Biodiversity hotspots are centers of biological diversity and particularly threatened by anthropogenic activities. Their true magnitude of species diversity and endemism, however, is still largely unknown as species diversity is traditionally assessed using morphological descriptions only, thereby ignoring cryptic species. This directly limits evidence-based monitoring and management strategies. Here we used molecular species delimitation methods to quantify cryptic diversity of the montane amphipods in the Irano-Anatolian and Caucasus biodiversity hotspots. Amphipods are ecosystem engineers in rivers and lakes. Species diversity was assessed by analysing two genetic markers (mitochondrial COI and nuclear 28S rDNA), compared with morphological assignments. Our results unambiguously demonstrate that species diversity and endemism is dramatically underestimated, with 42 genetically identified freshwater species in only five reported morphospecies. Over 90% of the newly recovered species cluster inside Gammarus komareki and G. lacustris; 69% of the recovered species comprise narrow range endemics. Amphipod biodiversity is drastically underestimated for the studied regions. Thus, the risk of biodiversity loss is significantly greater than currently inferred as most endangered species remain unrecognized and/or are only found locally. Integrative application of genetic assessments in monitoring programs will help to understand the true magnitude of biodiversity and accurately evaluate its threat status.

Tackling intraspecific genetic structure in distribution models better reflects species geographical range

Genetic diversity provides insight into heterogeneous demographic and adaptive history across organisms' distribution ranges. For this reason, decomposing single species into genetic units may represent a powerful tool to better understand biogeographical patterns as well as improve predictions of the effects of GCC (global climate change) on biodiversity loss. Using 279 georeferenced Iberian accessions, we used classes of three intraspecific genetic units of the annual plant Arabidopsis thaliana obtained from the genetic analyses of nuclear SNPs (single nucleotide polymorphisms), chloroplast SNPs, and the vernalization requirement for flowering. We used SDM (species distribution models), including climate, vegetation, and soil data, at the whole-species and genetic-unit levels. We compared model outputs for present environmental conditions and with a particularly severe GCC scenario. SDM accuracy was high for genetic units with smaller distribution ranges. Kernel density plots identified the environmental variables underpinning potential distribution ranges of genetic units. Combinations of environmental variables accounted for potential distribution ranges of genetic units, which shrank dramatically with GCC at almost all levels. Only two genetic clusters increased their potential distribution ranges with GCC. The application of SDM to intraspecific genetic units provides a detailed picture on the biogeographical patterns of distinct genetic groups based on different genetic criteria. Our approach also allowed us to pinpoint the genetic changes, in terms of genetic background and physiological requirements for flowering, that Iberian A. thaliana may experience with a GCC scenario applying SDM to intraspecific genetic units.

The eastern migratory caribou: the role of genetic introgression in ecotype evolution

Understanding the evolutionary history of contemporary animal groups is essential for conservation and management of endangered species like caribou (Rangifer tarandus). In central Canada, the ranges of two caribou subspecies (barren-ground/woodland caribou) and two woodland caribou ecotypes (boreal/eastern migratory) overlap. Our objectives were to reconstruct the evolutionary history of the eastern migratory ecotype and to assess the potential role of introgression in ecotype evolution. STRUCTURE analyses identified five higher order groups (i.e. three boreal caribou populations, eastern migratory ecotype and barren-ground). The evolutionary history of the eastern migratory ecotype was best explained by an early genetic introgression from barren-ground into a woodland caribou lineage during the Late Pleistocene and subsequent divergence of the eastern migratory ecotype during the Holocene. These results are consistent with the retreat of the Laurentide ice sheet and the colonization of the Hudson Bay coastal areas subsequent to the establishment of forest tundra vegetation approximately 7000 years ago. This historical reconstruction of the eastern migratory ecotype further supports its current classification as a conservation unit, specifically a Designatable Unit, under Canada's Species at Risk Act. These findings have implications for other sub-specific contact zones for caribou and other North American species in conservation unit delineation.

Adaptive introgression as a resource for management and genetic conservation in a changing climate

Current rates of climate change require organisms to respond through migration, phenotypic plasticity, or genetic changes via adaptation. We focused on questions regarding species' and populations' ability to respond to climate change through adaptation. Specifically, the role adaptive introgression, movement of genetic material from the genome of 1 species into the genome of another through repeated interbreeding, may play in increasing species' ability to respond to a changing climate. Such interspecific gene flow may mediate extinction risk or consequences of limited adaptive potential that result from standing genetic variation and mutation alone, enabling a quicker demographic recovery in response to changing environments. Despite the near dismissal of the potential benefits of hybridization by conservation practitioners, we examined a number of case studies across different taxa that suggest gene flow between sympatric or parapatric sister species or within species that exhibit strong ecotypic differentiation may represent an underutilized management option to conserve evolutionary potential in a changing environment. This will be particularly true where advanced-generation hybrids exhibit adaptive traits outside the parental phenotypic range, a phenomenon known as transgressive segregation. The ideas presented in this essay are meant to provoke discussion regarding how we maintain evolutionary potential, the conservation value of natural hybrid zones, and consideration of their important role in adaptation to climate.

Diverse seed banks favour adaptation of microalgal populations to future climate conditions

Selection of suitable genotypes from diverse seed banks may help phytoplankton populations to cope with environmental changes. This study examines whether the high genotypic diversity found in the Baltic cyst pool of the toxic dinoflagellate Alexandrium ostenfeldii is coupled to phenotypic variability that could aid short-term adaptation. Growth rates, cellular toxicities and bioluminescence of 34 genetically different clones isolated from cyst beds of four Baltic bloom sites were determined in batch culture experiments along temperature and salinity gradients covering present and future conditions in the Baltic Sea. For all parameters a significant effect of genotype on the response to temperature and salinity changes was identified. General or site-specific effects of the two factors remained minor. Clones thriving at future conditions were different from the best performing at present conditions, suggesting that genotypic shifts may be expected in the future. Increased proportions of highly potent saxitoxin were observed as a plastic response to temperature increase, indicating a potential for higher toxicity of future blooms. The observed standing variation in Baltic seed banks of A. ostenfeldii suggests that the population is likely to persist under environmental change.

Genetic differentiation of the regional Plutella xylostella populations across the Taiwan Strait based on identification of microsatellite markers

Movement of individuals through events, such as storms or crop transportation, may affect survival and distribution of insect pests, as well as population genetic structure at a regional scale. Understanding what factors contribute to gene flow in pest populations remains very important for sustainable pest management. The diamondback moth (Plutella xylostella) is an insect pest well known for its capacity of moving over short to long distances. Here, we used newly isolated microsatellite markers to analyze the genetic structure of nine populations across the Taiwan Strait of China (Taiwan and Fujian). A total of 12,152 simple sequence repeats (SSRs) were initially identified from the P. xylostella transcriptome (~94 Mb), with an average of 129 SSRs per Mb. Nine SSRs were validated to be polymorphic markers, and eight were used for this population genetic study. Our results showed that the P. xylostella populations could be divided into distinct two clusters, which is likely due to the year-round airflows in this region. A pattern of isolation by distance among the local populations within Fujian was found, and may be related to vegetable transportation. Considering the complexity of the P. xylostella population genetic structure from local and regional to global levels, we propose that developing ecologically sound strategies for managing this pest will require knowledge of the link between behavioral and population ecology and its genetic structure.

Genetic variation of Lymnaea stagnalis tolerance to copper: A test of selection hypotheses and its relevance for ecological risk assessment

The use of standardized monospecific testing to assess the ecological risk of chemicals implicitly relies on the strong assumption that intraspecific variation in sensitivity is negligible or irrelevant in this context. In this study, we investigated genetic variation in copper sensitivity of the freshwater snail Lymnaea stagnalis, using lineages stemming from eight natural populations or strains found to be genetically differentiated at neutral markers. Copper-induced mortality varied widely among populations, as did the estimated daily death rate and time to 50% mortality (LT50). Population genetic divergence in copper sensitivity was compared to neutral differentiation using the QST-FST approach. No evidence for homogenizing selection could be detected. This result demonstrates that species-level extrapolations from single population studies are highly unreliable. The study provides a simple example of how evolutionary principles could be incorporated into ecotoxicity testing in order to refine ecological risk assessment.