Speciation in the White-breasted Nuthatch (Sitta carolinensis): a multilocus perspective

Inferring the evolutionary and ecological processes that have shaped contemporary species distributions using the geographic distribution of gene lineages is the principal goal of phylogeographic research. Researchers in the field have recognized that inferences made from a single gene, often mitochondrial, can be informative regarding the pattern of diversification but lack conclusive information regarding the evolutionary mechanisms that led to the observed patterns. Here, we use a multilocus (20 loci) data set to explore the evolutionary history of the White‐breasted Nuthatch (Sitta carolinensis). A previous single‐locus study found S. carolinensis is comprised of four reciprocally monophyletic clades geographically restricted to the pine and oak forests of: (i) eastern North America, (ii) southern Rocky Mountain and Mexican Mountain ranges, (iii) Eastern Sierra Nevada and Northern Rocky Mountains and (iv) Pacific slope of North America. The diversification of the clades was attributed to the fragmentation of North American pine and oak woodlands in the Pliocene with subsequent divergences owing to the Pleistocene glacial cycles. Principal component, clustering and species tree analyses of the multilocus data resolved the same four groups or lineages found in the single‐locus study. Coalescent analyses and hypothesis testing of nested isolation and migration models indicate that isolation and not gene flow has been the major evolutionary mechanism responsible for shaping genetic variation, and all the divergence events within S. carolinensis have occurred in response to the Pleistocene glacial cycles.     

Mitochondrial DNA evolution in the Anaxyrus boreas species group

The Anaxyrus boreas species group currently comprises four species in western North America including the broadly distributed A. boreas, and three localized species, Anaxyrus nelsoni, Anaxyrusexsul and Anaxyrus canorus. Phylogenetic analyses of the mtDNA 12S rDNA, cytochrome oxidase I, control region, and restriction sites data, identified three major haplotype clades. The Northwest clade (NW) includes both subspecies of A. boreas and divergent minor clades in the middle Rocky Mountains, coastal, and central regions of the west and Pacific Northwest. The Southwest (SW) clade includes A. exsul, A. nelsoni, and minor clades in southern California. Anaxyrus canorus, previously identified as paraphyletic, has populations in both the NW and SW major clades. The Eastern major clade (E) includes three divergent lineages from southern Utah, the southern Rocky Mountains, and north of the Great Basin at the border of Utah and Nevada. These results identify new genetic variation in the eastern portion of the toad's range and are consistent with previous regional studies from the west coast. Low levels of control region sequence divergence between major clades (2.2-4.7% uncorrected pair-wise distances) are consistent with Pleistocene divergence and suggest that the phylogeographic history of the group was heavily influenced by dynamic Pleistocene glacial and climatic changes, and especially pluvial changes, in western North America. Results reported here may impact conservation plans in that the current taxonomy does not reflect the diversity in the group.     

Molecular evidence for Pleistocene glacial cycles driving diversification of a North American desert spider, Agelenopsis aperta

The influence of historical climatic vs. geological changes on species diversification patterns was investigated in a widely distributed North American desert spider, Agelenopsis aperta (Araneae: Agelenidae), with particular reference to Pleistocene glacial cycles and earlier patterns of mountain building. Levels of sequence divergence obtained from the mitochondrial gene, cytochrome oxidase I, dated to the Pleistocene, eliminating Rocky Mountain orogeny as a cause of diversification, as orogeny ended 4 million years ago. The results of phylogenetic and network analyses showed the presence of three geographically defined clades, which were consistent with the presence of at least three glacial refugia: (i) east of the Rocky Mountains; (ii) between the Rocky Mountains and Sierra Nevadas; and (iii) west of the Sierra Nevadas. In addition, populations within the Rocky Mountains exhibited significantly lower genetic diversity than populations east of the Rocky Mountains and the haplotypes found within the Rockies were a subset of eastern haplotypes. These patterns suggest that a post-Pleistocene range expansion occurred out of an eastern glacial refugium into the Rocky Mountains. Examination of phylogeographical studies of other North American desert taxa indicated that mountain building explained diversification patterns more effectively for some taxa but Pleistocene climate change was more important for others, including A. aperta.     

Geological barriers and restricted gene flow in the holarctic skipper Hesperia comma (Hesperiidae)

