DIGITALNA ARHIVA ŠUMARSKOG LISTA
prilagođeno pretraživanje po punom tekstu




ŠUMARSKI LIST 1-2/2018 str. 21     <-- 21 -->        PDF

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
RAZNOLIKOST I STRUKTURIRANOST HRVATSKIH KONTINENTALNIH I ALPSKO-DINARSKIH POPULACIJA BIJELE JOHE (Alnus incana /L./ Moench subsp. incana): GEOGRAFSKA I OKOLIŠNA IZOLACIJA KAO UZROK FENOTIPSKE DIVERGENCIJE
Igor POLJAK, Marilena IDŽOJTIĆ, Irena ŠAPIĆ, Patrik KORIJAN, Joso VUKELIĆ
Summary
We studied the morphological variation and its correlation to the environment in a boreal tree species, Alnus incana (L.) Moench subsp. incana, across the Croatian continental and Alpine-Dinaric biogeographical regions. A total of seven grey alder populations from environmentally divergent habitats were included in the study. We combined descriptive and multivariate statistical methods by using morphological leaf traits and environmental features to examine the diversity and structure of grey alder populations. High phenotypic variation, and two distinct morphotypes were found: small-leaf in the Alpine-Dinaric region, and large-leaf in the continental region. We identified a biogeographical structuring of populations with a high level of among-tree variation within the populations. Multivariate phenotypic and environmental analysis confirmed the existence of small- and large-leaf ecotypes. Additionally, a lack of significant among-population variation within biogeographical regions was observed. Furthermore, we found that phenotypic divergence of the studied grey alder populations can be explained as a result of a significant level of isolation both by distance and by environment.
Key words: phenotypic variation, leaf morphology, adaptation, ecological divergence, clinal variation, altitudinal gradient

ŠUMARSKI LIST 1-2/2018 str. 22     <-- 22 -->        PDF

Introduction
Uvod
Alnus incana (L.) Moench subsp. incana (hereafter Alnus incana), family Betulaceae, is a diploid (2n = 2x = 28), monoecious and wind-pollinated, relatively short-lived, small to medium-sized deciduous tree (Tallantire 1974; Mandák et al. 2016; Vít et al. 2017). It is a light-demanding, hygro-mesophilous species which tolerates low winter temperatures. The grey alder can be found on stream banks, lake shores, damp meadows and in wet mountain environments at elevations of up to 1800 m (Schütt and Lang 2014; Houston Durrant et al. 2016). Along with all other alder species, it is unique within the Betulaceae family since it can fix nitrogen in symbiotic root nodules with the actinomycetes bacteria from the genus Frankia (Benson et al. 2004; Vanden Heuvel 2011). In addition, it is an important pioneer species with the ability of fast regeneration and colonization of disturbed or damaged areas. Similar to the Norway spruce (Picea abies /L./ H.Karst), the range of the grey alder is divided into a northern (subsp. kolaensis /N.I.Orlova/Á.Löve et D.Löve) and a southern area, which meet in the Polish lowlands (Douda et al. 2014). The distribution of subsp. incana in the southern part of the range is linked with a patchy mountain occurrence (in the Alps, the northern Apennines, the Hercynian mountains, the Carpathians, the Bulgarian mountains, the Dinaric Alps, the Caucasus, and mountains in Turkey) and continues eastward across European Russia to western Siberia (Mandák et al. 2016). Grey alder is represented by two subspecies in North America: subsp. rugosa (Du Roi) R.T.Clausen and subsp. tenuifolia (Nutt.) Breitung.
In Croatia, A. incana inhabits north-western and north-eastern parts of the country at the southern border of its natural range of the Central European mountains (Vukelić and Rauš 1998; Vukelić 2012; Trinajstić 2008; Poljak et al. 2014). This part of the grey alder’s natural distribution is completely separated from its southern natural range in the Balkan Peninsula. Here, the grey alder appears in two biogeographical regions with contrasting climates: the continental region, along the course of the Drava river; and in the mountainous Alpine-Dinaric region in Gorski kotar, along the course of the Kupa river and its tributaries. This biogeographical divergence caused by strong climatic gradients and complex topography resulted in a clear floristic differentiation of grey alder stands in Croatia.
In the continental region of Croatia, the grey alder occurs mainly in riparian and floodplain forests along the main watercourse of the river Drava, where it forms smaller and isolated populations. These populations are included within the association Equiseto hyemali-Alnetum incanae Moor 1958 (Trinajstić 1964; Franjić et al. 1999). This association is characterized by the presence of hygrophilous species, such as: Fraxinus angustifolia Vahl subsp. oxycarpa (Willd.) Franco et Rocha Afonso, Ulmus laevis Pall., U. minor Mill., Prunus padus L., Equisetum hyemale L., and others. The former forest areas that were studied by Trinajstić (1964) are today no longer covered by forest. Riparian and floodplain forests close to urban centres are being dried out, and their habitats are transformed into residential areas and municipal infrastructure so that grey alder forests along the course of the Drava river are nowadays in strong regression (Vukelić 2012; Poljak et al. 2014).
Stands from the north-western Dinarides are defined within the association Lamio orvalae-Alnetum incanae Dakskobler 2010, where a geographic variant with Helleborus dumetorum is described (Vukelić et al. 2012, 2017). In general, the main characteristic of these stands is floristic heterogeneity and high richness of species, mainly as a consequence of the biogeographical position and floristic development of the western Dinarides (Bennett et al. 1991; Petit et al. 2003; Magri et al. 2006; Liepelt et al. 2009; Temunović et al. 2013) – one of the floristically richest forest areas in Europe. The following species are particularly prominent: Fraxinus excelsior L., Ulmus glabra Huds., Fagus sylvatica L., Acer pseudoplatanus L., Lamium orvala L., Scopolia carniolica Jacq., Helleborus dumetorum Waldst. et Kit., Lunaria rediviva L., Knautia drymeia Heuff. subsp. drymeia, and others. The occurrence and growth of grey alder forests in this region is conditioned by occasional flooding of the river Kupa and its tributaries over a length of about thirty kilometres. The river partially flows through a canyon and partially extends along to horizontally elongated terraces, on which forests were cut down in the past. In the last 50 years, land cultivation has been gradually abandoned, giving place to natural succession of vegetation.
To date, several genetic studies have been carried out on the grey alder. Microsatellite primers specific for the A. glutinosa (L.) Gaertn and A. incana, developed by Lepais and Bacles (2011), were used to evaluate the population genetic structure and diversity of 65 populations, and to reconstruct the historical pattern of postglacial biogeographical range expansion of the boreal tree species A. incana in Europe (Mandák et al. 2016). Furthermore, genetic diversity and clonal structure of 24 populations covering the Central European part of the species range from the Scandinavian Peninsula to the Balkans were analysed utilizing four nuclear microsatellite markers (Dering et al. 2016). Additionally, several studies have addressed the morphological variability of grey alder populations (Krauze-Michalska and Boratyńska 2013; Poljak et al. 2014). Krauze-Michalska and Boratyńska (2013) revealed a high level of morphological variation of leaves and significant differences between the Scandinavian and all other European populations of A. incana. By contrast, in our previous research of grey alder populations in Croatia (Poljak et al. 2014), along the upper course of the Drava river, the absence of inter-population variability was observed. However, no research of the grey alder’s morphologic variation on a larger geographical scale in Croatia has been performed thus far.

