Bat species assemblage predicts coronavirus prevalence.

Anthropogenic disturbances and the subsequent loss of biodiversity are altering species abundances and communities. Since species vary in their pathogen competence, spatio-temporal changes in host assemblages may lead to changes in disease dynamics. We explore how longitudinal changes in bat species assemblages affect the disease dynamics of coronaviruses (CoVs) in more than 2300 cave-dwelling bats captured over two years from five caves in Ghana. This reveals uneven CoV infection patterns between closely related species, with the alpha-CoV 229E-like and SARS-related beta-CoV 2b emerging as multi-host pathogens. Prevalence and infection likelihood for both phylogenetically distinct CoVs is influenced by the abundance of competent species and naïve subadults. Broadly, bat species vary in CoV competence, and highly competent species are more common in less diverse communities, leading to increased CoV prevalence in less diverse bat assemblages. In line with the One Health framework, our work supports the notion that biodiversity conservation may be the most proactive measure to prevent the spread of pathogens with zoonotic potential.

Environmental DNA highlights the influence of salinity and agricultural run-off on coastal fish assemblages in the Great Barrier Reef region.

Agricultural run-off in Australia's Mackay-Whitsunday region is a major source of nutrient and pesticide pollution to coastal and inshore ecosystems of the Great Barrier Reef. While the effects of run-off are well documented for the region's coral and seagrass habitats, the ecological impacts on estuaries, the direct recipients of run-off, are less known. This is particularly true for fish communities, which are shaped by the physico-chemical properties of coastal waterways that vary greatly in tropical regions. To address this knowledge gap, we used environmental DNA (eDNA) metabarcoding to examine fish assemblages at four locations (three estuaries and a harbour) subjected to varying levels of agricultural run-off during a wet and dry season. Pesticide and nutrient concentrations were markedly elevated during the sampled wet season with the influx of freshwater and agricultural run-off. Fish taxa richness significantly decreased in all three estuaries (F = 164.73, P = <0.001), along with pronounced changes in community composition (F = 46.68, P = 0.001) associated with environmental variables (largely salinity: 27.48% contribution to total variance). In contrast, the nearby Mackay Harbour exhibited a far more stable community structure, with no marked changes in fish assemblages observed between the sampled seasons. Among the four sampled locations, variation in fish community composition was more pronounced within the wet season (F = 2.5, P = 0.001). Notably, variation in the wet season was significantly correlated with agricultural contaminants (phosphorus: 6.25%, pesticides: 5.22%) alongside environmental variables (salinity: 5.61%, DOC: 5.57%). Historically contaminated and relatively unimpacted estuaries each demonstrated distinct fish communities, reflecting their associated catchment use. Our findings emphasise that while seasonal effects play a key role in shaping the community structure of fish in this region, agricultural contaminants are also important contributors in estuarine systems.

Sperm Storage in a Family-Living Lizard, the Tree Skink (Egernia striolata).

The ability to produce viable offspring without recently mating, either through sperm storage or parthenogenesis, can provide fitness advantages under a suite of challenging ecological scenarios. Using genetic analysis, we demonstrate that 3 wild-caught female Tree Skinks (Egernia striolata) reproduced in captivity with no access to males for over a year, and that this is best explained by sperm storage. To the best of our knowledge, this is the first time female sperm storage has been documented in any monogamous family-living reptile, including social Australian egerniine skinks (from the subfamily Egerniinae). Furthermore, by using paternal reconstruction of genotypes we show that captive-born offspring produced by the same females in the preceding year, presumably without sperm storage, were sired by different males. We qualitatively compared aspects of these females' mates and offspring between years. The parents of each litter were unrelated, but paternal and offspring genotypes from litters resulting from stored sperm were more heterozygous than those inferred to be from recent matings. Family-living egerniine skinks generally have low rates of multiple paternity, yet our study suggests that female sperm storage, potentially from outside social partners, offers the real possibility of benefits. Possible benefits include increasing genetic compatibility of mates and avoiding inbreeding depression via cryptic female choice. Sperm storage in Tree Skinks, a family-living lizard with a monogamous mating system, suggests that females may bet-hedge through extra-pair copulation with more heterozygous males, reinforcing the idea that females could have more control on reproductive outcomes than previously thought.

