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.     

Rethinking long‐term vegetation dynamics: multiple glacial refugia and local expansion of a species complex

Evidence is accumulating that some arcto‐boreal plant taxa persisted through the last glacial maximum (LGM) in Alaska and adjacent Canada. However, the spatial patterns of glacial persistence and associated postglacial colonization remain largely unknown. In this study, we investigated the LGM refugia of an alder (Alnus) species complex (n = 3 taxa) and assess the spatiotemporal dynamics of Alnus in this vast region. Specifically, we conducted high‐throughput DNA sequencing (ddRADseq) on Alnus foliar samples collected from a dense population network to investigate patterns of genetic structure and infer the presence of glacial lineages. Species distribution modeling (SDM) was used to investigate the probability and possible locations of glacial persistence. These analyses were integrated and then compared with fossil pollen data to identify the locations of refugial populations and spatial patterns of postglacial colonization. Our genetic analyses revealed two glacial lineages with separate geographic origins for each Alnus taxon, suggesting that the genus persisted in multiple LGM refugia. Non‐overlapping hindcast distributions based on SDMs further support the presence of multiple, spatially distinct refugia. These ddRADseq and SDM results, in conjunction with reassessment of fossil pollen records, suggest that Alnus expanded from several population nuclei that existed during the LGM and coalesced during the Holocene to form its present range. These results challenge the unidirectional model for postglacial vegetation expansion, implying that climate buffering associated with landscape heterogeneity and adaptation to millennial‐scale environmental variability played important roles in driving late‐Quaternary population dynamics.     

Spatial scale affects novel and disappeared climate change projections in Alaska

The formation of novel and disappeared climates between the last glacial maximum (LGM) and the present is important to consider to understand the expansion and contraction of species niches and distributions, as well as the formation and loss of communities and ecological interactions over time. Our choice in climate data resolution has the potential to complicate predictions of the ecological impacts of climate change, since climate varies from local to global scales and this spatial variation is reflected in climate data. To address this issue, we downscaled LGM and modern (1975–2005) 30‐year averaged climate data to 60‐m resolution for the entire state of Alaska for 10 different climate variables, and then upsampled each variable to coarser resolutions (60 m to 12 km). We modeled the distributions of novel and disappeared climates to evaluate the locations and fractional area of novel and disappeared climates for each of our climate variables and resolutions. Generally, novel and disappeared climates were located in southern Alaska, although there were cases where some disappeared climates existed within coastal and interior Alaska. Climate resolution affected the fractional area of novel and disappeared climates in three patterns: As the spatial resolution of climate became coarser, the fractional area of novel and disappeared climates (a) increased, (b) decreased, or (c) had no explainable relationship. Overall, we found the use of coarser climate data increased the fractional area of novel and disappeared climates due to decreased environmental variability and removal of climate extremes. Our results reinforce the importance of downscaling coarse climate data and suggest that studies analyzing the effects of climate change on ecosystems may overestimate or underestimate their conclusions when utilizing coarse climate data.     

Southern montane populations did not contribute to the recolonization of West Siberian Plain by Siberian larch (Larix sibirica): a range‐wide analysis of cytoplasmic markers

While many species were confined to southern latitudes during the last glaciations, there has lately been mounting evidence that some of the most cold‐tolerant species were actually able to survive close to the ice sheets. The contribution of these higher latitude outposts to the main recolonization thrust remains, however, untested. In the present study, we use the first range‐wide survey of genetic diversity at cytoplasmic markers in Siberian larch (Larix sibirica; four mitochondrial (mt) DNA loci and five chloroplast (cp) DNA SSR loci) to (i) assess the relative contributions of southern and central areas to the current L. sibirica distribution range; and (ii) date the last major population expansion in both L. sibirica and adjacent Larix species. The geographic distribution of cpDNA variation was uninformative, but that of mitotypes clearly indicates that the southernmost populations, located in Mongolia and the Tien‐Shan and Sayan Mountain ranges, had a very limited contribution to the current populations of the central and northern parts of the range. It also suggests that the contribution of the high latitude cryptic refugia was geographically limited and that most of the current West Siberian Plain larch populations likely originated in the foothills of the Sayan Mountains. Interestingly, the main population expansion detected through Approximate Bayesian Computation (ABC) in all four larch species investigated here pre‐dates the LGM, with a mode in a range of 220 000–1 340 000 years BP. Hence, L. sibirica, like other major conifer species of the boreal forest, was strongly affected by climatic events pre‐dating the Last Glacial Maximum.     

