Evolutionary Hotspots in the Mojave Desert

Genetic diversity within species provides the raw material for adaptation and evolution. Just as regions of high species diversity are conservation targets, identifying regions containing high genetic diversity and divergence within and among populations may be important to protect future evolutionary potential. When multiple co-distributed species show spatial overlap in high genetic diversity and divergence, these regions can be considered evolutionary hotspots. We mapped spatial population genetic structure for 17 animal species across the Mojave Desert, USA. We analyzed these in concurrence and located 10 regions of high genetic diversity, divergence or both among species. These were


Introduction
The genetic diversity within a species represents a fundamental level of biodiversity.It provides the basis for phenotypic variation and adaptation, and underlies the evolutionary processes of lineage diversification and speciation that contribute to the patterns of species-, community-and ecosystemlevel biodiversity evident today.To protect species-level biodiversity into the future, evolutionary processes must also be protected [1][2][3].By examining and mapping patterns of intraspecific genetic diversity across a landscape, regions with high evolutionary potential can be identified, where diversification and speciation may be more likely to occur [4][5][6].Herein, we characterize evolutionary potential in terms of interpopulation genetic divergence and intrapopulation genetic diversity.Geographic areas that support divergent populations may harbor greater evolutionary potential, as they can represent places where divergent lineages now meet (suture zones) or reflect abiotic drivers of adaptive variation (e.g., an ecotone of steep environmental gradients).High levels of intrapopulation genetic diversity may also provide a rich resource for evolutionary resilience, and can reflect large effective population sizes [7][8][9], the presence of a historical refuge with environmentally stable conditions [10], or a zone of mixing between previously isolated gene pools [11,12].Conversely, low or reduced genetic diversity is generally associated with reductions in fitness and survival [13,14].Comparative phylogeographic analyses have repeatedly shown that geographic barriers or past geological processes and climate alterations have left similar signatures on the genetic structure of populations across entire communities [15][16][17][18][19][20].When regions of high genetic divergence and diversity overlap for multiple species, these can be considered to be evolutionary hotspots.Protection of these hotspots may be important for preserving the process of lineage diversification and may impart resilience to changing environmental conditions and selective pressures [21,22].
Determining the locations of evolutionary hotspots may be particularly informative in regions where development pressures are growing, and land use planning efforts include conservation goals.In the Mojave Desert for example, widespread development of utility scale renewable energy facilities, transmission corridors, and associated infrastructure has been proposed and is being implemented, and land use planning strategies that balance renewable energy production with other land use and conservation goals are currently being devised.The Mojave Desert is the smallest of the four North American deserts, and is generally characterized by broad valleys with alluvial fan complexes, playas and dunes that are separated by plateaus and rugged mountain ranges.Major vegetation communities in the valleys are dominated by mixed shrublands while mountain slopes are dominated by woodlands [23][24][25].
The biogeographic history of the Mojave is complex.The modern desert environment in the Mojave is young; however, the first desert conditions originated approximately 5 million years ago (Ma) [26] and have fluctuated since that time with changing climatic conditions.Warm desert vegetation may only have become extensive through the region over the last 12,000 years [27].However, in a recent analysis of multiple arid-adapted small animal species in the Mojave and Sonoran Deserts, we found that many were present prior to this time, with divergence estimates for Mojave lineages ranging from 0.5 to 11.5 Ma [20], and indicating that genetic structure within desert adapted animal species has likely been influenced by expanding and contracting habitat and other isolating factors that have developed over this time period.In addition to climatic fluctuations, inundation and mountain uplift have influenced the Mojave Desert biota through time.In the late Miocene-late Pliocene (10-4 Ma) inundation events flooded the Gulf of California through the Colorado Desert region and contributed to the early Colorado River system, isolating the Mojave Desert from the Sonoran Desert [28,29].Uplift of the Sierra Nevada and Transverse Mountain Ranges, between 5-2 Ma, reversed the flow of the Mojave River to drain into the Mojave Basin and created rain shadows to the east.During the Pleistocene, between 1.2-0.5 Ma, the Mojave River slowly advanced, filling playas and basins.During glacial-interglacial cycles that followed, the Mojave oscillated between a lacustrine and dry system.During the last glacial maximum, the Mojave region was dominated by a fluvial-lacustrine system [30] and arid adapted taxa were likely isolated in small refugia [31].Piñon-juniper woodlands dominated the region [32], but eventually retreated as warming and drying trends intensified during the Holocene.
In modern times, human land use and modification in the arid lands of the western United States have further altered ecosystems, and have intensified in recent decades [33].In the Mojave Desert, urbanization, agriculture, transportation networks, military operations and energy production have already transformed sections of the landscape with impacts to wildlife [34,35].Further land use change is projected in this region, particularly utility-scale renewable energy development (USRED) on federal lands in response to federal and state goals and mandates for renewable energy production (e.g., the American Recovery and Reinvestment Act of 2009, CA Executive Order S-14-08 and CA Senate Bill X1-2).Increasing production of renewable energy will have benefits in curtailing greenhouse gas emissions, and promoting energy independence and economic growth, but development of USRED facilities, added transmission corridors, and their associated infrastructure may negatively impact local wildlife through habitat loss and fragmentation.Given the projected scale and rate of USRED development in the desert southwest, land use planning efforts must balance energy production and delivery goals with other land use and conservation priorities (e.g., California Desert Renewable Energy Conservation Plan).Additionally, agencies tasked with natural resource and wildlife management must understand potential impacts of planned development to local species and systems for mitigation and adaptive management purposes.
The focus of this study was to map patterns of intraspecific genetic divergence and diversity for multiple species throughout the Mojave Desert to locate evolutionary hotspots, where multiple species may have high evolutionary potential.We evaluated whether identified hotspots fell within protected lands, or overlapped with current and proposed USRED project sites and utility corridors.Resulting GIS-based maps may inform current landscape-based planning efforts, and help direct future investigations to particular regions of interest in the Mojave Desert and bordering ecoregions.

