Genetic diversity, structure, and demographic change in tanoak, Lithocarpus densiflorus (Fagaceae), the most susceptible species to sudden oak death in California

Knowledge of population genetic structure of tanoak (Lithocarpus densiflorus) is of interest to pathologists seeking natural variation in resistance to sudden oak death disease, to resource managers who need indications of conservation priorities in this species now threatened by the introduced pathogen (Phytophthora ramorum), and to biologists with interests in demographic processes that have shaped plant populations. We investigated population genetic structure using nuclear and chloroplast DNA (cpDNA) and inferred the effects of past population demographic processes and contemporary gene flow. Our cpDNA results revealed a strong pattern of differentiation of four regional groups (coastal California, southern Oregon, Klamath mountains, and Sierra Nevada). The chloroplast haplotype phylogeny suggests relatively deep divergence of Sierra Nevada and Klamath populations from those of coastal California and southern Oregon. A widespread coastal California haplotype may have resulted from multiple refugial sites during the Last Glacial Maximum or from rapid recolonization from few refugia. Analysis of nuclear microsatellites suggests two major groups: (1) central coastal California and (2) Sierra Nevada/Klamath/southern Oregon and an area of admixture in north coastal California. The low level of nuclear differentiation is likely to be due to pollen gene flow among populations during postglacial range expansion.     

Limited phylogeographic structure in the flightless ground beetle, Calathus ruficollis, in southern California

The California Floristic Province is home to more than 8000 species of beetles, yet their geographical patterns of supra- and infraspecific diversity remain largely unexplored. In this paper, we investigate the phylogeography and population demographics of a flightless ground beetle, Calathus ruficollis (Coleoptera: Carabidae), in southern California. We sampled 136 specimens from 25 localities divided into 10 populations using a fragment of the mitochondrial cytochrome oxidase I gene. We tested several hypotheses, including the association of geography with particular clades and populations, the degree of differentiation among regions, and the expansion of populations. Parsimony and Bayesian phylogenetic analyses along with nested clade analysis and amova indicate a deep split between the southern Sierra Nevada population and populations south and west. This split corresponds closely to the split between subspecies C. ruficollis ignicollis (southern Sierra Nevada) and C. ruficollis ruficollis . Populations otherwise exhibit limited geographical structure, though F st values indicate some local differentiation. Mismatch distributions and Fu's F s indicate range expansion of several populations, suggesting that some structure may have been obscured by recent exchange. The population of C. ruficollis on Santa Cruz Island, which might have been expected to be isolated, shares several haplotypes with mainland populations, appearing to represent multiple colonizations.     

Testing hypotheses of Pleistocene population history using coalescent simulations: phylogeography of the pygmy nuthatch (Sitta pygmaea)

In this paper, we use mitochondrial NADH dehydrogenase subunit 2 sequences to test Pleistocene refugial hypotheses for the pygmy nuthatch (Sitta pygmaea). Pygmy nuthatches are a common resident of long-needle pine forests in western North America and demonstrate a particular affinity with ponderosa pine (Pinus ponderosa). Palaeoecological and genetic data indicate that ponderosa pine was isolated in two Pleistocene refugia corresponding to areas in the southern Sierra Nevada in the west and southern Arizona and New Mexico in the east. We use coalescent simulations to test the hypothesis that pygmy nuthatches tracked the Pleistocene history of their preferred habitat and persisted in two refugia during the periods of glacial maxima. Coalescent simulation of population history does not support the hypothesis of two Pleistocene refugia for the pygmy nuthatch. Instead, our data are consistent with a single refuge model. Nucleotide diversity is greatest in the western populations of southern and coastal California. We suggest that the pygmy nuthatch expanded from a far western glacial refuge into its current distribution since the most recent glacial maximum.     

Lineage divergence and speciation in the Web‐toed Salamanders (Plethodontidae: Hydromantes) of the Sierra Nevada,California

