Cytogenetic Bioindication in Root Meristems for Vitality Assessment of Trees.

Our method describes how to collect forest tree root tips in the field, to store them for transfer to the lab, to pretreat root tips in order to arrest cells in metaphase, fix root tips to preserve specific morphological organizations, to stain fixed root tips by Feulgen's Reaction in order to increase contrast, and to prepare the root meristem for analyzing mitotic stages and chromosomal aberrations via light microscopy. We further describe how to classify chromosomal abnormalities and quantify them via aberration indices.

The drivers of intraspecific trait variation and their implications for future tree productivity and survival.

Forests are facing unprecedented levels of stress from pest and disease outbreaks, disturbance, fragmentation, development, and a changing climate. These selective agents act to alter forest composition from regional to cellular levels. Thus, a central challenge for understanding how forests will be impacted by future change is how to integrate across scales of biology. Phenotype, or an observable trait, is the product of an individual's genes (G) and the environment in which an organism lives (E). To date, researchers have detailed how environment drives variation in tree phenotypes over long time periods (e.g., long-term ecological research sites [LTERs]) and across large spatial scales (e.g., flux network). In parallel, researchers have discovered the genes and pathways that govern phenotypes, finding high degrees of genetic control and signatures of local adaptation in many plant traits. However, the research in these two areas remain largely independent of each other, hindering our ability to generate accurate predictions of plant response to environment, an increasingly urgent need given threats to forest systems. I present the importance of both genes and environment in determining tree responses to climate stress. I highlight why the difference between G versus E in driving variation is critical for our understanding of climate responses, then propose means of accelerating research that examines G and E simultaneously by leveraging existing long-term, large-scale phenotypic data sets from ecological networks and adding newly affordable sequence (-omics) data to both drill down to find the genes and alleles influencing phenotypes and scale up to find how patterns of demography and local adaptation may influence future response to change.

Phytophthora pluvialis maintenance, spore production and detached needle assays.

Phytophthora pluvialis is an oomycete that primarily infects Pinus radiata and Pseudotsuga menziesii causing the destructive foliar disease red needle cast (RNC). Recent observations show that P. pluvialis can also infect western hemlock inducing resinous cankers. High-throughput and reproducible infection assays are integral to find key information on tree health and oomycete pathogenicity. In this protocol, we describe the propagation and spore induction of P. pluvialis, followed by detached needle assays for verification and quantification of virulence of P. pluvialis in P. radiata needles. These needle assays can be employed for high-throughput screening of tree needles with diverse genetic backgrounds. In downstream analysis, Quantitative PCR (qPCR) was utilized to assess relative gene expression, as exemplified by candidate RxLR effector protein PpR01. Additional techniques like RNA sequencing, metabolomics, and proteomics can be combined with needle assays and can offer comprehensive insights into P. pluvialis infection mechanisms.

The super-pangenome of Populus unveils genomic facets for its adaptation and diversification in widespread forest trees.

Understanding the underlying mechanisms and links between genome evolution and adaptive innovations stands as a key goal in evolutionary studies. Poplars, among the world's most widely distributed and cultivated trees, exhibit extensive phenotypic diversity and environmental adaptability. In this study, we present a genus-level super-pangenome comprising 19 Populus genomes, revealing the likely pivotal role of private genes in facilitating local environmental and climate adaptation. Through the integration of pangenomes with transcriptomes, methylomes, and chromatin accessibility mapping, we unveil that the evolutionary trajectories of pangenes and duplicated genes are closely linked to local genomic landscapes of regulatory and epigenetic architectures, notably CG methylation in gene-body regions. Further comparative genomic analyses have enabled the identification of 142 202 structural variants across species that intersect with a significant number of genes and contribute substantially to both phenotypic and adaptive divergence. We have experimentally validated a ∼180-bp presence/absence variant affecting the expression of the CUC2 gene, crucial for leaf serration formation. Finally, we developed a user-friendly web-based tool encompassing the multi-omics resources associated with the Populus super-pangenome (http://www.populus-superpangenome.com). Together, the present pioneering super-pangenome resource in forest trees not only aids in the advancement of breeding efforts of this globally important tree genus but also offers valuable insights into potential avenues for comprehending tree biology.

Phytochrome-interacting factors play shared and distinct roles in regulating shade avoidance responses in Populus trees.

