ISSUE 03/2024 

Analyzing Population Diversity of Common Hornbeam (Carpinus betulus L.) in Native Habitats in Türkiye: A Focus on Leaf Characteristics

Article by Fahrettin Atar, Ebru Atar, Deniz Güney, İbrahim Turna | 01.10.2024 - 11:43

Schlüsselbegriffe: Plastizität, Anpassung, Vielfalt, Höhenlage, Gemeine Hainbuche

Available at https://doi.org/10.53203/fs.2403.1

Download as PDF

See below the issue 3/2024 as E-Paper or have a look at our E-Paper archive dating back to 1955.

Abstract

The primary aim of this investigation was to examine the morphological diversity of leaves in naturally occurring populations of common hornbeam (Carpinus betulus L.) in Türkiye. Twelve common hornbeam populations, distributed over four distinct watersheds in the Eastern Black Sea Region, were sampled from three altitude levels, ranging from sea level to 1200 meters. Nine morphological characteristics of the leaves were analyzed biometrically. We found a high degree of phenotypic heterogeneity both within and among the investigated populations. The eco-geographic principle's classification of populations was disclosed. Altitude appeared to be the primary factor influencing variation in the majority of leaf attributes from the Camlıhemşin, Çaykara, and Espiye watersheds. Compared to trees at lower altitude, trees at higher altitudes were characterized by smaller leaves. The populations in the Trabzon-Maçka watershed, however, did not show this trend. Largest leaves were observed for trees from the Trabzon-Maçka watershed region at the second altitude level. Overall, our findings support the notion that the Eastern Black Sea area is a hotspot of biodiversity and that the intricate evolutionary process leading to leaf morphological variability may be a consequence of adaptation and plasticity.

Zusammenfassung

Das Ziel dieser Arbeit ist, die morphologische Vielfalt von Blättern in natürlich vorkommenden Populationen der Gemeinen Hainbuche (Carpinus betulus L.) in der Türkei zu untersuchen. Zwölf Populationen der Gemeinen Hainbuche, die sich über vier verschiedene Einzugsgebiete in der östlichen Schwarzmeerregion verteilten, wurden innerhalb drei Höhenstufen zwischen Meereshöhe und 1200 Meter Seehöhe beprobt. Neun morphologische Merkmale der Blätter wurden biometrisch analysiert. Es wurde ein hoher Grad der phänotypischen Heterogenität sowohl innerhalb als auch zwischen den untersuchten Populationen festgestellt. Darüber hinaus wurde die Klassifizierung der Populationen nach dem ökogeografischen Prinzip aufgedeckt. Die Höhenlage war der primäre Faktor, der die Variation der meisten untersuchten Blattmerkmale in den Einzugsgebieten Camlıhemşin, Çaykara und Espiye signifikant beeinflusste. Im Vergleich zu Bäumen auf niedrigerer Seehöhe zeichneten sich Bäume in höheren Lagen durch kleinere Blätter aus. Die Populationen im Einzugsgebiet Trabzon-­Maçka zeigten diesen Trend nicht. Die größten Blätter in dieser Region wurden in der zweiten Höhenstufe beobachtet. Insgesamt weisen unsere Ergebnisse darauf hin, dass die östliche Schwarzmeerregion ein Hotspot für Biodiversität ist und der komplexe evolutionäre Prozess, der zur morphologischen Variabilität der Blätter führt, eine Folge von Anpassung und Plastizität ist.

1 Introduction

Plant adaptability and vitality features (reproduction, seeding, and growth) are affected by both biotic and abiotic influences in varied ecosystems (Matesanz et al. 2010; Varol et al. 2022). Changes in ecological variables, whether spatial or temporal, can impact the morphological features of plants (Herrera and Bazaga 2013; Koç et al. 2024; Özel et al. 2024). The development and adaptation of plant anatomical, physiological, and morphological characteristics are significantly influenced by the limitations associated with these variables. These limitations play a crucial role in shaping the diverse traits observed among plant populations (Nascimbene and Marini 2015). The eco-geographic principle's classification refers to the categorization of populations based on the geographic and ecological characteristics of their habitats. This principle considers factors such as altitude, latitude, climate, and soil types, which influence the phenotypic traits and genetic diversity within and among plant populations. Understanding these relationships helps in elucidating the adaptive mechanisms of species to different environmental conditions (Tekin et al. 2022).

The soil, temperature, and precipitation patterns of an ecosystem are all significantly influenced by altitude (Körner 2007; Dogan et al. 2023). The life history properties of plants and their evolutionary reactions could be influenced by a variety of climatic factors, encompassing variables such as rainfall, temperature and light (Leingärtner et al. 2014; Sevik et al. 2017; Ertürk et al. 2024). Species dwelling at higher altitudes exhibit morphological and physiological distinctions from those at lower altitudes, leading to differentiated characteristics based on the elevation changes (Körner 1999). Some researchers have reported that the morphological and anatomical characteristics of plants are influenced by altitude (Cordell et al. 1998; Özbucak et al. 2013). This situation indicates the adaptation ability of plants to factors such as water and nutrient deficiencies, temperature changes, etc. (Vitousek 1982; Aragon et al. 2012). 

The most significant organs of plants that vary in size, color, and shape are their leaves (Li et al. 2019). The capacity of a leaf to modify its morphology and physiology in response to environmental changes is known as phenotypic plasticity, and it is caused by contrasting environments (Schlichting 1986; Sultan 2000). Zotz et al. (2011) state that many plants undergo sudden and significant modifications in leaf morphology as a response to external influences. These plants are able to adapt to a changing environment due to the plasticity of their leaf morphology (Wells and Pigliucci 2000). Elevation has a significant impact on leaf morphology and physiology within a species (Körner 1999). With an increase in leaf surface area, plants exhibit a proportional increase in photosynthesis, transpiration, and respiration (Özen and Onay 2013). Therefore, it can be stated that the morphological characteristics of leaves are associated with the growth of plants. Moreover, specific leaf area (SLA), identified as an informative trait (Wright et al. 2004; Poorter et al. 2009), serves as an indicator for ecophysiological traits such as stress tolerance, leaf longevity, and growth rate (Weiher et al. 1999; Wright and Westoby 2002). A higher SLA signifies a reduction in leaf thickness, diminished leaf density, or a combination of the two (Scheepens et al. 2010). Poorter et al. (2009) highlighted the strong correlation between SLA and factors such as water availability, temperature and irradiance.

