Past climate change drives current genetic structure of an endangered freshwater mussel species

Historical‐to‐recent climate change and anthropogenic disturbance affect species distributions and genetic structure. The Rio Grande watershed of the United States and Mexico encompasses ecosystems that are intensively exploited, resulting in substantial degradation of aquatic habitats. While significant anthropogenic disturbances in the Rio Grande are recent, inhospitable conditions for freshwater organisms likely existed prior to such disturbances. A combination of anthropogenic and past climate factors may contribute to current distributions of aquatic fauna in the Rio Grande basin. We used mitochondrial DNA and 18 microsatellite loci to infer evolutionary history and genetic structure of an endangered freshwater mussel, Popenaias popeii, throughout the Rio Grande drainage. We estimated spatial connectivity and gene flow across extant populations of P. popeii and used ecological niche models (ENMs) and approximate Bayesian computation (ABC) to infer its evolutionary history during the Pleistocene. structure results recovered regional and local population clusters in the Rio Grande. ENMs predicted drastic reductions in suitable habitat during the last glacial maximum. ABC analyses suggested that regional population structure likely arose in this species during the mid‐to‐late Pleistocene and was followed by a late Pleistocene population bottleneck in New Mexico populations. The local population structure arose relatively recently, perhaps due to anthropogenic factors. Popenaias popeii, one of the few freshwater mussel species native to the Rio Grande basin, is a case study for understanding how both geological and anthropogenic factors shape current population genetic structure. Conservation strategies for this species should account for the fragmented nature of contemporary populations.     

Colonization from divergent ancestors: glaciation signatures on contemporary patterns of genomic variation in Collared Pikas (Ochotona collaris)

Identifying the genetic structure of a species and the factors that drive it is an important first step in modern population management, in part because populations evolving from separate ancestral sources may possess potentially different characteristics. This is especially true for climate‐sensitive species such as pikas, where the delimitation of distinct genetic units and the characterization of population responses to contemporary and historical environmental pressures are of particular interest. We combined a restriction site‐associated DNA sequencing (RADSeq) data set containing 4156 single nucleotide polymorphisms with ecological niche models (ENMs) of present and past habitat suitability to characterize population composition and evaluate the effects of historical range shifts, contemporary climates and landscape factors on gene flow in Collared Pikas, which are found in Alaska and adjacent regions of northwestern Canada and are the lesser‐studied of North America's two pika species. The results suggest that contemporary environmental factors contribute little to current population connectivity. Instead, genetic diversity is strongly shaped by the presence of three ancestral lineages isolated during the Pleistocene (~148 and 52 kya). Based on ENMs and genetic data, populations originating from a northern refugium experienced longer‐term stability, whereas both southern lineages underwent population expansion – contradicting the southern stability and northern expansion patterns seen in many other taxa. Current populations are comparable with respect to generally low diversity within populations and little‐to‐no recent admixture. The predominance of divergent histories structuring populations implies that if we are to understand and manage pika populations, we must specifically assess and accurately account for the forces underlying genetic similarity.     

Influence of late Quaternary climate change on present patterns of genetic variation in valley oak,Quercus lobata Née

Phylogeography and ecological niche models (ENMs) suggest that late Quaternary glacial cycles have played a prominent role in shaping present population genetic structure and diversity, but have not applied quantitative methods to dissect the relative contribution of past and present climate vs. other forces. We integrate multilocus phylogeography, climate‐based ENMs and multivariate statistical approaches to infer the effects of late Quaternary climate change on contemporary genetic variation of valley oak (Quercus lobata Née). ENMs indicated that valley oak maintained a stable distribution with local migration from the last interglacial period (~120 ka) to the Last Glacial Maximum (~21 ka, LGM) to the present compared with large‐scale range shifts for an eastern North American white oak (Quercus alba L.). Coast Range and Sierra Nevada foothill populations diverged in the late Pleistocene before the LGM [104 ka (28–1622)] and have occupied somewhat distinct climate niches, according to ENMs and coalescent analyses of divergence time. In accordance with neutral expectations for stable populations, nuclear microsatellite diversity positively correlated with niche stability from the LGM to present. Most strikingly, nuclear and chloroplast microsatellite variation significantly correlated with LGM climate, even after controlling for associations with geographic location and present climate using partial redundancy analyses. Variance partitioning showed that LGM climate uniquely explains a similar proportion of genetic variance as present climate (16% vs. 11–18%), and together, past and present climate explains more than geography (19%). Climate can influence local expansion–contraction dynamics, flowering phenology and thus gene flow, and/or impose selective pressures. These results highlight the lingering effect of past climate on genetic variation in species with stable distributions.     

