Water availability drives signatures of local adaptation in whitebark pine (Pinus albicaulis Engelm.) across fine spatial scales of the Lake Tahoe Basin,USA

Patterns of local adaptation at fine spatial scales are central to understanding how evolution proceeds, and are essential to the effective management of economically and ecologically important forest tree species. Here, we employ single and multilocus analyses of genetic data (n = 116 231 SNPs) to describe signatures of fine‐scale adaptation within eight whitebark pine (Pinus albicaulis Engelm.) populations across the local extent of the environmentally heterogeneous Lake Tahoe Basin, USA. We show that despite highly shared genetic variation (F_ST = 0.0069), there is strong evidence for adaptation to the rain shadow experienced across the eastern Sierra Nevada. Specifically, we build upon evidence from a common garden study and find that allele frequencies of loci associated with four phenotypes (mean = 236 SNPs), 18 environmental variables (mean = 99 SNPs), and those detected through genetic differentiation (n = 110 SNPs) exhibit significantly higher signals of selection (covariance of allele frequencies) than could be expected to arise, given the data. We also provide evidence that this covariance tracks environmental measures related to soil water availability through subtle allele frequency shifts across populations. Our results replicate empirical support for theoretical expectations of local adaptation for populations exhibiting strong gene flow and high selective pressures and suggest that ongoing adaptation of many P. albicaulis populations within the Lake Tahoe Basin will not be constrained by the lack of genetic variation. Even so, some populations exhibit low levels of heritability for the traits presumed to be related to fitness. These instances could be used to prioritize management to maintain adaptive potential. Overall, we suggest that established practices regarding whitebark pine conservation be maintained, with the additional context of fine‐scale adaptation.     

Conserving threatened riparian ecosystems in the American West: Precipitation gradients and river networks drive genetic connectivity and diversity in a foundation riparian tree (Populus angustifolia)

Gene flow is an evolutionary process that supports genetic connectivity and contributes to the capacity of species to adapt to environmental change. Yet, for most species, little is known about the specific environmental factors that influence genetic connectivity, or their effects on genetic diversity and differentiation. We used a landscape genetic approach to understand how geography and climate influence genetic connectivity in a foundation riparian tree (Populus angustifolia), and their relationships with specieswide patterns of genetic diversity and differentiation. Using multivariate restricted optimization in a reciprocal causal modelling framework, we quantified the relative contributions of riparian network connectivity, terrestrial upland resistance and climate gradients on genetic connectivity. We found that (i) all riparian corridors, regardless of river order, equally facilitated connectivity, while terrestrial uplands provided 2.5× more resistance to gene flow than riparian corridors. (ii) Cumulative differences in precipitation seasonality and precipitation of the warmest quarter were the primary climatic factors driving genetic differentiation; furthermore, maximum climate resistance was 45× greater than riparian resistance. (iii) Genetic diversity was positively correlated with connectivity (R² = 0.3744, p = .0019), illustrating the utility of resistance models for identifying landscape conditions that can support a species' ability to adapt to environmental change. From these results, we present a map highlighting key genetic connectivity corridors across P. angustifolia's range that if disrupted could have long‐term ecological and evolutionary consequences. Our findings provide recommendations for conservation and restoration management of threatened riparian ecosystems throughout the western USA and the high biodiversity they support.     

Detection of environmental and morphological adaptation despite high landscape genetic connectivity in a pest grasshopper (Phaulacridium vittatum)

Widespread species that exhibit both high gene flow and the capacity to occupy heterogeneous environments make excellent models for examining local selection processes along environmental gradients. Here we evaluate the influence of temperature and landscape variables on genetic connectivity and signatures of local adaptation in Phaulacridium vittatum, a widespread agricultural pest grasshopper, endemic to Australia. With sampling across a 900‐km latitudinal gradient, we genotyped 185 P. vittatum from 19 sites at 11,408 single nucleotide polymorphisms (SNPs) using ddRAD sequencing. Despite high gene flow across sites (pairwise F_ST = 0.0003–0.08), landscape genetic resistance modelling identified a positive nonlinear effect of mean annual temperature on genetic connectivity. Urban areas and water bodies had a greater influence on genetic distance among sites than pasture, agricultural areas and forest. Together, F_ST outlier tests and environmental association analysis (EAA) detected 242 unique SNPs under putative selection, with the highest numbers associated with latitude, mean annual temperature and body size. A combination of landscape genetic connectivity analysis together with EAA identified mean annual temperature as a key driver of both neutral gene flow and environmental selection processes. Gene annotation of putatively adaptive SNPs matched with gene functions for olfaction, metabolic detoxification and ultraviolet light shielding. Our results imply that this widespread agricultural pest has the potential to spread and adapt under shifting temperature regimes and land cover change.     