Leaf morphological characters such as length, width, and area generally exhibit a negative correlation with an increase in altitude (Körner et al. 1986). Paridari et al. (2013) observed that trees at higher elevations tend to have smaller leaf laminas compared to those at lower altitudes. In tropical highlands, as altitude increases, leaf mass per area (LMA) increases, while leaf size often decreases (Grubb 1977; Korner 1989). Therefore, studies investigating the relationships between altitude and plant morphology provide valuable insights into ecological research. (Paridari et al. 2013). Various outcomes have emerged in studies on plant morphometrics (Poljak et al. 2018; Adamidis et al. 2021). Morphological variability may be significantly impacted by variations in a number of variables, including temperature, light, precipitation, slope, depth of soil, and relative humidity. Forest trees commonly display adaptations to various environmental conditions (Savolainen et al. 2007). Additionally, genotype, environment, and their interactions all influence the phenotype of a tree (Falconer and Mackay 1996; Cobanoglu et al. 2023; Yigit et al. 2023). Moreover, morphological traits within a single tree can differ greatly, particularly in big trees (Hagemeier and Leuschner 2019). Additionally, morphological variability can be influenced by the kind of forest community, soil acidity and habitat, management type (Poljak et al. 2022), and other factors (Adamidis et al. 2021).

Including different taxon delimitations (Mac-Key 1988; Sękiewicz et al. 2016), population variability (Douaihy et al. 2012; Poljak et al. 2018),  morphological seed involucre variation (Boratyński et al. 2007; Atar 2022), physiological and morphological seed characterization (Drvodelić et al. 2015; Atar et al. 2017; Atar et al. 2020), cultivar characterization and selection (Ertan 2007; Poljak et al. 2016), and variation of macro- and micromorphological leaves (Bednorz 2006; Paridari et al. 2013; Amer et al. 2016; Güney et al. 2016; Aykut et al. 2017), morphological knowledge remains crucial in numerous domains within the field of plant science (Poljak et al. 2015; Atar and Turna 2018). Additionally, morphological characterization remains the formal process for cultivar registration and protection (Pereira-Lorenzo et al. 1996). Furthermore, the evolutionary adaptability of plants to various environments may be predicted using morphological plasticity (Nicotra et al. 2010).

The Carpinus genus, belonging to the Betulaceae family, comprises approximately 40 species indigenous to the Northern Hemisphere, ranging from Europe to eastern Asia, extending south to the Himalayas, and found in North and Central America (Krüssmann 1984; Furlow 1990; Hillier 1991; Suszka et al. 1996). In Türkiye, the oriental hornbeam (Carpinus orientalis Mill.) and common hornbeam (Carpinus betulus L.) species are naturally distributed. The common hornbeam is found in the Black Sea Region, which constitutes the entire Europe-Siberia floral area in Türkiye. Additionally, it locally spreads in the Amanos Mountains. The general geographical distribution of the common hornbeam includes all of Europe, Türkiye, the Caucasus, and Iran. It typically forms mixed forests with beech, chestnut, oak, and alder and occasionally occurs in pure stands (Yaltırık 1982; Anşin and Özkan 2006). The common hornbeam extends in the Black Sea forests from sea level up to an elevation of 1200–1300 m.

The aims of this study were to: 

(1) assess how the morphological variability of some common hornbeam leaf traits varies with altitude; 

(2) identify leaf variations within and between populations of the species; and 

(3) reveal how the morphological plasticity of the leaf traits varies with eco-geographic factors.

2 Materials and Methods

Atar_Fig1.jpg

Figure 1: The geographical locations of common hornbeam populations (Caudullo et al. 2017).

Abbildung 1: Die geografischen Standorte von Gemeine-Hainbuche-Populationen (Caudullo et al. 2017).

2.1 Selection of populations and the collection of leaf

Natural populations of common hornbeam were sampled for this study from four distinct watersheds in the Eastern Black Sea Region's forests. Leaf samples of common hornbeam were collected from populations that grew at varying altitudes in the watersheds of Rize-Çamlıhemşin, Trabzon-Çaykara, Trabzon-Maçka, and Giresun-Espiye in Türkiye (Figure 1).

The research encompassed 12 populations distributed across three distinct altitude zones, reaching up to 1200 meters above sea level (Table 1). Eco-geographic principle's classification was applied by analyzing the morphological traits of the leaves in relation to the altitude and geographic location of the sampled populations. This involved categorizing populations into distinct eco-geographic zones based on their altitude ranges (0-400 m, 400-800 m, 800-1200 m) and comparing the morphological variability within and between these zones. A total of 134 trees were sampled for leaves: 30 from Maçka, 46 from Çaykara, 30 from Çamlıhemşin, and 31 from Espiye. Leaf samples were collected from the outer (light-exposed) section of the tree crowns by climbing each tree, resulting in a total collection of 6850 leaves (50 leaves × 137 trees). Within the scope of the study, the heights of trees selected from each population vary between 15-20 m and their diameters range from 18-25 cm. The populations were selected from pure stands of hornbeam and mixed stands where hornbeam predominates, with beech, oak, and chestnut also present in the mix.

Atar_Tab1.jpg

Table 1: Information on the sites, where populations of common hornbeam were sampled.

Tabelle 1: Informationen der Standorte, auf denen Populationen der Gemeinen Hainbuche untersucht wurden.

2.2 Measurements of Leaf Characteristics

In the study, morphological characteristics such as leaf length, leaf width, leaf area, specific leaf area, petiole length, number of leaf veins, leaf vein angle, leaf shape index, and leaf length to petiole length ratio were measured (Table 2). Leaf vein angle has been obtained by measuring 3 different parts of the leaf, including the bottom, middle, and upper parts of the leaf. The ratio of leaf length to the corresponding leaf width was used to calculate the leaf shape index (Jeong et al. 2011; Wang and Zhang 2012). 

