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.

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 (F_ST = 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.

As old as the hills: Pliocene palaeogeographical processes influence patterns of genetic structure in the widespread, common shrub Banksia sessilis

The impact of Quaternary glaciation on the development of phylogeographic structure in plant species is well documented. In unglaciated landscapes, phylogeographic patterns tend to reflect processes relating to persistence and stochasticity, yet other factors, associated with the palaeogeographical history of the landscape, including geomorphological events, can also have a significant influence. The unglaciated landscape of south-western Western Australia is an ideal location to observe these ancient drivers of lineage diversification, with tectonic activity associated with the Darling Fault in the late Pliocene attributed to patterns of deep phylogeographic divergence in a widespread tree from this region. Interestingly, other species within this region have not shown this pattern and this palaeogeographical boundary therefore presents an opportunity to examine age and historical distribution of plant species endemic to this region. In this study, we assess patterns of genetic diversity and structure across 28 populations of the widespread shrub Banksia sessilis using three cpDNA markers and nine nuclear microsatellite markers. Sixteen cpDNA haplotypes were identified, comprising two major chloroplast DNA lineages that are estimated to have diverged in the Pliocene, approximately 3.3 million years ago. This timing coincides with major geomorphological processes in the landscape, including the separation of the Darling Plateau from the adjacent Swan Coastal Plain, as well as eustatic changes on the Swan Coastal Plain that are likely to have resulted in the physical isolation of historical plant lineages. Chloroplast lineages were broadly aligned with populations associated with older lateritic soils of the Darling Plateau and Geraldton sandplains or the younger sandy soils associated with the Swan Coastal Plain and Southern Coastline. This structural pattern of lateritic versus non-lateritic division was not observed in the nuclear microsatellite data that identified three genetic clades that roughly corresponded to populations in the North, South, and Central portions of the distributions.

Differential exposure and susceptibility to threats based on evolutionary history: how OCBIL theory informs flora conservation

Optimal conservation approaches have been proposed to differ for biota with contrasting evolutionary histories. Natural selection filters the distribution of plant traits over evolutionary time, with the current expression of traits mediating susceptibility to contemporary and often novel threats. We use old, climatically buffered, infertile landscape (OCBIL) theory to compile predictions regarding differences in exposure and susceptibility to key threats between OCBIL and young, often disturbed, fertile landscape (YODFEL) flora. Based on literature and existing data from the Southwest Australian Floristic Region (SWAFR), we evaluate evidence in support of our predictions, finding strong theoretical and empirical support for the proposition that exposure and/or impact of many threats differs between OCBILs and YODFELs. OCBILs have more exposure to land clearance from mining, whereas many YODFELs have greater exposure to land clearance from agriculture, and urban and industrial land uses, and greater overall levels of habitat loss and fragmentation. OCBIL flora are more susceptible to pathogens and extremes of fire interval than YODFEL flora, but conversely may have a greater capacity to persist in smaller populations if small populations featured in the evolutionary history of the species prior to anthropogenic fragmentation, and have substantial resistance to weed invasion. We argue that consideration of evolutionary history has an important role in informing conservation management.

Extensive Genetic Connectivity and Historical Persistence Are Features of Two Widespread Tree Species in the Ancient Pilbara Region of Western Australia

Phylogeographic studies can be used as a tool to understand the evolutionary history of a landscape, including the major drivers of species distributions and diversity. Extensive research has been conducted on phylogeographic patterns of species found in northern hemisphere landscapes that were affected by glaciations, yet the body of literature for older, unaffected landscapes is still underrepresented. The Pilbara region of north-western Australia is an ancient and vast landscape that is topographically complex, consisting of plateaus, gorges, valleys, and ranges, and experiences extreme meteorological phenomena including seasonal cyclonic activity. These features are expected to influence patterns of genetic structuring throughout the landscape either by promoting or restricting the movement of pollen and seed. Whilst a growing body of literature exists for the fauna endemic to this region, less is known about the forces shaping the evolution of plant taxa. In this study we investigate the phylogeography of two iconic Pilbara tree species, the Hamersley Bloodwood (Corymbia hamersleyana) and Western Gidgee (Acacia pruinocarpa), by assessing patterns of variation and structure in several chloroplast DNA regions and nuclear microsatellite loci developed for each species. Gene flow was found to be extensive in both taxa and there was evidence of long-distance seed dispersal across the region (pollen to seed ratios of 6.67 and 2.96 for C. hamersleyana and A. pruinocarpa, respectively), which may result from flooding and strong wind gusts associated with extreme cyclonic activity. Both species possessed high levels of cpDNA genetic diversity in comparison to those from formerly glaciated landscapes (C. hamersleyana = 14 haplotypes, A. pruinocarpa = 37 haplotypes) and showed evidence of deep lineage diversification occurring from the late Miocene, a time of intensifying aridity in this landscape that appears to be a critical driver of evolution in Pilbara taxa. In contrast to another study, we did not find evidence for topographic features acting as refugia for the widely sampled C. hamersleyana.

Platysace (Apiaceae) of south-western Australia: silent story tellers of an ancient human landscape

High gene flow and a population structure corresponding to human rather than geographical drivers are likely to be genetic patterns of human-dispersed plant taxa. We examined variation in geographical structure and gene flow estimates based on three non-coding regions of plastid DNA in three south-west Australian members of the Platysace genus to identify whether a human influence on dispersion of utilized taxa was detectable. Edible tubers of Platysace deflexa and Platysace trachymenioides have been harvested historically by Noongar traditional owners, whereas Platysace effusa has no known cultural significance. We found differences between utilized and non-utilized taxa, particularly when considered against the generally complex phylogeographical patterning in south-west Australian plant taxa. Platysace effusa showed a pattern of high population divergence, low gene flow and multiple refugia, consistent with a long evolutionary history, past climatic oscillations and persistence in a highly fragmented landscape. In contrast, higher gene flow estimates, less divergence between populations and common haplotypes in P. deflexa and in P. trachymenioides over the south-eastern part of its range are consistent with anthropogenic influences. This study contributes to the understanding of human influences on south-west Australian plant taxa that have been present since the late Pleistocene, but to date have received little scientific attention.

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.

Pervasive admixture between eucalypt species has consequences for conservation and assisted migration

Conservation management often uses information on genetic population structure to assess the importance of local provenancing for ecological restoration and reintroduction programs. For species that do not exhibit complete reproductive isolation, the estimation of population genetic parameters may be influenced by the extent of admixture. Therefore, to avoid perverse outcomes for conservation, genetically informed management strategies must determine whether hybridization between species is relevant, and the extent to which observed population genetic patterns are shaped by interspecific versus intraspecific gene flow. We used genotyping by sequencing to identify over 2,400 informative single nucleotide polymorphisms across 18 populations of Eucalyptus regnans F. Muell., a foundation tree species of montane forests in south-eastern Australia. We used these data to determine the extent of hybridization with another species, Eucalyptus obliqua L'Hér., and investigate how admixture influences genetic diversity parameters, by estimating metrics of genetic diversity and examining population genetic structure in datasets with and without admixed individuals. We found hybrid individuals at all sites and two highly introgressed populations. Hybrid individuals were not distributed evenly across environmental gradients, with logistic regression identifying hybrids as being associated with temperature. Removal of hybrids resulted in increases in genetic differentiation (F ST), expected heterozygosity, observed heterozygosity and the inbreeding coefficient, and different patterns of isolation by distance. After removal of hybrids and introgressed populations, mountain ash showed very little population genetic structure, with a small effect of isolation by distance, and very low global F ST(0.03). Our study shows that, in plants, decisions around provenancing of individuals for restoration depend on knowledge of whether hybridization is influencing population genetic structure. For species in which most genetic variation is held within populations, there may be little benefit in planning conservation strategies around environmental adaptation of seed sources. The possibility for adaptive introgression may also be relevant when species regularly hybridize.

Landscape genomic prediction for restoration of a Eucalyptus foundation species under climate change

As species face rapid environmental change, we can build resilient populations through restoration projects that incorporate predicted future climates into seed sourcing decisions. Eucalyptus melliodora is a foundation species of a critically endangered community in Australia that is a target for restoration. We examined genomic and phenotypic variation to make empirical based recommendations for seed sourcing. We examined isolation by distance and isolation by environment, determining high levels of gene flow extending for 500 km and correlations with climate and soil variables. Growth experiments revealed extensive phenotypic variation both within and among sampling sites, but no site-specific differentiation in phenotypic plasticity. Model predictions suggest that seed can be sourced broadly across the landscape, providing ample diversity for adaptation to environmental change. Application of our landscape genomic model to E. melliodora restoration projects can identify genomic variation suitable for predicted future climates, thereby increasing the long term probability of successful restoration.

Yellow box, or Eucalyptus melliodora, is an emblematic Australian tree that is essential to many native ecosystems. Some of these environments are now critically endangered, and replanting yellow box trees is one of the first steps to try to restore them.

However, it can be difficult for reforestation practitioners to decide where to collect the seeds they will use to replant an area. They have to select seeds that carry the genetic information that gives the trees the best chances of surviving now and in the future. This is a complex task because climate change creates fast-changing environments.

Here, Supple et al. collect genetic material from 275 E. melliodora trees at 36 different sites. Genetic analyses show that the yellow box trees from these sites exchange their genetic material and do not form isolated populations. This means that the seeds do not need to be sourced from near the reforestation site, but can be collected from areas much further away. This results in higher quality seeds for reforestation because seeds from more sites will include more genetic diversity.

Supple et al. then use information about the local climate, such as temperature and rain levels, at the sites where the samples were gathered to create a model that describes the relationship between genetic, geographical, and environmental factors. This helps identify which sites produce the seeds that will grow best under particular conditions. In addition, the seedlings from these sites are grown in the laboratory to see how well they fare in different types of environments. It therefore becomes possible to match a reforestation site with the seeds that will thrive in the future climate of the area.

Sharing this knowledge with conservationists will help to guide their replanting strategies for E. melliodora. The method can also be applied to other plant species and restoration projects, so it could ultimately shape resilient ecosystems that can cope with climate change.