Abstract
Understanding adaptive evolution and survival risks in understory herbs is crucial for the effective conservation of biodiversity. How environmental gradients shape species local adaptation patterns is not well understood, nor is how populations of understory herbs respond to a changing climate. In this study, we conducted population genomic analyses of Adenocaulon himalaicum (Asteraceae) with a pan-East Asian distribution, representing a good model for dominant understory herbs to elucidate adaptation mechanisms in heterogeneous forest ecosystems. Based on 34,398 putatively neutral single nucleotide polymorphisms (SNPs) across 27 populations, we identified three genetic lineages accompanied by high levels of genetic differentiation between populations. Our isolation by environment results (IBE) indicated a significant effect of environmental gradients on genomic variation of A. himalaicum (r = 0.18, p = 0.03). To decompose the relative contributions of climate, geography and population structure in explaining genetic variance, our partial RDA found that the prominent contribution of environmental effects (climatic and soil variables) explained 29% and 36% of the neutral and adaptive genetic variation, respectively. Using two genotype-environment association (GEA) methods, we identified 13 SNPs as candidates for core climate-related adaptation loci, with two of these loci further validated by qRT-PCR experiments. Projections of spatiotemporal genomic vulnerability under different future climate scenarios revealed that populations in the southeastern edge of the Himalayas, near the Sichuan Basin, the southernmost region of Northeast China and the northern Korean Peninsula, as well as northern Japan, were identified as the most vulnerable and should be prioritised for conservation. Therefore, our current study provides the genomic foundations for conservation and management strategies to elucidate how these understory herbs cope with future climate changes.