Determination of the specific leaf area: first of all, the area of each leaf was measured. Then the leaves were kept in a drying kiln for 24 hours at 65˚C. After that, oven-dry weights were measured by precision scales. After these processes, specific leaf area values were determined by proportioning the dry weight of the leaf area (Kardel et al. 2010; Scheepens et al. 2010).

Leaf samples were scanned at a resolution of 1200 dpi using a scanner (Hewlett-Packard Scanjet G4010). Subsequently, leaf-related measurements were conducted utilizing the Image Analysis Software (ImageJ). This software has been employed in various scientific studies to analyze variations in leaf characteristics (Bayramzadeh et al. 2008; Atar 2022).

 

Atar_Tab2.jpg

Table 2: List of the leaf morphological characteristics studied.

Tabelle 2: Liste der untersuchten blattmorphologischen Merkmale.

2.3 Statistical analysis

Descriptive statistical parameters were computed for each population to assess the extent of variation in each characteristic. Analysis of variance (ANOVA) was employed to identify whether there were statistically significant differences in terms of measured morphological characters both among populations and within populations. Using principal component analysis, hierarchical cluster analysis, K-means clustering method and discriminant analysis, the overall population diversity of leaf morphological characteristics was investigated. Pearson’s correlation coefficient was utilized to assess the correlation between pairs of morphological characteristics and altitude. The morphological plasticity index (PI) was calculated following the formula from Ceplick (1995): PI = (Xmax - Xmin) / Xmax; where Xmin and Xmax are the minimum and maximum mean value measurement of each morphological characteristics. The morphological plasticity index takes values between 0 and 1. The statistical package programs "R v.3.4.3" and "Windows SPSS Software 23.0" were used to analyze the data.

3 Results

The summary of the results from the descriptive statistical analysis (including arithmetic means and coefficients of variation), analysis of variance, and Duncan tests for the leaf traits in the sampled populations is presented in Table 3. It was determined that the values of leaf length, leaf width, and leaf area decrease with an increase in elevation in populations from each watershed. Generally, the highest values were detected at lower altitudes. However, in the Maçka watershed, an exception was observed where the leaf width did not show a significant decrease with increasing altitude, and in some cases, even showed slightly higher values at mid-altitudes (400-800 m) compared to the lowest altitude (0-400 m). In contrast, SLA and petiole length values increased with altitude in each population. The number of leaf veins and vein angle, however, did not exhibit a systematic change associated with elevation. For the morphological characteristics under study, the overall coefficients of variation varied from 11.1% (V) to 34.7% (LL/PL). Higher altitude populations showed the least amount of morphological variation for the bulk of the leaf characteristics that were studied. 

Atar_Tab3.jpg

Table 3: The results of analysis of variance, Duncan test and descriptive statistics concerning the values of leaf characteristics among populations. For abbreviations please see Table 1 and 2.

Tabelle 3: Ergebnisse der Varianzanalyse, des Duncan-Tests und der deskriptiven Statistiken bezüglich der Werte von Blattmerkmalen zwischen den Populationen. Für Abkürzungen siehe Tabelle 1 und 2.

The analysis of variance revealed significant differences (P<0.01) both within and between populations across all measured morphological characters. Groups formed by populations for each morphological character were identified using the Duncan test. Accordingly, LL, V, and LA resulted in 10, 6, and 5 groups, respectively. Furthermore, there were 8 groups among populations for SLA and PL and 7 groups for LW, LVA, LSI, and LL/PL (Table 3). 

The averages of populations located at the same altitude zones were taken to determine the variations in the measured characters across three altitude zones. The mean values, standard deviations, minimum and maximum values of first altitude (mean of CH1, C1, M1, E1 populations), second altitude (mean of CH2, C2, M2, E2 populations) and third altitude (mean of CH3, C3, M3, E3 populations) are presented in Figure 2. It has been observed that LL, LW, LA, LSI, and LL/PL values decrease with increasing altitude. On the contrary, SLA and PL values increased with altitude. No significant variation related to elevation was noted in the V and LVA values. The level of significance for all morphological characters measured across different altitudes was determined to be P<0.001. Additionally, based on the Duncan test results, each altitude formed a separate group for LL, SLA, LVA, LSI, and LL/PL characters. For LW, LA, PL, and V characters, two distinct groups were identified among altitudes.

Atar_Fig2.jpg

Figure 2: Variation of morphological characteristics of common hornbeam leaves depending on altitudes.

Abbildung 2: Variation der morphologischen Merkmale von Blättern der gemeinen Hainbuche in Abhängigkeit von der Höhenlage.

The results revealed significant variations in the plasticity index (PI) of each characteristic (Figure 3). Overall, the plasticity index in the 1st altitude zone was generally low when considering all parameters. LA, PL, SLA, and LL/PL parameters exhibited notably high plasticity indices. The plasticity index of LA was maximal at all three altitudes. In the 1st and 3rd altitudes, V, and in the 2nd altitude, LSI had the minimum plasticity index. In Figure 3, relationships between each morphological characteristic were depicted as a correlation matrix. Accordingly, it was determined that many measured leaf characters have significant correlations with each other at a 95% confidence level.

Atar_Fig3.jpg

Figure 3: A, Plasticity index depending on altitudes and B, Pearson correlation coefficients between pairs of morphological characters.

Abbildung 3: Plastizitätsindex in Abhängigkeit von den Höhenlagen und Pearson-Korrelationskoeffizienten zwischen Paaren von morphologischen Merkmalen.

Five principal components, elucidating approximately 96.71% of the overall variability, were derived as outlined in Table 4. PC I accounted for only one-third of the total variance (35.40%), followed by the subsequent five PCs (20.06%, 16.18%, 12.55%, and 12.52%, respectively). The majority of leaf-size parameters, including LL, LW, LA, and LL/PL, exhibited the strongest correlations with PC I, aligning with expectations. Conversely, PL demonstrated a negative correlation with PC II, while the only positively correlated variable with PC II was LSI. Additionally, LVA exhibited a high correlation with PC III, V was significantly correlated with PC IV, and the sole variable positively correlated with PC V was SLA.

Atar_Tab4.jpg

Table 4: Results of the PCA on morphological characteristics. 

