Quantifying intraspecific variation in host resistance and tolerance to a lethal pathogen.

Testing for intraspecific variation for host tolerance or resistance in wild populations is important for informing conservation decisions about captive breeding, translocation, and disease treatment. Here, we test the importance of tolerance and resistance in multiple populations of boreal toads (Anaxyrus boreas boreas) against Batrachochytrium dendrobatidis (Bd), the amphibian fungal pathogen responsible for the greatest host biodiversity loss due to disease. Boreal toads have severely declined in Colorado (CO) due to Bd, but toad populations challenged with Bd in western Wyoming (WY) appear to be less affected. We used a common garden infection experiment to expose post-metamorphic toads sourced from four populations (2 in CO and 2 in WY) to Bd and monitored changes in mass, pathogen burden and survival for 8 weeks. We used a multi-state modelling approach to estimate weekly survival and transition probabilities between infected and cleared states, reflecting a dynamic infection process that traditional approaches fail to capture. We found that WY boreal toads are more tolerant to Bd infection with higher survival probabilities than those in CO when infected with identical pathogen burdens. WY toads also appeared more resistant to Bd with a higher probability of infection clearance and an average of 5 days longer to reach peak infection burdens. Our results demonstrate strong intraspecific differences in tolerance and resistance that likely contribute to why population declines vary regionally across this species. Our multi-state framework allowed us to gain inference on typically hidden disease processes when testing for host tolerance or resistance. Our findings demonstrate that describing an entire host species as 'tolerant' or 'resistant' (or lack thereof) is unwise without testing for intraspecific variation.

Best practices for genetic and genomic data archiving.

Genetic and genomic data are collected for a vast array of scientific and applied purposes. Despite mandates for public archiving, data are typically used only by the generating authors. The reuse of genetic and genomic datasets remains uncommon because it is difficult, if not impossible, due to non-standard archiving practices and lack of contextual metadata. But as the new field of macrogenetics is demonstrating, if genetic data and their metadata were more accessible and FAIR (findable, accessible, interoperable and reusable) compliant, they could be reused for many additional purposes. We discuss the main challenges with existing genetic and genomic data archives, and suggest best practices for archiving genetic and genomic data. Recognizing that this is a longstanding issue due to little formal data management training within the fields of ecology and evolution, we highlight steps that research institutions and publishers could take to improve data archiving.

Thermal limits along tropical elevational gradients: Poison frog tadpoles show plasticity but maintain divergence across elevation.

Temperature is arguably one of the most critical environmental factors impacting organisms at molecular, organismal, and ecological levels. Temperature variation across elevation may cause divergent selection in physiological critical thermal limits (CTMAX and CTMIN). Generally, high elevation populations are predicted to withstand lower environmental temperatures than low elevation populations. Organisms can also exhibit phenotypic plasticity when temperature varies, although theory and empirical evidence suggest that tropical ectotherms have relatively limited ability to acclimate. To study the effect of temperature variation along elevational transects on thermal limits, we measured CTMAX and CTMIN of 934 tadpoles of a poison frog species, Epipedobates anthonyi, along two elevational gradients (200-1700 m asl) in southwestern Ecuador to investigate their thermal tolerance across elevation. We also tested if tadpoles could plastically shift their critical thermal limits in response to exposure to different temperatures representing the range of temperatures they experience in nature (20 °C, 24 °C, and 28 °C). Overall, we found that CTMAX did not change across elevation. In contrast, CTMIN was lower at higher elevations, suggesting that elevational variation in temperature influences this thermal trait. Moreover, all populations shifted their CTMAX and CTMIN according to treatment temperatures, demonstrating an acclimation response. Overall, trends in CTMIN among high, mid, and low elevation populations were maintained despite plastic responses to treatment temperature. These results demonstrate that, for tadpoles of E. anthonyi across tropical elevational gradients, temperature acts as a selective force for CTMIN, even when populations show acclimation abilities in both, CTMAX and CTMIN. Our findings advance our understanding on how environmental variation affects organisms' evolutionary trajectories and their abilities to persist in a changing climate in a tropical biodiversity hotspot.

Unraveling Adaptive Evolutionary Divergence at Microgeographic Scales.

