Abstract
Leaf variations of Betula medwediewii in Türkiye with respect to its distribution pattern were studied because this plant taxon is a Euxine element and a distinct relict species in Türkiye. In this context, the aim of this study is to investigate variations in leaf morphological characteristics within and among B. medwediewii populations in natural distribution areas in Türkiye. Thus, the morphological leaf characteristics of four populations of B. medwediewii growing at different elevations spanning from 1472 to 2065 m.a.s.l in two cities (Artvin and Rize) and four boroughs (Arhavi, Borçka, Murgul, and Çamlıhemşin) of Türkiye were analyzed. The four populations were selected based on their natural distribution and 1,200 leaf samples belonging to 40 individuals were measured with ImageJ. According to the results of correlation analysis, statistically significant relationships were determined among morphological leaf characteristics. Mean values for petiole length (1.16 cm), leaf width (5.11 cm), leaf length (7.96 cm), length of lamina (6.79 cm), leaf area (26.62 cm²), leaf vein angle (54.21°), and number of primary leaf veins (20.03) were determined in all populations. Because B. medwediewii contributes to biodiversity and can preserve the ecological stability of the Turkish forest area, its conservation is crucial.
Download PDF
Full Article
Leaf Morphological Variation in Betula medwediewii Regel Populations from Türkiye
Arzu Ergül Bozkurt *
Leaf variations of Betula medwediewii in Türkiye with respect to its distribution pattern were studied because this plant taxon is a Euxine element and a distinct relict species in Türkiye. In this context, the aim of this study is to investigate variations in leaf morphological characteristics within and among B. medwediewii populations in natural distribution areas in Türkiye. Thus, the morphological leaf characteristics of four populations of B. medwediewii growing at different elevations spanning from 1472 to 2065 m.a.s.l in two cities (Artvin and Rize) and four boroughs (Arhavi, Borçka, Murgul, and Çamlıhemşin) of Türkiye were analyzed. The four populations were selected based on their natural distribution and 1,200 leaf samples belonging to 40 individuals were measured with ImageJ. According to the results of correlation analysis, statistically significant relationships were determined among morphological leaf characteristics. Mean values for petiole length (1.16 cm), leaf width (5.11 cm), leaf length (7.96 cm), length of lamina (6.79 cm), leaf area (26.62 cm²), leaf vein angle (54.21°), and number of primary leaf veins (20.03) were determined in all populations. Because B. medwediewii contributes to biodiversity and can preserve the ecological stability of the Turkish forest area, its conservation is crucial.
DOI: 10.15376/biores.19.3.4871-4885
Keywords: Betula medwediewii; Leaf morphology; Relict species
Contact information: Artvin Çoruh University, Faculty of Forestry, Department of Forest Engineering, Artvin, Türkiye; *Corresponding author: arzu.ergulbozkurt@artvin.edu.tr
INTRODUCTION
Betula medwediewii Regel (family Betulaceae), which is known as Transcaucasian birch, is naturally distributed in four provinces: Rize and Artvin cities in north-eastern Türkiye, Adjara in south-western Georgia, and on Mount Jvari in western Georgia. This species usually occurs in subalpine mixed forests of Picea orientalis (L.) and in Rhododendron scrubs at an elevation of 1300 to 2160 m.a.s.l. It is a Euxine species, which means Western Black Sea Subregion (Euxine) of the Euro-Siberian Flora Area, and a distinct relict species in Türkiye (Browicz 1972; Davis 1988). This plant taxon is also rare in Georgia, where populations are very small, fragmented, and unstable. Its presence in Iran, Azerbaijan, and Armenia has not yet been confirmed (Shaw et al. 2014). On the Red List of Endangered Species, B. medwediewii is categorized as DD (Data Deficient) based on IUCN (International union for the conservation of nature) criteria (Wilson et al. 2018). According to Davis (1988) Artvin-Hatila valley is one of the natural distribution areas of B. medwediewii in Türkiye. However, during extensive field work between 2019 and 2020 years, this species was not encountered in the Artvin-Hatila valley.
Morphological differences in leaves can be associated with geographical region, climatic condition, and elevation (Gratani et al. 2012). Artvin is surrounded by high mountains and forests. In this region, the Eastern Black Sea oceanic rainfall regime is seen in the Arhavi district, and Arhavi is 30 km from Artvin. The elevation of the area ranges between 70 and 3000 asl. Camili is in the Borçka district, 30 km from Artvin, with an elevation ranging from 350 to 3440 asl, and it is under the influence of the Eastern Black Sea climate. This research area borders Georgia. Murgul is 48 km away from Artvin at an elevation of 250 to 3977 asl, and with abundant rainfall. Çamlıhemşin is 61 km away from Rize at an elevation of 295 to 3200 asl (Buser and Cvetic 1972; Akman 1999; Batan and Özdemir 2013; Yener et al. 2017), and the climate is characterized by abundant rain in each season.
Various research studies have indicated that a broad variety of plants’ ability to grow and survive in competition depends on leaf variation (Raschke 1960; Parkhurst and Loucks 1972; Givnish 1979; Hinckley et al. 1989; Gurevitch 1992; Calagari et al. 2006). Many studies on leaf morphology have shown that the shapes of leaves are of great importance for plant growth and are highly influenced by environmental conditions (Coleman et al. 1994; Pyakurel and Wang 2013). The main aim of this study was to analyze the leaf morphological characteristics of four populations of B. medwediewii growing at different elevations with different climatic conditions in two cities (Artvin and Rize) and four boroughs (Arhavi, Borçka, Murgul, and Çamlıhemşin) of Türkiye. The results are expected to show whether elevational and regional variations in these traits indicate future adaptive potential. At the same time, the results of the study will be instructive in terms of individuals that will outperform in the productive generation.
EXPERIMENTAL
Study Area
The four populations of B. medwediewii in Artvin-Arhavi (AA), Artvin-Borçka (Camili) (ABC), Artvin-Murgul (AMU), and Rize-Çamlıhemşin (RÇA) distributed in Türkiye were selected based on their natural distribution.
Fig. 1. Geographical locations of B. medwediewii populations (AAR (Artvin-Arhavi), ABC (Artvin-Borçka (Camili), AMU (Artvin-Murgul), and RÇA (Rize-Çamlıhemşin))
Leaf samples were taken from a limited number of random individuals, in the other natural distribution areas (Artvin-Arhavi, Artvin-Borçka (Camili), Artvin-Murgul, and Rize-Çamlıhemşin). B. medwediewii is in A8 according to the grid system adopted in the Flora of Turkey (Browicz 1972; Davis 1988) Türkiye is located between 26°- 45° east longitudes and 36°- 42° north latitudes (Fig. 1).
The climate types of the study area were determined according to the Thornthwaite Climate Classification (Thornthwaite 1948), which provides detailed results by considering the water storage capacity of the soil (Akman 2011). Long-term meteorological data were evaluated (2008 to 2020 year), and climate diagrams were created by using Thornthwaite’s system.
The city of Artvin is a semi-moist area with medium temperatures (mesothermal); water deficiency is moderate and mainly occurs in summer between July and September, whereas January, February, March, and December have excessive rainfall. The city of Rize is a very moist area with medium temperatures (mesothermal) with little or no water deficit and an ocean climate (Fig. 2).
Fig. 2. Climate diagrams of Artvin and Rize according to the Thornthwaite method
Table 1. Information on the Four Betula medwediewii Populations (Artvin-Arhavi (AAR), Artvin-Borçka (Camili) (ABC), Artvin-Murgul (AMU), Rize-Çamlıhemşin (RÇA) North (N), Northeast (NE), Southwest (SW))
Morphological Leaf Measurements
The morphological leaf characteristics petiole length (PL) (cm), leaf width (LW) (cm), leaf length (LL) (cm), length of lamina (LLA) (cm), leaf area (LA) (cm2), leaf vein angle (LVA) (°), and number of primary leaf veins (NLV) were determined using ImageJ (Image Analysis Software). From the four selected populations, a total of 1,200 leaf samples, belonging to 40 individuals, were measured (Table 1). Leaf vein angle was obtained by measuring the middle part of each leaf (Fig. 3).
In total, 40 individuals and 1,200 leaf samples (300 from AA, 300 from ABC, 300 from AMU, and 300 from RÇA) were selected (Table 1) in 2019 and 2020. The crown was visually divided into four different parts. Mature leaves were collected in the middle of the upper crown of the tree. Solely fully developed leaves of trees, without signs of damage from pest, disease or wilting were measured (Kozlov et al. 2021). Leaves were subjected to a preservation treatment that included pressing and drying (Tomaszewski and Górzkowska 2016). After that, leaf images were taken using a scanner and a ruler. Leaf traits were measured using ImageJ analysis software package (Schneider et al. 2012), calibrated by the ruler for accuracy (Mhamdi et al. 2013; Yildirim and Turna 2021) (Fig. 3).
Fig. 3. (a) General view of Betula medwediewii, (b) leaf samples, (c) measurement of leaf length, (d) leaf width, (e) length of petiole, (f) leaf vein angle, (g) length of lamina, and (h) leaf area, using ImageJ
Statistical Analysis
Data were analyzed using the statistical software package SPSS 20.0 software package (IBM, Armonk, NY, USA). Data were submitted to ANOVA (Analysis of Variance), P < 0.05, Duncan’s test, and cluster analysis. A variance analysis (one-way ANOVA) was used to determine variations among and within populations. Duncan’s test was performed to determine the groups that were found in terms of populations for all leaf traits. A cluster analysis using the Euclidean method was performed to examine leaf variation within and between populations. Pearson’s correlation was carried out to test for correlations between leaf morphological characteristics and altitude.
RESULTS
Variations in leaf morphology were observed among the populations. Figure 4 shows the mean values and standard deviations for the measured parameters. There were statistically significant differences (P < 0.01) among populations for all the measured characteristics of the leaves.
Based on the Duncan’s test, there were three groups in terms of petiole length, leaf length, leaf width, leaf area, and leaf vein angle; two groups in terms of length of lamina; and four groups in terms of the number of primary leaf veins. Among the four populations, the highest values of petiole length, leaf length, and length of lamina were obtained in the population of Rize-Çamlıhemşin, while the lowest values were obtained in the population of Artvin-Arhavi. The Artvin-Murgul population had the highest values of leaf width, leaf area, and leaf vein angle, while the lowest value was obtained in the populations of Rize-Çamlıhemşin and Artvin-Arhavi. The Artvin-Borçka (Camili) population had the highest number of primary leaf veins, while the lowest value was obtained in the Artvin-Arhavi population (Fig. 4).
Fig. 4. ANOVA and Duncan test results of the average petiole length (PL), leaf length (LL), length of lamina (LLA), leaf width (LW), leaf area (LA), leaf vein angle (LVA), and number of primary leaf veins (NLV for the four Betula medwediewii populations
There were two groups based on PL, LA, and NLV and three groups based on LVA in the Artvin-Arhavi population. There were six groups in terms of LW, LA, and NLV, five groups in terms of PL, LL, and LLA, and four groups in terms of LVA in the Artvin-Borçka (Camili) population. There were three groups based on LVA and NLV in the Artvin-Murgul population. There were three groups based on LL, LLA, and LVA and two groups based on LA in the Rize-Çamlıhemşin population (Table 2).
As a result of the variance analysis of the Artvin-Arhavi population, there were significant differences in terms of petiole length, leaf area, leaf vein angle, and number of primary leaf veins (α < 0.05), while there was no significant difference in leaf length, length of lamina, and leaf width. In the Artvin-Borcka (Camili) population, there was a significant difference based on all leaf characteristics. In the Artvin-Murgul population, there was a significant difference in terms of leaf vein angle and number of primary leaf veins. As a result of the ANOVA in the Rize-Çamlıhemşin populations, there was a significant difference in leaf length, length of lamina, leaf area, and leaf vein angle, while there was no significant difference in terms of petiole length, leaf width, and number of primary leaf veins (Table 2).
Table 2. Average Values and ANOVA Results for Petiole Length (PL), Leaf Width (LW), Leaf Length (LL), Length of lamina (LLA), Leaf Area (LA), Leaf Vein Angle (LVA) and Number of Primary Leaf Veins (NLV) Values within Populations
* Significance level (P) < 0.05 statistically difference
As a result of the study, the mean values according to the measurements made in all populations were petiole length (1.16 cm), leaf width (5.11 cm), leaf length (7.96 cm), length of lamina (6.79 cm), leaf area (26.62 cm²), leaf vein angle (54.21°), and number of primary leaf veins (20.03).
Table 3. Correlations Among Leaf Variables and Altitude
* Significance level (P) < 0.05 statistically difference (2-tailed)
Fig. 5. Dendrogram obtained by hierarchical cluster analysis including all leaf characteristics. Group one is Artvin-Arhavi, group two is Artvin-Murgul, and group three are Artvin-Camili and Rize-Çamlıhemşin
According to the ANOVA test, no significant correlation between annual temperature (while the annual temperature average value is 13.1 °C in Artvin-Arhavi, Borçka (Camili) and Murgul regions, it is 15.5 °C in Rize Çamlıhemşin region), annual precipitation (while annual precipitation total value is 685.5 mm in Artvin-Arhavi, Borçka (Camili) and Murgul regions, it is 2307.3 mm in Rize Çamlıhemşin region) and leaf traits was noted in this study; this could be because the leaf samples came from the same area (the Eastern Black Sea region). Regarding the variance analysis for all populations, there were meaningful results based on leaf length, length of lamina, leaf area, and altitude. The results of the correlation analysis for the measured seven leaf characteristics and altitude are provided in Table 3. According to the results of the correlation analysis, there were statistically significant positive and negative correlations between leaf characteristics, while the correlations for PL, LL, LW, LA, LVA, and NLV were not statistically significant. Both LLA and LL had a positive correlation, and LA and altitude had a negative correlation. Hierarchical clustering analysis was performed to determine how the populations were graphically grouped according to the measured leaf parameters (Fig. 5).
The significance of the groupings formed because of the hierarchical cluster analysis was tested via separation analysis, resulting in two groups. According to all measured leaf characteristics, the population from Artvin-Arhavi was in one group, while the populations from Artvin-Murgul, Rize-Çamlıhemşin, and Artvin-Borcka (Camili) formed another group.
DISCUSSION
The present study showed that morphological traits and geoclimatic parameters (combination of geographical and climatic conditions) were closely related to the distribution of Betula medwediewii populations, as this species occurred in a restricted area of Türkiye (Artvin and Rize, only in the Eastern Black Sea Region).
When the regions of Artvin-Arhavi (1772 to 2065 asl), Artvin-Borcka, Camili (1753 to 1792 asl), Artvin-Murgul (1472 to 1498 asl), and Rize-Çamlıhemşin (1818 to 1827 asl) were evaluated in terms of the natural distribution area of Betula medwediewii, this species was found at different elevations. However, this situation did not result in statistically significant differences in terms of leaf characteristics. The population from Artvin-Borcka, in the Camili region, provided the best results in terms of many morphological features (Table 2), most likely because of the other ecological features (microclimate, soil) in this region.
Regarding leaf morphological variation in relation to habitat conditions, Betula L. often shows significant differences in leaf morphological characteristics, such as leaf size and shape (Dancik et al. 1974; Sharik and Barnes 1979; Joel et al. 1994; Aspelmeier and Leuschner 2006; Uribe-Salas et al. 2008; Pyakurel and Wang 2013). Here, noticeable differences were also found in morphological leaf traits for B. medwediewii. The results of present study are in line with the findings of Kundu and Tigerstedt (1997), Bruschi et al. (2003), and Calagari et al. (2006) for other species (Azadirachta indica A. Juss., Quercus petraea (Matt.) Liebl., Quercus pubescens Willd., and Populus euphratica Oliv., respectively). Kleinschmit et al. (1995) and Bruschi et al. (2000) measured several morphological leaf parameters in Quercus petraea (Matt.) Liebl and found differences among the populations. Hiura et al. (1996) found differences in leaf size among various provenances of Fagus crenata Blume.
Information on leaf morphological traits, widely available in flora, could be used to strengthen predictive models of species distribution and vegetation function. In addition, leaf morphological trait variations could be exploited in models to assess the consequences of climate change for species and ecosystems (Wang et al. 2022). Through examining leaf morphological traits, it will be possible to obtain information about the climate, flora, and vegetation, which make up the ecosystem, using less time and less cost. Populations may persist in their current areas and withstand environmental changes if they have adaptive capacity (Gratani et al. 2012). However, this adaptation capacity of B. medwediewii seems quite low, considering its limited distribution areas. For that reason, B. medwediewii occasionally grows in arboreta and a few large gardens in Europe. Further research, including long-term monitoring, phenology studies, and seed germination capability studies, is needed to establish the competitive capacity of B. medwediewii. This plant taxon was assessed as category DD according to the IUCN criteria (Wilson et al. 2018). In this context, determining the efficiency of in situ conservation and development of propagation techniques of this taxon is very important.
The essential oil of leaves of B. medwediewii has antifungal properties against plant pathogenic fungi (Demirci et al. 2000). Taking into this account, the cultivation of this species can generate economic income. In this context, it is important to know the natural distribution areas of this species. To determine individuals that will outperform in the productive generation, leaf characteristics of individuals showing optimum and maximum development need to be examined. In addition, it is recommended that B. medwediewii individuals, which are thought to show superior performance in generative progeny, should be obtained from Artvin-Borçka and Artvin-Murgul regions, which have the healthiest and optimum leaf characteristics (according to the highest and optimum measurement values) (Table 2). Individuals from Rize-Çamlıhemşin can also be selected, but individuals in the Artvin-Arhavi region have very low values in terms of both optimum and maximum value of leaf characteristics. In this sense, determining the efficiency of in situ conservation and production of this taxon is important.
Betula medwediewii is categorized as vulnerable (VU) under criterion B1ab (iii,v) in the “Red List of the Endemic Plants of the Caucasus”, based on a provisional assessment undertaken during a Global Tree Specialist Group workshop in Tbilisi in 2005. Betula medwediewii nearly meets the threshold for qualifying as Vulnerable under criterion B1ab (Eristavi et al. 2014; Shaw et al. 2014). This taxon is heavily overgrazed by wild and domestic animals, especially in the Artvin-Borçka (Camili) region. Grazing causes damage to both the trees and the seeds, thereby reducing natural reproduction. In consequence, this taxon must be protected in situ (closing some areas to prevent grazing), or ex situ (greenhouse propagation of seedlings for reforestation strategies).
Betula medwediewii, which is a kind of relict plant taxon, has a narrow natural distribution in Türkiye. Although a very extensive field study was conducted in the Artvin-Hatila valley, which is one of the places specified as its natural distribution area in the flora of Turkey (Browicz 1972; Davis 1988), this plant taxon has not been encountered between 2019 and 2020. In addition, in other natural distribution areas of this plant taxon, a limited number of individuals were encountered. While plant individuals were determined, the distribution of plant taxa were too small in the research areas, and this plant taxon did not form a stand. The raised decline natural distribution of B. medwediewii, is concerned that this plant taxon would disappear in the natural distribution areas, in the future.
This manuscript analyzed leaf morphological characteristics of four populations, which have only restricted distribution in Türkiye, growing at different elevations and climatic conditions. Through researching the leaf characteristics of a particular location, ecological information about this location can be obtained. Thus, leaf morphology studies can facilitate the development of more economic and time-efficient cultivation practices. Depending on the plant taxa, the obtained data can be used to select places for plant production, determining the appropriate ecological conditions (e.g., climate, elevation, and aspect groups) and thereby reducing the mistake rate.
CONCLUSIONS
Determining variation within and among populations of Betula medwediewii was completed according to morphological leaf traits.
- The population from Artvin-Borcka, in the Camili region, provided the best results in terms of many morphological features. The small habitat range of this endemic and rare species often leads to a high risk of extinction (Rodrigues et al. 2006; Mhamdi et al. 2013; Yildirim and Turna 2021). For that reason, biodiversity conservation and the preservation of variation are important for the conservation of genetic resources. In situ and ex situ conservation activities should be applied to this plant taxon.
- Significant relationships between population size or rarity and measures of genetic variation have been found in several rare plant species (Karron et al. 1988; Billington 1991; Treuren et al. 1991; Oostermeijer et al. 1994). In this sense, the possibility of selecting individuals with superior performance in generative progeny is an important issue.
- In terms of leaf morphology, the individuals who have had the healthiest and optimum leaf characteristics, will likely be individuals that will outperform in the productive generation.
- As a result of the present study and many studies (Joel et al. 1994; Bruschi et al. 2003; Calagari et al. 2006; Aspelmeier and Leuschner 2006; Uribe-Salas et al. 2008; Pyakurel and Wang 2013; Yildirim and Turna 2021), it has been seen that elevation and climate types (especially microclimatic characteristics) affect the morphological characteristics of the leaves, so this should be taken into account especially when the using activities of leaf morphology as a bioindicator against environmental factors are carried out.
ACKNOWLEDGEMENTS
The author would like to express sincerest thanks to the staff of the Forest Enterprises of Artvin (Borçka and Murgul) and Rize (Çamlıhemşin) for their kind assistance with field studies. Dr. Nebahat Çimen was acknowledged for assistance with field studies.
REFERENCES CITED
Akman, Y. (1999). İklim ve Biyoiklim [Climate and Bioclimate], Kariyer Matbacılık, Ankara, Turkey.
Akman, Y. (2011). İklim ve Biyoiklim (Biyoiklim Metodları ve Türkiye İklimleri) [Climate and Bioclimate (Bioclimate Methods and Climates of Turkey)], Palme Yayın-Dağıtım, Ankara, Turkey.
Aspelmeier, S., and Leuschner, C. (2006). “Genotypic variation in drought response of silver birch (Betula pendula Roth: Leaf and root morphology and carbon partitioning,” Trees 20, 42-52. DOI: 10.1007/s00468-005-0011-9
Batan, N., and Özdemir, T. (2013). “Contributions to the moss flora of the Caucasian part (Artvin Province) of Turkey,” Turk. J. Bot. 37(2), 375-388. DOI: 10.3906/bot-1201-49
Billington, H. L. (1991). “Effect of population size on genetic variation in a dioecious conifer,” Conserv. Biol. 5, 115-119.
Browicz, K. (1972). Flora of Turkey and the East Aegean Islands, P. H. Davis (Ed.), Vol. 7, University Press, Edinburgh, UK, pp. 688-691.
Bruschi, P., Vendramin, G. G., Bussotti, F., and Grossoni, P. (2000). “Morphological and molecular differentiation between Quercus petraea (Matt.) Liebl. and Quercus pubescens Willd. (Fagaceae) in Northern and Central Italy,” Ann. Bot-London 85, 325-333. DOI: 10.1006/anbo.1999.1046
Bruschi, P., Grossoni, P., and Bussotti, F. (2003). “Within and among tree variation in leaf morphology of Quercus petraea (Matt.) Liebl. Natural populations,” Trees 17, 164-172. DOI: 10.1007/s00468-002-0218-y.
Buser, S., and Cvetic, S. (1972). Türkiyede’ki Murgul Bakır Yatağı Yöresinin Jeolojisi, Ljubljana Jeoloji Enstitüsü ve Nükleer ve Diğer Mineral Hammaddeleri Etüt ve Aramaları [Geology of the Murgul Copper Deposit Region in Turkey, Ljubljana Geological Institute and Survey and Exploration of Nuclear and Other Mineral Raw Materials], Jeoloji ve MadenEnstitüsü, Yayınları, Turkey.
Calagari, M., Modirrahmati, A. R., and Asadi, F. (2006). “Morphological variation in leaf traits of Populus euphratica Oliv. natural populations,” Int. J. Agri. Biol. 8(6), 754-758. DOI: 1560–8530/2006/08–6–754–758
Coleman, J. S., McConnaughay, K. D. M., and Ackerly, D. D. (1994). “Interpreting phenotypic variation in plants,” Trends. Ecol. Evol. (Personal Edition) 9, 187-191.
Dancik, B. P., and Barnes, B. V. (1974). “Leaf diversity in yellow birch (Betula alleghaniensis),” Can. J. Botany 52, 2407-2414. DOI: 10.1139/b74-312.
Davis, P. H. (1988). Flora of Turkey and the East Aegean Islands, Vol. 7, University Press, Edinburgh, UK.
Demirci, F., Demirci, B., Baser, K. H. C., and Guven, K. (2000). “The composition and antifungal bioassay of the essential oils of different Betula species growing in Türkiye,” Chem.