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

Integrating evolutionary genomics of forest trees to inform future tree breeding amid rapid climate change

Global climate change is leading to rapid and drastic shifts in environmental conditions, posing threats to biodiversity and nearly all life forms worldwide. Forest trees serve as foundational components of terrestrial ecosystems and play a crucial and leading role in combating and mitigating the adverse effects of extreme climate events, despite their own vulnerability to these threats. Therefore, understanding and monitoring how natural forests respond to rapid climate change is a key priority for biodiversity conservation. Recent progress in evolutionary genomics, driven primarily by cutting-edge multi-omics technologies, offers powerful new tools to address several key issues. These include precise delineation of species and evolutionary units, inference of past evolutionary histories and demographic fluctuations, identification of environmentally adaptive variants, and measurement of genetic load levels. As the urgency to deal with more extreme environmental stresses grows, understanding the genomics of evolutionary history, local adaptation, future responses to climate change, and conservation and restoration of natural forest trees will be critical for research at the nexus of global change, population genomics, and conservation biology. In this review, we explore the application of evolutionary genomics to assess the effects of global climate change using multi-omics approaches and discuss the outlook for breeding of climate-adapted trees.

Genomic and Epigenomic Influences on Resilience across Scales: Lessons from the Responses of Fish to Environmental Stressors

Understanding the factors that influence the resilience of biological systems to environmental change is a pressing concern in the face of increasing human impacts on ecosystems and the organisms that inhabit them. However, most considerations of biological resilience have focused at the community and ecosystem levels, whereas here we discuss how including consideration of processes occurring at lower levels of biological organization may provide insights into factors that influence resilience at higher levels. Specifically, we explore how processes at the genomic and epigenomic levels may cascade up to influence resilience at higher levels. We ask how the concepts of "resistance," or the capacity of a system to minimize change in response to a disturbance, and "recovery," or the ability of a system to return to its original state following a disturbance and avoid tipping points and resulting regime shifts, map to these lower levels of biological organization. Overall, we suggest that substantial changes at these lower levels may be required to support resilience at higher levels, using selected examples of genomic and epigenomic responses of fish to climate-change-related stressors such as high temperature and hypoxia at the levels of the genome, epigenome, and organism.

Plant genome evolution in the genus Eucalyptus is driven by structural rearrangements that promote sequence divergence

Genomes have a highly organized architecture (nonrandom organization of functional and nonfunctional genetic elements within chromosomes) that is essential for many biological functions, particularly gene expression and reproduction. Despite the need to conserve genome architecture, a high level of structural variation has been observed within species. As species separate and diverge, genome architecture also diverges, becoming increasingly poorly conserved as divergence time increases. However, within plant genomes, the processes of genome architecture divergence are not well described. Here we use long-read sequencing and de novo assembly of 33 phylogenetically diverse, wild and naturally evolving Eucalyptus species, covering 1-50 million years of diverging genome evolution to measure genome architectural conservation and describe architectural divergence. The investigation of these genomes revealed that following lineage divergence, genome architecture is highly fragmented by rearrangements. As genomes continue to diverge, the accumulation of mutations and the subsequent divergence beyond recognition of rearrangements become the primary driver of genome divergence. The loss of syntenic regions also contribute to genome divergence but at a slower pace than that of rearrangements. We hypothesize that duplications and translocations are potentially the greatest contributors to Eucalyptus genome divergence.

Genomics for monitoring and understanding species responses to global climate change

All life forms across the globe are experiencing drastic changes in environmental conditions as a result of global climate change. These environmental changes are happening rapidly, incur substantial socioeconomic costs, pose threats to biodiversity and diminish a species' potential to adapt to future environments. Understanding and monitoring how organisms respond to human-driven climate change is therefore a major priority for the conservation of biodiversity in a rapidly changing environment. Recent developments in genomic, transcriptomic and epigenomic technologies are enabling unprecedented insights into the evolutionary processes and molecular bases of adaptation. This Review summarizes methods that apply and integrate omics tools to experimentally investigate, monitor and predict how species and communities in the wild cope with global climate change, which is by genetically adapting to new environmental conditions, through range shifts or through phenotypic plasticity. We identify advantages and limitations of each method and discuss future research avenues that would improve our understanding of species' evolutionary responses to global climate change, highlighting the need for holistic, multi-omics approaches to ecosystem monitoring during global climate change.

Exploring the role of polymorphic interspecies structural variants in reproductive isolation and adaptive divergence in Eucalyptus