Microgeographical adaptation corresponds to elevational distributions of congeneric montane grasshoppers.

Local adaptation can occur at small spatial scales relative to the dispersal capacity of species. Alpine ecosystems have sharp environmental clines that offer an opportunity to investigate the effects of fine-scale shifts in species' niche breadth on adaptive genetic processes. Here we examine two grasshopper species endemic to the Australian Alps (Kosciuscola spp.) that differ in elevational niche breadth: one broader, K. usitatus (1400-2200 m), and one narrower, K. tristis (1600-2000 m). We examine signatures of selection with respect to environmental and morphological variables in two mountain regions using FST outlier tests and environmental association analyses (EAAs) applied to single nucleotide polymorphism (SNP) data (K. usitatus: 9017 SNPs, n = 130; K. tristis: 7363 SNPs, n = 135). Stronger genetic structure was found in the more narrowly distributed K. tristis, which showed almost twice the number of SNPs under putative selection (10.8%) compared with K. usitatus (5.3%). When examining SNPs in common across species (n = 3058), 260 SNPs (8.5%) were outliers shared across species, and these were mostly associated with elevation, a proxy for temperature, suggesting parallel adaptive processes in response to climatic drivers. Additive polygenic scores (an estimate of the cumulative signal of selection across all candidate loci) were nonlinearly and positively correlated with elevation in both species. However, a steeper correlation in K. tristis indicated a stronger signal of spatially varying selection towards higher elevations. Our study illustrates that the niche breadth of co-occurring and related species distributed along the same environmental cline is associated with differences in patterns of microgeographical adaptation.

The Influence of Environmental Variation on the Genetic Structure of a Poison Frog Distributed Across Continuous Amazonian Rainforest.

Biogeographic barriers such as rivers have been shown to shape spatial patterns of biodiversity in the Amazon basin, yet relatively little is known about the distribution of genetic variation across continuous rainforest. Here, we characterize the genetic structure of the brilliant-thighed poison frog (Allobates femoralis) across an 880-km-long transect along the Purus-Madeira interfluve south of the Amazon river, based on 64 individuals genotyped at 7609 single-nucleotide polymorphism (SNP) loci. A population tree and clustering analyses revealed 4 distinct genetic groups, one of which was strongly divergent. These genetic groups were concomitant with femoral spot coloration differences, which was intermediate within a zone of admixture between two of the groups. The location of these genetic groups did not consistently correspond to current ecological transitions between major forest types. A multimodel approach to quantify the relative influence of isolation-by-geographic distance (IBD) and isolation-by-environmental resistance (IBR) nevertheless revealed that, in addition to a strong signal of IBD, spatial genetic differentiation was explained by IBR primarily linked to dry season intensity (r2 = 8.4%) and canopy cover (r2 = 6.4%). We show significant phylogenetic divergence in the absence of obvious biogeographical barriers and that finer-scaled measures of genetic structure are associated with environmental variables also known to predict the density of A. femoralis.

Conservation prioritization can resolve the flagship species conundrum.

Conservation strategies based on charismatic flagship species, such as tigers, lions, and elephants, successfully attract funding from individuals and corporate donors. However, critics of this species-focused approach argue it wastes resources and often does not benefit broader biodiversity. If true, then the best way of raising conservation funds excludes the best way of spending it. Here we show that this conundrum can be resolved, and that the flagship species approach does not impede cost-effective conservation. Through a tailored prioritization approach, we identify places containing flagship species while also maximizing global biodiversity representation (based on 19,616 terrestrial and freshwater species). We then compare these results to scenarios that only maximized biodiversity representation, and demonstrate that our flagship-based approach achieves 79-89% of our objective. This provides strong evidence that prudently selected flagships can both raise funds for conservation and help target where these resources are best spent to conserve biodiversity.

Social Barriers in Ecological Landscapes: The Social Resistance Hypothesis.

