Population and landscape genomics provide insights into the adaptive genetic variation and future climate-induced vulnerability of the endangered tree species Phoebe bournei

Abstract

Elucidating the genomic underpinnings of adaptive variation is highly important for the conservation, landscape application, and management of ornamental trees against the backdrop of global climate change. However, research on the genetic mechanisms underlying climate adaptation in Phoebe bournei-a near-threatened subtropical tree species endemic to China, which is endowed with exceptionally high ornamental and ecological value-remains scarce. Whole-genome resequencing was conducted on 362 individuals from 27 natural populations across the geographical range of the species. Genome-environment association analyses were employed to identify 1556 climate-associated variants and 167 candidate genes associated with temperature and precipitation variables. Through functional annotation and expression profiling, pivotal genes, including TRX-M4 and FBD1, were identified as integral to drought and heat stress responses, with adaptive alleles displaying distinct geographic frequency distributions and significant phenotypic differentiation. Divergent evolutionary trajectories were deduced among populations, with southeastern populations distinguished by elevated genetic diversity and strong signatures of local adaptation. Nevertheless, projections derived from the Risk of Non-Adaptedness and gradient forest models suggest that these southeastern populations will face substantial genomic offset under future climate scenarios, signaling heightened vulnerability and the need for prioritized conservation and management. This study provides the first genome-wide perspective into the adaptive evolution of P. bournei and offers a robust foundation for its conservation and climate-resilient management.

References

[1] Peng S, Ellison AM, Davis CC . Climate change intensifies plant-pollinator mismatch and increases secondary extinction risk for plants in northern latitudes. Proc Natl Acad Sci USA. 2025; 122: e2506265122

[2] Román-Palacios C, Wiens JJ . Recent responses to climate change reveal the drivers of species extinction and survival. Proc Natl Acad Sci USA. 2020; 117: 4211-7

[3] Gugger PF, Liang CT, Sork VL, et al. Applying landscape genomic tools to forest management and restoration of Hawaiian koa (Acacia koa) in a changing environment . Evol Appl. 2018; 11: 231-42

[4] Fitzpatrick MC, Keller SR . Ecological genomics meets community-level modelling of biodiversity: mapping the genomic landscape of current and future environmental adaptation. Ecol Lett. 2015; 18: 1-16

[5] Ding X, Xiao JH, Li L, et al. Congruent species delimitation of two controversial gold-thread nanmu tree species based on morphological and restriction site-associated DNA sequencing data. J Syst Evol. 2018; 57: 234-46

[6] Han X, Zhang J, Han S, et al. The chromosome-scale genome of Phoebe bournei reveals contrasting fates of terpene synthase (TPS)-a and TPS-b subfamilies . Plant Commun. 2022; 3: 100410

[7] Jiao L, Lu Y, Zhang M, et al. Ancient plastid genomes solve the tree species mystery of the imperial wood “Nanmu” in the Forbidden City, the largest existing wooden palace complex in the world. Plants People Planet. 2022; 4: 696-709

[8] Heilmayr R, Dudney J, Moore FC . Drought sensitivity in Mesic forests heightens their vulnerability to climate change. Science. 2023; 382: 1171-7

[9] Xie N, Li B, Yu J, et al. Transcriptomic and proteomic analyses uncover the drought adaption landscape of Phoebe zhennan . BMC Plant Biol. 2022; 22: 95

[10] Wang D, Li J, Tang T . Determining the optimal density of Phoebe bournei plantations based on dynamic programming under close-to-nature management measures . Sustainability. 2022; 14: 847

[11] Yang ZJ, Wu XH, Grossnickle SC, et al. Formula fertilization promotes Phoebe bournei robust seedling cultivation . Forests. 2020; 11: 781

[12] Friis G, Cotterill N, Barsoum N, et al. Genetic consequences of tree planting versus natural colonisation: implications for afforestation programmes in the United Kingdom. Evol Appl. 2025; 18: e70146

[13] Hua F, Bruijnzeel LA, Meli P, et al. The biodiversity and ecosystem service contributions and trade-offs of forest restoration approaches. Science. 2022; 376: 839-44