Peripatric speciation and the importance of founder effects have long been controversial, and multilocus sequence data and coalescent methods now allow hypotheses of peripatric speciation to be tested in a rigorous manner. Using a multilocus phylogeographical data set for two species of salamanders (genus Hydromantes) from the Sierra Nevada of California, hypotheses of recent divergence by peripatric speciation and older, allopatric divergence were tested. Phylogeographical analysis revealed two divergent lineages within Hydromantes platycephalus, which were estimated to have diverged in the Pliocene. By contrast, a low‐elevation species, Hydromantes brunus, diverged from within the northern lineage of H. platycephalus much more recently (mid‐Pleistocene), during a time of major climatic change in the Sierra Nevada. Multilocus species tree estimation and coalescent estimates of divergence time, migration rate, and growth rate reject a scenario of ancient speciation of H. brunus with subsequent gene flow and introgression from H. platycephalus, instead supporting a more recent divergence with population expansion. Although the small, peripheral distribution of H. brunus suggests the possibility of peripatric speciation, the estimated founding population size of the species was too large to have allowed founder effects to be important in its divergence. These results provide evidence for both recent speciation, most likely tied to the climatic changes of the Pleistocene, and older lineage divergence, possibly due to geological events, and add to evidence that Pleistocene glacial cycles were an important driver of diversification in the Sierra Nevada.     

Concordant molecular and phenotypic data delineate new taxonomy and conservation priorities for the endangered mountain yellow-legged frog

The mountain yellow-legged frog Rana muscosa sensu lato , once abundant in the Sierra Nevada of California and Nevada, and the disjunct Transverse Ranges of southern California, has declined precipitously throughout its range, even though most of its habitat is protected. The species is now extinct in Nevada and reduced to tiny remnants in southern California, where as a distinct population segment, it is classified as Endangered. Introduced predators (trout), air pollution and an infectious disease (chytridiomycosis) threaten remaining populations. A Bayesian analysis of 1901 base pairs of mitochondrial DNA confirms the presence of two deeply divergent clades that come into near contact in the Sierra Nevada. Morphological studies of museum specimens and analysis of acoustic data show that the two major mtDNA clades are readily differentiated phenotypically. Accordingly, we recognize two species, Rana sierrae , in the northern and central Sierra Nevada, and R. muscosa , in the southern Sierra Nevada and southern California. Existing data indicate no range overlap. These results have important implications for the conservation of these two species as they illuminate a profound mismatch between the current delineation of the distinct population segments (southern California vs. Sierra Nevada) and actual species boundaries. For example, our study finds that remnant populations of R. muscosa exist in both the southern Sierra Nevada and the mountains of southern California, which may broaden options for management. In addition, despite the fact that only the southern California populations are listed as Endangered, surveys conducted since 1995 at 225 historic (1899–1994) localities from museum collections show that 93.3% ( n =146) of R. sierrae populations and 95.2% ( n =79) of R. muscosa populations are extinct. Evidence presented here underscores the need for revision of protected population status to include both species throughout their ranges.     

Multilocus phylogeography of the flightless darkling beetle Nyctoporis carinata (Coleoptera: Tenebrionidae) in the California Floristic Province: deciphering an evolutionary mosaic

The California Floristic Province (CFP) is considered a global biodiversity hotspot because of its confluence of high species diversity across a wide range of threatened habitats. To understand how biodiversity hotspots such as the CFP maintain and generate diversity, we conducted a phylogeographic analysis of the flightless darkling beetle, Nyctoporis carinata, using multiple genetic markers. Analyses of both nuclear and mitochondrial loci revealed an east–west genetic break through the Transverse Ranges and high genetic diversity and isolation of the southern Sierra Nevada Mountains. Overall, the results obtained suggest that this species has a deep evolutionary history whose current distribution resulted from migration out of a glacial refugium in the southern Sierra Nevada via the Transverse Ranges. This finding is discussed in light of similar genetic patterns found in other taxa to develop a foundation for understanding the biodiversity patterns of this dynamic area. © 2010 The Linnean Society of London, Biological Journal of the Linnean Society, 2010, 99 , 424–444.     

Effects of historical demography and ecological context on spatial patterns of genetic diversity within foxtail pine (Pinus balfouriana; Pinaceae) stands located in the Klamath Mountains, California