Study Region
This study was conducted in the Mojave Desert, using the spatial definition from the U.S. Environmental Protection Agency (EPA) Mojave Basin and Range Level III Ecoregion [36].The Mojave ranges in elevation from 146 m below sea level (Death Valley, CA, USA) to over 3,633 m (Mount Charleston, NV, USA).Precipitation in the Mojave ranges from 100 to 350 mm per year, and shows a strong longitudinal trend, with the western portions of the Mojave receiving more precipitation in the winter, and the eastern portions receiving increased precipitation in late summer.Temperatures range from below 0 °C in the winter months to over 54 °C in the summer months, and vary diurnally and geographically.We buffered the EPA Mojave Basin and Range Ecoregion by 10 km to ensure that our analyses covered all of the Mojave Desert, resulting in a study area of approximately 150,000 km 2 (Figure 1).
Genetic data consisted mainly of mitochondrial sequence data, which were available for all species, with nuclear sequence data or microsatellite data available for a subset of six species (Table 1).Alignments for all sequence datasets were obtained via CLUSTAL-W [54], as implemented in MEGA 4.0 [55].Sequence data were analyzed in ARLEQUIN 3.5 [56].For consistency, all sequence-based analyses were conducted assuming a Tamura and Nei [57] model of nucleotide evolution.For microsatellite data, allelic richness and heterozygosity were estimated in MSA [58].Whenever possible we incorporated data from multiple loci, as there are limitations to inferring population structure and diversity from mtDNA or any single locus [59,60].We also relied upon multiple taxa rather than multiple loci to cross validate observed patterns.This approach is consistent with our goal to uncover regional diversity patterns.

