Differentiation and Taxonomic Identification of Roburoid Oaks in the Caucasian and Crimean Regions Using Nuclear Microsatellite Markers

The inter- and intra-specific structure of genetic variability was studied using 18 microsatellite loci (nSSRs) in closely related roburoid oaks in the Crimean-Caucasian region.The seven most widespread Quercus taxa in the region were studied in 29 morphologically pure populations from different parts of the North Caucasus, Transcaucasia, Crimea and northeastern Europe. Most taxa were studied using nSSR markers for the first time. Among the 492 trees studied, Bayesian clustering method implemented in STRUCTURE identified clusters corresponding to the pedunculate oak Quercus robur, the Hartwiss oak Q. hartwissiana, the Caucasian oak Q. macranthera, the downy oak Q. pubescens and three subspecies of sessile oak: Q. petraea ssp. petraea, Q. petraea ssp. iberica, Q. petraea ssp. medwediewii. Geographic structure was identified within Q. robur, Q. pubescens and Q. p. ssp. petraea. The 18 nSSR loci used are efficient in the taxonomic assignment of individuals, and identifying hybrids. The close relationship between the “long-pedunculate” roburoid oaks (Q. robur and Q. hartwissiana) is shown, with a greater difference from other species. For one of the subspecies of sessile oak, widespread in the North Caucasus and Crimea Q. petraea ssp. medwediewii (syn. Q. calcarea), or limestone oak, significant differences from other taxa were found, reaching the species level. The assumption of a possible hybrid origin of this taxon as a result of hybridization of Q. petraea and Q. pubescens is not confirmed by genetic analysis. The other two subspecies of Q. petraea (Q. p. ssp. petraea and Georgian oak Q. p. ssp. iberica) are differentiated to a lesser extent and are related to each other, which confirms the legitimacy of distinguishing two geographically isolated taxa at the rank of subspecies. The highest variability was observed in Q. pubescens (He = 0.777). In Q. p. ssp. medwediewii variability was lower than in other widespread taxa (He = 0.652), and was approximately at the level of variability of Q. hartwissiana (He = 0.633) and Q. macranthera (He = 0.659). Clear differentiation of taxa by nuclear markers shows the limited introgression in closely related oak species in the Caucasus and Crimea. The identified genetic clusters can be used as reference groups for further population genetic studies of oaks in the Crimean-Caucasian region.