The density and dispersion of individuals, nonequilibrium demographics, and habitat fragmentation all affect the magnitude and extent of spatial genetic structure within forest tree populations. Here, we investigate the link between historical demography and spatial genetic structure within ecologically contrasting stands of foxtail pine (Pinus balfouriana) in the Klamath Mountains of northern California. We defined two stand types a priori, based largely on differences in foxtail pine density and basal area, and for each type we sampled two stands. Population expansions, likely from Pleistocene bottlenecks, were detected in three of the four stands. The magnitude and extent of spatial autocorrelation among genotypes at five nuclear microsatellites differed dramatically among stands, with those having lower foxtail pine density exhibiting strong patterns of isolation by distance. Moran's I statistics were 7-fold higher for the first distance class (<25 m) in these stands relative to those observed in stands with higher foxtail pine density (I(25) = 0.14 vs. 0.02). We conclude that differences in spatial genetic structure between stand types are due to differences in ecological attributes that affected expansion from inferred bottlenecks.     

Demographic responses to climate‐driven variation in habitat quality across the annual cycle of a migratory bird species

The demography and dynamics of migratory bird populations depend on patterns of movement and habitat quality across the annual cycle. We leveraged archival GPS‐tagging data, climate data, remote‐sensed vegetation data, and bird‐banding data to better understand the dynamics of black‐headed grosbeak (Pheucticus melanocephalus) populations in two breeding regions, the coast and Central Valley of California (Coastal California) and the Sierra Nevada mountain range (Sierra Nevada), over 28 years (1992–2019). Drought conditions across the annual cycle and rainfall timing on the molting grounds influenced seasonal habitat characteristics, including vegetation greenness and phenology (maturity dates). We developed a novel integrated population model with population state informed by adult capture data, recruitment rates informed by age‐specific capture data and climate covariates, and survival rates informed by adult capture–mark–recapture data and climate covariates. Population size was relatively variable among years for Coastal California, where numbers of recruits and survivors were positively correlated, and years of population increase were largely driven by recruitment. In the Sierra Nevada, population size was more consistent and showed stronger evidence of population regulation (numbers of recruits and survivors negatively correlated). Neither region showed evidence of long‐term population trend. We found only weak support for most climate–demographic rate relationships. However, recruitment rates for the Coastal California region were higher when rainfall was relatively early on the molting grounds and when wintering grounds were relatively cool and wet. We suggest that our approach of integrating movement, climate, and demographic data within a novel modeling framework can provide a useful method for better understanding the dynamics of broadly distributed migratory species.     

Phylogeography of the Northern Alligator Lizard (Squamata,Anguidae): Hidden diversity in a western endemic

Western North America includes the California Floristic Province and the Pacific Northwest, biologically diverse regions highlighted by a complex topography, geology, climate and history. A number of animals span these regions and show distinctive patterns of dispersal, vicariance and lineage diversification. Examining phylogeographic patterns in the fauna of this area aids in our understanding of the forces that have contributed to the generation and maintenance of regional biodiversity. Here, we investigate the biogeography and population structure of the Northern Alligator Lizard (Elgaria coerulea), a wide‐ranging anguid endemic to western North America. We sequenced two mtDNA fragments (ND2 and ND4) for 181 individuals across the range of the species and analysed these data with phylogenetic approaches to infer population and biogeographic history, and date major divergences within the taxon. We further used Bayesian clustering methods to assess major patterns of population structure and performed ecological niche modelling (ENM) to aid in our interpretation of geographic structure and diversification of E. coerulea lineages. Our phylogeographic examination of E. coerulea uncovered surprising diversity and structure, recovering 10 major lineages, each with substantial geographic substructure. While some divergences within the species are relatively old (Pliocene, 5.3–2.6 mya), most intraspecific variation appears to be of more recent origin (Pleistocene, 2.6 mya‐11,700 ya). Current diversity appears to have arisen in the Sierra Nevada Mountains and spread west and north since the Pliocene. Finally, our ENMs suggest that much of the Coast Ranges in California provided ideal habitat during the Last Glacial Maxima (LGM) that has since contracted dramatically and shifted northwards, whereas significant portions of the Sierra Nevada were unsuitable during the LGM and have since become more suitable. Interestingly, E. coerulea shares a number of genetic boundaries with other sympatric taxa, suggesting common historical events and geomorphological features have shaped the biota of this region.     

Phylogeographic evidence for two mesic refugia in a biodiversity hotspot

