Comparative tests of the species‐genetic diversity correlation at neutral and nonneutral loci in four species of stream insect

A fundamental question linking population genetics and community ecology is how adaptive processes (e.g., natural selection) and neutral processes (e.g., drift‐migration equilibrium) underpin the species‐genetic diversity correlation (SGDC). Here, we combine genome scans and outlier loci detection with community analysis to separately test for neutral and nonneutral SGDCs in four species of stream insect. We sampled 60 localities in Japan and examined the relationships among population AFLP band richness (Br), taxon richness of the total community (S) and of the trophic guild (S_tr), and 15 habitat parameters that could potentially drive adaptation and influence richness. Neutral Br was positively correlated with S only in the dominant species of these communities, suggesting Br may be constrained when intraspecific competition is pronounced. Nonneutral Br was correlated with S_tr in a species restricted to high elevations where habitat heterogeneity was highest. Community distance and genetic distance (β‐SGDC) was correlated in two of the four species at both neutral and nonneutral loci. Distance‐based redundancy analysis found geographic isolation and elevation to drive divergence of both communities and populations. This suggests that both neutral and adaptive divergence occurred through the shared influences of geographic isolation and local adaptation at the two levels of diversity.     

Floral integration,phenotypic covariance structure and pollinator variation in bumblebee‐pollinated Helleborus foetidus

By analysing patterns of phenotypic integration and multivariate covariance structure of five metric floral traits in nine Iberian populations of bumblebee‐pollinated Helleborus foetidus (Ranunculaceae), this paper attempts to test the general hypothesis that pollinators enhance floral integration and selectively modify phenotypic correlations between functionally linked floral traits. The five floral traits examined exhibited significant phenotypic integration at all populations, and both the magnitude and the pattern of integration differed widely among populations. Variation in extent and pattern of integration was neither distance‐dependent nor significantly related to between‐population variation in taxonomical composition and morphological diversity of the pollinator assemblage. Patterns of floral integration were closer to expectations derived from consideration of developmental affinities between floral whorls than to expectations based on a pollinator‐mediated adaptive hypothesis. Taken together, results of this study suggest that between‐population differences in magnitude and pattern of floral integration in H. foetidus are probably best explained as a consequence of random genetic sampling in the characteristically small and ephemeral populations of this species, rather than reflecting the selective action of current pollinators.     

Population genomics reveals multiple drivers of population differentiation in a sex‐role‐reversed pipefish

A major goal of molecular ecology is to identify the causes of genetic and phenotypic differentiation among populations. Population genomics is suitably poised to tackle these key questions by diagnosing the evolutionary mechanisms driving divergence in nature. Here, we set out to investigate the evolutionary processes underlying population differentiation in the Gulf pipefish, Syngnathus scovelli. We sampled approximately 50 fish from each of 12 populations distributed from the Gulf coast of Texas to the Atlantic coast of Florida and performed restriction‐site‐associated DNA sequencing to identify SNPs throughout the genome. After imposing quality and stringency filters, we selected a panel of 6348 SNPs present in all 12 populations, 1753 of which were not physically linked. We identified a genome‐wide pattern of isolation by distance, in addition to a more substantial genetic break separating populations in the Gulf of Mexico from those in the Atlantic. We also used several divergence outlier approaches and tests for genotype–environment correlations to identify 400 SNPs putatively involved in local adaptation. Patterns of phenotypic differentiation and variation diverged from the overall genomic pattern, suggesting that selection, phenotypic plasticity or demographic factors may be shaping phenotypes in distinct populations. Overall, our results suggest that population divergence is driven by a variety of factors in S. scovelli, including neutral processes and selection on multiple traits.     

Outlier SNP markers reveal fine‐scale genetic structuring across European hake populations (Merluccius merluccius)