Structural variations (SVs) play a significant role in speciation and adaptation in many species, yet few studies have explored the prevalence and impact of different categories of SVs. We conducted a comparative analysis of long-read assembled reference genomes of closely related Eucalyptus species to identify candidate SVs potentially influencing speciation and adaptation. Interspecies SVs can be either fixed differences or polymorphic in one or both species. To describe SV patterns, we employed short-read whole-genome sequencing on over 600 individuals of Eucalyptus melliodora and Eucalyptus sideroxylon, along with recent high-quality genome assemblies. We aligned reads and genotyped interspecies SVs predicted between species reference genomes. Our results revealed that 49,756 of 58,025 and 39,536 of 47,064 interspecies SVs could be typed with short reads in E. melliodora and E. sideroxylon, respectively. Focusing on inversions and translocations, symmetric SVs that are readily genotyped within both populations, 24 were found to be structural divergences, 2,623 structural polymorphisms, and 928 shared structural polymorphisms. We assessed the functional significance of fixed interspecies SVs by examining differences in estimated recombination rates and genetic differentiation between species, revealing a complex history of natural selection. Shared structural polymorphisms displayed enrichment of potentially adaptive genes. Understanding how different classes of genetic mutations contribute to genetic diversity and reproductive barriers is essential for understanding how organisms enhance fitness, adapt to changing environments, and diversify. Our findings reveal the prevalence of interspecies SVs and elucidate their role in genetic differentiation, adaptive evolution, and species divergence within and between populations.

Range-wide differential adaptation and genomic offset in critically endangered Asian rosewoods

In the billion-dollar global illegal wildlife trade, rosewoods have been the world's most trafficked wild product since 2005. Dalbergia cochinchinensis and Dalbergia oliveri are the most sought-after rosewoods in the Greater Mekong Subregion. They are exposed to significant genetic risks and the lack of knowledge on their adaptability limits the effectiveness of conservation efforts. Here, we present genome assemblies and range-wide genomic scans of adaptive variation, together with predictions of genomic offset to climate change. Adaptive genomic variation was differentially associated with temperature and precipitation-related variables between the species, although their natural ranges overlap. The findings are consistent with differences in pioneering ability and in drought tolerance. We predict their genomic offsets will increase over time and with increasing carbon emission pathway but at a faster pace in D. cochinchinensis than in D. oliveri. These results and the distinct gene-environment association in the eastern coastal edge of Vietnam suggest species-specific conservation actions: germplasm representation across the range in D. cochinchinensis and focused on hotspots of genomic offset in D. oliveri. We translated our genomic models into a seed source matching application, seedeR, to rapidly inform restoration efforts. Our ecological genomic research uncovering contrasting selection forces acting in sympatric rosewoods is of relevance to conserving tropical trees globally and combating risks from climate change.

Genetic assessment of the value of restoration planting within an endangered eucalypt woodland

Assessment of woodland restoration often focusses on stand demographics, but genetic factors likely influence long-term stand viability. We examined the genetic composition of Yellow Box (Eucalyptus melliodora) trees in endangered Box-Gum Grassy Woodland in SE Australia, some 30 years after planting with seeds of reportedly local provenance. Using DArT sequencing for 1406 SNPs, we compared genetic diversity and population structure of planted E. melliodora trees with remnant bushland trees, paddock trees and natural recruits. Genetic patterns imply that natural stands and paddock trees had historically high gene flow (among group pairwise FST = 0.04-0.10). Genetic diversity was highest among relictual paddock trees (He = 0.17), while diversity of revegetated trees was identical to natural bushland trees (He = 0.14). Bayesian clustering placed the revegetated trees into six genetic groups with four corresponding to genotypes from paddock trees, indicating that revegetated stands are mainly of genetically diverse, local provenance. Natural recruits were largely derived from paddock trees with some contribution from planted trees. A few trees have likely hybridised with other local eucalypt species which are unlikely to compromise stand integrity. We show that paddock trees have high genetic diversity and capture historic genetic variety and provide important foci for natural recruitment of genetically diverse and outcrossed seedlings.

Forecasting trait responses in novel environments to aid seed provenancing under climate change

Revegetation projects face the major challenge of sourcing optimal plant material. This is often done with limited information about plant performance and increasingly requires factoring resilience to climate change. Functional traits can be used as quantitative indices of plant performance and guide seed provenancing, but trait values expected under novel conditions are often unknown. To support climate-resilient provenancing efforts, we develop a trait prediction model that integrates the effect of genetic variation with fine-scale temperature variation. We train our model on multiple field plantings of Arabidopsis thaliana and predict two relevant fitness traits-days-to-bolting and fecundity-across the species' European range. Prediction accuracy was high for days-to-bolting and moderate for fecundity, with the majority of trait variation explained by temperature differences between plantings. Projection under future climate predicted a decline in fecundity, although this response was heterogeneous across the range. In response, we identified novel genotypes that could be introduced to genetically offset the fitness decay. Our study highlights the value of predictive models to aid seed provenancing and improve the success of revegetation projects.

The genomic and epigenetic footprint of local adaptation to variable climates in kiwifruit