Genome-wide SNPs detect fine-scale genetic structure in threatened populations of squirrel glider Petaurus norfolcensis

Australian arboreal mammals are experiencing significant population declines, particularly due to land clearing and resulting habitat fragmentation. The squirrel glider, Petaurus norfolcensis, is a threatened species in New South Wales, with a stronghold population in the Lake Macquarie Local Government Area (LGA) where fragmentation due to urbanization is an ongoing problem for the species conservation. Here we report on the use of squirrel glider mitochondrial (385 bp cytochrome b gene, 70 individuals) and nuclear DNA (6,834 SNPs, 87 individuals) markers to assess their population genetic structure and connectivity across 14 locations sampled in the Lake Macquarie LGA. The mitochondrial DNA sequences detected evidence of a historical genetic bottleneck, while the genome-wide SNPs detected significant population structure in the Lake Macquarie squirrel glider populations at scales as fine as one kilometer. There was no evidence of inbreeding within patches, however there were clear effects of habitat fragmentation and biogeographical barriers on gene flow. A least cost path analysis identified thin linear corridors that have high priority for conservation. These areas should be protected to avoid further isolation of squirrel glider populations and the loss of genetic diversity through genetic drift.

     

Integrative approaches to guide conservation decisions: Using genomics to define conservation units and functional corridors

Climate change and increasing habitat loss greatly impact species survival, requiring range shifts, phenotypic plasticity and/or evolutionary change for long‐term persistence, which may not readily occur unaided in threatened species. Therefore, defining conservation actions requires a detailed assessment of evolutionary factors. Existing genetic diversity needs to be thoroughly evaluated and spatially mapped to define conservation units (CUs) in an evolutionary context, and we address that here. We also propose a multidisciplinary approach to determine corridors and functional connectivity between CUs by including genetic diversity in the modelling while controlling for isolation by distance and phylogeographic history. We evaluate our approach on a Near Threatened Iberian endemic rodent by analysing genotyping‐by‐sequencing (GBS) genomic data from 107 Cabrera voles (Microtus cabrerae), screening the entire species distribution to define categories of CUs and their connectivity: We defined six management units (MUs) which can be grouped into four evolutionarily significant units (ESUs) and three (putatively) adaptive units (AUs). We demonstrate that the three different categories of CU can be objectively defined using genomic data, and their characteristics and connectivity can inform conservation decision‐making. In particular, we show that connectivity of the Cabrera vole is very limited in eastern Iberia and that the pre‐Pyrenean and part of the Betic geographic nuclei contribute the most to the species genetic diversity. We argue that a multidisciplinary framework for CU definition is essential and that this framework needs a strong evolutionary basis.     

How do low dispersal species establish large range sizes? The case of the water beetle Graphoderus bilineatus

Species’ dispersal abilities have been considered a major driving force in establishment and maintenance of large range sizes. However, recent studies question the general validity of this relationship because the relationship between dispersal ability and range size might in some cases be less important than species phylogeny or local spatial attributes. In this study we used the water beetle Graphoderus bilineatus a philopatric species of conservation concern in Europe as a model to explain large range size and to support effective conservation measures for such species that also have limited dispersal. We recorded the presence/absence of G. bilineatus and measured 14 habitat and 20 landscape variables at 228 localities in Estonia, Poland and Sweden within the core range of the species. Using information theory and average multivariate logistic regression models we determined that presence of G. bilineatus depended on landscape connectivity, distance to a possible source habitat, and stability of the site; however, specificity of habitat characteristics was not vital for the species. We reason that the large range of G. bilineatus is best explained by the historical combination of lakes, river systems and wetlands which used to be highly connected throughout the central plains of Europe. Our data suggest that a broad habitat niche can prevent landscape elements from becoming barriers for species like G. bilineatus. Therefore, we question the usefulness of site protection as conservation measures for G. bilineatus and similar philopatric species. Instead, conservation actions should be focused at the landscape level to ensure a long‐term viability of such species across their range.     

