Genomic determinants, architecture, and constraints in drought-related traits in Corymbia calophylla.

BACKGROUND: Drought adaptation is critical to many tree species persisting under climate change, however our knowledge of the genetic basis for trees to adapt to drought is limited. This knowledge gap impedes our fundamental understanding of drought response and application to forest production and conservation. To improve our understanding of the genomic determinants, architecture, and trait constraints, we assembled a reference genome and detected ~ 6.5 M variants in 432 phenotyped individuals for the foundational tree Corymbia calophylla. RESULTS: We found 273 genomic variants determining traits with moderate heritability (h2SNP = 0.26-0.64). Significant variants were predominantly in gene regulatory elements distributed among several haplotype blocks across all chromosomes. Furthermore, traits were constrained by frequent epistatic and pleiotropic interactions. CONCLUSIONS: Our results on the genetic basis for drought traits in Corymbia calophylla have several implications for the ability to adapt to climate change: (1) drought related traits are controlled by complex genomic architectures with large haplotypes, epistatic, and pleiotropic interactions; (2) the most significant variants determining drought related traits occurred in regulatory regions; and (3) models incorporating epistatic interactions increase trait predictions. Our findings indicate that despite moderate heritability drought traits are likely constrained by complex genomic architecture potentially limiting trees response to climate change.

Urban greening with shrubs can supercharge invertebrate abundance and diversity.

In urban areas, diverse and complex habitats for biodiversity are often lacking. This lack of diversity not only compromises essential ecological processes, such as pollination and nutrient cycling, but also diminishes the resilience of urban ecosystems to pests and diseases. To enhance urban biodiversity, a possible solution is to integrate shrubs alongside trees, thereby increasing the overall amount of vegetation, structural complexity and the associated resource diversity. Here, using a common garden experiment involving a variety of trees and shrubs planted alone and in combination, we evaluate how canopy-associated invertebrate assemblages are influenced by vegetation type. In particular, we test whether the presence of shrubs, alone or with trees, results in increased abundance and taxonomic richness of invertebrates, compared to trees on their own. We found that the overall abundance of invertebrates, and that of specific functional groups (e.g., herbivores, pollinators, detritivores), was higher on shrubs, compared to trees, and when trees and shrubs were planted in combination (relative to trees on their own). Our results suggest that planting shrub and tree species with wide and dense crowns can increase the associated abundance and taxonomic and functional group richness of invertebrate communities. Overall, our findings indicate that urban planning would benefit from incorporating shrubs alongside urban trees to maximise invertebrate abundance, diversity and function in urban landscapes.

Germination temperature sensitivity differs between co-occurring tree species and climate origins resulting in contrasting vulnerability to global warming.

Climate change is shifting temperatures from historical patterns, globally impacting forest composition and resilience. Seed germination is temperature-sensitive, making the persistence of populations and colonization of available habitats vulnerable to warming. This study assessed germination response to temperature in foundation trees in south-western Australia's Mediterranean-type climate forests (Eucalyptus marginata (jarrah) and Corymbia calophylla (marri)) to estimate the thermal niche and vulnerability among populations. Seeds from the species' entire distribution were collected from 12 co-occurring populations. Germination thermal niche was investigated using a thermal gradient plate (5-40°C). Five constant temperatures between 9 and 33°C were used to test how the germination niche (1) differs between species, (2) varies among populations, and (3) relates to the climate of origin. Germination response differed among species; jarrah had a lower optimal temperature and thermal limit than marri (T o 15.3°C, 21.2°C; ED50 23.4°C, 31°C, respectively). The thermal limit for germination differed among populations within both species, yet only marri showed evidence for adaptation to thermal origins. While marri has the capacity for germination at higher thermal temperatures, jarrah is more vulnerable to global warming exceeding safety margins. This discrepancy is predicted to alter species distributions and forest composition in the future.

Aridity drives clinal patterns in leaf traits and responsiveness to precipitation in a broadly distributed Australian tree species.