AbstractStriking examples of local adaptation at fine geographic scales are increasingly being documented in natural populations. However, the relative contributions made by natural selection, phenotype-dependent dispersal (when individuals disperse with respect to a habitat preference), and mate preference in generating and maintaining microgeographic adaptation and divergence are not well studied. Here, we develop quantitative genetics models and individual-based simulations (IBSs) to uncover the evolutionary forces that possibly drive microgeographic divergence. We also perform Bayesian estimation of the parameters in our IBS using empirical data on habitat-specific variation in bill morphology in the island scrub-jay (Aphelocoma insularis) to apply our models to a natural system. We find that natural selection and phenotype-dependent dispersal can generate the patterns of divergence we observe in the island scrub-jay. However, mate preference for a mate with similar bill morphology, even though observed in the species, does not play a significant role in driving divergence. Our modeling approach provides insights into phenotypic evolution occurring over small spatial scales relative to dispersal ranges, suggesting that adaptive divergence at microgeographic scales may be common across a wider range of taxa than previously thought. Our quantitative genetic models help to inform future theoretical and empirical work to determine how selection, habitat preference, and mate preference contribute to local adaptation and microgeographic divergence.

Effects of hunting on mating, relatedness, and genetic diversity in a puma population.

Hunting mortality can affect population abundance, demography, patterns of dispersal and philopatry, breeding, and genetic diversity. We investigated the effects of hunting on the reproduction and genetic diversity in a puma population in western Colorado, USA. We genotyped over 11,000 single nucleotide polymorphisms (SNPs), using double-digest, restriction site-associated DNA sequencing (ddRADseq) in 291 tissue samples collected as part of a study on the effects of hunting on puma population abundance and demography in Colorado from 2004 to 2014. The study was designed with a reference period (years 1-5), during which hunting was suspended, followed by a treatment period (years 6-10), in which hunting was reinstated. Our objectives were to examine the effects of hunting on: (1) paternity and male reproductive success; (2) the relatedness between pumas within the population, and (3) genetic diversity. We found that hunting reduced the average age of male breeders. The number of unique fathers siring litters increased each year without hunting and decreased each year during the hunting period. Mated pairs were generally unrelated during both time periods, and females were more closely related than males. Hunting was also associated with increased relatedness among males and decreased relatedness among females in the population. Finally, genetic diversity increased during the period without hunting and decreased each year when hunting was present. This study demonstrates the utility of merging demographic data with large-scale genomic datasets in order to better understand the consequences of management actions. Specifically, we believe that this study highlights the need for long-term experimental research in which hunting mortality is manipulated, including at least one non-harvested control population, as part of a broader adaptive, zone management scheme.

Conservation genomics of an endangered montane amphibian reveals low population structure, low genomic diversity and selection pressure from disease.

Wildlife diseases are a major global threat to biodiversity. Boreal toads (Anaxyrus [Bufo] boreas) are a state-endangered species in the southern Rocky Mountains of Colorado and New Mexico, and a species of concern in Wyoming, largely due to lethal skin infections caused by the amphibian chytrid fungus Batrachochytrium dendrobatidis (Bd). We performed conservation and landscape genomic analyses using single nucleotide polymorphisms from double-digest, restriction site-associated DNA sequencing in combination with the development of the first boreal toad (and first North American toad) reference genome to investigate population structure, genomic diversity, landscape connectivity and adaptive divergence. Genomic diversity (π = 0.00034-0.00040) and effective population sizes (Ne = 8.9-38.4) were low, likely due to post-Pleistocene founder effects and Bd-related population crashes over the last three decades. Population structure was also low, likely due to formerly high connectivity among a higher density of geographically proximate populations. Boreal toad gene flow was facilitated by low precipitation, cold minimum temperatures, less tree canopy, low heat load and less urbanization. We found >8X more putatively adaptive loci related to Bd intensity than to all other environmental factors combined, and evidence for genes under selection related to immune response, heart development and regulation and skin function. These data suggest boreal toads in habitats with Bd have experienced stronger selection pressure from disease than from other, broad-scale environmental variations. These findings can be used by managers to conserve and recover the species through actions including reintroduction and supplementation of populations that have declined due to Bd.

Observation of Collider Muon Neutrinos with the SND@LHC Experiment.

We report the direct observation of muon neutrino interactions with the SND@LHC detector at the Large Hadron Collider. A dataset of proton-proton collisions at sqrt[s]=13.6 TeV collected by SND@LHC in 2022 is used, corresponding to an integrated luminosity of 36.8 fb^{-1}. The search is based on information from the active electronic components of the SND@LHC detector, which covers the pseudorapidity region of 7.2<η<8.4, inaccessible to the other experiments at the collider. Muon neutrino candidates are identified through their charged-current interaction topology, with a track propagating through the entire length of the muon detector. After selection cuts, 8 ν_{µ} interaction candidate events remain with an estimated background of 0.086 events, yielding a significance of about 7 standard deviations for the observed ν_{µ} signal.

