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

Nonlinear effects of temperature on mosquito parasite infection across a large geographic climate gradient

Temperature drives ectothermic host - parasite interactions, making them particularly sensitive to climatic variation and change. To isolate the role of temperature, lab-based studies are increasingly used to assess and forecast disease risk under current and future climate conditions. However, in the field, the effects of temperature on parasitism may be mediated by other sources of variation, including local adaptation. To address the key knowledge gaps of how temperature influences host - parasite interactions and whether thermal responses measured in controlled experiments capture infection across temperature gradients in nature, we paired an extensive field survey of parasitism-by the ciliate Lambornella clarki on its tree hole mosquito host, Aedes sierrensis -with laboratory experiments describing parasitism thermal performance curves (TPCs) for six host populations from varying climates. We also investigated the mechanisms underlying the thermal biology of the host - parasite interaction by separately measuring TPCs for infection, host immunity, and parasite growth rates. Along the west coast of North America, across an 1100 km climate gradient spanning 12°C mean rainy season temperature variation, we found that parasitism peaked at intermediately cold temperatures, and was consistent both between field seasons and with the lab experiment results. The experiments produced no evidence of host intraspecific variation in temperature sensitivity to parasitism. Importantly, parasitism peaked at temperatures below the thermal optimum for free-living L. clarki due to the balance of temperature effects on parasite growth and reproduction against the strength of the host melanization immune response. The results suggest that nonlinear responses to temperature drive parasitism in nature, and that simple lab and field studies can accurately capture the thermal biology of multilayered host - parasite interactions. Data and code for this submission are provided on Dryad: http://datadryad.org/stash/share/CfZkk4LsJzljetJJnFZMDMrjuciTXMxrkrc95I2J3tA .

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.

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.

Adaptive plasticity in plant traits increases time to hydraulic failure under drought in a foundation tree

The viability of forest trees, in response to climate change-associated drought, will depend on their capacity to survive through genetic adaptation and phenotypic plasticity in drought tolerance traits. Genotypes with enhanced plasticity for drought tolerance (adaptive plasticity) will have a greater ability to persist and delay the onset of hydraulic failure. By examining populations from different climate-origins grown under contrasting soil water availability, we tested for genotype (G), environment (E) and genotype-by-environment (G × E) effects on traits that determine the time it takes for saplings to desiccate from stomatal closure to 88% loss of stem hydraulic conductance (time to hydraulic failure, THF). Specifically, we hypothesized that: (i) THF is dependent on a G × E interaction, with longer THF for warm, dry climate populations in response to chronic water deficit treatment compared with cool, wet populations, and (ii) hydraulic and allometric traits explain the observed patterns in THF. Corymbia calophylla saplings from two populations originating from contrasting climates (warm-dry or cool-wet) were grown under well-watered and chronic soil water deficit treatments in large containers. Hydraulic and allometric traits were measured and then saplings were dried-down to critical levels of drought stress to estimate THF. Significant plasticity was detected in the warm-dry population in response to water-deficit, with enhanced drought tolerance compared with the cool-wet population. Projected leaf area and total plant water storage showed treatment variation, and minimum conductance showed significant population differences driving longer THF in trees from warm-dry origins grown in water-limited conditions. Our findings contribute information on intraspecific variation in key drought traits, including hydraulic and allometric determinants of THF. It highlights the need to quantify adaptive capacity in populations of forest trees in climate change-type drought to improve predictions of forest die-back.

Expansion of Ash Dieback towards the scattered Fraxinus excelsior range of the Italian peninsula

Hymenoscyphus fraxineus, causal agent of Ash Dieback, has posed a threat to Fraxinus excelsior (common ash) in Europe since the 1990s. In south-western Europe, optimal climatic conditions for H. fraxineus become scattered and host density decreases, reducing disease spread rates. To date, the Ash Dieback agent has not been reported from southern and most of central Italy, where native F. excelsior is present as small fragmented populations. This study examines the expansion of Ash Dieback into central Italy, and it considers the consequences of further local spread with regards to the loss of F. excelsior genetic resource. Symptomatic F. excelsior were sampled from sixteen sites in northern and central Italy during 2020. Specimens were analyzed with a culturomics and a quantitative PCR approach. A bibliographic search of F. excelsior floristic reports was conducted for the creation of a detailed range map. The combined use of both techniques confirmed the presence of H. fraxineus in all the sites of central Italy where host plants were symptomatic. These new records represent the southern limit of the current known distribution of this pathogen in Italy, and together with Montenegro, in Europe. The characterization of the F. excelsior scattered range suggests that further spread of Ash Dieback across southern Italy is a realistic scenario. This presents a threat not just to the southern European proveniences of F. excelsior, but to the species as a whole, should Ash Dieback lead to the loss of warm climate adapted genetic material, which may become increasingly valuable under climate change.