Distance,elevation and environment as drivers of diversity and divergence in bumble bees across latitude and altitude

Identifying drivers of dispersal limitation and genetic differentiation is a key goal in biogeography. We examine patterns of population connectivity and genetic diversity using restriction site‐associated DNA sequencing (RADseq) in two bumble bee species, Bombus vosnesenskii and Bombus bifarius, across latitude and altitude in mountain ranges from California, Oregon and Washington, U.S.A. Bombus vosnesenskii, which occurs across a broader elevational range at most latitudes, exhibits little population structure while B. bifarius, which occupies a relatively narrow higher elevation niche across most latitudes, exhibits much stronger population differentiation, although gene flow in both species is best explained by isolation with environmental niche resistance. A relationship between elevational habitat breadth and genetic diversity is also apparent, with B. vosnesenskii exhibiting relatively consistent levels of genetic diversity across its range, while B. bifarius has reduced genetic diversity at low latitudes, where it is restricted to high‐elevation habitat. The results of this study highlight the importance of the intersect between elevational range and habitat suitability in influencing population connectivity and suggest that future climate warming will have a fragmenting effect even on populations that are presently well connected, as they track their thermal niches upward in montane systems.     

The shaping of genetic variation in edge‐of‐range populations under past and future climate change

With rates of climate change exceeding the rate at which many species are able to shift their range or adapt, it is important to understand how future changes are likely to affect biodiversity at all levels of organisation. Understanding past responses and extent of niche conservatism in climatic tolerance can help predict future consequences. We use an integrated approach to determine the genetic consequences of past and future climate changes on a bat species, Plecotus austriacus. Glacial refugia predicted by palaeo‐modelling match those identified from analyses of extant genetic diversity and model‐based inference of demographic history. Former refugial populations currently contain disproportionately high genetic diversity, but niche conservatism, shifts in suitable areas and barriers to migration mean that these hotspots of genetic diversity are under threat from future climate change. Evidence of population decline despite recent northward migration highlights the need to conserve leading‐edge populations for spearheading future range shifts.     

Noninvasive sampling reveals population genetic structure in the Royle’s pika,Ochotona roylei,in the western Himalaya

Understanding population genetic structure of climate‐sensitive herbivore species is important as it provides useful insights on how shifts in environmental conditions can alter their distribution and abundance. Herbivore responses to the environment can have a strong indirect cascading effect on community structure. This is particularly important for Royle's pika (Lagomorpha: Ochotona roylei), a herbivorous talus‐dwelling species in alpine ecosystem, which forms a major prey base for many carnivores in the Himalayan arc. In this study, we used seven polymorphic microsatellite loci to detect evidence for recent changes in genetic diversity and population structure in Royle's pika across five locations sampled between 8 and 160 km apart in the western Himalaya. Using four clustering approaches, we found the presence of significant contemporary genetic structure in Royle's pika populations. The detected genetic structure could be primarily attributed to the landscape features in alpine habitat (e.g., wide lowland valleys, rivers) that may act as semipermeable barriers to gene flow and distribution of food plants, which are key determinants in spatial distribution of herbivores. Pika showed low inbreeding coefficients (F_IS) and a high level of pairwise relatedness for individuals within 1 km suggesting low dispersal abilities of talus‐dwelling pikas. We have found evidence of a recent population bottleneck, possibly due to effects of environmental disturbances (e.g., snow melting patterns or thermal stress). Our results reveal significant evidence of isolation by distance in genetic differentiation (F_ST range = 0.04–0.19). This is the first population genetics study on Royle's pika, which helps to address evolutionary consequences of climate change which are expected to significantly affect the distribution and population dynamics in this talus‐dwelling species.     