Linking critical thermal maximum to mortality from thermal stress in a cold-water frog.

Estimates of organismal thermal tolerance are frequently used to assess physiological risk from warming, yet the assumption that these estimates are predictive of mortality has been called into question. We tested this assumption in the cold-water-specialist frog, Ascaphus montanus. For seven populations, we used dynamic experimental assays to measure tadpole critical thermal maximum (CTmax) and measured mortality from chronic thermal stress for 3 days at different temperatures. We tested the relationship between previously estimated population CTmax and observed mortality, as well as the strength of CTmax as a predictor of mortality compared to local stream temperatures capturing varying timescales. Populations with higher CTmax experienced significantly less mortality in the warmest temperature treatment (25°C). We also found that population CTmax outperformed stream temperature metrics as the top predictor of observed mortality. These results demonstrate a clear link between CTmax and mortality from thermal stress, contributing evidence that CTmax is a relevant metric for physiological vulnerability assessments.

Where to draw the line? Expanding the delineation of conservation units to highly mobile taxa.

Conservation units (CUs) are an essential tool for maximizing evolutionary potential and prioritizing areas across a species' range for protection when implementing conservation and management measures. However, current workflows for identifying CUs on the basis of neutral and adaptive genomic variation largely ignore information contained in patterns of isolation by distance (IBD), frequently the primary signal of population structure in highly mobile taxa, such as birds, bats, and marine organisms with pelagic larval stages. While individuals located on either end of a species' distribution may exhibit clear genetic, phenotypic, and ecological differences, IBD produces subtle changes in allele frequencies across space, making it difficult to draw clear boundaries for conservation purposes in the absence of discrete population structure. Here, we highlight potential pitfalls that arise when applying common methods for delineating CUs to continuously distributed organisms and review existing methods for detecting subtle breakpoints in patterns of IBD that can indicate barriers to gene flow in highly mobile taxa. In addition, we propose a new framework for identifying CUs in all organisms, including those characterized by continuous genomic differentiation, and suggest several possible ways to harness the information contained in patterns of IBD to guide conservation and management decisions.

Genetic diversity goals and targets have improved, but remain insufficient for clear implementation of the post-2020 global biodiversity framework.

Genetic diversity among and within populations of all species is necessary for people and nature to survive and thrive in a changing world. Over the past three years, commitments for conserving genetic diversity have become more ambitious and specific under the Convention on Biological Diversity's (CBD) draft post-2020 global biodiversity framework (GBF). This Perspective article comments on how goals and targets of the GBF have evolved, the improvements that are still needed, lessons learned from this process, and connections between goals and targets and the actions and reporting that will be needed to maintain, protect, manage and monitor genetic diversity. It is possible and necessary that the GBF strives to maintain genetic diversity within and among populations of all species, to restore genetic connectivity, and to develop national genetic conservation strategies, and to report on these using proposed, feasible indicators.

Habitat connectivity and host relatedness influence virus spread across an urbanising landscape in a fragmentation-sensitive carnivore.

Spatially heterogeneous landscape factors such as urbanisation can have substantial effects on the severity and spread of wildlife diseases. However, research linking patterns of pathogen transmission to landscape features remains rare. Using a combination of phylogeographic and machine learning approaches, we tested the influence of landscape and host factors on feline immunodeficiency virus (FIVLru) genetic variation and spread among bobcats (Lynx rufus) sampled from coastal southern California. We found evidence for increased rates of FIVLru lineage spread through areas of higher vegetation density. Furthermore, single-nucleotide polymorphism (SNP) variation among FIVLru sequences was associated with host genetic distances and geographic location, with FIVLru genetic discontinuities precisely correlating with known urban barriers to host dispersal. An effect of forest land cover on FIVLru SNP variation was likely attributable to host population structure and differences in forest land cover between different populations. Taken together, these results suggest that the spread of FIVLru is constrained by large-scale urban barriers to host movement. Although urbanisation at fine spatial scales did not appear to directly influence virus transmission or spread, we found evidence that viruses transmit and spread more quickly through areas containing higher proportions of natural habitat. These multiple lines of evidence demonstrate how urbanisation can change patterns of contact-dependent pathogen transmission and provide insights into how continued urban development may influence the incidence and management of wildlife disease.

Search for Dark Matter Axions with CAST-CAPP.