Shallow population structure is generally reported for most marine fish and explained as a consequence of high dispersal, connectivity and large population size. Targeted gene analyses and more recently genome‐wide studies have challenged such view, suggesting that adaptive divergence might occur even when neutral markers provide genetic homogeneity across populations. Here, 381 SNPs located in transcribed regions were used to assess large‐ and fine‐scale population structure in the European hake (Merluccius merluccius), a widely distributed demersal species of high priority for the European fishery. Analysis of 850 individuals from 19 locations across the entire distribution range showed evidence for several outlier loci, with significantly higher resolving power. While 299 putatively neutral SNPs confirmed the genetic break between basins (F_CT = 0.016) and weak differentiation within basins, outlier loci revealed a dramatic divergence between Atlantic and Mediterranean populations (F_CT range 0.275–0.705) and fine‐scale significant population structure. Outlier loci separated North Sea and Northern Portugal populations from all other Atlantic samples and revealed a strong differentiation among Western, Central and Eastern Mediterranean geographical samples. Significant correlation of allele frequencies at outlier loci with seawater surface temperature and salinity supported the hypothesis that populations might be adapted to local conditions. Such evidence highlights the importance of integrating information from neutral and adaptive evolutionary patterns towards a better assessment of genetic diversity. Accordingly, the generated outlier SNP data could be used for tackling illegal practices in hake fishing and commercialization as well as to develop explicit spatial models for defining management units and stock boundaries.     

Machine‐learning‐based detection of adaptive divergence of the stream mayfly Ephemera strigata populations

Adaptive divergence is a key mechanism shaping the genetic variation of natural populations. A central question linking ecology with evolutionary biology is how spatial environmental heterogeneity can lead to adaptive divergence among local populations within a species. In this study, using a genome scan approach to detect candidate loci under selection, we examined adaptive divergence of the stream mayfly Ephemera strigata in the Natori River Basin in northeastern Japan. We applied a new machine‐learning method (i.e., random forest) besides traditional distance‐based redundancy analysis (dbRDA) to examine relationships between environmental factors and adaptive divergence at non‐neutral loci. Spatial autocorrelation analysis based on neutral loci was employed to examine the dispersal ability of this species. We conclude the following: (a) E. strigata show altitudinal adaptive divergence among the populations in the Natori River Basin; (b) random forest showed higher resolution for detecting adaptive divergence than traditional statistical analysis; and (c) separating all markers into neutral and non‐neutral loci could provide full insight into parameters such as genetic diversity, local adaptation, and dispersal ability.     

Lack of parallel genetic patterns underlying the repeated ecological divergence of beach and stream‐spawning kokanee salmon

Recent progress in methods for detecting adaptive population divergence in situ shows promise for elucidating the conditions under which selection acts to generate intraspecific diversity. Rapid ecological diversification is common in fishes; however, the role of phenotypic plasticity and adaptation to local environments is poorly understood. It is now possible to investigate genetic patterns to make inferences regarding phenotypic traits under selection and possible mechanisms underlying ecotype divergence, particularly where similar novel phenotypes have arisen in multiple independent populations. Here, we employed a bottom‐up approach to test for signatures of directional selection associated with divergence of beach‐ and stream‐spawning kokanee, the obligate freshwater form of sockeye salmon (Oncorhynchus nerka). Beach‐ and stream‐spawners co‐exist in many post‐glacial lakes and exhibit distinct reproductive behaviours, life‐history traits and spawning habitat preferences. Replicate ecotype pairs across five lakes in British Columbia, Canada were genotyped at 57 expressed sequence tag‐linked and anonymous microsatellite loci identified in a previous genome scan. Fifteen loci exhibited signatures of directional selection (high F_ST outliers), four of which were identified in multiple lakes. However, the lack of parallel genetic patterns across all lakes may be a result of: 1) an inability to detect loci truly under selection; 2) alternative genetic pathways underlying ecotype divergence in this system; and/or 3) phenotypic plasticity playing a formative role in driving kokanee spawning habitat differences. Gene annotations for detected outliers suggest pathogen resistance and energy metabolism as potential mechanisms contributing to the divergence of beach‐ and stream‐spawning kokanee, but further study is required.     

High degree of genetic differentiation in marine three‐spined sticklebacks (Gasterosteus aculeatus)

Populations of widespread marine organisms are typically characterized by a low degree of genetic differentiation in neutral genetic markers, but much less is known about differentiation in genes whose functional roles are associated with specific selection regimes. To uncover possible adaptive population divergence and heterogeneous genomic differentiation in marine three‐spined sticklebacks (Gasterosteus aculeatus), we used a candidate gene‐based genome‐scan approach to analyse variability in 138 microsatellite loci located within/close to (<6 kb) functionally important genes in samples collected from ten geographic locations. The degree of genetic differentiation in markers classified as neutral or under balancing selection—as determined with several outlier detection methods—was low (F_ST = 0.033 or 0.011, respectively), whereas average F_ST for directionally selected markers was significantly higher (F_ST = 0.097). Clustering analyses provided support for genomic and geographic heterogeneity in selection: six genetic clusters were identified based on allele frequency differences in the directionally selected loci, whereas four were identified with the neutral loci. Allelic variation in several loci exhibited significant associations with environmental variables, supporting the conjecture that temperature and salinity, but not optic conditions, are important drivers of adaptive divergence among populations. In general, these results suggest that in spite of the high degree of physical connectivity and gene flow as inferred from neutral marker genes, marine stickleback populations are strongly genetically structured in loci associated with functionally relevant genes.     

