The complex of beech taxa in western Eurasia (Fagus sylvatica L. s.l.) represents one of the most widespread broad-leaved trees in the temperate forests of the northern hemisphere. Two subspecies are generally recognized, F. sylvatica ssp. sylvatica in Europe and F. sylvatica ssp. orientalis in southwestern Asia, but several strongly differentiated regional populations were identified within this complex. To reveal geographical and genomic patterns of hybridization and introgression, genetic variation at 6789 single-nucleotide polymorphisms (SNP) resulting from ddRAD sequencing was studied in 288 specimens representing 207 populations distributed across the whole range. Two clusters resulting from a Bayesian analysis of population Structure well coincided with the subspecies, while increasing the number of clusters allowed to distinguish divergent lineages within subspecies. Two zones of potential hybridization between such lineages were analyzed: western Balkans at the contact of both subspecies and the Alpine area at the contact of populations coming from different glacial refugia. Geographical distribution of hybrid indices in the Balkan hybrid zone indicates asymmetric gene flow preferentially directed from F. orientalis to F. sylvatica. Eight SNPs showed a significant shift of genomic cline centers, seven of them toward excess F. orientalis ancestry, suggesting that introgression irregularly affects different regions of the beech genome. On the other hand, no genomic cline outliers were found in the Alpine hybrid zone. Implications for assisted gene flow are shortly discussed.
