Genetic diversity and differentiation among populations of the pedunculate oak (Quercus robur) at the eastern margin of its range based on a new set of 95 SNP loci

Abstract

The aim of this study is to determine the genetic variability and differentiation among populations of pedunculate oak (Quercus robur) at the eastern margin of the species distribution range. We applied new set of 95 geographically informative nuclear SNP (single nucleotide polymorphism) loci developed using Double Digest Restriction Site Associated DNA (ddRAD), a new generation DNA sequencing technology. The study area is located in the basins of the Volga River and its two main tributaries, Kama and Belaya. Despite the strong reduction and fragmentation of pedunculate oak forests over the last several centuries and comparatively small sizes of present stands, a relatively similar genetic diversity of seven populations was observed (on average, allelic diversity υ_а = 1.472–1.603, observed heterozygosity H_O = 0.305–0.358, expected heterozygosity H_E = 0.272–0.343). The genetic distances between populations vary between 0.102 and 0.170 (on average, d₀ = 0.128). We have identified the absence of a spatial genetic structure and the existence of genetically peculiar populations in the territory studied. The genetic distances of Gregorius between populations vary between 0.102 and 0.170 (on average, d₀ = 0.128). Genetic differentiation among populations was statistically significant with a mean of the fixation index F_ST of 0.075. Populations from the Privolzhye Upplands and from northern parts of the study area make the most contribution to overall inter-population differentiation (d₀ = 0.137). No significant differences exist among the samples at lowlands along the Volga river (d₀ = 0.110). It was concluded that the multiplicity and different ages of local populations and the complexity of migration routes of the pedunculate oak from Pleistocene refugia and secondary Holocene shelters can cause this phenomenon.

Keywords

Population / Quercus robur / Gene pool / Single nucleotide polymorphism

References

[1] Aruova LB, Ospanova ZN, Aruov BB, Alibekova NT, Shashpan Z, Kyrgizbaev AT. Cyclic tests of joints of glued wooden structures. IOP Conf Ser Mat Sci Eng, 2020 829 1

[2] Blanc-Jolivet C, Bakhtina S, Yanbaev R, Yanbaev Y, Mader M, Guichoux E, Degen B. Development of new SNPs loci on Quercus robur and Quercus petraea for genetic studies covering the whole species’ distribution range. Conserv Genet Resour, 2020, 1: 1-4.

[3] Buschbom J, Yanbaev Y, Degen B. Efficient long-distance gene flow into an isolated relict oak stand. J Hered, 2011, 102(4): 464-472.

[4] Degen B, Yanbaev R, Yanbaev Y. Genetic differentiation of Quercus robur in the South-Ural. Silvae Genet, 2019, 68(1): 111-115.

[5] Ducousso A, Bordacs S (2004) Technical guidelines for genetic conservation and use: pedunculate and sessile oaks. Bioversity International, 6 p

[6] Dumolin S, Demesure B, Petit RJ. Inheritance of chloroplast and mitochondrial genomes in pedunculate oak investigated with an efficient PCR method. Theor Appl Genet, 1995, 91(8): 1253-1256.

[7] Evanno G, Regnaut S, Goudet J. Detecting the number of clusters of individuals using the software STRUCTURE: a simulation study. Mol Ecol, 2005, 14(8): 2611-2620.

[8] Gorchakovsky PL. Plants of European broad-leaved forests on the Eastern limit of their range, 1968, Sverdlovsk: Institute of Ecology of Plants and Animals 207 p (In Russian)

[9] Gregorius HR. A unique genetic distance. Biom J, 1984, 26(1): 13-18.

[10] Gregorius HR. The relationship between the concepts of genetic diversity and differentiation. Theor Appl Genet, 1987, 74(3): 397-401.

[11] Ivetić V, Devetaković J, Nonić M, Stanković D, Šijačić-Nikolić M. Genetic diversity and forest reproductive material-from seed source selection to planting. iForest, 2016, 9(5): 801-812.

[12] König AO, Ziegenhagen B, Van Dam BC, Csaikl UM, Coart E, Degen B, Burg K, de Vries SMG, Petit RJ. Chloroplast DNA variation of oaks in western Central Europe and genetic consequences of human influences. For Ecol Manag, 2002, 156(1–3): 147-166.

[13] Kopelman NM, Mayzel J, Jakobsson M, Rosenberg NA, Mayrose I. Clumpak: a program for identifying clustering modes and packaging population structure inferences across K. Mol Ecol Resour, 2015, 15(5): 1179-1191.

[14] Kozharinov AV, Borisov PV. Distribution of oak forests on the territory of the Eastern Europe for the last 13 thousand years. Lesovedenie, 2012, 5: 22-28. (In Russian)

