Local adaptation and demographic history of vulnerable conifer Cephalotaxus oliveri in southern China

Abstract

Having a comprehensive understanding of genetic differentiation, responses to environmental change and demographic history is critical for genetic improvement and conservation efforts. Forest trees are an excellent resource for understanding population differentiation and adaptive genetic variation due to their ability to adapt to different climates and environments. Cephalotaxus oliveri is a relict conifer endemic to China. In this study, we generated transcriptome data and identified 17 728 high-quality single-nucleotide polymorphisms (SNPs) from 18 populations. We found significant negative correlations between expression diversity and nucleotide diversity within and among populations, suggesting that gene expression and nucleotide diversity have a reciprocal relationship when the species adapts to the environment. The analyses of population structure showed that C. oliveri displayed a striking genetic structure with four groups. BayeScEnv and RDA methods detected the signatures of local adaptation, and identified that 738 outlier SNPs were associated with precipitation, temperature and soil conditions across heterogeneous environmental conditions. Approximate Bayesian computation analyses showed that the first and second divergence occurred in the late Miocene (c. 10.075 million years ago [Ma]) and the middle Pleistocene transition (c. 0.815 Ma), respectively. Ecological niche modeling of C. oliveri revealed signs of westward expansion after the last glacial maximum, while it was predicted to experience significant range contractions in future climate change scenarios. Geographical factors and environmental factors in southern China have played a critical role in establishing the current genetic diversity and population structure of C. oliveri. This study provides an important reference for forest resource management and conservation for C. oliveri.

Keywords

Cephalotaxus oliveri / demographic history / ecological niche modeling / genetic differentiation / local adaptation / population transcriptome

References