Patterns of genetic variation within a species may be a consequence of historical factors, such as past fragmentation, as well as current barriers to gene flow. Using sequence data from the mitochondrial cytochrome oxidase subunit II region (COII) and the nuclear gene wingless, we conducted a phylogeographical study of the holarctic skipper Hesperia comma to elucidate patterns of genetic diversity and to infer historical and contemporary processes maintaining genetic variation. One hundred and fifty-one individuals were sampled from throughout North America and Eurasia, focusing on California and adjacent regions in the western United States where morphological diversity is highest compared to the rest of the range. Analyses of sequence data obtained from both genes revealed a well-supported division between the Old and New World. Within western North America, wingless shows little geographical structure, while a hierarchical analysis of genetic diversity of COII sequences indicates three major clades: a western clade in Oregon and Northern California, an eastern clade including the Great Basin, Rocky Mountains and British Columbia, and a third clade in southern California. The Sierra Nevada and the Transverse Ranges appear to be the major barriers to gene flow for H. comma in the western United States. Relatively reduced haplotype diversity in Eurasia compared to North America suggests that populations on the two continents have been affected by different historical processes.     

Extensive chloroplast haplotype variation indicates Pleistocene hybridization and radiation of North American Arabis drummondii, A. x divaricarpa, and A. holboellii (Brassicaceae)

Arabis drummondii, A. holboellii and their hybrid A. x divaricarpa are widespread perennials of open habitats in North America. A phylogenetic analysis based on noncoding chloroplast DNA sequences (trnL intron and trnL/F intergenic spacer) resolved A. drummondii as a monophyletic taxon, but found A. holboellii to bear chloroplast haplotypes from highly diverged evolutionary lineages. This raised the question of a possible polyphyletic origin of A. holboellii. Arabis x divaricarpa was found to be of recent and polytopic origin, a result consistent with its presumed hybrid origin. One hundred and three chloroplast haplotypes were detected within 719 Arabis accessions investigated. The majority of chloroplast-types were estimated to have arisen prior to the Wisconsin glaciation. Phylogeographical analysis using nested clade analysis, suggested for A. holboellii (i). past fragmentation events, partitioning genetic variation in several instances between the Sierra Nevada, the Southern Rocky Mountains and the Colorado Plateau on the one hand and the Central to Northern Rockies of the United States and adjacent Cascades on the other; and for both parental species (ii). recolonization of major areas formerly covered by the Wisconsin glaciation by three haplotypes; and (iii). restricted gene flow indicating isolation by distance in areas south of the last glacial maximum. Arabis x divaricarpa is closely codistributed with its parental species and resampled their haplotypes. The highest genetic diversity was found in the Rocky Mountains from Idaho and Montana south to Utah and Colorado. This area was further hypothesized to have played a major role in the origin of both parental species and probably represented an important glacial refugium. However, evidence for glacial refugia was also found in arctic and boreal regions of Alaska and near the Great Lakes. In comparison to nuclear ribosomal internal transcribed spacer data, chloroplast DNA divergence was very high and evidently predated the origin of A. drummondii and possibly A. holboellii. Divergence of major chloroplast lineages dates back to the middle of the Pleistocene at least. Extensive hybridization is the most likely evolutionary factor working on A. holboellii to explain the revealed discrepancy in nuclear DNA and chloroplast DNA diversification.     

Differential regional response of the bushy‐tailed woodrat (Neotoma cinerea) to late Quaternary climate change