Tabelle 4: Ergebnisse der Hauptkomponentenanalyse (PCA) der morphologischen Merkmale.

K-means clustering and hierarchical cluster analysis were used to determine the degree of dissimilarity or similarity among the investigated populations (Figure 4). Discriminant analysis was used to evaluate the significance of the groups that emerged from the cluster analysis. Discriminant analysis revealed that splitting into two groups was statistically significant (P<0.05). Accordingly, in terms of all measured leaf characters, the Çaykara-3, Maçka-1, Espiye-2, and Espiye-3 populations were in the same group, and the other populations formed the other group. Although statistically significant separation into two groups was identified, cluster analysis suggests the potential existence of six distinct groups based on leaf characters. The spatial distribution of these six different groups is illustrated in Figure 4. Populations C3 and E2 are grouped together, while M1 and E3 form another distinct group. The remaining populations do not belong to any specific group. As depicted in Figure 4, populations from neighboring basins tend to form clusters with each other.

Atar_Fig4.jpg

Figure 4: A, Dendrogram and B, map of population groups related to hierarchical cluster analysis.

Abbildung 4: Dendrogramm und Karte der Populationengruppen im Zusammenhang mit der hierarchischen Clusteranalyse.

4 Discussion

High diversity in morphological characteristics was demonstrated by the research done, and variations within and between populations were confirmed for all the parameters that were evaluated. Our findings showed the non-homogeneity of populations concerning leaf characteristics. Variation in most of these characteristics appeared to be related mostly to altitude. Moreover, all measured leaf characteristics provided evidence that the species has evolved plastic strategies based on varying altitudes. Additionally, variations and differences among populations at the same altitude may be due to different gene and genotype frequencies among and within populations. Di Pierro et al. (2017) investigated the genetic basis of local adaptation across altitudinal gradients in 18 natural Alpine Norway spruce populations and reported that most of the genetic variability was found within populations, but there was little but significant variation between geographic groups.

The leaf morphological characteristics of the majority of woody plants are influenced by a variety of abiotic variables along environmental gradients. Several studies have reported substantial phenotypic variations as adaptive responses to changes in altitude (Premoli et al. 2007; Bresson et al. 2011; Varol et al. 2022). In our study, it was observed that the leaf length, leaf width, and leaf area values in populations from each watershed decrease with increasing altitude, and generally, the highest values are obtained in populations located at lower altitudes. Atar et al. (2014) reported significant differences in the morphological characteristics of Carpinus orientalis leaves across three different altitudes. They indicated that leaf width ranged from 2.11 to 2.20 cm, leaf length from 4.24 to 4.38 cm, leaf area from 6.39 to 6.66 cm2, and specific leaf area from 157.56 to 170.15 cm2 g-1. Paridari et al. (2013) obtained similar findings in their morphometric study on common hornbeam leaves. Their research revealed that trees at higher altitudes exhibited smaller leaf lamina compared to those at lower altitudes. In the study conducted by Atar (2022), 13 morphological characteristics were analyzed in common hornbeam seed involucres, and it was reported that populations exhibit high phenotypic variability both among themselves and within their own groups. Additionally, trees at higher altitudes were characterized by smaller involucres compared to those at lower altitudes. Moreover, a study investigating how altitude affects the biometric characteristics of hornbeam tree indicated important effects of elevation variations on the characteristics under investigation (Kiaei et al. 2019). Changes in abiotic factors such as soil moisture, air temperature, and atmospheric CO2 concentrations along altitude gradients likely lead to strong differential selection, resulting in the local adaptation of Carpinus betulus populations (Panek and Waring 1995; Beerling et al. 1996; Halbritter et al. 2018).  

It has been reported for a variety of plant species that as altitude increases, leaf length, width, and area typically decrease. The morphological traits and morphological plasticity changes of 17 populations of Primula nivalis at varying elevations were investigated by Abdulsalam and Li (2018). As a result of their study, it was reported that plant height, leaf length and petiole indicated a significant negative correlation with altitude and as the altitude increased, plant height and leaf length decreased. In the study investigating the morphological characters of the mericarp of 24 species in the Selineae tribe, it was stated that the length of mericarp varies from 2.7 to 11.4 mm, while the width of mericarp varies from 1.4 to 9.3 mm (Lee et al. 2018). Yousefezadeh et al. (2010) observed significant variations in parameters such as leaf width, leaf vein angle, and number of leaf pairs in Parrotia persica based on altitude. Moreover, Korner et al. (1986) reported that the length, width, and area of leaves often decrease with an increase in altitude. Research on Alnus incana conducted by Poljak et al. (2018) indicated an inverse relationship between petiole length and leaf blade area values concerning altitude.

A high level of specific leaf area (SLA) indicates a wider ecological niche for a given species and rising production (photosynthesis) (Ishii et al. 2018). In our study, SLA and PL values increased in each population with the altitude's rise. Another study conducted in Hawaiin, examining the variations in leaf morphology of Metrosideros polymorpha species based on environmental factors, also reported similar results. It was observed that as the elevation increases from 70 m to 2350 m, leaf length per area decreases and leaf mass per area increases (Joel et al. 1994). Additionally, the increase in leaf mass per area with elevation that it was observed could reflect decreasing nutrient availability with increasing elevation (Vitousek et al. 1992). Öztürk et al. (2015) analyzed leaf area and leaf area index (LAI) parameters as indicators of leaf development to determine the impact of air and soil temperature on Carpinus betulus leaves. The study reported strong positive correlations between air and soil temperatures and leaf area as well as leaf area index. In an extensive investigation of Catalan Quercus ilex forests by Ogaya and Peñuelas (2007), increased leaf mass per area was noted in drier conditions, as well as in sunnier locations with lower air temperatures. Similarly, Kolozsvari et al. (2020) evaluated certain morphometric and physiological indices on two woody plant species, Fagus sylvatica and Carpinus betulus, and characterized their plant behavior in different forest regions, elevations, and light intensities. The study concluded that hornbeam leaves exhibited higher specific leaf weight compared to beech leaves, with the highest value obtained from samples taken from trees growing at forest edges where light intensity was higher than within the forest.   Stiegel and Mantilla-Contreras (2018) conducted an experimental study aiming to analyze the environmental effects, including microclimate and soil moisture, prevalent in the upper canopies of forest stands, on leaf traits of juvenile Fagus sylvatica, Acer pseudoplatanus, and Carpinus betulus. Following the treatments, it was observed that the leaves of Carpinus betulus exhibited a tendency towards decreased SLA, attributed to the treatment involving warmer and drier microclimates. The average SLA values recorded for Carpinus betulus were 155 cm2 g-1. Furthermore, the environmental impact on SLA was found to be more pronounced under conditions of lower soil moisture. Indeed, it has been highlighted that leaves possessing elevated leaf mass per area exhibit a structure that reduces both transpirational water loss and photosynthetic rates (Hultine and Marshall 2000). In another study investigating the variation in SLA of mature Carpinus betulus trees in a closed forest, based on leaf position and growth period, it was reported that the SLA values were lowest in the upper part of the crowns, with the highest SLA values obtained in measurements taken in May (Huzulák and Eliáš 1975). 