A full understanding of adaptive genetic variation at the genomic level will help address questions of how organisms adapt to diverse climates. Actinidia eriantha is a shade-tolerant species, widely distributed in the southern tropical region of China, occurring in spatially heterogeneous environments. In the present study we combined population genomic, epigenomic, and environmental association analyses to infer population genetic structure and positive selection across a climatic gradient, and to assess genomic offset to climatic change for A. eriantha. The population structure is strongly shaped by geography and influenced by restricted gene flow resulting from isolation by distance due to habitat fragmentation. In total, we identified 102 outlier loci and annotated 455 candidate genes associated with the genomic basis of climate adaptation, which were enriched in functional categories related to development processes and stress response; both temperature and precipitation are important factors driving adaptive variation. In addition to single-nucleotide polymorphisms (SNPs), a total of 27 single-methylation variants (SMVs) had significant correlation with at least one of four climatic variables and 16 SMVs were located in or adjacent to genes, several of which were predicted to be involved in plant response to abiotic or biotic stress. Gradient forest analysis indicated that the central/east populations were predicted to be at higher risk of future population maladaptation under climate change. Our results demonstrate that local climate factors impose strong selection pressures and lead to local adaptation. Such information adds to our understanding of adaptive mechanisms to variable climates revealed by both population genome and epigenome analysis.

Understanding local plant extinctions before it is too late: bridging evolutionary genomics with global ecology

Understanding evolutionary genomic and population processes within a species range is key to anticipating the extinction of plant species before it is too late. However, most models of biodiversity risk under global change do not account for the genetic variation and local adaptation of different populations. Population diversity is critical to understanding extinction because different populations may be more or less susceptible to global change and, if lost, would reduce the total diversity within a species. Two new modeling frameworks advance our understanding of extinction from a population and evolutionary angle: Rapid climate change-driven disruptions in population adaptation are predicted from associations between genomes and local climates. Furthermore, losses of population diversity from global land-use transformations are estimated by scaling relationships of species' genomic diversity with habitat area. Overall, these global eco-evolutionary methods advance the predictability - and possibly the preventability - of the ongoing extinction of plant species.

Interspecies genome divergence is predominantly due to frequent small scale rearrangements in Eucalyptus

Synteny, the ordering of sequences within homologous chromosomes, must be maintained within the genomes of sexually reproducing species for the sharing of alleles and production of viable, reproducing offspring. However, when the genomes of closely related species are compared, a loss of synteny is often observed. Unequal homologous recombination is the primary mechanism behind synteny loss, occurring more often in transposon rich regions, and resulting in the formation of chromosomal rearrangements. To examine patterns of synteny among three closely related, interbreeding, and wild Eucalyptus species, we assembled their genomes using long-read DNA sequencing and de novo assembly. We identify syntenic and rearranged regions between these genomes and estimate that ~48% of our genomes remain syntenic while ~36% is rearranged. We observed that rearrangements highly fragment microsynteny. Our results suggest that synteny between these species is primarily lost through small-scale rearrangements, not through sequence loss, gain, or sequence divergence. Further examination of identified rearrangements suggests that rearrangements may be altering the phenotypes of Eucalyptus species. Our study also underscores that the use of single reference genomes in genomic variation studies could lead to reference bias, especially given the scale at which we show potentially adaptive loci have highly diverged, deleted, duplicated and/or rearranged. This study provides an unbiased framework to look at potential speciation and adaptive loci among a rapidly radiating foundation species of woodland trees that are free from selective breeding seen in most crop species.

Multiscale spatial analysis of headwater vulnerability in South-Central Chile reveals a high threat due to deforestation and climate change

Headwaters represent an essential component of hydrological, ecological, and socioeconomical systems, by providing constant water streams to the complete basin. However, despite the high importance of headwaters, there is a lack of vulnerability assessments worldwide. Identifying headwaters and their vulnerability in a spatially explicit manner can enable restauration and conservation programs. In this study, we assess the vulnerability of headwaters in South-Central Chile (38.4 to 43.2°S) considering multiple degradation factors related to climate change and land cover change. We analyzed 2292 headwaters, characterizing multiple factors at five spatial scales by using remote sensing data related to Land Use and Cover Change (LUCC), human disturbances, vegetation cover, climate change, potential water demand, and physiography. We then generated an index of vulnerability by integrating all the analyzed variables, which allowed us to map the spatial distribution of headwater vulnerability. Finally, to estimate the main drivers of degradation, we performed a Principal Components Analysis with an Agglomerative Hierarchical Clustering, that allowed us to group headwaters according to the analyzed factors. The largest proportion of most vulnerable headwaters are located in the north of our study area with 48.1 %, 62.1 %, and 28.1 % of headwaters classified as highly vulnerable at 0, 10, and 30 m scale, respectively. The largest proportion of headwaters are affected by Climate Change (63.66 %) and LUCC (23.02 %) on average across all scales. However, we identified three clusters, in which the northern cluster is mainly affected by LUCC, while the Andean and Coastal clusters are mainly affected by climate change. Our results and methods present an informative picture of the current state of headwater vulnerability, identifying spatial patterns and drivers at multiple scales. We believe that the approach developed in this study could be useful for new studies in other zones of the world and can also promote Chilean headwater conservation.

