Aspectual Fertility Variation and Its Effect on Gene Diversity of Seeds in Natural Stands of Taurus Cedar (Cedrus libani A. Rich.).
There are many environmental (i.e., soil properties, climate, altitude, aspect, and slope) and biological (i.e., genetic structure of population or individual plant) factors on forestry practices. Aspect defining the compass direction that a slope faces is one of the most important environmental factors because of its easy application and cheap determination for these practices such as silvicultural, watershed, and other forestry operations for managers because of its effect on soil and forest formation and habitat. Aspect can have a strong influence on temperature. This is because of the angle of the sun in the northern and southern hemispheres which is less than 90 degrees or directly overhead. For instance, aspect is a criterion in selection and establishment of seed sources and selection of gene conservation areas  and their managements. However, aspectual fertility variation and gene diversity have not been investigated or reported in the species or other plant species until now, while the environmental factors were related to fertility variation and gene diversity based on reproductive output (e.g., [2-5]). Fertility is defined as an individual ability to give progeny (i.e., reproductive success). Estimation of fertility variation among genotypes is one of the important tools used in different purposes such as gene conservation, seed production programs, managing forest genetic resources, and evolutional and genetic management of populations for plant breeding [2, 6-8] based on sustainable forestry. It was known that there were many environmental and genetical factors on reproductive characters and variation among population and within population. However, while many studies (most of them on estimation practices) have been carried out on fertility variation (e.g., [4, 9-14]), environmental studies on fertility variation and its related parameters (i.e., effective number of parents and gene diversity) are very limited [2, 3].
The main objective of the present study is to evaluate the fertility variation in selected aspectual natural stands of Taurus cedar basis of cone production over a period of three consecutive years and to estimate the effective number of parents and gene diversity in the stands. The impacts of aspect on fertility variation and gene diversity were also discussed for future forestry practices of the species.
2. Materials and Methods
Five different natural stands called also aspectual populations located at latitude between 37[degrees]49' and 37[degrees]53' N, longitude between 31[degrees] 15' and 31[degrees]20' E, and average elevation of 1600 m of Taurus cedar were sampled aspectually (flat, north, south, west, and east aspects) in Isparta region at southern part of Turkey. Numbers of mature cone were counted from the same numbered thirty trees selected phenotypically during three consecutive years (2013-2015).
Cone fertility variation ([[PSI].sub.C]) was estimated based on individual cone production as 
[mathematical expression not reproducible] (1)
where N is the census number; [Con.sub.i] is the cone fertility of the ith individual. C[V.sub.C] is the coefficient of variation in total fertility. In this paper, the fertility of ith individual was estimated by the proportion of cone production in the population. Therefore, the cone fertility represents the total contribution as zygotic parents.
The effective number of parents ([N.sub.p]) was estimated as 
[N.sub.p] = N/[[PSI].sub.C]. (2)
The relative effective number of parents ([N.sub.r]) was estimated relatively compared to the census number as [N.sub.r] = [N.sub.p]/N.
Gene diversity (GD) was estimated based on effective number of parents ([N.sub.p]) as 
GD = 1 - 0.5/[N.sub.p]. (3)
One-way ANOVA was performed in order to elucidate the differences between years as well as among stands:
[Y.sub.ij] = [mu] + [P.sub.i] + [e.sub.ij], (4)
where [Y.sub.ij] is the observation from the jth tree of the ith population/population, [mu] is overall mean, [P.sub.i] is the random effect of the ith population, and [e.sub.ij] is random error.
Phenotypic correlations ([r.sub.p])  among years for cone production were calculated using SPSS statistical package program as
[r.sub.p] = CO[V.sub.f(x,y)]/[square root of ([[sigma].sup.2.sub.f(x))] [square root of ([[sigma].sup.2.sub.f(y))]. (5)
CO[V.sub.f(x,y)] is the covariance between the characters x and y; [[sigma].sup.2.sub.f(x)] and [[sigma].sup.2.sub.f(y)] are the phenotypic variances of character x and character y, respectively.
3. Results and Discussion
3.1. Cone Production. There were large differences among aspects and years and also within aspect and year for cone production (Table 1 and Figure 1).
The averages of cone number were 19.4, 47.2, and 75.5 for the years. Third year was good seed year by average of 75.5 cones per tree of the species (Table 1). The results' well accordance with that good seed year of the species in natural stands was once in three years . The average of cone number was the highest in flat (23.5) for 2013, in south (92.1) for 2014, and in flat (95.7) for 2015, while it was lowest in south (16.3), in east (18.2), and in north (39.4) for the years, respectively (Table 1). The most abundant six trees (20% of total tree) produced 36.1%, 44.2%, 50.4%, 49.1%, and 38.3% of total cone production in the flat, south, west, east, and north populations for 2013, respectively, while it was 60.5%, 56%, 69.2%, 58%, and 41.4% for 2014 and 39.2%, 50.1%, 35.3%, 50.5%, and 52.6% for 2015. Besides, 23 individual trees did not produce any cone during observation in combined population and years; it was the highest in east population by 11 trees in polled years. The differences were also higher within aspect and year (Figure 1). For instance, coefficient of variation of cone production was the highest (105.0% and 175.2%) in west for 2013 and 2014. It was the lowest (52.7%) in the same aspect for 2015. It was calculated as 79.4%, 155.3%, and 81.0% for the years in polled populations. South aspect produced the highest cone (67.4) by the highest coefficient of variation of cone production (93.5%) for polled years. They were the lowest in north aspect (26.1% and 63.6%) (Table 1). These variations showed importance of stand/aspect and individual selection for higher cone production. Statistically significant differences (p [less than or equal to] 0.05) were found among stands for cone productions of 2014 and 2015 based on results of analysis of variance, while it was opposite for 2013. Large differences in fertility were reported among trees in natural populations of different forest tree species (e.g., [17, 18]) and among clones in seed orchards [8, 17, 19, 20]. Differences in gamete contribution could be genetic  and could be influenced by environmental factors  such as water regime of aspects  and its effect on soil characters. Environmental factors were more important than genetic constitution for the reproductive variation in Pinus sylvestris (e.g., [4, 5]). It was also found that age, elevation, and crown closure were important factors in seed yield of Pinus brutia . Altitude, age, and environmental variation, mainly in soil properties, could be effective on the observed variation in fruiting and seed set within each population in the natural forest .
There were positive and significant (p [less than or equal to] 0.01) correlations amongyears for cone production in the stands (Table 2). The results had well accordance with other forest tree species carried out on Pinus brutia .
3.2. Fertility Variation, Effective Number of Parents, and Gene Diversity. Fertility variations changed for the aspect and years (Table 3 and Figure 2). Fertility variation was 1.55 for 2013, 3.05 for 2014, 1.64 for 2015, and 1.78 for polled years, while it was 1.24, 2.16, 1.40, and 1.42 for flat and 2.01, 3.83, 1.27, and 1.47 for west aspect (Table 3 and Figure 2). Data may be much affected by west population in 2014. These results showed that seed collecting years and/or stands could be important factors to transmit gene diversity from current to next generations as also emphasized  based on sustainable forestry practices. It was suggested that the sibling coefficient ([PSI]) of natural stands as a heuristic rule of thumb could be set to three ([PSI] = 3) and that of seed orchards could be set to two ([PSI] = 2) . Estimated fertility variations in the present study could be acceptable level for typical natural populations except of west population of 2013 (3.83) and for 2014 (3.05) (Table 3). These results emphasized importance of aspect and seed collection year.
The production of cones, flowers, pollen, fruits, and seeds has been used to estimate fertility and fertility variation in many plant species [18, 26, 27]. However, data collection on cone production was easier, cheaper, and more accurate than that of strobilus count. The tree keeps the cone in longer period than strobili in a year. So, data collection period is longer in cone than in strobilus counts as also emphasized . In the present study, cone data was collected from only one population for each aspect. Therefore, it is needed to collect more data from large population on fertility variation to draw accurate conclusion.
The effective number of parents ([N.sub.p]) was mirrored to the fertility variation (Table 3). The effective number of parents was 84.3 (56.2% of census number) for combined year populations, while it was between 16.3 (54.2%, south) and 21.5 (71.8%, north). The results showed that the populations which had low effective number of parents (north and south) could be balanced to the ideal populations by traditional or genetic forest tending (e.g., removing of unproductive tree). Silvicultural practices such as natural regeneration are applied in good flowering year to obtain higher gene diversity for next generation. As presented in Table 3, gene diversity was similar (0.99) for the years in polled stands. However, it was ranged from 0.93 to 0.98 for years and stands (Table 3). Besides, while 2015 seemed to be good cone production year than 2013 (Table 1), fertility variation was similar (1.64 and 1.55) for the years (Table 3). The result implied that the gene diversity in good flowering year could be similar to the gene diversity in moderate or poor flowering years, which may be important for genetic researchers and natural regenerations of the species as also emphasized . In contrast to years, stands had large fertility variation within year. For instance, it was between 1.24 (flat) and 2.01 (west) for 2013 (Table 3). It emphasized importance of aspect in selection of gene conservation areas and seed sources. North aspect could be taken into consideration in establishment and selection of seed sources and gene conservation areas based on fertility variation and gene diversity.
Aspect is one of the most important environmental factors for forestry practices because of its easy application for managers. It is used in many purposes in forestry from selection of seed sources and nursery practice to natural regeneration or forest establishment. Results of the study also emphasized importance of aspect in forestry practices such as preparing of silvicultural planning aspect level and management of gene conservation areas in watershed.
The authors declare that there are no competing interests.
 S. Urgenc, Forest Tree Improvement, Istanbul University Press, Istanbul, Turkey, 1982.
 N. Bilir, K.-S. Kang, and D. Lindgren, "Fertility variation in six populations of Brutian pine (Pinus brutia Ten.) over altitudinal ranges," Euphytica, vol. 141, no. 1-2, pp. 163-168, 2005.
 N. Bilir, "Fertility variation in wild rose (Rosa canina) over habitat classes," International Journal of Agriculture and Biology, vol. 13, no. 1, pp. 110-114, 2011.
 N. Bilir, F. Prescher, S. Ayan, and D. Lindgren, "Growth characters and number of strobili in clonal seed orchards of Pinus sylvestris," Euphytica, vol. 152, no. 2, pp. 293-301, 2006.
 N. Bilir, F. Prescher, D. Lindgren, and J. Kroon, "Variation in cone and seed characters in clonal seed orchards of Pinus sylvestris," New Forests, vol. 36, no. 2, pp. 187-199, 2008.
 A. D. Bila, Fertility variation and its effects on gene diversity in forest tree populations [Ph.D. thesis], Swedish University of Agricultural Science, Umea, Sweden, Acta Universitatis Agriculturae Sueciae, Silvestria, 166, 2000.
 Y. A. El-Kassaby, "Evaluation of the tree-improvement delivery system: factors affecting genetic potential," Tree Physiology, vol. 15, no. 7-8, pp. 545-550, 1995.
 K.-S. Kang, A. D. Bila, A. M. Harju, and D. Lindgren, "Estimation of fertility variation in forest tree populations," Forestry, vol. 76, no. 3, pp. 329-344, 2003.
 K. S. Kang and D. Lindgren, "Fertility variation and its effect on the relatedness of seeds in Pinus densiflora, Pinus thunbergii and Pinus koraiensis clonal seed orchards," Silvae Genetica, vol. 47, no. 4, pp. 196-201, 1998.
 A. D. Bila, D. Lindgren, and T. J. Mullin, "Fertility variation and its effect on diversity over generations in a teak plantation (Tectona grandis L.f.)," Silvae Genetica, vol. 48, no. 3-4, pp. 109-114, 1999.
 N. Bilir, K.-S. Kang, and D. Lindgren, "Fertility variation and effective number in the seed production areas of Pinus radiata and Pinus pinaster," Silvae Genetica, vol. 52, no. 2, pp. 75-77, 2003.
 R. Kamalakannan, M. Varghese, N. Bilir, and D. Lindgren, "Conversion of a progeny trial of Eucalyptus tereticornis to a seedling seed orchard considering gain and fertility," in Proceedings of the Conference of Low Input Breeding and Conservation of Forest Tree Species, Antalya, Turkey, October 2006.
 R. Kamalakannan, M. Varghese, and D. Lindgren, "Fertility variation and its implications on relatedness in seed crops in seedling seed orchards of Eucalyptus camaldulensis and E. tereticornis," Silvae Genetica, vol. 56, no. 6, pp. 253-259, 2007.
 M. Varghese, A. Nicodemus, B. Nagarajan, and D. Lindgren, "Impact of fertility variation on gene diversity and drift in two clonal seed orchards of teak (Tectona grandisLinn. f.)," New Forests, vol. 31, no. 3, pp. 497-512, 2006.
 D. S. Falconer, Introduction to Quantitative Genetics, Longman Scientific Technical, Longman Group, 1989.
 M. Boydak, "Regeneration of Lebanon cedar (Cedrus libani A. Rich.) on karstic lands in Turkey," Forest Ecology and Management, vol. 178, no. 3, pp. 231-243, 2003.
 A. D. Bila and D. Lindgren, "Fertility variation in Milletia sthuhlmanii, Brachystegia spiciformis, Brachystegia bohemii and Leucaena leucocphala and its effect on relatedness in seeds," Forest Genetics, vol. 5, pp. 119-129, 1998.
 K. L. Shea, "Effects of population structure and cone production on outcrossing rates in Engelmann spruce and subalpine fir," Evolution, vol. 41, no. 1, pp. 124-136, 1987.
 N. Bilir, K. S. Kang, and H. Ozturk, "Fertility variation and gene diversity in clonal seed orchards of Pinus brutia, Pinus nigra and Pinus sylvestris in Turkey," Silvae Genetica, vol. 51, no. 2-3, pp. 112-115, 2002.
 S. Keskin, "Clonal variation in flowering and cone characteristics in a Pinus brutia seed orchard," SAFRI Technical Bulletin, vol. 9, 1999.
 G. Eriksson, A. Jonsson, and D. Lindgren, "Flowering in a clonal trial of Picea abies (Karst.)," Studia Forestalia Sueccica, Stockholm, Sweden, vol. 110, pp. 4-45, 1973.
 T. Hedegart, "Breeding systems, variation and genetic improvement of teak," in Tropical Trees: Variation, Breeding and Conservation Linnean Society Symposium Series, vol. 2, Academic Press, New York, NY, USA, 1976.
 R. Salve and B. Allen-Diaz, "Variations in soil moisture content in a rangeland catchment," Journal of Range Management, vol. 54, no. 1, pp. 44-51, 2001.
 U. Eler, "Seed yield in Calabrian Cluster pine (Pinus brutia Ten.) by age," Technical Bulletin 225, 1990.
 K. S. Kang, Genetic gain and gene diversity of seed orchard crops [Ph.D. thesis], Swedish University of Agricultural Science, Umea, Sweden, Acta Universitatis Agriculturae Sueciae, Silvestria 187, 2001.
 K. Roeder, B. Devlin, and B. G. Lindsay, "Application of maximum likelihood methods to population genetic data for the estimation of individual fertilities," Biometrics, vol. 45, no. 2, pp. 363-379, 1989.
 O. Savolainen, K. Karkkainen, A. Harju, T. Nikkanen, and M. Rusanen, "Fertility variation in Pinus sylvestris: a test of sexual allocation theory," American Journal of Botany, vol. 80, no. 9, pp. 1016-1020, 1993.
 N. Bilir and K. S. Kang, "Estimation of fertility variation by strobili and cone productions in Taurus cedar (Cedrus libani A. Rich.) populations," in Proceedings of the IUFRO Forest Tree Breeding Conference, Prague, Czech Republic, August 2014.
Nilufer Yazici and Nebi Bilir
Forestry Faculty, Suleyman Demirel University, Isparta, TR-32260, Turkey
Correspondence should be addressed to Nebi Bilir; email@example.com
Received 27 October 2016; Accepted 19 December 2016; Published 9 January 2017
Academic Editor: Pablo M. Vergara
Caption: Figure 1: Averages of cone production for aspectual populations and years.
Caption: Figure 2: Fertility variations for aspectual populations and years.
Table 1: Averages ([bar.x]) and coefficient of variation (CV%) for number of cones for the aspects and years. Year 2013 2014 Aspects [bar.x] CV% [bar.x] CV% Flat 23.5 54.0 74.1 120.2 South 16.3 78.5 92.1 113.8 West 20.3 105.0 33.0 175.2 East 16.8 92.9 18.2 117.6 North 20.2 60.4 18.8 66.5 Polled 19.4 79.4 47.2 155.3 Year 2015 Total Aspects [bar.x] CV% [bar.x] CV% Flat 95.7 64.1 64.5 66.0 South 93.6 86.1 67.4 93.5 West 85.0 52.7 46.2 69.3 East 65.7 89.2 33.5 76.4 North 39.4 97.0 26.1 63.6 Polled 75.9 81.0 47.5 88.6 Table 2: Correlation coefficients for cone production among years. [r.sub.p] * 2014 2015 Total 2013 0.570 0.408 0.867 2014 -- 0.292 0.894 2015 -- 0.489 * The correlations are significant at the 0.01 level. Table 3: Fertility variation ([PSI]), effective number of parents ([N.sub.p]), relative effective number of parents ([N.sub.r]), and gene diversity (GD) for the aspects and years. Year 2013 Aspects [PSI] [N.sub.p] [N.sub.r] GD Flat 1.24 23.4 80.7 0.98 South 1.54 18.7 64.8 0.97 West 2.01 14.5 49.9 0.97 East 1.71 16.3 58.4 0.97 North 1.31 22.2 76.7 0.98 Polled 1.55 92.5 64.7 0.99 Year 2014 Aspects [PSI] [N.sub.p] [N.sub.r] GD Flat 2.16 12.5 46.4 0.96 South 2.25 13.3 44.4 0.96 West 3.83 7.57 26.1 0.93 East 1.71 12.9 58.4 0.96 North 1.27 21.2 78.5 0.98 Polled 3.05 44.2 32.8 0.99 Year 2015 Aspects [PSI] [N.sub.p] [N.sub.r] GD Flat 1.40 21.5 71.6 0.98 South 1.72 17.5 58.3 0.97 West 1.27 23.7 78.9 0.98 East 1.71 17.0 58.5 0.97 North 1.91 15.7 52.5 0.97 Polled 1.64 90.8 60.9 0.99 Year Total Aspects [PSI] [N.sub.p] [N.sub.r] GD Flat 1.42 21.1 70.3 0.98 South 1.84 16.3 54.2 0.97 West 1.47 20.5 68.2 0.98 East 1.57 19.1 63.8 0.97 North 1.39 21.5 71.8 0.98 Polled 1.78 84.3 56.2 0.99
|Printer friendly Cite/link Email Feedback|
|Title Annotation:||Research Article|
|Author:||Yazici, Nilufer; Bilir, Nebi|
|Publication:||International Journal of Genomics|
|Date:||Jan 1, 2017|
|Previous Article:||Differential Analysis of Genetic, Epigenetic, and Cytogenetic Abnormalities in AML.|
|Next Article:||Highly Variable Genomic Landscape of Endogenous Retroviruses in the C57BL/6J Inbred Strain, Depending on Individual Mouse, Gender, Organ Type, and...|