Population assignment reveals low migratory connectivity in a weakly structured songbird

Understanding migratory connectivity is essential for determining the drivers behind population dynamics and for implementing effective conservation strategies for migratory species. Genetic markers provide a means to describe migratory connectivity; however, they can be uninformative for species with weak population genetic structure, which has limited their application. Here, we demonstrated a genomic approach to describing migratory connectivity in the prothonotary warbler, Protonotaria citrea, a Neotropical songbird of conservation concern. Using 26,189 single nucleotide polymorphisms (SNPs), we revealed regional genetic structure between the Mississippi River Valley and the Atlantic Seaboard with overall weak genetic differentiation among populations (F_ST = 0.0055; 95% CI: 0.0051–0.0059). Genetic variation had a stronger association with geographic rather than environmental factors, with each explaining 14.5% and 8.2% of genetic variation, respectively. By varying the numbers of genomic markers used in population assignment models with individuals of known provenance, we identified a maximum assignment accuracy (89.7% to site, 94.3% to region) using a subset of 600 highly differentiated SNPs. We then assigned samples from nonbreeding sites to breeding region and found low migratory connectivity. Our results highlight the importance of filtering markers for informative loci in models of population assignment. Quantifying migratory connectivity for weakly structured species will be useful for expanding studies to a wider range of migratory species across taxonomic groups and may contribute to a deeper understanding of the evolution of migratory strategies.     

Conservation genomics of anadromous Atlantic salmon across its North American range: outlier loci identify the same patterns of population structure as neutral loci

Anadromous Atlantic salmon (Salmo salar) is a species of major conservation and management concern in North America, where population abundance has been declining over the past 30 years. Effective conservation actions require the delineation of conservation units to appropriately reflect the spatial scale of intraspecific variation and local adaptation. Towards this goal, we used the most comprehensive genetic and genomic database for Atlantic salmon to date, covering the entire North American range of the species. The database included microsatellite data from 9142 individuals from 149 sampling locations and data from a medium‐density SNP array providing genotypes for >3000 SNPs for 50 sampling locations. We used neutral and putatively selected loci to integrate adaptive information in the definition of conservation units. Bayesian clustering with the microsatellite data set and with neutral SNPs identified regional groupings largely consistent with previously published regional assessments. The use of outlier SNPs did not result in major differences in the regional groupings, suggesting that neutral markers can reflect the geographic scale of local adaptation despite not being under selection. We also performed assignment tests to compare power obtained from microsatellites, neutral SNPs and outlier SNPs. Using SNP data substantially improved power compared to microsatellites, and an assignment success of 97% to the population of origin and of 100% to the region of origin was achieved when all SNP loci were used. Using outlier SNPs only resulted in minor improvements to assignment success to the population of origin but improved regional assignment. We discuss the implications of these new genetic resources for the conservation and management of Atlantic salmon in North America.     

Conservation of genetic diversity hotspots of the high‐valued relic yellowhorn (Xanthoceras sorbifolium) considering climate change predictions

Genetic structure and major climate factors may contribute to the distribution of genetic diversity of a highly valued oil tree species Xanthoceras sorbifolium (yellowhorn). Long‐term over utilization along with climate change is affecting the viability of yellowhorn wild populations. To preserve the species known and unknown valuable gene pools, the identification of genetic diversity “hotspots” is a prerequisite for their consideration as in situ conservation high priority. Chloroplast DNA (cpDNA) diversity was high among 38 natural populations (H_d = 0.717, K = 4.616, Tajmas’ D = ?0.22) and characterized by high genetic divergence (F_ST = 0.765) and relatively low gene flow (N_m = 0.03), indicating populations isolation reflecting the species’ habitat fragmentation and inbreeding depression. Six out of the studied 38 populations are defined as genetic diversity “hotspots.” The number and geographic direction of cpDNA mutation steps supported the species southwest to northeast migration history. Climatic factors such as extreme minimum temperature over 30 years indicated that the identified genetic “hotspots” are expected to experience 5°C temperature increase in next following 50 years. The results identified vulnerable genetic diversity “hotspots” and provided fundamental information for the species’ future conservation and breeding activities under the anticipated climate change. More specifically, the role of breeding as a component of a gene resource management strategy aimed at fulfilling both utilization and conservation goals.     