Aridity shapes species distributions and plant growth and function worldwide. Yet, plant traits often show complex relationships with aridity, challenging our understanding of aridity as a driver of evolutionary adaptation. We grew nine genotypes of Eucalyptus camaldulensis subsp. camaldulensis sourced from an aridity gradient together in the field for ~650 days under low and high precipitation treatments. Eucalyptus camaldulesis is considered a phreatophyte (deep-rooted species that utilizes groundwater), so we hypothesized that genotypes from more arid environments would show lower aboveground productivity, higher leaf gas-exchange rates, and greater tolerance/avoidance of dry surface soils (indicated by lower responsiveness) than genotypes from less arid environments. Aridity predicted genotype responses to precipitation, with more arid genotypes showing lower responsiveness to reduced precipitation and dry surface conditions than less arid genotypes. Under low precipitation, genotype net photosynthesis and stomatal conductance increased with home-climate aridity. Across treatments, genotype intrinsic water-use efficiency and osmotic potential declined with increasing aridity while photosynthetic capacity (Rubisco carboxylation and RuBP regeneration) increased with aridity. The observed clinal patterns indicate that E. camaldulensis genotypes from extremely arid environments possess a unique strategy defined by lower responsiveness to dry surface soils, low water-use efficiency, and high photosynthetic capacity. This strategy could be underpinned by deep rooting and could be adaptive under arid conditions where heat avoidance is critical and water demand is high.

Contributions of phenotypic integration, plasticity and genetic adaptation to adaptive capacity relating to drought in Banksia marginata (Proteaceae).

The frequency and intensity of drought events are predicted to increase because of climate change, threatening biodiversity and terrestrial ecosystems in many parts of the world. Drought has already led to declines in functionally important tree species, which are documented in dieback events, shifts in species distributions, local extinctions, and compromised ecosystem function. Understanding whether tree species possess the capacity to adapt to future drought conditions is a major conservation challenge. In this study, we assess the capacity of a functionally important plant species from south-eastern Australia (Banksia marginata, Proteaceae) to adapt to water-limited environments. A water-manipulated common garden experiment was used to test for phenotypic plasticity and genetic adaptation in seedlings sourced from seven provenances of contrasting climate-origins (wet and dry). We found evidence of local adaptation relating to plant growth investment strategies with populations from drier climate-origins showing greater growth in well-watered conditions. The results also revealed that environment drives variation in physiological (stomatal conductance, predawn and midday water potential) and structural traits (wood density, leaf dry matter content). Finally, these results indicate that traits are coordinated to optimize conservation of water under water-limited conditions and that trait coordination (phenotypic integration) does not constrain phenotypic plasticity. Overall, this study provides evidence for adaptive capacity relating to drought conditions in B. marginata, and a basis for predicting the response to climate change in this functionally important plant species.

Soil chemistry and fungal communities are associated with dieback in an Endangered Australian shrub.

Background and aims: Field surveys across known populations of the Endangered Persoonia hirsuta (Proteaceae) in 2019 suggested the soil environment may be associated with dieback in this species. To explore how characteristics of the soil environment (e.g., pathogens, nutrients, soil microbes) relate to dieback, a soil bioassay (Experiment 1) was conducted using field soils from two dieback effected P. hirsuta populations. Additionally, a nitrogen addition experiment (Experiment 2) was conducted to explore how the addition of soil nitrogen impacts dieback. Methods: The field soils were baited for pathogens, and soil physiochemical and microbial community characteristics were assessed and related to dieback among plants in the field and nursery-grown plants inoculated with the same field soils. Roots from inoculated plants were harvested to confirm the presence of soil pathogens and root-associated endophytes. Using these isolates, a dual culture antagonism assay was performed to examine competition among these microbes and identify candidate pathogens or pathogen antagonists. Results: Dieback among plants in the field and Experiment 1 was associated with soil physiochemical properties (nitrogen and potassium), and soil microbes were identified as significant indicators of healthy and dieback-affected plants. Plants in Experiment 2 exhibited greater dieback when treated with elevated nitrogen. Additionally, post-harvest culturing identified fungi and other soil pathogens, some of which exhibited antagonistic behavior. Conclusion: This study identified candidate fungi and soil physiochemical properties associated with observed dieback and dieback resistance in an Endangered shrub and provides groundwork for further exploring what drives dieback and how it can be managed to promote the conservation of wild populations. Supplementary Information: The online version contains supplementary material available at 10.1007/s11104-022-05724-7.

