Data from: Influence of late Quaternary climate change on present patterns of genetic variation in valley oak, Quercus lobata Née

Phylogeography and ecological niche models (ENMs) suggest that late Quaternary glacial cycles have played a prominent role in shaping present population genetic structure and diversity, but have not applied quantitative methods to dissect the relative contribution of past and present climate vs. other forces. We integrate multilocus phylogeography, climate-based ENMs and multivariate statistical approaches to infer the effects of late Quaternary climate change on contemporary genetic variation of valley oak (Quercus lobata Née). ENMs indicated that valley oak maintained a stable distribution with local migration from the last interglacial period (~120 ka) to the Last Glacial Maximum (~21 ka, LGM) to the present compared with large-scale range shifts for an eastern North American white oak (Quercus alba L.). Coast Range and Sierra Nevada foothill populations diverged in the late Pleistocene before the LGM [104 ka (28–1622)] and have occupied somewhat distinct climate niches, according to ENMs and coalescent analyses of divergence time. In accordance with neutral expectations for stable populations, nuclear microsatellite diversity positively correlated with niche stability from the LGM to present. Most strikingly, nuclear and chloroplast microsatellite variation significantly correlated with LGM climate, even after controlling for associations with geographic location and present climate using partial redundancy analyses. Variance partitioning showed that LGM climate uniquely explains a similar proportion of genetic variance as present climate (16% vs. 11–18%), and together, past and present climate explains more than geography (19%). Climate can influence local expansion–contraction dynamics, flowering phenology and thus gene flow, and/or impose selective pressures. These results highlight the lingering effect of past climate on genetic variation in species with stable distributions.Aligned, phased elongation factor 1-alpha DNA sequencesAligned, phased elongation factor 1-alpha (ef1a) DNA sequences in NEXUS format.Qlobata_ef1a_forDryad.txtAligned, phased glyceraldehyde 3-phosphate dehydrogenase DNA sequencesAligned, phased glyceraldehyde 3-phosphate dehydrogenase (g3pdh) DNA sequences in NEXUS formatQlobata_g3pdh_forDryad.txtMicrosatellite (SSR) dataSheet 1 contains genotypes for 6 chloroplast microsatellite (cpSSR) loci and Sheet 2 contains genotypes for 6 nuclear microsatellite (nSSR) loci for each valley oak (Quercus lobata) sample. GPS coordinates, sample code and population codes are also given. Sheet 3 is a readme file indicating which additional articles to cite when using the data.Gugger_etal SSR forDryad.xlsxOccurrence points for ecological niche modelsOccurrence points used to construct ecological niche models for Quercus lobata (Sheet 1) and Quercus alba (Sheet 2).Gugger_etal ENM forDryad.xls