Seascape genomics reveals adaptive divergence in a connected and commercially important mollusc,the greenlip abalone (Haliotis laevigata), along a longitudinal environmental gradient

Populations of broadcast spawning marine organisms often have large sizes and are exposed to reduced genetic drift. Under such scenarios, strong selection associated with spatial environmental heterogeneity is expected to drive localized adaptive divergence, even in the face of connectivity. We tested this hypothesis using a seascape genomics approach in the commercially important greenlip abalone (Haliotis laevigata). We assessed how its population structure has been influenced by environmental heterogeneity along a zonal coastal boundary in southern Australia linked by strong oceanographic connectivity. Our data sets include 9,109 filtered SNPs for 371 abalones from 13 localities and environmental mapping across ~800 km. Genotype–environment association analyses and outlier tests defined 8,786 putatively neutral and 323 candidate adaptive loci. From a neutral perspective, the species is better represented by a metapopulation with very low differentiation (global F_ST = 0.0081) and weak isolation by distance following a stepping‐stone model. For the candidate adaptive loci, however, model‐based and model‐free approaches indicated five divergent population clusters. After controlling for spatial distance, the distribution of putatively adaptive variation was strongly correlated to selection linked to minimum sea surface temperature and oxygen concentration. Around 80 candidates were annotated to genes with functions related to high temperature and/or low oxygen tolerance, including genes that influence the resilience of abalone species found in other biogeographic regions. Our study includes a documented example about the uptake of genomic information in fisheries management and supports the hypothesis of adaptive divergence due to coastal environmental heterogeneity in a connected metapopulation of a broadcast spawner.     

Subspecies delineation amid phenotypic,geographic and genetic discordance in a songbird

Understanding the processes that drive divergence within and among species is a long‐standing goal in evolutionary biology. Traditional approaches to assessing differentiation rely on phenotypes to identify intra‐ and interspecific variation, but many species express subtle morphological gradients in which boundaries among forms are unclear. This intraspecific variation may be driven by differential adaptation to local conditions and may thereby reflect the evolutionary potential within a species. Here, we combine genetic and morphological data to evaluate intraspecific variation within the Nelson's (Ammodramus nelsoni) and salt marsh (Ammodramus caudacutus) sparrow complex, a group with populations that span considerable geographic distributions and a habitat gradient. We evaluated genetic structure among and within five putative subspecies of A. nelsoni and A. caudacutus using a reduced‐representation sequencing approach to generate a panel of 1929 SNPs among 69 individuals. Although we detected morphological differences among some groups, individuals sorted along a continuous phenotypic gradient. In contrast, the genetic data identified three distinct clusters corresponding to populations that inhabit coastal salt marsh, interior freshwater marsh and coastal brackish–water marsh habitats. These patterns support the current species‐level recognition but do not match the subspecies‐level taxonomy within each species—a finding which may have important conservation implications. We identified loci exhibiting patterns of elevated divergence among and within these species, indicating a role for local selective pressures in driving patterns of differentiation across the complex. We conclude that this evidence for adaptive variation among subspecies warrants the consideration of evolutionary potential and genetic novelty when identifying conservation units for this group.     

Extreme genetic similarity does not predict non‐breeding distribution of two closely related warblers

