Data from: Genome-wide admixture and ecological niche modeling reveal the maintenance of species boundaries despite long history of interspecific gene flow

1.The maintenance of species boundaries despite interspecific gene flow has been a continuous source of interest in evolutionary biology. Many hybridizing species have porous genomes with regions impermeable to introgression, conferring reproductive barriers between species. 2.We used ecological niche modeling to study the glacial and postglacial recolonization patterns between the widely hybridizing spruce species Picea glauca and P. engelmannii in western North America. 3.Genome-wide estimates of admixture based on a panel of 311 candidate gene single nucleotide polymorphisms (SNP) from 290 genes were used to assess levels of admixture and introgression, and to identify loci putatively involved in adaptive differences or reproductive barriers between species. 4.Our paleoclimatic modeling suggests these two closely related species have a long history of hybridization and introgression, dating to at least 21,000 years ago, yet species integrity is maintained by a combination of strong environmental selection and reduced current interspecific gene flow. 5.20 loci showed evidence of divergent selection, including six loci that were both Fst outliers and associated with climatic gradients, and fourteen loci that were either outliers or showed associations with climate. These included genes responsible for carbohydrate metabolism, signal transduction, and transcription factors.Genome-wide admixture in the Picea glauca x P. engelmannii hybrid zone_311SNPsThis file contains Illumina genomic data for 311 candidate gene single nucleotide polymorphisms (SNPs) tested in 700 individuals of P.glauca, P.engelmannii and their hybrids in British Columbia (Canada). The file has a format similar to the "Structure" (Pritchard et al.2000) program, in which the first column is the genotype ID, second column is the population ID and all other columns are the SNP IDs. Please note that each SNP is described by two columns, each containing one allele. Missing values are denoted by "-9", according to the Structure format requirements.SNPs_311.txtGenome-wide admixture in the Picea glauca x P. engelmannii hybrid zone_Small data setSNPs_86.txtClimate_dataThis file contains the location (latitude, longitude and elevation) of the 771 sampled individuals (white spruce, Engelmann spruce and their hybrids). It also contains the climatic variables used in the environmental association analysis. 21 different climatic variables are included for each of the sampled individuals