The roles of divergent and parallel molecular evolution contributing to thermal adaptive strategies in trees.

Local adaptation is a driver of biological diversity, and species may develop analogous (parallel evolution) or alternative (divergent evolution) solutions to similar ecological challenges. We expect these adaptive solutions would culminate in both phenotypic and genotypic signals. Using two Eucalyptus species (Eucalyptus grandis and Eucalyptus tereticornis) with overlapping distributions grown under contrasting 'local' temperature conditions to investigate the independent contribution of adaptation and plasticity at molecular, physiological and morphological levels. The link between gene expression and traits markedly differed between species. Divergent evolution was the dominant pattern driving adaptation (91% of all significant genes); but overlapping gene (homologous) responses were dependent on the determining factor (plastic, adaptive or genotype by environment interaction). Ninety-eight percent of the plastic homologs were similarly regulated, while 50% of the adaptive homologs and 100% of the interaction homologs were antagonistical. Parallel evolution for the adaptive effect in homologous genes was greater than expected but not in favour of divergent evolution. Heat shock proteins for E. grandis were almost entirely driven by adaptation, and plasticity in E. tereticornis. These results suggest divergent molecular evolutionary solutions dominated the adaptive mechanisms among species, even in similar ecological circumstances. Suggesting that tree species with overlapping distributions are unlikely to equally persist in the future.

Signatures of natural selection in a foundation tree along Mediterranean climatic gradients.

Temperature and precipitation regimes are rapidly changing, resulting in forest dieback and extinction events, particularly in Mediterranean-type climates (MTC). Forest management that enhance forests' resilience is urgently required, however adaptation to climates in heterogeneous landscapes with multiple selection pressures is complex. For widespread trees in MTC 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. We sampled 28 populations across the geographic distribution of Eucalyptus marginata (jarrah), in South-west Western Australia, and obtained 13,534 independent single nucleotide polymorphic (SNP) markers across the genome. Three genotype-association analyses that employ different ways of correcting population structure were used to identify putatively adapted SNPs associated with independent climate variables. While overall levels of population differentiation were low (FST  = 0.04), environmental association analyses found a total of 2336 unique SNPs associated with temperature and precipitation variables, with 1440 SNPs annotated to genic regions. Considerable allelic turnover was identified for SNPs associated with temperature seasonality and mean precipitation of the warmest quarter, suggesting that both temperature and precipitation are important factors in adaptation. SNPs with similar gene functions had analogous allelic turnover along climate gradients, while SNPs among temperature and precipitation variables had uncorrelated patterns of adaptation. These contrasting patterns provide evidence that there may be standing genomic variation adapted to current climate gradients, providing the basis for adaptive management strategies to bolster forest resilience in the future.

Major biogeographic barriers in eastern Australia have shaped the population structure of widely distributed Eucalyptus moluccana and its putative subspecies.