Projected climate change causes loss and redistribution of genetic diversity in a model metapopulation of a medium‐good disperser

Climate change causes species ranges to shift geographically as individuals colonise new suitable temperature zones or fail to reproduce where climate conditions fall below tolerance levels. Little is known about the potential loss of genetic diversity in such dynamic ranges. We investigated the level and distribution of neutral genetic diversity in shifting metapopulations during three scenarios of temperature increase projected for this century and at various degrees of weather variability. We used an individual‐based and spatially explicit metapopulation model in which temperature zones were simulated to move across a fragmented landscape following different climate change scenarios. Although the connectivity between habitat patches allowed the species, modelled after the middle spotted woodpecker Dendrocopos medius, to move along with the shifting temperature range, existing neutral genetic diversity was lost under all three temperature increase scenarios. This was independent of the loss of individuals. The explanation for this effect is that only a part of the original genetic variation moved into the newly colonised habitat. Under increased weather variability the number of individuals and the number of alleles per locus were persistently lower. However, the pattern of changes in allele distributions under temperature zone shifts was the same under all weather variability levels. Genetic differentiation between populations had a tendency to increase at metapopulation range margins, but decreased again when population sizes increased in time. Increased weather variability led to increased variation around the mean genetic differentiation across the metapopulation. Our results illustrate the usefulness of more realistic models for studying the effects of climate change on metapopulations. They indicate that biodiversity monitoring indices based on species occurrence and abundance are not a good proxy for the trend in the level of genetic diversity. Further, the results underline the importance of conserving areas where species have existed for a long time as modern refugia for genetic diversity.     

Good‐bye to tropical alpine plant giants under warmer climates? Loss of range and genetic diversity in Lobelia rhynchopetalum

The main aim of this paper is to address consequences of climate warming on loss of habitat and genetic diversity in the enigmatic tropical alpine giant rosette plants using the Ethiopian endemic Lobelia rhynchopetalum as a model. We modeled the habitat suitability of L. rhynchopetalum and assessed how its range is affected under two climate models and four emission scenarios. We used three statistical algorithms calibrated to represent two different complexity levels of the response. We analyzed genetic diversity using amplified fragment length polymorphisms and assessed the impact of the projected range loss. Under all model and scenario combinations and consistent across algorithms and complexity levels, this afro‐alpine flagship species faces massive range reduction. Only 3.4% of its habitat seems to remain suitable on average by 2,080, resulting in loss of 82% (CI 75%–87%) of its genetic diversity. The remaining suitable habitat is projected to be fragmented among and reduced to four mountain peaks, further deteriorating the probability of long‐term sustainability of viable populations. Because of the similar morphological and physiological traits developed through convergent evolution by tropical alpine giant rosette plants in response to diurnal freeze‐thaw cycles, they most likely respond to climate change in a similar way as our study species. We conclude that specialized high‐alpine giant rosette plants, such as L. rhynchopetalum, are likely to face very high risk of extinction following climate warming.     

Genetic variability and the effect of habitat fragmentation in spotted suslik Spermophilus suslicus populations from two different regions

Endangered species worldwide exist in remnant populations, often within fragmented landscapes. Although assessment of genetic diversity in fragmented habitats is very important for conservation purposes, it is usually impossible to evaluate the amount of diversity that has actually been lost. Here, we compared population structure and levels of genetic diversity within populations of spotted suslik Spermophilus suslicus, inhabiting two different parts of the species range characterized by different levels of habitat connectivity. We used microsatellites to analyze 10 critically endangered populations located at the western part of the range, where suslik habitat have been severely devastated due to agriculture industrialization. Their genetic composition was compared with four populations from the eastern part of the range where the species still occupies habitat with reasonable levels of connectivity. In the western region, we detected extreme population structure (F _ST = 0.20) and levels of genetic diversity (Allelic richness ranged from 1.45 to 3.07) characteristic for highly endangered populations. Alternatively, in the eastern region we found significantly higher allelic richness (from 5.09 to 5.81) and insignificant population structure (F _ST = 0.03). As we identified a strong correlation between genetic and geographic distance and a lack of private alleles in the western region, we conclude that extreme population structure and lower genetic diversity is due to recent habitat loss. Results from this study provide guidelines for conservation and management of this highly endangered species.     