Phylogeography and niche modelling of the relict plant Amborella trichopoda (Amborellaceae) reveal multiple Pleistocene refugia in New Caledonia

Amborella trichopoda Baill. (Amborellaceae, Amborellales), the sole living member of the sister group to all other extant angiosperms, is endemic to New Caledonia. We addressed the intraspecific phylogeography of Amborella by investigating whether its present population genetic structure could be related to its current and past habitats. We found moderate range‐wide genetic diversity based on nuclear microsatellite data and detected four well‐differentiated, geographically distinct genetic groups using Bayesian clustering analyses. We modelled the ecological niche of Amborella based on the current climatic and environmental conditions. The predictive ability of the model was very good throughout the Central East mainland zone, but Amborella was predicted in the northern part of the island where this plant has not been reported. Furthermore, no significant barrier was detected based on habitat suitability that could explain the genetic differentiation across the area. Conversely, we found that the main genetic clusters could be related to the distribution of the suitable habitat at the last glacial maximum (LGM, c. 21 000 years BP), when Amborella experienced a dramatic 96.5% reduction in suitable area. At least two lineages survived in distinct putative refugia located in the Massif des Lèvres and in the vicinity of Mount Aoupinié. Our findings finally confirmed the importance of LGM rainforest refugia in shaping the current intra‐ and interspecific diversity in New Caledonian plants and revealed the possibility of an as yet unreported refugium. The combination of niche modelling and population genetics thereby offered novel insight into the biogeographical history of an emblematic taxon.     

Post‐glacial distribution of a Mongolian mayfly inferred from population genetic analysis

相似文献   

The population demography of Betula maximowicziana,a cool‐temperate tree species in Japan,in relation to the last glacial period: its admixture‐like genetic structure is the result of simple population splitting not admixing

Conservation of the local genetic variation and evolutionary integrity of economically and ecologically important trees is a key aspect of studies involving forest genetics, and a population demographic history of the target species provides valuable information for this purpose. Here, the genetic structure of 48 populations of Betula maximowicziana was assessed using 12 expressed sequence tag–simple sequence repeat (EST‐SSR) markers. Genetic diversity was lower in northern populations than southern ones and structure analysis revealed three groups: northern and southern clusters and an admixed group. Eleven more genomic‐SSR loci were added and the demographic history of these three groups was inferred by approximate Bayesian computation (ABC). The ABC revealed that a simple split scenario was much more likely than isolation with admixture, suggesting that the admixture‐like structure detected in this species was due to ancestral polymorphisms. The ABC analysis suggested that the population growth and divergence of the three groups occurred 96 800 (95% CI, 20 500–599 000) and 28 300 (95% CI, 8700–98 400) years ago, respectively. We need to be aware of several sources of uncertainty in the inference such as assumptions about the generation time, overlapping of generations, confidence intervals of the estimated parameters and the assumed model in the ABC. However, the results of the ABC together with the model‐based maps of reconstructed past species distribution and palaeoecological data suggested that the modern genetic structure of B. maximowicziana originated prior to the last glacial maximum (LGM) and that some populations survived in the northern range even during the LGM.     

Climate hindcasts: exploring the disjunct distribution of Diopatra biscayensis

The ecosystem engineer onuphid polychaete Diopatra biscayensis has a continuous population in the Bay of Biscay from the Cantabria coast in Spain to southern Brittany in France. A group of disjunct populations also are found in the English Channel, separated from the Biscay population by more than 400 coastal kilometers. It remains unclear whether D. biscayensis is native to the Bay of Biscay; it is also debated whether the disjunct populations in the English Channel are relics of a formerly continuous population, or the product of recent introductions through aquaculture. Here, we use climate hindcasts to explore hypotheses about the D. biscayensis historical distribution in Europe. If D. biscayensis is native, its range would have been restricted to southern Iberia and the Mediterranean during the Last Glacial Maximum (21,000 BP). However, the species is completely absent from both regions today, further supporting its interpretation as a non‐native species. If it was historically present in Europe, the climate hindcasts are congruent with range contraction in the Last Glacial Maximum (21,000 BP), expansion in the Mid‐Holocene Warm Period (6000 BP), and contraction again in the past 1000 years (850–1850), prior to the first reports of D. biscayensis on the Spanish and French Atlantic coasts. However, the simulations do not support there being climatic refugia along the English Channel coast that would account for the existence of relic populations. Taken together, the evidence suggests that D. biscayensis has been introduced to the Bay of Biscay, and that disjunct populations in the English Channel are the result of recent transport through human activities, perhaps aquaculture.     