[14] Beale CM, Lennon JJ . Incorporating uncertainty in predictive species distribution modelling. Philos Trans R Soc Lond Ser B Biol Sci. 2012; 367: 247-58

[15] Guisan A, Thuiller W . Predicting species distribution: offering more than simple habitat models. Ecol Lett. 2005; 8: 993-1009

[16] Lind BM, Candido-Ribeiro R, Singh P, et al. How useful are genomic data for predicting maladaptation to future climate? Glob Chang Biol. 2024; 30: e17227

[17] Thuiller W, Albert C, Araújo MB, et al. Predicting global change impacts on plant species’ distributions: future challenges. Perspect Plant Ecol Evol Syst. 2008; 9: 137-52

[18] Feng L, Du FK . Landscape genomics in tree conservation under a changing environment. Front Plant Sci. 2022; 13: 822217

[19] Rellstab C, Zoller S, Walthert L, et al. Signatures of local adaptation in candidate genes of oaks (Quercus spp.) with respect to present and future climatic conditions . Mol Ecol. 2016; 25: 5907-24

[20] van Tiel N, Fopp F, Brun P, et al. Regional uniqueness of tree species composition and response to forest loss and climate change. Nat Commun. 2024; 15: 4375

[21] Wang P, Liu J, Liu B, et al. Evidence for a role of chloroplastic m-type thioredoxins in the biogenesis of photosystem II in Arabidopsis. Plant Physiol. 2013; 163: 1710-28

[22] Yanagisawa S, Schmidt RJ . Diversity and similarity among recognition sequences of Dof transcription factors. Plant J. 1999; 17: 209-14

[23] Xu K, Zhao Y, Zhao Y, et al. Soybean F-box-like protein GmFBL144 interacts with small heat shock protein and negatively regulates plant drought stress tolerance. Front Plant Sci. 2022; 13: 823529

[24] Li Q, Wang W, Wang W, et al. Wheat F-box protein gene TaFBA1 is involved in plant tolerance to heat stress. Front Plant Sci. 2018; 9: 521

[25] Wang X, Li Z, Shi Y, et al. Strigolactones promote plant freezing tolerance by releasing the WRKY41-mediated inhibition of CBF/DREB1 expression. EMBO J . 2023; 42: e112999

[26] Zhou QY, Tian AG, Zou HF, et al. Soybean WRKY-type transcription factor genes, GmWRKY13, GmWRKY21, and GmWRKY54, confer differential tolerance to abiotic stresses in transgenic Arabidopsis plants. Plant Biotechnol J. 2008; 6: 486-503

[27] Wyatt GE, Hamrick JL, Trapnell DW . The role of anthropogenic dispersal in shaping the distribution and genetic composition of a widespread north American tree species. Ecol Evol. 2021; 11: 11515-32

[28] Tian S, Kou Y, Zhang Z, et al. Phylogeography of Eomecon chionantha in subtropical China: the dual roles of the Nanling Mountains as a glacial refugium and a dispersal corridor . BMC Evol Biol. 2018; 18: 20

[29] Shi MM, Michalski SG, Welk E, et al. Phylogeography of a widespread Asian subtropical tree: genetic east-west differentiation and climate envelope modelling suggest multiple glacial refugia. J Biogeogr. 2014; 41: 1710-20

[30] Xu J, Deng M, Jiang X-L, et al. Phylogeography of Quercus glauca (Fagaceae), a dominant tree of east Asian subtropical evergreen forests, based on three chloroplast DNA interspace sequences . Tree Genet Genomes. 2014; 11: 805

[31] Cornille A, Giraud T, Smulders MJ, et al. The domestication and evolutionary ecology of apples. Trends Genet. 2013; 30: 57-65

[32] Miller AJ, Gross BL . From forest to field: perennial fruit crop domestication. Am J Bot. 2011; 98: 1389-414

[33] VanWallendael A, Soltani A, Emery NC, et al. A molecular view of plant local adaptation: incorporating stress-response networks. Annu Rev Plant Biol. 2019; 70: 559-83

