Abstract

Introduction Estimated future climate scenarios can be used to predict where hotspots of endemism may occur over the next century, but life history, ecological and genetic traits will be important in informing the varying responses within myriad taxa. Essential to predicting the consequences of climate change to individual species will be an understanding of the factors that drive genetic structure within and among populations. Here, I review the factors that influence the genetic structure of plant species in California, but are applicable elsewhere; existing levels of genetic variation, life history and ecological characteristics will affect the ability of an individual taxon to persist in the presence of anthropogenic change.

Factors influencing the distribution of genetic variation Persistence in the face of climate change is likely determined by life history characteristics: dispersal ability, generation time, reproductive ability, degree of habitat specialization, plant–insect interactions, existing genetic diversity and availability of habitat or migration corridors. Existing levels of genetic diversity in plant populations vary based on a number of evolutionary scenarios that include endemism, expansion since the last glacial maximum, breeding system and current range sizes.

Regional priorities and examples A number of well-documented examples are provided from the California Floristic Province. Some predictions can be made for the responses of plant taxa to rapid environmental changes based on geographic position, evolutionary history, existing genetic variation, and ecological amplitude.

Conclusions, Solutions and Recommendations The prediction of how species will respond to climate change will require a synthesis drawing from population genetics, geography, palaeontology and ecology. The important integration of the historical factors that have shaped the distribution and existing genetic structure of California’s plant taxa will enable us to predict and prioritize the conservation of species and areas most likely to be impacted by rapid climate change, human disturbance and invasive species.

INTRODUCTION

The extent to which high levels of genetic variation within and among populations are associated with the ability of a species to persist in the face of both short- and long-term disturbance regimes has long been deliberated (Lowe and Allendorf, 2010; Davies et al., 2016). Understanding the distribution of genetic variation on the landscape and the life history or ecological parameters that influence this variation will be important in identifying the taxa most vulnerable to anthropogenic change (Peter and Slatkin, 2013; Brandvain et al., 2014). High levels of heterozygosity have been repeatedly shown to confer resistance to environmental change (Hanski et al., 2006; Willi et al., 2006; Bonin et al., 2007). Model simulations suggest that species with low heterozygosity but large ranges will suffer a spatial redistribution of heterozygosity, depending on the level and length of the disturbance, whereas taxa with limited distributions with severe and frequent disturbance will suffer dramatic declines in heterozygosity regardless of initial levels of heterozygosity (Davies et al., 2016).

Here I focus primarily on molecularly-based population genetic structure and the life history characteristics, historical landscapes and climatic events that have shaped, and will shape, the distribution of genetic variation within and among California plant populations. California, specifically the California Floristic Province (CFP), is a recognized biodiversity hotspot (Myers et al., 2000) with high levels of endemism (Loarie et al., 2008; Kraft et al., 2010; Burge et al., 2016) that has been, and will continue to be, impacted by anthropogenic change. Various models have predicted different outcomes in climatic shifts and at this juncture, because it is difficult to predict how whole plant communities may respond, it will be important to focus on the adaptive ability of individual species and conserving regions that contain high levels of diversity. Genetic diversity as measured by neutral markers does not necessarily provide a measure of adaptive ability, but the study of neutral genetic markers can be useful in the interpretation of past landscapes, refugia and gene flow (Holderegger et al., 2006). An examination of fitness and heterozygosity via a meta-analysis of multiple datasets substantiates that a loss of heterozygosity as measured by neutral markers does indeed reflect a negative impact on fitness (Reed and Frankham, 2003). However, genes under selection are preferred for interpreting how populations may respond to climate changes as neutral genetic markers can lead to incorrect interpretations, particularly in highly fragmented populations (Hoffman and Willi, 2008). High levels of genetic diversity within a taxon dispersed across a diverse array of climate zones and habitat types should confer resilience to varying temperatures. Taxa with low or fragmented patterns of genetic diversity, limited geographic ranges or high habitat specificity are those most in danger of population extirpation or extinction loss due to anthropogenic change. Based on a review of the distribution of genetic variation of California plant taxa, examples of existing and potential refugia under climate change are also identified and discussed. Although the examples here are drawn primarily from the California Floristic Province, the conclusions and recommendations are broadly applicable, particularly in areas with high topographic diversity.

Molecularly-based population genetic data have been very useful in determining past historical events and will be useful in estimating the ability of taxa to respond to contemporary climate change. Extant species that have survived via glacial refugia are likely to show genetic signatures of within-refugium genetic drift and selection or, conversely, harbour historical genetic diversity not found throughout the range of a taxon; however, combinations of these scenarios may also occur. Even where phenotypic variation is present, measuring genetic variability due to recent divergence can be a problem, particularly in places like California, where speciation and evolutionary diversification have occurred recently (Boykin et al., 2005; Burge et al., 2011). Detection of patterns of recent genetic divergence now can be achieved, however, through the use of markers more suitable for measuring population-level variation, such as microsatellites, restriction site-associated DNA markers (RADseq) (Davey and Blaxter 2010) and whole-genome sequences (Brandvain et al., 2014).

Varying levels of diversity among refugia can help infer historical ecological patterns (Gugger et al., 2010) and may be used to identify those areas most valuable for conservation. Depending on the temporal scale, dispersal events, bottlenecks followed by population growth and gene flow among populations will all leave signatures of genetic variation (Slatkin and Hudson, 1991). The efficacy of contemporary refugia in conserving taxa will depend on the capacity of refugia to support enough genetic diversity for population viability and to maintain evolutionary processes (Ovaskainen, 2002). The topographical and latitudinal variation present in California will require the regional identification of areas that are spatially appropriate in size and proximity to the current ranges of the maximum number of species. Those areas that have been identified as containing high levels of species diversity and high genetic diversity within and among populations should be prioritized for conservation.

Connectivity associated with gene flow is further influenced by population demographics that contribute to patterns of genetic variation on physical landscapes; however, evolutionary genetic models more strongly support dynamics among populations than demographic models in Clarkia xantiana spp. xantiana (Moeller et al., 2011; Davies et al., 2016). Population differentiation is easier to resolve in species with low rates of dispersal due to high genetic structuring (high FST values), whereas high levels of gene flow (high Nm) can increase the difficulty of interpretation and will be reflected by low FST values. Divergent taxa residing within the same region with similar ecological requirements, life histories and behaviours often share similar patterns of genetic diversity (Bermingham and Moritz, 1998). However, loci under selection can alter expected patterns based on demographic patterns or the historical distribution of a population. Thus, precise identification of refugia and vicariant events requires data from a number of species, as differences in life history characteristics and ecological requirements affect lineages differently in time and space. Ultimately, understanding the genetic structures of populations and how they have been influenced by vicariant and dispersal events is important for understanding the potential of fragmented populations to adapt to anthropogenic change. The recognition of regional patterns of genetic diversity offers the opportunity not only to conserve regions with the highest biological diversity but also to conserve evolutionary processes in the face of climate change.

Range shifts of species in response to climate change are a serious concern (Warren et al., 2001; Zacherl et al., 2003). Mean annual temperatures in California increased by 1 °C between 1950 and 2000 (LaDochy et al., 2007) and are projected to rise 3·8–5·8 °C by 2100 under high emission scenarios by 2050 (Cayan et al., 2008). Vulnerability to climate change cannot be predicted by geographic range and habitat specificity, based on a lack of correlation of Climate Change Vulnerability Index scores and rarity type (Anacker et al., 2013); however, geographic range, connectivity and habitat specificity all influence genetic variation. The velocity of temperature change is estimated globally to be 0·42 km year−1 under some IPCC scenarios, a rate that most plants cannot meet, with the highest velocities in flooded grasslands and deserts and the lowest in temperate montane regions (Loarie et al., 2008). An argument has been made to prioritize conservation in regions with high levels of evolutionary diversity to capitalize on genetic and functional diversity (Cadotte and Davies, 2010). Measures of evolutionary diversity need to be developed to maximize ecosystem health and minimize losses to biological diversity. Unfortunately, genetic diversity is often ignored in climate change decision-making (Cadotte and Davies, 2010). Montane areas will likely provide refuge for many species and the creation of large desert reserves may reduce losses of xeric taxa (Loarie et al., 2008). Small reserves are problematic in that they provide only temporary protection. The question remains as to how to mitigate for the loss of habitat in California for 4844 native plant species (Baldwin et al., 2012) in over 435 vegetation alliances (Sawyer et al., 2009), in a rapidly changing climate on a landscape constantly disturbed by humans. Whether we can effectively estimate the ability of a species to migrate based on past migrational events further complicated by human landscape disturbance remains to be seen. Regardless, a combination of ecological and genetic data will provide the best foundation for conserving the ability of the California biota to respond to rapid environmental change (Eckert, 2011; Gugger et al., 2011). In the following sections I review the factors that have shaped the extant genetic structure of plant species in California and argue that existing levels of genetic variation, combined with life history and ecological characteristics, affect the ability of an individual taxon to persist in the presence of anthropogenic change. Of particular importance are dispersal ability, generation time, reproductive ability, degree of habitat specialization, plant–insect interactions and the availability of habitat or migration corridors.

FACTORS INFLUENCING THE DISTRIBUTION OF GENETIC VARIATION

Habitat specificity

Little is known about the interaction of range shifts, habitat specificity and genetic diversity. Dispersal and habitat limitation studies in Mimulus leptaleus and Mimuluslaciniatus in the California Sierra Nevada found that habitat limitation was more important than dispersal ability in limiting range expansion (Sexton and Dickman, 2016). However, in some cases there may be a lack of overlap between a species’ future latitudinal or elevational range and the habitat needed to support it (Bellard et al., 2012). Genetic structuring among populations of serpentine endemics can be significant and, as with many edaphic endemics, the level of genetic variation contained within each outcrop or island can have important conservation implications (Wolf et al., 2000). If there is little variation contained in a small population and the population becomes isolated due to the extirpation of adjacent populations, it may succumb to processes such as genetic drift and inbreeding depression. In a resource-limited world, populations that are genetically depauperate and spatially isolated may not be compelling for other than ex situ conservation. In contrast, populations with high levels of variation in both small and large populations are more viable for conservation (Lesica and Allendorf, 1995).

The California Coast Ranges were formed from a variety of geological substrates, which further complicates the inference of the evolutionary history of many California plant taxa. The late Cenozoic formation of much of the Coast Ranges has resulted in low levels of genetic differentiation among plant species, as found in Ceanothus (Burge et al., 2011) and Collinsia (Baldwin et al., 2011). Despite this recent evolutionary divergence, among-population genetic variation can be high within these species if they are endemic to specific habitats. Increasingly dry summer periods in California beginning about 15 million years ago (MYA) (Baldwin and Sanderson, 1998) may have resulted in intense selection via the restriction of some species populations to the fine-grained, high moisture capacity serpentine soils which often contain seeps (Raven and Axelrod, 1978). Without experimentation for adaptive variation, it will be difficult to prioritize the conservation of populations. For example, 32 Coast Range populations of Calystegia collina, a serpentine endemic, were found to have low gene flow among populations in Napa, Lake and Sonoma counties. Gene flow is thought to be primarily mediated by pollinators, which, should they decline, will result in even lower gene flow (Preston and Dempster, 2012). The maintenance of gene flow and genetic diversity will be important in maximizing the ability of populations to respond to evolutionary processes.

Plant–insect interactions