Phylogeographic studies of flora in species-rich south-western Australia point to complex evolutionary histories, reflecting patterns of persistence and resilience to climatic changes during the Pleistocene. We asked whether coastal areas of the mid-west and south, as well as granite outcrops and inland ranges, have acted as major refugia within this region during Pleistocene climatic fluctuations by analysing phylogeographic patterns in the shrub Calothamnus quadrifidus R.Br. (Myrtaceae). We determined variation in chloroplast DNA data for 41 populations across the geographic range. Relationships and major clades were resolved using parsimony and Bayesian analyses. We tested for demographic and spatial expansion of the major clades and estimated clade divergence dates using an uncorrelated, lognormal relaxed clock based on two conservative chloroplast mutation rates. Two distinct phylogeographic clades were identified showing divergence during the Pleistocene, consistent with other phylogeographic studies of south-west Australian flora, emphasising the impact of climatic oscillations in driving divergence in this landscape. The southern clade was more diverse, having higher haplotype diversity and greater genetic structure, while the northern clade showed evidence of fluctuation in population size. Regions of high haplotype diversity with adjacent areas of low diversity observed in each clade indicated the locations of two coastal refugia: one on the south coast and another along the mid-west coast. This is the first evidence for major Pleistocene refugia using chloroplast genetic data in a common, widespread species from this region.     

Phylogeography and systematics of the Peromyscus eremicus species group and the historical biogeography of North American warm regional deserts

Phylogeographic relationships among 26 populations from throughout the geographic range of the Peromyscus eremicus species group are described based on sequence data for a 699-bp fragment of the mitochondrial DNA COIII gene. Distance, maximum-likelihood, and maximum-parsimony analyses of phylogenetic trees generated under four separate character-weighting strategies and representing five alternative biogeographic hypotheses revealed the existence of a cryptic species (Peromyscus fraterculus, previously included under P. eremicus) on the Baja California Peninsula and adjacent southwestern California and two distinct forms of P. eremicus, one from the Mojave, Sonoran, and northwestern Chihuahuan regional deserts (West) and one from the remainder of the Chihuahuan Desert (East). Distinctiveness of P. fraterculus is supported by previous morphometric and allozyme analyses, including comparisons with neighboring P. eremicus and parapatric P. eva, with which P. fraterculus shares a sister taxon relationship. Divergence of the eva + fraterculus, West + East eremicus, and P. merriami haplotype lineages likely occurred in the late Neogene (3 Ma), in response to northern extension of the Sea of Cortéz and elevation of the Sierra Madre Occidental; divergence of eva from fraterculus is concordant with the existence of a trans-Peninsular seaway during the Pleistocene (1 Ma); and divergence of West from East eremicus occurred during the Pleistocene pluvial-interpluvial cycles, but well before the Wisconsinan glacial interval. The sequence of divergence within the eremicus species group and causal association of geological events of the Neogene and Holocene provide a working hypothesis against which phylogeographic patterns among other arid-adapted species of the warm regional deserts of North America may be compared.     

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.     

Geographic partitioning of chloroplast DNA variation in the genusDatisca (Datiscaceae)

Datisca (Datiscaceae) is a ditypic genus with an intercontinentally disjunct distribution. Chloroplast DNA restriction site data was obtained from 23 populations and four 10–20 year old herbarium specimens ofD. glomerata and three populations ofD. cannabina from throughout their geographic ranges in western North America and southwest-central Asia, respectively. InD. glomerata, plastome diversity is partitioned geographically. All populations from southern California have a common plastome, while most populations north of this region share a relatively divergent plastome (0.49% sequence divergence). Likewise, these plastomes are highly divergent (0.87% mean sequence divergence) from those found inD. cannabina. Biogeographic processes dating to the Pleistocene and Late Miocene may be responsible for these intra- and interspecific patterns of chloroplast DNA divergence.     

Long-distance dispersal vs vicariance: the origin and genetic diversity of alpine plants in the Spanish Sierra Nevada

Here, we investigated the origin and genetic diversity of four alpine plant species co-occurring in the Spanish Sierra Nevada and other high mountains in south-western Europe by analysis of amplified fragment length polymorphisms (AFLPs). In Kernera saxatilis, Silene rupestris and Gentiana alpina we found intraspecific phylogroups corresponding to mountain regions as predicted by the vicariance hypothesis. Moreover, genetic distances between Sierra Nevada and Pyrenees populations were always higher than those between populations from the Pyrenees and the south-western Alps/Massif Central. This suggests successive disruption of gene exchange between mountain ranges as postglacial climatic warming proceeded from south to north. In Papaver alpinum, our data indicate that a central Pyrenean population arose via long-distance dispersal from the Sierra Nevada, and that vicariant separation events between the Sierra Nevada and the Pyrenees and between the Pyrenees and the south-western Alps occurred simultaneously. Overall, Sierra Nevada populations of all species investigated here preserve unexpectedly high (or not exceptionally reduced) genetic diversity. This testifies to the important influence of long-term isolation, i.e. vicariance, on genetic diversity through fostering the accumulation of new mutations and/or the fixation of ancestral ones.     