ŠUMARSKI LIST 1-2/2018 str. 23     <-- 23 -->        PDF

In this study, we combined descriptive and multivariate statistical methods using the morphological leaf traits and environmental features to examine the diversity of grey alder populations distributed across the Croatian continental and Alpine-Dinaric biogeographical regions. Our main objectives were: 1) to reveal the levels of population variation and divergence within and among natural populations of A. incana; and 2) to test correlations between geographic, environmental and morphological variation.
Material and methods
Materijal i metode
Plant material and leaf characters studied – Biljni materijal i istraživane značajke lista
The material for morphometric analysis was collected in July 2016 in seven natural populations of A. incana in Croatia (Figure 1). The study encompassed four populations from the Alpine-Dinaric region (P01 - Gerovo, P02 – Turke, P03 – Čedanj; P04 – Lukovdol), and three populations from the continental region (P05 – Varaždin, P06 – Slatina; P07 – Donji Miholjac). Each of the studied populations was represented by 20 trees and each tree by 20 healthy and undamaged leaves of short shoots from the external, sunlit part of the crown. The minimum distance between individual trees was at least 50 m, to avoid the sampling of close relatives. Vouchers for the populations studied were deposited in the herbarium at the Faculty of Forestry of the University of Zagreb (DEND).
The leaves were scanned and measured using the WinFolia program (WinFolia TM2001). A total of nine characteristics were measured (Figure 2). The accuracy of measurements was 0.1 mm, and the following morphological characteristics

ŠUMARSKI LIST 1-2/2018 str. 24     <-- 24 -->        PDF

were measured on each leaf: leaf blade area (LA); leaf blade length (LL); maximum leaf blade width (MLW); leaf blade length, measured from the leaf base to the point of maximum leaf width (PMLW); leaf blade width at 50% of leaf blade length (LW1); leaf blade width at 90% of leaf blade length (LW2); angle closed by the main leaf vein and the line defined by the leaf blade base and a point on the leaf margin, at 10% of leaf blade length (LA1); angle closed by the main leaf vein and the line defined by leaf blade base and a point on the leaf margin, at 25% of leaf blade length (LA2); and petiole length (PL). Finally, a total of 2800 leaves, dried in a herbarium, were measured, and 25 200 simple data values were obtained.
Statistical analyses – Statističke analize
The measured morphological characteristics were described by standard descriptive statistical parameters: arithmetic mean, standard deviation, coefficient of variation (CV%), and the percentiles 0% (minimum), 5%, 25% (lower quartile), 50% (median), 75% (upper quartile), 95% and 100% (maximum). To assess the possibility of conducting multivariate statistical analyses and parametric tests, the symmetry, unimodality and homoscedasticity of data were verified (Sokal and Rohlf 2012). Assumptions of normality were checked using the Shapiro–Wilk test, and the assumption of homogeneity of variance using Levene’s test. A hierarchical analysis of variance (ANOVA) was performed to examine the partition of phenotypic variation between the continental and the Alpine-Dinaric regions, among populations within the regions, and within populations.
The relationship between average values of morphological leaf characteristics and geographical longitude, latitude, and altitude (e.g. Krauze-Michalska and Boratyńska 2013) and the influence of climatic variables on leaf characteristics were tested using Spearman’s coefficient (Sokal and Rohlf 2012). Climate data were obtained from the WorldClim database with a spatial resolution close to a square km (Hijmans et al. 2005). To evaluate the correlation between multicharacter differences among populations, a Mantel test (Mantel 1967) was performed on the matrices of Euclidean distances. First, the correlations among all 19 WorldClim bioclimatic variables and topographic variables for all presence points were calculated to exclude the highly correlated ones, whilst keeping the variables useful in predicting the distribution limits of trees, such as climatic averages and extremes (Zimmermann et al. 2009; Temunović et al. 2012). We computed and tested the correlations between: (1) the matrix of the geographical distances between pairs of populations and the matrix of morphological differences among populations – “isolation by distance” (Wright 1943); and (2) the matrix of environmental distances and the matrix of morphological differences among populations – “isolation by environmental distance” (Mendez et al. 2010). The significance level was assessed after 10,000 permutations, and the Mantel test was performed with the R package “Vegan” (Oksanen et al. 2017).
Multivariate statistical methods were used to identify the population differentiation (McGarigal et al. 2000). Pearson’s correlation coefficient was used to identify interactions between leaf traits and to detect potential redundant variables, i.e. highly correlated variables were excluded from the analyses. The conducted cluster analysis resulted in a hierarchical tree, where the unweighted pair-group method with arithmetic mean (UPGMA) was used to join the clusters, and the Euclidean distance to define the distance between the studied objects. The K-means method was applied to detect phenotypic structure and define the number of K-groups that best explained the morphological variation of populations (e.g Douaihy et al. 2012; Boratyński et al. 2013; Sobierajska et al. 2016). In addition, the biogeographical structure of the studied populations was further

ŠUMARSKI LIST 1-2/2018 str. 25     <-- 25 -->        PDF