Detailed knowledge of migratory connectivity can facilitate effective conservation of Neotropical migrants by helping biologists understand where and when populations may be most limited. We studied the migratory behavior and non‐breeding distribution of two closely related species of conservation concern, the Golden‐winged Warbler (Vermivora chrysoptera) and Blue‐winged Warbler (Vermivora cyanoptera). Although both species have undergone dynamic range shifts and population changes attributed to habitat loss and social interactions promoting competition and hybridization, full life‐cycle conservation planning has been limited by a lack of information about their non‐breeding ecology. Because recent work has demonstrated that the two species are nearly identical genetically, we predicted that individuals from a single breeding population would have similar migratory timing and overwintering locations. In 2015, we placed light‐level geolocators on 25 males of both species and hybrids in an area of breeding sympatry at the Fort Drum Military Installation in Jefferson and Lewis counties, New York. Despite extreme genetic similarity, non‐breeding locations and duration of migration differed among genotypes. Golden‐winged Warblers (N = 2) overwintered > 1900 km southeast of the nearest Blue‐winged Warbler (N = 3) and spent nearly twice as many days in migration; hybrids (N = 2) had intermediate wintering distributions and migratory timing. Spring migration departure dates were staggered based on distance from the breeding area, and all birds arrived at the breeding site within 8 days of each other. Our results show that Golden‐winged Warblers and Blue‐winged Warblers in our study area retain species‐specific non‐breeding locations despite extreme genetic similarity, and suggest that non‐breeding locations and migratory timing vary along a genetic gradient. If the migratory period is limiting for these species, our results also suggest that Golden‐winged Warblers in our study population may be more vulnerable to population decline than Blue‐winged Warblers because they spend almost twice as many days migrating.     

Exploring genomic variation associated with drought stress in Picea mariana populations

Predicted increases in drought and heat stress will likely induce shifts in species bioclimatic envelopes. Genetic variants adapted to water limitation may prove pivotal for species response under scenarios of increasing drought. In this study, we aimed to explore this hypothesis by investigating genetic variation in 16 populations of black spruce (Picea mariana) in relation to climate variables in Alaska. A total of 520 single nucleotide polymorphisms (SNPs) were genotyped for 158 trees sampled from areas of contrasting climate regimes. We used multivariate and univariate genotype‐by‐environment approaches along with available gene annotations to investigate the relationship between climate and genetic variation among sampled populations. Nine SNPs were identified as having a significant association with climate, of which five were related to drought stress response. Outlier SNPs with respect to the overall environment were significantly overrepresented for several biological functions relevant for coping with variable hydric regimes, including osmotic stress response. This genomic imprint is consistent with local adaptation of black spruce to drought stress. These results suggest that natural selection acting on standing variation prompts local adaptation in forest stands facing water limitation. Improved understanding of possible adaptive responses could inform our projections about future forest dynamics and help prioritize populations that harbor valuable genetic diversity for conservation.     

Identifying climate refugia and its potential impact on Tibetan brown bear (Ursus arctos pruinosus) in Sanjiangyuan National Park,China

Climate change has direct impacts on wildlife and future biodiversity protection efforts. Vulnerability assessment and habitat connectivity analyses are necessary for drafting effective conservation strategies for threatened species such as the Tibetan brown bear (Ursus arctos pruinosus). We used the maximum entropy (MaxEnt) model to assess the current (1950–2000) and future (2041–2060) habitat suitability by combining bioclimatic and environmental variables, and identified potential climate refugia for Tibetan brown bears in Sanjiangyuan National Park, China. Next, we selected Circuit model to simulate potential migration paths based on current and future climatically suitable habitat. Results indicate a total area of potential suitable habitat under the current climate scenario of approximately 31,649.46 km², of which 28,778.29 km² would be unsuitable by the 2050s. Potentially suitable habitat under the future climate scenario was projected to cover an area of 23,738.6 km². Climate refugia occupied 2,871.17 km², primarily in the midwestern and northeastern regions of Yangtze River Zone, as well as the northern region of Yellow River Zone. The altitude of climate refugia ranged from 4,307 to 5,524 m, with 52.93% lying at altitudes between 4,300 and 4,600 m. Refugia were mainly distributed on bare rock, alpine steppe, and alpine meadow. Corridors linking areas of potentially suitable brown bear habitat and a substantial portion of paths with low‐resistance value were distributed in climate refugia. We recommend various actions to ameliorate the impact of climate change on brown bears, such as protecting climatically suitable habitat, establishing habitat corridors, restructuring conservation areas, and strengthening monitoring efforts.     

Spatial population genomics of the brown rat (Rattus norvegicus) in New York City