Gene flow connects coastal populations of a habitat specialist,the Clapper Rail Rallus crepitans

Examining population genetic structure can reveal patterns of reproductive isolation or population mixing and inform conservation management. Some avian species are predicted to exhibit minimal genetic differentiation among populations as a result of the species high mobility, with habitat specialists tending to show greater fine‐scale genetic structure. To explore the relationship between habitat specialization and gene flow, we investigated the genetic structure of a saltmarsh specialist with high potential mobility across a wide geographical range of fragmented habitat. Little variation among mitochondrial sequences (620 bp from ND2) was observed among 149 individual Clapper Rails Rallus crepitans sampled along the Atlantic coast of the USA, with the majority of individuals at all sampling sites sharing a single haplotype. Genotyping of nine microsatellite loci across 136 individuals revealed moderate genetic diversity, no evidence of bottlenecks and a weak pattern of genetic differentiation that increased with geographical distance. Multivariate analyses, Bayesian clustering and an AMOVA all suggested a lack of genetic structuring across the Atlantic coast of the USA, with all individuals grouped into a single interbreeding population. Spatial autocorrelation analyses showed evidence of weak female philopatry and a lack of male philopatry. We conclude that high gene flow connecting populations of this habitat specialist may result from the interaction of ecological and behavioural factors that promote dispersal and limit natal philopatry and breeding‐site fidelity. As climate change threatens saltmarshes, the genetic diversity and population connectivity of Clapper Rails may promote resilience of their populations. This finding helps inform about potential fates of other similarly behaving saltmarsh specialists on the Atlantic coast.     

A macroecological approach to evolutionary rescue and adaptation to climate change

Despite the widespread use of ecological niche models (ENMs) for predicting the responses of species to climate change, these models do not explicitly incorporate any population‐level mechanism. On the other hand, mechanistic models adding population processes (e.g. biotic interactions, dispersal and adaptive potential to abiotic conditions) are much more complex and difficult to parameterize, especially if the goal is to predict range shifts for many species simultaneously. In particular, the adaptive potential (based on genetic adaptations, phenotypic plasticity and behavioral adjustments for physiological responses) of local populations has been a less studied mechanism affecting species’ responses to climatic change so far. Here, we discuss and apply an alternative macroecological framework to evaluate the potential role of evolutionary rescue under climate change based on ENMs. We begin by reviewing eco‐evolutionary models that evaluate the maximum sustainable evolutionary rate under a scenario of environmental change, showing how they can be used to understand the impact of temperature change on a Neotropical anuran species, the Schneider's toad Rhinella diptycha. Then we show how to evaluate spatial patterns of species’ geographic range shift using such models, by estimating evolutionary rates at the trailing edge of species distribution estimated by ENMs and by recalculating the relative amount of total range loss under climate change. We show how different models can reduce the expected range loss predicted for the studied species by potential ecophysiological adaptations in some regions of the trailing edge predicted by ENMs. For general applications, we believe that parameters for large numbers of species and populations can be obtained from macroecological generalizations (e.g. allometric equations and ecogeographical rules), so our framework coupling ENMs with eco‐evolutionary models can be applied to achieve a more accurate picture of potential impacts from climate change and other threats to biodiversity.     

Identifying refugia from climate change using coupled ecological and genetic data in a transitional Mediterranean‐temperate tree species

Populations occurring in areas of overlap between the current and future distribution of a species are particularly important because they can represent “refugia from climate change”. We coupled ecological and range‐wide genetic variation data to detect such areas and to evaluate the impacts of habitat suitability changes on the genetic diversity of the transitional Mediterranean‐temperate tree Fraxinus angustifolia. We sampled and genotyped 38 natural populations comprising 1006 individuals from across Europe. We found the highest genetic diversity in western and northern Mediterranean populations, as well as a significant west to east decline in genetic diversity. Areas of potential refugia that correspond to approximately 70% of the suitable habitat may support the persistence of more than 90% of the total number of alleles in the future. Moreover, based on correlations between Bayesian genetic assignment and climate, climate change may favour the westward spread of the Black Sea gene pool in the long term. Overall, our results suggest that the northerly core areas of the current distribution contain the most important part of the genetic variation for this species and may serve as in situ macrorefugia from ongoing climate change. However, rear‐edge populations of the southern Mediterranean may be exposed to a potential loss of unique genetic diversity owing to habitat suitability changes unless populations can persist in microrefugia that have facilitated such persistence in the past.     

