New insights into adaptation and population structure of cork oak using genotyping by sequencing

Species respond to global climatic changes in a local context. Understanding this process, including its speed and intensity, is paramount due to the pace at which such changes are currently occurring. Tree species are particularly interesting to study in this regard due to their long generation times, sedentarism, and ecological and economic importance. Quercus suber L. is an evergreen forest tree species of the Fagaceae family with an essentially Western Mediterranean distribution. Despite frequent assessments of the species' evolutionary history, large-scale genetic studies have mostly relied on plastidial markers, whereas nuclear markers have been used on studies with locally focused sampling strategies. In this work, Genotyping by sequencing is used to derive 1,996 single nucleotide polymorphism markers to assess the species' evolutionary history from a nuclear DNA perspective, gain insights into how local adaptation is shaping the species' genetic background, and to forecast how Q. suber may respond to global climatic changes from a genetic perspective. Results reveal (a) an essentially unstructured species, where (b) a balance between gene flow and local adaptation keeps the species' gene pool somewhat homogeneous across its distribution, but still allowing (c) variation clines for the individuals to cope with local conditions. Risk of Non-Adaptedness (RONA) analyses suggest that for the considered variables and most sampled locations, (d) the cork oak should not require large shifts in allele frequencies to survive the predicted climatic changes. Future directions include integrating these results with ecological niche modeling perspectives, improving the RONA methodology, and expanding its use to other species. With the implementation presented in this work, the RONA can now also be easily assessed for other organisms.