Genomic structure and diversity of oak populations in British Parklands

Abstract

The two predominant oak species in Britain are Quercus robur (English or pedunculate oak) and Q. petraea (sessile oak). We sequenced the whole genomes of 386 oak trees from four British parkland sites and found over 50 million nuclear single nucleotide polymorphisms (SNPs), allowing us to identify 360 Q. robur, ten Q. petraea and 16 hybrid individuals using clustering methods. Comparing Q. robur and Q. petraea trees from Attingham Park, we found that the nuclear genomes of the two species are largely undifferentiated but identified 81 coding regions exhibiting strong interspecific differentiation. The nuclear genomes of our 360 Q. robur individuals showed no clear differentiation among the four parkland sites. Scans for selective sweeps in Q. robur highlighted regions containing genes with putative involvement in stress tolerance, one of which was moderately differentiated from Q. petraea. Reconstructions of past effective population sizes suggested a long population size decline in both Q. robur and Q. petraea over the Pleistocene, but population growth after the last glacial maximum. We assembled the whole chloroplast genomes of 287 Q. robur, 8 Q. petraea and 14 hybrid trees. In a phylogenetic network, these fell into five major haplotypes, which were shared among species but differed in frequency among parkland sites. We matched our chloroplast genome haplotypes to restriction enzyme fragment haplotypes identified in older studies that had surveyed ancient woodlands in Britain and much of Europe. This suggested that the parkland populations in our study derive from local seed sources.