[34] White NJ, Beckerman AP, Snook RR, et al. Experimental evolution of local adaptation under unidimensional and multidimensional selection. Curr Biol. 2022; 32: 1310-8.e4

[35] López-Goldar X, Agrawal AA . Ecological interactions, environmental gradients, and gene flow in local adaptation. Trends Plant Sci. 2021; 26: 796-809

[36] Kawecki TJ . Adaptation to marginal habitats. Annu Rev Ecol Evol Syst. 2008; 39: 321-42

[37] Bernatchez L, Ferchaud AL, Berger CS, et al. Genomics for monitoring and understanding species responses to global climate change. Nat Rev Genet. 2024; 25: 165-83

[38] Sang Y, Long Z, Dan X, et al. Genomic insights into local adaptation and future climate-induced vulnerability of a keystone forest tree in East Asia. Nat Commun. 2022; 13: 6541

[39] Savolainen O, Lascoux M, Merilä J . Ecological genomics of local adaptation. Nat Rev Genet. 2013; 14: 807-20

[40] Marand AP, Eveland AL, Kaufmann K, et al. Cis-regulatory elements in plant development, adaptation, and evolution. Annu Rev Plant Biol. 2023; 74: 111-37

[41] Fagny M, Austerlitz F . Polygenic adaptation: integrating population genetics and gene regulatory networks. Trends Genet. 2021; 37: 631-8

[42] Höllinger I, Pennings PS, Hermisson J . Polygenic adaptation: from sweeps to subtle frequency shifts. PLoS Genet. 2019; 15: e1008035

[43] de Miguel M, Rodríguez-Quilón I, Heuertz M, et al. Polygenic adaptation and negative selection across traits, years and environments in a long-lived plant species (Pinus pinaster Ait., Pinaceae) . Mol Ecol. 2022; 31: 2089-105

[44] Chen WJ, Chang SH, Hudson ME, et al. Contribution of transcriptional regulation to natural variations in Arabidopsis. Genome Biol. 2005; 6: R32

[45] Dubin MJ, Zhang P, Meng D, et al. DNA methylation in Arabidopsis has a genetic basis and shows evidence of local adaptation. ELife. 2015; 4: e05255

[46] Kirschner P, Perez MF, Záveská E, et al. Congruent evolutionary responses of European steppe biota to late Quaternary climate change. Nat Commun. 2022; 13: 1921

[47] Song M, Dodson J, Lu F, et al. Central China as LGM plant refugia: insights from biome reconstruction for palaeoclimate information. Sci Total Environ. 2024; 942: 173783

[48] Tierney JE, Zhu J, King J, et al. Glacial cooling and climate sensitivity revisited. Nature. 2020; 584: 569-73

[49] Qiu YX, Fu CX, Comes HP . Plant molecular phylogeography in China and adjacent regions: tracing the genetic imprints of quaternary climate and environmental change in the world’s most diverse temperate flora. Mol Phylogenet Evol. 2011; 59: 225-44

[50] Lan Y. Study on the geographical distribution changes of Chinese Nanmu in historical periods. J Chin Hist. 1995; 4: 19-32

[51] Kaufman DS, Broadman E . Revisiting the Holocene global temperature conundrum. Nature. 2023; 614: 425-35

[52] Uchiyama K, Ujino-Ihara T, Nakao K, et al. Climate-associated genetic variation and projected genetic offsets for Cryptomeria japonica D. Don under future climate scenarios . Evol Appl. 2025; 18: e70077

[53] Gougherty AV, Keller SR, Fitzpatrick MC . Maladaptation, migration and extirpation fuel climate change risk in a forest tree species. Nat Clim Chang. 2021; 11: 166-71

[54] Liu J, Zhao G, Zheng Y, et al. Genetic diversity and adaptive evolutionary history of Sapindus in China: insights from whole-genome resequencing of 100 representative individuals . Plant Biotechnol J. 2025; 23: 2485-500