We have investigated the impact of recognized biogeographic barriers on genetic differentiation of grey box (Eucalyptus moluccana), a common and widespread tree species of the family Myrtaceae in eastern Australian woodlands, and its previously proposed four subspecies moluccana, pedicellata, queenslandica, and crassifolia. A range of phylogeographic analyses were conducted to examine the population genetic differentiation and subspecies genetic structure in E. moluccana in relation to biogeographic barriers. Slow evolving markers uncovering long term processes (chloroplast DNA) were used to generate a haplotype network and infer phylogeographic barriers. Additionally, fast evolving, hypervariable markers (microsatellites) were used to estimate demographic processes and genetic structure among five geographic regions (29 populations) across the entire distribution of E. moluccana. Morphological features of seedlings, such as leaf and stem traits, were assessed to evaluate population clusters and test differentiation of the putative subspecies. Haplotype network analysis revealed twenty chloroplast haplotypes with a main haplotype in a central position shared by individuals belonging to the regions containing the four putative subspecies. Microsatellite analysis detected the genetic structure between Queensland (QLD) and New South Wales (NSW) populations, consistent with the McPherson Range barrier, an east-west spur of the Great Dividing Range. The substructure was detected within QLD and NSW in line with other barriers in eastern Australia. The morphological analyses supported differentiation between QLD and NSW populations, with no difference within QLD, yet some differentiation within NSW populations. Our molecular and morphological analyses provide evidence that several geographic barriers in eastern Australia, including the Burdekin Gap and the McPherson Range have contributed to the genetic structure of E. moluccana. Genetic differentiation among E. moluccana populations supports the recognition of some but not all the four previously proposed subspecies, with crassifolia being the most differentiated.

Regarding the F-word: The effects of data filtering on inferred genotype-environment associations.

Genotype-environment association (GEA) methods have become part of the standard landscape genomics toolkit, yet, we know little about how to best filter genotype-by-sequencing data to provide robust inferences for environmental adaptation. In many cases, default filtering thresholds for minor allele frequency and missing data are applied regardless of sample size, having unknown impacts on the results, negatively affecting management strategies. Here, we investigate the effects of filtering on GEA results and the potential implications for assessment of adaptation to environment. We use empirical and simulated data sets derived from two widespread tree species to assess the effects of filtering on GEA outputs. Critically, we find that the level of filtering of missing data and minor allele frequency affect the identification of true positives. Even slight adjustments to these thresholds can change the rate of true positive detection. Using conservative thresholds for missing data and minor allele frequency substantially reduces the size of the data set, lessening the power to detect adaptive variants (i.e., simulated true positives) with strong and weak strengths of selection. Regardless, strength of selection was a good predictor for GEA detection, but even some SNPs under strong selection went undetected. False positive rates varied depending on the species and GEA method, and filtering significantly impacted the predictions of adaptive capacity in downstream analyses. We make several recommendations regarding filtering for GEA methods. Ultimately, there is no filtering panacea, but some choices are better than others, depending on the study system, availability of genomic resources, and desired objectives.

Vulnerability to xylem cavitation of Hakea species (Proteaceae) from a range of biomes and life histories predicted by climatic niche.

BACKGROUND AND AIMS: Extreme drought conditions across the globe are impacting biodiversity, with serious implications for the persistence of native species. However, quantitative data on physiological tolerance are not available for diverse flora to inform conservation management. We quantified physiological resistance to cavitation in the diverse Hakea genus (Proteaceae) to test predictions based on climatic origin, life history and functional traits. METHODS: We sampled terminal branches of replicate plants of 16 species in a common garden. Xylem cavitation was induced in branches under varying water potentials (tension) in a centrifuge, and the tension generating 50 % loss of conductivity (stem P50) was characterized as a metric for cavitation resistance. 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). KEY RESULTS: There was significant variation in stem P50 among species, which was negatively associated with the species climate origin (rainfall and aridity). Cavitation resistance 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. CONCLUSIONS: This study highlights climate (rainfall and aridity), rather than life history and functional traits, as the key predictor of variation in cavitation resistance (stem P50). Rainfall for species origin was the best predictor of cavitation resistance, explaining variation in stem P50, which appears to be a major determinant of species distribution. This study also indicates that stem P50 is an adaptive trait, genetically determined, and hence reliable and robust for predicting species vulnerability to climate change. Our findings will contribute to future prediction of species vulnerability to drought and adaptive management under climate change.