Insects are important conduits of gene flow in many plant species, maintaining within- and among-population levels of diversity. Insects are particularly susceptible to temperature and the availability of food resources, primarily the plants with which they evolved (Pelini et al., 2009). However, because insects and plants disperse and reproduce at varying rates their respective responses to climate change could become a problem for either group (Schweiger et al., 2008). Butterfly species have demonstrated migration to higher elevations or latitudes of 240 km/30 years (Parmesan et al., 1999), whereas tree migration responds to climate change at 20–40 km/100 years (Davis and Shaw, 2001). The context of the plant–animal interaction also needs to be taken into consideration, specifically the strength of the relationship or the availability of alternative host plants or pollinators. Differential phenological shifts among plants and their pollinators decrease the diet breadth of pollinators via the loss of floral resources and are predicted to result in increased inbreeding and eventual population declines in the plant hosts (Memmott et al., 2007). For example, the butterfly species Erynnis propertius, an important generalist pollinator, occurs on the Pacific coast from Baja California to British Columbia, where the host plants, species of Quercus, occur (Guppy and Shepard, 2001; Pelini et al., 2009). Southern and central California populations of E. propertius feed on Quercusagrifolia and higher-elevation populations feed on Quercuskelloggii. Host plant experiments have found that E. propertius have zero survival on other Quercus species. Southwestern Oregon populations are dependent on Quercusgarryana; however, overlaps between Q. agrifolia and Q. garryana are very limited. If E. propertius has reduced fitness on alternative host plants and the host plants migrate at much slower rates than E. propertius, its survival could be prevented if it is unable to switch hosts (Pelini et al., 2009). The loss of generalist pollinators such as E. propertius at a large scale has the potential to be devastating to plant taxa dependent on outcrossing. Specialist pollinators are expected to be the most vulnerable to extinction when plant–pollinator phenologies are decoupled, but further experimentation is necessary to make more accurate predictions (Forrest, 2015).

Comparative phylogeographic studies among Yucca brevifolia, two obligate pollinators (Tegeticula antithetica and Tegeticulasynthetica) and two species of parasitizing yucca moths (Prodoxus weethumpi and Prodoxussordidus) provides evidence for the effects of climate change on species distributions and genetic structure (Smith et al., 2011). Coalescent-based DNA sequence analyses and distribution modelling of Y. brevifolia trees since the Last Glacial Maximum (LGM) indicate that demographic changes in all species were present prior to climatic change during the Holocene (Smith et al., 2011). For the insect species some divergence is apparent, with average FST values of 0·12; however, Y. brevifolia variation is quite different, with high global FST values of 0·87, which strongly correlate with geographic distance, which is structured by pollinator association. The differences in FST values between Y. brevifolia and their associated insects likely reflect differences in dispersal ability. Demographic and genetic analyses indicate that all of the species had concerted population growth in the late Pleistocene beginning about 0·2 MYA, were little affected by the LGM, and reached their current distribution about 50 thousand years ago (KYA). These populations are projected to contract with increased CO2 levels and increased temperatures; however, some new habitats may become available (Dole et al., 2003). Decoupling plant–insect interactions via phenological changes or changes in distribution has the potential to become severe enough to result in extinction events of insect and plant taxa or both (Gilman et al., 2012).

Connectivity and disturbance

Populations residing in refugia from the Pleistocene glaciation cycles harbour much genetic diversity and are often found in cooler and more northern microclimates. However, low latitudinal or elevational portions and disjunct populations of a species’ range can also be important reservoirs of adaptive genetic variation (Hampe and Petit, 2005). Although these latter populations may not be currently connected, they likely persisted under both cold and warm conditions and are important sources of genetic distinctiveness, i.e. they may contain genetic diversity not found in other populations.

Increased human settlement can result in less connectivity and increases in regional genetic structure via increased population differentiation. Due to rapid environmental change, this increase in population differentiation is likely to be based on random effects such as genetic drift rather than selection due to sudden declines in population size. The resulting non-adaptive variation will lead to greater susceptibility of a population to further environmental change. The disruption of previously contiguous habitats can result in inbreeding depression and the artificial introduction of genotypes from other locations may be maladaptive (Ellstrand and Elam, 1993; Elam et al., 1998). Alternatively, artificially introduced gene flow can have negative consequences due to differential local adaptations and lead to outbreeding depression (Templeton, 1986; Sexton et al., 2011; Halbur et al., 2014). Generally, if gene flow levels have been historically high, the introduction from similar populations is unlikely to be harmful (Elam et al., 1998). Additionally, gene flow can be beneficial where range limits among populations are similar, and for populations at risk gene flow may be recommended as a conservation option (Sexton et al., 2011).

The extreme loss of mesic and wetland habitats throughout California, but particularly in the Central Valley, has isolated populations and removed opportunities for adaptive radiation and gene flow. For example, formerly widespread but habitat-specific taxa such as those that occur in vernal pools would have historically experienced occasional recolonization events due to periodic widespread flooding that would have resulted in gene flow among existing populations. Chloropyron palmatum, a federal endangered species and vernal pool endemic, is found only in five genetically distinct populations (FST = 0·23) in the Central Valley. Seed dispersal historically occurred via flooding, which, due to development, now rarely occurs and is thought to have reduced gene flow among populations (Ayres et al., 2015). Similarly, Allium munzii, which has suffered 80–90 % loss of its coastal sage scrub habitat, primarily due to urbanization, has very strong genetic structure and low gene flow among the remaining 18 small, isolated populations (Mashayekhi and Columbus, 2015). Species such as C. palmatum and A. munzii are unlikely to survive climate change without ex situ conservation. Unfortunately, assisted migration as a remedy for climate change has both risks and benefits, and translocating genotypes requires caution (Frankham, 2010; Vitt et al., 2010).

Climate alone does not predict the distribution of plant species, and even at large scales there may be other factors that contribute to the range of a species (Serra-Diaz et al., 2013). In addition to biotic interactions and traits such as dispersal ability, disturbance can be a major factor, changing light and soil conditions and influencing the demographics of a population (Serra-Diaz et al., 2015). Interacting with climate change, disturbance regimes and land management can hinder or benefit recruitment (Boulangeat et al., 2014). Diverse, intact landscapes, particularly in mountainous areas, may provide suitable habitat for many taxa, due to varying microrefugial climatic conditions, such as found on north-facing slopes (Vanderwel and Purves, 2014). In general, human-caused disturbances such as increased fire frequency will likely influence species distributions more than climate change, but understanding the interaction between disturbance and climate change will be important for the persistence of not just individual species but also of ecosystems (Penman et al., 2014; Syphard and Keeley, 2015).

Human activities, such as introducing non-native species and disturbing, fragmenting or otherwise degrading natural ecosystems, promote hybridization of previously allopatric congeners (Vilà et al., 2000). When introgression accompanies hybridization, dilution of the native gene pool can occur (Abbott, 1992; Johnson JR et al., 2010; Johnson MG et al., 2016), and can result in extinction of the native species (Anttila et al., 1998). Often underrated, human-mediated hybridization is one of the leading causes of biodiversity loss via homogenization of formerly distinct populations or taxa (Muhlfeld et al., 2009). Through introgressive hybridization, populations become less differentiated, which can limit evolutionary flexibility as locally adapted genes become lost along with the populations that once contained them (Rieseberg, 1991; Ellstrand and Elam, 1993). Conversely, hybridization between previously isolated populations may, in some cases, confer genetic resilience via new gene combinations that may be necessary to survive climate change. For example, populations isolated due to human disturbance may benefit from hybridization from populations in more southern latitudes.

Dispersal ability

Life history traits strongly influence dispersal distance (Stevens et al., 2012; Tamme et al., 2014) and consequently influence the distribution of genetic variation on the landscape. Seed and fruit mass and morphology will greatly influence dispersal distance depending on the vector (Thomson et al., 2010), but plant height has generally been found to be more important than seed mass in predicting dispersal distance (Moles et al., 2004; Tamme et al., 2014). Wind- and vertebrate-dispersed species, not surprisingly, disperse further than insect-dispersed species and generally have higher levels of gene flow and genetic similarity among populations (Hamrick et al., 1979), conferring some resilience to climate change depending on habitat availability. It is important to consider, however, how heterogeneous environments also influence dispersal ability. The highest levels of genetic diversity of Notholithocarpus densiflorus (four cpDNA haplotypes) are in southern Oregon, the Klamath Ranges and the northern Sierra Nevada, with one haplotype in the Coast Ranges from Humboldt to Santa Barbara counties (Nettel et al., 2009). The distribution of these haplotypes is likely a reflection of a complex terrain combined with heavy fruits that limit seed dispersal. In contrast to the cpDNA data, the nuclear genome reflects coastal groupings south of San Francisco and in an interior northern Sierra Nevada, Klamath and Southern Oregon group, which indicates that genetic mixing is more effective through pollen movement than through seed (Nettel et al., 2009). The occurrence of this species in complex topography combined with genetic mixing via wind pollination may facilitate its survival under climate change.

Species lifespan and reproductive rates

Species that have long lives have a different vulnerability to extinction than species that have short lives. Regardless of the type of lifespan, critical to survival and recolonization is the availability of migration corridors and source populations at a number of spatial scales. Tree species are long-lived and often have high effective population sizes and high rates of reproduction; along with delayed seed production, they often show genetic signatures of expansion and contraction and lag behind climate change. Other considerations affecting response to climate include a number of factors that influence reproduction rate, such as demographics within a population, variation of seed production within the range, sensitivity of germination to climate and biotic factors, and the availability of suitable microsites across the range of the species (Kroiss and Hillerislambers, 2015). Testing the hypothesis that climate change affects reproductive outcomes, seedling to mature tree size classes were used as a surrogate for age distribution for 46 tree species from northern California through the Pacific Northwest. Seedlings of 33 species occurred in ranges with a mean annual temperature of 0·120  °C lower than the mature trees, a difference consistent with climate change (Monleon and Lintz, 2015). Interestingly, Pinus lambertiana, Pinus jeffreyi, Calocedrus decurrens and Abies concolor, common components of montane regions of California, showed statistically significant shifts towards warmer areas.

Short-lived taxa are expected to respond more rapidly to climate change than the long-lived tree species. Encelia farinosa is a short-lived perennial that occurs in the Sonoran, Mojave and Peninsular Deserts and shows genetic signatures that reflect migration from refugial populations (Fehlberg and Ranker, 2009). A sample of cpDNA sequences reveals that 80·64 % of the variation is explained by differences within locations, consistent with past fragmentation and range expansion. Fossil records from middens for at least 13·8 KYA suggest that Encelia expanded west and north following the LGM. In general, short-lived species with effective regional dispersal, such as E. farinosa, will likely survive in microclimatic areas in adjoining mountainous regions in the face of rapid environmental change.

Two short-lived desert species, Cryptantha flava and Carrichtera annua, were found to be able to buffer stochastic changes under climate change due to elasticity in seed bank emergence (Salguero-Gómez et al., 2012), and this is consistent with model predictions (Koons et al., 2009). Long-lived species may have lower annual reproductive rates but their capacity to withstand incremental changes to their environment may be buffered by wider demographic reaction norms, i.e. greater phenotypic plasticity within different demographic stages (Morris et al., 2008; Koons et al., 2009; Chevin et al., 2013). A number of taxa indicate a non-linear, species-specific response to climate change that will require experimental assessment on an individual basis (Morris et al., 2008; Doak and Morris, 2010).

Existing genetic diversity