A box on the river: The phylogenetics and phylogeography of Eucalyptus baueriana (Eucalyptus sect. Adnataria ser. Heterophloiae)

We present a phylogeographic study of the tree species Eucalyptus baueriana Schauer, which occurs in disjunct areas on the near coastal plains and ranges of the south-east Australian mainland. DArTseq data are used to build a phylogeny including E. baueriana and closely related taxa to test its monophyly, test the genetic distinctness of the three subspecies of E. baueriana, and investigate relationships between its disjunct populations. Additionally, we use population structure analysis to investigate the genetic distinctness of populations, and MaxEnt to investigate the environmental factors potentially influencing the species' distribution. We show E. baueriana is monophyletic and most closely related to three other Blue Box eucalypt species: E. conica H.Deane & Maiden, E. dalveenica T.L.Collins, R.L.Andrew & J.J.Bruhl and E. magnificata L.A.S.Johnson & K.D.Hill, with some evidence for genetic introgression between these taxa. Within E. baueriana, the deepest genetic breaks do not correspond with the subspecies classification as the two geographically restricted subspecies, together with samples of the more widespread E. baueriana subsp. baueriana from west of the Gippsland lowlands, form a south-western clade with that is sister to other populations of subsp. baueriana. The oldest genetic break in the species occurs in far eastern Gippsland (Victoria), corresponding to one of the shortest geographic disjunctions in the species' distribution. Genetic breaks in other species have been observed in this region which is broadly referred to as the southern transition zone. Both total annual rainfall and the seasonality of this rainfall are hypothesised to affect the species' distribution; gaps in its distribution are in areas of higher rainfall that support closed forest and in regions with more winter dominated rainfall.

Genetic diversity loss in the Anthropocene

Anthropogenic habitat loss and climate change are reducing species' geographic ranges, increasing extinction risk and losses of species' genetic diversity. Although preserving genetic diversity is key to maintaining species' adaptability, we lack predictive tools and global estimates of genetic diversity loss across ecosystems. We introduce a mathematical framework that bridges biodiversity theory and population genetics to understand the loss of naturally occurring DNA mutations with decreasing habitat. By analyzing genomic variation of 10,095 georeferenced individuals from 20 plant and animal species, we show that genome-wide diversity follows a mutations-area relationship power law with geographic area, which can predict genetic diversity loss from local population extinctions. We estimate that more than 10% of genetic diversity may already be lost for many threatened and nonthreatened species, surpassing the United Nations' post-2020 targets for genetic preservation.

Climate Adaptation, Drought Susceptibility, and Genomic-Informed Predictions of Future Climate Refugia for the Australian Forest Tree Eucalyptus globulus

Understanding the capacity of forest tree species to adapt to climate change is of increasing importance for managing forest genetic resources. Through a genomics approach, we modelled spatial variation in climate adaptation within the Australian temperate forest tree Eucalyptus globulus, identified putative climate drivers of this genomic variation, and predicted locations of future climate refugia and populations at-risk of future maladaptation. Using 812,158 SNPs across 130 individuals from 30 populations (i.e., localities) spanning the species’ natural range, a gradientForest algorithm found 1177 SNPs associated with locality variation in home-site climate (climate-SNPs), putatively linking them to climate adaptation. Very few climate-SNPs were associated with population-level variation in drought susceptibility, signalling the multi-faceted nature and complexity of climate adaptation. Redundancy analysis (RDA) showed 24% of the climate-SNP variation could be explained by annual precipitation, isothermality, and maximum temperature of the warmest month. Spatial predictions of the RDA climate vectors associated with climate-SNPs allowed mapping of genomically informed climate selective surfaces across the species’ range under contemporary and projected future climates. These surfaces suggest over 50% of the current distribution of E. globulus will be outside the modelled adaptive range by 2070 and at risk of climate maladaptation. Such surfaces present a new integrated approach for natural resource managers to capture adaptive genetic variation and plan translocations in the face of 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.

Using genomics to guide seed-sourcing at the right taxonomical level for ecological restoration projects: The complex case of Carex bigelowii s.lat. in Norway