The CAST-CAPP axion haloscope, operating at CERN inside the CAST dipole magnet, has searched for axions in the 19.74 µeV to 22.47 µeV mass range. The detection concept follows the Sikivie haloscope principle, where Dark Matter axions convert into photons within a resonator immersed in a magnetic field. The CAST-CAPP resonator is an array of four individual rectangular cavities inserted in a strong dipole magnet, phase-matched to maximize the detection sensitivity. Here we report on the data acquired for 4124 h from 2019 to 2021. Each cavity is equipped with a fast frequency tuning mechanism of 10 MHz/ min between 4.774 GHz and 5.434 GHz. In the present work, we exclude axion-photon couplings for virialized galactic axions down to gaγγ = 8 × 10-14 GeV-1 at the 90% confidence level. The here implemented phase-matching technique also allows for future large-scale upgrades.

Apathogenic proxies for transmission dynamics of a fatal virus.

Identifying drivers of transmission-especially of emerging pathogens-is a formidable challenge for proactive disease management efforts. While close social interactions can be associated with microbial sharing between individuals, and thereby imply dynamics important for transmission, such associations can be obscured by the influences of factors such as shared diets or environments. Directly-transmitted viral agents, specifically those that are rapidly evolving such as many RNA viruses, can allow for high-resolution inference of transmission, and therefore hold promise for elucidating not only which individuals transmit to each other, but also drivers of those transmission events. Here, we tested a novel approach in the Florida panther, which is affected by several directly-transmitted feline retroviruses. We first inferred the transmission network for an apathogenic, directly-transmitted retrovirus, feline immunodeficiency virus (FIV), and then used exponential random graph models to determine drivers structuring this network. We then evaluated the utility of these drivers in predicting transmission of the analogously transmitted, pathogenic agent, feline leukemia virus (FeLV), and compared FIV-based predictions of outbreak dynamics against empirical FeLV outbreak data. FIV transmission was primarily driven by panther age class and distances between panther home range centroids. FIV-based modeling predicted FeLV dynamics similarly to common modeling approaches, but with evidence that FIV-based predictions captured the spatial structuring of the observed FeLV outbreak. While FIV-based predictions of FeLV transmission performed only marginally better than standard approaches, our results highlight the value of proactively identifying drivers of transmission-even based on analogously-transmitted, apathogenic agents-in order to predict transmission of emerging infectious agents. The identification of underlying drivers of transmission, such as through our workflow here, therefore holds promise for improving predictions of pathogen transmission in novel host populations, and could provide new strategies for proactive pathogen management in human and animal systems.

Genomics-informed delineation of conservation units in a desert amphibian.

Delineating conservation units (CUs, e.g., evolutionarily significant units, ESUs, and management units, MUs) is critical to the recovery of declining species because CUs inform both listing status and management actions. Genomic data have strengths and limitations in informing CU delineation and related management questions in natural systems. We illustrate the value of using genomic data in combination with landscape, dispersal and occupancy data to inform CU delineation in Nevada populations of the Great Basin Distinct Population Segment of the Columbia spotted frog (Rana luteiventris). R. luteiventris occupies naturally fragmented aquatic habitats in this xeric region, but beaver removal, climate change and other factors have put many of these populations at high risk of extirpation without management intervention. We addressed three objectives: (i) assessing support for ESUs within Nevada; (ii) evaluating and revising, if warranted, the current delineation of MUs; and (iii) evaluating genetic diversity, effective population size, adaptive differentiation and functional connectivity to inform ongoing management actions. We found little support for ESUs within Nevada but did identify potential revisions to MUs based on unique landscape drivers of connectivity that distinguish these desert populations from those in the northern portion of the species range. Effective sizes were uniformly small, with low genetic diversity and weak signatures of adaptive differentiation. Our findings suggest that management actions, including translocations and genetic rescue, might be warranted. Our study illustrates how a carefully planned genetic study, designed to address priority management goals that include CU delineation, can provide multiple insights to inform conservation action.

Global genetic diversity status and trends: towards a suite of Essential Biodiversity Variables (EBVs) for genetic composition.