Spatial, climate and ploidy factors drive genomic diversity and resilience in the widespread grass Themeda triandra.

Global climate change poses a significant threat to natural communities around the world, with many plant species showing signs of climate stress. Grassland ecosystems are not an exception, with climate change compounding contemporary pressures such as habitat loss and fragmentation. In this study, we assess the climate resilience of Themeda triandra, a foundational species and the most widespread plant in Australia, by assessing the relative contributions of spatial, environmental and ploidy factors to contemporary genomic variation. Reduced-representation genome sequencing on 472 samples from 52 locations was used to test how the distribution of genomic variation, including ploidy polymorphism, supports adaptation to hotter and drier climates. We explicitly quantified isolation by distance (IBD) and isolation by environment (IBE) and predicted genomic vulnerability of populations to future climates based on expected deviation from current genomic composition. We found that a majority (54%) of genomic variation could be attributed to IBD, while an additional 22% (27% when including ploidy information) could be explained by two temperature and two precipitation climate variables demonstrating IBE. Ploidy polymorphisms were common within populations (31/52 populations), indicating that ploidy mixing is characteristic of T. triandra populations. Genomic vulnerabilities were found to be heterogeneously distributed throughout the landscape, and our analysis suggested that ploidy polymorphism, along with other factors linked to polyploidy, reduced vulnerability to future climates by 60% (0.25-0.10). Our data suggests that polyploidy may facilitate adaptation to hotter climates and highlight the importance of incorporating ploidy in adaptive management strategies to promote the resilience of this and other foundation species.

Leaf trait variation is similar among genotypes of Eucalyptus camaldulensis from differing climates and arises in plastic responses to the seasons rather than water availability.

We used a widely distributed tree Eucalyptus camaldulensis subsp. camaldulensis to partition intraspecific variation in leaf functional traits to genotypic variation and phenotypic plasticity. We examined if genotypic variation is related to the climate of genotype provenance and whether phenotypic plasticity maintains performance in a changing environment. Ten genotypes from different climates were grown in a common garden under watering treatments reproducing the wettest and driest edges of the subspecies' distribution. We measured functional traits reflecting leaf metabolism and associated with growth (respiration rate, nitrogen and phosphorus concentrations, and leaf mass per area) and performance proxies (aboveground biomass and growth rate) each season over a year. Genotypic variation contributed substantially to the variation in aboveground biomass but much less in growth rate and leaf traits. Phenotypic plasticity was a large source of the variation in leaf traits and performance proxies and was greater among sampling dates than between watering treatments. The variation in leaf traits was weakly correlated to performance proxies, and both were unrelated to the climate of genotype provenance. Intraspecific variation in leaf traits arises similarly among genotypes in response to seasonal environmental variation, instead of long-term water availability or climate of genotype provenance.

The fate of carbon in a mature forest under carbon dioxide enrichment.

Atmospheric carbon dioxide enrichment (eCO2) can enhance plant carbon uptake and growth1-5, thereby providing an important negative feedback to climate change by slowing the rate of increase of the atmospheric CO2 concentration6. Although evidence gathered from young aggrading forests has generally indicated a strong CO2 fertilization effect on biomass growth3-5, it is unclear whether mature forests respond to eCO2 in a similar way. In mature trees and forest stands7-10, photosynthetic uptake has been found to increase under eCO2 without any apparent accompanying growth response, leaving the fate of additional carbon fixed under eCO2 unclear4,5,7-11. Here using data from the first ecosystem-scale Free-Air CO2 Enrichment (FACE) experiment in a mature forest, we constructed a comprehensive ecosystem carbon budget to track the fate of carbon as the forest responded to four years of eCO2 exposure. We show that, although the eCO2 treatment of +150 parts per million (+38 per cent) above ambient levels induced a 12 per cent (+247 grams of carbon per square metre per year) increase in carbon uptake through gross primary production, this additional carbon uptake did not lead to increased carbon sequestration at the ecosystem level. Instead, the majority of the extra carbon was emitted back into the atmosphere via several respiratory fluxes, with increased soil respiration alone accounting for half of the total uptake surplus. Our results call into question the predominant thinking that the capacity of forests to act as carbon sinks will be generally enhanced under eCO2, and challenge the efficacy of climate mitigation strategies that rely on ubiquitous CO2 fertilization as a driver of increased carbon sinks in global forests.