analysed using the Monmonier’s maximum difference algorithm, implemented in Barrier 2.2 software (Manni et al. 2004). A canonical discriminant analysis was performed to calculate the discriminatory power of characters among biogeographical regions. The proportion of correctly classified individuals into biogeographical regions was determined using classificatory discriminant analyses (Oksanen et al. 2017). To examine the level of environmental and phenotypic divergence between the studied biogeographical regions, principal component analysis (PCA) was used. The input data in multivariate statistical methods were previously standardized, i.e. standardization of characters to zero mean and unit standard deviation was performed prior to each multivariate analysis.
The above statistical analyses were conducted using the STATISTICA 8.0 statistical program (StatSoft, Inc. 2001), and “MorphoTools” R scripts in R v.3.2.2 (R Core Team, 2016) following the manual of Koutecký (2015).
Results
Rezultati
The results of the conducted descriptive statistical analysis are shown by biogeographical region (Figure 3). The highest mean values for all analysed characteristics were observed in the populations from the continental region. Coefficients of variation among populations for particular morphological traits were generally similar (data not shown). However,

ŠUMARSKI LIST 1-2/2018 str. 26     <-- 26 -->        PDF

in the populations from the Alpine-Dinaric region, somewhat higher coefficients of variations were observed. The most variable characteristics with the CV above 25% were leaf blade surface (LA = 32.65%) and leaf blade width at 90% of leaf blade length (LW2 = 30.30%). In addition, a high degree of variability was also noted for petiole length

ŠUMARSKI LIST 1-2/2018 str. 28     <-- 28 -->        PDF


ŠUMARSKI LIST 1-2/2018 str. 27     <-- 27 -->        PDF

with the total CV of 24.27%. As opposed to that, the lowest value of the coefficient of variation was found for the angle closed by the main leaf vein and the line defined by the leaf blade base and a point on the leaf margin, at 10% (LA1 = 10.46%) and 25% (LA2 = 7.82%) of leaf blade length.
The results of the hierarchical analysis of variance (ANOVA) for each characteristic are shown in Table 1. The studied regions differ significantly for all studied variables. Likewise, statistically significant differences among trees within populations were observed. The AMOVA analysis showed that most of the morphological diversity was attributable to the differences between regions, confirming the geographical structuring of populations. However, a highly significant percentage of variation was explained by the differences among individuals within populations for the variables: leaf blade length, measured from the leaf base to the point of maximum leaf width (PMLW); leaf blade width at 90% (LW2); petiole length (PL); and leaf angles LA1 and LA2. Among populations within regions, there were significant differences in four of the nine analysed traits (LA, MLW, LW1, PL), where only a small proportion of the total variation was explained.
As expected, strong correlations between leaf morphological traits were observed. Almost all measured leaf traits correlated with each other at a statistically significant level. The highest correlations, considering all populations, were found between leaf blade area (LA); leaf blade length (LL); maximum leaf blade width (MLW); and leaf blade width at 50% of leaf blade length (LW1). We also observed strong and significant correlations between the angles closed by the main leaf vein and the lines defined by leaf blade base and points on the leaf margin, at 10 % (LA1) and 25% (LA2) of leaf blade length. Non-significant correlations were found between petiole length (PL) and the above-mentioned leaf traits (LA1 and LA2).
Using Spearman’s correlation coefficient (Table 2), a highly positive correlation was found between geographical longitude and the six measured variables (LA, LL, MLW, LW1, LA1, LA2). On the other hand, all measured variables were highly negatively correlated with the altitude and total precipitation of the warmest quarter of the year. Likewise, almost all leaf traits were significantly negatively correlated with the distance-to-water, annual mean precipitation, precipitation of the coldest quarter and isothermality. By contrast, almost all leaf traits were positively correlated with the annual mean temperature, maximum temperature of the warmest month and precipitation seasonality.
The Mantel test identified significant correlations between the morphological, geographical, and environmental distance matrices. Correlations were higher between morphological and geographical distance matrices (r = 0.91, p = 0.0018), and slightly smaller but still very strong between the morphological and environmental distance matrices (r = 0.80, p = 0.0029).
The structure of the seven grey alder populations was inferred by the K-means clustering method. The most probable division was detected at K = 2, and the estimated population structure is shown in Figure 1. If the proportion of a certain population was equal to or higher than 0.75, it was assumed that the population belonged to one cluster, and if it was lower than 0.75, it was assumed that the population had a mixed origin. The populations from the Alpine-Dinaric region grouped together into cluster A, and the populations from continental Croatia grouped together into cluster B. Only one population, P05 (Varaždin), was of mixed origin, with the dominant proportion from cluster B. The results obtained with the K-means clustering method were congruent with the hierarchical tree (Figure 4C), where the unweighted pair-group method with arithmetic mean (UPGMA) was used to join the clusters.

ŠUMARSKI LIST 1-2/2018 str. 34     <-- 34 -->        PDF

unutar-populacijska raznolikost. Kontinentalne populacije su imale veće vrijednosti svih mjerenih značajki lista u odnosu na populacije iz alpsko-dinarske regije.
Analizom varijance, utvrđene su statistički značajne razlike između istraživanih regija te između stabla unutar populacija, dok su se populacije unutar regija razlikovale samo za četiri istraživane značajke (tablica 1). Unutarpopulacijska varijabilnost veća je od međupopulacijske varijabilnosti, što je odlika mnogih vrsta drveća. Multivarijatne analize potvrđuju postojanje jasne razlike između kontinentalnih i alpsko-dinarskih populacija bijele johe (slika 1, slika 4, tablica 3), što upućuje na mogućnost da istraživane populacije predstavljaju dva različita ekotipa.
Morfološka udaljenost između populacija bila je korelirana s geografskom i ekološkom udaljenošću, što potvrđuje uzorak „isolation by distance“ i „isolation by environmental distance“. Drugim riječima, utvrđeno je da geografska udaljenost i ekološka divergencija imaju važan utjecaj na oblikovanje morfološke strukture populacija te da je ona posljedica njihove moguće lokalne prilagodbe na okolišne uvjete i njihove geografske udaljenosti (tablica 2).
Dobivene spoznaje o varijabilnosti populacija mogu doprinijeti razvoju učinkovitijih planova očuvanja i gospodarenja sivom johom u Hrvatskoj. Kako bi se potvrdili dobiveni zaključci o varijabilnosti populacija sive johe dobiveni morfometrijskim metodama, istraživanja je svakako potrebno proširiti i na molekularno-biološke metode.
Ključne riječi: fenotipska varijabilnost, morfologija lista, adaptabilnost, ekološka divergencija, klinalna varijabilnost, visinski gradijent

