Close Home
iForest - Biogeosciences and Forestry
vol. 4, pp. 186-189
Copyright © 2011 by the Italian Society of Silviculture and Forest Ecology
doi: 10.3832/ifor0580-004

Collection: IUFRO RG 7.01 2010 - Antalya (Turkey)
“Adaptation of Forest Ecosystems to Air Pollution and Climate Change”
Guest Editors: Elena Paoletti, Yusuf Serengil

Technical Reports

Population genetic structure of Platanus orientalis L. in Bulgaria

M. GruevaCorresponding author, P. Zhelev


Eastern plane (Platanus orientalis L.) occurs in Europe only in the eastern part of Mediterranean region, including some remnants in Calabria and Sicily (southern Italy - [23]). Its natural range spreads eastward through Caucasus and Asia Minor to Iran ([12], [15]). P. orientalis is a fast growing species having both economical and ornamental importance. A limitation for a wider cultivation of the species in Europe is its frost sensitivity and therefore in the other parts of central and western Europe hybrid plane (London plane, Platanus × acerifolia Willd.) is preferred ([18], [14]).

In Bulgaria P. orientalis appears southwards of the 42nd latitudinal parallel and occupies mainly river banks. In contrast to other riparian species, it is very tolerant to soil conditions and can grow on very limited soil layers and on gravel along the rivers, thus depending mostly on soil humidity. Even though it is rather abundant along the rivers in the southern part of the Balkan peninsula and in the Aegean islands, the species is considered rare in Bulgaria and was included in the Red Data Book of Bulgaria ([24]). The habitats of P. orientalis (Palaearctic code 44.711) are included in the Appendix 1 of the resolution no. 4/1996 of the Permanent Committee of the Bern convention. Therefore, this species requires measures for conservation and sustainable management.

Populations of this species in Bulgaria are located in three regions over four river basins. Since river basins are separated by relatively high mountain ridges and some of the populations are also rather distant geographically, it is of particular interest to study the genetic variation and the differentiation between and within populations of P. orientalis. Both biochemical and morphological quantitative markers could be useful for such a study, since the molecular markers and quantitative traits tend to express different levels and patterns of variation. Although the use of both types of traits for population structure and differentiation stu-dies is controversial ([10], [6]), there are some examples of successful application of such an approach ([9]).

The objective of the present study was to assess the genetic variation and population structure of Platanus orientalis in Bulgaria, where it reaches the northernmost limit of its natural distribution. The knowledge of the distribution of genetic variation within and among populations could be useful for better understanding the species evolution and for designing gene conservation strategies.

Material and methods 

Plant material

The material for the study was collected from nine populations, representing the whole area of distribution of the species in Bulgaria (Tab. 1). Fifty to sixty randomly chosen trees per population were sampled. Dormant winter buds, sometimes starting to flush, were used for the analysis.

Tab. 1 - Bulgarian populations of Platanus orientalis from which samplings were collected.

Isozyme analysis

Enzymes extraction was done after grinding the bud tissue in Tris-HCl extraction buffer pH 7.3. Ten milligrams of Polyvinil-polypyrrholidone (PVPP-40) were added as a stabilizing agent in the plastic vessels, where the extraction took place. Before extraction, 15 mg dithiothreitol (DTT), 500 mg saccharose, 150 mg Polyvinil-pyrrholidone (PVP-40) and 5 mg Na2EDTA were dissolved in 15 ml extraction buffer, and 100 μl ß-mercaptoethanol were also added to the solution.

Standard 12 % starch gel electrophoresis was applied to separate the isoenzyme variants, in two buffer systems: Lithium-borate - Tris-citrate pH 8.1 discontinuous buffer system ([1]) and Tris-citrate pH 7.0 continuous buffer system ([19]). The enzyme systems and loci scored are listed in Tab. 2.

Tab. 2 - Investigated enzyme systems. (A): Lithium borate (pH 8.1) Tris Citrate (pH 8.1) discontinuous buffer system ([1]). (TC): Tris Citrate (pH 7.0) continuous buffer system ([19]).