Testing the ‘Pleistocene species pump’ in alpine habitats: lineage diversification of flightless ground beetles (Coleoptera: Carabidae: Nebria) in relation to altitudinal zonation

Alpine species often have similar demographic responses to Pleistocene climate changes, but exhibit different spatial patterns of genetic diversity. Using a comparative phylogeographical approach, we examined the factors influencing lineage formation in three alpine carabid beetles of the genus Nebria Latreille inhabiting the California Sierra Nevada. These flightless beetles differ in altitudinal zonation and habitat preferences, but overlap spatially, have limited dispersal capacities and share life history characteristics. Species distribution modelling predicted decreasing population connectivity in relation to increasing altitudinal preferences. Diversity patterns at the cytochrome oxidase subunit I gene revealed north–south genetic structure and recent population growth in all three species. The high‐elevation‐restricted species, Nebria ingens Horn, exhibited a deep phylogeographical split, morphological divergence and evidence of limited, unidirectional gene flow towards the south. This was supported by additional data from three nuclear genes and isolation with migration analysis. Nebria spatulata Van Dyke, inhabiting an intermediate altitudinal range, exhibited fixed morphological differences between northern and southern populations, but showed limited structure. The broadly distributed Nebria ovipennis LeConte showed less structure and lacked morphological variation. Diversification of these Nebria species supports the role of altitudinal zonation in lineage formation and is consistent with the Pleistocene species pump model. © 2012 The Linnean Society of London, Biological Journal of the Linnean Society, 2012, ?? , ??–??.     

Multilocus phylogeographic assessment of the California Mountain Kingsnake (Lampropeltis zonata) suggests alternative patterns of diversification for the California Floristic Province

Phylogeographic inference can determine the timing of population divergence, historical demographic processes, patterns of migration, and when extended to multiple species, the history of communities. Single‐locus analyses can mislead interpretations of the evolutionary history of taxa and comparative analyses. It is therefore important to revisit previous single‐locus phylogeographic studies, particularly those that have been used to propose general patterns for regional biotas and the processes responsible for generating inferred patterns. Here, we employ a multilocus statistical approach to re‐examine the phylogeography of Lampropeltis zonata. Using nonparametic and Bayesian species delimitation, we determined that there are two well‐supported species within L. zonata. Ecological niche modelling supports the delimitation of these taxa, suggesting that the two species inhabit distinct climatic environments. Gene flow between the two taxa is low and appears to occur unidirectionally. Further, our data suggest that gene flow was mediated by females, a rare pattern in snakes. In contrast to previous analyses, we determined that the divergence between the two lineages occurred in the late Pliocene (c. 2.07 Ma). Spatially and temporally, the divergence of these lineages is associated with the inundation of central California by the Monterey Bay. The effective population sizes of the two species appear to have been unaffected by Pleistocene glaciation. Our increased sampling of loci for L. zonata, combined with previously published multilocus analyses of other sympatric species, suggests that previous conclusions reached by comparative phylogeographic studies conducted within the California Floristic Province should be reassessed.     

Geographic differentiation of monoterpenes from Abies procera and Abies magnifica

Cortical oleoresins from 352 Abies magnifica and A. procera trees, collected in 35 localities, were analysed for composition of their monoterpene fractions. The monoterpenes were composed primarily of -pinene, β-pinene, limonene and β-phellandrene, with camphene and 3-carene appearing only sporadically in larger amounts. Limonene, β-phellandrene, -pinene and 3-cerene were used to study the transitional populations in northern California and southern and central Oregon and the transitional-like populations in southern Sierra Nevada. The populations segregated latitudinally into three related clusters—above 44° (A.procera), between 44° and 40° (transitional), and below 40° (A.magnifica). A.magnifica from southern Sierra Nevada and A.magnifica from Mtn. Shasta differed in a number of parameters, with southern Sierra Nevada populations being chemically much closer to typical A. magnifica from central and northern Sierra Nevada. While monoterpene and morphological data did not allow a clear Interpretation of the status of the transitional and southern populations, the paleobotanical evidence favored the recent evolution over introgression. The desirability of additional studies was indicated.     