ŠUMARSKI LIST 1-2/2018 str. 29     <-- 29 -->        PDF

The results of the discriminant analysis suggest that the differentiation between the studied biogeographical regions is significant. Although with some overlap, the analysis revealed a clear separation between the Alpine-Dinaric and continental populations (Figure 4D). The proportion of correctly classified individuals into two studied biogeographical regions “Alpine-Dinaric” and “continental” is shown in Figure 1B. The overall classification success was 87.1% of correctly classified individuals. Within the Alpine-Dinaric group, 88.7% of the individuals were included in the correct population group. A lower rate of correctly classified individuals was observed in the continental region (85.0%).
Principal component analysis (PCA) has clearly supported a separation of the continental and Alpine-Dinaric populations along the altitudinal, temperature and precipitation gradients (Figure 4B, Table 3). The results show that continental populations are characterised by low elevations and warmer habitats with higher temperatures and lower precipitation. By contrast, high-elevation populations from the Alpine-Dinaric region are distributed within cooler habitats with lower winter temperatures and higher precipitation. In addition, PCA analysis of all individuals has supported the divergence of two morphologically distinct groups (Figure 4A).
Discussion
Rasprava
The results clearly demonstrate a high phenotypic diversity of grey alder populations in Croatia. In general, populations from the continental region had larger and wider leaves than populations from high altitudes in the mountainous Alpine-Dinaric region. In addition, large-leaf populations showed slightly lower variation possibly due to an overall reduction in the number of individuals, and fragmentation and isolation of populations due to human impact (Vukelić 2012; Poljak et al. 2014). On the other hand, environmental heterogeneity over very short distances within the Alpine-Dinaric region could result in higher phenotypic intra-population variability. A significant variability of morphological characteristics of the grey alder leaves was also reported by Krauze-Michalska and Boratyńska (2013), and Poljak et al. (2014). Nevertheless, in our previous study of grey alder populations, along the upper course of the river Drava, somewhat higher within-population variability was recorded. Furthermore, our results indicated that a natural hybridization has occurred between the common and grey alder in those populations, but at relatively low rates. However, these events may influence the diversity and structure of populations (Barton 2001; Poljak et al. 2017), and increase within-population morphological variability. In the current study, the above-mentioned influence was strongly reduced by examining the plant material for the analysis, i.e. hybrid individuals were subsequently excluded from the study.
The results of the hierarchical analysis of variance (ANOVA) were in line with the expectations of high morphological variation within populations and low differentiation between populations, as observed in alders (Krauze-Michalska and Boratyńska 2013; Poljak et al. 2014) and other wind-pollinated floodplain tree species (Jarni et al. 2011; Zebec et al. 2010, 2014). The relatively high level of among-tree variation within the populations is probably a result of both phenotypic adaptation to specific micro-environmental conditions experienced by each tree, and genetic differentiation among individual trees (Brus et al. 2011; Poljak et al. 2015). However, the AMOVA analysis showed that a great proportion of the total variation was attributable to the differences between regions, confirming the geographical structuring of populations. Moreover, most of the population pairs from both biogeographical regions had non-significant pairwise values (data not shown). The lack of ­significant among-population variation within biogeographical regions could be explained with the population longitudinal distribution along the rivers and floodplains with no barriers to gene flow, where free pollen and seed dispersion between the populations occurs (Temunović et al. 2012). More specifically, this can confirm the findings of our previous study of grey alder populations in Croatia (Poljak et al. 2014), along the upper course of the Drava ­river, where absence of inter-population variability was also observed.
The results of the hierarchical analysis of variance have been confirmed by multivariate statistical methods, suggesting the existence of a clear divergence between the populations from two biogeographical regions. Such differences in variability were also found in other woody species that occur in different biogeographical regions with contrasting climates (Škvorc et al. 2005; Temunović et al. 2012; Poljak et al. 2015; Zebec et al. 2016).
We found that the divergence among populations tended to follow an altitudinal cline along which populations from lower altitudes had larger and wider leaves than populations from high altitudes. In general, it is well known that leaf morphological traits, such as leaf length and width, are negatively correlated with altitude (Körner et al. 1986, 1989; Hovenden and Vander Schoor 2003; Bresson et al. 2011; Paridari et al. 2013). However, many environmental features accompany altitudinal changes to which plant populations adapt: those physically tied to meters above sea level, such as atmospheric pressure, temperature and clear-sky turbidity; and those that are not generally altitude specific, such as moisture, hours of sunshine, wind, season length, geology and even human land use (Körner 2007). In fact,

ŠUMARSKI LIST 1-2/2018 str. 30     <-- 30 -->        PDF