There is a growing demand for ecological restoration using suitable seeds following international standards or national legal demands for local seed-sourcing. However, before selecting the appropriate geographic origin of seeds, it is vital to explore taxonomic complexity related to the focal taxa. We used ddRAD-seq to screen genomic diversity within Carex bigelowii s.lat. focussing on Norway. This species complex is considered a candidate for seeding, but presents considerable morphological, ecological, and genetic variation. The genetic structure of 132 individuals of C. bigelowii s.lat., including Carex nigra as an outgroup, was explored using ordinations, clustering analyses, and a genetic barrier algorithm. Two highly divergent clusters were evident, supporting the recognition of two taxonomic units "C. dacica" and C. bigelowii "subsp. bigelowii". Previously defined seed-sourcing regions for C. bigelowii s.lat. did not consider the known taxonomic complexity, and therefore interpreted the overall genetic structure as seed-sourcing regions, not taxa. We estimated genetic neighborhood sizes within each taxon to be 100-150 km and 300 km, respectively, indicating species-specific delimitations of local seed-sourcing regions. Frequent hybrids, local genetic distinctiveness, and suggested ecotypes add complexity to the discussed seed-sourcing regions. Our results show how genomic screening of diversity and structure in a species complex can alleviate the taxonomic impediment, inform practical questions, and legal requirements related to seed-sourcing, and together with traditional taxonomic work provide necessary information for a sound management of biodiversity.

Chromosome-specific polymorphic SSR markers in tropical eucalypt species using low coverage whole genome sequences: systematic characterization and validation

Eucalyptus is one of the major plantation species with wide variety of industrial uses. Polymorphic and informative simple sequence repeats (SSRs) have broad range of applications in genetic analysis. In this study, two individuals of Eucalyptus tereticornis (ET217 and ET86), one individual each from E. camaldulensis (EC17) and E. grandis (EG9) were subjected to whole genome resequencing. Low coverage (10×) genome sequencing was used to find polymorphic SSRs between the individuals. Average number of SSR loci identified was 95,513 and the density of SSRs per Mb was from 157.39 in EG9 to 155.08 in EC17. Among all the SSRs detected, the most abundant repeat motifs were di-nucleotide (59.6%–62.5%), followed by tri- (23.7%–27.2%), tetra- (5.2%–5.6%), penta- (5.0%–5.3%), and hexa-nucleotide (2.7%–2.9%). The predominant SSR motif units were AG/CT and AAG/TTC. Computational genome analysis predicted the SSR length variations between the individuals and identified the gene functions of SSR containing sequences. Selected subset of polymorphic markers was validated in a full-sib family of eucalypts. Additionally, genome-wide characterization of single nucleotide polymorphisms, InDels and transcriptional regulators were carried out. These variations will find their utility in genome-wide association studies as well as understanding of molecular mechanisms involved in key economic traits. The genomic resources generated in this study would provide an impetus to integrate genomics in marker-trait associations and breeding of tropical eucalypts.

Achievements and Challenges of Genomics-Assisted Breeding in Forest Trees: From Marker-Assisted Selection to Genome Editing

Forest tree breeding efforts have focused mainly on improving traits of economic importance, selecting trees suited to new environments or generating trees that are more resilient to biotic and abiotic stressors. This review describes various methods of forest tree selection assisted by genomics and the main technological challenges and achievements in research at the genomic level. Due to the long rotation time of a forest plantation and the resulting long generation times necessary to complete a breeding cycle, the use of advanced techniques with traditional breeding have been necessary, allowing the use of more precise methods for determining the genetic architecture of traits of interest, such as genome-wide association studies (GWASs) and genomic selection (GS). In this sense, main factors that determine the accuracy of genomic prediction models are also addressed. In turn, the introduction of genome editing opens the door to new possibilities in forest trees and especially clustered regularly interspaced short palindromic repeats and CRISPR-associated protein 9 (CRISPR/Cas9). It is a highly efficient and effective genome editing technique that has been used to effectively implement targetable changes at specific places in the genome of a forest tree. In this sense, forest trees still lack a transformation method and an inefficient number of genotypes for CRISPR/Cas9. This challenge could be addressed with the use of the newly developing technique GRF-GIF with speed breeding.

OMICs, Epigenetics, and Genome Editing Techniques for Food and Nutritional Security

The incredible success of crop breeding and agricultural innovation in the last century greatly contributed to the Green Revolution, which significantly increased yields and ensures food security, despite the population explosion. However, new challenges such as rapid climate change, deteriorating soil, and the accumulation of pollutants require much faster responses and more effective solutions that cannot be achieved through traditional breeding. Further prospects for increasing the efficiency of agriculture are undoubtedly associated with the inclusion in the breeding strategy of new knowledge obtained using high-throughput technologies and new tools in the future to ensure the design of new plant genomes and predict the desired phenotype. This article provides an overview of the current state of research in these areas, as well as the study of soil and plant microbiomes, and the prospective use of their potential in a new field of microbiome engineering. In terms of genomic and phenomic predictions, we also propose an integrated approach that combines high-density genotyping and high-throughput phenotyping techniques, which can improve the prediction accuracy of quantitative traits in crop species.

Genomic Approaches for Conservation Management in Australia under Climate Change