Mapping Divergence Landscapes
We used a phylogeographic approach to map patterns of genetic divergence estimated from mtDNA (and when available nuclear) sequences within species.For each species, genetic divergence landscapes were mapped in a Geographic Information System (GIS) based on pairwise genetic divergence among collection locations defined by the original study authors [6].Sequence divergence was calculated as D A , Nei's net number of nucleotide differences between each pair of collection locations, where the average number of differences between individuals among locations is corrected for the average number of differences within each location [61].For datasets where multiple loci were sequenced, pairwise distances were averaged across loci after normalization.We tested for a signal of increasing genetic differentiation with increasing geographic distance (or Isolation by Distance; IBD) among sampling locations using reduced major axis regression as implemented in IBDWS 3.16 [62].For species where statistically significant IBD was detected, we used the regression residuals rather than raw genetic divergence values to create divergence landscapes, in order to highlight locations where divergence was high while accounting for variation in the geographic distance among sampling points.Genetic distances or residual distances were visualized as genetic landscapes in ArcGIS 9.3, using the Genetic Landscapes Toolbox [63].Divergence landscapes were created using the Single Species Divergence Tool.Pairwise genetic divergence values (or residuals) were mapped to the geographic midpoints between collection locations, as described and employed in previous visualization methods [64,65].A surface was interpolated from the midpoints using Inverse Distance Weighted interpolation (power = 2, variable search radius with 12 points, grid cell size 1 km 2 ).To avoid extrapolating beyond the original collection locations, individual species surfaces were clipped to the spatial extent of collection locations.
Finally, we masked interpolated genetic landscapes with suitable habitat derived from habitat suitability models created or obtained for each species.Functionally, this allowed us to focus on locations where genetic divergence was high across suitable (and potentially occupied) habitat, highlighting regions of potential lineage and population re-contact (or suture zones) rather than current geographic barriers.For all but one species, (O.canadensis), habitat suitability was statistically modeled from presence data using generalized additive modeling, maximum entropy (MaxEnt) [66], or both.Full modeling methods and results are available elsewhere [67][68][69], and so are only generally summarized here.Species occurrence records were obtained from online public databases for museum collections (VertNet, HerpNET, and MaNIS), and other regional data sources.For each species, conceptual models of potential habitat requirements were developed to identify the environmental covariates thought to influence their geographic distribution.These covariate layers were represented in a GIS as spatial raster layers, with a 1km 2 grid cell size.Model performance was evaluated using the receiver operator curve statistic (AUC) [70,71], Boyce Index (BI) [72], Akaike Information Criterion (AIC) [73,74] and Bayesian Information Criterion (BIC) [75] test statistics, and environmental response curves were visually assessed.Models that ranked high based on these criteria were then averaged into a single final habitat suitability map for each species.Finally, we reclassified the averaged habitat suitability maps into binary representations of suitable/unsuitable habitat using the 5th percentile of habitat suitability scores from the occurrence localities [76].For O. canadensis, a statistical habitat suitability model encompassing the study region was not available at the time of analysis.An expert opinion potential habitat model [77], rasterized at a 1 km 2 grid cell size, was used to mask genetic landscapes instead.
To highlight areas of congruence among species, the genetic divergence landscapes for all taxa were averaged into a single multi-species genetic landscape.To assure that each species received equal weighting in the multi-species genetic landscape, we re-scaled individual species divergence landscapes by dividing each cell value by the maximum.Because the spatial coverage of available genetic data varied among species, the number of species represented in each 1 km 2 grid cell of the multi-species genetic landscape also varied.To assess multi-species concordance, we clipped the multi-species genetic landscape extent to areas with coverage for three or more species.We considered the most divergent cells in the multi-species genetic landscape to be greater than 1.5 standard deviations from the mean cell value [6].Single species genetic divergence landscapes were then visually examined in relation to the average landscape to determine which species contributed to divergent areas.To further assess variation across individual species landscapes that contributed to the average landscape, we also calculated the coefficient of variation among single species landscapes, calculated as the standard deviation divided by the mean.

Creating Diversity Landscapes
Diversity landscapes were created from genetic diversity estimates within sampling locations.To maximize the number of individuals grouped per sampling location, we grouped sampling locations that were within 20 km of one another that were not divided by potential barriers (e.g., rivers, mountains, highways, etc.).New population point locations were assigned as the midpoints of grouped sampling locations.All remaining singleton sampling locations were excluded from further analyses.Genetic diversity was calculated in two ways.For sequence data, we calculated -sequence diversity‖ as the average sequence divergence among individuals within a population as π i , under the Tamura and Nei model of nucleotide evolution [57].This statistic is less biased by differences in sample size among populations than others such as the number of alleles or segregating sites [78].For species for which the majority of sampling locations contained three or more individuals, we also calculated -gene diversity‖ for sequence datasets [79].For haplotypic data this statistic is defined as the probability that two randomly chosen haplotypes are different in the sample, and is equivalent to expected heterozygosity for diploid data.For microsatellite datasets we calculated diversity as expected heterozygosity (H E ).Gene diversity and H E are different from sequence diversity in that they focus on the number of gene copies, but do not contain information on the relative divergence of gene copies.For species where mtDNA sequences and multiple microsatellite loci were available, we chose to model gene diversity from microsatellites (H E ) as multiple loci should give a better representation of levels of genomic diversity within a species.
Diversity landscapes were created using IDW interpolation in the Genetic Landscapes GIS Toolbox.Diversity values are specific to their population or collection location and so were mapped to their population locations, rather than to the midpoints between locations as done in divergence mapping.As with divergence landscapes, individual species surfaces were clipped to the extent of the collection locations, and suitable habitat was used to bound the analyses.Finally, diversity landscapes (sequence diversity, or gene diversity/H E ) were averaged among species and the coefficient of variation was also calculated.Single species genetic diversity landscapes were visually examined to determine which species contributed to t he most divergent areas on the average landscapes.

