Understanding species response to climate change is essential for the conservation and utilization of species resources under rapid climate change in the future. In this study, the present and future suitable distribution range of Keteleeria davidiana, a tertiary relict gymnosperm, was predicted based on the maximum entropy model (MaxEnt). A total of 158 occurrence records were collected after removing the duplicated records. Six low-correlation climate variables were used to predict species distributions. The three key climate factors that affect the distribution of K. davidiana were temperature seasonality (34.96%), mean temperature of the coldest quarter (28.30%) and precipitation seasonality (13.58%). The most suitable zone of the temperature seasonality for K. davidiana was between 377.4 and 843.4. The highly suitable area was located in the mountainous regions of central and southeast China, which accounted for 13.39% of the whole study area. With climate warming in the future, the highly suitable distribution area of K. davidiana was estimated to decrease by 35% (SSP1-2.6 scenario) or 85% (SSP5-8.5 scenario). This study has provided a sufficient scientific basis for the future in situ and ex situ conservation of K. davidiana.

Quercus section Cyclobalanopsis is a dominant woody lineage in East Asian evergreen broadleaved forests. Regardless of its significant ecological and economic importance, little was known on the genome of this unique oak group. Quercus glauca, also known as ring-cup oak, is the most widespread tree in section Cyclobalanopsis. In this study, a high-quality haplotype-resolved reference genome of Q. glauca was generated from PacBio CCS and Hi-C reads. The genome size, contig N50, and scaffold N50 of Q. glauca are 902.8 Mb, 7.6 Mb, and 69.3 Mb, respectively, for haplotype1, and 913.2 Mb, 7.2 Mb, and 71.5 Mb, respectively, for haplotype2. The LTR Assembly Index of the Q. glauca genome was more than 22. A total of 37,460 and 38,312 protein-coding genes were predicted in haplotype1 and haplotype2, respectively. Homologous chromosomes of Q. glauca had excellent gene pair collinearity. A total of 905 R-genes were detected in Q. glauca. The divergence between Q. glauca and Q. suber occurred about 44.6 million years ago, agreeable to the earliest fossil records of section Cyclobalanopsis found in East Asia. The global climate change before the late Miocene and local climate change since the Pliocene were the main causes of the effective population size change of Q. glauca. The high-quality genome assembly of the most widespread species of section Cyclobalanopsis can provide essential genomic resources for unraveling the evolution mystery of the main oak lineages, and to promote the application of genomic data on the interspecific introgression, local adaptation, and speciation studies of oaks.