Conservation genetics has informed threatened species management for several decades. With the advent of advanced DNA sequencing technologies in recent years, it is now possible to monitor and manage threatened populations with even greater precision. Climate change presents a number of threats and challenges, but new genomics data and analytical approaches provide opportunities to identify critical evolutionary processes of relevance to genetic management under climate change. Here, we discuss the applications of such approaches for threatened species management in Australia in the context of climate change, identifying methods of facilitating viability and resilience in the face of extreme environmental stress. Using genomic approaches, conservation management practices such as translocation, targeted gene flow, and gene-editing can now be performed with the express intention of facilitating adaptation to current and projected climate change scenarios in vulnerable species, thus reducing extinction risk and ensuring the protection of our unique biodiversity for future generations. We discuss the current barriers to implementing conservation genomic projects and the efforts being made to overcome them, including communication between researchers and managers to improve the relevance and applicability of genomic studies. We present novel approaches for facilitating adaptive capacity and accelerating natural selection in species to encourage resilience in the face of climate change.

A decision framework for evidence-based climate adaptation interventions

Climate change is outpacing the ability of species and populations to naturally adapt warranting active interventions to boost climate resilience. In their review paper, Gaitán-Espitia and Hobday discuss how such interventions may, however, counter natural evolutionary processes and adaptive capacity if not underpinned by background knowledge from genes through to ecosystems. They present a robust decision framework to guide implementation of climate adaptation interventions to avoid unintended evolutionary outcomes.

Divergent lineages in a semi-arid mallee species, Eucalyptus behriana, correspond to a major geographic break in southeastern Australia

AIM

To infer relationships between populations of the semi-arid, mallee eucalypt, Eucalyptus behriana, to build hypotheses regarding evolution of major disjunctions in the species' distribution and to expand understanding of the biogeographical history of southeastern Australia.

LOCATION

Southeastern Australia.

TAXON

Eucalyptus behriana (Myrtaceae, Angiospermae).

METHODS

We developed a large dataset of anonymous genomic loci for 97 samples from 11 populations of E. behriana using double digest restriction site-associated DNA sequencing (ddRAD-seq), to determine genetic relationships between the populations. These relationships, along with species distribution models, were used to construct hypotheses regarding environmental processes that have driven fragmentation of the species' distribution.

RESULTS

Greatest genetic divergence was between populations on either side of the Lower Murray Basin. Populations west of the Basin showed greater genetic divergence between one another than the eastern populations. The most genetically distinct population in the east (Long Forest) was separated from others by the Great Dividing Range. A close relationship was found between the outlying northernmost population (near West Wyalong) and those in the Victorian Goldfields despite a large disjunction between them.

CONCLUSIONS

Patterns of genetic variation are consistent with a history of vicariant differentiation of disjunct populations. We infer that an early disjunction to develop in the species distribution was that across the Lower Murray Basin, an important biogeographical barrier separating many dry sclerophyll plant taxa in southeastern Australia. Additionally, our results suggest that the western populations fragmented earlier than the eastern ones. Fragmentation, both west and east of the Murray Basin, is likely tied to climatic changes associated with glacial-interglacial cycles although it remains possible that major geological events including uplift of the Mount Lofty Ranges and basalt flows in the Newer Volcanics Province also played a role.

Genomic Prediction of (Mal)Adaptation Across Current and Future Climatic Landscapes

Signals of local adaptation have been found in many plants and animals, highlighting the heterogeneity in the distribution of adaptive genetic variation throughout species ranges. In the coming decades, global climate change is expected to induce shifts in the selective pressures that shape this adaptive variation. These changes in selective pressures will likely result in varying degrees of local climate maladaptation and spatial reshuffling of the underlying distributions of adaptive alleles. There is a growing interest in using population genomic data to help predict future disruptions to locally adaptive gene-environment associations. One motivation behind such work is to better understand how the effects of changing climate on populations’ short-term fitness could vary spatially across species ranges. Here we review the current use of genomic data to predict the disruption of local adaptation across current and future climates. After assessing goals and motivationsunderlying the approach, we review the main steps and associated statistical methods currently in use and explore our current understanding of the limits and future potential of using genomics to predict climate change (mal)adaptation.

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.

Genomic evidence of introgression and adaptation in a model subtropical tree species, Eucalyptus grandis

The genetic consequences of adaptation to changing environments can be deciphered using population genomics, which may help predict species' responses to global climate change. Towards this, we used genome-wide SNP marker analysis to determine population structure and patterns of genetic differentiation in terms of neutral and adaptive genetic variation in the natural range of Eucalyptus grandis, a widely cultivated subtropical and temperate species, serving as genomic reference for the genus. We analysed introgression patterns at subchromosomal resolution using a modified ancestry mapping approach and identified provenances with extensive interspecific introgression in response to increased aridity. Furthermore, we describe potentially adaptive genetic variation as explained by environment-associated SNP markers, which also led to the discovery of what is likely a large structural variant. Finally, we show that genes linked to these markers are enriched for biotic and abiotic stress responses.