Niche shifting in response to warming climate after the last glacial maximum: inference from genetic data and niche assessments in the chisel‐toothed kangaroo rat (Dipodomys microps)

During Pleistocene glacial‐interglacial cycles, the geographic range is often assumed to have shifted as a species tracks its climatic niche. Alternatively, the geographic range would not necessarily shift if a species can adapt in situ to a changing environment. The potential for a species to persist in place might increase with the diversity of habitat types that a species exploits. We evaluate evidence for either range shift or range stability between the last glacial maximum (LGM) and present time in the chisel‐toothed kangaroo rat (Dipodomys microps), an endemic of the Great Basin and Mojave deserts. We modeled how the species’ range would have changed if the climatic niche of the species remained conserved between the LGM and present time. The climatic models imply that if D. microps inhabited the same climatic niche during the LGM as it does today, the species would have persisted primarily within the warm Mojave Desert and expanded northwards into the cold Great Basin only after the LGM. Contrary to the climatic models, the mitochondrial DNA assessment revealed signals of population persistence within the current distribution of the species throughout at least the latest glacial‐interglacial cycle. We concluded that D. microps did not track its climatic niche during late Pleistocene oscillations, but rather met the challenge of a changing environment by shifting its niche and retaining large portions of its distribution. We speculate that this kind of response to fluctuating climate was possible because of ‘niche drifting’, an alteration of the species’ realized niche due to plasticity in various biological characters. Our study provides an example of an approach to reconstruct species’ responses to past climatic changes that can be used to evaluate whether and to what extent taxa have capacity to shift their niches in response to the changing environment – information becoming increasingly important to predicting biotic responses to future environmental changes.     

Quaternary climate instability is correlated with patterns of population genetic variability in Bombus huntii

Climate oscillations have left a significant impact on the patterns of genetic diversity observed in numerous taxa. In this study, we examine the effect of Quaternary climate instability on population genetic variability of a bumble bee pollinator species, Bombus huntii in western North America. Pleistocene and contemporary B. huntii habitat suitability (HS) was estimated with an environmental niche model (ENM) by associating 1,035 locality records with 10 bioclimatic variables. To estimate genetic variability, we genotyped 380 individuals from 33 localities at 13 microsatellite loci. Bayesian inference was used to examine population structure with and without a priori specification of geographic locality. We compared isolation by distance (IBD) and isolation by resistance (IBR) models to examine population differentiation within and among the Bayesian inferred genetic clusters. Furthermore, we tested for the effect of environmental niche stability (ENS) on population genetic diversity with linear regression. As predicted, high‐latitude B. huntii habitats exhibit low ENS when compared to low‐latitude habitats. Two major genetic clusters of B. huntii inhabit western North America: (a) a north genetic cluster predominantly distributed north of 28°N and (b) a south genetic cluster distributed south of 28°N. In the south genetic cluser, both IBD and IBR models are significant. However, in the north genetic cluster, IBD is significant but not IBR. Furthermore, the IBR models suggest that low‐latitude montane populations are surrounded by habitat with low HS, possibly limiting dispersal, and ultimately gene flow between populations. Finally, we detected high genetic diversity across populations in regions that have been climatically unstable since the last glacial maximum (LGM), and low genetic diversity across populations in regions that have been climatically stable since the LGM. Understanding how species have responded to climate change has the potential to inform management and conservation decisions of both ecological and economic concerns.     

Emerging climate‐driven disturbance processes: widespread mortality associated with snow‐to‐rain transitions across 10° of latitude and half the range of a climate‐threatened conifer

