Abstract: Understanding the effects of the demographic dynamics and environmental heterogeneity on the genomic variation of forest species is important, not only for uncovering the evolutionary history of the species, but also for predicting their ability to adapt to climate change. In this study, we combined a common garden experiment with range-wide population genomics analyses to infer the demographic history and characterize patterns of local adaptation in a subtropical oak species, Quercus acutissima (Carruthers). We scanned approximately 8% of the oak genome using a balanced representation of both genic and non-genic regions and identified a total of 55 361 single nucleotide polymorphisms (SNPs) in 167 trees. Genomic diversity analyses revealed an east–west split in the species distribution range. Coalescent-based model simulations inferred a late Pleistocene divergence in Q. acutissima between the east and west groups as well as subsequent preglaciation population expansion events. Consistent with observed genetic differentiation, morphological traits also showed east–west differentiation and the biomass allocation in seedlings was significantly associated with precipitation. Environment was found to have a significant and stronger impact on the non-neutral than the neutral SNPs, and also significantly associated with the phenotypic differentiation, suggesting that, apart from the geography, environment had played a role in determining non-neutral and phenotypic variation. Our approach, which combined a common garden experiment with landscape genomics data, validated the hypothesis of local adaptation of this long-lived oak tree of subtropical China. Our study joins the small number of studies that have combined genotypic and phenotypic data to detect patterns of local adaptation.

Key words: demographic history, isolation by distance, isolation by environment, local adaptation, phenotype variation, RAD sequencing