Phylogeography,population genetics and distribution modelling reveal vulnerability of Scirpus longii (Cyperaceae) and the Atlantic Coastal Plain Flora to climate change

A proactive approach to conservation must be predictive, anticipating how habitats will change and which species are likely to decline or prosper. We use composite species distribution modelling to identify suitable habitats for 18 members of the North American Atlantic Coastal Plain Flora (ACPF) since the Last Glacial Maximum and project these into the future. We then use Scirpus longii (Cyperaceae), a globally imperiled ACPF sedge with many of the characteristics of extinction vulnerability, as a case study. We integrate phylogeographical and population genetic analyses and species distribution modelling to develop a broad view of its current condition and prognosis for conservation. We use genotyping‐by‐sequencing to characterize the genomes of 142 S. longii individuals from 20 populations distributed throughout its range (New Jersey to Nova Scotia). We measure the distribution of genetic diversity in the species and reconstruct its phylogeographical history using the snapp and rase models. Extant populations of S. longii originated from a single refugium south of the Laurentide ice sheet around 25 ka. The genetic diversity of S. longii is exceedingly low, populations exhibit little genetic structure and the species is slightly inbred. Projected climate scenarios indicate that nearly half of extant populations of S. longii will be exposed to unsuitable climate by 2070. Similar changes in suitable habitat will occur for many other northern ACPF species—centres of diversity will shift northward and Nova Scotia may become the last refuges for those species not extinguished.     

Combining ecological niche modelling and morphology to assess the range‐wide population genetic structure of bobcats (Lynx rufus)

Despite a broad distribution, general habitat requirements, and a large dispersal potential, bobcats (Lynx rufus) exhibit a genetic division that longitudinally transects central North America. We investigated (1) whether the climate of the Last Glacial Maximum (LGM; 21 kya) isolated bobcats into refugia and also whether the current climate influences gene flow between the segregate populations and (2) whether the geographical patterns in cranial morphology reflect population identity. We created ecological niche models (ENMs) to evaluate climatic suitability and to estimate distributions of the disparate populations under both historical (LGM) and contemporary conditions. We used two‐dimensional geometric morphometric methods to evaluate variations in the cranium and mandible. These variations were then regressed across geographical variables to assess morphological differences throughout the range of the bobcat. ENMs projected onto LGM climate provided evidence of refugia during the LGM via increased suitability in the north‐west and south‐east portions of this species' range. Contemporarily, our models suggest that the Great Plains may be restricting bobcat migration and gene flow, effectively maintaining disparate populations. Morphological analyses identified a significant linear trend in shape variation across latitudinal and longitudinal gradients rather than distinct morphological divergence between lineages. Similar shape variations, however, did converge in approximate locations of assumed refugia. The findings of the present study provide a robust assessment of the biogeographical considerations for the population genetic structure of bobcats.     

Pleistocene climate changes shaped the divergence and demography of Asian populations of the great tit Parus major: evidence from phylogeographic analysis and ecological niche models