Repeated extreme heatwaves result in higher leaf thermal tolerances and greater safety margins

The frequency and severity of heatwave events are increasing, exposing species to conditions beyond their physiological limits. Species respond to heatwaves in different ways, however it remains unclear if plants have the adaptive capacity to successfully respond to hotter and more frequent heatwaves. We exposed eight tree populations from two climate regions grown under cool and warm temperatures to repeated heatwave events of moderate (40°C) and extreme (46°C) severity to assess adaptive capacity to heatwaves. Leaf damage and maximum quantum efficiency of photosystem II (F_v /F_m ) were significantly impacted by heatwave severity and growth temperatures, respectively; populations from a warm-origin avoided damage under moderate heatwaves compared to those from a cool-origin, indicating a degree of local adaptation. We found that plasticity to heatwave severity and repeated heatwaves contributed to enhanced thermal tolerance and lower leaf temperatures_, leading to greater thermal safety margins (thermal tolerance minus leaf temperature) in a second heatwave. Notably, while we show that adaptation and physiological plasticity are important factors affecting plant adaptive capacity to thermal stress, plasticity of thermal tolerances and thermal safety margins provides the opportunity for trees to persist among fluctuating heatwave exposures.

Aridity drives coordinated trait shifts but not decreased trait variance across the geographic range of eight Australian trees

Large intraspecific functional trait variation strongly impacts many aspects of communities and ecosystems, and is the medium upon which evolution works. Yet intraspecific trait variation is inconsistent and hard to predict across traits, species and locations. We measured within-species variation in leaf mass per area (LMA), leaf dry matter content (LDMC), branch wood density (WD), and allocation to stem area vs leaf area in branches (branch Huber value (HV)) across the aridity range of seven Australian eucalypts and a co-occurring Acacia species to explore how traits and their variances change with aridity. Within species, we found consistent increases in LMA, LDMC and WD and HV with increasing aridity, resulting in consistent trait coordination across leaves and branches. However, this coordination only emerged across sites with large climate differences. Unlike trait means, patterns of trait variance with aridity were mixed across populations and species. Only LDMC showed constrained trait variation in more xeric species and drier populations that could indicate limits to plasticity or heritable trait variation. Our results highlight that climate can drive consistent within-species trait patterns, but that patterns might often be obscured by the complex nature of morphological traits, sampling incomplete species ranges or sampling confounded stress gradients.

Assessing adaptive and plastic responses in growth and functional traits in a 10-year-old common garden experiment with pedunculate oak (Quercus robur L.) suggests that directional selection can drive climatic adaptation

Understanding how tree species will respond to a future climate requires reliable and quantitative estimates of intra-specific variation under current climate conditions. We studied three 10-year-old common garden experiments established across a rainfall and drought gradient planted with nearly 10,000 pedunculate oak (Quercus robur L.) trees from ten provenances with known family structure. We aimed at disentangling adaptive and plastic responses for growth (height and diameter at breast height) as well as for leaf and wood functional traits related to adaptation to dry environments. We used restricted maximum likelihood approaches to assess additive genetic variation expressed as narrow-sense heritability (h2), quantitative trait differentiation among provenances (QST), and genotype-by-environment interactions (GxE). We found strong and significant patterns of local adaptation in growth in all three common gardens, suggesting that transfer of seed material should not exceed a climatic distance of approximately 1°C under current climatic conditions, while transfer along precipitation gradients seems to be less stringent. Moreover, heritability reached 0.64 for tree height and 0.67 for dbh at the dry margin of the testing spectrum, suggesting significant additive genetic variation of potential use for future selection and tree breeding. GxE interactions in growth were significant and explained less phenotypic variation than origin of seed source (4% versus 10%). Functional trait variation among provenances was partly related to drought regimes at provenances origins but had moderate explanatory power for growth. We conclude that directional selection, either naturally or through breeding, is the most likely and feasible outcome for pedunculate oak to adapt to warmer and drier climate conditions in the future.

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.

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.