Plants adjust their growth and development in response to changing light caused by canopy shade. The molecular mechanisms underlying shade avoidance responses have been widely studied in Arabidopsis and annual crop species, yet the shade avoidance signalling in woody perennial trees remains poorly understood. Here, we first showed that PtophyB1/2 photoreceptors serve conserved roles in attenuating the shade avoidance syndrome (SAS) in poplars. Next, we conducted a systematic identification and characterization of eight PtoPIF genes in Populus tomentosa. Knocking out different PtoPIFs led to attenuated shade responses to varying extents, whereas overexpression of PtoPIFs, particularly PtoPIF3.1 and PtoPIF3.2, led to constitutive SAS phenotypes under normal light and enhanced SAS responses under simulated shade. Notably, our results revealed that distinct from Arabidopsis PIF4 and PIF5, which are major regulators of SAS, the Populus homologues PtoPIF4.1 and PtoPIF4.2 seem to play a minor role in controlling shade responses. Moreover, we showed that PtoPIF3.1/3.2 could directly activate the expression of the auxin biosynthetic gene PtoYUC8 in response to shade, suggesting a conserved PIF-YUC-auxin pathway in modulating SAS in tree. Overall, our study provides insights into shared and divergent functions of PtoPIF members in regulating various aspects of the SAS in Populus.

The genome assembly of Island Oak (Quercus tomentella), a relictual island tree species.

Island oak (Quercus tomentella) is a rare relictual island tree species that exists only on six islands off the coast of California and Mexico, but was once widespread throughout mainland California. Currently, this species is endangered by threats such as non-native plants, grazing animals, and human removal. Efforts for conservation and restoration of island oak currently underway could benefit from information about its range-wide genetic structure and evolutionary history. Here we present a high-quality genome assembly for Q. tomentella, assembled using PacBio HiFi and Omni-C sequencing, developed as part of the California Conservation Genomics Project (CCGP). The resulting assembly has a length of 781 Mb, with a contig N50 of 22.0 Mb and a scaffold N50 of 63.4 Mb. This genome assembly will provide a resource for genomics-informed conservation of this rare oak species. Additionally, this reference genome will be the first one available for a species in Quercus section Protobalanus, a unique oak clade present only in western North America.

Low-frequency somatic mutations are heritable in tropical trees Dicorynia guianensis and Sextonia rubra.

Somatic mutations potentially play a role in plant evolution, but common expectations pertaining to plant somatic mutations remain insufficiently tested. Unlike in most animals, the plant germline is assumed to be set aside late in development, leading to the expectation that plants accumulate somatic mutations along growth. Therefore, several predictions were made on the fate of somatic mutations: mutations have generally low frequency in plant tissues; mutations at high frequency have a higher chance of intergenerational transmission; branching topology of the tree dictates mutation distribution; and exposure to UV (ultraviolet) radiation increases mutagenesis. To provide insights into mutation accumulation and transmission in plants, we produced two high-quality reference genomes and a unique dataset of 60 high-coverage whole-genome sequences of two tropical tree species, Dicorynia guianensis (Fabaceae) and Sextonia rubra (Lauraceae). We identified 15,066 de novo somatic mutations in D. guianensis and 3,208 in S. rubra, surprisingly almost all found at low frequency. We demonstrate that 1) low-frequency mutations can be transmitted to the next generation; 2) mutation phylogenies deviate from the branching topology of the tree; and 3) mutation rates and mutation spectra are not demonstrably affected by differences in UV exposure. Altogether, our results suggest far more complex links between plant growth, aging, UV exposure, and mutation rates than commonly thought.

Transcriptome profile of pecan scab resistant and susceptible trees from a pecan provenance collection.