[15] Lepoittevin C, Bodénès C, Chancerel E, Villate L, Lang T, Lesur I, Boury C, Ehrenmann F, Zelenica D, Boland A, Besse C, Garnier-Géré P, Plomion C, Kremer A. Single-nucleotide polymorphism discovery and validation in high-density SNP array for genetic analysis in European white oaks. Mol Ecol Resour, 2015, 15(6): 1446-1459.

[16] Leroy T, Rougemont Q, Dupouey JL, Bodénès C, Lalanne C, Belser C, Labadie K, Le Provost G, Aury JM, Kremer A, Plomion C. Massive postglacial gene flow between European white oaks uncovered genes underlying species barriers. New Phytol, 2020, 226(4): 1183-1197.

[17] Leummens HJL. Finlayson C, Milton G, Prentice R, Davidson N. Volga river basin (Russia). The wetland book, 2016, Dordrecht: Springer 1 11

[18] McKenney D, Pedlar J, O’Neill G. Climate change and forest seed zones: past trends, future prospects and challenges to ponder. Forest Chron, 2009, 85(2): 258-266.

[19] Moracho E, Moreno G, Jordano P, Hampe A. Unusually limited pollen dispersal and connectivity of pedunculate oak (Quercus robur) refugial populations at the species’ southern range margin. Mol Ecol, 2016, 25(14): 3319-3331.

[20] Neyshtadt MI. History of forests and paleogeography of USSR in Holocene, 1957, Moscow: Izd-vo AN SSSR 402 (In Russian)

[21] Peterson BK, Weber JN, Kay EH, Fisher HS, Hoekstra HE. Double digest RADseq: an inexpensive method for de novo SNP discovery and genotyping in model and non-model species. PLoS One, 2012 7 5

[22] Petit RJ, Brewer S, Bordács S, Burg K, Cheddadi R, Coart E, Cottrell J, Csaikl UM, van Dam B, Deans JD, Espinel S, Fineschi S, Finkeldey R, Glaz I, Goicoechea PG, Jensen JS, König AO, Lowe AJ, Kremer A. Identification of refugia and post-glacial colonisation routes of European white oaks based on chloroplast DNA and fossil pollen evidence. For Ecol Manag, 2002, 156(1–3): 49-74.

[23] Petit RJ, Csaikl UM, Bordács S, Burg K, Coart E, Cottrell J, van Dam B, Deans JD, Dumolin-Lapègue S, Fineschi S, Finkeldey R, Gillies A, Glaz I, Goicoechea PG, Jensen JS, König AO, Lowe AJ, Madsen SF, Mátyás G, Munro RC, Olalde M, Pemonge MH, Popescu F, Slade D, Tabbener H, Taurchini D, de Vries SGM, Ziegenhagen B, Kremer A. Chloroplast DNA variation in European white oaks: phylogeography and patterns of diversity based on data from over 2600 populations. For Ecol Manag, 2002, 156(1–3): 5-26.

[24] Pohjanmies T, Elshibli S, Pulkkinen P, Rusanen M, Vakkari P, Korpelainen H, Roslin T. Fragmentation-related patterns of genetic differentiation in pedunculate oak (Quercus robur) at two hierarchical scales. Silva Fenn, 2016 50 2 1510

[25] Pritchard JK, Stephens M, Donnelly P. Inference of population structure using multilocus genotype data. Genetics, 2000, 155(2): 945-959.

[26] Semerikov LF. Population structure of woody plants (on the example of oak in the European part of the USSR and the Caucasus), 1986, Moscow: Nauka 140 p (In Russian)

[27] Temunović M, Garnier-Géré P, Morić M, Franjić J, Ivanković M, Bogdan S, Hampe A. Candidate gene SNP variation in floodplain populations of pedunculate oak (Quercus robur L.) near the species’ southern range margin: weak differentiation yet distinct associations with water availability. Mol Ecol, 2020, 29(13): 2359-2378.

[28] Tsaralunga VV, Furmenkova ES, Kryukova AA. External signs of pathology of the pedunculate oak, 2015, Voronezh: Voronezh State Technical University 231 p (In Russian)

[29] Weir BS, Cockerham CC. Estimating F-statistics for the analysis of population structure. Evolution, 1984, 38: 1358-1370.

[30] Willing EM, Dreyer C, Van Oosterhout C. Estimates of genetic differentiation measured by F ST do not necessarily require large sample sizes when using many SNP markers. PLoS One, 2012 7 8

[31] Zaurbekov N, Aidosov A, Zaurbekova N, Aidosov G, Zaurbekova G, Zaurbekov I. Emission spread from mass and energy exchange in the atmospheric surface layer: two-dimensional simulation. Energ Source A, 2018, 40(23): 2832-2841.