There are extreme examples of small isolated populations with low levels of diversity within and among populations (Syring et al., 2007). The long-lived Pinus torreyana, known from two populations with about 9000 individuals, has very low levels of genetic diversity. An analysis of 59 isozyme loci identified all as homozygous, suggested to be a result of reduction to less than 50 individuals during the xerothermic period 8·5–3·5 KYA (Ledig and Conkle, 1983). Pinus balfouriana occurs at 1830–3400 m in the southern Sierra Nevada and the Klamath Ranges, separated by 500 km. There is some morphological divergence but little allozyme divergence among the populations (Mastrogiuseppe and Mastrogiuseppe, 1980; Oline et al., 2000). An analysis of cpDNA, mtDNA and nucDNA genes or regions indicates divergence during the Sherwin glacial maximum event during the Pleistocene, extirpating populations in the central and northern Sierra Nevada, but the genetic data are not consistent with the Holocene xerotherm (Eckert and Carstens, 2008). Pinus balfouriana shows greater genetic diversity than other range-restricted California conifers and, although at great risk in the White Mountains due to climate change, it should survive in some populations in the Rocky Mountains. Another long-lived species, Pinus longaeva, once widespread throughout the Great Basin during the LGM, is now limited to alpine habitat, and the White Mountains represents the western extent of the species distribution. Despite the fact that they grow in isolated alpine habitats, they retain high levels of genetic diversity because of the life history characteristics of longevity, large population size and high seed set, and no evidence of genetic bottlenecks (Lee et al., 2002). However, due to their limited number of populations in a habitat that is predicted not to exist in the next century, their survival in situ is dubious (Dirnböck et al., 2010; Chen et al., 2011).

Hesperocyparis macrocarpa consists of small, isolated populations but has fairly high levels of genetic diversity, indicating a former historic, widespread distribution (Kafton, 1976). Although widely planted as an ornamental, the remaining native populations, although protected, are in urban areas along the coast, which are prone to human disturbance. Further placing this species at risk in situ will be the reduction of coastal fog projected to occur as a result of climate change (Farjon, 2013). Hesperocyparis forbesii, an endemic in the southern California Peninsular Ranges and northwest Baja California, has widely separated native populations; however, genetic diversity as measured by allozymes was not associated with geographic distance (Truesdale and McClenagan, 1998). A close relative of H. forbesii, Hesperocyparisguadalupensis, is found only on Guadalupe Island off the Pacific coast of Baja California. Paternally inherited cpDNA markers compared genetic diversity among populations of H. forbesii and H. guadalupensis and identified higher genetic diversity in H. guadalupensis than in its mainland congener. Contrary to expectations, H. guadalupensis shows no evidence of a genetic bottleneck and may reflect the environmental stability of the island habitat over a long evolutionary history. In contrast, other species support a pattern indicating that the mainland and more northern Peninsular Ranges experienced glacial and interglacial periods during the Pleistocene, which may have led to less predictable interactions within and among populations (Escobar García et al., 2012).

REGIONAL PRIORITIES AND EXAMPLES

It would be ideal to conserve as much ecological, species and genetic diversity as possible; however, climate change and the addition of 14 million people to California’s population by 2060 (State of California, 2014) will result in significant losses. Conservation priorities may have to be focused on areas that contain the greatest number of endemic species or high biodiversity, e.g. the central Coast Ranges, the Klamath–Siskiyou region, the eastern Transverse Ranges, the Tehachapi Mountains, the Santa Lucia Mountains and the eastern Mojave Desert. Areas with high levels of endemism that include species with the greatest range restrictions are on the central coast, in the Sierra Nevada and in the San Bernardino Mountains (Loarie et al., 2008; Kraft et al., 2010). The Mojave Desert and Great Basin have the highest level of neoendemics, or species that are range-restricted, in part because of their recent evolutionary emergence. Further complicating the diversity of California, however, is the presence of a number of palaeoendemic taxa, or those that were once widespread and are now range-restricted. Issues to consider in conservation include local, regional and taxonomic rarity, ecological function and threats to biodiversity at multiple levels.

Klamath–Siskiyou region

The International Union for Conservation of Nature (IUCN) has identified the Klamath–Siskiyou ecoregion as an ‘area of global botanical significance’ (DellaSala et al., 1999; Villa-Lobos, 2003). The biodiversity of this region is partially a result of species from the north reaching their southern limit and species from the south reaching their northern limit of distribution. Geological complexity, especially in those areas that provide edaphic soils, further contributes to the botanical diversity of the area. This area escaped oceanic submergence, major glaciation effects and volcanic activity throughout its history. A fairly stable climate and geographic complexity likely buffered plant community change and population loss for over 100 million years (Whittaker, 1961; Villa-Lobos, 2003). The characteristics that will continue to make the Klamath–Siskiyou ecoregion an important refugial area are topographic diversity, an array of soil types and high levels of biological diversity.

Studies of genetic variation within and among populations of plant species in northern California are primarily focused on conifer taxa due to their economic importance. The Klamath Ranges contain forests that are considered roughly equivalent to the Paleogene circumboreal forests (Sawyer and Thornburgh, 1977). Currently, the Klamath–Siskiyou region contains more than 281 endemic plant species and 3500 plant species (Sawyer, 2006). Pleistocene glaciations throughout California were significant forces affecting intraspecific, interspecific and intergeneric evolutionary patterns in plants and animals, particularly those with low dispersal ability. In contrast, the Klamath–Siskiyou region experienced minimal glaciation and served as a refugium for many plant taxa (Smith and Sawyer, 1988; Soltis et al., 1997; Sawyer, 2006).

Pleistocene glaciations are increasingly recognized as important in shaping the genetic structure of California plant species, as has been demonstrated in Pinus albicaulis, P. lambertiana (Liston et al., 2007) and Pinusmonticola (Steinhoff et al., 1990). Pinus albicaulis is likely a post-diversification migrant from Europe across Beringia, which was unglaciated during the Pleistocene and LGM, and is a known refugium (Brubaker et al., 2005). Pinus albicaulis shows little differentiation between populations, even though it has a discontinuous distribution in alpine habitat (Jorgensen and Hamrick, 1997). Both P. albicaulis and P. lambertiana demonstrate divergence into northern and southern clades, with secondary contact zones in the Klamath–Siskiyou Mountains (Eckert and Hall, 2006). Pinus lambertiana occurs from Oregon south to northern Baja California and individuals from Oregon, and the Klamath and the north coast ranges of California share a common cpDNA haplotype (‘N’), with plants from the Sierra Nevada and the Transverse and Peninsular Ranges having a different haplotype (‘S’) (Gernandt et al., 2005). There is a contact zone between the two haplotypes in northeastern California in a region about 150 km from Mount Lassen. Recent secondary contact suggests connections from northern and southern refugia following Holocene glaciations. The southern Cascades/Klamath–Siskiyou region is proposed to be a northern glacial refugium for P. lambertiana (Eckert and Hall, 2006; Syring et al., 2007). As climates warm, those conifer species with high genetic diversity, ranges that span the largest latitudes, and the broadest ecological amplitude may fare well, such as P.lambertiana, but the outlook is more tenuous for high-elevation species such as P. monticola, P. longaeva, and P. albicaulis. Pinus albicaulis populations were recently found to have a range of seedlings 90·0 m lower in altitude but 91·8 km north in latitude compared with the range of mature trees; mean annual temperature of the populations is similar and consistent with a gradual shift north due to climate change (Monleon and Lintz, 2015).

Patterns of genetic diversity in California conifers indicate a complex array of palaeoendemics, neoendemics that are variously uncommon or widespread. Generally, California conifers generally maintain high levels of genetic diversity, and many of these species, such as Abies concolor, Pinus ponderosa and Pseudotsuga menziesii, are likely to persist, but there are many notable exceptions. Neotoma spp. midden, palynological and genetic data support two refugia for P.ponderosa during the Pleistocene, the southern Sierra Nevada and the southern Rocky Mountains (Betancourt et al., 1990; Latta and Mitton, 1999). Due to high levels of genetic diversity and a large range, P. ponderosa populations are expected to contract but likely survive climate change, whereas even with high levels of diversity those species with more limited ranges are more precarious (Jump and Peñuelas, 2005).

Picea breweriana is currently found only in the Klamath Ranges but was formerly widespread in Idaho, Nevada, southern California and central Oregon in the Pliocene and Miocene. It now occurs in roughly the north–south extent of the Klamath region at elevations of 560–2300 m in small, disjunct populations. Isozyme loci show a lack of correspondence between genetic and geographic distance that is likely a reflection of some genetic drift (Ledig et al., 2005). Levels of heterozygosity and migration between populations are lower than for other conifers (Hamrick et al., 1992) but also decrease with latitudinal increases, which are estimated to be due to postglacial dispersal (Hamrick and Godt, 1996; Ledig, 2000; Ledig et al., 2005). The response of species such as P. breweriana is equivocal; because it is long-lived it may not be able to disperse rapidly enough to respond to warming climates, but due to the microclimates available within the complex topography of the Klamath Ranges it may survive.

It will be difficult to predict, even with large-scale experimentation or measures of adaptive genetic variation, how species with large ranges but low genetic diversity will fare in the face of climate change. For example, Cornus nuttallii is endemic to North America with a distribution from the southern California mountains to southwestern British Columbia. Extensive sampling from throughout the range of the species indicates that there is no isolation by distance and there is generally low diversity range-wide for both the nuclear and cpDNA genomes. A decline in diversity with latitude suggests there was long-distance dispersal from a southern refugium following glaciation. The presence of two chloroplast haplotypes in a north–south division in southern Oregon, low genetic distance and a lack of private alleles also support a long bottleneck during the LGM that resulted in a north–south division in southern Oregon (Keir et al., 2011). Cornus nuttallii has been significantly affected by the introduced Discula destructiva (dogwood anthracnose), particularly in moist sites, which occur in the northern range of the species and result in more mature trees at lower elevations and latitudes and indicates that factors other than climate change are more important in affecting distribution (Monleon and Lintz, 2015). The scenario for C. nuttallii is consistent with the prediction that taxa isolated by the environment will result in more locally adapted genotypes that may be able to migrate long distances but ultimately are poorly adapted to different environments. In contrast, species that are genetically isolated by distance may be able to respond more rapidly than those that require specific habitat types, but gene flow may be limited by dispersal ability (Sexton et al., 2014).

Alpine habitats

Alpine habitats in California are generally considered to be those at 3500 m or higher, with a growing season that is limited to the period between the melt of winter snowpack and first snowfall, depending on the latitude, and with intense summer drought (Sharsmith, 1940). A virtually contiguous alpine habitat exists from Tuolumne to Inyo counties in the Sierra Nevada, in addition to habitats on Mount Shasta in northern California, on peaks in southwestern California (e.g. San Gorgonio Mountain) and in the Great Basin on Mount Patterson and the White/Inyo Mountains (Rundel, 2011). Additional alpine habitats at slightly lower elevations are found on Mount Eddy and Thompson Peak in the Klamath–Siskiyou Mountains, and in the Panamint Mountains east of the White/Inyo Mountains in southern California. Species that occur in these habitats can be divided into those that are endemic to subalpine and alpine habitats and those that have broader distributions. California has 31 endemic alpine species (Rundel, 2011).

Significant changes are projected to occur in the subalpine and alpine flora and fauna of California as climate change progresses (Hayhoe et al., 2004; Moritz et al., 2008). Of the plants found in the high Sierra Nevada and Transverse Ranges, those with the smallest ranges are projected to be the most heavily impacted by climate change (Loarie et al., 2008). Alpine organisms often experience population contraction during glacial periods and dramatic population expansion during interglacial periods (DeChaine and Martin, 2004). Depending on dispersal ability, generation time and habitat specificity, genetic diversity is variable but often alpine organisms leave strong genetic signatures. Neoendemics in general are expected to exhibit a variety of patterns due to isolation, colonization of extreme environments, and strong selection that can result in rapid speciation (Schoville et al., 2010).