Pecan scab is a devastating disease that causes damage to pecan (Carya illinoinensis (Wangenh.) K. Koch) fruit and leaves. The disease is caused by the fungus Venturia effusa (G. Winter) and the main management practice for controlling the disease is by application of fungicides at 2-to-3-week intervals throughout the growing season. Besides disease-related yield loss, application of fungicides can result in considerable cost and increases the likelihood of fungicide resistance developing in the pathogen. Resistant cultivars are available for pecan growers; although, in several cases resistance has been overcome as the pathogen adapts to infect resistant hosts. Despite the importance of host resistance in scab management, there is little information regarding the molecular basis of genetic resistance to pecan scab.The purpose of this study was to elucidate mechanisms of natural pecan scab resistance by analyzing transcripts that are differentially expressed in pecan leaf samples from scab resistant and susceptible trees. The leaf samples were collected from trees in a provenance collection orchard that represents the natural range of pecan in the US and Mexico. Trees in the orchard have been exposed to natural scab infections since planting in 1989, and scab ratings were collected over three seasons. Based on this data, ten susceptible trees and ten resistant trees were selected for analysis. RNA-seq data was collected and analyzed for diseased and non-diseased parts of susceptible trees as well as for resistant trees. A total of 313 genes were found to be differentially expressed when comparing resistant and susceptible trees without disease. For susceptible samples showing scab symptoms, 1,454 genes were identified as differentially expressed compared to non-diseased susceptible samples. Many genes involved in pathogen recognition, defense responses, and signal transduction were up-regulated in diseased samples of susceptible trees, whereas differentially expressed genes in pecan scab resistant samples were generally down-regulated compared to non-diseased susceptible samples.Our results provide the first account of candidate genes involved in resistance/susceptibility to pecan scab under natural conditions in a pecan orchard. This information can be used to aid pecan breeding programs and development of biotechnology-based approaches for generating pecan cultivars with more durable scab resistance.

Multiomics analyses reveal the mechanisms of the responses of subalpine treeline trees to phenology and winter low-temperature stress.

Withstanding extreme cold stress is a prerequisite for alpine treeline trees to persist and survive. However, the underlying mechanism by which treeline trees sense phenological changes and survive hard winters has not been fully elucidated. Here, we investigated the physiology, transcriptome, and metabolome of the subalpine treeline species Larix chinensis to identify the molecular mechanism of phenological and cold resistance. Calcium and antioxidant enzyme activities (e.g., superoxide dismutase and glutathione peroxidase) are essential for coping with winter cold stress in L. chinensis. Transcriptome analysis revealed that circadian rhythm and phytohormone signalling transduction played important roles in regulating L. chinensis phenological changes and cold stress responses. The variations in the transcriptome identified were accompanied by the specific accumulation of flavones, flavonols, and monosaccharides. The flavonoid biosynthesis and phenylpropanoid biosynthesis pathways played important roles in the adaptation of L. chinensis to the extreme winter environment, and flavone and flavonol biosynthesis was an important pathway involved in bud burst. In addition, temperature and photoperiod had synergistic influences on the formation and release of bud dormancy. Thus, our findings provided new insights into the mechanism of subalpine treeline formation.

High-quality genome resource of Lasiodiplodia pseudotheobromae associated with die-back on Eucalyptus trees.

OBJECTIVES: Lasiodiplodia pseudotheobromae is an important fungal pathogen associated with die-back, canker and shoot blight in many plant hosts with a wide geographic distribution. The aim of our study was to provide high-quality genome assemblies and sequence annotation resources of L. pseudotheobromae, to facilitate future studies on the systematics, population genetics and genomics of the fungal pathogen L. pseudotheobromae. DATA DESCRIPTION: High-quality genomes of five L. pseudotheobromae isolates were sequenced based on Oxford Nanopore technology (ONT) and Illumina HiSeq sequencing platform. The total size of each assembly ranged from 43 Mb to 43.86 Mb and over 11,000 protein-coding genes were predicted from each genome. The proteins of predicted genes were annotated using multiple public databases, among the annotated protein-coding genes, more than 4,300 genes were predicted as potential virulence genes by the Pathogen Host Interactions (PHI) database. Moreover, the genome comparative analysis among L. pseudotheobromae and other closely related species revealed that 7,408 gene clusters were shared among them and 152 gene clusters unique to L. pseudotheobromae. This genome and associated datasets provided here will serve as a useful resource for further analyses of this fungal pathogen species.

Adaptive potential of maritime pine under contrasting environments.