[55] Bourzac K. Microbiologists at COP28 push for a seat at the climate-policy table. Nature. 2023

[56] Knight CG, Nicolitch O, Griffiths RI, et al. Soil microbiomes show consistent and predictable responses to extreme events. Nature. 2024; 636: 690-6

[57] Xie PJ, Ke YT, Kuo LY . Modified CTAB protocols for high-molecular-weight DNA extractions from ferns. Appl Plant Sci. 2023; 11: e11526

[58] Wang K, Li M, Hakonarson H . ANNOVAR: functional annotation of genetic variants from high-throughput sequencing data. Nucleic Acids Res. 2010; 38: e164

[59] Purcell S, Neale B, Todd-Brown K, et al. PLINK: a tool set for whole-genome association and population-based linkage analyses. Am J Hum Genet. 2007; 81: 559-75

[60] Alexander DH, Novembre J, Lange K . Fast model-based estimation of ancestry in unrelated individuals. Genome Res. 2009; 19: 1655-64

[61] Ringnér M . What is principal component analysis? Nat Biotechnol. 2008; 26: 303-4

[62] Minh BQ, Schmidt HA, Chernomor O, et al. IQ-TREE 2: new models and efficient methods for phylogenetic inference in the genomic era. Mol Biol Evol. 2020; 37: 1530-4

[63] Korunes KL, Samuk K . pixy: unbiased estimation of nucleotide diversity and divergence in the presence of missing data. Mol Ecol Resour. 2021; 21: 1359-68

[64] Danecek P, Auton A, Abecasis G, et al. The variant call format and VCFtools. Bioinf. 2011; 27: 2156-8

[65] Zhang C, Dong SS, Xu JY, et al. PopLDdecay: a fast and effective tool for linkage disequilibrium decay analysis based on variant call format files. Bioinf. 2019; 35: 1786-8

[66] Pickrell JK, Pritchard JK . Inference of population splits and mixtures from genome-wide allele frequency data. PLoS Genet. 2012; 8: e1002967

[67] Patton AH, Margres MJ, Stahlke AR, et al. Contemporary demographic reconstruction methods are robust to genome assembly quality: a case study in Tasmanian devils. Mol Biol Evol. 2019; 36: 2906-21

[68] Caye K, Jumentier B, Lepeule J, et al. LFMM 2: fast and accurate inference of gene-environment associations in genome-wide studies. Mol Biol Evol. 2019; 36: 852-60

[69] Gain C, François O . LEA 3: factor models in population genetics and ecological genomics with R. Mol Ecol Resour. 2021; 21: 2738-48

[70] Dixon P . VEGAN, a package of R functions for community ecology. J Veg Sci. 2003; 14: 927-30

[71] Yuan S, Shi Y, Zhou BF, et al. Genomic vulnerability to climate change in Quercus acutissima, a dominant tree species in East Asian deciduous forests . Mol Ecol. 2023; 32: 1639-55

[72] Orsini L, Vanoverbeke J, Swillen I, et al. Drivers of population genetic differentiation in the wild: isolation by dispersal limitation, isolation by adaptation and isolation by colonization. Mol Ecol. 2013; 22: 5983-99

[73] Mangiola S, Doyle MA, Papenfuss AT . Interfacing Seurat with the R tidy universe. Bioinf. 2021; 37: 4100-7

[74] Wickham H, Sievert C . ggplot2: Elegant Graphics for Data Analysis . Vol. 10. New York: Springer; 2009

[75] Yu J, Yin K, Liu Y, et al. Co-expression network analysis reveals PbTGA4 and PbAPRR2 as core transcription factors of drought response in an important timber species Phoebe bournei . Front Plant Sci. 2023; 14: 1297235

[76] Gautier M, Klassmann A, Vitalis R . REHH 2.0: a reimplementation of the R package REHH to detect positive selection from haplotype structure. Mol Ecol Resour. 2017; 17: 78-90

[77] Dong SS, He WM, Ji JJ, et al. LDBlockShow: a fast and convenient tool for visualizing linkage disequilibrium and haplotype blocks based on variant call format files. Brief Bioinform. 2021; 22: bbaa227