Aim To reconstruct the regional biogeographical history of the bushy‐tailed woodrat, Neotoma cinerea (Rodentia: Cricetidae), across its distribution using multiple sources of information, including genetic data, ecological niche models and the palaeorecord. Location Western North America. Methods We analysed complete cytochrome b gene (1143 bp) sequences from 182 specimens of N. cinerea using Bayesian and coalescent methods to infer phylogenetic relationships, time of major divergences, and recent demographic trends. For comparison, we developed clade‐specific ecological niche models for groups of interest and analysed spatial trends of N. cinerea in the palaeorecord as well as temporal frequency trends across strata of individual palaeomiddens. Results We found two largely allopatric clades within N. cinerea, with several regionally distinct subclades showing contrasting recent population dynamics. Niche models showed consistent habitat at the Last Glacial Maximum (LGM) and modern times in the Rocky Mountains and northern United States, while the Great Basin may have been markedly less suitable at the LGM than today. The palaeorecord showed great spatio‐temporal variability in the presence of N. cinerea, but documents broad‐scale patterns of occupancy and regional population trends. Main conclusions The Quaternary dynamics and evolutionary history of N. cinerea appear to have been shaped by both vicariant events associated with geographical barriers and the availability of suitable habitat through time. Divergence of the two major clades dates to the Pliocene–Pleistocene transition, with clades separated by the Green and Colorado rivers and northern Rocky Mountain Pleistocene glaciations. We found largely concordant genetic, niche model and palaeorecord patterns suggesting long‐term population stability in the Rocky Mountains, while extant clades in the Great Basin and far north appear to have expanded or re‐expanded into these regions relatively recently. Furthermore, disjunct haplotype distributions, regional demographic history, and historical distribution of suitable habitat suggest that the Great Basin has been a particularly dynamic region.     

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.     

The roles of Neogene geology and late Pleistocene lake levels in shaping the genetic structure of the Lahontan redside shiner Richardsonius egregius (Teleostei: Cyprinidae)

Shifting drainage patterns in western North America, shaped by geological activity and changing global climates, have influenced the evolution of many aquatic taxa. We investigated the role of late Pleistocene high stands in pluvial Lake Lahontan on the genetic structure of Richardsonius egregius, a minnow endemic to the Lahontan Basin of the western Great Basin. We used the mitochondrial cytochrome b gene to generate a phylogeny and assess intraspecific genetic diversity, to estimate divergence times between clades, and to evaluate whether gene flow currently occurs. The results obtained show that R. egregius exhibits genetic divergence between eastern and western Lahontan Basin populations. Divergence time estimates show that intraspecific genetic diversification began in the Pliocene or early Pleistocene, before the pluvial lake high stands associated with the last glacial maximum. These results imply that the fluctuating water levels in pluvial Lake Lahontan had a minimal effect on shaping the genetic architecture of R. egregius. Coalescent analyses using the immigration with migration model show that contemporary gene flow between eastern and western Lahontan Basin populations does not occur. © 2011 The Linnean Society of London, Biological Journal of the Linnean Society, 2011, 104 , 163–176.     

Phylogeography of Douglas‐fir based on mitochondrial and chloroplast DNA sequences: testing hypotheses from the fossil record

The integration of fossil and molecular data can provide a synthetic understanding of the ecological and evolutionary history of an organism. We analysed range‐wide maternally inherited mitochondrial DNA and paternally inherited chloroplast DNA sequence data with coalescent simulations and traditional population genetic methods to test hypotheses of population divergence generated from the fossil record of Douglas‐fir (Pseudotsuga menziesii), an ecologically and economically important western North American conifer. Specifically, we tested (i) the hypothesis that the Pliocene orogeny of the Cascades and Sierra Nevada caused the divergence of coastal and Rocky Mountain Douglas‐fir varieties; and (ii) the hypothesis that multiple glacial refugia existed on the coast and in the Rocky Mountains. We found that Douglas‐fir varieties diverged about 2.11 Ma (4.37 Ma–755 ka), which could be consistent with a Pliocene divergence. Rocky Mountain Douglas‐fir probably resided in three or more glacial refugia. More variable molecular markers would be required to detect the two coastal refugia suggested in the fossil record. Comparison of mitochondrial DNA and chloroplast DNA variation revealed that gene flow via pollen linked populations isolated from seed exchange. Postglacial colonization of Canada from coastal and Rocky Mountain refugia near the ice margin at the Last Glacial Maximum produced a wide hybrid zone among varieties that formed almost exclusively by pollen exchange and chloroplast DNA introgression, not seed exchange. Postglacial migration rates were 50–165 m/year, insufficient to track projected 21st century warming in some regions. Although fossil and genetic data largely agree, each provides unique insights.     

Patterns of morphological, biochemical, and molecular evolution in the Oeneis chryxus complex (Lepidoptera: Satyridae): a test of historical biogeographical hypotheses.