Different scales and frequencies of glaciations developed in Europe and Asia during the Pleistocene. Because species’ responses to climate change are influenced by interactive factors including ecology and local topography, the pattern and tempo of species diversification may vary significantly across regions. The great tit Parus major is a widespread Eurasian passerine with a range that encircles the central Asian desert and high‐altitude areas of the Tibetan Plateau. A number of genetic studies have assessed the effect of paleo‐climate changes on the distribution of the European population. However, none have comprehensively addressed how paleo‐climate change affected the distribution of the great tit in China, an apparent hotspot of P. major subspecific diversity. Here, we describe likely paleo‐climatic effects on P. major populations in China based on a combination of phylogeography and ecological niche models (ENMs). We sequenced three mitochondrial DNA markers from 28 populations (213 individuals), and downloaded 112 sequences from outside its Chinese range. As the first step in clarifying the intra‐specific relationships among haplotypes, we attempted to clarify the divergence and demography of populations in China. Phylogeographic analysis revealed that P. major is comprised of five highly divergent clades with geographic breaks corresponding to steep mountains and dry deserts. A previously undescribed monophyletic clade with high genetic diversity, stable niches and a long and independent evolutionary history was detected in the mountainous areas of southwest China. The estimated times at which these clades diverged was traced back to the Early‐Middle Pleistocene (2.19–0.61 mya). Contrary to the post‐LGM (the Last Glacial Maximum) expansion of European populations, demographic history indicates that Asian populations expanded before the LGM after which they remained relatively stable or grew slowly through the LGM. ENMs support this conclusion and predict a similar distribution in the present and the LGM. Our genetic and ecological results demonstrate that Pleistocene climate changes shaped the divergence and demography of P. major in China.     

Global warming threatens the persistence of Mediterranean brown trout

Current climate change exacerbates the environmental restrictions on temperate species inhabiting low latitude edges of their geographical ranges. We examined how temperature variations due to current and future climate change are likely to affect populations’ persistence of stream‐dwelling brown trout Salmo trutta at the vulnerable southern periphery of its range. Analysis of 33 years of air temperature data (1975–2007) by time‐series models indicated a significant upward trend and a pronounced shift in air temperature around 1986‐1987. This warming is associated with an ongoing population decline of brown trout, most likely caused by a loss of suitable thermal habitat in lower latitudes since the 1980s. Population decrease may not be attributed to physical habitat modification or angler pressure, as carrying capacity remained stable and populations were not overexploited. We developed regional temperature models, which predicted that unsuitable thermal habitat for brown trout increased by 93% when comparing climate conditions between 1975–1986 and 1993–2004. Predictions from climate envelope models showed that current climate change may be rendering unsuitable 12% of suitable thermal habitat each decade, resulting in an overall population decrease in the lower reaches of around 6% per year. Furthermore, brown trout catches markedly decreased 20% per year. Projections of thermal habitat loss under the ecologically friendly B2 SRES scenario showed that brown trout may lose half of their current suitable habitat within the study area by 2040 and become almost extinct by 2100. In parallel to the upstream movement of brown trout thermal habitat, warm water species are increasing their relative abundance in salmonid waters. Empirical evidence was provided of how current climate change threatens some of the most healthy native brown trout populations in Southern Europe and how forthcoming climate change is expected to further decrease the conservation status of the species.     

Montane refugia predict population genetic structure in the Large‐blotched Ensatina salamander

Understanding the biotic consequences of Pleistocene range shifts and fragmentation remains a fundamental goal in historical biogeography and evolutionary biology. Here, we combine species distribution models (SDM) from the present and two late Quaternary time periods with multilocus genetic data (mitochondrial DNA and microsatellites) to evaluate the effect of climate‐induced habitat shifts on population genetic structure in the Large‐blotched Ensatina (Ensatina eschscholtzii klauberi), a plethodontid salamander endemic to middle and high‐elevation conifer forest in the Transverse and Peninsular Ranges of southern California and northern Baja California. A composite SDM representing the range through time predicts two disjunct refugia, one in southern California encompassing the core of the species range and the other in the Sierra San Pedro Mártir of northern Baja California at the southern limit of the species range. Based on our spatial model, we would expect a pattern of high connectivity among populations within the northern refugium and, conversely, a pattern of isolation due to long‐term persistence of the Sierra San Pedro Mártir population. Our genetic results are consistent with these predictions based on the hypothetical refugia in that (i) historical measures of population connectivity among stable areas are correlated with gene flow estimates; and (ii) there is strong geographical structure between separate refugia. These results provide evidence for the role of recent climatic change in shaping patterns of population persistence and connectivity within the Transverse and Peninsular Ranges, an evolutionary hotspot.     

Climate structures genetic variation across a species' elevation range: a test of range limits hypotheses