The loci were numbered according to their position from the anode, i.e., the fastest one was numbered 1 and so forth, and the alleles were numbered according to their relative mobility relative to the most common allele (allele 100).

Data analysis

Allele frequencies were calculated from diploid genotypes. Genetic diversity within populations was characterized by the following parameters: mean number of alleles (N a); effective number of alleles (N e = 1/Σ i p i 2, where p i is the frequency of the i-th allele) for a locus and harmonic mean of all loci; percent of polymorphic loci, applying 0.05 criterion, which means that a locus is considered polymorphic if the frequency of the predominant allele does not exceed 0.95; expected heterozygosity H e = 1 - Σ i (1 - p i 2 ), where p i is the frequency of the i-th allele ([13]); observed (direct-count) heterozygosity (H o); inbreeding coefficient F = 1 - H o /H e. The software Biosys-1 ([22]) was used. Shannon index of diversity (I’) was used as an additional criterion. Genetic distances between population pairs were computed according to Nei ([13]). The differentiation was tested by means of F-statistics ([26]) following the method of Weir & Cockerham ([25]) and using the software GENEPOP v.3.1c. ([17]). While F ST was used to measure differentiation over populations, F IS and F IT were used to test the deviation from Hardy-Weinberg expectations within populations and in the total population, respectively. Multidimensional scaling (Principal Coordinate Analysis - [7]) was applied for better interpretation of genetic distances by means of the software SYN-TAX ([16]).

Results and discussion 

Polymorphism and diversity

Nine of the enzyme systems studied were polymorphic - LAP, GOT, GDH, PGI, PGM, MNR, ADH, MDH, and SKDH. The allele frequencies are available from authors upon request. The predominant allele was the same in all populations, which corresponds to the so-called “minor polymorphism” ([11]). The polymorphism and diversity parameters (Tab. 3) indicate that the level of intrapopulation genetic diversity in Platanus orientalis is well within the range of the figures reported for cross-pollinating, long-lived broadleaved species ([8]). The mean allele number per locus varied within the narrow range 2.2 to 2.3. Maximum number of alleles (5) was found in the Got-2 locus. The effective allele number (or the reciprocal of homozygosity - [2]) ranged from 1.294 (Asenovgrad) to 1.406 (Slavyanka). The Shannon’s index of diversity, borrowed from the information theory and numerical ecology ([3]) is sometimes used as an additional measure of diversity. The Got-2 locus had the highest value, with the highest allele number (Tab. 4).

Tab. 3 - Characteristics of genetic diversity and polymorphism. (N): mean number of individuals studied; (N a): mean number of alleles/locus; (N e): effective allele number (harmonic mean); (P): percent of polymorphic loci (0.05 criterion); (H o): observed heterozygosity; (H e): expected heterozygosity; (SE): standard error; (F): inbreeding coefficient (F = 1 - H o/H e); (*): significantly different from zero (p≤0.05).
Tab. 4 - Shannon index of diversity (I’).

The percent of polymorphic loci ranged from 53.8 (Melnik) to 76.9 % (Topolovo - 0.05 criterion). These results are somewhat higher than the average values summarized in Hamrick et al. ([8]) for tree species (59.5 %). Both observed and expected heterozygosities were high and the differences between them for the respective populations were low, as indicated by the inbreeding coefficient (Tab. 3). The highest heterozygosity was found in the population of Slavyanka (0.289) and the lowest one was in the population of Melnik (0.229).

Genetic differentiation among populations

Usually the genetic differentiation is measured by comparing the pairwise genetic distances ([13]) and the average indicator of the differentiation F ST (or G ST). The overall genetic differentiation (F ST) was 0.077 (Tab. 5), fairly similar to values summarized by Hamrick et al. ([8]) for species with similar life-history characteristics. High level of inbreeding was found at MNR locus, which differed from the remaining loci. The populations were most differentiated at locus Lap-2, which could be due to occurrence of silent alleles, a relatively frequent phenomenon in LAP isozymes, but undetectable when diploid tissue is analyzed.