our results showed that the leaf measured characteristics that varied in the grey alder populations were significantly related to climate features. Populations of the grey alder from high altitudes in the mountainous Alpine-Dinaric region occupy habitats with lower temperatures accompanied by high irradiance and larger temperature oscillations. On the other hand, low altitude populations from the continental region are exposed to higher temperatures, lower radiation and higher precipitation variability. In addition, growing season tends to increase toward warmer low-elevation sites. Furthermore, common-garden experiments have shown that populations from high elevations have lower growth rates than populations from high altitudes (Vitasse et al. 2009, 2013; Montesinos-Navarro et al. 2011). Similar trends were observed in our study (authors’ personal observations). This pattern was likely due to differential selection pressure among the biogeographical regions, i.e. populations from low altitudes tend to be adapted to the warmer climate under which selection has favoured a high allocation to growth and competitive ability, whereas populations from high altitudes display lower growth rates and greater cold-tolerance (Körner 2003).
Surprisingly, low- vs high-elevation populations significantly increased leaf area at decreasing annual precipitation and precipitation of the warmest and coldest quarter. Similar results were also obtained by Meier and Leuschner (2008) for Fagus sylvatica stands in central Germany. Authors concluded that stand leaf area of the beech along this precipitation gradient is not a simple function of water availability, but is controlled by several abiotic factors including spring temperature and possibly also nitrogen supply, which both tend to increase toward drier sites, thus overlaying any negative effect of water shortage on leaf development. Since the grey alder is a hygro-mesophilous species which mostly occupied riversides, we assume that differences in precipitation among biogeographical regions did not influence population divergence. Similar conclusions were reported by Krauze-Michalska and Boratyńska (2013). Nevertheless, our results revealed that leaves of grey alder populations decrease in almost all measured leaf traits with increasing distance-to-water. This variation trend can be explained by the specific differences in the geomorphology of the terrain among biogeographical regions. The majority of grey alder populations in the Alpine-Dinaric region inhabit a drier riparian part, on elongated terraces and wet slopes along the course of the river Kupa and its tributaries, which are generally 2 to 6 m above the average water level. By contrast, continental populations of the grey alder, along the river Drava, grow along an active watercourse with intensive shifting of pebbly, and sometimes sandy sediments, on which aerated and humus soils are formed. The presence of different species in these stands also reflects different biogeographical characteristics (Trinajstić 1964; Franjić et al. 1999; Vukelić et al. 2012, 2017).
In general, tree species have high levels of morphological variability (Poljak et al. 2012; McKown et al. 2014; Zebec et al. 2015), and phenotypic differences among populations are often a result of the environmental distances between populations. Nevertheless, our results suggest that phenotypic divergence of the studied grey alder populations is the result of a significant level of isolation both by distance and by the environment. The isolation by distance assumes that phenotypic differentiation between populations increases with increasing geographic distance because of limited pollen and seed dispersal (Wright 1943, 1946; Orsini et al. 2013). On the other hand, gene flow among ecologically divergent habitats is reduced because of reduced success of the establishment of immigrants from different environments, as a result of local genetic adaptation (Nosil and Crespi 2004; Noisl et al. 2005, 2008, 2009; Orsini et al. 2013; DeWoody et al. 2015). However, we cannot exclude the possibility that adaptive differentiation and historical migration processes acted in combination to produce the observed pattern of morphological variation in the studied populations (Temunović et al. 2012; DeWoody et al. 2015). Although the main refugium during the last glacial period was located in Central Europe (Huntley and Birks 1983; Douda et al. 2014; Mandák et al. 2016), A. incana probably survived in situ in river canyons and protected valleys along the north-western Dinaric Alps, while north-eastern populations from the continental biogeographical region could represent newly colonised populations from the refugium in Central Europe.
Conclusions
Zaključci
The study revealed the existence of two morphologically and, to a large extent, geographically and environmentally distinct and well-defined groups of grey alder populations in Croatia. These patterns confirm that multiple evolutionary processes influence the morphological diversity and structure of the populations (DeWoody et al. 2015). In ­general, two distinct morphotypes were found: small-leaf in the mountainous Alpine-Dinaric region, and large-leaf in the continental region. These ecologically divergent ­habitats may have led to the general subdivision into two distinct ecotypes. Adaptation to local environmental conditions in the continental region probably resulted in a higher growth rate and large-leaf morphotypes which increase the species’ competitive ability and long-term success, whereas in the Alpine-Dinaric region natural selection has probably favoured greater cold-tolerance (Körner 2003; ­Vitasse et al. 2009). We also believe that our results could be valuable baseline data for the development of more efficient management plans for this boreal tree species.

ŠUMARSKI LIST 1-2/2018 str. 31     <-- 31 -->        PDF