Historical and Contemporary DNA Indicate Fisher Decline and Isolation Occurred Prior to the European Settlement of California

Establishing if species contractions were the result of natural phenomena or human induced landscape changes is essential for managing natural populations. Fishers (Martes pennanti) in California occur in two geographically and genetically isolated populations in the northwestern mountains and southern Sierra Nevada. Their isolation is hypothesized to have resulted from a decline in abundance and distribution associated with European settlement in the 1800s. However, there is little evidence to establish that fisher occupied the area between the two extant populations at that time. We analyzed 10 microsatellite loci from 275 contemporary and 21 historical fisher samples (1880–1920) to evaluate the demographic history of fisher in California. We did not find any evidence of a recent (post-European) bottleneck in the northwestern population. In the southern Sierra Nevada, genetic subdivision within the population strongly influenced bottleneck tests. After accounting for genetic subdivision, we found a bottleneck signal only in the northern and central portions of the southern Sierra Nevada, indicating that the southernmost tip of these mountains may have acted as a refugium for fisher during the anthropogenic changes of the late 19^th and early 20^th centuries. Using a coalescent-based Bayesian analysis, we detected a 90% decline in effective population size and dated the time of decline to over a thousand years ago. We hypothesize that fisher distribution in California contracted to the two current population areas pre-European settlement, and that portions of the southern Sierra Nevada subsequently experienced another more recent bottleneck post-European settlement.     

Pseudotsuga menziesii chemical races of California and Oregon

Cortical oleoresin from more than 300 Douglas-fir trees sampled in over 30 localities throughout California and southern Oregon was analyzed for its monoterpene hydrocarbon composition. The results confirmed the existence of a distinct chemical race found in its pure unhybridized form throughout the Sierra Nevada as well as near Lompoc in the Purissima Hills of coastal California. Chemical features characterizing this race included near-absence of sabinene and presence of larger amounts of - and β-pinenes in the cortical essential oil. Elsewhere in California and southernmost Oregon, the Sierra Nevada race has extensively intergraded with the coastal Douglas-fir race from western Oregon, Washington, and British Columbia (var. menziesii), a race producing large amounts of sabinene and lesser amounts of - and β-pinenes. On the basis of chemical composition, the Sierra Nevada race was judged to be closer to the inland Douglas-fir races (var. glauca) than to the coastal race. Evolution of the Sierra Nevada race is explained by the emergence of the Great Basin arid regions during the PliocenePleistocene epochs which disjoined Western Douglas-fir populations from their eastern inland counterparts.     

Phylogeography of the mountain chickadee (Poecile gambeli): diversification, introgression, and expansion in response to Quaternary climate change

Since the late 1990s, molecular techniques have fuelled debate about the role of Pleistocene glacial cycles in structuring contemporary avian diversity in North America. The debate is still heated; however, there is widespread agreement that the Pleistocene glacial cycles forced the repeated contraction, fragmentation, and expansion of the North American biota. These demographic processes should leave genetic 'footprints' in modern descendants, suggesting that detailed population genetic studies of contemporary species provide the key to elucidating the impact of the late Quaternary (late Pleistocene-Holocene). We present an analysis of mitochondrial DNA (mtDNA) variation in the mountain chickadee (Poecile gambeli) in an attempt to examine the genetic evidence of the impact of the late Quaternary glacial cycles. Phylogenetic analyses reveal two strongly supported clades of P. gambeli: an Eastern Clade (Rocky Mountains and Great Basin) and a Western Clade (Sierra Nevada and Cascades). Post-glacial introgression is apparent between these two clades in the Mono Lake region of Central California. Within the Eastern Clade there is evidence of isolation-by-distance in the Rocky Mountain populations, and of limited gene flow into and around the Great Basin. Coalescent analysis of genetic variation in the Western Clade indicates that northern (Sierra Nevada/Cascades) and southern (Transverse/Peninsular Ranges) populations have been isolated and evolving independently for nearly 60,000 years.