Biodiversity underlies ecosystem resilience, ecosystem function, sustainable economies, and human well-being. Understanding how biodiversity sustains ecosystems under anthropogenic stressors and global environmental change will require new ways of deriving and applying biodiversity data. A major challenge is that biodiversity data and knowledge are scattered, biased, collected with numerous methods, and stored in inconsistent ways. The Group on Earth Observations Biodiversity Observation Network (GEO BON) has developed the Essential Biodiversity Variables (EBVs) as fundamental metrics to help aggregate, harmonize, and interpret biodiversity observation data from diverse sources. Mapping and analyzing EBVs can help to evaluate how aspects of biodiversity are distributed geographically and how they change over time. EBVs are also intended to serve as inputs and validation to forecast the status and trends of biodiversity, and to support policy and decision making. Here, we assess the feasibility of implementing Genetic Composition EBVs (Genetic EBVs), which are metrics of within-species genetic variation. We review and bring together numerous areas of the field of genetics and evaluate how each contributes to global and regional genetic biodiversity monitoring with respect to theory, sampling logistics, metadata, archiving, data aggregation, modeling, and technological advances. We propose four Genetic EBVs: (i) Genetic Diversity; (ii) Genetic Differentiation; (iii) Inbreeding; and (iv) Effective Population Size (Ne ). We rank Genetic EBVs according to their relevance, sensitivity to change, generalizability, scalability, feasibility and data availability. We outline the workflow for generating genetic data underlying the Genetic EBVs, and review advances and needs in archiving genetic composition data and metadata. We discuss how Genetic EBVs can be operationalized by visualizing EBVs in space and time across species and by forecasting Genetic EBVs beyond current observations using various modeling approaches. Our review then explores challenges of aggregation, standardization, and costs of operationalizing the Genetic EBVs, as well as future directions and opportunities to maximize their uptake globally in research and policy. The collection, annotation, and availability of genetic data has made major advances in the past decade, each of which contributes to the practical and standardized framework for large-scale genetic observation reporting. Rapid advances in DNA sequencing technology present new opportunities, but also challenges for operationalizing Genetic EBVs for biodiversity monitoring regionally and globally. With these advances, genetic composition monitoring is starting to be integrated into global conservation policy, which can help support the foundation of all biodiversity and species' long-term persistence in the face of environmental change. We conclude with a summary of concrete steps for researchers and policy makers for advancing operationalization of Genetic EBVs. The technical and analytical foundations of Genetic EBVs are well developed, and conservation practitioners should anticipate their increasing application as efforts emerge to scale up genetic biodiversity monitoring regionally and globally.

Habitat-linked genetic variation supports microgeographic adaptive divergence in an island-endemic bird species.

We investigated the potential mechanisms driving habitat-linked genetic divergence within a bird species endemic to a single 250-km2 island. The island scrub-jay (Aphelocoma insularis) exhibits microgeographic divergence in bill morphology across pine-oak ecotones on Santa Cruz Island, California (USA), similar to adaptive differences described in mainland congeners over much larger geographic scales. To test whether individuals exhibit genetic differentiation related to habitat type and divergence in bill length, we genotyped over 3000 single nucleotide polymorphisms in 123 adult island scrub-jay males from across Santa Cruz Island using restriction site-associated DNA sequencing. Neutral landscape genomic analyses revealed that genome-wide genetic differentiation was primarily related to geographic distance and differences in habitat composition. We also found 168 putatively adaptive loci associated with habitat type using multivariate redundancy analysis while controlling for spatial effects. Finally, two genome-wide association analyses revealed a polygenic basis to variation in bill length with multiple loci detected in or near genes known to affect bill morphology in other birds. Our findings support the hypothesis that divergent selection at microgeographic scales can cause adaptive divergence in the presence of ongoing gene flow.

Hunting alters viral transmission and evolution in a large carnivore.

Hunting can fundamentally alter wildlife population dynamics but the consequences of hunting on pathogen transmission and evolution remain poorly understood. Here, we present a study that leverages a unique landscape-scale quasi-experiment coupled with pathogen-transmission tracing, network simulation and phylodynamics to provide insights into how hunting shapes feline immunodeficiency virus (FIV) dynamics in puma (Puma concolor). We show that removing hunting pressure enhances the role of males in transmission, increases the viral population growth rate and increases the role of evolutionary forces on the pathogen compared to when hunting was reinstated. Changes in transmission observed with the removal of hunting could be linked to short-term social changes while the male puma population increased. These findings are supported through comparison with a region with stable hunting management over the same time period. This study shows that routine wildlife management can have impacts on pathogen transmission and evolution not previously considered.

Parasites as conservation tools.