Acknowledgements
Zahvale
This research was carried out with the financial support of the Ministry of Agriculture of the Republic of Croatia.
References
Literatura
Barton, N.H., 2001: The role of hybridization in evolution, Mol Ecol, 10 (3): 551–568.
Bennett, K.D., P.C. Tzedakis, K.J. Willis, 1991: Quaternary refugia of north European trees, J. Biogeogr, 18: 103–115.
Benson, D.R., B.D. Vanden Heuvel, D. Potter, 2004: Actinorhizal symbioses: diversity and biogeography, In: Plant microbiology (eds. M. Gillings, A. Holmes), Garland Science/BIOS Scientific Publishers, pp. 97–128, Oxford.
Bohonak, A.J., 1999: Gene flow, and population structure, Q Rev Biol, 74: 21–45.
Boratyński, A., A.K. Jasińska, K. Marcysiak, M. Mazur, A.M. Romo, K. Boratyńska, K. Sobierajska, G. Iszkuło, 2013: Morphological differentiation supports the genetic pattern of the geographic structure of Juniperus thurifera (Cupressaceae), Plant Syst Evol, 299: 773–784.
Bresson, C.C., Y. Vitasse, A. Kremer, S. Delzon, 2011: To what extent is altitudinal variation of functional traits driven by genetic adaptation in European oak and beech? Tree Physiol, 31 (11): 1164–1174.
Brus, R., D. Ballian, P. Zhelev, M. Pandža, M. Bobinac, J. Acevski, Y. Raftoyannis, K. Jarni, 2011: Absence of geographical structure of morphological variation in Juniperus oxycedrus L. subsp. oxycedrus in the Balkan Peninsula, Eur J Forest Res, 130: 657–670.
DeWoody, J., H. Trewin, G. Taylor, 2015: Genetic and morphological differentiation in Populus nigra L.: isolation by colonization or isolation by adaptation? Mol Ecol, 24: 2641–2655.
Douda, J., J. Doudová, A. Drašnarová, P. Kuneš, V. Hadincová, K. Krak, P. Zákravský, B. Mandák, 2014: Migration patterns of subgenus Alnus in Europe since the Last Glacial Maximum: A systematic review, PLoS ONE 9: e88709.
Douaihy, B., K. Sobierajska, A.K. Jasińska, K. Boratyńska, T. Ok, A. Romo, N. Machon, Y. Didukh, M.B. Dagher-Kharrat, A. Boratyński, 2012: Morphological versus molecular markers to describe variability in Juniperus excelsa subsp. excelsa (Cupressaceae), AOB Plants 2012: plr003.
Franjić, J., I. Trinajstić, Ž. Škvorc, M. Presečan, I. Samardžić 1999: A contribution to the knowledge of the distribution of Equisetum hyemale L. (Equisetaceae) in Croatia, Nat Croat, 8 (2): 95–100.
Hijmans, R.J., S.E. Cameron, J.L. Parra, P.G. Jones, A. Jarvis, 2005: Very high resolution interpolated climate surfaces for global land areas, Int J Climatol, 25: 1965–1978.
Houston Durrant, T., D. de Rigo, G. Caudullo, 2016: Alnus incana in Europe: distribution, habitat, usage and threats. In: J. San-Miguel Ayanz, D. de Rigo, G. Caudullo, T. Houston Durrant, A. Mauri, (Eds.), European Atlas of Forest Tree Species. Publ. Off. EU, 66-67 pp., Luxembourg.
Hovenden, M.J., J.K. Vander Schoor, 2003: Nature vs nurture in the leaf morphology of Southern beech, Nothofagus cunninghamii (Nothofagaceae), New Phytol, 161: 585–594.
Huntley, B., H. Birks, 1983: An atlas of past and present pollen maps for Europe: 0–13000 years ago. Cambridge University Press., Cambridge.
Jarni, K., M. Westergren, H. Kraigher, R. Brus, 2011: Morphological variability of Fraxinus angustifolia Vahl in the north-western Balkans, Acta Soc Bot Pol, 80 (3): 245–252.
Körner, C., 2003: Alpine plant life: functional plant ecology of high mountain ecosystems. Springer-Verlag, Berlin, Heidelberg, and New York.
Körner, C., 2007: The use of ‘altitude’ in ecological research, Trends Ecol Evol, 22(11): 569–574.
Körner, C., M. Neumayer, S. Menendez-Riedl, A. Smeets-Scheel, 1989: Functional morphology of mountain plants, Flora, 182: 353–383.
Körner, C., P. Bannister, A.F. Mark, 1986: Altitudinal variation in stomatal conductance, nitrogen content and leaf anatomy in different plant life forms in New Zealand, Oecologia, 69: 577–588.
Koutecký, P., 2015: MorphoTools: a set of R functions for morphometric analysis, Plant Syst Evol, 301: 1115–1121.
Krauze-Michalska, E., K. Boratyńska, 2013: European geography of Alnus incana leaf variation, Plant Biosyst, 147 (3): 601-610.
Lepais, O., C.F.E. Bacles, 2011: De novo discovery and multiplexed amplification of microsatellite markers for black alder (Alnus glutinosa) and related species using SSR-enriched shotgun pyrosequencing, J Hered, 102: 627–632.
Liepelt, S., R. Cheddadi, J.L. de Beaulieu, B. Fady, D. Gömöry, D. Hussendörfer, M. Konnert, T. Litt, R. Longauer, R. Terhürne-Berson, B. Ziegenhagen, 2009: Postglacial range expansion and its genetic imprints in Abies alba (Mill.) - A synthesis from palaeobotanic and genetic data, Rev Palaeobot Palyno, 153: 139–149.
Magri, D., G.G. Vendramin, B. Comps, I. Dupanloup, T. Geburek, D. Gömöry, M. Latałowa, T. Litt, L. Paule, J.M. Roure, I. Tantau, W.O. van der Knaap, R. Petit, J.L. de Beaulieu, 2006: A new scenario for the quaternary history of European beech populations: palaeobotanical evidence and genetic consequences, New Phytol, 171:199–221.
Mandák, B., A. Havrdová, K. Krak, V. Hadincová, P. Vít, P. Zákravský, J. Douda, 2016: Recent similarity in distribution ranges does not mean a similar postglacial history: A phylogeographical study of the boreal tree species Alnus incana based on microsatellite and chloroplast DNA variation, New Phytol, 210: 1395–1407.
Manni, F., E. Guérard, E. Heyer, 2004: Geographic patterns of (genetic, morphologic, linguistic) variation: how barriers can be detected by ‘Monmonier’s algorithm’, Hum Biol, 76: 173–190.
Mantel, N., 1967: The detection of disease clustering and a generalized regression approach, Cancer Res, 27: 209–220.
McGarigal, K., S. Cushman, S. Stafford, 2000: Multivariate statistics for wildlife and ecology research, Springer Verlag, 283 pp., New York.
McKown, A.D., R.D. Guy, J. Klápště, A. Geraldes, M. Friedmann, Q.C.B. Cronk, Y.A. El-Kassaby, S.D. Mansfield, C.J. Douglas, 2014: Geographical and environmental gradients shape phenotypic trait variation and genetic structure in Populus trichocarpa, New Phytol, 201: 1263–1276.

ŠUMARSKI LIST 1-2/2018 str. 32     <-- 32 -->        PDF