Perceptions of Similarity Can Mislead Provenancing Strategies—An Example from Five Co-Distributed Acacia Species

Ecological restoration requires balancing levels of genetic diversity to achieve present-day establishment as well as long-term sustainability. Assumptions based on distributional, taxonomic or functional generalizations are often made when deciding how to source plant material for restoration. We investigate this assumption and ask whether species-specific data is required to optimize provenancing strategies. We use population genetic and environmental data from five congeneric and largely co-distributed species of Acacia to specifically ask how different species-specific genetic provenancing strategies are based on empirical data and how well a simple, standardized collection strategy would work when applied to the same species. We find substantial variability in terms of patterns of genetic diversity and differentiation across the landscape among these five co-distributed Acacia species. This variation translates into substantial differences in genetic provenancing recommendations among species (ranging from 100% to less than 1% of observed genetic variation across species) that could not have been accurately predicted a priori based on simple observation or overall distributional patterns. Furthermore, when a common provenancing strategy was applied to each species, the recommended collection areas and the evolutionary representativeness of such artificially standardized areas were substantially different (smaller) from those identified based on environmental and genetic data. We recommend the implementation of the increasingly accessible array of evolutionary-based methodologies and information to optimize restoration efforts.

Contrasting patterns of local adaptation along climatic gradients between a sympatric parasitic and autotrophic tree species

Sympatric tree species are subject to similar climatic drivers, posing a question as to whether they display comparable adaptive responses. However, no study has explicitly examined local adaptation of co-occurring parasitic and autotrophic plant species to the abiotic environment. Here we test the hypotheses that a generalist parasitic tree would display a weaker signal of selection and that genomic variation would associate with fewer climatic variables (particularly precipitation) but have similar spatial patterns to a sympatric autotrophic tree species. To test these hypotheses, we collected samples from 17 sites across the range of two tree species, the hemiparasite Nuytsia floribunda (n = 264) and sympatric autotroph Melaleuca rhaphiophylla (n = 272). We obtained 5,531 high-quality genome-wide single nucleotide polymorphisms (SNPs) for M. rhaphiophylla and 6,727 SNPs for N. floribunda using DArTseq genome scan technology. Population differentiation and environmental association approaches were used to identify signals of selection. Generalized dissimilarly modelling was used to detect climatic and spatial patterns of local adaptation across climatic gradients. Overall, 322 SNPs were identified as putatively adaptive for the autotroph, while only 57 SNPs were identified for the parasitic species. We found genomic variation to associate with different sets of bioclimatic variables for each species, with precipitation relatively less important for the parasite. Spatial patterns of predicted adaptive variability were different and indicate that co-occurring species with disparate life history traits may not respond equally to selective pressures (i.e., temperature and precipitation). Together, these findings provide insight into local adaptation of sympatric parasitic and autotrophic tree species to abiotic environments.

Combining genotype, phenotype, and environmental data to delineate site-adjusted provenance strategies for ecological restoration

Despite the importance of climate-adjusted provenancing to mitigate the effects of environmental change, climatic considerations alone are insufficient when restoring highly degraded sites. Here we propose a comprehensive landscape genomic approach to assist the restoration of moderately disturbed and highly degraded sites. To illustrate it we employ genomic data sets comprising thousands of single nucleotide polymorphisms from two plant species suitable for the restoration of iron-rich Amazonian Savannas. We first use a subset of neutral loci to assess genetic structure and determine the genetic neighbourhood size. We then identify genotype-phenotype-environment associations, map adaptive genetic variation, and predict adaptive genotypes for restoration sites. Whereas local provenances were found optimal to restore a moderately disturbed site, a mixture of genotypes seemed the most promising strategy to recover a highly degraded mining site. We discuss how our results can help define site-adjusted provenancing strategies, and argue that our methods can be more broadly applied to assist other restoration initiatives.

Range-wide neutral and adaptive genetic structure of an endemic herb from Amazonian Savannas

Conserving genetic diversity in rare and narrowly distributed endemic species is essential to maintain their evolutionary potential and minimize extinction risk under future environmental change. In this study we assess neutral and adaptive genetic structure and genetic diversity in Brasilianthus carajensis (Melastomataceae), an endemic herb from Amazonian Savannas. Using RAD sequencing we identified a total of 9365 SNPs in 150 individuals collected across the species' entire distribution range. Relying on assumption-free genetic clustering methods and environmental association tests we then compared neutral with adaptive genetic structure. We found three neutral and six adaptive genetic clusters, which could be considered management units (MU) and adaptive units (AU), respectively. Pairwise genetic differentiation (F _ST) ranged between 0.024 and 0.048, and even though effective population sizes were below 100, no significant inbreeding was found in any inferred cluster. Nearly 10 % of all analysed sequences contained loci associated with temperature and precipitation, from which only 25 sequences contained annotated proteins, with some of them being very relevant for physiological processes in plants. Our findings provide a detailed insight into genetic diversity, neutral and adaptive genetic structure in a rare endemic herb, which can help guide conservation and management actions to avoid the loss of unique genetic variation.