Evolution of population structure in an estuarine‐dependent marine fish

Restriction site‐associated DNA (RAD) sequencing was used to characterize neutral and adaptive genetic variation among geographic samples of red drum, Sciaenops ocellatus, an estuarine‐dependent fish found in coastal waters along the southeastern coast of the United States (Atlantic) and the northern Gulf of Mexico (Gulf). Analyses of neutral and outlier loci revealed three genetically distinct regional clusters: one in the Atlantic and two in the northern Gulf. Divergence in neutral loci indicated gradual genetic change and followed a linear pattern of isolation by distance. Divergence in outlier loci was at least an order of magnitude greater than divergence in neutral loci, and divergence between the regions in the Gulf was twice that of divergence between other regions. Discordance in patterns of genetic divergence between outlier and neutral loci is consistent with the hypothesis that the former reflects adaptive responses to environmental factors that vary on regional scales, while the latter largely reflects drift processes. Differences in basic habitat, initiated by glacial retreat and perpetuated by contemporary oceanic and atmospheric forces interacting with the geomorphology of the northern Gulf, followed by selection, appear to have led to reduced gene flow among red drum across the northern Gulf, reinforcing differences accrued during isolation and resulting in continued divergence across the genome. This same dynamic also may pertain to other coastal or nearshore fishes (18 species in 14 families) where genetically or morphologically defined sister taxa occur in the three regions.     

BRINGING HABITAT INFORMATION INTO STATISTICAL TESTS OF LOCAL ADAPTATION IN QUANTITATIVE TRAITS: A CASE STUDY OF NINE‐SPINED STICKLEBACKS

Detection of footprints of historical natural selection on quantitative traits in cross‐sectional data sets is challenging, especially when the number of populations to be compared is small and the populations are subject to strong random genetic drift. We extend a recent Bayesian multivariate approach to differentiate between selective and neutral causes of population differentiation by the inclusion of habitat information. The extended framework allows one to test for signals of selection in two ways: by comparing the patterns of population differentiation in quantitative traits and in neutral loci, and by comparing the similarity of habitats and phenotypes. We illustrate the framework using data on variation of eight morphological and behavioral traits among four populations of nine‐spined sticklebacks (Pungitius pungitius). In spite of the strong signal of genetic drift in the study system (average F_ST = 0.35 in neutral markers), strong footprints of adaptive population differentiation were uncovered both in morphological and behavioral traits. The results give quantitative support for earlier qualitative assessments, which have attributed the observed differentiation to adaptive divergence in response to differing ecological conditions in pond and marine habitats.     

Calls,colours, shape,and genes: a multi‐trait approach to the study of geographic variation in the Amazonian frog Allobates femoralis

Evolutionary divergence in behavioural traits related to mating may represent the initial stage of speciation. Direct selective forces are usually invoked to explain divergence in mate‐recognition traits, often neglecting a role for neutral processes or concomitant differentiation in ecological traits. We adopted a multi‐trait approach to obtain a deeper understanding of the mechanisms behind allopatric divergence in the Amazonian frog, Allobates femoralis. We tested the null hypothesis that geographic distance between populations correlates with genetic and phenotypic divergence, and compared divergence between mate‐recognition (acoustic) and ecological (coloration, body‐shape) traits. We quantified geographic variation in 39 phenotypic traits and a mitochondrial DNA marker among 125 individuals representing eight populations. Geographic variation in acoustic traits was pronounced and tracked the spatial genetic variation, which appeared to be neutral. Thus, the evolution of acoustic traits tracked the shared history of the populations, which is unexpected for pan‐Amazonian taxa or for mate‐recognition traits. Divergence in coloration appeared uncorrelated with genetic distance, and might be partly attributed to local selective pressures, and perhaps to Batesian mimicry. Divergence in body‐shape traits was low. The results obtained depict a complex evolutionary scenario and emphasize the importance of considering multiple traits when disentangling the forces behind allopatric divergence. ©2009 The Linnean Society of London, Biological Journal of the Linnean Society, 2009, 98 , 826–838.     