Across animal societies, individuals invest time and energy in social interactions. The social landscape that emerges from these interactions can then generate barriers that limit the ability of individuals to disperse to, and reproduce in, groups or populations. Therefore, social barriers can contribute to the difference between the physical capacity for movement through the habitat and subsequent gene flow. We call this contributing effect 'social resistance'. We propose that social resistance can act as an agent of selection on key life-history strategies and promote the evolution of social strategies that facilitate effective dispersal. By linking landscape genetics and social behaviour, the social resistance hypothesis generates predictions integrating dispersal, connectivity, and life-history evolution.

The evolution of polymorphism in the warning coloration of the Amazonian poison frog Adelphobates galactonotus.

While intraspecific variation in aposematic signals can be selected for by different predatory responses, their evolution is also contingent on other processes shaping genetic variation. We evaluate the relative contributions of selection, geographic isolation, and random genetic drift to the evolution of aposematic color polymorphism in the poison frog Adelphobates galactonotus, distributed throughout eastern Brazilian Amazonia. Dorsal coloration was measured for 111 individuals and genetic data were obtained from 220 individuals at two mitochondrial genes (mtDNA) and 7963 Single Nucleotide Polymorphisms (SNPs). Four color categories were described (brown, blue, yellow, orange) and our models of frog and bird visual systems indicated that each color was distinguishable for these taxa. Using outlier and correlative analyses we found no compelling genetic evidence for color being under divergent selection. A time-calibrated mtDNA tree suggests that the present distribution of dorsal coloration resulted from processes occurring during the Pleistocene. Separate phylogenies based on SNPs and mtDNA resolved the same well supported clades, each containing different colored populations. Ancestral character state analysis provided some evidence for evolutionary transitions in color type. Genetic structure was more strongly associated with geographic features, than color category, suggesting that the distribution of color is explained by localized processes. Evidence for geographic isolation together with estimates of low effective population size implicates drift as playing a key role in color diversification. Our results highlight the relevance of considering the neutral processes involved with the evolution of traits with important fitness consequences.

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 FST  = 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, FST 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.

Morphological Variation Tracks Environmental Gradients in an Agricultural Pest, Phaulacridium vittatum (Orthoptera: Acrididae).

Invertebrate pests often show high morphological variation and wide environmental tolerances. Knowledge of how phenotypic variation is associated with environmental heterogeneity can elucidate the processes underpinning these patterns. Here we examine morphological variation and relative abundance along environmental gradients in a widespread agricultural pest, native to Australia, the wingless grasshopper Phaulacridium vittatum (Sjöstedt). We test for correlations between body size, wing presence, and stripe polymorphism with environmental variables. Using multiple regression and mixed-effects modeling, body size and stripe polymorphism were positively associated with solar radiation, and wing presence was positively associated with foliage projective cover (FPC). There were no associations between body size or morphological traits with relative abundance. However, relative abundance was positively associated with latitude, soil moisture, and wind speed, but was negatively associated with FPC. Therefore, sites with low relative abundance and high forest cover were more likely to contain winged individuals. Overall, our results suggest that environmental and climatic conditions strongly influence the relative abundance and the distribution of morphotypes in P. vittatum, which is likely to affect dispersal and fitness in different landscapes. This knowledge is useful for informing how environmental change might influence the future spread and impact of this agricultural pest.

Different environmental gradients affect different measures of snake ß-diversity in the Amazon rainforests.

Mechanisms generating and maintaining biodiversity at regional scales may be evaluated by quantifying ß-diversity along environmental gradients. Differences in assemblages result in biotic complementarities and redundancies among sites, which may be quantified through multi-dimensional approaches incorporating taxonomic ß-diversity (TBD), functional ß-diversity (FBD) and phylogenetic ß-diversity (PBD). Here we test the hypothesis that snake TBD, FBD and PBD are influenced by environmental gradients, independently of geographic distance. The gradients tested are expected to affect snake assemblages indirectly, such as clay content in the soil determining primary production and height above the nearest drainage determining prey availability, or directly, such as percentage of tree cover determining availability of resting and nesting sites, and climate (temperature and precipitation) causing physiological filtering. We sampled snakes in 21 sampling plots, each covering five km2, distributed over 880 km in the central-southern Amazon Basin. We used dissimilarities between sampling sites to quantify TBD, FBD and PBD, which were response variables in multiple-linear-regression and redundancy analysis models. We show that patterns of snake community composition based on TBD, FBD and PBD are associated with environmental heterogeneity in the Amazon. Despite positive correlations between all ß-diversity measures, TBD responded to different environmental gradients compared to FBD and PBD. Our findings suggest that multi-dimensional approaches are more informative for ecological studies and conservation actions compared to a single diversity measure.