Plant functional traits differ in adaptability and are predicted to be differentially affected by climate change.

Climate change is testing the resilience of forests worldwide pushing physiological tolerance to climatic extremes. Plant functional traits have been shown to be adapted to climate and have evolved patterns of trait correlations (similar patterns of distribution) and coordinations (mechanistic trade-off). We predicted that traits would differentiate between populations associated with climatic gradients, suggestive of adaptive variation, and correlated traits would adapt to future climate scenarios in similar ways.We measured genetically determined trait variation and described patterns of correlation for seven traits: photochemical reflectance index (PRI), normalized difference vegetation index (NDVI), leaf size (LS), specific leaf area (SLA), δ13C (integrated water-use efficiency, WUE), nitrogen concentration (NCONC), and wood density (WD). All measures were conducted in an experimental plantation on 960 trees sourced from 12 populations of a key forest canopy species in southwestern Australia.Significant differences were found between populations for all traits. Narrow-sense heritability was significant for five traits (0.15-0.21), indicating that natural selection can drive differentiation; however, SLA (0.08) and PRI (0.11) were not significantly heritable. Generalized additive models predicted trait values across the landscape for current and future climatic conditions (>90% variance). The percent change differed markedly among traits between current and future predictions (differing as little as 1.5% (δ13C) or as much as 30% (PRI)). Some trait correlations were predicted to break down in the future (SLA:NCONC, δ13C:PRI, and NCONC:WD).Synthesis: Our results suggest that traits have contrasting genotypic patterns and will be subjected to different climate selection pressures, which may lower the working optimum for functional traits. Further, traits are independently associated with different climate factors, indicating that some trait correlations may be disrupted in the future. Genetic constraints and trait correlations may limit the ability for functional traits to adapt to climate change.

Adaptive variation for growth and resistance to a novel pathogen along climatic gradients in a foundation tree.

Natural ecosystems are under pressure from increasing abiotic and biotic stressors, including climate change and novel pathogens, which are putting species at risk of local extinction, and altering community structure, composition and function. Here, we aim to assess adaptive variation in growth and fungal disease resistance within a foundation tree, Corymbia calophylla to determine local adaptation, trait heritability and genetic constraints in adapting to future environments. Two experimental planting sites were established in regions of contrasting rainfall with seed families from 18 populations capturing a wide range of climate origins (~4,000 individuals at each site). Every individual was measured in 2015 and 2016 for growth (height, basal diameter) and disease resistance to a recently introduced leaf blight pathogen (Quambalaria pitereka). Narrow-sense heritability was estimated along with trait covariation. Trait variation was regressed against climate-of-origin, and multivariate models were used to develop predictive maps of growth and disease resistance. Growth and blight resistance traits differed significantly among populations, and these differences were consistent between experimental sites and sampling years. Growth and blight resistance were heritable, and comparisons between trait differentiation (Q ST) and genetic differentiation (F ST) revealed that population differences in height and blight resistance traits are due to divergent natural selection. Traits were significantly correlated with climate-of-origin, with cool and wet populations showing the highest levels of growth and blight resistance. These results provide evidence that plants have adaptive growth strategies and pathogen defence strategies. Indeed, the presence of standing genetic variation and trait heritability of growth and blight resistance provide capacity to respond to novel, external pressures. The integration of genetic variation into adaptive management strategies, such as assisted gene migration and seed sourcing, may be used to provide greater resilience for natural ecosystems to both biotic and abiotic stressors.

Desiccation time during drought is highly predictable across species of Eucalyptus from contrasting climates.