Gene flow may influence the formation of species range limits, and yet little is known about the patterns of gene flow with respect to environmental gradients or proximity to range limits. With rapid environmental change, it is especially important to understand patterns of gene flow to inform conservation efforts. Here we investigate the species range of the selfing, annual plant, Mimulus laciniatus, in the California Sierra Nevada. We assessed genetic variation, gene flow, and population abundance across the entire elevation‐based climate range. Contrary to expectations, within‐population plant density increased towards both climate limits. Mean genetic diversity of edge populations was equivalent to central populations; however, all edge populations exhibited less genetic diversity than neighbouring interior populations. Genetic differentiation was fairly consistent and moderate among all populations, and no directional signals of contemporary gene flow were detected between central and peripheral elevations. Elevation‐driven gene flow (isolation by environment), but not isolation by distance, was found across the species range. These findings were the same towards high‐ and low‐elevation range limits and were inconsistent with two common centre‐edge hypotheses invoked for the formation of species range limits: (i) decreasing habitat quality and population size; (ii) swamping gene flow from large, central populations. This pattern demonstrates that climate, but not centre‐edge dynamics, is an important range‐wide factor structuring M. laciniatus populations. To our knowledge, this is the first empirical study to relate environmental patterns of gene flow to range limits hypotheses. Similar investigations across a wide variety of taxa and life histories are needed.     

Predicting impacts of climate change on habitat connectivity of Kalopanax septemlobus in South Korea

Understanding the drivers of habitat distribution patterns and assessing habitat connectivity are crucial for conservation in the face of climate change. In this study, we examined a sparsely distributed tree species, Kalopanax septemlobus (Araliaceae), which has been heavily disturbed by human use in temperate forests of South Korea. We used maximum entropy distribution modeling (MaxEnt) to identify the climatic and topographic factors driving the distribution of the species. Then, we constructed habitat models under current and projected climate conditions for the year 2050 and evaluated changes in the extent and connectivity of the K. septemlobus habitat. Annual mean temperature and terrain slope were the two most important predictors of species distribution. Our models predicted the range shift of K. septemlobus toward higher elevations under medium-low and high emissions scenarios for 2050, with dramatic reductions in suitable habitat (51% and 85%, respectively). In addition, connectivity analysis indicated that climate change is expected to reduce future levels of habitat connectivity. Even under the Representative Construction Pathway (RCP) 4.5 medium-low warming scenario, the projected climate conditions will decrease habitat connectivity by 78%. Overall, suitable habitats for K. septemlobus populations will likely become more isolated depending on the severity of global warming. The approach presented here can be used to efficiently assess species and habitat vulnerability to climate change.     

Cryptic habitats and cryptic diversity: unexpected patterns of connectivity and phylogeographical breaks in a Mediterranean endemic marine cave mysid

The marine cave‐dwelling mysid Hemimysis margalefi is distributed over the whole Mediterranean Sea, which contrasts with the poor dispersal capabilities of this brooding species. In addition, underwater marine caves are a highly fragmented habitat which further promotes strong genetic structuring, therefore providing highly informative data on the levels of marine population connectivity across biogeographical regions. This study investigates how habitat and geography have shaped the connectivity network of this poor disperser over the entire Mediterranean Sea through the use of several mitochondrial and nuclear markers. Five deeply divergent lineages were observed among H. margalefi populations resulting from deep phylogeographical breaks, some dating back to the Oligo‐Miocene. Whether looking at the intralineage or interlineage levels, H. margalefi populations present a high genetic diversity and population structuring. This study suggests that the five distinct lineages observed in H. margalefi actually correspond to as many separate cryptic taxa. The nominal species, H. margalefi sensu stricto, corresponds to the westernmost lineage here surveyed from the Alboran Sea to southeastern Italy. Typical genetic breaks such as the Almeria‐Oran Front or the Siculo‐Tunisian Strait do not appear to be influential on the studied loci in H. margalefi sensu stricto. Instead, population structuring appears more complex and subtle than usually found for model species with a pelagic dispersal phase. The remaining four cryptic taxa are all found in the eastern basin, but incomplete lineage sorting is suspected and speciation might still be in process. Present‐day population structure of the different H. margalefi cryptic species appears to result from past vicariance events started in the Oligo‐Miocene and maintained by present‐day coastal topography, water circulation and habitat fragmentation.