Surveys of allozyme allele frequency and mitochondrial DNA (mtDNA) sequence variation were employed to test historical biogeographical hypotheses on the origin and unique distribution of the synchronized biennial, high-altitude butterflies of the Oeneis chryxus complex in western North America. Populations of O. c. stanislaus and O. ivallda from the central and northern Sierra Nevada are indistinguishable by use of allozyme allele frequency data, possessed nearly identical mtDNA cytochrome oxidase subunit 1 (COI) haplotypes, and were found to be relatively distantly related to O. c. chryxus from the Snake Range in eastern Nevada. However, individuals of O. ivallda from Piute Pass in the southern Sierra Nevada are more variable, with some individuals sharing mtDNA characteristics with O. c. chryxus. We find little support for the hypothesis proposed by W. Hovanitz in 1940 that O. c. stanislaus invaded the central Sierra Nevada from across the Great Basin and displaced O. ivallda, but cannot reject the hypothesis that ancestral Oeneis dispersed across the Great Basin to California. This result is congruent with hypotheses of dispersal across the Great Basin for the origin of some Sierran alpine organisms.     

Population genomic evidence for multiple Pliocene refugia in a montane‐restricted harvestman (Arachnida,Opiliones, Sclerobunus robustus) from the southwestern United States

The integration of ecological niche modelling into phylogeographic analyses has allowed for the identification and testing of potential refugia under a hypothesis‐based framework, where the expected patterns of higher genetic diversity in refugial populations and evidence of range expansion of nonrefugial populations are corroborated with empirical data. In this study, we focus on a montane‐restricted cryophilic harvestman, Sclerobunus robustus, distributed throughout the heterogeneous Southern Rocky Mountains and Intermontane Plateau of southwestern North America. We identified hypothetical refugia using ecological niche models (ENMs) across three time periods, corroborated these refugia with population genetic methods using double‐digest RAD‐seq data and conducted population‐level phylogenetic and divergence dating analyses. ENMs identify two large temporally persistent regions in the mid‐latitude highlands. Genetic patterns support these two hypothesized refugia with higher genetic diversity within refugial populations and evidence for range expansion in populations found outside hypothesized refugia. Phylogenetic analyses identify five to six genetically divergent, geographically cohesive clades of S. robustus. Divergence dating analyses suggest that these separate refugia date to the Pliocene and that divergence between clades pre‐dates the late Pleistocene glacial cycles, while diversification within clades was likely driven by these cycles. Population genetic analyses reveal effects of both isolation by distance (IBD) and isolation by environment (IBE), with IBD more important in the continuous mountainous portion of the distribution, while IBE was stronger in the populations inhabiting the isolated sky islands of the south. Using model‐based coalescent approaches, we find support for postdivergence migration between clades from separate refugia.     

Postglacial colonisation patterns and the role of isolation and expansion in driving diversification in a passerine bird

Pleistocene glacial cycles play a major role in diversification and speciation, although the relative importance of isolation and expansion in driving diversification remains debated. We analysed mitochondrial DNA sequence data from 15 great reed warbler (Acrocephalus arundinaceus) populations distributed over the vast Eurasian breeding range of the species, and revealed unexpected postglacial expansion patterns from two glacial refugia. There were 58 different haplotypes forming two major clades, A and B. Clade A dominated in Western Europe with declining frequencies towards Eastern Europe and the Middle East, but showed a surprising increase in frequency in Western and Central Asia. Clade B dominated in the Middle East, with declining frequencies towards north in Central and Eastern Europe and was absent from Western Europe and Central Asia. A parsimonious explanation for these patterns is independent postglacial expansions from two isolated refugia, and mismatch distribution analyses confirmed this suggestion. Gene flow analyses showed that clade A colonised both Europe and Asia from a refugium in Europe, and that clade B expanded much later and colonised parts of Europe from a refugium in the Middle East. Great reed warblers in the eastern parts of the range have slightly paler plumage than western birds (sometimes treated as separate subspecies; A. a. zarudnyi and A. a. arundinaceus, respectively) and our results suggest that the plumage diversification took place during the easterly expansion of clade A. This supports the postglacial expansion hypothesis proposing that postglacial expansions drive diversification in comparatively short time periods. However, there is no indication of any (strong) reproductive isolation between clades and our data show that the refugia populations became separated during the last glaciation. This is in line with the Pleistocene speciation hypothesis invoking that much longer periods of time in isolation are needed for speciation to occur.     