Colonization of the high Sierra Nevada is suggested to have occurred from north to south, based on a greater number of endemic species in the southern Sierra Nevada (Raven and Axelrod, 1978). North-to-south colonization of high montane taxa more generally may explain why Mount Lassen shares greater floristic affinity with Mount Shasta than the closer northern Sierra Nevada (Gillett et al., 1995). Similarly, the Klamath Mountains are the southern limit of a number of boreal alpine species (Howell, 1944). Generally, these patterns support a northern origin for the alpine floras of the Cascade Range and Sierra Nevada. Species with their southern limit of distribution in the northern or central high Sierra Nevada are Carex whitneyi, Podistera nevadensis, Claytonia megarhiza, Thalictrum alpinum, Galium grayanum and Salix nivalis (Rundel, 2011). As a result of this colonization pattern and generally decreasing levels of genetic diversity from north to south, southern endemic species and populations will be most vulnerable. Nine of the Sierra Nevada alpine species are restricted to elevations of 3500 m and 36 species are endemic to the Sierra Nevada at ≥2700 m (Rundel, 2011).

Geographically, the alpine flora has diverse origins. The largest group of species (34·3 %) are those that are widespread in the mountainous western USA, 20·5 % occur in the Intermountain Great Basin, 15·8 % in the Sierra Nevada and into the Cascade Range, and 13·6 % are widely distributed in boreal or arctic–alpine regions (Rundel, 2011). A majority of both the Sweetwater (94 %) and White Mountain (90 %) alpine floras are shared with the Sierra Nevada (Scott, 1995). These floras should be considered at great risk. However, many species present there may survive in the northern Sierra Nevada. Models to predict how populations will respond to climate change will need to include changes in population size, demographics and connectivity during changing warm and cold phases, with the intensity of gene flow dependent on individual species characteristics.

Models examining the effects of climate change on European alpine species for temperature increases of 2·5–4  °C found significant changes in population-level genetic structure that became more dramatic as temperatures increased (Jay et al., 2012). Although these models are based on neutral genetic variation, adaptive alleles are expected to follow the same patterns (Thibert-Plante and Hendry, 2010; Bierne et al., 2011). Additionally, invasions by lower-elevation or more southern populations has the potential to increase genetic variation in these populations, provided there are opportunities for gene flow (Verhoeven et al., 2011).

Although there are molecularly-based population genetic studies of alpine plant species in the European Alps and the Spanish Sierra Nevada, beyond Abronia alpina there are no range-wide population genetic studies of California alpine taxa. Abronia alpina is a rare alpine species endemic to the southern Sierra Nevada that is obligately outcrossing, making it dependent on pollinators to maintain genetic diversity. Although genetic diversity is currently high, a loss or upward elevational migration of pollinators will minimally result in a loss of genetic diversity and possibly extinction (Jabis et al., 2011). It is critical to the persistence of California’s alpine plant taxa that assessments of population genetic diversity and connectivity and projections to habitat change be conducted.

Central Valley wetlands

The vernal pools of California have suffered devastating losses in the past 150 years, mostly due to habitat conversion. Formerly located throughout the Central Valley, less than 10 % of vernal pools are estimated to remain (Holland, 1976). Almost 200 plant species are either endemic or associated with vernal pools, 55 % of which are endemic to California (Holland, 1976). The presence of vernal pools in flat, easily accessible land has made them prone to development pressure and alterations to surrounding areas can negatively affect their hydrological properties (Cheatham, 1976). In addition to their direct loss, reduced connectivity among vernal pools has resulted in the extirpation and loss of gene flow in a number of taxa, both plant and animal.

Conservation of vernal pool taxa requires the assessment of gene flow within and among populations to determine their degree of genetic isolation (Sloop et al., 2011) and their ability to persist under altered hydrological and climate change. An assessment of genetic variation in Neostapfia colusana, a federal threatened species endemic to vernal pools, found high levels of heterozygosity (HO = 0·68, HE = 0·71), but an FST value of 0·268 (P < 0·0001) between northern and southern populations indicates that gene flow is very limited among the regions and likely due to the loss of geographically intermediate populations (Sloop et al., 2011). Some populations had lower than expected heterozygosity that was likely due to geographic isolation and increased inbreeding. Similarly, Tuctoria greenei, a federal endangered vernal pool endemic species evaluated for population-level genetic variation, found high heterozygosity (HO = 0·77, HE = 0·79), also with regional subdivision between northern and southern populations, although gene flow was not as restricted as in N. colusana (FST = 0·11, P = 0·0001) (Gordon et al., 2012).

The Limnanthaceae is a small clade endemic to western North America, with the exception of Floerkea proserpinacoides, which occurs in pockets throughout North America (Kishore et al., 2004). The genus Limnanthes consists of 18 taxa that occur primarily in wetland habitats in the Central Valley, 12 of which are listed as threatened, rare or endangered. Limnanthes floccosa ssp. californica is a federal and state-listed endangered species endemic to vernal pools along the eastern edge of the Sacramento Valley and currently known from only 14 geographically isolated populations. The species has fairly low genetic diversity, with HO = 0·10 and HE = 0·19, and high among-population genetic structure, as measured by pairwise FST values of 0·12 – 0·79 (Sloop et al., 2011). Of particular concern for this species is its high rate of inbreeding and barriers to gene flow. Pollinator declines and reduced aquatic dispersal due to habitat loss will likely further erode genetic diversity in this species.

Quercus lobata, endemic to California and the western USA, occurs in remnant populations in the Central Valley and the surrounding valleys of the Coast Ranges, Sierra Nevada and Transverse Ranges of California. The Mediterranean climate that became established by the mid-Pliocene in west central California and continued through the Pleistocene glaciations affected distributions of Quercus and associated species (Griffin, 1988). Genetic data support strong geographic structuring of populations within the species, determined to be about 200 km (Grivet et al., 2006); however, the genetic connectivity of current populations supports long-distance colonization (Sork et al., 2010). Combined with palynological data, these results suggest the Neogene was an important period for Quercus woodlands as they expanded and contracted with glaciations. It will be important to maintain the remaining connectivity in Q. lobata to protect its ability to adapt, as southern populations are more likely to be lost.

Taxa such as Carex, a worldwide genus with its centres of diversity in Northern Hemisphere temperate regions, are likely to fare better in the face of warming climates; however, more isolated populations or species are likely to be in peril. Only 30 known plant species have a bipolar distribution, six of which are found in the genus Carex (Moore and Chater, 1971). Diversification via long-distance dispersal in the late Miocene through the Pliocene is hypothesized for bipolar Carex species (Escudero et al., 2010). Molecular studies using chloroplast and nuclear DNA of the subgenera Vignea and Carex, both of which occur at high latitudes in the Northern and Southern Hemispheres, suggest that periodic glaciations events may have played a role in speciation (Martín-Bravo and Escudero, 2012). This pattern is consistent with Carex as whole (2000 species) as it has an inverse latitudinal gradient of species richness (Hillebrand, 2004), which will likely become more pronounced as climates warm.

Coast, Transverse and Peninsular Ranges

Taxa with populations that have close proximity to the confluence of two or more divergent floras often have high levels of genetic diversity, reflecting a complex evolutionary history. It will be essential to maintain landscape connectivity in these regions to maintain migrational corridors and gene flow in the context of anthropogenic change. The Tehachapi Mountains are at the confluence of major regional biotas and historically have provided migrational corridors between the central coast, Transverse Ranges, the Sierra Nevada and the Mojave Desert. The Tehachapi Mountains and surrounding areas and have long been recognized as containing high biological diversity and evolutionary divergence (Stebbins and Major, 1965; Patton and Smith, 1990). Examples of non-plant taxa with clade breaks across the Tehachapi Mountains are numerous (Schierenbeck, 2014); however, surveys of the population genetics or high-resolution phylogenetic analyses of plants in this area are few. High genetic diversity with evidence of northern expansion into the Sierra Nevada or western expansion into the Coast Range is found in the Tehachapi Mountains for some taxa. The area has long been recognized as an evolutionary hotspot, with the proposal of the area for genetic connectivity between the coast and Transverse Range (Peabody and Savage, 1958) and later as an area for high hybridization between xeric and mesic lineages (Remington, 1968). High levels of relictual plant species, or palaeoendemics with very restricted distributions, provide evidence of range expansion of cooler-adapted species into southern latitudes, with concentrations in the San Jacinto and Santa Rosa Mountains (Stebbins and Major,1965; Raven and Axelrod, 1978; Baldwin, 2014). The San Gorgonio Pass in the San Bernardino Mountains has a very strong cline of habitat types and is the location of many divergent clades (Patton et al., 2008). Conservation efforts have improved in recent years but much of the area remains under threat. The Tejon Ranch, recently placed under conservation easement, is 109 265 ha at the juncture of the Tehachapi Mountains, the San Joaquin Valley, the Transverse Ranges and the Mojave Desert and is home to 20 state-listed and federally listed species and 60 other rare taxa. This area will continue to serve as an important reserve for genetic diversity and cooler microclimates will provide refuge for many species unable to migrate with climate change.

The Transverse Ranges and Los Angeles Basin Pliocene embayment break is an important phylogeographic break and often separates central and southern California for a number of clades (Calsbeek et al., 2003; Chatzimanolis and Caterino, 2007; Polihronakis and Caterino, 2010). Species distributions were affected by expansion and contraction of mesic habitats during the Pleistocene. Expansion in southern clades happened from southern refugia and relationships among genetic divergence and geographic distance indicate a gradual movement north. Unfortunately, the Transverse Ranges have only about 30 % of their original habitat intact, including the Ventana Wilderness Area in the Santa Lucia Range and the Ventura region. There are very high levels of fragmentation and isolation of contiguous habitat. Although much of the land falls within the Los Padres National Forest, the area is heavily impacted by logging, air pollution, disturbance to aquatic habitat, overgrazing, heavy recreational use and fire suppression (David Olson, TNC, pers. comm.).

Genetic distances among taxa in the Transverse and Peninsular Ranges are generally low and reflect late Pleistocene divergences (Calsbeek et al., 2003). Populations of species with northern populations may persist but not those with limited southern California distributions, limited dispersal or strong sensitivity to human disturbance. There are few plant studies of this region, but amphibian studies indicate that it will be important to identify those areas that will be suitable for reproduction and ongoing gene flow in this evolutionary hotspot (Devitt et al., 2013).

Taxa with limited distributions in heavily populated areas are particularly at risk. For example, Dirca occidentalis, a California endemic found only in the San Francisco Bay Area, is known from four genetically distinct populations hypothesized to have resulted from climatic events within the last 20 000 years (Graves and Schrader, 2008). Species such as D. occidentalis, dependent on mesic habitat in areas with dense human populations, are in jeopardy.

CONCLUSIONS

The prediction of how species will respond to climate change will require a synthesis drawing from population genetics, geography, palaeontology and multiple ecological parameters (Table 1). The important integration of these fields will enable us to predict and prioritize conservation areas during a time of rapid climate change, human disturbance and invasive species. Phylogeographic methods, ecological niche modelling, fossils and palaeoclimatic reconstructions all contribute to the inference of refugia during the LGM, which are now reflected as areas with high genetic diversity or as contact zones. These refugia will continue to be significant sources of genetic variation, but because climate is now warming, not cooling, populations at lower elevations and latitudes will also provide important reservoirs of genetic diversity. Although earlier interpretations of the Pleistocene fossil record indicated that responses to climate change were idiosyncratic (Graham et al., 1996), recent evidence shows a number of patterns among groups with similar life history characteristics. Vagility, isolation of the habitat, reproductive output and generation time all affect the ability of species to respond to climate change.

Table 1

Factors influencing the genetic structure of plant populations and individual species’ responses to climate change. Note that ability to respond to climate change can be dependent on the interaction of these effects. Positive to negative should be considered a continuum