Tab. 5 - F-statistics and putative migration. (Nm): number of migrants per generation calculated as Nm = 0.25 (1 - F ST)/ ST.

The putative number of migrants calculated by using the F ST values varied from 2 to 16 for the different loci (mean = 2.97). The average number of migrant calculated by the Private allele method ([20]) was 7.49. The differences could be due to the specific approaches for calculation or to violations of the assumptions, but in both cases this number was below 10. This fact indicates that there were mechanisms that make the isolation ineffective. This could be due to the light pollen grains of the species, and also to the fact that its fruits could be transported by the water stream at long distances. We hypothesize that even a small number of migrants is sufficient to overpower the differences among populations that could appear solely by genetic drift.

Locus MNR again showed different patterns of variation. It was the least differentiating locus, and with high inbreeding. As its visualization and interpretation was the easiest one, we can exclude possible misinterpretation of the electrophoregram as an explanation; we thus hypothesize that there were some mechanisms of selection most affecting this locus.

El-Kassaby & Yanchuk ([5]) and Davidson & El-Kassaby ([4]) modified the approach used by Slatkin ([21]) for assessing the genetic uniqueness of a population by removing it from the dataset and recalculating the overall level of differentiation among remaining populations. The procedure is repeated n times, where n is the number of populations (see the references above for details). Obviously, when a genetically most different population is out of the data set, the overall differentiation will be the lowest of all cases. In our study this population was Ivaylovgrad (results not shown).

The genetic distances were of magnitude 0.005 to 0.085, and the most different population again appeared to be Ivaylovgrad, which was also the most distant and isolated one (Tab. 6). Ivaylovgrad showed the highest genetic distances toward the other populations, while highest genetic distances among them were 0.078 (Ivaylovgrad - Goce Delchev) and 0.069 (Ivaylovgrad - Slavyanka).

Tab. 6 - Genetic distances between the population pairs.

The Multidimensional Scaling (Fig. 1) confirms these patterns of variation and differentiation and clearly shows the outstanding position of Ivaylovgrad, followed by Topolovo, Goce Delchev and Kresna.

Fig. 1 - Multidimensional scaling of the populations on the basis of the genetic distances.

The analysis of the genetic variation and population structure of Platanus orientalis in Bulgaria showed a good level of genetic diversity, both within and among populations. This fact facilitates the selection of populations for conservation. Although the decision requires deeper studies and applying different approaches, this study demonstrated the usefulness of genetic information in designing gene conservation strategies.