Phylogeography and species biogeography of montane Great Basin stoneflies

Sky islands are ideal systems for determining the effects of climatic oscillations on species distributions and genetic structure. Our study focused on montane stonefly populations in the Great Basin of western North America. We used niche-based distribution modelling, phylogeography and traditional species-based biogeography to test several hypotheses as follows: (i) genetic differentiation among Doroneuria baumanni populations will be independent of hydrologic connectivity (headwater model); (ii) Sky islands were colonized when habitat was more continuous and populations likely experienced multiple expansions and contractions; (iii) Colonization events were coincident with the late Pleistocene and Holocene; and (iv) Shared topography and climate history will result in concordant patterns of genetic differentiation in D. baumanni and occurrences of 32 stonefly species across the region. Overall, Φ(ST) 's and coalescent-based estimates of migration were consistent with the headwater model. Maximum likelihood and Bayesian gene trees identified three major nonoverlapping east-west clades. Distribution modelling indicated more suitable habitat in the Great Basin during the Last Glacial Maximum than at present, but none during the last interglacial period. Demographic analyses showed evidence of population expansion in one of the three major east-west clades. Intra-clade divergence times (60,000-183,000ybp) were well within the late Pleistocene while among-clade divergence times (499.000-719,000ybp) were deeper. Genetic differentiation in D. baumanni and distributions of stonefly species were significantly concordant. These results imply that climatic oscillations have played major roles in shaping the genetic structure and distributions of Great Basin stoneflies, but that divergence among clades occurred much earlier than our late Pleistocence/early Holocene predictions.     

Dispersal differences predict population genetic structure in Mormon crickets

Research investigating the geographical context of speciation has primarily focused on abiotic factors such as the role of Pleistocene glacial cycles, or geotectonic events. Few study systems allow a direct comparison of how biological differences, such as dispersal behaviour, affect population genetic structure of organisms that were subdivided during the Pleistocene. Mormon crickets exist in solitary and gregarious 'phases', which broadly correspond with an east-west mtDNA division across the Rocky Mountains. Gregarious individuals form bands that can move up to 2 km daily. This study assessed whether population genetic structure results mainly from deep Pleistocene vicariance or if we can also detect more recent genetic patterns due to phase and dispersal differences superimposed on the older, deeper divisions. We found that separation in refugia was a more important influence on genetic divergence than phase, with the Rockies acting as a barrier that separated Mormon cricket populations into eastern and western refugia during Pleistocene glacial cycles. However, patterns of isolation by distance differ between eastern and western clades for both mitochondrial and nuclear DNA, with greater divergence within the eastern, solitary clade. An mtDNA haplotype mismatch distribution is compatible with historical population expansion in the western clade but not in the eastern clade. A persistent (and possibly sex-biased) difference in dispersal ability has most likely influenced the greater population genetic structure seen in the eastern clade, emphasizing the importance of the interaction of Quaternary climate fluctuations and geography with biotic factors in producing the patterns of genetic subdivision observed today.     

Genotyping‐by‐sequencing and ecological niche modeling illuminate phylogeography,admixture, and Pleistocene range dynamics in quaking aspen (Populus tremuloides)

Populus tremuloides is the widest‐ranging tree species in North America and an ecologically important component of mesic forest ecosystems displaced by the Pleistocene glaciations. Using phylogeographic analyses of genome‐wide SNPs (34,796 SNPs, 183 individuals) and ecological niche modeling, we inferred population structure, ploidy levels, admixture, and Pleistocene range dynamics of P. tremuloides, and tested several historical biogeographical hypotheses. We found three genetic lineages located mainly in coastal–Cascades (cluster 1), east‐slope Cascades–Sierra Nevadas–Northern Rockies (cluster 2), and U.S. Rocky Mountains through southern Canadian (cluster 3) regions of the P. tremuloides range, with tree graph relationships of the form ((cluster 1, cluster 2), cluster 3). Populations consisted mainly of diploids (86%) but also small numbers of triploids (12%) and tetraploids (1%), and ploidy did not adversely affect our genetic inferences. The main vector of admixture was from cluster 3 into cluster 2, with the admixture zone trending northwest through the Rocky Mountains along a recognized phenotypic cline (Utah to Idaho). Clusters 1 and 2 provided strong support for the “stable‐edge hypothesis” that unglaciated southwestern populations persisted in situ since the last glaciation. By contrast, despite a lack of clinal genetic variation, cluster 3 exhibited “trailing‐edge” dynamics from niche suitability predictions signifying complete northward postglacial expansion. Results were also consistent with the “inland dispersal hypothesis” predicting postglacial assembly of Pacific Northwestern forest ecosystems, but rejected the hypothesis that Pacific‐coastal populations were colonized during outburst flooding from glacial Lake Missoula. Overall, congruent patterns between our phylogeographic and ecological niche modeling results and fossil pollen data demonstrate complex mixtures of stable‐edge, refugial locations, and postglacial expansion within P. tremuloides. These findings confirm and refine previous genetic studies, while strongly supporting a distinct Pacific‐coastal genetic lineage of quaking aspen.     