FactorPositiveNegative
Genetic factors
 Effective population sizeLarge, but influenced by range sizeSmall
 Locally adapted genotypesDepends on limits of adaptation and habitat availabilityGenerally, limits ability to migrate
FSTVariable. Dependent on number of populations, higher FST values can indicate lack of gene flow, but also higher total within-species diversity
Ecological factors
 PollinationWind/generalistsSpecialists
 DispersalWidespreadLimited
Demographic factors
 Reproductive rateHigh, continuousLow, occasional
 Demographic reaction normsWideNarrow
 DemographyMany age-groupsLack of consistent reproduction. Dependent on human-induced disturbance regime
Geographic factors
 Spatial isolationVariable, dependent on extent of isolation and life history factors associated with gene flowRemoval of geographically intermediate populations
 RangeLargeSmall
 EndemismVariable, dependent on geographic extent endemism and proximity to other populations
 RarityDependent on extent and reason for rarityDisturbance, hybridization with close relatives
 Habitat availabilityWidespreadLimited, or isolated
 Migration corridorsWidely available, intactNot available, heavily disturbed
FactorPositiveNegative
Genetic factors
 Effective population sizeLarge, but influenced by range sizeSmall
 Locally adapted genotypesDepends on limits of adaptation and habitat availabilityGenerally, limits ability to migrate
FSTVariable. Dependent on number of populations, higher FST values can indicate lack of gene flow, but also higher total within-species diversity
Ecological factors
 PollinationWind/generalistsSpecialists
 DispersalWidespreadLimited
Demographic factors
 Reproductive rateHigh, continuousLow, occasional
 Demographic reaction normsWideNarrow
 DemographyMany age-groupsLack of consistent reproduction. Dependent on human-induced disturbance regime
Geographic factors
 Spatial isolationVariable, dependent on extent of isolation and life history factors associated with gene flowRemoval of geographically intermediate populations
 RangeLargeSmall
 EndemismVariable, dependent on geographic extent endemism and proximity to other populations
 RarityDependent on extent and reason for rarityDisturbance, hybridization with close relatives
 Habitat availabilityWidespreadLimited, or isolated
 Migration corridorsWidely available, intactNot available, heavily disturbed
Table 1

Factors influencing the genetic structure of plant populations and individual species’ responses to climate change. Note that ability to respond to climate change can be dependent on the interaction of these effects. Positive to negative should be considered a continuum

FactorPositiveNegative
Genetic factors
 Effective population sizeLarge, but influenced by range sizeSmall
 Locally adapted genotypesDepends on limits of adaptation and habitat availabilityGenerally, limits ability to migrate
FSTVariable. Dependent on number of populations, higher FST values can indicate lack of gene flow, but also higher total within-species diversity
Ecological factors
 PollinationWind/generalistsSpecialists
 DispersalWidespreadLimited
Demographic factors
 Reproductive rateHigh, continuousLow, occasional
 Demographic reaction normsWideNarrow
 DemographyMany age-groupsLack of consistent reproduction. Dependent on human-induced disturbance regime
Geographic factors
 Spatial isolationVariable, dependent on extent of isolation and life history factors associated with gene flowRemoval of geographically intermediate populations
 RangeLargeSmall
 EndemismVariable, dependent on geographic extent endemism and proximity to other populations
 RarityDependent on extent and reason for rarityDisturbance, hybridization with close relatives
 Habitat availabilityWidespreadLimited, or isolated
 Migration corridorsWidely available, intactNot available, heavily disturbed
FactorPositiveNegative
Genetic factors
 Effective population sizeLarge, but influenced by range sizeSmall
 Locally adapted genotypesDepends on limits of adaptation and habitat availabilityGenerally, limits ability to migrate
FSTVariable. Dependent on number of populations, higher FST values can indicate lack of gene flow, but also higher total within-species diversity
Ecological factors
 PollinationWind/generalistsSpecialists
 DispersalWidespreadLimited
Demographic factors
 Reproductive rateHigh, continuousLow, occasional
 Demographic reaction normsWideNarrow
 DemographyMany age-groupsLack of consistent reproduction. Dependent on human-induced disturbance regime
Geographic factors
 Spatial isolationVariable, dependent on extent of isolation and life history factors associated with gene flowRemoval of geographically intermediate populations
 RangeLargeSmall
 EndemismVariable, dependent on geographic extent endemism and proximity to other populations
 RarityDependent on extent and reason for rarityDisturbance, hybridization with close relatives
 Habitat availabilityWidespreadLimited, or isolated
 Migration corridorsWidely available, intactNot available, heavily disturbed

The risk of extinction can be affected by genetic diversity, range size, demography, body size and life history characteristics, especially fecundity. The most vulnerable species are endemic species because of their habitat specificity and rarity, and the most vulnerable ecosystems are those with high levels of spatial isolation, e.g. alpine habitat and serpentine islands. Spatial isolation can mean different evolutionary dynamics in different populations of the same species in terms of gene flow, genetic drift, and strong selection from climate change and invasive species. Ecosystem function is dependent on biological diversity as many ecosystem processes are sensitive to declines in biodiversity. Local declines in biodiversity can be more devastating than global declines as the benefits that organisms provide for local ecosystems can be lost before a global extinction.

POTENTIAL SOLUTIONS

Issues to consider in mitigating climate change include local, regional and taxonomic rarity, ecological function and current threats to biodiversity. The specific nature of species–environment interactions sometimes appears idiosyncratic because of our lack of specific understanding of an ecosystem or the ecological requirements of particular species. Experimental examination of the biotic factors that control genetic diversity and their interaction with climate change should focus on areas that contain the greatest number of endemic species or large amounts of biodiversity, e.g. the central Coast Ranges, the Klamath–Siskiyou region, the eastern Transverse Ranges, the Tehachapi Mountains, the Santa Lucia Mountains and the eastern Mojave Desert. Areas with high levels of endemism that include species with the greatest range restrictions are on the central coast, in the Sierra Nevada and in the San Bernardino Mountains (Kraft et al., 2010). The Mojave Desert and Great Basin have the highest level of neoendemics. High levels of endemism, however, do not necessarily coincide with high levels of genetic diversity, which will likely affect the ability of a species to persist in a rapidly changing climate. Due to its extant biodiversity, latitudinal position, topographical complex and minimal development, and the presence of many taxa with splits between northern and southern clades, the Klamath–Siskiyou region will remain an important refugial region.

Far northern California represents a confluence of the volcanic Cascades, the geologically ancient and complex Klamath–Siskiyou region, the northern Sierra Nevada and the Warner Mountains, separated from the Sierra Nevada by the Modoc Plateau. The species that have hybridization zones or reach the northern or southern limit of their distribution provide varying degrees of divergence. Population genetic studies of plant taxa in the Warner Mountains are virtually nil outside of conifers, but there have been a number of vertebrate and invertebrate studies that establish this region as a contact zone. Because divergence patterns of plant and animal follow similar patterns (Calsbeek et al., 2003), this region should reveal interesting population structures within and among plant taxa. Although not currently considered a priority area, this region has been poorly explored and merits consideration. For example, this region has been documented for the evolutionary divergence of northern and southern subspecies of the bird species Picoides albolarvatus (Alexander and Burns, 2006) and Pacific Northwest and California lineages of Strix occidentalis (Barrowclough et al., 2005). Habitat and climatic differences correspond to a split between the subspecies bird taxa Aphelocoma californica woodhouseii and Aphelocoma californicacalifornica over a distance of about 500 km, and similar scenarios exist between bird taxa Baeolophus inornatus and Baeolophusridgwayi and mammal taxa Pica nuttalli and Picahudsonia (Delaney et al., 2008).

RECOMMENDATIONS

The general long-term goal of conservation genetics is to conserve historical genetic variation that may be critical to the long-term evolutionary survival of a species. Usually this is accomplished through the conservation of genetically distinct populations. In addition, quantitative trait locus analysis can identify traits of functional diversity or those that are important in ecological interactions. In addition, molecular clock data, which inform divergence times combined with the knowledge of vicariant events, can be critical to the survival of a species (Shaffer et al., 2004). For example, highly restricted palaeoendemic taxa with little genetic variation, such as Hesperocyparis forbesii, will be difficult for large-scale conservation outside of its current range. Conservation strategies for recently emerged taxa will vary depending on the time since emergence, genetic variation, existence of sister taxa, and range. As populations become demographically at risk, some reintroductions may be necessary and an understanding of the genetic basis of adaptive variation will need to be accomplished with genetically and ecologically similar populations.

Genetically divergent populations occurring in fragmented habitats need designation as evolutionarily significant management units. Recommendations to ensure the survival of species include the use of genetically appropriate choices for plant materials in ecosystem restoration to maintain biological diversity, and include provision of the genetic diversity necessary for evolution to continue to occur. Appropriate genetic linkages among fragmented populations must be established, evolutionary interactions conserved, surrounding biological diversity maintained, and changes to ecosystem function minimized (Rogers and Montalvo, 2004).

Genotype specificity can be very important in site restoration. Species without quantitative or molecular differentiation suggest that propagule transfer should have a lower risk of maladaptation. Every effort should be made to assess projected climatic shifts and whether the reproductive and phenotypic traits are suitable to future scenarios. The use of native genotypes in restoration is important in the avoidance of outbreeding depression, which can result in the decline of a native population (Erickson, 2008; Aitken and Whitlock, 2013). Local ecotypes are adapted to site conditions, and even though the use of natives is often more expensive there is greater restoration success with their use (Lesica and Allendorf, 1999). However, there have been recent recommendations that local genotypes could be supplemented with non-local genotypes as bet-hedging against unforeseen environmental change (Doak and Morris, 2010).

Phylogeographic studies can be very useful in defining evolutionarily significant units for conservation purposes (Roderick, 1996). Concordance among divergent major clades provides the basis needed for major conservation efforts in particular regions (Avise et al., 1987). Many regions of high diversity are in need of further support in terms of conservation efforts that will strengthen legal protection, provide and strengthen migratory corridors, and increase the size of parks and preserves to protect biological diversity. Population genetic data can be informative for management and should be utilized to determine whether genetic diversity is being lost (Funk et al., 2010).

LITERATURE CITED

Abbott
RJ.
1992
.
Plant invasions, interspecific hybridization and the evolution of new plant taxa
.
Trends in Ecology and Evolution
7
:
401
405
.

Aitken
SN
Whitlock
MC.
2013
.
Assisted gene flow to facilitate local adaptation to climate change
.
Annual Review of Ecology, Evolution, and Systematics
44
:
367
388
.

Alexander
MP
Burns
KJ.
2006
.
Intraspecific phylogeography and adaptive divergence in the white-headed woodpecker
.
Condor
108
:
489
508
.

Anacker
BL
Gogol-Prokurat
M
Leidhom
K
Schoenig
S.
2013
.
Climate change vulnerability assessment of rare plants in California
.
Madroño
60
:
193
210
.

Anttila
CK
Daehler
CC
Rank
NE
Strong
DR.
1998
.
Greater male fitness of a rare invader (Spartina alterniflora, Poaceae) threatens a common native (Spartina foliosa) with hybridization
.
American Journal of Botany
85
:
1597
1601
.

Avise
JC
Arnold
J
Ball
RM
, et al.
1987
.
Intraspecific phylogeography: the mitochondrial DNA bridge between population genetics and systematics
.
Annual Review of Ecology and Systematics
18
:
489
522
.

Ayres
DR
Fleishman
E
Launer
A
Lee
AK
Zippin
D.
2015
.
Genetic structure and demography of Chloropyron palmatum, an endangered annual plant
.
Madroño
62
:
139
149
.