Climate change is causing rapid changes to forest disturbance regimes worldwide. While the consequences of climate change for existing disturbance processes, like fires, are relatively well studied, emerging drivers of disturbance such as snow loss and subsequent mortality are much less documented. As the climate warms, a transition from winter snow to rain in high latitudes will cause significant changes in environmental conditions such as soil temperatures, historically buffered by snow cover. The Pacific coast of North America is an excellent test case, as mean winter temperatures are currently at the snow–rain threshold and have been warming for approximately 100 years post‐Little Ice Age. Increased mortality in a widespread tree species in the region has been linked to warmer winters and snow loss. Here, we present the first high‐resolution range map of this climate‐sensitive species, Callitropsis nootkatensis (yellow‐cedar), and document the magnitude and location of observed mortality across Canada and the United States. Snow cover loss related mortality spans approximately 10° latitude (half the native range of the species) and 7% of the overall species range and appears linked to this snow–rain transition across its range. Mortality is commonly >70% of basal area in affected areas, and more common where mean winter temperatures is at or above the snow–rain threshold (>0 °C mean winter temperature). Approximately 50% of areas with a currently suitable climate for the species (相似文献   

Waterscape genetics of the yellow perch ( Perca flavescens): patterns across large connected ecosystems and isolated relict populations

Comparisons of a species' genetic diversity and divergence patterns across large connected populations vs. isolated relict areas provide important data for understanding potential response to global warming, habitat alterations and other perturbations. Aquatic taxa offer ideal case studies for interpreting these patterns, because their dispersal and gene flow often are constrained through narrow connectivity channels that have changed over geological time and/or from contemporary anthropogenic perturbations. Our research objective is to better understand the interplay between historic influences and modern‐day factors (fishery exploitation, stocking supplementation and habitat loss) in shaping population genetic patterns of the yellow perch Perca flavescens (Percidae: Teleostei) across its native North American range. We employ a modified landscape genetics approach, analysing sequences from the entire mitochondrial DNA control region and 15 nuclear DNA microsatellite loci of 664 spawning adults from 24 populations. Results support that perch from primary glacial refugium areas (Missourian, Mississippian and Atlantic) founded contemporary northern populations. Genetic diversity today is highest in southern (never glaciated) populations and also is appreciable in northern areas that were founded from multiple refugia. Divergence is greater among isolated populations, both north and south; the southern Gulf Coast relict populations are the most divergent, reflecting their long history of isolation. Understanding the influence of past and current waterway connections on the genetic structure of yellow perch populations may help us to assess the roles of ongoing climate change and habitat disruptions towards conserving aquatic biodiversity.     

Postglacial range shift and demographic expansion of the marine intertidal snail Batillaria attramentaria

To address the impacts of past climate changes, particularly since the last glacial period, on the history of the distribution and demography of marine species, we investigated the evolutionary and demographic responses of the intertidal batillariid gastropod, Batillaria attramentaria, to these changes, using the snail as a model species in the northwest Pacific. We applied phylogeographic and divergence population genetic approaches to mitochondrial COI sequences from B. attramentaria. To cover much of its distributional range, 197 individuals collected throughout Korea and 507 publically available sequences (mostly from Japan) were used. Finally, a Bayesian skyline plot (BSP) method was applied to reconstruct the demographic history of this species. We found four differentiated geographic groups around Korea, confirming the presence of two distinct, geographically subdivided haplogroups on the Japanese coastlines along the bifurcated routes of the warm Tsushima and Kuroshio Currents. These two haplogroups were estimated to have begun to split approximately 400,000 years ago. Population divergence analysis supported the hypothesis that the Yellow Sea was populated by a northward range expansion of a small fraction of founders that split from a southern ancestral population since the last glacial maximum (LGM: 26,000–19,000 years ago), when the southern area became re‐submerged. BSP analyses on six geographically and genetically defined groups in Korea and Japan consistently demonstrated that each group has exponentially increased approximately since the LGM. This study resolved the phylogeography of B. attramentaria as a series of events connected over space and time; while paleoceanographic conditions determining the connectivity of neighboring seas in East Asia are responsible for the vicariance of this species, the postglacial sea‐level rise and warming temperatures have played a crucial role in rapid range shifts and broad demographic expansions of its populations.     

Harbor seal use of glacier ice and terrestrial haul‐outs in the Kenai Fjords,Alaska

Harbor seals, Phoca vitulina, use diverse haul‐out substrates including ice calved by tidewater glaciers. Numbers of seals at glacial and terrestrial haul‐outs on the southeastern Kenai Peninsula, Alaska, were assessed using aerial, vessel, and video surveys. Mean annual abundance at glacial and terrestrial haul‐outs differed temporally. From 2004 to 2011, numbers of seals counted during the molt increased 5.4%/yr at glacial haul‐outs and 9%/yr at terrestrial haul‐outs while numbers of pups increased 5.0%/yr at glacial sites and 1.5%/yr at terrestrial sites. Numbers of seals without pups counted during pupping increased 7.96%/yr at glacial sites and 5.1%/yr at terrestrial sites. Results indicate that pupping and molting locations are not equivalent and population monitoring during the molt does not necessarily reflect habitat association of pupping seals. Ratios of pups to total seals counted during pupping and the subsequent molt were used to contrast habitat use. Low proportions of pups at terrestrial haul‐outs, relative to most glacial haul‐outs, indicate an overall preference for pupping in glacial haul‐outs. High proportions of pups at most glacial sites (during pupping and molting) suggest reduced use of tidewater glacier habitats by nonbreeders and molting seals. Results suggest more seals associate with glacial haul‐outs than currently estimated.