Nucleotide variation at 12-16 nuclear loci was studied in three spruce species from the Qinghai-Tibetan Plateau (QTP), Picea likiangensis, P. wilsonii, and P. purpurea, and one species from the Tian Shan mountain range, P. schrenkiana. Silent nucleotide diversity was limited in P. schrenkiana and high in the three species from the QTP, with values higher than in boreal spruce species, despite their much more restricted distributions compared with that of the boreal species. In contrast to European boreal species that have experienced severe bottlenecks in the past, coalescent-based analysis suggests that DNA polymorphism in the species from the QTP and adjacent areas is compatible with the standard neutral model (P. likiangensis, P. wilsonii, and P. schrenkiana) or with population growth (P. purpurea). In order to test if P. purpurea is a diploid hybrid of P. likiangensis and P. wilsonii, we used a combination of approaches, including model-based inference of population structure, isolation-with-migration models, and recent theoretical results on the effect of introgression on the geographic distribution of diversity. In contrast to the three other species, each of which was predominantly assigned to a single cluster in the Structure analysis, P. purpurea individuals were scattered over the three main clusters and not, as we had expected, confined to the P. likiangensis and P. wilsonii clusters. Furthermore, the contribution of P. schrenkiana was by far the largest one. In agreement with this, the divergence between P. purpurea and P. schrenkiana was lower than the divergence of either P. likiangensis or P. wilsonii from P. schrenkiana. These results, together with previous ones showing that P. purpurea and P. wilsonii share the same haplotypes at both chloroplast and mitochondrial markers, suggest that P. purpurea has a complex origin, possibly involving additional species.