BACKGROUND: Predicting the adaptability of forest tree populations under future climates requires a better knowledge of both the adaptive significance and evolvability of measurable key traits. Phenotypic plasticity, standing genetic variation and degree of phenotypic integration shape the actual and future population genetic structure, but empirical estimations in forest tree species are still extremely scarce. We analysed 11 maritime pine populations covering the distribution range of the species (119 families and 8 trees/family, ca. 1300 trees) in a common garden experiment planted at two sites with contrasting productivity. We used plant height as a surrogate of fitness and measured five traits (mean and plasticity of carbon isotope discrimination, specific leaf area, needle biomass, Phenology growth index) related to four different strategies (acquisitive economics, photosynthetic organ size, growth allocation and avoidance of water stress). RESULTS: Estimated values of additive genetic variation would allow adaptation of the populations to future environmental conditions. Overall phenotypic integration and selection gradients were higher at the high productivity site, while phenotypic integration within populations was higher at the low productivity site. Response to selection was related mainly to photosynthetic organ size and drought-avoidance mechanisms rather than to water use efficiency. Phenotypic plasticity of water use efficiency could be maladaptive, resulting from selection for height growth. CONCLUSIONS: Contrary to the expectations in a drought tolerant species, our study suggests that variation in traits related to photosynthetic organ size and acquisitive investment of resources drive phenotypic selection across and within maritime pine populations. Both genetic variation and evolvability of key adaptive traits were considerably high, including plasticity of water use efficiency. These characteristics would enable a relatively fast micro-evolution of populations in response to the ongoing climate changes. Moreover, differentiation among populations in the studied traits would increase under the expected more productive future Atlantic conditions.

Genomic evaluation for breeding and genetic management in Cordia africana, a multipurpose tropical tree species.

BACKGROUND: Planting tested forest reproductive material is crucial to ensure the increased resilience of intensively managed productive stands for timber and wood product markets under climate change scenarios. Single-step Genomic Best Linear Unbiased Prediction (ssGBLUP) analysis is a cost-effective option for using genomic tools to enhance the accuracy of predicted breeding values and genetic parameter estimation in forest tree species. Here, we tested the efficiency of ssGBLUP in a tropical multipurpose tree species, Cordia africana, by partial population genotyping. A total of 8070 trees from three breeding seedling orchards (BSOs) were phenotyped for height. We genotyped 6.1% of the phenotyped individuals with 4373 single nucleotide polymorphisms. The results of ssGBLUP were compared with pedigree-based best linear unbiased prediction (ABLUP) and genomic best linear unbiased prediction (GBLUP), based on genetic parameters, theoretical accuracy of breeding values, selection candidate ranking, genetic gain, and predictive accuracy and prediction bias. RESULTS: Genotyping a subset of the study population provided insights into the level of relatedness in BSOs, allowing better genetic management. Due to the inbreeding detected within the genotyped provenances, we estimated genetic parameters both with and without accounting for inbreeding. The ssGBLUP model showed improved performance in terms of additive genetic variance and theoretical breeding value accuracy. Similarly, ssGBLUP showed improved predictive accuracy and lower bias than the pedigree-based relationship matrix (ABLUP). CONCLUSIONS: This study of C. africana, a species in decline due to deforestation and selective logging, revealed inbreeding depression. The provenance exhibiting the highest level of inbreeding had the poorest overall performance. The use of different relationship matrices and accounting for inbreeding did not substantially affect the ranking of candidate individuals. This is the first study of this approach in a tropical multipurpose tree species, and the analysed BSOs represent the primary effort to breed C. africana.

The shared ancestry between the C9orf72 hexanucleotide repeat expansion and intermediate-length alleles using haplotype sharing trees and HAPTK.

The C9orf72 hexanucleotide repeat expansion (HRE) is a common genetic cause of amyotrophic lateral sclerosis (ALS) and frontotemporal dementia (FTD). The inheritance is autosomal dominant, but a high proportion of subjects with the mutation are simplex cases. One possible explanation is de novo expansions of unstable intermediate-length alleles (IAs). Using haplotype sharing trees (HSTs) with the haplotype analysis tool kit (HAPTK), we derived majority-based ancestral haplotypes of HRE samples and discovered that IAs containing ≥18-20 repeats share large haplotypes in common with the HRE. Using HSTs of HRE and IA samples, we demonstrate that the longer IA haplotypes are largely indistinguishable from HRE haplotypes and that several ≥18-20 IA haplotypes share over 5 Mb (>600 markers) haplotypes in common with the HRE haplotypes. These analysis tools allow physical understanding of the haplotype blocks shared with the majority-based ancestral haplotype. Our results demonstrate that the haplotypes with longer IAs belong to the same pool of haplotypes as the HRE and suggest that longer IAs represent potential premutation alleles.

Ancient tree genomes for old questions.

