In all, 10 nuclear loci were re-sequenced in four spruce species. Three of the species are boreal species with very large natural ranges: Picea mariana and P. glauca are North American, and P. abies, is Eurasian. The fourth species, P. breweriana, is a Tertiary relict from Northern California, with a very small natural range. Although the boreal species population sizes have fluctuated through the Ice Ages, P. breweriana is believed to have had a rather stable population size through the Quaternary. Indeed, the average Tajima's D was close to zero in this species and negative in the three boreal ones. Reflecting differences in current population sizes, nucleotide diversity was an order of magnitude lower in P. breweriana than in the boreal species. This is in contrast to the similar and high levels of heterozygosity observed in previous studies at allozyme loci across species. As the species have very different histories and effective population sizes, selection at allozyme loci rather than demography appears to be a better explanation for this discrepancy. Parameters of Isolation-with-Migration (IM) models were also estimated for pairs of species. Shared polymorphisms were extensive and fixed polymorphisms few. Divergence times were much shorter than those previously reported. There was also evidence of historical gene flow between P. abies and P. glauca. The latter was more closely related to P. abies than to its sympatric relative P. mariana. This last result suggests that North American and Eurasian species might have been geographically much closer in the recent past than they are today.