Landscape characteristics influence morphological and genetic differentiation in a widespread raptor (Buteo jamaicensis)

Landscape-scale population genetic structure in vagile vertebrates was commonly considered to be a contradiction in terms whereas recent studies have demonstrated behaviour and habitat associated structure in several such species. We investigate whether landscape features influence morphological and genetic differentiation in a widespread, mobile raptor. To accurately describe genetic differentiation associated with regional landscape factors, we first investigated subspecies relationships at a continental scale. We used 17 microsatellite loci and five morphological measurements to investigate differentiation between eastern and western subspecies of red-tailed hawks ( Buteo jamaicensis ) and to identify patterns between differentiation and habitat within western North America. Bayesian and frequency-based analyses of microsatellite data revealed clear distinctions between B. j. borealis (eastern) and B. j. calurus (western) samples. Furthermore, hawks sampled in Texas were stouter than those collected from the Rocky Mountains and farther west. Among western samples, birds from the Great Basin, Rocky Mountains, and Washington were significantly different in morphology than those from Oregon and California. We identified a pattern of isolation by distance among western breeding sites around the Sierra Nevada. Given the long-range dispersal capabilities of raptors, this pattern suggests that population-specific habitat preferences, corresponding with habitat breaks between eastern and western slopes of the Sierra Nevada, and/or regionally variable population densities limit migration between the Mediterranean habitat of central California and the xeric habitats of southern California and interior west. We suggest habitat preferences and regionally disparate population densities may play a role in shaping genetic structure in vagile avian taxa.     

Cryptic Neogene vicariance and Quaternary dispersal of the red-spotted toad (Bufo punctatus): insights on the evolution of North American warm desert biotas

We define the geographical distributions of mitochondrial DNA (mtDNA) lineages embedded within a broadly distributed, arid-dwelling toad, Bufo punctatus. These patterns were evaluated as they relate to hypothesized vicariant events leading to the formation of desert biotas within western North America. We assessed mtDNA sequence variation among 191 samples from 82 sites located throughout much of the species' range. Parsimony-based haplotype networks of major identified lineages were used in nested clade analysis (NCA) to further elucidate and evaluate shallow phylogeographic patterns potentially associated with Quaternary (Pleistocene-Holocene) vicariance and dispersal. Phylogenetic analyses provided strong support for three monophyletic lineages (clades) within B. punctatus. The geographical distributions of the clades showed little overlap and corresponded to the general boundaries of the Peninsular Desert, and two continental desert regions, Eastern (Chihuahuan Desert-Colorado Plateau) and Western (Mojave-Sonoran deserts), geographically separated along the Rocky Mountains and Sierra Madre Occidental. The observed divergence levels and congruence with postulated events in earth history implicate a late Neogene (latest Miocene-early Pliocene) time frame for separation of the major mtDNA lineages. Evaluation of nucleotide and haplotype diversity and interpretations from NCA reveal that populations on the Colorado Plateau resulted from a recent, likely post-Pleistocene, range expansion from the Chihuahuan Desert. Dispersal across historical barriers separating major continental clades appear to be recent, resulting in secondary contacts in at least two areas. Given the observed contact between major clades, we speculated as to why the observed deep phylogeographic structure has not been eroded during the multiple previous interglacials of the Pleistocene.     

Identification of source‐sink dynamics in mountain lions of the Great Basin