Soil and forest structure predicts large-scale patterns of occurrence and local abundance of a widespread Amazonian frog.

The distribution of biodiversity within the Amazon basin is often structured by sharp environmental boundaries, such as large rivers. The Amazon region is also characterized by subtle environmental clines, but how they might affect the distributions and abundance of organisms has so far received less attention. Here, we test whether soil and forest characteristics are associated with the occurrence and relative abundance of the forest-floor dwelling Aromobatid frog, Allobates femoralis. We applied a structured sampling regime along an 880 km long transect through forest of different density. High detection probabilities were estimated for A. femoralis in each of the sampling modules. Using generalized linear mixed-effects models and simple linear regressions that take detectability into account, we show that A. femoralis is more abundant in open forests than in dense forests. The presence and relative abundance of A. femoralis is also positively associated with clay-rich soils, which are poorly drained and therefore likely support the standing water bodies required for reproduction. Taken together, we demonstrate that relatively easy-to-measure environmental features can explain the distribution and abundance of a widespread species at different spatial scales. Such proxies are of clear value to ecologists and conservation managers working in large inaccessible areas such as the Amazon basin.

Monogamous sperm storage and permanent worker sterility in a long-lived ambrosia beetle.

The lifetime monogamy hypothesis claims that the evolution of permanently unmated worker castes always requires maximal full-sibling relatedness to be established first. The long-lived diploid ambrosia beetle Austroplatypus incompertus (Schedl) is known to be highly social, but whether it has lifetime sterile castes has remained unclear. Here we show that the gallery systems of this beetle inside the heartwood of live Eucalyptus trees are always inhabited by a single core family, consisting of a lifetime-inseminated mother, permanently unmated daughter workers, and immatures that are always full siblings to each other and their adult caretakers. Overall sex ratios are even. Males always disperse and only survive as stored sperm, but female offspring either disperse to mate and found their own colony or assume unmated worker roles, probably surviving for many years without any reproductive potential because tarsal loss precludes later dispersal. A well-supported Platypodinae phylogeny has allowed us to infer that parental monogamy evolved before a lifetime-unmated worker caste emerged, confirming the prediction that monogamy and full-sibling relatedness are necessary conditions for the evolution of such workers. The initially very challenging but ultimately long-term stable nesting habitat in live trees appears to have provided the crucial benefit/cost factor for maintaining selection for permanently sterile workers after strict monogamy and lifetime sperm storage had become established in this curculionid coleopteran lineage.

The search for loci under selection: trends, biases and progress.

Detecting genetic variants under selection using FST outlier analysis (OA) and environmental association analyses (EAAs) are popular approaches that provide insight into the genetic basis of local adaptation. Despite the frequent use of OA and EAA approaches and their increasing attractiveness for detecting signatures of selection, their application to field-based empirical data have not been synthesized. Here, we review 66 empirical studies that use Single Nucleotide Polymorphisms (SNPs) in OA and EAA. We report trends and biases across biological systems, sequencing methods, approaches, parameters, environmental variables and their influence on detecting signatures of selection. We found striking variability in both the use and reporting of environmental data and statistical parameters. For example, linkage disequilibrium among SNPs and numbers of unique SNP associations identified with EAA were rarely reported. The proportion of putatively adaptive SNPs detected varied widely among studies, and decreased with the number of SNPs analysed. We found that genomic sampling effort had a greater impact than biological sampling effort on the proportion of identified SNPs under selection. OA identified a higher proportion of outliers when more individuals were sampled, but this was not the case for EAA. To facilitate repeatability, interpretation and synthesis of studies detecting selection, we recommend that future studies consistently report geographical coordinates, environmental data, model parameters, linkage disequilibrium, and measures of genetic structure. Identifying standards for how OA and EAA studies are designed and reported will aid future transparency and comparability of SNP-based selection studies and help to progress landscape and evolutionary genomics.

