Local adaptation is a major driver of biological diversity, and related species may develop analogous (parallel evolution) or alternative (divergent evolution) solutions to similar ecological challenges. We expect these adaptive solutions between closely related organisms would culminate in both phenotypic and genotypic signals. In this study, we employ a reciprocal transplant, glasshouse experiment with two Eucalyptus species ( E. grandis and E. tereticornis) with large, overlapping distributions grown under contrasting ‘local’ temperature conditions (tropic and temperate) to investigate the independent contribution of adaptation, plasticity, and their interaction at molecular, physiological and morphological levels. We find key traits differ in their response. The link between gene expression and traits markedly differed between species. Divergent evolution was the dominant pattern driving adaptation as unique gene responses (91% of all significant genes) was the greatest factor driving differentiation; but overlapping gene (homologous) responses were dependent on the determining factor (plastic, adaptive, or genotype by environment interaction). 98% of the plastic homologs were similarly regulated, while 50% of the adaptive homologs and 100% of the interaction homologs were antagonistically regulated. Therefore, parallel evolution for the adaptive effect in homologous genes was greater than expected but not in favour of divergent evolution. Further, heat shock proteins for E. grandis were almost entirely driven by adaptive responses, while plasticity drove the response in E. tereticornis. These results suggest divergent molecular evolutionary solutions dominated the adaptive mechanisms among species, even in similar ecological circumstances. Thus, trees with overlapping distributions are unlikely to equally persist in the future, suggesting that management of future forests to changing temperature conditions must be species specific.

Temperature and precipitation regimes are rapidly changing, resulting in forest dieback and local extinction events, particularly in Mediterranean-type climates. Strategic forest management approaches that enhance forests’ resilience to future climates are urgently required, however adaptation to climates in heterogeneous landscapes with multiple selection pressures may be complex. For widespread trees in Mediterranean-type climates we hypothesized that patterns of local adaptation are associated with climate; precipitation is a stronger factor of adaptation than temperature; functionally related genes show similar signatures of adaptation; and adaptive variants are independently sorting across the landscape. To test our hypotheses, we sampled 28 populations across the geographic and climatic distribution of Eucalyptus marginata (jarrah), in south-west Western Australia, and obtained 13,534 independent single nucleotide polymorphic (SNP) markers across the genome. While overall levels of population differentiation were low (FST=0.04), environmental association analyses found a total of 2,336 unique SNPs potentially associated with five climate variables of temperature and precipitation. Allelic turnover was identified for SNPs associated with temperate seasonality and mean precipitation of the warmest quarter (39.2% and 36.9% deviance explained, respectively), suggesting that both temperature and precipitation are important factors in adaptation. SNPs within similarly function genes, according to gene ontology enrichment analysis, had analogous allelic turnover along climate gradients, while SNPs among temperature and precipitation variables had orthogonal patterns of adaptation. These contrasting patterns of adaptation provide evidence that there may be standing genomic variation adapted to changing climates, providing the substrate needed to promote adaptive management strategies to bolster forest resilience in the future.

Extreme drought conditions across the globe are impacting biodiversity with serious implications for the persistence of native species. However, quantitative data on drought tolerance is not available for diverse flora to inform conservation management. We quantified physiological drought tolerance in the diverse Hakea genus (Proteaceae) to test predictions based on climatic-origin, life history and functional traits. We sampled terminal branches of replicate plants of 16 species in a common garden. Xylem cavitation was induced in branches under varying water potential (tension) in a centrifuge and the tension generating 50% loss of conductivity (stem P50) was characterized as a metric for drought tolerance. The same branches were used to estimate plant functional traits, including wood density, specific leaf area, and Huber value (sap flow area to leaf area ratio). There was significant variation in stem P50 among species, which was negatively associated with the species climate-origin (rainfall and aridity). Drought tolerance did not differ among life histories; however, a drought avoidance strategy with terete leaf form and greater Huber value may be important for species to colonize and persist in the arid biome. Our findings will contribute to future prediction of species vulnerability to drought and adaptive management under climate change.