GENETIC AND PHENOTYPIC DIVERGENCE BETWEEN LOW‐ AND HIGH‐ALTITUDE POPULATIONS OF TWO RECENTLY DIVERGED CINNAMON TEAL SUBSPECIES

Spatial variation in the environment can lead to divergent selection between populations occupying different parts of a species’ range, and ultimately lead to population divergence. The colonization of new areas can thus facilitate divergence in beneficial traits, yet with little differentiation at neutral genetic markers. We investigated genetic and phenotypic patterns of divergence between low‐ and high‐altitude populations of cinnamon teal inhabiting normoxic and hypoxic regions in the Andes and adjacent lowlands of South America. Cinnamon teal showed strong divergence in body size (PC1; P_ST= 0.56) and exhibited significant frequency differences in a single nonsynonymous α‐hemoglobin amino acid polymorphism (Asn/Ser‐α9; F_ST= 0.60) between environmental extremes, despite considerable admixture of mtDNA and intron loci (F_ST= 0.004–0.168). Inferences of strong population segregation were further supported by the observation of few mismatched individuals in either environmental extreme. Coalescent analyses indicated that the highlands were most likely colonized from lowland regions but following divergence, gene flow has been asymmetric from the highlands into the lowlands. Multiple selection pressures associated with high‐altitude habitats, including cold and hypoxia, have likely shaped morphological and genetic divergence within South American cinnamon teal populations.     

Sexual selection and population divergence II. Divergence in different sexual traits and signal modalities in field crickets (Teleogryllus oceanicus)

Sexual selection can target many different types of traits. However, the relative influence of different sexually selected traits during evolutionary divergence is poorly understood. We used the field cricket Teleogryllus oceanicus to quantify and compare how five traits from each of three sexual signal modalities and components diverge among allopatric populations: male advertisement song, cuticular hydrocarbon (CHC) profiles and forewing morphology. Population divergence was unexpectedly consistent: we estimated the among‐population (genetic) variance‐covariance matrix, D , for all 15 traits, and D_max explained nearly two‐thirds of its variation. CHC and wing traits were most tightly integrated, whereas song varied more independently. We modeled the dependence of among‐population trait divergence on genetic distance estimated from neutral markers to test for signatures of selection versus neutral divergence. For all three sexual trait types, phenotypic variation among populations was largely explained by a neutral model of divergence. Our findings illustrate how phenotypic integration across different types of sexual traits might impose constraints on the evolution of mating isolation and divergence via sexual selection.     

Pleiotropic quantitative trait loci contribute to population divergence in traits associated with life-history variation in Mimulus guttatus

Evolutionary biologists seek to understand the genetic basis for multivariate phenotypic divergence. We constructed an F2 mapping population (N = 539) between two distinct populations of Mimulus guttatus. We measured 20 floral, vegetative, and life-history characters on parents and F1 and F2 hybrids in a common garden experiment. We employed multitrait composite interval mapping to determine the number, effect, and degree of pleiotropy in quantitative trait loci (QTL) affecting divergence in floral, vegetative, and life-history characters. We detected 16 QTL affecting floral traits; 7 affecting vegetative traits; and 5 affecting selected floral, vegetative, and life-history traits. Floral and vegetative traits are clearly polygenic. We detected a few major QTL, with all remaining QTL of small effect. Most detected QTL are pleiotropic, implying that the evolutionary shift between these annual and perennial populations is constrained. We also compared the genetic architecture controlling floral trait divergence both within (our intraspecific study) and between species, on the basis of a previously published analysis of M. guttatus and M. nasutus. Eleven of our 16 floral QTL map to approximately the same location in the interspecific map based on shared, collinear markers, implying that there may be a shared genetic basis for floral divergence within and among species of Mimulus.     

Genome‐wide association analyses reveal polygenic genomic architecture underlying divergent shell morphology in Spanish Littorina saxatilis ecotypes