Baldwin
BG.
2014
.
Origins of plant diversity in the California Floristic Province
.
Annual Review of Ecology, Evolution, and Systematics
45
:
347
369
.

Baldwin
BG
Sanderson
MJ.
1998
.
Age and rate of diversification of the Hawaiian silversword alliance (Compositae)
.
Proceedings of the National Academy of Sciences of the USA
95
:
9402
9406
.

Baldwin
BG
Kalisz
S
Armbruster
WS.
2011
.
The need for a detailed phylogenetic understanding of Collinsia Nutt. (Plantaginaceae) has grown with increasing attention to the genus as a subject of diverse studies in mating-system
.
American Journal of Botany
98
:
731
753
.

Baldwin
BG
Goldman
DH
Keil
DJ
Patterson
RW
Rosatti
TJ
Wilken
DH.
eds.
2012
.
The Jepson manual: vascular plants of California
, 2nd edn.
Berkeley
:
University of California Press
.

Barrowclough
GF
Groth
JG
Mertz
LA
Gutiérrez
RJ.
2005
.
Genetic structure, introgression, and a narrow hybrid zone between northern and southern spotted owls (Strix occidentalis)
.
Molecular Ecology
14
:
1109
1120
.

Bellard
C
Bertelsmeier
C
Leadley
P
Thuiller
W
Courchamp
F.
2012
.
Impacts of climate change on the future of biodiversity
.
Ecology Letters
15
:
365
377
.

Bermingham
E
Moritz
C.
1998
.
Comparative phylogeography: concepts and applications
.
Molecular Ecology
7
:
367
369
.

Betancourt
JL
Van Devender
TR
Martin
PS.
1990
. Synthesis and prospectus. In:
Betancourt
JL
Van Devender
TR
Martin
PS
, eds.
Packrat middens. The last 40,000 years of biotic change
.
Tucson
:
University of Arizona Press
,
435
447
.

Bierne
N
Welch
J
Loire
E
Bonhomme
F
David
P.
2011
.
The coupling hypothesis: why genome scans may fail to map local adaptation genes
.
Molecular Ecology
20
:
2044
2072
.

Bonin
A
Ehrich
D
Manel
S.
2007
.
Statistical analysis of amplified fragment length polymorphism data: a toolbox for molecular ecologists and evolutionist
.
Molecular Ecology
16
:
3737
3758
.

Boulangeat
I
Georges
D
Dentant
C
Bonet
R
, et al.
2014
.
Anticipating the spatio-temporal response of plant diversity and vegetation structure to climate and land use change in a protected area
.
Ecography
37
:
1
10
.

Boykin
LM
Vasey
MC
Parker
VT
Patterson
R.
2005
.
Two lineages of Arctostaphylos (Ericaceae) identified using the internal transcribed spacer (ITS) region of the nuclear genome
.
Madroño
52
:
139
147
.

Brandvain
Y
Kenney
AM
Flagel
L
Coop
G
Sweigart
AL.
2014
.
Speciation and introgression between Mimulus nasutus and Mimulus guttatus
.
PLoS Genetics
10
:
e1004410
. doi:10.1371/journal.pgen.1004410.

Brubaker
LB
Anderson
PM
Edwards
ME
Lozhkin
AV.
2005
.
Beringia as a glacial refugium for boreal trees and shrubs: new perspectives from mapped pollen data
.
Journal of Biogeography
32
:
833
848
.

Burge
DO
Erwin
DM
Islam
MB
, et al.
2011
.
Diversification of Ceanothus (Rhamnaceae) in the California Floristic Province
.
International Journal of Plant Sciences
172
:
1137
1164
.

Burge
DO
Thorne
JH
Harrison
SP
, et al.
2016
.
Plant diversity and endemism in the California Floristic Province
.
Madroño
63
:
3
206
.

Cadotte
MW
Davies
J.
2010
.
Rarest of the rare: advances in combining evolutionary distinctiveness and scarcity to inform conservation at biogeographical scales
.
Diversity and Distributions
16
:
376
385
.

Calsbeek
R
Thompson
JN
Richardson
J.
2003
.
Patterns of molecular evolution and diversification in a biodiversity hotspot: the California Floristic Province
.
Molecular Ecology
12
:
1021
1029
.

Cayan
DR
Maurer
EP
Dettinger
MD
Tyree
M
Hayhoe
K.
2008
.
Climate change scenarios for the California region
.
Climatic Change
87
:
S21
S42
.

Chatzimanolis
S
Caterino
MS.
2007
.
Toward a better understanding of the ′Transverse Range Break′: lineage diversification in southern California
.
Evolution
61
:
2127
2141
.

Cheatham
NH.
1976
. Conservation of vernal pools. In:
Jain
S.
ed.
Vernal pools: their ecology and conservation. Institute of Ecology Publication No. 9
.
Davis
:
University of California
,
86
89
.

Chen
I-C
Hill
JK
Ohlemüller
R
Roy
DB
Thomas
CD.
2011
.
Rapid range shifts of species associated with high levels of climate warming
.
Science
333
:
1024
1026
.

Chevin
L-M
Collins
S
Lefevre
F.
2013
.
Phenotypic plasticity and evolutionary demographic responses to climate change: taking theory out to the field
.
Functional Ecology
27
:
966
979
.

Davey
JW
Blaxter
ML.
2010
.
RADseq: next generation population genetics
.
Briefings in Functional Genomics
9
:
416
423
.

Davies
ID
Cary
GJ
Landguth
EL
Lindenmayer
DB
Banks
SC.
2016
.
Implications of recurrent disturbance for genetic diversity
.
Ecology and Evolution
6
:
1181
1196
.

Davis
MB
Shaw
RG.
2001
.
Range shifts and adaptive responses to quaternary climate change
.
Science
292
:
673
679
.

DeChaine
EG
Martin
AP.
2004
.
Historic cycles of fragmentation and expansion in Parnassius smintheus (Papilionidae) inferred using mitochondrial DNA
.
Evolution
58
:
113
127
.

Delaney
KS
Zafar
S
Wayne
RK.
2008
.
Genetic divergence and differentiation within the western scrub-jay (Aphelocoma californica)
Auk
125
:
839
849
.

DellaSala
DA
Ried
SB
Frest
TJ
Strittholt
JR
Olson
DM.
1999
.
A global perspective on the biodiversity of the Klamath-Siskiyou Ecoregion
.
Natural Areas Journal
19
:
300
319
.

Devitt
TJ
Devitt
C
Hollingsworth
BD
McGuire
JA
Moritz
C.
2013
.
Montane refugia predict population genetic structure in the large-blotched Ensatina salamander
.
Molecular Ecology
22
:
1650
1665
.

Dirnböck
T
Essl
F
Rabitsch
W.
2010
.
Disproportional risk for habitat loss of high-altitude endemic species under climate change
.
Global Change Biology
17
:
990
996
.

Doak
DF
Morris
WF.
2010
.
Demographic compensation and tipping points in climate-induced range shifts
.
Nature
467
:
959
962
.

Dole
KP
Loik
ME
Sloan
LC.
2003
.
The relative importance of climate change and the physiological effects of CO2 on freezing tolerance for the future distribution of Yucca brevifolia
.
Global and Planetary Change
36
:
137
146
.

Eckert
AJ.
2011
.
Seeing the forest for the trees: statistical phylogeography in a changing world
.
New Phytologist
189
:
894
897
.

Eckert
AJ
Carstens
BC.
2008
.
Does gene flow destroy phylogenetic signal? The performance of three methods for estimating species phylogenies in the presence of gene flow
.
Molecular Phylogenetics and Evolution
49
:
832
842
.

Eckert
AJ
Hall
BD.
2006
.
Phylogeny, historical biogeography, and patterns of diversification for Pinus (Pinaceae): phylogenetic tests of fossil-based hypotheses
.
Molecular Phylogenetics and Evolution
40
:
166
182
.

Elam
DR
Lewiner
K
Keeler-Wolf
T
Flint
SA.
1998
. Vernal pool assessment, preliminary report. Prepared by the California Natural Diversity Database. State of California, Sacramento, CA, USA: Resources Agency, Department of Fish and Game.

Ellstrand
NC
Elam
DR.
1993
.
Population genetic consequences of small population size: implications for plant conservation
.
Annual Review of Ecology and Systematics
24
:
217
242
.

Erickson
VJ.
2008
.
Developing native plant germplasm for national forests and grasslands in the Pacific Northwest
.
Native Plant Journal
9
:
3255
3266
.

Escobar García
P
Winkler
M
Flatscher
R
, et al. .
2012
.
Extensive range persistence in peripheral and interior refugia characterizes Pleistocene range dynamics in a widespread alpine plant species (Senecio carniolicus, Asteraceae)
.
Molecular Ecology
21
:
1255
1270
.

Escudero
M
Hipp
AL
Luceño
M.
2010
.
Karyotype stability and predictors of chromosome number variation in sedges: a study in Carex section Spriostachyae (Cyperaceae)
.
Molecular Phylogenetics and Evolution
57
:
353
363
.

Farjon
A.
2013
. Cupressus macrocarpa. In: The IUCN Red List of Threatened Species 2013: e.T30375A2793139. http://dx.doi.org/10.2305/IUCN.UK.2013-1.RLTS.T30375A2793139.en.

Fehlberg
SD
Ranker
TA.
2009
.
Evolutionary history and phylogeography of Encelia farinosa (Asteraceae) from the Sonoran, Mojave, and Peninsular Deserts
.
Molecular Phylogenetics and Evolution
50
:
326
335
.

Forrest
JRK.
2015
.
Plant-pollinator interactions and phenological change: what can we learn about climate impacts from experiments and observations?
Oikos
124
:
4
13
.

Frankham
R.
2010
.
Challenges and opportunities of genetic approaches to biological conservation
.
Biological Conservation
143
:
1919
1927
.

Funk
WC
Forsman
ED
Johnson
M
Mullins
TD
Haig
SM.
2010
.
Evidence for recent population bottlenecks in northern spotted owls (Strix occidentalis caurina)
.
Conservation Genetics
11
:
1013
1021
.

Gernandt
DS
Geada López
G
Ortiz García
S
Liston
A.
2005
.
Phylogeny and classification of Pinus
.
Taxon
54
:
29
42
.

Gillett
GW
Howell
JT
Leschke
H
, revised by
Oswald
VH
Showers
DW
Showers
M.
1995
.
A flora of Lassen Volcanic National Park, California
.
Sacramento, CA
:
California Native Plant Society
.

Gilman
RT
Fabina
NS
Abbott
KC
Rafferty
NE.
2012
.
Evolution of plant-pollinator mutualisms in response to climate change
.
Evolutionary Applications
5
:
2
16
.

Gordon
SP
Sloop
CM
Davis
HG
Cushman
JH.
2012
.
Population genetic diversity and structure of two rare vernal pool grasses in central California
.
Conservation Genetics
13
:
117
130
.

Graham
RW
Lundelius
EL
Jr
Graham
MA.
1996
.
Spatial response of mammals to late-Quaternary environmental fluctuations
.
Science
272
:
1601
1606
.

Graves
WR
Schrader
JA.
2008
.
At the interface of phylogenetics and population genetics, the phylogeography of Dirca occidentalis (Thymelaeaceae)
.
American Journal of Botany
95
:
1454
1465
.

Griffin
JR.
1988
. Oak woodland. In:
Barbour
MG
Major
J.
eds.
Terrestrial vegetation of California. Special Publication 9
, 2nd edn.
Sacramento
:
California Native Plant Society
,
383
415
.

Grivet
DM-F
Petit
RJ
Sork
VL.
2006
.
Contrasting patterns of historical colonization in white oaks (Quercus spp.) in California and Europe
.
Molecular Ecology
15
:
4085
4093
.