Under the weather?-The direct effects of climate warming on a threatened desert lizard are mediated by their activity phase and burrow system.

For ectotherms such as lizards, the importance of behavioural thermoregulation in avoiding thermal extremes is well-established and is increasingly acknowledged in modern studies of climate warming and its impacts. Less appreciated and understood are the buffering roles of retreat sites and activity phase, in part because of logistical challenges of studying below-ground activity. Burrowing and nocturnal activity are key behavioural adaptations that have enabled a diverse range of reptiles to survive extreme environmental temperatures within hot desert regions. Yet, the direct impact of recent global warming on activity potential has been hypothesised to have caused extinctions in desert lizards, including the Australian arid zone skink Liopholis kintorei. We test the relevance of this hypothesis through a detailed characterisation of the above- and below-ground thermal and hydric microclimates available to, and used by, L. kintorei. We integrate operative temperatures with observed body temperatures to construct daily activity budgets, including the inference of subterranean behaviour. We then assess the likelihood that contemporary and future local extinctions in this species, and those of similar burrowing habits, could be explained by the direct effects of warming on its activity budget and exposure to thermal extremes. We found that L. kintorei spent only 4% of its time active on the surface, primarily at dusk, and that overall potential surface activity will be increased, not restricted, with climate warming. The burrow system provides an exceptional buffer to current and future maximum extremes of temperature (≈40°C reduction from potential surface temperatures), and desiccation (burrows near 100% humidity). Therefore, any climate warming impacts on this species are likely to be indirect. Our findings reflect the general buffering capacity of underground microclimates, therefore, our conclusions for L. kintorei are more generally applicable to nocturnal and crepuscular ectotherms, and highlight the need to consider the buffering properties of retreat sites and activity phase when forecasting climate change impacts.

Geographic vocal variation and perceptual discrimination abilities in male Australian sea lions.

Vocal characteristics can vary among and within populations. In species with geographic variation in the structure of vocalizations, individuals may have the ability to discriminate between calls from local and non-local individuals. The ability to distinguish differences in acoustic signals is likely to have a significant influence on the outcome of social interactions between individuals, including potentially mate selection and breeding success. Pinnipeds (seals, fur seals, sea lions and walruses) are highly vocal yet the Australian sea lion (Neophoca cinerea) is the only eared seal known to show geographic vocal variation in male barks. Barks are produced in many social interactions and encode sufficient information for both individual and colony identity to be discriminable. Yet until now, whether males could themselves discriminate these bark differences was unclear. We performed playback experiments in four breeding colonies to investigate whether males can discriminate local from non-local barks. Overall, males responded more strongly to barks from their own colony compared to barks from other colonies regardless of whether those other colonies were close or distant. Competition for females is high in Australian sea lions, but mating periods are asynchronous across colonies. The ability to correctly assess whether a male is from the same colony, thus representing a potential competitor for mates, or merely a visitor from elsewhere, may influence how males interact with others. Given the high cost of fighting, the ability to discern competitors may influence the nature of male-male interactions and ultimately influence how they allocate reproductive effort.

Influence of adaptive capacity on the outcome of climate change vulnerability assessment.

Climate change vulnerability assessment (CCVA) has become a mainstay conservation decision support tool. CCVAs are recommended to incorporate three elements of vulnerability - exposure, sensitivity and adaptive capacity - yet, lack of data frequently leads to the latter being excluded. Further, weighted or unweighted scoring schemes, based on expert opinion, may be applied. Comparisons of these approaches are rare. In a CCVA for 17 Australian lizard species, we show that membership within three vulnerability categories (low, medium and high) generally remained similar regardless of the framework or scoring scheme. There was one exception however, where, under the warm/dry scenario for 2070, including adaptive capacity lead to five fewer species being classified as highly vulnerable. Two species, Eulamprus leuraensis and E. kosciuskoi, were consistently ranked the most vulnerable, primarily due to projected losses in climatically suitable habitat, narrow thermal tolerance and specialist habitat requirements. Our findings provide relevant information for prioritizing target species for conservation and choosing appropriate conservation actions. We conclude that for the species included in this study, the framework and scoring scheme used had little impact on the identification of the most vulnerable species. We caution, however, that this outcome may not apply to other taxa or regions.