Parasite success typically depends on a close relationship with one or more hosts; therefore, attributes of parasitic infection have the potential to provide indirect details of host natural history and are biologically relevant to animal conservation. Characterization of parasite infections has been useful in delineating host populations and has served as a proxy for assessment of environmental quality. In other cases, the utility of parasites is just being explored, for example, as indicators of host connectivity. Innovative studies of parasite biology can provide information to manage major conservation threats by using parasite assemblage, prevalence, or genetic data to provide insights into the host. Overexploitation, habitat loss and fragmentation, invasive species, and climate change are major threats to animal conservation, and all of these can be informed by parasites.

Los Parásitos como Herramienta de Conservación Resumen El éxito de los parásitos depende típicamente de la relación cercana con uno o más hospederos; por lo tanto, las características de la infección parasitaria tienen potencial para proporcionar detalles indirectos de la historia natural del hospedero y son biológicamente relevantes para la conservación animal. La caracterización de las infecciones parasitarias ha sido útil para definir a las poblaciones hospederas y ha servido como sustituto para la evaluación de la calidad ambiental. Los estudios innovadores de la biología de parásitos pueden proporcionar información para manejar las principales amenazas a la conservación mediante la información proporcionada por el conjunto de parásitos, su prevalencia o genética que proporciona conocimiento sobre el hospedero. La sobreexplotación, la pérdida del hábitat y la fragmentación, las especies invasoras y el cambio climático son las principales amenazas para la conservación animal y a todas pueden ser informadas mediante los parásitos.

Contrasting environmental drivers of genetic and phenotypic divergence in an Andean poison frog (Epipedobates anthonyi).

Phenotypic and genetic divergence are shaped by the homogenizing effects of gene flow and the differentiating processes of genetic drift and local adaptation. Herein, we examined the mechanisms that underlie phenotypic (size and color) and genetic divergence in 35 populations (535 individuals) of the poison frog Epipedobates anthonyi along four elevational gradients (0-1800 m asl) in the Ecuadorian Andes. We found phenotypic divergence in size and color despite relatively low genetic divergence at neutral microsatellite loci. Genetic and phenotypic divergence were both explained by landscape resistance between sites (isolation-by-resistance, IBR), likely due to a cold and dry mountain ridge between the northern and southern elevational transects that limits dispersal and separates two color morphs. Moreover, environmental differences among sites also explained genetic and phenotypic divergence, suggesting isolation-by-environment (IBE). When northern and southern transects were analyzed separately, genetic divergence was predicted either by distance (isolation-by-distance, IBD; northern) or environmental resistance between sites (IBR; southern). In contrast, phenotypic divergence was primarily explained by environmental differences among sites, supporting the IBE hypothesis. These results indicate that although distance and geographic barriers are important drivers of population divergence, environmental variation has a two-fold effect on population divergence. On the one hand, landscape resistance between sites reduces gene flow (IBR), while on the other hand, environmental differences among sites exert divergent selective pressures on phenotypic traits (IBE). Our work highlights the importance of studying both genetic and phenotypic divergence to better understand the processes of population divergence and speciation along ecological gradients.

Adaptive divergence in bill morphology and other thermoregulatory traits is facilitated by restricted gene flow in song sparrows on the California Channel Islands.

Disentangling the effects of neutral and adaptive processes in maintaining phenotypic variation across environmental gradients is challenging in natural populations. Song sparrows (Melospiza melodia) on the California Channel Islands occupy a pronounced east-west climate gradient within a small spatial scale, providing a unique opportunity to examine the interaction of genetic isolation (reduced gene flow) and the environment (selection) in driving variation. We used reduced representation genomic libraries to infer the role of neutral processes (drift and restricted gene flow) and divergent selection in driving variation in thermoregulatory traits with an emphasis on the mechanisms that maintain bill divergence among islands. Analyses of 22,029 neutral SNPs confirm distinct population structure by island with restricted gene flow and relatively large effective population sizes, suggesting bill differences are probably not a product of genetic drift. Instead, we found strong support for local adaptation using 3294 SNPs in differentiation-based and environmental association analyses coupled with genome-wide association tests. Specifically, we identified several putatively adaptive and candidate loci in or near genes involved in bill development pathways (e.g., BMP, CaM, Wnt), confirming the highly complex and polygenic architecture underlying bill morphology. Furthermore, we found divergence in genes associated with other thermoregulatory traits (i.e., feather structure, plumage colour, and physiology). Collectively, these results suggest strong divergent selection across an island archipelago results in genomic changes in a suite of traits associated with climate adaptation over small spatial scales. Future research should move beyond studying univariate traits to better understand multidimensional responses to complex environmental conditions.