The results revealed significant differences in the plasticity index for each characteristic. LA, PL, SLA, and LL/PL parameters exhibited notably high plasticity indices. In research conducted by Paridari et al. (2013) on Carpinus betulus species, characteristics such as vein count, tooth count, and stomatal sizes demonstrated minimal plasticity, whereas stomatal density and leaf area exhibited the greatest degree of plasticity. The highest plasticity was found in the leaf-apex length of Alnus subcordata by Akbarian et al. (2011) and the leaf-petiole length of Castanea sativa by Zarafshar et al. (2010).

In this study, although the diversity of the examined leaf characteristics in populations investigated in the Çamlıhemşin, Çaykara, and Espiye basins was mostly associated with altitude, this was not the case for populations in the Trabzon-Maçka watershed. In this region, trees at the second altitude exhibited the largest leaf sizes. Atar et al. (2014) conducted a study investigating the elevation-dependent variation of certain leaf characteristics in Carpinus orientalis in the Trabzon-Maçka basin, which encompasses three elevation zones, similar to our study. Consistent with our study results, this research also found that parameters such as leaf length, leaf area, and specific leaf area exhibited the highest values in populations belonging to the second altitude zone. This situation is likely dependent on habitat variables related to management type and soil characteristics. Common hornbeam forests are managed in two ways: as high forests and as coppice forests. Hornbeam forests are managed as high forests in the majority of the populations under investigation, and they are usually found in mixed stands with chestnut, beech, and oak. However, the hornbeam trees located in the first altitude zone of the Maçka watershed may be continually cut and regenerated, resembling a coppice management system. Therefore, the lower leaf sizes of trees at the first altitude can be explained by this phenomenon. Poljak et al. (2022) reported similar findings, indicating that chestnut stands managed as coppice yield smaller fruits compared to chestnut stands managed as high forests. Furthermore, they highlighted that factors such as soil acidity, habitat type, and forest community type can influence fruit size (Poljak et al. 2012). The washing caused by excessive rainfall on the slopes facing north in the Eastern Black Sea region leads to lime removal from the soil. Consequently, the soil exhibits an acidic nature. In addition to the kind of management, habitat modifications, like soil acidity, can significantly affect leaf size in the Maçka watershed.

5 Conclusions

In conclusion, genetic diversity plays a pivotal role in forest ecosystem adaptation and stability overall, especially in the face of environmental stress and inadequate management that pose long-term threats. Therefore, it is essential to delineate the genetic makeup of C. betulus and other natural forest resources, particularly to facilitate the implementation of both in-situ and ex-situ conservation efforts. This research highlights the significance of altitude, along with its associated factors like rainfall and temperature, as a key driver of morphological variations in C. betulus leaves. Additionally, the leaves of C. betulus exhibit heightened adaptive plasticity and stability amidst fluctuating environmental conditions. To comprehensively grasp the altitude adaptation of this species, more extensive studies incorporating ecophysiological characteristics are warranted. 

Acknowledgements

This study includes a part of the master's thesis prepared by Fahrettin Atar at Karadeniz Technical University, Institute of Natural Sciences.

Funding

A part of this study was supported by TUBITAK (The Scientific and Technological Research Council of Türkiye, Project Number: TOVAG-112O033). 

Conflict of interest

The authors have no relevant financial or non-financial interests to disclose.

References

Abdusalam A, Li Q (2018) Morphological plasticity and adaptation level of distylous Primula nivalis in a heterogeneous alpine environment. Plant Diversity 40(6): 284−291.

Adamidis GC, Varsamis G, Tsiripidis I, Dimitrakopoulos PG, Papageorgiou AC (2021) Patterns of leaf morphological traits of beech (Fagus sylvatica L.) along an altitudinal gradient. Forests 12(10): 1297.

Akbarian MR, Tabari M, Akbarinia M, Zarafshar M, Meave JA, Yousefzadeh H, Satarian A (2011) Effected of elevation gradiant on leaf and stomata morphology of Caucasian alder (Alnus subcordata) in Hyrcanian forests (Iran). Journal of Folia Oecologica 38: 1−7.

Amer W, Hamdy R, Mahdy R (2016) Leaf micromorphological features of the genus Cordia L.(Boraginaceae) from Egypt. Bangladesh Journal of Plant Taxonomy 23(2): 119-131.

Anşin R, Özkan ZC (2006) Spermatophytha, Woody Plant. Karadeniz Technical University, Faculty of Forestry, Publication No: 167, Faculty Publication No: 19, Trabzon.

Aragón G, Martínez I, García A (2012) Loss of epiphytic diversity along a latitudinal gradient in southern Europe. Science of the total environment 426: 188-195.

Atar E, Atar F, Güney D, Turna İ, Seyis E (2014) Doğu Gürgeni’nde (Carpinus orientalis Miller) Yaprak Karakterlerine Ait Bazı Morfolojik Özelliklerin Yükseltiye Bağlı Olarak Değişimi. II. Ulusal Akdeniz Orman Ve Çevre Sempozyumu, 22-24 Ekim, Isparta, Türkiye.