Ashton GC, Braden AWH (1961). Serum ß-globulin polymorphism in mice. Australian Journal of Biological Sciences 14: 248-253.
::Online::Google Scholar::
Crow JF, Kimura M (1970). An introduction to population genetics theory. Harper & Row, New York, USA, pp. 591.
::Google Scholar::
Crow JF (2001). Shannon’s brief foray into genetics. Genetics 159 (3): 915-917.
::Online::Google Scholar::
Davidson R, El-Kassaby YA (1997). Genetic diversity and gene conservation of Pacific silver fir (Abies amabilis) on Vancouver Island, British Columbia. Forest Genetics 4 (2): 85-98.
::Google Scholar::
El-Kassaby YA, Yanchuk AD (1995). Genetic variation of Pacific yew in British Columbia and its conservation. In: “Population genetics and genetic conservation of forest trees” (Baradat P, Adams WT, Mller-Starck G eds). SPB Academic Publishing bv, Hague, The Netherland, pp. 227-235.
::Google Scholar::
Felsenstein J (1986). Population differences in quantitative characters and gene frequencies: a comment on the papers by Lewontin and Rogers. American Naturalist 127 (5): 731-732.
::CrossRef::Google Scholar::
Gower JC (1966). Some distance properties of latent root and vector methods used in multivariate analysis. Biometrika 53: 325-338.
::Online::Google Scholar::
Hamrick JL, Godt MJW, Sherman-Broyles SL (1992). Factors influencing levels of genetic diversity in woody plant species. New Forests 6: 95-124.
::CrossRef::Google Scholar::
Lascoux M, Thorsén J, Gullberg U (1996). Population structure of a riparian willow species, Salix viminalis L. Genetical Research 68: 45-54.
::CrossRef::Google Scholar::
Lewontin RC (1984). Detecting population differences in quantitative characters as opposed to gene frequencies. American Naturalist 123 (1): 115-124.
::CrossRef::Google Scholar::
Lewontin RC (1985). Population genetics. Annual Review of Genetics 19: 81-102.
::CrossRef::Google Scholar::
Mossadegh A (1979). Mini-monograph on Platanus orientalis L. in Iran. Technical consultation on fast-growing plantation of broadleaved trees for Mediterranean and temperate zones. FAO-FGB-79-8/2, FAO, Rome, pp. 19.
::Google Scholar::
Nei M (1978). Estimation of average heterozygosity and genetic distance from a small number of individuals. Genetics 89 (3): 583-590.
::Online::Google Scholar::
Panetsos K, Scaltsoyannes A, Alizoti P (1994). Vegetative propagation of Platanus orientalis × P. occidentalis F1 hybrids by stem cuttings. Forest Genetics 1 (3): 125-130.
::Google Scholar::
Panetsos K (1984). Genetics and breeding of the genus Platanus. Reunion AGRIMED, Antibes (France) 29-30 October 1984.
::Google Scholar::
Podani J (1993). SYN-TAX 5.0: Computer programs for multivariate data analysis in ecology and systematics. Abstracta Botanica 17: 289-302.
::Google Scholar::
Raymond M, Rousset F (1995). GENEPOP (version 1.2): population genetics software for exact tests and ecumenicism. Journal of Heredity 86: 248-249.
::Online::Google Scholar::
Santamour FS (1970). Hybrid vigor in seedlings of re-created London plane. Nature 225: 1159-1160.
::CrossRef::Google Scholar::
Shaw CR, Prasad R (1970). Starch gel electrophoresis of enzymes - a compilation of recipes. Biochemical Genetics 4: 297-320.
::CrossRef::Google Scholar::
Slatkin M, Barton N (1989). A comparison of three indirect methods for estimating average levels of gene flow. Evolution 43: 1349-1368.
::CrossRef::Google Scholar::
Slatkin M (1985). Rare alleles as indicators of gene flow. Evolution 39: 53-65.
::CrossRef::Google Scholar::
Swofford DL, Selander RB (1989). BIOSYS-1: a computer program for the analysis of allelic variation in population genetics and biochemical systematics. Illinois Nat. Hist. Survey, Illinois, USA, pp. 43.
::Google Scholar::
Tutin TG (1964). Platanaceae. In: “Flora Europaea 1” (Tutin TG et al. eds). Cambridge University Press, Cambridge, UK, pp. 384.
::Google Scholar::
Velcev V (1984). Red data book of Bulgaria 1 plants. Bulgarian Academy of Sciences, Sofia, Bulgaria, pp. 447. [in Bulgarian].
::Google Scholar::
Weir BS, Cockerham CC (1984). Estimating F-statistics for the analysis of population structure. Evolution 38: 1358-1370.
::CrossRef::Google Scholar::
Wright S (1965). The interpretation of population structure by F-statistics with special regard to systems of mating. Evolution 19: 395-420.
::CrossRef::Google Scholar::


Grueva M, Zhelev P (2011).
Population genetic structure of Platanus orientalis L. in Bulgaria
iForest - Biogeosciences and Forestry 4: 186-189. - doi: 10.3832/ifor0580-004
First Previous Next Last
© iForest

Download Reference

Paper ID# ifor0580-004
Title Population genetic structure of Platanus orientalis L. in Bulgaria
Authors Grueva M, Zhelev P
Close Download