Most foundational work on the evolution and migration of plant species relies on genomic data from contemporary samples. Ancient plant samples can give us access to allele sequences and distributions on the landscape dating back to the mid Holocene or earlier (Gugerli et al., 2005). Nuclear DNA from ancient wood, however, has been mostly inaccessible until now. In a From the Cover article in this issue of Molecular Ecology, Wagner et al. (2023) present the first resequenced nuclear genomes from ancient oak wood, including two samples dated to the 15th century and one that dates to more than 3500 years ago. These ancient tree genomes open the possibility for investigating species adaptation, migration, divergence, and hybridisation in the deep past. They pave the way for what we hope will be a new era in the use of paleogenomics to study Holocene tree histories.

Efficient genomics-based 'end-to-end' selective tree breeding framework.

Since their initiation in the 1950s, worldwide selective tree breeding programs followed the recurrent selection scheme of repeated cycles of selection, breeding (mating), and testing phases and essentially remained unchanged to accelerate this process or address environmental contingencies and concerns. Here, we introduce an "end-to-end" selective tree breeding framework that: (1) leverages strategically preselected GWAS-based sequence data capturing trait architecture information, (2) generates unprecedented resolution of genealogical relationships among tested individuals, and (3) leads to the elimination of the breeding phase through the utilization of readily available wind-pollinated (OP) families. Individuals' breeding values generated from multi-trait multi-site analysis were also used in an optimum contribution selection protocol to effectively manage genetic gain/co-ancestry trade-offs and traits' correlated response to selection. The proof-of-concept study involved a 40-year-old spruce OP testing population growing on three sites in British Columbia, Canada, clearly demonstrating our method's superiority in capturing most of the available genetic gains in a substantially reduced timeline relative to the traditional approach. The proposed framework is expected to increase the efficiency of existing selective breeding programs, accelerate the start of new programs for ecologically and environmentally important tree species, and address climate-change caused biotic and abiotic stress concerns more effectively.

Tracking population structure and phenology through time using ancient genomes from waterlogged white oak wood.

Whole genome characterizations of crop plants based on ancient DNA have provided unique keys for a better understanding of the evolutionary origins of modern cultivars, the pace and mode of selection underlying their adaptation to new environments and the production of phenotypes of interest. Although forests are among the most biologically rich ecosystems on earth and represent a fundamental resource for human societies, no ancient genome sequences have been generated for trees. This contrasts with the generation of multiple ancient reference genomes for important crops. Here, we sequenced the first ancient tree genomes using two white oak wood remains from Germany dating to the Last Little Ice Age (15th century CE, 7.3× and 4.0×) and one from France dating to the Bronze Age (1700 BCE, 3.4×). We assessed the underlying species and identified one medieval remains as a hybrid between two common oak species (Quercus robur and Q. petraea) and the other two remains as Q. robur. We found that diversity at the global genome level had not changed over time. However, exploratory analyses suggested that a reduction of diversity took place at different time periods. Finally, we determined the timing of leaf unfolding for ancient trees for the first time. The study extends the application of ancient wood beyond the classical proxies of dendroclimatology, dendrochronology, dendroarchaeology and dendroecology, thereby enhancing resolution of inferences on the responses of forest ecosystems to past environmental changes, epidemics and silvicultural practices.

Application of new breeding techniques in fruit trees.

Climate change and rapid adaption of invasive pathogens pose a constant pressure on the fruit industry to develop improved varieties. Aiming to accelerate the development of better-adapted cultivars, new breeding techniques have emerged as a promising alternative to meet the demand of a growing global population. Accelerated breeding, cisgenesis, and CRISPR/Cas genome editing hold significant potential for crop trait improvement and have proven to be useful in several plant species. This review focuses on the successful application of these technologies in fruit trees to confer pathogen resistance and tolerance to abiotic stress and improve quality traits. In addition, we review the optimization and diversification of CRISPR/Cas genome editing tools applied to fruit trees, such as multiplexing, CRISPR/Cas-mediated base editing and site-specific recombination systems. Advances in protoplast regeneration and delivery techniques, including the use of nanoparticles and viral-derived replicons, are described for the obtention of exogenous DNA-free fruit tree species. The regulatory landscape and broader social acceptability for cisgenesis and CRISPR/Cas genome editing are also discussed. Altogether, this review provides an overview of the versatility of applications for fruit crop improvement, as well as current challenges that deserve attention for further optimization and potential implementation of new breeding techniques.