[78] Udvardi MK, Czechowski T, Scheible WR . Eleven golden rules of quantitative RT-PCR. Plant Cell. 2008; 20: 1736-7

[79] Fu N, Wang L, Sun Q, et al. Genome-wide identification of the bHLH transcription factor family and the regulatory roles of PbbHLH74 in response to drought stress in Phoebe bournei . Int J Biol Macromol. 2024; 283: 137760

[80] Wild J, Kopecký M, Macek M, et al. Climate at ecologically relevant scales: a new temperature and soil moisture logger for long-term microclimate measurement. Agric For Meteorol . 2019; 268: 40-7

[81] Feng LY, Smith SJ, Braun C, et al. The generation of gridded emissions data for CMIP6. Geosci Model Dev. 2020; 13: 461-82

[82] Barnosky AD, Matzke N, Tomiya S, et al. Has the earth’s sixth mass extinction already arrived? Nature. 2011; 471: 51-7

[83] Cahill AE, Aiello-Lammens ME, Fisher-Reid MC, et al. How does climate change cause extinction? Proc Biol Sci . 2013; 280: 20121890

[84] Eccles R, Zhang H, Hamilton D . A review of the effects of climate change on riverine flooding in subtropical and tropical regions. J Water Clim Change. 2019; 10: 687-707

[85] El-Barougy RF, Dakhil MA, Halmy MWA, et al. Potential extinction risk of Juniperus phoenicea under global climate change: towards conservation planning . Glob Ecol Conserv. 2023; 46: e02541

[86] England PR, Osler GHR, Woodworth LM, et al. Effects of intense versus diffuse population bottlenecks on microsatellite genetic diversity and evolutionary potential. Conserv Genet. 2003; 4: 595-604

[87] Johnson O, Ribas CC, Aleixo A, et al. Amazonian birds in more dynamic habitats have less population genetic structure and higher gene flow. Mol Ecol. 2023; 32: 2186-205

[88] Laine AL, Tylianakis JM . The coevolutionary consequences of biodiversity change. Trends Ecol Evol. 2024; 39: 745-56

[89] Leigh DM, Hendry AP, Vázquez-Domínguez E, et al. Estimated six per cent loss of genetic variation in wild populations since the industrial revolution. Evol Appl. 2019; 12: 1505-12

[90] Li S, Wei F . A revision of the name of the valuable timber nanmu. Guihaia. 1988; 8: 297-300

[91] Li T, Min X . Dormancy characteristics and germination requirements of Phoebe bournei seed . Sci Hortic. 2020; 260: 108903

[92] Norén K, Angerbjörn A . Genetic perspectives on northern population cycles: bridging the gap between theory and empirical studies. Biol Rev Camb Philos Soc. 2014; 89: 493-510

[93] Rivers M . The global tree assessment-red listing the world’s trees. BG J. 2017; 14: 16-9

[94] Shin YJ, Midgley GF, Archer ERM, et al. Actions to halt biodiversity loss generally benefit the climate. Glob Chang Biol. 2022; 28: 2846-74

[95] Song MY, Wang L, Zhang YT, et al. Temporospatial pattern of flavonoid metabolites and potential regulatory pathway of PbMYB211-coordinated kaempferol-3-O-rhamnoside biosynthesis in Phoebe bourne . Plant Physiol Biochem. 2023; 202: 107913

[96] Su S, Jin N, Nie R, et al. Effects of monochromatic light on growth and development of Phoebe bournei seedlings tend to be "extreme" . Plant Physiol Biochem. 2025; 227: 110116

[97] Urban MC . Climate change: accelerating extinction risk from climate change. Sci. 2015; 348: 571-3

[98] Walther GR, Post E, Convey P, et al. Ecological responses to recent climate change. Nature. 2002; 416: 389-95

[99] Yu H, Sui X, Sun M, et al. Relative importance of ecological, evolutionary and anthropogenic pressures on extinction risk in Chinese angiosperm genera. Front Ecol Evol. 2022; 10: 844509