Atar F, Bayraktar A, Atar E, Turna İ (2017) Morphological diversity of common hornbeam (Carpinus betulus L.) seeds in the Eastern Black Sea Region of Turkey, International Forestry and Environment Symposium, Trabzon, Turkey, 7 - 10 Nov, pp. 203.

Atar F (2022) Population Variability of Common Hornbeam (Carpinus betulus L.) in North-Eastern Part of Turkey According to the Involucre Morphology. Šumarski list 146(11-12): 507-516.

Atar F, Turna, I (2018) Fruit and seedling diversity among sweet chestnut (Castanea sativa Mill.) populations in Turkey. Sumarski List 142(11-12): 611−619.

Atar F, Bayraktar A, Yıldırım N, Turna İ, Güney D (2020) Fruit and seed diversity of Smilax excelsa in the Black Sea Region, Turkey. Turkish Journal of Forestry Research. 7(1): 1-8.

Aykut Y, Emel U, Tekin BM (2017) Morphological variability of evergreen oaks (Quercus) in Turkey. Bangladesh Journal of Plant Taxonomy 24(1): 39-47.

Bayramzadeh V, Funada R, Kubo T (2008) Relationships between vessel element anatomy and physiological as well as morphological traits of leaves in Fagus crenata seedlings originating from different provenances. Trees 22: 217-224.

Bednorz L (2006) Morphological variability of leaves of Sorbus torminalis (L.) Crantz in Poland. Acta Soc Bot Pol 75(3): 233-243.

Beerling DJ, Heath J, Woodward FI, Mansfield TA (1996) Interactions in trees: observations and mechanisms. New Phytologist 134: 235-242.

Boratyński A, Boratyńska K, Mazur M, Marcysiak K (2007) Seed involucre variation in Carpinus betulus (Corylaceae) in Poland. Acta Biologica Cracoviensia Sereies Botanica 49(1): 103-111.

Bresson CC, Vitasse Y, Kremer A, Delzon S (2011) To what extent is altitudinal variation of functional traits driven by genetic adaptation in European oak and beech?. Tree physiology 31(11): 1164-1174.

Caudullo G, Welk E, San-Miguel-Ayanz J (2017) Chorological maps for the main European woody species. Data in brief 12: 662–666. 

Ceplick GP (1995) Genotypic variation and plasticity of clonal growth in relation to nutrient availability in Amphibromus scabrivalvis. J. Ecol 83: 459e468.

Cobanoglu H, Sevik H, Koç İ (2023) Do annual rings really reveal Cd, Ni, and Zn pollution in the air related to traffic density? An example of the Cedar tree. Water, air, & soil pollution 234(2): 65.

Cordell S, Goldstein G, Mueller-Dombois D, Webb D, Vitousek PM (1998) Physiological and morphological variation in Metrosideros polymorpha, a dominant Hawaiian tree species, along an altitudinal gradient: the role of phenotypic plasticity. Oecologia 113: 188-196.

Di Pierro, EA, Mosca E, González-Martínez SC, Binelli G, Neale DB, La Porta N (2017) Adaptive variation in natural Alpine populations of Norway spruce (Picea abies [L.] Karst) at regional scale: Landscape features and altitudinal gradient effects. Forest Ecology and Management 405: 350-359.

Dogan S, Kilicoglu C, Akinci H, Sevik H, Cetin M (2023) Determining the suitable settlement areas in Alanya with GIS-based site selection analyses. Environmental science and pollution research 30(11): 29180-29189.

Douaihy B, Sobierajska K, Jasińska AK, Boratyńska K, Ok T, Romo A, NMachon N, Didukh Y, Dagher-Kharrat MD, Boratyński A (2012) Morphological versus molecular markers to describe variability in Juniperus excelsa subsp. excelsa (Cupressaceae). AoB Plants pls013.

Drvodelić D, Jemrić, T, Oršanić M, Paulić V (2015) Fruits size of wild apple (Malus sylvestris/L./Mill.): impact on morphological and physiological properties of seeds. Sumarki List 139(3-4): 145-153.

Ertan E (2007) Variability in leaf and fruit morphology and in fruit composition of chestnuts (Castanea sativa Mill.) in the Nazilli region of Turkey. Genet. Resour. Crop Ev. 54: 691-699.

Ertürk N, Arıcak B, Şevik H, Yiğit N (2024) Possible Change in Distribution Areas of Abies in Kastamonu due to Global Climate Change. Kastamonu University Journal of Forestry Faculty 24(1): 81-91.

Falconer DS, Mackay TFC (1996) Introduction to quantitative genetics, 4th Edn. London: Pearson Education Ltd.

Furlow JJ (1990). The genera of Betulaceae in the southeastern United States. Journal of the Arnold Arboretum 71(1): 1-67.

Grubb PJ (1977) Control of forest growth and distribution on wet tropical mountains with specific reference to mineral nutrition. Annual Review of Ecology and Systematics 8: 83-107.

Güney D, Turna H, Turna İ, Kulaç Ş, Atar F, Fılız E (2016) Variations within and among populations depending on some leaf characteristics of oriental beech Fagus orientalis Lipsky. Biological Diversity and Conservation 9/2 (S1): 1-9.

Hagemeier M, Leuschner C (2019) Functional crown architecture of five temperate broadleaf tree species: Vertical gradients in leaf morphology, leaf angle and leaf area density. Forests 10: 265.

Halbritter AH, Fior S, Keller I, Billeter R, Edwards PJ, Holderegger R, Karrenberg S, Pluess AR, Widmer A, Alexander JM (2018) Trait differentiation and adaptation of plants along elevation gradients. J. Evolution. Biol. 31(6): 784-800.

Herrera CM, Bazaga P (2013) Epigenetic correlates of plant phenotypic plasticity: DNA methylation differs between prickly and nonprickly leaves in heterophyllous Ilex aquifolium (Aquifoliaceae) trees. Bot. J. Linn. Soc. 171: 441e452.

Hillier J (1991) The Hillier manual of trees and shrubs (6th edn). Melksham, Wiltshire, UK: Redwood Press. 704 p.

Hultine KR, Marshall JD (2000) Altitude trends in conifer leaf morphology and stable carbon isotope composition. Oecologia 123: 32-40.