Assessing Protected Status and Vulnerability
Average genetic diversity and divergence surfaces were assessed in combination with land status boundaries and USRED boundaries to determine if regions of high genetic divergence and diversity fell within lands we considered protected or at-risk, or within USRED project footprints or transmission corridors.A land status layer was compiled for the study area using a regional land management/ownership map categorized into three status levels: -Protected‖, -Uncertain‖ or -At-Risk‖.The land management and ownership maps were obtained from the Bureau of Land Management (BLM) for California, Nevada, Arizona and Utah.We considered areas where the administrative status of land parcels was generally mandated to protect or preserve natural resources and ecosystems present as -Protected‖.Land parcels that we categorized as -At-Risk‖ had administrative mandates that allow land uses that are incompatible with preserving habitat for our focal species-even if portions of those parcels are currently undisturbed or occupied by our focal species.
Finally, the -Uncertain‖ status designates land parcels for which the effect of administrative status is currently unknown.Lands that we categorized as -Protected‖ included: BLM-Areas of Critical Environmental Concern and National Conservation Areas, U.S. Fish and Wildlife Service refuges, National Park Service lands, State wildlife departments, and State Parks and Wilderness Areas.Lands having -Uncertain‖ status included: U.S. Department of Energy, State lands, U.S. Forest Service lands, and Wilderness Study Areas.And finally, lands that we categorized as -At-Risk‖ were: BLM Unclassified lands (lands that do not have special designations other than multiple use), Bureau of Reclamation, Department of Defense, local municipality and county government, Off-Highway Vehicle areas, and private property.
Shapefiles of current and proposed USRED facilities mapped by BLM in California, Nevada and parts of Arizona and Utah were provided by the Renewable Energy Project Managers for the Southern Nevada and the California Desert Districts of the BLM, or were obtained from the Solar Energy Development Programmatic Environmental Impact Statement [80].Footprints of pending USRED identified by BLM were compiled only for BLM lands, while development on private or state owned land are not represented.These data are not readily available, and were not included in our analysis.We also compiled areas that are designated for transmission corridors under the California Desert Conservation Area Plan of 1980 [81], and the West-wide Designation of Energy Corridors [82].All files were converted into raster format at a 1 km 2 grid cell size and overlaid on average genetic landscapes, and to regions of highest divergence and diversity (categorized as 1.5 standard deviations from mean values).Because USRED project sites and attendant transmission lines will likely have different impacts on wildlife, we examined these layers separately and in combination.

Evolutionary Hotspots
The average divergence landscape encompassed 110,089 km 2 after clipping to the extent with three or more species overlapping, representing 73% of the study area.Data were lacking in the most northern and western portions of the study region (Figure 2A).The most divergent areas (categorized as greater than 1.5 standard deviations above the mean) encompassed 4,864 km 2 or 4.4% of the analyzed area.We expected highest areas of divergence across habitat transition zones and our results generally concur with this, as hotspots generally ring the ecoregion.We identified 10 locations with highest levels of genetic divergence (Table 2; Figure 3 labeled regions A-J).These hotspots varied in the number and composition of contributing individual species (Table 2).Concordance was greatest in the Colorado River hotspot, with 10 species showing high divergence in this region, while the Ivanpah Valley and Virgin Mountains hotspots reflect high divergence in only one species each, G. agassizii and A. punctatus, respectively.However, the Ivanpah hotspot was also a region of high gene diversity in several other species examined.While relatively few species had high sequence divergence or diversity in the Virgin Mountains, this region has been previously identified as a zone of secondary contact and hybridization in amphibians [37], a group that is only represented by one species in our dataset.Across all species, variation among datasets was generally highest in the central portion of the study region.(Figure 2A inset).