Gugger
PF
Sugita
S
Cavender-Bares
J.
2010
.
Phylogeography of Douglas-fir based on mitochondrial and chloroplast DNA sequences: testing hypotheses from the fossil record
.
Molecular Ecology
19
:
1877
1897
.

Gugger
PF
González-Rodríguez
A
Rodríguez-Correa
H
Sugita
S
Cavender-Bares
J.
2011
.
Southward Pleistocene migration of Douglas-fir into Mexico: phylogeography, ecological niche modeling, and conservation of ‘rear edge’ populations
.
New Phytologist
189
:
1185
1199
.

Guppy
CS
Shepard
JH.
2001
.
Butterflies of British Columbia: including Western Alberta, Southern Yukon, the Alaska Panhandle, Washington, Northern Oregon, Northern Idaho, and northwestern Montana
.
Vancouver, BC
:
UBC Press
.

Halbur
MM
Sloop
CM
Zanis
MJ
Emery
NC.
2014
.
The population biology of mitigation: impacts of habitat creation on an endangered plant species
.
Conservation Genetics
15
:
679
695
.

Hampe
A
Petit
RJ.
2005
.
Conserving biodiversity under climate change: the rear edge matters
.
Ecology Letters
8
:
461
467
.

Hamrick
JL
Linhart
YB
Mitton
JB.
1979
.
Relationships between life history characteristics and electrophoretically detectable genetic variation in plants
.
Annual Review of Ecology and Systematics
10
:
173
200
.

Hamrick
JL
Godt
MJW
Sherman-Broyles
SL.
1992
.
Factors influencing levels of genetic diversity in woody plant species
.
New Forests
6
:
95
124
.

Hamrick
JL
Godt
MJW.
1996
. Conservation genetics of endemic plant species. In:
Avise
JC
Hamrick
JL.
eds.
Conservation genetics: case histories from nature
.
London
:
Chapman & Hall
,
281
304
.

Hanski
I
Saastamoinen
M
Ovaskainen
O.
2006
.
Dispersal-related life-history trade-offs in a butterfly metapopulation
.
Journal of Animal Ecology
75
:
91
100
.

Hayhoe
K
Cayan
D
Field
CB
, et al.
2004
.
Emissions pathways, climate change, and impacts on California
.
Proceedings of the National Academy of Sciences of the USA
101
:
12422
12427
.

Hillebrand
H.
2004
.
On the generality of the latitudinal diversity gradient
.
American Naturalist
163
:
192
211
.

Hoffman
AA
Willi
Y.
2008
.
Detecting genetic responses to environmental change
.
Nature Reviews Genetics
9
:
421
432
.

Holderegger
R
Kamm
U
Gugerli
F.
2006
.
Adaptive vs. neutral genetic diversity: implications for landscape genetics
.
Landscape Ecology
21
:
797
807
.

Holland
R.F.
1976
. The vegetation of vernal pools: a survey. In:
Jain
S.
ed.
Vernal pools: their ecology and conservation. Institute of Ecology Publication No. 9
.
Davis
:
University of California
,
11
15
.

Howell
JT.
1944
.
Certain plants of the Marble Mountains in California with remarks on the boreal flora of the Klamath area
.
Wasmann Collector
6
:
13
20
.

Jabis
MD
Ayers
TJ
Allan
GJ.
2011
.
Pollinator-mediated gene flow fosters genetic variability in a narrow alpine endemic Abronia alpina (Nytaginaceae)
.
American Journal of Botany
98
:
1583
1594
.

Jay
F
Manel
S
Alvarez
N
, et al.
2012
.
Forecasting changes in population genetic structure of alpine plants in response to global warming
.
Molecular Ecology
21
:
2354
2368
.

Johnson
JR
Fitzpatrick
BM
Shaffer
HB.
2010
.
Retention of low-fitness genotypes over six decades of admixture between native and introduced tiger salamanders
.
BMC Evolutionary Biology
10
:
147
.

Johnson
MG
Lang
K
Manos
P
Golet
GH
Schierenbeck
KA.
2016
.
Evidence for genetic pollution of a California native tree, Platanus racemosa, via recent, ongoing introgressive hybridization with an introduced ornamental species
.
Conservation Genetics
17
:
593
. doi:10.1007/s10592-015-0808-z.

Jorgensen
SM
Hamrick
JL.
1997
.
Biogeography and population genetics of whitebark pine, Pinus albicaulis
.
Canadian Journal of Forest Research
27
:
1574
1585
.

Jump
AS
Peñuelas
J.
2005
.
Running to stand still: adaptation and the response of plants to rapid climate change
.
Ecology Letters
8
:
1010
1020
.

Kafton
DL.
1976
. Isozyme variability and reproductive phenology in Monterrey cypress. PhD Thesis, University of California, Berkeley.

Keir
KR
Bemmels
JB
Aitken
SN.
2011
.
Low genetic diversity, moderate local adaptation, and phylogeographic insights in Cornus nuttallii (Cornaceae)
.
American Journal of Botany
98
:
1327
1336
.

Kishore
VK
Velasco
P
Shintani
DK
, et al.
2004
.
Conserved simple sequence repeats for the Limnanthaceae (Brassicales)
.
Theoretical and Applied Genetics
108
:
450
457
.

Koons
DN
Pavard
S
Baudisch
A
Metcalf
CJE.
2009
.
Is life-history buffering or lability adaptive in stochastic environments?
Oikos
118
:
972
980
.

Kraft
NJB
Baldwin
BG
Ackerly
DD.
2010
.
Range size, taxon age and hotspots of neoendemism in the California flora
.
Diversity and Distributions
16
:
403
413
.

Kroiss
SJ
Hillerislambers
J.
2015
.
Recruitment limitation of long-lived conifers: implications for climate change responses
.
Ecology
96
:
1286
1297
.

LaDochy
S
Medina
R
Patzert
W.
2007
.
Recent California climate variability: spatial and temporal patterns in temperature trends
.
Climate Research
33
:
159
169
.

Latta
RG
Mitton
JB.
1999
.
Historical separation and present gene flow through a zone of secondary contact in ponderosa pine
.
Evolution
53
:
769
776
.

Ledig
FT.
2000
.
Founder effects and the genetic structure of Coulter pine
.
Journal of Heredity
91
:
307
315
.

Ledig
FT
Conkle
MT.
1983
.
Gene diversity and genetic structure in a narrow endemic, Torrey pine (Pinus torreyana parry ex carr.)
.
Evolution
37
:
79
85
.

Ledig
FT
Hodgkiss
PD
Johnson
DR.
2005
.
Genic diversity, genetic structure, and mating system and Brewer spruce (Pinaceae), a relict of the Arcto-Tertiary forest
.
American Journal of Botany
92
:
1975
1986
.

Lee
SW
Ledig
FT
Johnson
DR.
2002
.
Genetic variation at allozyme and RAPD markers in Pinus longaeva (Pinaceae) of the White Mountains, California
.
American Journal of Botany
89
:
566
577
.

Lesica
P
Allendorf
FW.
1995
.
When are peripheral populations valuable for conservation?
Conservation Biology
9
:
753
760
.

Lesica
P
Allendorf
FW.
1999
.
Ecological genetics and the restoration of plant communities: mix or match?
Restoration Ecology
7
:
42
50
.

Liston
A
Parker-Defeniks
M
Syring
JV
Willyard
A
Cronn
R.
2007
.
Interspecific phylogenetic analysis enhances intraspecific phylogeographical inference: a case study in Pinus lambertiana
.
Molecular Ecology
16
:
3926
3937
.

Loarie
SR
Carter
BE
Hayhoe
K
, et al.
2008
.
Climate change and the future of California’s endemic flora
.
PLoS One
3
:
e2502
. doi:10.1371/journal.pone.0002502.e2502.

Lowe
WH
Allendorf
FW.
2010
.
What can genetics tell us about population connectivity?
Molecular Ecology
19
:
3038
3051
.

Martín-Bravo
S
Escudero
M.
2012
. Biogeography of flowering plants: a case study in mignonettes (Resedaceae) and sedges (Carex, Cyperaceae). In:
Stevens
L.
ed.
Global advances in biogeography
.
Rijeka, Croatia
:
Intech
,
257
290
.

Mashayekhi
S
Columbus
JT.
2015
.
Genetic diversity of Allium munzii (Amaryllidaceae), a rare southern California species and implication for its conservation
.
Biochemical Systematics and Ecology
59
:
91
99
.

Mastrogiuseppe
RJ
Mastrogiuseppe
JD.
1980
.
A study of Pinus balfouriana Grev. & Balf. (Pinaceae)
.
Systematic Botany
5
:
86
104
.

Memmott
J
Craze
PG
Waser
NM
Price
MV.
2007
.
Global warming and the disruption of plant-pollinator interactions
.
Ecology Letters
10
:
710
717
.

Moeller
DA
Geber
MA
Tiffin
P.
2011
.
Population genetics and the evolution of geographic range limits in an annual plant
.
American Naturalist
178
:
S44
–S61.

Moles
AT
Falster
DS
Leishman
MR
Westoby
M.
2004
.
Small-seeded species produce more seeds per square metre of canopy per year, but not per individual per lifetime
.
Journal of Ecology
92
:
384
396
.

Monleon
VJ
Lintz
HE.
2015
.
Evidence of tree species’ range shifts in a complex landscape
.
PLoS One
10
:
e0118069
. doi:10.1371/journal. pone.0118069.

Moore
DM
Chater
OA.
1971
.
Studies on bipolar species. I. Carex
.
Botaniska Notiser
124
:
317
334
.

Moritz
C
Patton
JL
Conroy
CJ
Parra
JL
White
GC
Beissinger
SR.
2008
.
Impact of a century of climate change on small-mammal communities in Yosemite National Park, USA
.
Science
322
:
261
264
.

Morris
WF
Pfister
CA
Tuljapurkar
S
, et al.
2008
.
Longevity can buffer plant and animal populations against changing climatic variability
.
Ecology
89
:
19
25
.

Muhlfeld
CC
McMahon
TE
Boyer
MC
Gresswell
RE.
2009
.
Local-habitat, watershed, and biotic factors in the spread of hybridization between native westslope cutthroat trout and introduced rainbow trout
.
Transactions of the American Fisheries Society
138
:
1036
1051
.

Myers
N
Mittermeier
RA
Mittermeier
CG
da Fonseca
GA
Kent
J.
2000
.
Biodiversity hotspots for conservation priorities
.
Nature
403
:
853
858
.

Nettel
A
Dodd
RS
Afzal-Rafii
Z.
2009
.
Genetic diversity, structure, and demographic change in tanoak, Lithocarpus densiflorus (Fagaceae), the most susceptible species to sudden oak death in California
.
American Journal of Botany
96
:
2224
2233
.

Oline
DK
Mitton
JB
Grant
MC.
2000
.
Population and subspecific genetic differentiation in the foxtail pine (Pinus balfouriana)
.
Evolution
54
:
1813
1819
.

Ovaskainen
O.
2002
.
Long-term persistence of species and the SLOSS problem
.
Journal of Theoretical Biology
218
:
419
433
.

Parmesan
C
Ryrholm
M
Stefanescu
C
, et al.
1999
.
Poleward shifts in geographical ranges of butterfly species associated with regional warming
.
Nature
399
:
579
583
.

Patton
JL
Smith
MF.
1990
.
The evolutionary dynamics of the pocket-gopher Thomonys bottae, with emphasis on California populations. University of California Publications in Zoology, Volume 123
.
Oakland
:
University of California Press
.