Loss of the accessory chromosome converts a pathogenic tree-root fungus into a mutualistic endophyte.

Some fungal accessory chromosomes (ACs) may contribute to virulence in plants. However, the mechanisms by which ACs determine specific traits associated with lifestyle transitions along a symbiotic continuum are not clear. Here we delineated the genetic divergence in two sympatric but considerably variable isolates (16B and 16W) of the poplar-associated fungus Stagonosporopsis rhizophilae. We identified a âˆ¼0.6-Mb horizontally acquired AC in 16W that resulted in a mildly parasitic lifestyle in plants. Complete deletion of the AC (Δ16W) significantly altered the fungal phenotype. Specifically, Δ16W was morphologically more similar to 16B, showed enhanced melanization, and established beneficial interactions with poplar plants, thereby acting as a dark septate endophyte. RNA sequencing (RNA-seq) analysis showed that AC loss induced the upregulation of genes related to root colonization and biosynthesis of indole acetic acid and melanin. We observed that the AC maintained a more open status of chromatin across the genome, indicating an impressive remodeling of cis-regulatory elements upon AC loss, which potentially enhanced symbiotic effectiveness. We demonstrated that the symbiotic capacities were non-host-specific through comparable experiments on Triticum- and Arabidopsis-fungus associations. Furthermore, the three isolates generated symbiotic interactions with a nonvascular liverwort. In summary, our study suggests that the AC is a suppressor of symbiosis and provides insights into the underlying mechanisms of mutualism with vascular plants in the absence of traits encoded by the AC. We speculate that AC-situated effectors and other potential secreted molecules may have evolved to specifically target vascular plants and promote mild virulence.

Wood of trees: Cellular structure, molecular formation, and genetic engineering.

Wood is an invaluable asset to human society due to its renewable nature, making it suitable for both sustainable energy production and material manufacturing. Additionally, wood derived from forest trees plays a crucial role in sequestering a significant portion of the carbon dioxide fixed during photosynthesis by terrestrial plants. Nevertheless, with the expansion of the global population and ongoing industrialization, forest coverage has been substantially decreased, resulting in significant challenges for wood production and supply. Wood production practices have changed away from natural forests toward plantation forests. Thus, understanding the underlying genetic mechanisms of wood formation is the foundation for developing high-quality, fast-growing plantation trees. Breeding ideal forest trees for wood production using genetic technologies has attracted the interest of many. Tremendous studies have been carried out in recent years on the molecular, genetic, and cell-biological mechanisms of wood formation, and considerable progress and findings have been achieved. These studies and findings indicate enormous possibilities and prospects for tree improvement. This review will outline and assess the cellular and molecular mechanisms of wood formation, as well as studies on genetically improving forest trees, and address future development prospects.

Genome sequences and population genomics provide insights into the demographic history, inbreeding, and mutation load of two 'living fossil' tree species of Dipteronia.

'Living fossils', that is, ancient lineages of low taxonomic diversity, represent an exceptional evolutionary heritage, yet we know little about how demographic history and deleterious mutation load have affected their long-term survival and extinction risk. We performed whole-genome sequencing and population genomic analyses on Dipteronia sinensis and D. dyeriana, two East Asian Tertiary relict trees. We found large-scale genome reorganizations and identified species-specific genes under positive selection that are likely involved in adaptation. Our demographic analyses suggest that the wider-ranged D. sinensis repeatedly recovered from population bottlenecks over late Tertiary/Quaternary periods of adverse climate conditions, while the population size of the narrow-ranged D. dyeriana steadily decreased since the late Miocene, especially after the Last Glacial Maximum (LGM). We conclude that the efficient purging of deleterious mutations in D. sinensis facilitated its survival and repeated demographic recovery. By contrast, in D. dyeriana, increased genetic drift and reduced selection efficacy, due to recent severe population bottlenecks and a likely preponderance of vegetative propagation, resulted in fixation of strongly deleterious mutations, reduced fitness, and continuous population decline, with likely detrimental consequences for the species' future viability and adaptive potential. Overall, our findings highlight the significant impact of demographic history on levels of accumulation and purging of putatively deleterious mutations that likely determine the long-term survival and extinction risk of Tertiary relict trees.