Abstract
Intraspecific genetic variance and gene flow can support the adaptive evolution of species challenged by climate shifts or novel environmental conditions. Less well understood is how genome organization and gene flow interact in closely related species during evolutionary divergence and differentiation. Here we conducted genomic footprint analyses to determine how three species of Pterocarya (P . stenoptera, P. hupehensis, and P. macroptera), which are sympatric but occupy different elevational niches, adapted to the heterogeneous environment of the Qinling-Daba Mountains, China. We identified candidate genes for environmental adaptation (i.e., PIEZO1, WRKY39, VDAC3, CBL1, and RAF), and also identified regions of gene introgression between P. hupehensis and P. macroptera that show lower genetic load and higher genetic diversity than the rest of their genomes. The same introgressed regions are notably situated in areas of minimal genetic divergence yet they are characterized by elevated recombination rates. We also identified candidate genes within these introgressed regions related to environmental adaptation (TPLC2, CYCH;1, LUH, bHLH112, GLX1, TLP-3, and ABC1). Our findings have thus clarified the important role of gene flow in ecological adaptation and revealed genomic signatures of past introgression. Together, these findings provide a stronger theoretical basis for understanding the ecological adaptation and conservation of Quaternary relict woody plants in East Asia.