Human commensal species such as rodent pests are often widely distributed across cities and threaten both infrastructure and public health. Spatially explicit population genomic methods provide insights into movements for cryptic pests that drive evolutionary connectivity across multiple spatial scales. We examined spatial patterns of neutral genomewide variation in brown rats (Rattus norvegicus) across Manhattan, New York City (NYC), using 262 samples and 61,401 SNPs to understand (i) relatedness among nearby individuals and the extent of spatial genetic structure in a discrete urban landscape; (ii) the geographic origin of NYC rats, using a large, previously published data set of global rat genotypes; and (iii) heterogeneity in gene flow across the city, particularly deviations from isolation by distance. We found that rats separated by ≤200 m exhibit strong spatial autocorrelation (r = .3, p = .001) and the effects of localized genetic drift extend to a range of 1,400 m. Across Manhattan, rats exhibited a homogeneous population origin from rats that likely invaded from Great Britain. While traditional approaches identified a single evolutionary cluster with clinal structure across Manhattan, recently developed methods (e.g., fineSTRUCTURE, sPCA, EEMS) provided evidence of reduced dispersal across the island's less residential Midtown region resulting in fine‐scale genetic structuring (F_ST = 0.01) and two evolutionary clusters (Uptown and Downtown Manhattan). Thus, while some urban populations of human commensals may appear to be continuously distributed, landscape heterogeneity within cities can drive differences in habitat quality and dispersal, with implications for the spatial distribution of genomic variation, population management and the study of widely distributed pests.     

Habitat predictors of genetic diversity for two sympatric wetland‐breeding amphibian species

Population genetic diversity is widely accepted as important to the conservation and management of wildlife. However, habitat features may differentially affect evolutionary processes that facilitate population genetic diversity among sympatric species. We measured genetic diversity for two pond‐breeding amphibian species (Dwarf salamanders, Eurycea quadridigitata; and Southern Leopard frogs, Lithobates sphenocephalus) to understand how habitat characteristics and spatial scale affect genetic diversity across a landscape. Samples were collected from wetlands on a longleaf pine reserve in Georgia. We genotyped microsatellite loci for both species to assess population structures and determine which habitat features were most closely associated with observed heterozygosity and rarefied allelic richness. Both species exhibited significant population genetic structure; however, structure in Southern Leopard frogs was driven primarily by one outlier site. Dwarf salamander allelic richness was greater at sites with less surrounding road area within 0.5 km and more wetland area within 1.0 and 2.5 km, and heterozygosity was greater at sites with more wetland area within 0.5 km. In contrast, neither measure of Southern Leopard frog genetic diversity was associated with any habitat features at any scale we evaluated. Genetic diversity in the Dwarf salamander was strongly associated with land cover variables up to 2.5 km away from breeding wetlands, and/or results suggest that minimizing roads in wetland buffers may be beneficial to the maintenance of population genetic diversity. This study suggests that patterns of genetic differentiation and genetic diversity have associations with different habitat features across different spatial scales for two syntopic pond‐breeding amphibian species.     

The challenge of separating signatures of local adaptation from those of isolation by distance and colonization history: The case of two white pines

Accurately detecting signatures of local adaptation using genetic‐environment associations (GEAs) requires controlling for neutral patterns of population structure to reduce the risk of false positives. However, a high degree of collinearity between climatic gradients and neutral population structure can greatly reduce power, and the performance of GEA methods in such case is rarely evaluated in empirical studies. In this study, we attempted to disentangle the effects of local adaptation and isolation by environment (IBE) from those of isolation by distance (IBD) and isolation by colonization from glacial refugia (IBC) using range‐wide samples in two white pine species. For this, SNPs from 168 genes, including 52 candidate genes for growth and phenology, were genotyped in 133 and 61 populations of Pinus strobus and P. monticola, respectively. For P. strobus and using all 153 SNPs, climate (IBE) did not significantly explained among‐population variation when controlling for IBD and IBC in redundancy analyses (RDAs). However, 26 SNPs were significantly associated with climate in single‐locus GEA analyses (Bayenv2 and LFMM), suggesting that local adaptation took place in the presence of high gene flow. For P. monticola, we found no evidence of IBE using RDAs and weaker signatures of local adaptation using GEA and F_ST outlier tests, consistent with adaptation via phenotypic plasticity. In both species, the majority of the explained among‐population variation (69 to 96%) could not be partitioned between the effects of IBE, IBD, and IBC. GEA methods can account differently for this confounded variation, and this could explain the small overlap of SNPs detected between Bayenv2 and LFMM. Our study illustrates the inherent difficulty of taking into account neutral structure in natural populations and the importance of sampling designs that maximize climatic variation, while minimizing collinearity between climatic gradients and neutral structure.     