Gene flow between diverging populations experiencing dissimilar ecological conditions can theoretically constrain adaptive evolution. To minimize the effect of gene flow, alleles underlying traits essential for local adaptation are predicted to be located in linked genome regions with reduced recombination. Local reduction in gene flow caused by selection is expected to produce elevated divergence in these regions. The highly divergent crab‐adapted and wave‐adapted ecotypes of the marine snail Littorina saxatilis present a model system to test these predictions. We used genome‐wide association (GWA) analysis of geometric morphometric shell traits associated with microgeographic divergence between the two L. saxatilis ecotypes within three separate sampling sites. A total of 477 snails that had individual geometric morphometric data and individual genotypes at 4,066 single nucleotide polymorphisms (SNPs) were analyzed using GWA methods that corrected for population structure among the three sites. This approach allowed dissection of the genomic architecture of shell shape divergence between ecotypes across a wide geographic range, spanning two glacial lineages. GWA revealed 216 quantitative trait loci (QTL) with shell size or shape differences between ecotypes, with most loci explaining a small proportion of phenotypic variation. We found that QTL were evenly distributed across 17 linkage groups, and exhibited elevated interchromosomal linkage, suggesting a genome‐wide response to divergent selection on shell shape between the two ecotypes. Shell shape trait‐associated loci showed partial overlap with previously identified outlier loci under divergent selection between the two ecotypes, supporting the hypothesis of diversifying selection on these genomic regions. These results suggest that divergence in shell shape between the crab‐adapted and wave‐adapted ecotypes is produced predominantly by a polygenic genomic architecture with positive linkage disequilibrium among loci of small effect.     

Adaptive differentiation in floral traits in the presence of high gene flow in scarlet gilia (Ipomopsis aggregata)

Plant–pollinator interactions are thought to be major drivers of floral trait diversity. However, the relative importance of divergent pollinator‐mediated selection vs. neutral processes in floral character evolution has rarely been explored. We tested for adaptive floral trait evolution by comparing differentiation at neutral genetic loci to differentiation at quantitative floral traits in a putative Ipomopsis aggregata hybrid zone. Typical I. aggregata subsp. candida displays slender white tubular flowers that are typical of flowers pollinated by hawkmoths, and subsp. collina displays robust red tubular flowers typical of flowers pollinated by hummingbirds; yet, hybrid flower morphs are abundant across the East Slope of the Colorado Rockies. We estimated genetic differentiation (F_ST) for nuclear and chloroplast microsatellite loci and used a half‐sib design to calculate quantitative trait divergence (Q_ST) from collection sites across the morphological hybrid zone. We found little evidence for population structure and estimated mean F_ST to be 0.032. Q_ST values for several floral traits including corolla tube length and width, colour, and nectar volume were large and significantly greater than mean F_ST. We performed multivariate comparisons of neutral loci to genetic correlations within and between populations and found a strong signal for divergent selection, suggesting that specific combinations of floral display and reward traits may be the targets of selection. Our results show little support for historical subspecies categories, yet floral traits are more diverged than expected due to drift alone. Non‐neutral divergence for multivariate quantitative traits suggests that selection by pollinators is maintaining a correlation between display and reward traits.     

Divergence at neutral and non-neutral loci in Drosophila buzzatii populations and their hybrids

The impact of intraspecific hybridisation on fitness and morphological traits depends on the history of natural selection and genetic drift, which may have led to differently coadapted gene-complexes in the parental populations. The divergence at neutral and non-neutral loci between populations can be evaluated by estimating F_ST and Q_ST respectively, and hence give an estimate of drift and selection in the populations. Here we investigate (1) whether divergence between populations in quantitative traits (wing size and shape) can be attributed to selection or drift alone, (2) The impact of intraspecific hybridisation on estimators for divergence at neutral (F_ST) and non-neutral loci (Q_ST) in hybrids, (3) If measurement of shape is more informative than size in order to detect divergence in quantitative traits between populations. The aims were addressed by performing two hybridisations between three populations of Drosophila buzzatii, one between populations from Argentina and the Canary Islands (separated for 200 years), and the other between populations from Argentina and Australia (separated for 80 years). We observed the highest divergence at neutral loci between the Argentinean and Canary Island populations, but highest morphological divergence between the Argentinean and Australian populations, indicating that natural selection is acting on the wings. Divergence based on Q_ST measures in the hybrids was sensitive towards increased phenotypic variance (σ²p) within groups and should be used with care when σ²p of populations differ. Our results indicate that measures of shape give a better estimate of divergence at the underlying quantitative traits loci than measures of size.     

Selection on outlier loci and their association with adaptive phenotypes in Littorina saxatilis contact zones