Huzulák J, Eliáš P (1975) Within-crown pattern of ecophysiological features in leaves of Acer campestre and Carpinus betulus. Folia Geobot. Phytotaxon. 10: 337–350.

Ishii HR, Horikawa SI, Noguchi Y, Azuma W (2018) Variation of intra-crown leaf plasticity of Fagus crenata across its geographical range in Japan. For Ecol Manag 429:437-448.

Jeong N, Moon JK, Kim HS, Kim CG, Jeong SC (2011) Fine genetic mapping of the genomic region controlling leaflet shape and number of seeds per pod in the soybean. Theor Appl Genet 122:865-874. 

Joel G, Aplet G, Vitousek PM (1994) Leaf morphology along environmental gradients in Hawaiian Metrosideros polymorpha. Biotropica 17-22.

Kardel F, Wuyts K, Babanezhad M, Vitharana UWA, Wuytack T, Potters G, Samson R (2010) Assessing urban habitat quality based on specific leaf area and stomatal characteristics of Plantago lanceolata L. Environ Pollut 158:788-794.

Kiaei M, Moosavi V, Ebadi SE (2019) Effects of altitude on density and biometric properties of hornbeam wood (Carpinus betulus). For. Syst. 28(2): e011.

Koç İ, Canturk U, Isinkaralar K, Ozel HB, Sevik H (2024) Assessment of metals (Ni, Ba) deposition in plant types and their organs at Mersin City, Türkiye. Environmental Monitoring and Assessment 196(3): 282.

Kolozsvari AG, Sala F, Nutrition P (2020) Morphometric and physiological analysis of Fagus Sylvatica and Carpinus Betulus leaves. Res. J. Agric. Sci. 52: 66–71.

Körner C (1989) The nutritional status of plants from high elevation: a worldwide comparison. Oecologia 81: 379-391.

Körner C (1999) Alpine plant life: functional plant ecology of high mountain ecosystems. Springer. Berlin Heidelberg, New York. p. 344.

Körner C (2007) The use of altitude in ecological research. Trends in Ecology & Evolution 22: 567-574.

Körner C, Bannister P, Mark AF (1986) Altitudinal variation in stomatal conductance, nitrogen content and leaf anatomy in different plant life forms in New Zealand. Oecologia 69: 577-588.

Krüssmann, G (1984) Manual of Cultivated Broad-Leaved Trees and Shrubs, volume 1, A-D, Timber Press Inc., Portland, Oregon.

Lee C, Kim J, Darshetkar AM, Choudhary RK, Park SH, Lee J, & Choi S (2018) Mericarp morphology of the tribe Selineae (Apiaceae, Apioideae) and its taxonomic implications in Korea. Bangladesh Journal of Plant Taxonomy 25(2): 175-186.

Leingärtner A, Hoiss B, Krauss J, Steffan-Dewenter I (2014) Combined effects of extreme climatic events and elevation on nutritional quality and herbivory of alpine plants. PloS one 9(4): e93881.

Li PC, Skvortsov AK (1999) Betulaceae. In:Wu, Z.Y., Raven, P., Hong, D. (Eds.), Flora of China. Science Press/Missouri Botanical Garden Press, Beijing/St. Louis 4: 289-300.

MacKey J (1988) A plant breeder’s aspect on the taxonomy of cultivated plants. Biologisches Zentralblatt 107: 369-379.

Matesanz S, Gianoli E, Valladares F (2010) Global change and the evolution of phenotypic plasticity in plants. Ann. N.Y. Acad. Sci. 1206: 35e55.

Nascimbene J, Marini L (2015) Epiphytic lichen diversity along elevational gradients: biological traits reveal a complex response to water and energy. J. Biogeogr. 42: 1222e1232.

Nicotra AB, Atkin OK, Bonser SP, Davidson AM, Finnegan EJ, Mathesius U, Poot P, Purugganan MD, Richards CL, Valladares F, van Kleunen M (2010) Plant phenotypic plasticity in a changing climate. Trends Plant Sci. 15: 684e692.

Ogaya R, Peñuelas J (2007) Leaf mass per area ratio in Quercus ilex leaves under a wide range of climatic conditions. The importance of low temperatures. J Acta Oecologica 31: 168-173.

Özbucak TB, Akçin ÖE, Ertürk Ö (2013) The change in ecological, anatomical and antimicrobiological properties of the medicinal plant Tilia rubra dc. subsp. caucasica (rupr.) v. engler along an elevational gradient. Pakistan Journal of Botany 45(5): 1735-1742.

Özel HB, Şevik H, Yıldız Y, Çobanoğlu H (2024) Effects of Silver Nanoparticles on Germination and Seedling Characteristics of Oriental Beech (Fagus orientalis) Seeds. BioResources 19(2): 2135-2148.

Özen HÇ, Onay A (2013) Bitki fizyolojisi. 2. Basım, Nobel Akademik Yayıncılık, s. 349.

Öztürk M, Bolat İ, Ergün A (2015) Influence of air–soil temperature on leaf expansion and LAI of Carpinus betulus trees in a temperate urban forest patch. Agricultural and Forest Meteorology 200: 185-191.

Panek JA, Waring RH (1995) Carbon isotope variation in Douglas-fir foliage: improving the δ13C–climate relationship. Tree Physiology 15: 657.

Paridari IC, Jalali SG, Sonboli A, Zarafshar M, Bruschi P (2013) Leaf macro-and micro-morphological altitudinal variability of Carpinus betulus in the Hyrcanian forest (Iran). J. Forestry Res. 24(2): 301-307.

Pereira-Lorenzo S, Fernandez-Lopez J, Moreno-González J (1996) Variability and grouping of Northwestern Spanish chestnut cultivars. I. Morphological traits. Journal of the American Society for Horticultural Science 121(2): 183-189.

Premoli AC, Raffaele E, Mathiasen P (2007) Morphological and phenological differences in Nothofagus pumilio from contrasting elevations: evidence from a common garden. Austral Ecology 32(5): 515-523.

Poljak I, Kajba D, Ljubić I, Idžojtić M (2015) Morphological variability of leaves of Sorbus domestica L. in Croatia. Acta Soc. Bot. Pol. 84(2): 249-259.