Do plant populations on distinct inselbergs talk to each other? A case study of genetic connectivity of a bromeliad species in an Ocbil landscape

Here, we explore the historical and contemporaneous patterns of connectivity among Encholirium horridum populations located on granitic inselbergs in an Ocbil landscape within the Brazilian Atlantic Forest, using both nuclear and chloroplast microsatellite markers. Beyond to assess the E. horridum population genetic structure, we built species distribution models across four periods (current conditions, mid‐Holocene, Last Glacial Maximum [LGM], and Last Interglacial) and inferred putative dispersal corridors using a least‐cost path analysis to elucidate biogeographic patterns. Overall, high and significant genetic divergence was estimated among populations for both nuclear and plastid DNA (Φ_ST(n) = 0.463 and Φ_ST(plastid) = 0.961, respectively, p < .001). For nuclear genome, almost total absence of genetic admixture among populations and very low migration rates were evident, corroborating with the very low estimates of immigration and emigration rates observed among E. horridum populations. Based on the cpDNA results, putative dispersal routes in Sugar Loaf Land across cycles of climatic fluctuations in the Quaternary period revealed that the populations’ connectivity changed little during those events. Genetic analyses highlighted the low genetic connectivity and long‐term persistence of populations, and the founder effect and genetic drift seemed to have been very important processes that shaped the current diversity and genetic structure observed in both genomes. The genetic singularity of each population clearly shows the need for in situ conservation of all of them.     

Single nucleotide polymorphisms unravel hierarchical divergence and signatures of selection among Alaskan sockeye salmon (Oncorhynchus nerka) populations

Background  

Disentangling the roles of geography and ecology driving population divergence and distinguishing adaptive from neutral evolution at the molecular level have been common goals among evolutionary and conservation biologists. Using single nucleotide polymorphism (SNP) multilocus genotypes for 31 sockeye salmon (Oncorhynchus nerka) populations from the Kvichak River, Alaska, we assessed the relative roles of geography (discrete boundaries or continuous distance) and ecology (spawning habitat and timing) driving genetic divergence in this species at varying spatial scales within the drainage. We also evaluated two outlier detection methods to characterize candidate SNPs responding to environmental selection, emphasizing which mechanism(s) may maintain the genetic variation of outlier loci.     

Phylogeography and population genetics of pine butterflies: Sky islands increase genetic divergence

The sky islands of southeastern Arizona (AZ) mark a major transition zone between tropical and temperate biota and are considered a neglected biodiversity hotspot. Dispersal ability and host plant specificity are thought to impact the history and diversity of insect populations across the sky islands. We aimed to investigate the population structure and phylogeography of two pine‐feeding pierid butterflies, the pine white (Neophasia menapia) and the Mexican pine white (Neophasia terlooii), restricted to these “islands” at this transition zone. Given their dependence on pines as the larval hosts, we hypothesized that habitat connectivity affects population structure and is at least in part responsible for their allopatry. We sampled DNA from freshly collected butterflies from 17 sites in the sky islands and adjacent high‐elevation habitats and sequenced these samples using ddRADSeq. Up to 15,399 SNPs were discovered and analyzed in population genetic and phylogenetic contexts with Stacks and pyRAD pipelines. Low genetic differentiation in N. menapia suggests that it is panmictic. Conversely, there is strong evidence for population structure within N. terlooii. Each sky island likely contains a population of N. terlooii, and clustering is hierarchical, with populations on proximal mountains being more related to each other. The N. menapia habitat, which is largely contiguous, facilitates panmixia, while the N. terlooii habitat, restricted to the higher elevations on each sky island, creates distinct population structure. Phylogenetic results corroborate those from population genetic analyses. The historical climate‐driven fluxes in forest habitat connectivity have implications for understanding the biodiversity of fragmented habitats.