Catastrophic failure of the water transport pathway in trees is a principal mechanism of mortality during extreme drought. To be able to predict the probability of mortality at an individual and landscape scale we need knowledge of the time for plants to reach critical levels of hydraulic failure. We grew plants of eight species of Eucalyptus originating from contrasting climates before allowing a subset to dehydrate. We tested whether a trait-based model of time to plant desiccation tcrit , from stomatal closure gs90 to a critical level of hydraulic dysfunction Ψcrit is consistent with observed dry-down times. Plant desiccation time varied among species, ranging from 96.2 to 332 h at a vapour-pressure deficit of 1 kPa, and was highly predictable using the tcrit model in conjunction with a leaf shedding function. Plant desiccation time was longest in species with high cavitation resistance, strong vulnerability segmentation, wide stomatal-hydraulic safety, and a high ratio of total plant water content to leaf area. Knowledge of tcrit in combination with water-use traits that influence stomatal closure could significantly increase our ability to predict the timing of drought-induced mortality at tree and forest scales.

Standing genomic variation within coding and regulatory regions contributes to the adaptive capacity to climate in a foundation tree species.

Global climate is rapidly changing, and the ability for tree species to adapt is dependent on standing genomic variation; however, the distribution and abundance of functional and adaptive variants are poorly understood in natural systems. We test key hypotheses regarding the genetics of adaptive variation in a foundation tree: genomic variation is associated with climate, and genomic variation is more likely to be associated with temperature than precipitation or aridity. To test these hypotheses, we used 9,593 independent, genomic single-nucleotide polymorphisms (SNPs) from 270 individuals sampled from Corymbia calophylla's entire distribution in south-western Western Australia, spanning orthogonal temperature and precipitation gradients. Environmental association analyses returned 537 unique SNPs putatively adaptive to climate. We identified SNPs associated with climatic variation (i.e., temperature [458], precipitation [75] and aridity [78]) across the landscape. Of these, 78 SNPs were nonsynonymous (NS), while 26 SNPs were found within gene regulatory regions. The NS and regulatory candidate SNPs associated with temperature explained more deviance (27.35%) than precipitation (5.93%) and aridity (4.77%), suggesting that temperature provides stronger adaptive signals than precipitation. Genes associated with adaptive variants include functions important in stress responses to temperature and precipitation. Patterns of allelic turnover of NS and regulatory SNPs show small patterns of change through climate space with the exception of an aldehyde dehydrogenase gene variant with 80% allelic turnover with temperature. Together, these findings provide evidence for the presence of adaptive variation to climate in a foundation species and provide critical information to guide adaptive management practices.

Drought tolerance traits do not vary across sites differing in water availability in Banksia serrata (Proteaceae).

Interspecific variation in plant hydraulic traits plays a major role in shaping species distributions across climates, yet variation within species is poorly understood. Here we report on intraspecific variation of hydraulic traits in Banksia serrata (L.f.) sampled from three sites characterised by contrasting climates (warm-wet, warm-dry and cool-wet). Hydraulic characteristics including vulnerability to embolism, hydraulic conductance, pressure-volume traits and key morphological traits were measured. Vulnerability to embolism in leaf and stem, defined by the water potential inducing 50 and 88% loss of hydraulic conductivity (P50 and P88 respectively), did not differ across sites. However, plants from the warm-dry environment exhibited higher stem conductivity (Ks) than the cool-wet environment. Leaf turgor loss point (TLP) did not vary among sites, but warm-dry site plants showed lower leaf capacitance (C*FT) and higher modulus of elasticity (ε) than the other two sites. Plants from the cool-wet site had lower specific leaf area (SLA) and plants from the warm-dry site had lower sapwood density (WD). Overall, key hydraulic traits were generally conserved across populations despite differences in mean site water availability, and the safety-efficiency trade-off was absent in this species. These results suggest that B. serrata has limited ability to adjust hydraulic architecture in response to environmental change and thus may be susceptible to climate change-type drought stress.