Poljak I, Idžojtić M, Šapić I, Korijan P, Vukelić P (2018) Diversity and structure of Croatian continental and Alpine-Dinaric populations of grey alder (Alnus incana /L./ Moench subsp. incana): Isolation by distance and environment explains phenotypic divergence. Sumarski List 142(1-2): 35-48.

Poljak I, Idžojtić M, Zebec M, Perković N (2012) The variability of European sweet chestnut (Castanea sativa Mill.) in the region of northwest Croatia according to morphology of fruits. Sumarski List 136(9-10): 479-489.

Poljak I, Vahčić N, Gačić M, Idžojtić M (2016) Morphological characterization and chemical composition of fruits of the traditional Croatian chestnut variety ‘Lovran Marron’. Food Technol. Biotechnol 54(2): 189-199.

Poljak I, Vahčić N, Liber Z, Šatović Z, Idžojtić M (2022) Morphological and chemical variation of wild sweet chestnut (Castanea sativa Mill.) populations. Forests 13(1): 55.

Poorter H, Niinemets U, Poorter L, Wright IJ, Villar R (2009) Causes and consequences of variation in leaf mass per area (LMA): a meta-analysis. New Phytol 182:565-588.

Savolainen O, Pyhäjärvi T, Knürr T (2007) Gene flow and local adaptation in trees. Annu. Rev. Ecol. Evol. Syst. 38: 595-619.

Scheepens JF, Frei ES, Stöcklin J (2010) Genotypic and environmental variation in specific leaf area in a widespread Alpine plant after transplantation to different altitudes. Oecologia 164:141-150.

Schlichting CD (1986) The evolution of phenotypic plasticity in plants. Ann Rev Ecol Syst 17:667-693.

Sękiewicz K, Boratyńska K, Dagher-Kharrat MB, Ok T, Boratyński A (2016) Taxonomic differentiation of Cupressus sempervirens and C. atlantica based on morphometric evidence. Syst Biodivers 14(5): 494-508.

Sevik H, Cetin M, Kapucu O, Aricak B, Canturk U (2017) Effects of light on morphologic and stomatal characteristics of Turkish Fir needles (Abies nordmanniana subsp. bornmulleriana Mattf.). Fresenius Environmental Bulletin 26(11): 6579-6587.

Stiegel S, Mantilla-Contreras J (2018) Experimental Study of Environmental Effects: Leaf Traits of Juvenile Fagus sylvatica, Acer pseudoplatanus, and Carpinus betulus Are Comparable to Leaves of Mature Trees in Upper Canopies. Int. J. Ecol. 1–8.

Sultan SE (2000) Phenotypic plasticity for plant development, function and life history. Trends Plant Sci 5:537-542.

Suszka B, Muller C, Bonnet-Masimbert M (1998) Seeds of forest broadleaves: from harvest to sowing. Paris: Institut National de la Recerche Agronomique. 294 p.

Tekin O, Cetin M, Varol T, Ozel HB, Sevik H, Zeren Cetin I (2022) Altitudinal migration of species of Fir (Abies spp.) in adaptation to climate change. Water, air, & soil pollution 233(9): 385.

Varol T, Canturk U, Cetin M, Ozel HB, Sevik H, Zeren Cetin, I (2022) Identifying the suitable habitats for Anatolian boxwood (Buxus sempervirens L.) for the future regarding the climate change. Theoretical and Applied Climatology 150(1): 637-647.

Varol T, Cetin M, Ozel HB, Sevik H, Zeren Cetin I (2022) The effects of climate change scenarios on Carpinus betulus and Carpinus orientalis in Europe. Water, Air, & Soil Pollution, 233(2), 45.

Vitousek P (1982) Nutrient cycling and nutrient use effenciency. Amer. Nat. 119: 553-572.

Vitousek PM, Aplet G, Turner D, Lockwood JJ (1992) The Mauna Loa environmental matrix: foliar and soil nutrients. Oecologia 89(3): 372-382.

Wang Z, Zhang L (2012) Leaf shape alters the coefficients of leaf area estimation models for Saussurea stoliczkai in central Tibet. Photosynthetica 50:337-342.

Weiher E, Keddy PA (1999) Assembly rules as general constraints on community composition. In: Weiher E, Keddy PA (eds) Ecological assembly rules: perspectives, advances, retreats. Cambridge University Press, Cambridge, pp 251-271.

Wells CL, Pigliucci M (2000) Adaptive phenotypic plasticity: the case of heterophylly in aquatic plants. Perspect Plant Ecol Evol Syst 3:1-18.

Wright IJ, Westoby M (2002) Leaves at low versus high rainfall: coordination of structure, lifespan and physiology. New Phytol 155:403-416.

Wright IJ, Reich PB, Westoby M, Ackerly DD, Baruch Z, Bongers F, Cavender-Bares J, Chapin T, Cornelissen JH, Diemer M (2004) The worldwide leaf economics spectrum. Nature 428:821.

Yaltırık F (1982) Flora of Turkey and Aest Eagen Island, University Press, Edinburgh.

Yigit N, Öztürk A, Sevik H, Özel HB, Ramadan Kshkush FE, Işık B (2023) Clonal Variation Based on Some Morphological and Micromorphological Characteristics in the Boyabat (Sinop/Turkey) Black Pine (Pinus nigra subsp. pallasiana (Lamb.) Holmboe) Seed Orchard. BioResources 18(3): 4850-4865.

Yousefzadeh H, Tabari M, Akbarinia M, Akbarian MR, Bussotti F (2010) Morphological plasticity of Parrotia persica leaves in eastern Hyrcanian forests (Iran) is related to altitude. Nordic Journal of Botany 28: 344-349.

Zarafshar M, Akbarinia M, Bruschi P, Hosseiny SM, Yousefzadeh H, Taieby M, Sattarian A (2010) Phenotypic variation in chestnut (Castanea sativa Mill.) natural populations in Hyrcanian forest (north of Iran), revealed by leaf morphometrics. Folia Oecologica 37: 113-121.

Zotz G, Wilhelm K, Becker A (2011) Heteroblasty-A Review. Bot Rev 77:109-151.