Genome-wide signatures of environmental adaptation in European aspen (Populus tremula) under current and future climate conditions

Future climate change has been predicted to disrupt local adaptation in many perennial plants, such as forest trees, but the magnitude and location of these effects are thus far poorly understood. Here, we assess local adaptation to current climate in European aspen (Populus tremula) by using environmental association analyses to identify genetic variants associated with two representative climate variables describing current day variation in temperature and precipitation. We also analysed patterns of genetic differentiation between southern and northern populations and observe that regions of high genetic differentiation are enriched for SNPs that are significantly associated with climate. Using variants associated with climate, we examined patterns of isolation by distance and environment and used spatial modelling to predict the geographic distribution of genomic variation in response to two scenarios of future climate change. We show that climate conditions at a northern reference site will correspond to climate conditions experienced by current day populations located 4-8 latitude degrees further south. By assessing the relationship between phenotypic traits and vegetative fitness, we also demonstrate that southern populations harbour genetic variation that likely would be adaptive further north under both climate change scenarios. Current day populations at the lagging edge of the distribution in Sweden can therefore serve as sources for introducing adaptive alleles onto northern populations, but the likelihood of this largely depends on naturally occurring levels of gene flow.

Reframing conservation physiology to be more inclusive, integrative, relevant and forward-looking: reflections and a horizon scan

Applying physiological tools, knowledge and concepts to understand conservation problems (i.e. conservation physiology) has become commonplace and confers an ability to understand mechanistic processes, develop predictive models and identify cause-and-effect relationships. Conservation physiology is making contributions to conservation solutions; the number of 'success stories' is growing, but there remain unexplored opportunities for which conservation physiology shows immense promise and has the potential to contribute to major advances in protecting and restoring biodiversity. Here, we consider how conservation physiology has evolved with a focus on reframing the discipline to be more inclusive and integrative. Using a 'horizon scan', we further explore ways in which conservation physiology can be more relevant to pressing conservation issues of today (e.g. addressing the Sustainable Development Goals; delivering science to support the UN Decade on Ecosystem Restoration), as well as more forward-looking to inform emerging issues and policies for tomorrow. Our horizon scan provides evidence that, as the discipline of conservation physiology continues to mature, it provides a wealth of opportunities to promote integration, inclusivity and forward-thinking goals that contribute to achieving conservation gains. To advance environmental management and ecosystem restoration, we need to ensure that the underlying science (such as that generated by conservation physiology) is relevant with accompanying messaging that is straightforward and accessible to end users.

Landscape drivers of genomic diversity and divergence in woodland Eucalyptus

Spatial genetic patterns are influenced by numerous factors, and they can vary even among coexisting, closely related species due to differences in dispersal and selection. Eucalyptus (L'Héritier 1789; the "eucalypts") are foundation tree species that provide essential habitat and modulate ecosystem services throughout Australia. Here we present a study of landscape genomic variation in two woodland eucalypt species, using whole-genome sequencing of 388 individuals of Eucalyptus albens and Eucalyptus sideroxylon. We found exceptionally high genetic diversity (π ≈ 0.05) and low genome-wide, interspecific differentiation (FST  = 0.15) and intraspecific differentiation between localities (FST  ≈ 0.01-0.02). We found no support for strong, discrete population structure, but found substantial support for isolation by geographic distance (IBD) in both species. Using generalized dissimilarity modelling, we identified additional isolation by environment (IBE). Eucalyptus albens showed moderate IBD, and environmental variables have a small but significant amount of additional predictive power (i.e. IBE). Eucalyptus sideroxylon showed much stronger IBD and moderate IBE. These results highlight the vast adaptive potential of these species and set the stage for testing evolutionary hypotheses of interspecific adaptive differentiation across environments.

The potential of genomics for restoring ecosystems and biodiversity

Billions of hectares of natural ecosystems have been degraded through human actions. The global community has agreed on targets to halt and reverse these declines, and the restoration sector faces the important but arduous task of implementing programmes to meet these objectives. Existing and emerging genomics tools offer the potential to improve the odds of achieving these targets. These tools include population genomics that can improve seed sourcing, meta-omics that can improve assessment and monitoring of restoration outcomes, and genome editing that can generate novel genotypes for restoring challenging environments. We identify barriers to adopting these tools in a restoration context and emphasize that regulatory and ethical frameworks are required to guide their use.

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.