A fundamental issue in speciation research is to evaluate phenotypic variation and the genomics driving the evolution of reproductive isolation between sister taxa. Above all, hybrid zones are excellent study systems for researchers to examine the association of genetic differentiation, phenotypic variation and the strength of selection. We investigated two contact zones in the marine gastropod Littorina saxatilis and utilized landmark‐based geometric morphometric analysis together with amplified fragment length polymorphism (AFLP) markers to assess phenotypic and genomic divergence between ecotypes under divergent selection. From genetic markers, we calculated the cline width, linkage disequilibrium and the average effective selection on a locus. Additionally, we conducted an association analysis linking the outlier loci and phenotypic variation between ecotypes and show that a proportion of outlier loci are associated with key adaptive phenotypic traits.     

Parallels,nonparallels, and plasticity in population differentiation of threespine stickleback within a lake

The frequent occurrence of parallel phenotypic divergence in similar habitats is often evoked when emphasizing the role of ecology in adaptive radiation and speciation. However, because phenotypic plasticity can contribute to the observed pattern of divergence, confirmation of divergence at loci underlying phenotypic traits is important for confirming adaptive divergence. In the present study, we examine parallel morphological, neutral, and potentially adaptive genetic divergence of threespine stickleback inhabiting different habitats within a lake. Three genetic clusters best explained the neutral genetic structure within the lake; however, morphological differences were only weakly connected to genetic clusters and there was considerable phenotypic variation within clusters. Among the factors that could contribute to the observed pattern of morphological and genetic divergence are phenotypic plasticity, selective mortality of hybrids, and habitat choice based on morphology. Several loci are identified as outliers indicating divergent selection between the morphs and some parallels in morphological and adaptive genetic divergence are found in stickleback spawning at two lava sites. However, neutral genetic structure indicates considerable genetic connectivity among the two lava sites, and the parallels in morphology may therefore represent selective distribution of phenotypes rather than parallel divergence. © 2009 The Linnean Society of London, Biological Journal of the Linnean Society, 2009, 98 , 803–813.     

Environmental effects on molecular and phenotypic variation in populations of Eruca sativa across a steep climatic gradient

In Israel Eruca sativa has a geographically narrow distribution across a steep climatic gradient that ranges from mesic Mediterranean to hot desert environments. These conditions offer an opportunity to study the influence of the environment on intraspecific genetic variation. For this, we combined an analysis of neutral genetic markers with a phenotypic evaluation in common‐garden experiments, and environmental characterization of populations that included climatic and edaphic parameters, as well as geographic distribution. A Bayesian clustering of individuals from nine representative populations based on amplified fragment length polymorphism (AFLP) divided the populations into a southern and a northern geographic cluster, with one admixed population at the geographic border between them. Linear mixed models, with cluster added as a grouping factor, revealed no clear effects of environment or geography on genetic distances, but this may be due to a strong association of geography and environment with genetic clusters. However, environmental factors accounted for part of the phenotypic variation observed in the common‐garden experiments. In addition, candidate loci for selection were identified by association with environmental parameters and by two outlier methods. One locus, identified by all three methods, also showed an association with trichome density and herbivore damage, in net‐house and field experiments, respectively. Accordingly, we propose that because trichomes are directly linked to defense against both herbivores and excess radiation, they could potentially be related to adaptive variation in these populations. These results demonstrate the value of combining environmental and phenotypic data with a detailed genetic survey when studying adaptation in plant populations.     

Signatures of diversifying selection at EST‐SSR loci and association with climate in natural Eucalyptus populations

Understanding the environmental parameters that drive adaptation among populations is important in predicting how species may respond to global climatic changes and how gene pools might be managed to conserve adaptive genetic diversity. Here, we used Bayesian F_ST outlier tests and allele–climate association analyses to reveal two Eucalyptus EST‐SSR loci as strong candidates for diversifying selection in natural populations of a southwestern Australian forest tree, Eucalyptus gomphocephala (Myrtaceae). The Eucalyptus homolog of a CONSTANS‐like gene was an F_ST outlier, and allelic variation showed significant latitudinal clinal associations with annual and winter solar radiation, potential evaporation, summer precipitation and aridity. A second F_ST outlier locus, homologous to quinone oxidoreductase, was significantly associated with measures of temperature range, high summer temperature and summer solar radiation, with important implications for predicting the effect of temperature on natural populations in the context of climate change. We complemented these data with investigations into neutral population genetic structure and diversity throughout the species range. This study provides an investigation into selection signatures at gene‐homologous EST‐SSRs in natural Eucalyptus populations, and contributes to our understanding of the relationship between climate and adaptive genetic variation, informing the conservation of both putatively neutral and adaptive components of genetic diversity.