Patton
JL
Huckaby
DG
Alvarez-Castaneda
STL.
2008
.
The evolutionary history and a systematic revision of woodrats of the Neotoma lepida group
.
Berkeley, CA, USA
:
University of California Press
. Published to California Scholarship Online in March 2012: doi:10.1525/california/9780520098664.001.0001.

Peabody
FE
Savage
JM.
1958
. Evolution of a Coast Range corridor in California and its effects on the origin and dispersion of living amphibians and reptiles. In:
Hubbs
CL.
ed.
Zoogeography
.
Washington, DC
:
American Association for the Advancement of Science
,
159
186
.

Pelini
SL
Dzurisin
JDK
Prior
KM
, et al.
2009
.
Translocation experiments with butterflies reveal limits to enhancement of poleward populations under climate change
.
Proceedings of the National Academy of Sciences of the USA
106
:
11160
11165
.

Penman
TD
Collins
L
Syphard
AD
Keeley
JE
Bradstock
RA.
2014
.
Influence of fuels, weather and the built environment on the exposure of property to wildfire
.
PLoS One:
9
:
e111414
. doi:10.1371/journal.pone.0111414.

Peter
BM
Slatkin
M.
2013
.
Detecting range expansions from genetic data
.
Evolution
67
:
3274
3289
.

Polihronakis
M
Caterino
MS.
2010
.
Contrasting patterns of phylogeographic relationships in sympatric sister species of ironclad beetles (Zopheridae: Phloeodes spp.) in California’s Transverse Ranges
.
BMC Evolutionary Biology
10
:
195
.

Preston
RE
Dempster
LT.
2012
. Convolulaceae In:
Baldwin
BG
Goldman
DH
Keil
DJ
Patterson
R
Rosatti
TJ
Wilken
DH.
eds.
The Jepson manual: vascular plants of California
, 2nd edn.
Berkeley
:
University of California Press
,
654
664
.

Raven
PH
Axelrod
DI.
1978
.
Origin and relationships of the California flora
.
Sacramento
:
California Native Plant Society
.

Reed
DH
Frankham
R.
2003
.
Correlation between fitness and genetic diversity
.
Conservation Biology
17
:
230
237
.

Remington
CL.
1968
. Suture-zones of hybrid interaction between recently joined biotas. In:
Dobzhansky
T
Hecht
MK
Steere
WC.
eds.
Evolutionary biology
.
New York
:
Plenum
,
321
428
.

Rieseberg
LH.
1991
. Hybridization in rare plants: insights from case studies in Helianthus and Cercocarpus. In:
Falk
DA
Holsinger
KE
, eds.
Conservation of rare plants: biology and genetics
.
New York
:
Oxford University Press
,
171
181
.

Roderick
GK.
1996
.
Geographic structure of insect populations: gene flow, phylogeography, and their uses
.
Annual Review of Entomology
41
:
325
352
.

Rogers
DL
Montalvo
AM.
2004
. Genetically appropriate choices for plant materials to maintain biological diversity. USDA Forest Service Rocky Mountain Region. Region, Lakewood, CO. http://www.fs.usda.gov/Internet/FSE_DOCUMENTS/fsbdev3_039080.pdf.

Rundel
PW.
2011
.
The diversity and biogeography of the alpine flora of the Sierra Nevada, California
.
Madroño
58
:
153
184
.

Salguero-Gómez
Siewert
W
Casper
BB
Tielbörger
K.
2012
.
A demographic approach to study effects of climate change in desert plants
.
Philosophical Transactions of Royal Society of London Series B: Biological Sciences
367
:
3100
3114
.

Sawyer
JO.
2006
.
Northwestern California
.
Berkeley
:
University of California Press
.

Sawyer
JO
Thornburgh
DA.
1977
. Montane and sub-alpine vegetation of the Klamath Mountains. In:
Barbour
MG
Major J
eds.
Terrestrial vegetation of California
.
New York
:
John Wiley and Sons
,
699
732
.

Sawyer
JO
Keeler-Wolf
T
Evens
JM.
2009
.
A manual of California vegetation
, 2nd edn.
Sacramento
:
California Native Plant Society
.

Schierenbeck
KA.
2014
.
Phylogeography of California
.
Berkeley
:
University of California Press
.

Schoville
SD
Lam
AW
Roderick
GK.
2010
.
A range-wide genetic bottleneck overwhelms contemporary landscape factors and local abundance in shaping genetic patterns of an alpine butterfly (Lepidoptera: Pieridae: Colias behrii)
.
Molecular Ecology
21
:
4242
4256
.

Schweiger
O
Settele
J
Kudrna
O
Klotz
S
Kühn
I.
2008
.
Climate change can cause spatial mismatch of trophically interacting species
.
Ecology
89
:
3472
3479
.

Scott
RW.
1995
.
The alpine flora of the Rocky Mountains: Volume 1. The Middle Rockies
.
Salt Lake City
:
University of Utah Press
.

Serra-Diaz
JM
Keenan
TF
Ninyerola
M
Sabat
S
Gracia
C
Lloret
F.
2013
.
Geographical patterns of congruence and incongruence between correlative species distribution models and a process-based ecophysiological growth model
.
Journal of Biogeography
40
:
1928
1938
.

Serra-Diaz
JM
Scheller
RM
Syphard
AD
Franklin
J.
2015
.
Disturbance and climate microrefugia mediate tree range shifts during climate change
.
Landscape Ecology
30
:
1039
1053
.

Sexton
JP
Dickman
EE.
2016
.
What can local and geographic population limits tell us about distributions?
American Journal of Botany
103
:
129
139
.

Sexton
JP
Strauss
SY
Rice
KJ.
2011
.
Gene flow increases fitness at the warm edge of a species' range
.
Proceedings of the National Academy of Sciences of the USA
108
:
11704
11709
.

Sexton
JP
Hangartner
SB
Hoffmann
AA.
2014
.
Genetic isolation by environment or distance: which pattern of gene flow is most common?
Evolution
68
:
1
15
.

Shaffer
HB
Pauly
GB
Oliver
JC
Trenham
PC.
2004
.
The molecular phylogenetics of endangerment: cryptic variation and historical phylogeography of the California tiger salamander, Ambystoma californiense
.
Molecular Ecology
13
:
3033
3049
.

Sharsmith
CW.
1940
. A contribution to the history of the alpine flora of the Sierra Nevada. PhD Thesis, University of California, Berkeley.

Slatkin
M
Hudson
RR.
1991
.
Pairwise comparisons of mitochondrial DNA sequences in stable and exponentially growing populations
.
Genetics
129
:
555
562
.

Sloop
CM
Pickens
C
Gordon
SP.
2011
.
Conservation genetics of Butte County meadowfoam (Limnanthes floccosa ssp. californica Arroyo), an endangered vernal pool endemic
.
Conservation Genetics
12
:
211
323
.

Smith
CI
Tank
S
Godsoe
W
, et al.
2011
.
Comparative phylogeography of a coevolved community: concerted population expansions in Joshua trees and four yucca moths
.
PLoS One
6
:
e25628
. doi:10.1371/journal.pone.0025628.

Smith
JP
Sawyer
JO
Jr.
1988
.
Endemic vascular plants of northwestern California and southwestern Oregon
.
Madroño
35
:
54
69
.

Soltis
DE
Gitzendanner
MA
Strenge
DD
Soltis
PS.
1997
.
Chloroplast DNA intraspecific phylogeography of plant from the Pacific Northwest of North America
.
Plant Systematics and Evolution
206
:
353
373
.

Sork
VL
Davis
FW
Westfall
R
, et al.
2010
.
Gene movement and genetic association with regional climate gradients in California valley oak (Quercus lobata Née) in the face of climate change
.
Molecular Ecology
19
:
3806
3823
.

State of California, Department of Finance
.
2014
. Report P-1 (age): state and county population projections by major age groups, 2010-2060. Sacramento, CA.

Stebbins
GL
Major
J.
1965
.
Endemism and speciation in the California flora
.
Ecological Monographs
35
:
1
35
.

Stevens
VM
Trochet
A
Van Dyck
H
Clobert
J
Baguette
M.
2012
.
How is dispersal integrated in life histories: a quantitative analysis using butterflies
.
Ecology Letters
15
:
74
86
.

Steinhoff
RJ
Joyce
DG
Fins
L.
1990
.
Isozyme variation in Pinus monticola
.
Canadian Journal of Forest Research
13
:
1122
1132
.

Syring
J
Farrell
K
Businský
R
Cronn
R
Liston
A.
2007
.
Widespread genealogical nonmonophyly in species of Pinus subgenus Strobus
.
Systematic Biology
56
:
163
181
.

Syphard
AD
Keeley
JE.
2015
.
Location, timing and extent of wildfire vary by cause of ignition
.
International Journal of Wildland Fire
24
:
37
47
.

Tamme
R
Gotzenberger
L
Zobel
M
, et al.
2014
.
Predicting species’ maximum dispersal distances from simple plant traits
.
Ecology
95
:
505
513
.

Templeton
AR.
1986
.
Further comments on the statistical analysis of DNA-DNA hybridization data
.
Molecular Biology and Evolution
3
:
290
295
.

Thibert-Plante
X
Hendry
AP.
2010
.
When can ecological speciation be detected with neutral loci?
Molecular Ecology
19
:
2301
2314
.

Thomson
FJ
Moles
AT
Auld
TD
Ramp
D
Ren
SQ
Kingsford
RT.
2010
.
Chasing the unknown: predicting seed dispersal mechanisms from plant traits
.
Journal of Ecology
98
:
1310
1318
.

Truesdale
HD
McClenaghan
LR
Jr.
1998
.
Population genetic structure of tecate cypress (Cupressaceae)
.
Southwestern Naturalist
43
:
363
373
.

Vanderwel
MC
Purves
DW.
2014
.
How do disturbances and environmental heterogeneity affect the pace of forest distribution shifts under climate change?
Ecography
37
:
10
20
.

Verhoeven
KJF
Macel
M
Wolfe
LM
Biere
A.
2011
.
Population admixture, plant invasions and the balance between local adaptation and inbreeding depression
.
Proceedings of the Royal Society of London Series B: Biological Sciences
278
:
2
8
.

Vilà
M
Weber
E
D’Antonio
CM.
2000
.
Conservation implications of invasion by plant hybridization
. Biological Invasions
2
:
207
217
.

Villa-Lobos
J.
2003
.
Klamath-Siskiyou: jewel of the Pacific Coast
.
Plant Talk
31
:
29
33
.

Vitt
P
Havens
K
Kramer
AT
Sollenberger
D
Yates
E.
2010
.
Assisted migration of plants: changes in latitudes, changes in attitudes
.
Biological Conservation
143
:
18
27
.

Warren
MS
Hill
JK
Thomas
JA
, et al.
2001
.
Rapid responses of British butterflies to opposing forces of climate and habitat change
.
Nature
414
:
65
69
.

Whittaker
RH.
1961
.
Vegetation history of the Pacific Coast States and the “central” significance of the Klamath Region
.
Madroño
16
:
5
23
.

Willi
Y
Van Buskirk
J
Schmid
B
Fischer
M.
2006
.
Genetic isolation of fragmented populations is exacerbated by drift and selection
.
Journal of Evolutionary Biology
20
:
534
542
.

Wolf
AT
Howe
RW
Hamrick
JL.
2000
.
Genetic diversity and population structure of the serpentine endemic Calystegia collina (Convolvulaceae) in Northern California
.
American Journal of Botany
87
:
1138
1146
.

Zacherl
D
Gaines
SD
Lonhart
SI.
2003
.
The limits to biogeographical distributions: insights from the northward range extension of the marine snail Kelletia kelletii (Forbes, 1852)
.
Journal of Biogeography
30
:
913
924
.