C
Sequence diversity (π) was averaged across 14 species (A.punctatus, C. occipitalis, C. penicillatus, G. agassizii, H. selenopoides, H. theologus, L. trivirgata, O. canadensis, P. longimembris, P. gilberti, S. magister, U. scoparia, X. vigilis, and X. mohavensis).The remaining three species in our dataset (C.bicinctores, D. dorsalis, and T. bottae) were excluded from diversity analyses because only single individuals were sampled at disparate locations that could not be combined.The average sequence diversity landscape covered 87,982 km 2 when clipped to include three or more species (Figure 2B).High sequence diversity hotspots encompassed 6,724 km 2 , or approximately 7.6% of the area.Sequence diversity hotspots were mainly concentrated along the western and southern boundaries of the Mojave Ecoregion and overlapped with regions of high divergence (Table 2, Figure 3).Individual species showed the greatest concordance at the Ord Mountain -Lucerne Valley hotspot with 5 species showing high sequence diversity (Table 2).Variation among species was generally greatest in the northeast portion of the study region, with patches throughout the central and southwest Mojave Desert (Figure 2B inset).
Gene diversity (heterozygosity) was averaged across 13 species (A.punctatus, C. occipitalis, C. penicillatus, G. agassizii, H. selenopoides, H. theologus, L. trivirgata, O. canadensis, P. longimembris, P. gilberti, S. magister, U. scoparia, X. vigilis, and X. mohavensis).Sceloporus magister was excluded because most sampling locations contained less than three individuals.The gene diversity average genetic landscape covered 71,468 km 2 when clipped to the extent of three or more species and showed some geographical concordance with sequence divergence and diversity (Figure 2C).
Highest diversity was concentrated along the western and southern ecoregional transition zones (Table 2, Figure 3).In addition, high average gene diversity was evident in the central Mojave near the Ivanpah Valley (Figure 4).The Indio Hills and Pluvial lakes hotspots represented the most species, with six species contributing to each (Table 2).Variation among datasets was generally greatest in the north central and northeastern Mojave Desert (Figure 2C   The complex climatic and geological history of the Mojave Desert likely contributed to the formation of evolutionary hotspot regions detected.Hotspots are clustered in locations along the western and southern regions of the Mojave Desert, where past inundation, uplift, and the cyclical formation of riverine and lacustrine systems may have isolated lineages [50].These areas may represent secondary contact zones or -suture zones‖ for previously isolated lineages across multiple species [83,84].Regions of high gene diversity and heterozygosity may also reflect large effective population sizes [7][8][9]85] or refugia [10].In addition to historical isolating events, current environmental conditions may contribute to diversity and divergence hotspots in regions with steep environmental gradients (e.g., steep ecotones between mountains and basins).
This study focused on patterns of genetic variation that are presumed to be selectively neutral.Clearly, the ultimate measure of evolutionary potential is genetic variation underlying traits that will be under selection in the future.Although higher intrapopulation genetic diversity measured at neutral loci can be associated with higher fitness and lower extinction rates [14,86], correlations between neutral and adaptive variation are not always strong [87,88].However, inferences of adaptive potential are strengthened when the genetic hotspots overlap spatially with zones of lineage recontact and/or steep environmental gradients.Hybridization between previously isolated and divergent lineages or populations can create novel gene combinations that can facilitate speciation and adaptive evolution in some cases [11,[89][90][91][92][93].Adaptive variation is also often concentrated across ecological transitions [94][95][96].Although not included in this study, previous genetic and morphological analysis of the Neotoma lepida group revealed zones of lineage recontact and hybridization between coastal and desert morphological groups in Kelso Valley (near hotspot B) and in Morongo Valley (hotspot E) [97].Likewise, in the Mojave Desert annual Linanthus parrayae, dimorphism in flower color occurs in local populations found within hotspot D near Lucerne Valley and Victorville, and along the base of the San Gabriel Mountains near Palmdale [98].Flowers are white throughout the majority of the species range, but blue flowers occur with variable frequencies in these mixed populations.Color dimorphism seems to be maintained by temporally variable selection associated with annual rainfall patterns and differential water use [99,100].The adaptive potential stored in ecotonal evolutionary hotspots is likely to become increasingly important as climatic conditions change in the future, and environmental gradients intensify, weaken or shift spatially [22,101,102].Other studies have also highlighted regions of the Mojave.Our broader analysis encompassing both the Mojave and Sonoran Deserts and including a subset of the species considered here also highlighted hotspots of high neutral diversity and divergence near Ord Mountain-Lucerne Valley, along the Mojave/Sonoran ecotone and across the Colorado River (corresponding to hotspots D, E, F and G in Figure 3) [20].Finally, in an independent analysis focusing on California plants, Kraft et al. [5] recognized the Mojave as an important region of evolutionary potential in California, as it contains some of the youngest neoendemic vascular plants.
Because our genetic landscapes are interpolations from point data, the resulting patterns are highly dependent on the number and dispersion of collection locations across the landscape.When individual species datasets vary in sampling location and density (as with those compiled for this study), undoubtedly, uncertainty in the spatial location of hotspots is introduced.With this in mind, we caution against using the resulting hotspots maps for pinpointing exact locations for conservation purposes, rather we feel their utility is greatest for more broadly identifying regions of evolutionary potential within the ecoregion.
Mapped USRED project footprints encompassed approximately 3,547 km 2 across the study region, while transmission corridors occupied approximately 14,483 km 2 .The area of overlap of average divergence and diversity landscapes with USRED project footprints ranged from 2,563 to 3,209 km 2 and overlap with transmission corridors ranged from 8,503 to 10,733 km 2 (Table 3).USRED footprints overlapped with 3-7% of the area designated as divergence or diversity hotspots, but this rose to 10-17% with the inclusion of transmission corridors (Table 3).Certain hotspot regions overlapped more than others.Hotspots B (Sierra-Tehachapi Transition Zone), D (Ord Mountain-Lucerne Valley), F (Pluvial Lakes) and G (Colorado River) showed the most overlap with project footprints, while Hotspots C (Antelope Valley-Mojave Transition) and I (Ivanpah Valley), showed the most overlap with transmission corridors.In total, 6 of the 10 identified divergence and diversity hotspot regions have the potential to be impacted by energy and infrastructure development (Figure 4).Because evolutionary hotspots tend to occur in ecological transition zones, they may not necessarily be included in national and state parks, or other protected lands aimed at preserving exemplars of ecoregions and geomorphic provinces.We found that greater than half of the total area identified as divergence and diversity hotspots fell outside of designated protected lands.While less than 10% of the total area identified as diversity and divergence hotspots overlapped with current and pending USRED project sites, four hotspots showed substantial spatial overlap.When transmission corridors were included, overlap with hotspots increased up to 17%.Mapped energy corridors have high overlap with two additional hotspots in our study area.Given the potential impact of USRED in these areas, these six evolutionary hotspots may deserve further investigation, both in terms of habitat use and fine scale genetic structure across these regions, and in terms of specific impacts of USRED on wildlife populations in these regions.
While we have evaluated the spatial overlap between hotspot regions and planned renewable development, we do not have data on the extent to which these developments will directly impact individuals and populations in these locations.Perhaps of greatest importance to maintaining the functionality of evolutionary hotspots is mitigating the potential reduction of population size and connectivity due to habitat loss and fragmentation within and surrounding development footprints.Habitat disturbance within USRED project footprints may vary depending on the number and size of wind turbines or solar arrays placed at a site, vegetation clearing, access roads, fencing, and other project requirements [103].Many of the mapped USRED projects currently slated for development within hotspot regions will generate wind energy.The documented effects of large scale wind energy farms on wildlife are varied.Research has largely focused on direct mortality of birds and bats from collisions [103,104], with impacts varying by life history and ecology of species [105,106], wind turbine type [106], and other specific site characteristics (e.g., region, layout design of the site, topology, weather, and lighting [107]).Direct effects on ground dwelling animals are largely unknown.Behavioral changes due to turbine noise have been documented in ground squirrel populations [108], but no differences were noted in individual growth rates and mortality in a population of desert tortoises on a wind farm when compared to other locations [109].The effects of solar development projects on wildlife have not been well documented in the scientific literature [but see 80,110], and will likely vary extensively with the type of system and specific site management practices, although it should be noted that sites developed in the Mojave appear to have high levels of impact (e.g., blading with removal of all topsoil and vegetation, fencing, etc.).Associated USRED infrastructure including transmission corridors and road networks can affect connectivity [111,112], and may have a greater impact on connectivity in the Mojave than production sites themselves.Transmission corridors and surrounding right-of-ways throughout the Mojave Desert could encompass up to four times as much area as project footprints alone.In addition these crisscross through both protected and unprotected lands, further fragmenting the Mojave Desert.Over the long term, fragmentation and isolation may lead to loss of genetic diversity and increased divergence among sites, however this may take many generations to be measurable [113,114].While transmission corridors and associated road networks may represent barriers to movement for some species, they may be permeable or even facilitate movement for others [115][116][117].Variable disturbance characteristics can influence permeability (e.g., corridor width, grading, vegetation removal, paving, lighting, fencing, culverts, berms and traffic volume (reviewed in [103,112,118]).Transmission corridors and roads may also act as conduits for exotic and invasive species [119,120].For example, power lines, roads and other linear right-of-ways provide nesting and perching sites for predatory birds, including the common raven, Corvus corax [121][122][123].Large raven populations subsidized by human development in the Mojave Desert and elsewhere pose a threat to juvenile desert tortoises and other sensitive prey species [124].Increased perch availability in xeric habitats has been studied elsewhere and shown to have a negative impact on certain prey species of lizards [125].
Future survey and genetic research efforts may better resolve patterns of genetic diversity in the Mojave Desert, and refine our assessment of evolutionary hotspots.Our initial work was opportunistic, relying mainly on previously conducted population genetic and phylogeographic studies.Increasing genetic sampling of additional species representing different ecotypes or species guilds, and gathering genomic data representative of both functional and selectively neutral diversity could provide greater resolution of population-level patterns across the landscape, and determine whether the hotspots identified here are indicative of a wider range of species, and how well these reflect patterns of adaptive genetic diversity.The lack of coverage in the northern Mojave represents a significant data gap.Because hotspots tend to occur at ecotones, the northern transition between the Mojave and Great Basin may also retain high genetic diversity.This region may be particularly important if climate change results in northward range shifts for some species [126,127].Finally, focused surveys aimed at determining population status, fine-scale habitat suitability, and movement corridors within identified hotspot regions could further assess the potential site-specific impacts of USRED and other development to wildlife populations in these regions.

Conclusions
Throughout the Mojave Desert, we identified 10 regions where intraspecific genetic divergence and diversity were high for up to 17 species.Because genetic diversity underlies differentiation and adaptation, preserving regions with the greatest levels of genetic diversity and differentiation may help to preserve the evolutionary potential of these species.Thus we recognize these regions as evolutionary hotspots.While our analyses are coarse, we can generally assess where gaps in protection exist and where overlap with development may occur at a landscape-scale.Six of the 10 identified evolutionary hotspots show substantial overlap with existing and proposed renewable energy development sites and infrastructure.These results provide a working hypothesis for identifying regions of evolutionary significance and thus conservation importance in the Mojave.The compiled landscape-level, GIS-based maps of genetic information could be applied in ongoing regional conservation planning initiatives such as the Desert Renewable Energy Conservation Plan [128], the Landscape Conservation Cooperatives [129], and the U.S. Department of the Interior, Bureau of Land Management's Landscape Approach to Managing Public Lands [130].

Figure 1 .
Figure 1.Mojave Desert study region including major geographic features.

Figure 2 .BFigure 2 .
Figure 2. (A): Average sequence divergence for 17 animal species; (B): Average sequence diversity for 14 animal species; (C): Average gene diversity for 13 animal species.The insets in each figure show the coefficient of variation among individual species layers included in each average. inset).

Figure 4 .
Figure 4.Ten regional hotspots of sequence divergence, sequence diversity and gene diversity overlaid with existing and pending utility scale renewable energy development (USRED) project site footprints and energy corridors.Impervious surfaces represent existing urban development.

Table 1 .
Species and genetic datasets examined.The number of sampling locations within the study region and mean number of individuals sampled per location are given for divergence and diversity estimates.Numbers for gene diversity estimates based on microsatellite markers are provided in parentheses.NI denotes not included.

Table 2 .
Identified divergence and diversity hotspots and individual species contributing to each hotspot.Species codes refer to the first two letters of the genus and species name (full names provided in the table footer).Geographic and geological factors that may have contributed to high levels of genetic divergence and diversity within each region are also listed.

Table 3 .
Overlap between divergence and diversity landscapes, land status, proposed utility scale renewable energy development (USRED), and transmission corridors.