Meier, I.C., C. Leuschner, 2008: Leaf size and leaf area index in Fagus sylvatica forests: competing effects of precipitation, temperature, and nitrogen availability, Ecosystems 11: 655–669.
Mendez, M., H.C. Rosenbaum, A. Subramaniam, C. Yackulic, P. Bordino, 2010: Isolation by environmental distance in mobile marine species: molecular ecology of franciscana dolphins at their southern range, Mol Ecol 19: 2212–2228.
Montesinos-Navarro, A., J. Wig, F.X. Pico, S.J. Tonsor, 2011: Arabidopsis thaliana populations show clinal variation in a climatic gradient associated with altitude, New Phytol, 189: 282–294.
Nosil, P., B.J. Crespi, 2004: Does gene flow constrain adaptive divergence or vice versa? A test using ecomorphology and sexual isolation in Timema cristinae walking-sticks, Evolution, 58: 102–112.
Nosil, P., D.J. Funk, D. Ortiz-Barrientos, 2009: Divergent selection and heterogeneous genomic divergence, Mol Ecol, 18: 375–402.
Nosil, P., S.P. Egan, D.J. Funk, 2008: Heterogeneous genomic differentiation between walking-stick ecotypes: ‘isolation by adaptation’ and multiple roles for divergent selection, Evolution, 62: 316–336.
Nosil, P., T.H. Vines, D.J. Funk, 2005: Perspective: reproductive isolation caused by natural selection against immigrants from divergent habitats, Evolution, 59: 705–719.
Oksanen, J., F.G. Blanchet, M. Friendly, R. Kindt, P. Legendre, D. McGlinn, P.R. Minchin, R.B. O’Hara, G.L. Simpson, P. Solymos, M.H.H. Stevens, E. Szoecs, H. Wagner, 2017: Package ‘vegan’: community ecology package. Version 2.4-2, URL https://cran.r-project.org, https://github.com/vegandevs/vegan.
Orsini, L., J. Vanoverbeke, I. Swillen, J. Mergeay, L. De Meester, 2013: Drivers of population genetic differentiation in the wild: isolation by dispersal limitation, isolation by adaptation and isolation by colonization, Mol Ecol, 22: 5983–5999.
Paridari, I.C., S.G. Jalali, A. Sonboli, M. Zarafshar, P. Bruschi, 2013: Leaf macro- and micro-morphological altitudinal variability of Carpinus betulus in the Hyrcanian forest (Iran), J Forestry Res, 24 (2): 301-307.
Petit R.J., I. Aguinagalde, J.L. de Beaulieu, C. Bittkau, S. Brewer, R. Cheddadi, R. Ennos, S. Fineschi, D. Grivet, M. Lascoux, A. Mohanty, G. Müller-Starck, B. Demesure-Musch, A. Palmé, J.P. Martin, S. Rendell, G.G. Vendramin, 2003: Glacial refugia: hotspots but not meltingpots of genetic diversity, Science, 300: 1563–1565.
Poljak, I., D. Kajba, I. Ljubić, M. Idžojtić, 2015: Morphological variability of leaves of Sorbus domestica L. in Croatia, Acta Soc Bot Pol, 84 (2): 249–259.
Poljak, I., M. Idžojtić, I. Šapić, J. Vukelić, M. Zebec, 2014. Population variability of grey (Alnus incana /L./ Moench) and black alder (A. glutinosa /L./ Gaertn.) in the Mura and Drava region according to the leaf morphology, Sumar List, 138 (1-2): 7–17.
Poljak, I., M. Idžojtić, Z. Šatović, M. Ježić, M. Ćurković-Perica, B. Simovski, J. Acevski, Z. Liber, 2017: Genetic diversity of the sweet chestnut (Castanea sativa Mill.) in Central Europe and the western part of the Balkan Peninsula and evidence of marron genotype introgression into wild populations, Tree Genet Genomes, 13:18.
Poljak, I., M. Idžojtić, M. Zebec, N. Perković, 2012: The variability of European sweet chestnut (Castanea sativa Mill.) in the region of northwest Croatia according to morphology of fruits, Sumar List, 136 (9-10): 479–489.
R Core Team, 2016. R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria. http://www.R-project.org/
Schütt, P., U.M. Lang, 2014: Alnus incana. Enzyklopädie der Holzgewächse: Handbuch und Atlas der Dendrologie, pp. 1–10.
Škvorc, Ž., J. Franjić, M. Idžojtić, 2005: Population structure of Quercus pubescens Willd. (Fagaceae) in Croatia according to morphology of leaves, Acta Bot Hung, 47: 183–-196.
Sobierajska, K., K. Boratyńska, A. Jasińska, M. Dering, T. Ok, B. Douaihy, M.B. Dagher-Kharrat, A. Romo, A. Boratyński, 2016: Effect of the Aegean Sea barrier between Europe and Asia on differentiation in Juniperus drupacea (Cupressaceae), Bot J Linn Soc, 180 (3): 365–385.
Sokal, R.R., F.J. Rohlf, 2012: Biometry: the principles and practice of statistics in biological research, 4th edition, W.H. Freeman and Co., 937 pp., New York.
StatSoft, Inc. 2001: STATISTICA (data analysis software system), version 8.0.
Tallantire, P.A., 1974: The palaeohistory of the grey alder (Alnus incana (L.) Moench) and black alder (A. glutinosa (L.) Gaertn.) in Fennoscandia, New Phytol, 73: 529–546.
Temunović, M., J. Franjić, Z. Šatović, M. Grgurev, N. Frascaria-Lacoste, J.F. Fernández-Manjarrés, 2012: Environmental heterogeneity explains the genetic structure of continental and Mediterranean populations of Fraxinus angustifolia Vahl, PLoS ONE, 7, e42764.
Temunović, M., N. Frascaria-Lacoste, J. Franjić, Z. Šatović, J.F. Fernández-Manjarrés, 2013: Identifying refugia from climate change using coupled ecological and genetic data in a transitional Mediterranean-temperate tree species, Mol Ecol, 22: 2128–2142.
Trinajstić, I., 1964: Vegetacija obalnog područja rijeke Drave u široj okolici Varaždina, Magistarski rad, Prirodoslovno-matematički fakultet Sveučilišta u Zagrebu.
Trinajstić, I., 2008: Biljne zajednice Republike Hrvatske, Akademija šumarskih znanosti Zagreb, 179 pp., Zagreb.
Vanden Heuvel, B.D., 2011: Alnus. In: Wild Crop Relatives: Genomic and Breeding Resources, Forest Trees (ed. C. Kole), Springer-Verlag Berlin Heidelberg, pp. 1–14.
Vít, P., J. Douda, K. Krak, A. Havrdová, B. Mandák, 2017: Two new polyploid species closely related to Alnus glutinosa in Europe and North Africa – An analysis based on morphometry, karyology, flow cytometry and microsatellites, Taxon, 66 (3): 567–583.
Vitasse, Y., A. Lenz, C. Kollas, C.F. Randin, G. Hoch, C. Korner, 2013: Genetic vs. non-genetic responses of leaf morphology and growth to elevation in temperate tree species, Funct Ecol, 28: 243–252.
Vitasse, Y., S. Delzon, C.C. Bresson, R. Michalet, A. Kremer, 2009: Altitudinal differentiation in growth and phenology among populations of temperate tree species in a common garden, Can J For Res, 39: 1259–1269.
Vukelić, J., 2012: Šumska vegetacija Hrvatske, Sveučilište u Zagrebu, 403 pp., Zagreb.
Vukelić, J., D. Baričević, I. Šapić, 2012: Phytocoenological chahacteristics of forests of grey alder (Alnus incana /L./ Moench) in Gorski kotar, Nat Croat, 21 (1): 49–64.
Vukelić, J., Đ. Rauš, 1998: Šumarska fitocenologija i šumske zajednice u Hrvatskoj, Šumarski fakultet Sveučilišta u Zagrebu, 310 pp., Zagreb.

ŠUMARSKI LIST 1-2/2018 str. 33     <-- 33 -->        PDF