Natural and anthropogenic boundaries have been shown to affect population dynamics and population structure for many species with movement patterns at the landscape level. Understanding population boundaries and movement rates in the field for species that are cryptic and occur at low densities is often extremely difficult and logistically prohibitive; however genetic techniques may offer insights that have previously been unattainable. We analysed thirteen microsatellite loci for 739 mountain lions (Puma concolor) using muscle tissue samples from individuals in the Great Basin throughout Nevada and the Sierra Nevada mountain range to test the hypothesis that heterogeneous hunting pressure results in source‐sink dynamics at the landscape scale. We used a combination of non‐spatial and spatial model‐based Bayesian clustering methods to identify genetic populations. We then used a recently developed Bayesian multilocus genotyping method to estimate asymmetrical rates of contemporary movement between those subpopulations and to identify source and sink populations. We identified two populations at the highest level of genetic structuring with a total of five subpopulations in the Great Basin of Nevada and the Sierra Nevada range. Our results suggest that source‐sink dynamics occur at landscape scales for wide‐ranging species, such as mountain lions, and that source populations may be those that are under relatively less hunting pressure and that occupy refugia.     

Near‐future forest vulnerability to drought and fire varies across the western United States

Recent prolonged droughts and catastrophic wildfires in the western United States have raised concerns about the potential for forest mortality to impact forest structure, forest ecosystem services, and the economic vitality of communities in the coming decades. We used the Community Land Model (CLM) to determine forest vulnerability to mortality from drought and fire by the year 2049. We modified CLM to represent 13 major forest types in the western United States and ran simulations at a 4‐km grid resolution, driven with climate projections from two general circulation models under one emissions scenario (RCP 8.5). We developed metrics of vulnerability to short‐term extreme and prolonged drought based on annual allocation to stem growth and net primary productivity. We calculated fire vulnerability based on changes in simulated future area burned relative to historical area burned. Simulated historical drought vulnerability was medium to high in areas with observations of recent drought‐related mortality. Comparisons of observed and simulated historical area burned indicate simulated future fire vulnerability could be underestimated by 3% in the Sierra Nevada and overestimated by 3% in the Rocky Mountains. Projections show that water‐limited forests in the Rocky Mountains, Southwest, and Great Basin regions will be the most vulnerable to future drought‐related mortality, and vulnerability to future fire will be highest in the Sierra Nevada and portions of the Rocky Mountains. High carbon‐density forests in the Pacific coast and western Cascades regions are projected to be the least vulnerable to either drought or fire. Importantly, differences in climate projections lead to only 1% of the domain with conflicting low and high vulnerability to fire and no area with conflicting drought vulnerability. Our drought vulnerability metrics could be incorporated as probabilistic mortality rates in earth system models, enabling more robust estimates of the feedbacks between the land and atmosphere over the 21st century.     

Pleistocene origin and population history of a neoendemic alpine butterfly

Alpine environments underwent dramatic transformation during glacial–interglacial cycles, with the consequence that geographical, ecological and demographic changes of alpine populations provided the opportunity for formation of neoendemic species. Several biogeographical models have been proposed to account for the unique history of alpine populations, with different expectations of genetic divergence and speciation. The expanding alpine archipelago model proposes that alpine populations expand spatially and demographically during glacial events, dispersing between mountain ranges. Under this model, alpine populations are unlikely to diverge in isolation due to substantial interpopulation gene flow. In contrast, the alpine archipelago refuge model proposes that gene flow during glacial phases is limited and populations expand demographically during interglacial phases, increasing genetic isolation and the likelihood of speciation. We assess these models by reconstructing the evolutionary history of Colias behrii, a morphologically and ecologically distinct alpine butterfly restricted to the California Sierra Nevada. C. behrii exhibits very low genetic diversity at mitochondrial and nuclear loci, limited population structure and evidence of population expansion. C. behrii and Rocky Mountain C. meadii share identical mitochondrial haplotypes, while in contrast, nuclear data indicate common ancestry between C. behrii and Cascades Range Colias pelidne. The conflict in gene genealogies may be a result of recent expansion in North American Colias, but an isolation with migration analysis indicates that genetic patterns in C. behrii might result from differential introgression following hybridization. Based on the timing of population expansion and gene flow between mountain ranges, the expanding alpine archipelago model is supported in C. behrii.