1	CW Ahrens, M Byrne, PD Rymer. 2019. Standing genomic variation within coding and regulatory regions contributes to the adaptive capacity to climate in a foundation tree species. Molecular Ecology 28: 2502–2516.
2	DH Alexander, K Lange. 2011. Enhancements to the ADMIXTURE algorithm for individual ancestry estimation. BMC Bioinformatics 12: 246.
3	ZS An, SC Clemens, J Shen, XK Qiang, ZD Jin, YB Sun, WL Prell, JJ Luo, SM Wang, H Xu, YJ Cai, WJ Zhou, XD Liu, WG Liu, ZG Shi, LB Yan, XY Xiao, H Chang, F Wu, L Ai, FY Lu. 2011. Glacial-interglacial Indian summer monsoon dynamics. Science 333: 719–723.
4	ZS An, JE Kutzbach, WL Prell, SC Porter. 2001. Evolution of Asian monsoons and phased uplift of the Himalayan Tibetan plateau since Late Miocene times. Nature 411: 62–66.
5	B Buchfink, C Xie, DH Huson. 2015. Fast and sensitive protein alignment using DIAMOND. Nature Methods 12: 59–60.
6	PU Clark, D Archer, D Pollard, JD Blum, JA Rial, V Brovkin, AC Mix, NG Pisias, M Roy. 2006. The middle Pleistocene transition: Characteristics, mechanisms, and implications for long-term changes in atmospheric pCO2. Quaternary Science Reviews 25: 3150–3184.
7	A Conesa, S Gotz, JM Garcia-Gomez, J Terol, M Talon, M Robles. 2005. Blast2GO: A universal tool for annotation, visualization and analysis in functional genomics research. Bioinformatics 21: 3674–3676.
8	JM Cornuet, V Ravigne, A Estoup. 2010. Inference on population history and model checking using DNA sequence and microsatellite data with the software DIYABC (v1.0). BMC Bioinformatics 11: 401.
9	P Danecek, A Auton, G Abecasis, CA Albers, E Banks, MA DePristo, RE Handsaker, G Lunter, GT Marth, ST Sherry, G McVean, R Durbin, G Genomes Project Anal. 2011. The variant call format and VCFtools. Bioinformatics 27: 2156–2158.
10	JW Davey, PA Hohenlohe, PD Etter, JQ Boone, JM Catchen, ML Blaxter. 2011. Genome-wide genetic marker discovery and genotyping using next-generation sequencing. Nature Reviews Genetics 12: 499–510.
11	TP Dawson, ST Jackson, JI House, IC Prentice, GM Mace. 2011. Beyond predictions: Biodiversity conservation in a changing climate. Science 332: 53–58.
12	PB deMenocal. 1995. Plio-Pleistocene African climate. Science (Washington, DC) 270: 53–59.
13	P de Villemereuil, OE Gaggiotti. 2015. A new F-ST-based method to uncover local adaptation using environmental variables. Methods in Ecology and Evolution 6: 1248–1258.
14	FK Du, M Hou, WT Wang, KS Mao, A Hampe. 2017. Phylogeography of Quercus aquifolioides provides novel insights into the Neogene history of a major global hotspot of plant diversity in south-west China. Journal of Biogeography 44: 294–307.
15	FK Du, TR Wang, YY Wang, S Ueno, G de Lafontaine. 2020. Contrasted patterns of local adaptation to climate change across the range of an evergreen oak, Quercus aquifolioides. Evolutionary Applications 13: 2377–2391.
16	AJ Eckert, JL Wegrzyn, B Pande, KD Jermstad, JM Lee, JD Liechty, BR Tearse, KV Krutovsky, DB Neale. 2009. Multilocus patterns of nucleotide diversity and divergence reveal positive selection at candidate genes related to cold hardiness in coastal Douglas fir (Pseudotsuga menziesii var. menziesii). Genetics 183: 289–298.
17	JR Etterson. 2004. Evolutionary potential of Chamaecrista fasciculata in relation to climate change. II. Genetic architecture of three populations reciprocally planted along an environmental gradient in the great plains. Evolution 58: 1459–1471.
18	L Excoffier, HEL Lischer. 2010. Arlequin suite ver 3.5: A new series of programs to perform population genetics analyses under Linux and Windows. Molecular Ecology Resources 10: 564–567.
19	LQ Fan, HL Zheng, RI Milne, L Zhang, KS Mao. 2018. Strong population bottleneck and repeated demographic expansions of Populus adenopoda (Salicaceae) in subtropical China. Annals of Botany 121: 665–679.
20	Y Fan, CL Zhang, WD Wu, W He, L Zhang, X Ma. 2017. Analysis of genetic diversity and structure pattern of Indigofera pseudotinctoria in Karst habitats of the Wushan Mountains using AFLP markers. Molecules 22: 1734.
21	RD Finn, J Clements, SR Eddy. 2011. HMMER web server: Interactive sequence similarity searching. Nucleic Acids Research 39: W29–W37.
22	M Foll, O Gaggiotti. 2006. Identifying the environmental factors that determine the genetic structure of populations. Genetics 174: 875–891.
23	C Fu, J Lu, ZY Li, RL Liu. 2017. Analysis of the characteristics of geographical distribution and main community of species Cephalotaxus oliveri in Jiangxi. Jiangxi Science 35: 16–22.
24	LG Fu, N Li, RR Mill. 1999. Cephalotaxaceae. In: Wu ZY, Raven PH, Hong DY eds. Flora of China. Beijing: Science Press; St. Louis: Missouri Botanical Garden Press. 4: 85–88.
25	LM Gao, M Moeller, XM Zhang, ML Hollingsworth, J Liu, RR Mill, M Gibby, DZ Li. 2007. High variation and strong phylogeographic pattern among cpDNA haplotypes in Taxus wallichiana (Taxaceae) in China and North Vietnam. Molecular Ecology 16: 4684–4698.
26	DA Garfield, DE Runcie, CC Babbitt, R Haygood, WJ Nielsen, GA Wray. 2013. The impact of gene expression variation on the robustness and evolvability of a developmental gene regulatory network. PLoS Biology 11: e1001696.
27	DG Gavin, MC Fitzpatrick, PF Gugger, KD Heath, F Rodriguez-Sanchez, SZ Dobrowski, A Hampe, FS Hu, MB Ashcroft, PJ Bartlein, JL Blois, BC Carstens, EB Davis, G de Lafontaine, ME Edwards, M Fernandez, PD Henne, EM Herring, ZA Holden, W-s Kong, JQ Liu, D Magri, NJ Matzke, MS McGlone, F Saltre, AL Stigall, Y-HE Tsai, JW Williams. 2014. Climate refugia: Joint inference from fossil records, species distribution models and phylogeography. New Phytologist 204: 37–54.
28	MJS Gibson, LC Moyle. 2020. Regional differences in the abiotic environment contribute to genomic divergence within a wild tomato species. Molecular Ecology 29: 2204–2217.
29	W Gong, C Chen, C Dobes, CX Fu, MA Koch. 2008. Phylogeography of a living fossil: Pleistocene glaciations forced Ginkgo biloba L. (Ginkgoaceae) into two refuge areas in China with limited subsequent postglacial expansion. Molecular Phylogenetics and Evolution 48: 1094–1105.
30	X Gong, SS Luan, KH Hung, CC Hwang, CJ Lin, YC Chiang, TY Chiang. 2011. Population structure of Nouelia insignis (Asteraceae), an endangered species in southwestern China, based on chloroplast DNA sequences: Recent demographic shrinking. Journal of Plant Research 124: 221–230.
31	MG Grabherr, BJ Haas, M Yassour, JZ Levin, DA Thompson, I Amit, X Adiconis, L Fan, R Raychowdhury, QD Zeng, ZH Chen, E Mauceli, N Hacohen, A Gnirke, N Rhind, F di Palma, BW Birren, C Nusbaum, K Lindblad-Toh, N Friedman, A Regev. 2011. Full-length transcriptome assembly from RNA-Seq data without a reference genome. Nature Biotechnology 29(644): U130.
32	ZL Gu, SA Rifkin, KP White, WH Li. 2004. Duplicate genes increase gene expression diversity within and between species. Nature Genetics 36: 577–579.
33	XJ Guo, JR Wang. 2017. Global identification, structural analysis and expression characterization of bHLH transcription factors in wheat. BMC Plant Biology 17: 90.
34	M Ha, ED Kim, ZJ Chen. 2009. Duplicate genes increase expression diversity in closely related species and allopolyploids. Proceedings of the National Academy of Sciences of the United States of America 106: 2295–2300.
35	BJ Haas, A Papanicolaou, M Yassour, M Grabherr, PD Blood, J Bowden, MB Couger, D Eccles, B Li, M Lieber, MD MacManes, M Ott, J Orvis, N Pochet, F Strozzi, N Weeks, R Westerman, T William, CN Dewey, R Henschel, RD Leduc, N Friedman, A Regev. 2013. De novo transcript sequence reconstruction from RNA-seq using the Trinity platform for reference generation and analysis. Nature Protocols 8: 1494–1512.
36	JAP Hamlin, ML Arnold. 2015. Neutral and selective processes drive population differentiation for Iris hexagona. Journal of Heredity 106: 628–636.
37	A Hampe, RJ Petit. 2005. Conserving biodiversity under climate change: The rear edge matters. Ecology Letters 8: 461–467.
38	PW Harrison, AE Wright, JE Mank. 2012. The evolution of gene expression and the transcriptome-phenotype relationship. Seminars in Cell & Developmental Biology 23: 222–229.
39	MJ Head, PL Gibbard. 2005. Early-Middle Pleistocene transitions: An overview and recommendation for the defining boundary. Geological Society, London, Special Publications 247: 1–18.
40	AA Hoffmann, CM Sgro. 2011. Climate change and evolutionary adaptation. Nature 470: 479–485.
41	M Jakoby, B Weisshaar, W Droge-Laser, J Vicente-Carbajosa, J Tiedemann, T Kroj, F Parcy, bZIP Research Group. 2002. bZIP transcription factors in Arabidopsis. Trends in Plant Science 7: 106–111.
42	YH Ji, CK Liu, JB Landis, M Deng, JH Chen. 2021. Plastome phylogenomics of Cephalotaxus (Cephalotaxaceae) and allied genera. Annals of Botany 127: 697–708.
43	HX Jia, GJ Liu, JB Li, J Zhang, P Sun, ST Zhao, X Zhou, MZ Lu, JJ Hu. 2020a. Genome resequencing reveals demographic history and genetic architecture of seed salinity tolerance in Populus euphratica. Journal of Experimental Botany 71: 4308–4320.
44	KH Jia, W Zhao, PA Maier, XG Hu, YQ Jin, SS Zhou, SQ Jiao, YA El-Kassaby, TL Wang, XR Wang, JF Mao. 2020b. Landscape genomics predicts climate change-related genetic offset for the widespread Platycladus orientalis (Cupressaceae). Evolutionary Applications 13: 665–676.
45	IK Jordan, YI Wolf, EV Koonin. 2004. Duplicated genes evolve slower than singletons despite the initial rate increase. BMC Evolutionary Biology 4: 22.
46	J Josse, F Husson. 2016. missMDA: A package for handling missing values in multivariate data analysis. Journal of Statistical Software 70: 1–31.
47	TJ Kawecki, D Ebert. 2004. Conceptual issues in local adaptation. Ecology Letters 7: 1225–1241.
48	SA Khan, MZ Li, SM Wang, HJ Yin. 2018. Revisiting the role of plant transcription factors in the battle against abiotic stress. International Journal of Molecular Sciences 19: 1634.
49	YX Kou, SM Cheng, S Tian, B Li, DM Fan, YJ Chen, DE Soltis, PS Soltis, ZY Zhang. 2016. The antiquity of Cyclocarya paliurus (Juglandaceae) provides new insights into the evolution of relict plants in subtropical China since the late Early Miocene. Journal of Biogeography 43: 351–360.
50	YX Kou, L Zhang, DM Fan, SM Cheng, DZ Li, RGJ Hodel, ZY Zhang. 2020. Evolutionary history of a relict conifer, Pseudotaxus chienii (Taxaceae), in southeast China during the late Neogene: Old lineage, young populations. Annals of Botany 125: 105–117.
51	Z Kvacek, WC Rember. 2000. Shared miocene conifers of the Clarkia flora and Europe. Acta Universitatis Carolinae Geologica 44: 75–85.
52	V Ledent, M Vervoort. 2001. The basic helix-loop-helix protein family: Comparative genomics and phylogenetic analysis. Genome Research 11: 754–770.
53	CR Lee, T Mitchell-Olds. 2011. Quantifying effects of environmental and geographical factors on patterns of genetic differentiation. Molecular Ecology 20: 4631–4642.
54	DW Li, YF Liu, CH Zhong, HW Huang. 2010. Morphological and cytotype variation of wild kiwifruit (Actinidia chinensis complex) along an altitudinal and longitudinal gradient in central-west China. Botanical Journal of the Linnean Society 164: 72–83.
55	H Li, B Handsaker, A Wysoker, T Fennell, J Ruan, N Homer, G Marth, G Abecasis, R Durbin, P Genome Project Data. 2009. The sequence alignment/map format and SAMtools. Bioinformatics 25: 2078–2079.
56	JL Li, RI Milne, DF Ru, JB Miao, WJ Tao, L Zhang, JJ Xu, JQ Liu, KS Mao. 2020. Allopatric divergence and hybridization within Cupressus chengiana (Cupressaceae), a threatened conifer in the northern Hengduan Mountains of western China. Molecular Ecology 29: 1250–1266.
57	X Li, MH Zhao, YJ Xu, Y Li, M Tigabu, XY Zhao. 2021a. Genetic diversity and population differentiation of Pinus koraiensis in China. Horticulturae 7: 104.
58	Y Li, LC Shi, NC Pei, SA Cushman, YT Si. 2021b. Transcriptomic responses to drought stress among natural populations provide insights into local adaptation of weeping forsythia. BMC Plant Biology 21: 273.
59	N Lin, JB Landis, YX Sun, XH Huang, X Zhang, Q Liu, HJ Zhang, H Sun, HC Wang, T Deng. 2021a. Demographic history and local adaptation of Myripnois dioica (Asteraceae) provide insight on plant evolution in northern China flora. Ecology and Evolution 11: 8000–8013.
60	YP Lin, T Mitchell-Olds, CR Lee. 2021b. The ecological, genetic and genomic architecture of local adaptation and population differentiation in Boechera stricta. Proceedings of the Royal Society B: Biological Sciences 288: 20202472.
61	HJ Liu, Z Wang, YL Zhang, MH Li, T Wang, YJ Su. 2023. Geographic isolation and environmental heterogeneity contribute to genetic differentiation in Cephalotaxus oliveri. Ecology and Evolution 13: e9869.
62	HJ Liu, YL Zhang, Z Wang, YJ Su, T Wang. 2021a. Development and application of EST-SSR markers in Cephalotaxus oliveri from transcriptome sequences. Frontiers in Genetics 12: 759557.
63	JJ Liu, R Sniezko, M Murray, N Wang, H Chen, A Zamany, RN Sturrock, D Savin, A Kegley. 2016. Genetic diversity and population structure of whitebark Pine (Pinus albicaulis Engelm.) in western North America. PLoS One 11: e0167986.
64	L Liu, Z Wang, YJ Su, T Wang. 2021b. Population transcriptomic sequencing reveals allopatric divergence and local adaptation in Pseudotaxus chienii (Taxaceae). BMC Genomics 22: 388.
65	W Liu, M Kang, H Tian, HW Huang. 2013. A range wide geographic pattern of genetic diversity and population structure of Castanea mollissima populations inferred from nuclear and chloroplast microsatellites. Tree Genetics & Genomes 9: 975–987.
66	J Lopez-Pujol, FM Zhang, HQ Sun, TS Ying, S Ge. 2011. Centres of plant endemism in China: Places for survival or for speciation? Journal of Biogeography 38: 1267–1280.
67	LM Lu, LF Mao, T Yang, JF Ye, B Liu, HL Li, M Sun, JT Miller, S Mathews, HH Hu, YT Niu, DX Peng, YH Chen, SA Smith, M Chen, KL Xiang, CT Le, V-C Dang, AM Lu, PS Soltis, DE Soltis, JH Li, ZD Chen. 2018. Evolutionary history of the angiosperm flora of China. Nature 554: 234–238.
68	SR Manchester, ZD Chen, AM Lu, K Uemura. 2009. Eastern Asian endemic seed plant genera and their paleogeographic history throughout the Northern Hemisphere. Journal of Systematics and Evolution 47: 1–42.
69	PG Meirmans. 2015. Seven common mistakes in population genetics and how to avoid them. Molecular Ecology 24: 3223–3231.
70	HH Meng, T Su, XY Gao, J Li, XL Jiang, H Sun, ZK Zhou. 2017. Warm-cold colonization: Response of oaks to uplift of the Himalaya-Hengduan Mountains. Molecular Ecology 26: 3276–3294.
71	YF Miao, M Herrmann, FL Wu, XL Yan, SL Yang. 2012. What controlled Mid-Late Miocene long-term aridification in Central Asia?—Global cooling or Tibetan Plateau uplift: A review. Earth-Science Reviews 112: 155–172.
72	MW Nachman, BA Payseur. 2012. Recombination rate variation and speciation: Theoretical predictions and empirical results from rabbits and mice. Philosophical Transactions of the Royal Society B: Biological Sciences 367: 409–421.
73	NS Nagalingum, CR Marshall, TB Quental, HS Rai, DP Little, S Mathews. 2011. Recent synchronous radiation of a living fossil. Science 334: 796–799.
74	L Ni, J Schinnerl, MF Bao, BJ Zhang, J Wu, XH Cai. 2016. Two key biogenetic intermediates of Cephalotaxus alkaloids from Cephalotaxus oliveri and C. lanceolata. Tetrahedron Letters 57: 5201–5204.
75	SH Niu, J Li, WH Bo, WF Yang, A Zuccolo, S Giacomello, X Chen, FX Han, JH Yang, YT Song, YM Nie, B Zhou, PY Wang, Q Zuo, H Zhang, JJ Ma, J Wang, LJ Wang, QY Zhu, HH Zhao, ZM Liu, XM Zhang, T Liu, SR Pei, ZM Li, Y Hu, YH Yang, WZ Li, YJ Zan, LH Zhou, JX Lin, TQ Yuan, W Li, Y Li, HR Wei, HX Wu. 2022. The Chinese pine genome and methylome unveil key features of conifer evolution. Cell 185: 204–217.
76	L Orsini, J Vanoverbeke, I Swillen, J Mergeay, L De Meester. 2013. Drivers of population genetic differentiation in the wild: Isolation by dispersal limitation, isolation by adaptation and isolation by colonization. Molecular Ecology 22: 5983–5999.
77	SJ Phillips, M Dudik. 2008. Modeling of species distributions with Maxent: New extensions and a comprehensive evaluation. Ecography 31: 161–175.
78	T Qiu, LL Jiang, SZ Li, YF Yang. 2017. Small-scale habitat-specific variation and adaptive divergence of photosynthetic pigments in different alkali soils in reed identified by common garden and genetic tests. Frontiers in Plant Science 7: 2016.
79	D Ren, FC Meng, H Liu, T Xiao, JJ Lu, LG Lin, XP Chen, QW Zhang. 2018. Novel biflavonoids from Cephalotaxus oliveri Mast. Phytochemistry Letters 24: 150–153.
80	O Savolainen, M Lascoux, J Merila. 2013. Ecological genomics of local adaptation. Nature Reviews Genetics 14: 807–820.
81	O Shaul. 2002. Magnesium transport and function in plants: The tip of the iceberg. BioMetals 15: 309–323.
82	GL Shi, ZY Zhou, ZM Xie. 2010. A new Cephalotaxus and associated epiphyllous fungi from the Oligocene of Guangxi, South China. Review of Palaeobotany and Palynology 161: 179–195.
83	KM Shih, CT Chang, JD Chung, YC Chiang, SY Hwang. 2018. Adaptive genetic divergence despite significant isolation-by-distance in populations of Taiwan cow-tail fir (Keteleeria davidiana var. formosana). Frontiers in Plant Science 9: 92.
84	KB Singh, RC Foley, L Onate-Sanchez. 2002. Transcription factors in plant defense and stress responses. Current Opinion in Plant Biology 5: 430–436.
85	M Song, T Xiao, QS Wu, HF Kwok, G Chan, LG Lin, XQ Zhang, WC Ye, SJ He, QW Zhang. 2021. Biflavonoids from the twigs and leaves of Cephalotaxus oliveri Mast. and their alpha-glucosidase inhibitory activity. Natural Product Research 36: 3085–3094.
86	P Sornaraj, S Luang, S Lopato, M Hrmova. 2016. Basic leucine zipper (bZIP) transcription factors involved in abiotic stresses: A molecular model of a wheat bZIP factor and implications of its structure in function. Biochimica Et Biophysica Acta-General Subjects 1860: 46–56.
87	J Stapley, J Reger, PGD Feulner, C Smadja, J Galindo, R Ekblom, C Bennison, AD Ball, AP Beckerman, J Slate. 2010. Adaptation genomics: The next generation. Trends in Ecology & Evolution 25: 705–712.
88	S Stojnic, EV Avramidou, B Fussi, M Westergren, S Orlovic, B Matovic, B Trudic, H Kraigher, FA Aravanopoulos, M Konnert. 2019. Assessment of genetic diversity and population genetic structure of Norway spruce (Picea abies (L.) Karsten) at its southern lineage in Europe. Implications for Conservation of Forest Genetic Resources. Forests 10: 258.
89	YQ Sun, W Zhao, CQ Xu, YL Xu, YA El-Kassaby, AR De La Torre, JF Mao. 2020. Genetic variation related to high elevation adaptation revealed by common garden experiments in Pinus yunnanensis. Frontiers in Genetics 10: 1405.
90	YX Sun, MJ Moore, LL Yue, T Feng, HJ Chu, ST Chen, YH Ji, HC Wang, JQ Li. 2014. Chloroplast phylogeography of the East Asian Arcto-Tertiary relict Tetracentron sinense (Trochodendraceae). Journal of Biogeography 41: 1721–1732.
91	F Tajima. 1989. Statistical method for testing the neutral mutation hypothesis by DNA polymorphism. Genetics 123: 585–595.
92	P Tiffin, J Ross-Ibarra. 2014. Advances and limits of using population genetics to understand local adaptation. Trends in Ecology & Evolution 29: 673–680.
93	AD Twyford, CA Kidner, RA Ennos. 2015. Maintenance of species boundaries in a Neotropical radiation of Begonia. Molecular Ecology 24: 4982–4993.
94	GP Wagner, G Booth, H Bagheri-Chaichian. 1997. A population genetic theory of canalization. Evolution 51: 329–347.
95	H Walther, Z Kvacek. 2007. Early Oligocene flora of Seifhennersdorf (Saxony). Sbornik Narodniho Muzea v Praze Rada B Prirodni Vedy 63: 85–174.
96	CB Wang, T Wang, YJ Su. 2014. Phylogeography of Cephalotaxus oliveri (Cephalotaxaceae) in relation to habitat heterogeneity, physical barriers and the uplift of the Yungui Plateau. Molecular Phylogenetics and Evolution 80: 205–216.
97	IJ Wang. 2013. Examining the full effects of landscape heterogeneity on spatial genetic variation: A multiple matrix regression approach for quantifying geographic and ecological isolation. Evolution 67: 3403–3411.
98	IJ Wang, GS Bradburd. 2014. Isolation by environment. Molecular Ecology 23: 5649–5662.
99	JF Wang, YZ Pan, X Gong, YC Chiang, C Kuroda. 2011. Chloroplast DNA variation and phylogeography of Ligularia tongolensis (Asteraceae), a species endemic to the Hengduan Mountains region of China. Journal of Systematics and Evolution 49: 108–119.
100	T Wang, Z Wang, F Xia, YJ Su. 2016a. Local adaptation to temperature and precipitation in naturally fragmented populations of Cephalotaxus oliveri, an endangered conifer endemic to China. Scientific Reports 6: 25031.
101	Z Wang, Q Deng, YJ Su. 2016b. Habitat soil properties of the natural population of Pseudotaxus chienii. Ecological Science 35: 208–213.
102	R Warren, J VanDerWal, J Price, JA Welbergen, I Atkinson, J Ramirez-Villegas, TJ Osborn, A Jarvis, LP Shoo, SE Williams, J Lowe. 2013. Quantifying the benefit of early climate change mitigation in avoiding biodiversity loss. Nature Climate Change 3: 678–682.
103	JL Wegrzyn, JD Liechty, KA Stevens, L-S Wu, CA Loopstra, H Vasquez-Gross, WM Dougherty, BY Lin, JJ Zieve, PJ Martinez-Garcia, C Holt, M Yandell, A Zimin, JA Yorke, M Crepeau, D Puiu, SL Salzberg, P de Jong, K Mockaitis, D Main, CH Langley, DB Neale. 2014. Unique features of the loblolly pine (Pinus taeda L.) megagenome revealed through sequence annotation. Genetics 196: 891–909.
104	A Whitehead, DL Crawford. 2006. Neutral and adaptive variation in gene expression. Proceedings of the National Academy of Sciences of the United States of America 103: 5425–5430.
105	S Wright. 1943. Isolation by distance. Genetics 28: 114–138.
106	HH Xia, BS Wang, W Zhao, J Pan, JF Mao, XR Wang. 2018. Combining mitochondrial and nuclear genome analyses to dissect the effects of colonization, environment, and geography on population structure in Pinus tabuliformis. Evolutionary Applications 11: 1931–1945.
107	S Xiao, ZQ Mu, CR Cheng, J Ding. 2019. Three new biflavonoids from the branches and leaves of Cephalotaxus oliveri and their antioxidant activity. Natural Product Research 33: 321–327.
108	YW Xing, RH Ree. 2017. Uplift-driven diversification in the Hengduan Mountains, a temperate biodiversity hotspot. Proceedings of the National Academy of Sciences of the United States of America 114: E3444–E3451.
109	XY Xiong, JB Gou, QG Liao, YL Li, Q Zhou, GQ Bi, C Li, R Du, XT Wang, TS Sun, LJ Guo, HF Liang, PJ Lu, YY Wu, ZH Zhang, D-K Ro, Y Shang, SW Huang, JB Yan. 2021. The Taxus genome provides insights into paclitaxel biosynthesis. Nature Plants 7: 1026–1036.
110	J Xu, M Deng, XL Jiang, M Westwood, YG Song, R Turkington. 2015a. Phylogeography of Quercus glauca (Fagaceae), a dominant tree of East Asian subtropical evergreen forests, based on three chloroplast DNA interspace sequences. Tree Genetics & Genomes 11: 805.
111	Q Xu, SL Xing, CY Zhu, W Liu, YY Fan, Q Wang, ZH Song, WH Yang, F Luo, F Shang, LF Kang, WL Chen, J Yan, JQ Li, T Sang. 2015b. Population transcriptomics reveals a potentially positive role of expression diversity in adaptation. Journal of Integrative Plant Biology 57: 284–299.
112	HF Yan, CY Zhang, FY Wang, CM Hu, XJ Ge, G Hao. 2012. Population expanding with the phalanx model and lineages split by environmental heterogeneity: A case study of Primula obconica in subtropical China. PLoS One 7: e41315.
113	J Yan, ZH Song, Q Xu, LF Kang, CY Zhu, SL Xing, W Liu, J Greimler, T Zuest, JQ Li, T Sang. 2017. Population transcriptomic characterization of the genetic and expression variation of a candidate progenitor of Miscanthus energy crops. Molecular Ecology 26: 5911–5922.
114	J Yang, SH Lee, ME Goddard, PM Visscher. 2011. GCTA: A tool for genome-wide complex trait analysis. American Journal of Human Genetics 88: 76–82.
115	ZH Yang, XD Lang, SF Li, JR Su. 2015. Niche and interspecific interaction of dominant tree population in Cephalotaxus oliveri communities. Forest Research 28: 473–478.
116	JB Yoder, J Stanton-Geddes, P Zhou, R Briskine, ND Young, P Tiffin. 2014. Genomic signature of adaptation to climate in Medicago truncatula. Genetics 196: 1263–1275.
117	JC Zachos, GR Dickens, RE Zeebe. 2008. An early Cenozoic perspective on greenhouse warming and carbon-cycle dynamics. Nature 451: 279–283.
118	JW Zhang, A D'Rozario, XQ Liang, ZK Zhou. 2019. Middle Miocene Cephalotaxus (Taxaceae) from Yunnan, Southwest China, and its implications to taxonomy and evolution of the genus. Palaeoworld 28: 381–402.
119	L Zhang, WQ Su, R Tao, WY Zhang, JJ Chen, PY Wu, CH Yan, Y Jia, RM Larkin, D Lavelle, MJ Truco, SR Chin-Wo, RW Michelmore, HH Kuang. 2017. RNA sequencing provides insights into the evolution of lettuce and the regulation of flavonoid biosynthesis. Nature Communications 8: 2264.
120	S Zhu, J Chen, J Zhao, HP Comes, P Li, C Fu, X Xie, R Lu, W Xu, Y Feng, W Ye, S Sakaguchi, Y Isagi, L Li, M Lascoux, Y Qiu. 2020. Genomic insights on the contribution of balancing selection and local adaptation to the long-term survival of a widespread living fossil tree, Cercidiphyllum japonicum. New Phytologist 228: 1674–1689.
121	A Zimin, KA Stevens, M Crepeau, A Holtz-Morris, M Koriabine, G Marcais, D Puiu, M Roberts, JL Wegrzyn, PJ de Jong, DB Neale, SL Salzberg, JA Yorke, CH Langley. 2014. Sequencing and assembly of the 22-Gb loblolly pine genome. Genetics 196: 875–890.