Vukelić, J., I. Šapić, A. Alegro, V. Šegota, I. Stankić, D. Baričević, 2017: Phytocoenological analysis of grey alder (Alnus incana L.) forests in the Dinarides of Croatia and their relationship with affiliated communities, Tuexenia, 37: 65–78.
WinFolia TM, 2001: Regent Instruments Inc., Quebec, Canada, version PRO 2005b.
Wright, S., 1943: Isolation by distance, Genetics, 28: 139–156.
Wright, S., 1946: Isolation by distance under diverse systems of mating, Genetics, 31: 39–59.
Zebec, M., M. Idžojtić, I. Poljak, 2014: Morphological variability of the field elm (Ulmus minor Mill. sensu latissimo) in continental Croatia, Sumar List, 138 (11-12): 563–572.
Zebec, M., M. Idžojtić, I. Poljak, I. Mihaldinec, 2010: The variability of field elm (Ulmus minor Mill. sensu latissimo) in Croatian Drava river valley according to the leaf morphology, Sumar List, 134 (11-12): 569–580.
Zebec, M., M. Idžojtić, I. Poljak, I. Modrić, 2015: Population variability of wych elm (Ulmus glabra Huds.) in the mountainous region of Croatia according to the leaf morphology, Sumar List 139 (9-10): 429–439.
Zebec, M., M. Idžojtić, Z. Šatović, I. Poljak, Z. Liber, 2016: Alive and kicking, or, living on borrowed time? – Microsatellite diversity in natural populations of the endangered Ulmus minor Mill. sensu latissimo from Croatia, Acta Bot Croat, 75 (1): 53–59.
Zimmermann, N.E., N.G. Yoccoz, T.C. Edwards, E.S. Meier, W. Thuiller, A. Guisan, D.R. Schmatz, P.B. Pearman, 2009: Climatic extremes improve predictions of spatial patterns of tree species, P Natl Acad Sci USA 106: 19723–19728.
Sažetak
Bijela, odnosno siva joha, Alnus incana (L.) Moench subsp. incana, listopadna je, jednodomna i anemofilna drvenasta vrsta iz porodice Betulaceae. Prirodno je rasprostranjena u sjevernoj, srednjoj i mjestimično južnoj Europi, zapadnom Sibiru i Turskoj, od nizinskoga pojasa do preko 1800 m nadmorske visine u alpskom području. Bijela joha je mezofilna do higrofilna vrsta koja dobro podnosi oštru planinsku i kontinentalnu klimu, a raste na karbonatnim, ali i silikatnim, pjeskovitim, šljunčanim, slabo humusnim tlima. Najčešća je u planinskim dolinama, uz potoke i rijeke gdje na vlažnim staništima s mnogim higrofilnim vrstama tvori guste, fragmentarne sastojine, često obnovljene iz panja. To su povremeno plavljeni lokaliteti uz vodotoke ili okolne padine kroz koje se postrano procjeđuje voda prema glavnom vodotoku.
Prisutnost bijele johe u Hrvatskoj vezana je za dva oštro ograničena područja: sjeveroistočno, u nizinskom dijelu Hrvatske (110 do 190 m nadmorske visine) uz tok rijeke Drave; i sjeverozapadno, u alpsko-dinarskom području uz rijeku Kupu i njene pritoke (220 do 550 m nadmorske visine).
U sjeveroistočnom području rasprostiru se sastojine definirane kao Equiseto hyemali-Alnetum incanae Moor 1958 (Trinajstić 1964; Franjić et al. 1999), u kojima se u odnosu na sastojine bijele johe u gorskome području kao razlikovne vrste izdvajaju: Ulmus laevis Pall., U. minor Mill., Prunus padus L., Equisetum hyemale L. i druge. Ove sastojine danas su rjeđe nego nekada, manje se koriste, a zbog nekadašnjeg načina sječe vrlo su često panjače. Nekadašnje šumske površine koje je Trinajstić istraživao (1964) danas više nisu pod šumom. Poplavne i vlažne šume u blizini urbanih središta isušuju se, a na njihovim su staništima izgrađena naselja i komunalna infrastruktura tako da su danas šume sive johe uz tok rijeke Drave u velikoj regresiji.
U alpsko-dinarskom području istraživane sastojine svrstane su u šumsku zajednicu Lamio orvalae-Alnetum incanae Dakskobler 2010, a opisane su kao geografska varijanta Helleborus dumetorum (Vukelić et al. 2012, 2017). Ona se ističe mnogim vrstama iz dinarske vegetacijske zone, poput Fraxinus excelsior L., Ulmus glabra Huds., Fagus sylvatica L., Acer pseudoplatanus L., Lamium orvala L., Lunaria rediviva L., Scopolia carniolica Jacq., Helleborus dumetorum Waldst. et Kit., Knautia drymeia Heuff. subsp. drymeia i druge. Siva joha ovdje raste u kanjonu Kupe i na obalnim terasama koje su uglavnom antropogenizirane, a nekoć su bile korištene u poljoprivrednoj proizvodnji.
Glavni ciljevi rada bili su utvrditi morfološku raznolikost i strukturiranost populacija bijele johe u Hrvatskoj te istražiti na koji način ekološki i geografski čimbenici utječu na varijabilnost populacija.
Materijal za morfometrijsku analizu sakupljen je iz sedam prirodnih populacija – tri iz kontinentalne i četiri iz alpsko-dinarske regije (slika 1). Svaka populacija predstavljena je s 20 stabala, a svako stablo s 20 zdravih i neoštećenih listova s kratkih fertilnih izbojaka iz vanjskog, osvijetljenog dijela krošnje. Listovi su sakupljeni početkom srpnja 2016. godine, odnosno sredinom vegetacijskog razdoblja, kada su s obzirom na dimenzije i oblik u potpunosti razvijeni. Listovi su skenirani i izmjereni pomoću programa WinFolia. Ukupno je izmjereno devet lisnih značajki (slika 2).
Deskriptivnim statističkim metodama utvrđena je visoka varijabilnost listova bijele johe, kao i jasan trend divergencije između kontinentalnih i alpsko-dinarskih populacija (slika 3). Najvarijabilnijim su se pokazale značajke: površina plojke, širina lista mjerena na 90 % dužine plojke i dužina peteljke. Za razliku od toga, najmanji stupanj varijabilnosti zabilježen je za varijable koje opisuju osnovu lisne plojke. Iako je trend variranja pojedinih značajki po populacijama bio sličan, u populacijama alpsko-dinarske regije zabilježena je nešto viša