Combining dispersal, landscape connectivity and habitat suitability to assess climate-induced changes in the distribution of Cunningham's skink, Egernia cunninghami.

The ability of species to track their climate niche is dependent on their dispersal potential and the connectivity of the landscape matrix linking current and future suitable habitat. However, studies modeling climate-driven range shifts rarely address the movement of species across landscapes realistically, often assuming "unlimited" or "no" dispersal. Here, we incorporate dispersal rate and landscape connectivity with a species distribution model (Maxent) to assess the extent to which the Cunningham's skink (Egernia cunninghami) may be capable of tracking spatial shifts in suitable habitat as climate changes. Our model was projected onto four contrasting, but equally plausible, scenarios describing futures that are (relative to now) hot/wet, warm/dry, hot/with similar precipitation and warm/wet, at six time horizons with decadal intervals (2020-2070) and at two spatial resolutions: 1 km and 250 m. The size of suitable habitat was projected to decline 23-63% at 1 km and 26-64% at 250 m, by 2070. Combining Maxent output with the dispersal rate of the species and connectivity of the intervening landscape matrix showed that most current populations in regions projected to become unsuitable in the medium to long term, will be unable to shift the distance necessary to reach suitable habitat. In particular, numerous populations currently inhabiting the trailing edge of the species' range are highly unlikely to be able to disperse fast enough to track climate change. Unless these populations are capable of adaptation they are likely to be extirpated. We note, however, that the core of the species distribution remains suitable across the broad spectrum of climate scenarios considered. Our findings highlight challenges faced by philopatric species and the importance of adaptation for the persistence of peripheral populations under climate change.

Contrasting Patterns of Gene Flow for Amazonian Snakes That Actively Forage and Those That Wait in Ambush.

Knowledge of genetic structure, geographic distance and environmental heterogeneity can be used to identify environmental features and natural history traits that influence dispersal and gene flow. Foraging mode is a trait that might predict dispersal capacity in snakes, because actively foragers typically have greater movement rates than ambush predators. Here, we test the hypothesis that 2 actively foraging snakes have higher levels of gene flow than 2 ambush predators. We evaluated these 4 co-distributed species of snakes in the Brazilian Amazon. Snakes were sampled along an 880 km transect from the central to the southwest of the Amazon basin, which covered a mosaic of vegetation types and seasonal differences in climate. We analyzed thousands of single nucleotide polymorphisms to compare patterns of neutral gene flow based on isolation by geographic distance (IBD) and environmental resistance (IBR). We show that IBD and IBR were only evident in ambush predators, implying lower levels of dispersal than the active foragers. Therefore, gene flow was high enough in the active foragers analyzed here to prevent any build-up of spatial genotypic structure with respect to geographic distance and environmental heterogeneity.

Cunningham's skinks show low genetic connectivity and signatures of divergent selection across its distribution.

Establishing corridors of connecting habitat has become a mainstay conservation strategy to maintain gene flow and facilitate climate-driven range shifts. Yet, little attention has been given to ascertaining the extent to which corridors will benefit philopatric species, which might exhibit localized adaptation. Measures of genetic connectivity and adaptive genetic variation across species' ranges can help fill this knowledge gap. Here, we characterized the spatial genetic structure of Cunningham's skink (Egernia cunninghami), a philopatric species distributed along Australia's Great Dividing Range, and assessed evidence of localized adaptation. Analysis of 4,274 SNPs from 94 individuals sampled at four localities spanning 500 km and 4° of latitude revealed strong genetic structuring at neutral loci (mean FST ± SD = 0.603 ± 0.237) among the localities. Putatively neutral SNPs and those under divergent selection yielded contrasting spatial patterns, with the latter identifying two genetically distinct clusters. Given low genetic connectivity of the four localities, we suggest that the natural movement rate of this species is insufficient to keep pace with spatial shifts to its climate envelope, irrespective of habitat availability. In addition, our finding of localized adaptation highlights the risk of outbreeding depression should the translocation of individuals be adopted as a conservation management strategy.