Next Article in Journal
Antifungal Activity of Cedrol from Cunninghamia lanceolate var. konishii against Phellinus noxius and Its Mechanism
Previous Article in Journal
Nematocidal Potential of Phenolic Acids: A Phytochemical Seed-Coating Approach to Soybean Cyst Nematode Management
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Plants
Volume 13
Issue 2
10.3390/plants13020320
plants-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Patterns of Leaf and Fruit Morphological Variation in Marginal Populations of Acer tataricum L. subsp. tataricum

by
Igor Poljak
1,
Antonio Vidaković
1,
Luka Benić
1,
Katarina Tumpa
1,
Marilena Idžojtić
1 and
Zlatko Šatović
2,3,*
1
Institute of Forest Genetics, Dendrology and Botany, Faculty of Forestry and Wood Technology, University of Zagreb, Svetošimunska cesta 23, HR-10000 Zagreb, Croatia
2
Department for Seed Science and Technology, Faculty of Agriculture, University of Zagreb, Svetošimunska cesta 25, HR-10000 Zagreb, Croatia
3
Centre of Excellence for Biodiversity and Molecular Plant Breeding, Svetošimunska cesta 25, HR-10000 Zagreb, Croatia
*
Author to whom correspondence should be addressed.
Plants 2024, 13(2), 320; https://doi.org/10.3390/plants13020320
Submission received: 4 December 2023 / Revised: 15 January 2024 / Accepted: 18 January 2024 / Published: 21 January 2024
(This article belongs to the Section Plant Development and Morphogenesis)
Browse Figures
Versions Notes

Abstract

:
Marginal populations are usually smaller and more isolated and grow in less favourable conditions than those at the distribution centre. The variability of these populations is of high importance, as it can support the adaptations needed for the conditions that they grow in. In this research, the morphological variability of eight Tatar maple (Acer tataricum L. subsp. tataricum) populations was analysed. Tatar maple is an insect-pollinated and wind-dispersed shrub/tree, whose northwestern distribution edge is in southeastern Europe. Morphometric methods were used to analyse the variability of the populations using leaf and fruit morphology. The research revealed significant differences between and within populations. Furthermore, differences in the distribution of the total variability were noted, which suggest that different evolutionarily factors affect different plant traits. Correlation analysis confirmed a weak dependency between the vegetative and generative traits. In addition, no evidence was found for the presence of isolation by environment (IBE). However, the Mantel test for isolation by distance (IBD) was significant for the leaf morphometric traits and non-significant for the fruit morphometric traits. Being the marginal leading-edge populations, they are younger and were less likely to have had time for adaptation to local environments, which would have resulted in the development of IBE. Overall, edge populations of Tatar maple were characterised by great morphological variability, which helps these populations in their response to the intensive selective pressures they face in their environment.

1. Introduction

A certain degree of morphological variability of vegetative and generative organs is characteristic of plant species [1]. This variability is, to some extent, conditioned by genetic variability [2] but also by environmental heterogeneity [1,3]. Many authors have determined the influence of various environmental conditions on leaf morphology, including temperature [4,5], precipitation [6,7] and insolation gradients [8,9]. Other studies confirm that not only are the leaves prone to environmentally induced morphological changes, but the fruits are as well [10,11]. However, most of the research on morphological variability to date has been carried out on either leaves or fruits, rarely both [12,13].
In wind-pollinated species, unhindered gene flow usually occurs, which has a homogenizing effect on a species’ genetic and morphological variability [14]. Therefore, differentiation between populations is expected to be low, as confirmed for Acer campestre L. [15], A. monspessulanum L. [16] and Pinus sylvestris L. [14]. This, however, does not imply a low diversity within populations; on the contrary, it is usually significantly higher than between populations [17]. Variability, both genetic and morphological, is crucial for adaptation, biodiversity and speciation [18] and is under the influence of several evolutionary processes, including gene flow, genetic drift and natural selection. There are two well-known patterns that explain divergence among populations: isolation by distance (IBD) and isolation by environment (IBE). The isolation-by-distance model indicates that differentiation between populations increases with geographical distance [19], arising from limited gene flow and the presence of genetic drift. The original model of isolation by distance presented by Wright [20] showed how population structure is generated by limits on parent-offspring dispersal distances and that this process is expected to be more pronounced as population size decreases. On the other hand, isolation by environment (IBE) explains genetic and morphological differentiation through environmental differences between populations, where genetic/morphological and environmental distances are positively correlated, independent of geographical distance [19,21]. However, these patterns are not exclusive, as confirmed in the study by Sexton et al. [22], indicating that differences among populations can follow both IBD and IBE patterns.
Populations that are particularly affected by limiting gene flow and other selective pressures are marginal populations [23,24]. Therefore, marginal populations are likely to have increased genetic differentiation [25,26]. As a result, they are characterised by the ability to adapt to suboptimal environments, which makes them important from the evolutionary point of view [23]. Furthermore, this adaptability marks them as highly important for the long-term conservation of genetic diversity of species [27]. With climate change in mind, marginal populations will represent potential gene pools for the expansion and shift of the species away from the drought in the south towards the newly formed, favorable habitats to the north [28]. They are usually characterised by greater demographic stochasticity and differentiation whilst simultaneously being the youngest and with the lowest regional diversity [25,27].
Although research into the plant species marginal populations has been conducted from a genetics point of view [28,29,30,31], such research from the morphological point of view is sparse [32]. In this study, therefore, we aimed at using the morphological characteristics to gain insight into the diversity of marginal populations. As a model species, the Tatar maple (Acer tataricum L. subsp. tataricum, Sapindaceae) was selected, which reaches its northern and northwestern ranges’ limits in southeastern Europe. It is a slow-growing, small tree that reaches a height of 5 to 7 m, often with multiple or forked trunks. Flowers are scented and white and pollinated by insects. The fruits are samaras, winged nutlets that ripen in late summer or early autumn. The species is predominantly found in the lowlands or the colline zone, with an overall range between 300 and 1700 m a.s.l. Although native to warm and dry climates, it is successfully cultivated across Europe and Asia, although the trees are susceptible to frost [33].
Leaf and fruit morphological traits of eight populations of Tatar maple were analysed, with the aim of investigating the following: (1) morphological intra- and interpopulation diversity and (2) the level of phenotypical variability in the context of isolation by distance (IBD) and isolation by environment (IBE).

2. Materials and Methods

2.1. Plant Material

In the research, eight populations of Tatar maple were sampled (Figure 1, Table S1). Samples were collected only from fully developed individuals 5 to 7 m tall and about 30 years old. All analysed populations grew in homogeneous pedological conditions. According to Bartha [34], Tatar maple fruits ripen in late August and early September, which is why we conducted sample collection in late September, to make sure that the samaras were fully ripened and that the leaf fall had not started. Leaf and fruit samples were collected simultaneously, 20 leaves and samaras per individual plant, with a total of 10 to 15 individuals per population. For sampling, well-lit branches were selected. The plant material was stored in airtight bags in the field and herbarised the same day upon arrival in the laboratory, following the protocol of plant collecting and documentation by Simpson [35]. Samples were herbarised between newspaper sheets in the herbarium of the Faculty of Forestry and Wood Technology in Zagreb (Herbarium DEND).

2.2. Leaf and Fruit Morphometric Analysis

After being herbarised and left to fully dry, the leaves were scanned with a MICROTEK ScanMaker 4800 (Microtek International, Inc., Hsinchu, Taiwan). The created scans were used in the WinFolia [36] morphometric analysis program by using the in-program parameters as follows: leaf area (LA); leaf length (LL); maximum leaf width (MLW); leaf length, measured from the leaf base to the point of maximum leaf width (PMLW); leaf blade width at 50% (LW1) and 90% of leaf blade length (LW2); and petiole length (PL). In addition, three leaf shape measurements were conducted, describing the leaf shape by angles enclosed by the main leaf vein (the centre of the leaf blade) and the line connecting the leaf blade base to a set point on the leaf margin, at 10% (LA1) and 25% (LA2) of the total leaf blade length, as well as the shape of the leaf blade itself (FC) (Figure 2).
Maple fruits, i.e., the mericarps, were also scanned and measured using the WinFolia program [36], with the following traits measured: mericarp area (MA); mericarp length (ML); maximum mericarp width (MMW); length of the mericarp, measured from the base to the point of maximum width (PMMW); width of the mericarp at 90% of the mericarp’s length (MW90); nut length (NL); and nut width (NW). For fruits, an interactive measurement was used for the angle enclosed by the wings (WA) (Figure 2).

2.3. Statistical Analysis

Descriptive statistics were calculated for each individual trait and for each population, with the goal of revealing the overall range of their variability [37]. In addition, Pearson correlation coefficients were calculated between all leaf and fruit traits. Analysis of variance was used to determine the variability between the studied populations, as well as between trees/shrubs within the populations. The “tree/shrub” factor was nested within the “population” factor. Descriptive statistics and analysis of variance were carried out using the STATISTICA software package Version 13 [38].
The Mantel test was used to evaluate the correlations between the morphometric, geographic, and environmental data [37]. In this study, four dissimilarity matrices were calculated in order to describe the differences between the analysed populations: (1) leaf morphometric differences as squared Mahalanobis distances between the pairs of populations; (2) fruit morphometric differences as squared Mahalanobis distances between the pairs of populations; (3) environmental distances as the Euclidian distances between the population means for the first three PCs of the principal component analysis; and (4) geographic distance from the latitude and longitude of the sampling site. The significance level was assessed after 10,000 permutations as implemented in NTSYS-pc Ver. 2.21L [39].

3. Results

3.1. Correlations

A total of 93 statistically significant correlations were observed between the tested pairs of leaf and fruit morphological traits, none of which were of a negative character. The strength of a relationship based on the r value was classified according to Sokal and Rohlf [37] as follows: none or very weak (r < 0.3); weak (0.3 < r < 0.5); moderate (0.5 < r < 0.7); and strong (r > 0.7). Out of the total 93 correlations, only 22 demonstrated r values larger than 0.7.
When leaf characteristics are considered (Table S2), 39 significant correlations were detected, among which 15 were strong. The traits that were correlated to all of the other characters were LA, MLW, LW1, LW2 and PL. However, the largest number of strong correlations was observed for LA (6) and MLW (5), while correlations were the weakest in traits related to leaf shape: FC, LA1 and LA2. A strong correlation among leaf shape traits was detected only between LA1 and LA2. No significant correlation was detected between LL and FC or LA1 and LA2 or between PMLW and the same three traits.
Observing the fruit morphological traits (Table S3), significant correlations were noted between all of the measured traits except WA. In total, seven strong correlations were detected, four of which refer to MA (ML, MMW, PMMW, MW90). In general, nut characteristics were less significantly correlated with other fruit traits.
Analysis between leaf and fruit traits showed 33 significant correlations (Table S4), all of which were classified as very weak or weak. No significant correlations were detected between MMW, NW and WA and any leaf traits. On the other hand, traits ML, PMMW and NL were significantly correlated to all of the leaf traits except LA1 and LA2, which showed no significant correlations with any of the fruit traits.

3.2. Leaf Phenotypic Traits

The results of the performed statistical analysis for leaf morphometric traits are shown in Table 1, by population and for all populations together. The average leaf length (LL) was 6.29 cm, whereas the average maximum blade width was 4.08 cm. Mean petiole length (PL) was 4.30 cm, and the overall mean leaf area, for all populations, was 18.38 cm2. The highest mean coefficient of variation was noted for the trait LW2 (37.40%), with the second highest noted for LA (36.45%). The lowest coefficients of variations were noted for both leaf angles measured, 4.44% and 6.50%, for LA1 and LA2, respectively.
The most variable population was P3, with the maximum values of coefficient of variation for seven out of ten measured traits (LA, LL, MLW, PMLW, LW1, PL, LA2). Furthermore, population P2 was characterised by the overall largest leaves, having the highest mean values of six leaf traits (LA, LL, MLW, LW1, LW2, PL). On the other hand, population P8 was characterised by the smallest leaves (LA, LL, MLW, PMLW, LW1, LW2, PL), as well as the narrowest leaf angles (LA1, LA2) and the most elongated leaf blade shape (FC). In addition, P8 was characterised by the most homogeneous values, i.e., the least variable leaf traits, with five traits having the minimum coefficient of variation values (LA, LL, MLW, LW1, LW2).

3.3. Fruit Phenotypic Traits

The results of the morphometric analysis of fruits have been presented on the overall sample, as well as on the individual population’s level (Table 2). The average mericarp area (MA), i.e., the area of the winged nutlet, on the overall level was 2.53 cm2, whereas the overall mean length (ML) was 3.18 cm. The average nut length (NL) came up to 1.23 cm, with a mean width (NW) of 0.63 cm. The angle enclosed by the two samaras, i.e., the two mericarps, had a mean value of 66.35°. The most variable trait was the single mericarp area (MA), with a coefficient of variation of 25.96%. The remaining fruit traits were characterised by prominently lower coefficients of variation, which ranged between 12.38% for mericarp length (ML) up to 19.88% for the angle enclosed by the two mericarps (WA).
Population P8 was characterised by having the lowest mean values for six of the measured traits: MA, ML, MMW, PMMW, MW90 and NL. In contrast, population P1 was defined by having the longest mericarp length (ML) and largest mericarp area (MA). In addition, maximum values were also found in population P6, for PMMW, NW and WA. The least variable population was P1, with lowest coefficient of variation values for four out of eight measured traits: MA, MMW, NW and WA. As the most variable population, P7 stands out, with four traits having the highest values of coefficient of variation (MA, ML, PMMW, NW).

3.4. Analysis of Variance—ANOVA

To determine the significance of differentiation by each individual trait, ANOVA analysis was performed. Both populations and individual maple shrubs/trees differed significantly when analysed leaf traits are considered (Table 3). For all analysed traits, between 50 and 75% of total variance could be ascribed to variance between leaves of the same individual, found as the error component. The smallest component of the variance could be ascribed to the interpopulation diversity, only 4.58–16.28% of total variance, with the overall lowest value of 4.58% found in LA1. The intrapopulation diversity had intermediate values, with a range of 16.49–36.64%.
Results of the performed ANOVA analysis for measured fruit traits can be found in detail in Table 4. The results of ANOVA analysis for fruit traits of Tatar maple showed significant differences between populations of Tatar maple, as well as between trees/shrubs of the same population, for all of the eight measured traits. The majority of measured traits had the largest proportion of variability assigned to the intrapopulation variability, around 50%. The proportions of variability assigned to the interpopulation component and the error component were similar, with only two traits, NW and PMNW, demonstrating notably larger proportions assigned to the error component.

3.5. Mantel Test

The results of the Mantel test are shown in Figure 3. Statistical significance was observed only between geographical and leaf morphological distances, which indicates that morphologically similar populations were also geographically closer. Linear regression of those distances was used to obtain a regression line, and the coefficient of determination accounted for R2 = 0.3421. Thus, 34.21% of leaf morphological variability among analysed populations could be explained by their geographic distance. On the other hand, a significant relationship between environmental distances between populations and their leaf (r = 0.407, p = 0.083) and fruit (r = 0.361, p = 0.140) morphological distances was not recorded. In addition, no statistically significant correlations between fruit morphological distances and geographical distances (r = 0.451, p = 0.098), nor between leaf and fruit morphological distances (r = 0.319, p = 0.149), were disclosed.

4. Discussion

Tatar maple is a poorly researched species, with no discernible morphological studies conducted so far. Even so, the data on leaf dimensions can be found in several written publications, with leaf size being reported to be between 6 and 10 cm long and 3 to 8 cm wide, with a petiole of about 2.5 cm [34,40,41]. These values are in accordance with those revealed in this study, albeit the petiole length of the populations presented here is much longer, on average 4.30 cm. This deviation in petiole length might be due to adjustments in phyllotaxy, directed towards better light capture efficiency [42], and due to the fact that plants tend to display longer petioles on more productive habitats [43]. Adjustments in the petiole length can also be stimulated for better nutrient transport efficiency and leaf support and protection [44]. Such adaptations in petiole length are enabled through differences in the total number of cells, individual cell size, or both [45]. In addition, the average mericarp length (ML) in this study was 3.18 cm, and it is slightly longer than the values reported by Bartha [34] and Idžojtić [46], with reported mericarp values of 2 to 3 cm.
Positive correlations between leaf/fruit size morphological traits, like leaf/fruit area, length and width, were noted in our research. Such a relation between named traits is expected and has been previously reported for woody plant species within Populus L. [47], Theobroma L. [13] and Alnus L. [48] species, as well as on a general level [49]. According to the latter, the length and width of leaves/fruits directly influence within-species variation of total leaf/fruit area with high accuracy. On the other hand, the relationship between leaf and fruit morphological traits is less studied. Fishler et al. [50] indicate that the leaf area is one of the main limiting factors in fruit growth, as greater photosynthetically active leaf area better supplies fruit growth, therefore determining its size. In our research, there was a very weak or weak correlation between leaf and fruit traits, few of which were statistically significant. Similar results were reported for Eucommia ulmoides Oliv. [51] and Malania oleifera Chun and S.K. Lee [52].
With regard to the coefficients of variability, they ranged from 4.44% to 36.45% for leaf traits. The high levels of variability found for leaf traits aren’t surprising, considering the importance of morphological adaptability of leaves for the survival of plants, since leaves are the key players in the interaction with the environment [53]. Leaves, as organs sensitive to changes in the environment, exhibit phenotypic variability as a response to stressors, which in turn can alter the primary morphogenesis [9]. Previous studies have even confirmed a global trend of leaf trait variation across individual species with broad and diverse biogeographical and environmental affinities [54,55]. This is true for numerous plant species, in accordance with various environmental factors, such as light availability [56,57,58], altitude [59,60,61] and relative humidity [62,63], and attests to the importance of leaf adaptations for survival and growth [53]. Unlike leaves, fruits had significantly lower coefficients of variation, within the range of 12.38% to 25.96%. The lower variability of fruit traits revealed in this research is not uncommon. Leaf traits have been continuously found to be more variable than those of fruits, for numerous species, such as Pterocarpus erinaceus Poir. [64], Staphylea japonica (Thunb.) Mabb. [65] and Malus zumi (Matsum.) Rehder [66]. This fact could be attributed to the pronounced influence of the environmental factors on the leaves as main assimilation organs [67], which favours greater variability.
The high levels of variability in leaf shape and size, as well as in the fruits of Tatar maple, were found on both levels, intra- and interpopulation. In both cases, intra-population variability was greater than inter-population variability. This result is generally consistent with that predicted by population genetics theory when the plant species are taken into consideration. In general, it is well known that woody species maintain more variation within populations than other types of species, such as annual and perennial herbaceous species, but have less variation among populations [17]. Woody species with large geographic ranges, outcrossing breeding systems and wind or animal-ingested seed dispersal have more genetic diversity within species and populations but less variation among populations than woody species with other combinations of traits. In addition, long-lived plant species with large ranges and high fecundities typically have large, stable populations. Such populations are resistant to chance fluctuations in gene or genotype frequencies and should maintain more variation than populations that experience large fluctuations in size.
When the results of the hierarchical variance analysis are considered, some differences in the distribution of the total variability can be noted when leaf and fruit variabilities are compared. For leaves, the greatest component of the total variability is assigned to the leaves’ variability within each tree, and, in the case of fruits, it is assigned to the variability of the individuals within populations. Furthermore, interpopulation variability was somewhat greater when measured fruit traits were compared to those measured on leaves. These different variance partitioning patterns suggest to us that different factors affect reproductive and vegetative traits. As already mentioned, reproductive traits can be considered more evolutionarily conserved than vegetative traits [68], and vegetative traits can be considered evolutionarily more “plastic” [69]. In general, leaf traits’ components of variability indicated greater influence of the environmental conditions and leaf placement on the shoots, as the error component was by far the most significant in the overall variability. In other words, leaves are the organs that respond to the heterogeneity of the environment stronger than fruits, possibly due to their significance in the overall adaptability and survival of the plants [70]. Fruit morphology, on the other hand, has been shown to be under stronger genetic influence, i.e., adaptive evolution, in numerous species, and as such is less likely to demonstrate high levels of variability [71,72,73]. Similar patterns of variability in leaf and fruit traits were found in other species as well: Castanea sativa Mill. [74,75,76,77], Ulmus laevis Pall [78] and Sorbus domestica L. [79].
In addition to lower levels of morphological variability, fruits analysed in this research showed no significant correlation to either isolation by distance (IBD) or isolation by environment (IBE). On the other hand, leaf morphology was found to have a significant correlation to IBD, with more distant populations demonstrating greater morphological differences. In marginal populations, both the specific environments and geographical distance can affect morphological variability, as those populations are often fragmented geographically and, due to varying environments, develop habitat-specific adaptations [80,81]. Such morphological differentiation of marginal populations has previously been noted in Juniperus thurifera L. [82] in northern Africa and Pinus sylvestris L. [83] in the western Mediterranean, with both species demonstrating significant morphological differences in their marginal populations, much like the populations analysed here. Being the marginal populations, they are younger [28] and were less likely to have had time for pronounced adaptation to local environments, which would have resulted in development of IBE. Additionally, though geographically isolated, they are found in relatively homogeneous pedological and climatological conditions, which could hinder ecological diversification. This lack of selection signature in leading-edge populations has been previously noted in younger, leading-edge populations of Biscutella laevigata L., which demonstrated a smaller degree of local adaptation than the trailing-edge populations of the same species, thus inhibiting the IBE [84].

5. Conclusions

Southeastern Europe represents the northwestern edge of Tatar maple natural distribution. In such marginal populations, both the specific environments and geographical distance can affect morphological variability, as those populations are often geographically and ecologically fragmented. Indeed, in this study, high levels of variability were found for leaf traits, attributed to their morphological adaptability and variability conditioned by environmental heterogeneity. Unlike leaves, fruits had significantly lower coefficients of variation. The high levels of morphological variability of Tatar maple were found at both intra- and interpopulation levels. In both cases, intra-population variability was greater than inter-population variability. Weak but significant correlations were found between some of the leaf and fruit characters. However, according to the Mantel test, the only significant correlation was found between geographical and leaf morphological distances , i.e., isolation-by-distance pattern. Overall, studied edge populations of Tatar maple are characterised by great morphological variability, which represents a coping mechanism that helps these populations in their response to the intensive selective pressures. In the future, such pressures will be even more frequent due to climate change, so northward expansion of the Tatar maple natural range could be expected.

Supplementary Materials

The following supporting information can be downloaded at https://www.mdpi.com/article/10.3390/plants13020320/s1: Table S1: Geographic coordinates, altitudes, bioclimatic variables and number of samples for eight studied Acer tataricum L. subsp. tataricum populations. Bioclimatic variables: BIO1 (Annual Mean Temperature); BIO2 (Mean Diurnal Range (Mean of monthly (max temp—min temp)); BIO3 (Isothermality (BIO2/BIO7) (×100)); BIO4 (Temperature Seasonality (standard deviation ×100)); BIO5 (Max Temperature of Warmest Month); BIO6 (Min Temperature of Coldest Month); BIO7 (Temperature Annual Range (BIO5-BIO6)); BIO8 (Mean Temperature of Wettest Quarter); BIO9 (Mean Temperature of Driest Quarter); BIO10 (Mean Temperature of Warmest Quarter); BIO11 (Mean Temperature of Coldest Quarter); BIO12 (Annual Precipitation); BIO13 (Precipitation of Wettest Month); BIO14 (Precipitation of Driest Month); BIO15 (Precipitation Seasonality (Coefficient of Variation)); BIO16 (Precipitation of Wettest Quarter); BIO17 (Precipitation of Driest Quarter); BIO18 (Precipitation of Warmest Quarter); BIO19 (Precipitation of Coldest Quarter). Populations: P1–Odransko polje; P2–Lipovljani; P3–Veliki Grđevac; P4–Mali Grđevac; P5–Grubišno Polje; P6–Virovitica; P7–Požega; P8–Županja; Table S2: Pearson’s correlation coefficients between analysed leaves’ traits. Leaves’ morphometric traits acronyms: leaf area (LA); leaf length (LL); maximum leaf width (MLW); leaf length, measured from the leaf base to the point of maximum leaf width (PMLW); leaf blade width at 50% (LW1) and 90% of leaf blade length (LW2); petiole length (PL); form coefficient (FC); angles enclosed by the main leaf vein (the centre of the leaf blade) and the line connecting the leaf blade base to a set point on the leaf margin, at 10% (LA10) and 25% (LA25). n.s. not significant at p > 0.05; * significant at 0.01 < p < 0.05; ** significant at 0.001 < p < 0.01; *** significant at p < 0.001; Table S3: Pearson’s correlation coefficients between analysed fruit traits. Fruit morphometric traits acronyms: mericarp area (MA); mericarp length (ML); maximum mericarp width (MMW); length of the mericarp, measured from the basis to the point of maximum width (PMMW); width of mericarp at 90% of mericarp’s length (MW90); nut length (NL); nut width (NW); angle enclosed by the wings (WA). n.s. not significant at p > 0.05; * significant at 0.01 < p < 0.05; ** significant at 0.001 < p < 0.01; *** significant at p < 0.001; Table S4: Pearson’s correlation coefficients between analysed leaves and fruit traits. Leaves’ morphometric traits acronyms: leaf area (LA); leaf length (LL); maximum leaf width (MLW); leaf length, measured from the leaf base to the point of maximum leaf width (PMLW); leaf blade width at 50% (LW1) and 90% of leaf blade length (LW2); petiole length (PL); form coefficient (FC); angles enclosed by the main leaf vein (the centre of the leaf blade) and the line connecting the leaf blade base to a set point on the leaf margin, at 10% (LA10) and 25% (LA25). Fruit morphometric traits acronyms: mericarp area (MA); mericarp length (ML); maximum mericarp width (MMW); length of the mericarp, measured from the basis to the point of maximum width (PMMW); width of the mericarp at 90% of the mericarp’s length (MW90); nut length (NL); nut width (NW); angle enclosed by the wings (WA). n.s. not significant at p > 0.05; * significant at 0.01 < p < 0.05; ** significant at 0.001 < p < 0.01; *** significant at p < 0.001.

Author Contributions

Conceptualisation, I.P. and Z.Š.; methodology, I.P., M.I. and Z.Š.; software, I.P. and Z.Š.; validation, I.P., M.I. and Z.Š.; formal analysis, I.P. and Z.Š.; investigation, A.V. and L.B.; resources, I.P.; data curation, I.P. and Z.Š.; writing—original draft preparation, I.P., A.V. and K.T.; writing—review and editing, M.I. and Z.Š.; visualisation, I.P., A.V. and Z.Š.; supervision, I.P., M.I. and Z.Š. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the financial support of the University of Zagreb.

Data Availability Statement

Data are contained within the article or Supplementary Material.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Matesanz, S.; Gianoli, E.; Valladares, F. Global change and the evolution of phenotypic plasticity in plants. Ann. N. Y. Acad. Sci. 2010, 1206, 35–55. [Google Scholar] [CrossRef] [PubMed]
  2. Khadivi-Khub, A.; Karimi, S.; Kameli, M. Morphological diversity of naturally grown Crataegus monogyna (Rosaceae, Maloideae) in Central Iran. Braz. J. Bot. 2015, 38, 921–936. [Google Scholar] [CrossRef]
  3. Vázquez, D.P.; Gianoli, E.; Morris, W.F.; Bozinovic, F. Ecological and evolutionary impacts of changing climatic variability. Biol. Rev. 2015, 92, 22–42. [Google Scholar] [CrossRef] [PubMed]
  4. Royer, D.L.; Meyerson, L.A.; Robertson, K.M.; Adams, J.M. Phenotypic Plasticity of Leaf Shape along a Temperature Gradient in Acer rubrum. PLoS ONE 2009, 4, e7653. [Google Scholar] [CrossRef]
  5. McKee, M.L.; Royer, D.L.; Poulos, H.M. Experimental evidence for species-dependent responses in leaf shape to temperature: Implications for paleoclimate inference. PLoS ONE 2019, 14, e0218884. [Google Scholar] [CrossRef]
  6. Leuschner, C.; Voß, S.; Foetzki, A.; Clases, Y. Variation in leaf area index and stand leaf mass of European beech across gradients of soil acidity and precipitation. Plant Ecol. 2006, 186, 247–258. [Google Scholar] [CrossRef]
  7. Souza, M.L.; Duarte, A.A.; Lovato, M.B.; Fagundes, M.; Valladares, F.; Lemos-Filho, J.P. Climatic factors shaping intraspecific leaf trait variation of a neotropical tree along a rainfall gradient. PLoS ONE 2018, 13, e0208512. [Google Scholar] [CrossRef]
  8. Ackerly, D.; Knight, C.; Weiss, S.; Barton, K.; Starmer, K. Leaf size, specific leaf area and microhabitat distribution of chaparral woody plants: Contrasting patterns in species level and community level analyses. Oecologia 2002, 130, 449–457. [Google Scholar] [CrossRef]
  9. Xu, F.; Guo, W.; Xu, W.; Wei, Y.; Wang, R. Leaf morphology correlates with water and light availability: What consequences for simple and compound leaves? Prog. Nat. Sci. 2009, 19, 1789–1798. [Google Scholar] [CrossRef]
  10. Soleimani, A.; Etemad, V.; Calagari, M.; Namiranian, M.; Shirvani, A. Influence of climatic factors on fruit morphological traits in Populus euphratica Oliv. Ann. For. Res. 2014, 57, 31–38. [Google Scholar] [CrossRef]
  11. Gao, S.; Wang, B.; Liu, F.; Zhao, J.; Yuan, J.; Xiao, S.; Masabni, J.; Zou, F.; Yuan, D. Variation in Fruit Morphology and Seed Oil Fatty Acid Composition of Camellia oleifera Collected from Diverse Regions in Southern China. Horticulturae 2022, 8, 818. [Google Scholar] [CrossRef]
  12. Ertan, E. Variability in leaf and fruit morphology and in fruit composition of chestnuts (Castanea sativa Mill.) in the Nazilli region of Turkey. Genet. Resour. Crop Evol. 2007, 54, 691–699. [Google Scholar] [CrossRef]
  13. Santos, R.C.; Pires, J.L.; Correa, R.X. Morphological characterization of leaf, flower, fruit and seed traits among Brazilian Theobroma L. species. Genet. Resour. Crop Evol. 2012, 59, 327–345. [Google Scholar] [CrossRef]
  14. Wójkiewicz, B.; Litkowiec, M.; Wachowiak, W. Contrasting patterns of genetic variation in core and peripheral populations of highly outcrossing and wind pollinated forest tree species. AoB Plants 2016, 8, plw054. [Google Scholar] [CrossRef] [PubMed]
  15. Kvesić, S.; Hodžić, M.M.; Ballian, D.; Gömöry, D.; Fussi, B. Genetic variation of a widespread subdominant tree species (Acer campestre L.) in Bosnia and Herzegovina. Tree Genet. Genomes 2020, 16, 82. [Google Scholar] [CrossRef]
  16. Motahari, B.; Shabanian, N.; Rahmani, M.-S.; Mohammad-Hasani, F. Genetic diversity and genetic structure of Acer monspessulanum L. across Zagros forests of Iran using molecular markers. Gene 2021, 769, 145245. [Google Scholar] [CrossRef]
  17. Hamrick, J.L.; Godt, M.J.W.; Sherman-Broyles, S.L. Factors influencing levels of genetic diversity in woody plant species. In Population Genetics of Forest Trees; Springer: Berlin/Heidelberg, Germany, 1992; pp. 95–124. [Google Scholar]
  18. Ramel, C. Biodiversity and intraspecific genetic variation. Pure Appl. Chem. 1998, 70, 2079–2084. [Google Scholar] [CrossRef]
  19. Morente-López, J.; García, C.; Lara-Romero, C.; García-Fernández, A.; Draper, D.; Iriondo, J.M. Geography and Environment Shape Landscape Genetics of Mediterranean Alpine Species Silene ciliata Poiret. (Caryophyllaceae). Front. Plant Sci. 2018, 9, 1698. [Google Scholar] [CrossRef]
  20. Wright, S. Isolation by distance. Genetics 1943, 28, 114–138. [Google Scholar] [CrossRef]
  21. Wang, I.J.; Bradburd, G.S. Isolation by environment. Mol. Ecol. 2014, 23, 5649–5662. [Google Scholar] [CrossRef]
  22. Sexton, J.P.; Hangartner, S.B.; Hoffmann, A. Genetic isolation by environment or distance: Which pattern of gene flow is most common? Evolution 2013, 68, 1–15. [Google Scholar] [CrossRef] [PubMed]
  23. Kawecki, T.J. Adaptation to Marginal Habitats. Annu. Rev. Ecol. Evol. Syst. 2008, 39, 321–342. [Google Scholar] [CrossRef]
  24. Sexton, J.P.; McIntyre, P.J.; Angert, A.L.; Rice, K.J. Evolution and Ecology of Species Range Limits. Annu. Rev. Ecol. Evol. Syst. 2009, 40, 415–436. [Google Scholar] [CrossRef]
  25. Eckert, C.G.; Samis, K.E.; Lougheed, S.C. Genetic variation across species’ geographical ranges: The central–marginal hypothesis and beyond. Mol. Ecol. 2008, 17, 1170–1188. [Google Scholar] [CrossRef] [PubMed]
  26. Langin, K.M.; Sillett, T.S.; Funk, W.C.; Morrison, S.A.; Ghalambor, C.K. Partial support for the central–marginal hypothesis within a population: Reduced genetic diversity but not increased differentiation at the range edge of an island endemic bird. Heredity 2017, 119, 8–15. [Google Scholar] [CrossRef] [PubMed]
  27. Hampe, A.; Petit, R.J. Conserving biodiversity under climate change: The rear edge matters. Ecol. Lett. 2005, 8, 461–467. [Google Scholar] [CrossRef]
  28. Vakkari, P.; Rusanen, M.; Heikkinen, J.; Huotari, T.; Kärkkäinen, K. Patterns of genetic variation in leading-edge populations of Quercus robur: Genetic patchiness due to family clusters. Tree Genet. Genomes 2020, 16, 73. [Google Scholar] [CrossRef]
  29. Morente-Lopez, J.; Lara-Romero, C.; Garica-Fernandez, A.; Rubio Teso, M.T.; Prieto-Benitez, S.; Iriondo, J.M. Gene flow effects on populations inhabiting marginal areas: Origin matters. J. Ecol. 2020, 109, 139–153. [Google Scholar] [CrossRef]
  30. Paica, I.C.; Banciu, C.; Maria, G.M.; Vladimirescu, M.; Manole, A. Genetic diversity in marginal populations of Nitraria schoberi L. from Romania. Forests 2022, 14, 882. [Google Scholar] [CrossRef]
  31. Kitamura, K.; Uchiyama, K.; Ueno, S.; Ishizuka, W.; Tsuyama, I.; Goto, S. Geographical Gradients of Genetic Diversity and Differentiation among the Southernmost Marginal Populations of Abies sachalinensis Revealed by EST-SSR Polymorphism. Forests 2020, 11, 233. [Google Scholar] [CrossRef]
  32. Farris, M.A.; Schaal, B.A. Morphological and genetic variation in ecologically central and marginal populations of Rumex acetosella L. (Polygonaceae). Am. J. Bot. 1983, 70, 246–255. [Google Scholar] [CrossRef]
  33. Bartha, D. Enzyklopadie der Holzgewachse; University of Sopron: Sopron, Hungary, 1997. [Google Scholar]
  34. Bartha, D. Acer tataricum L. In Enzyklopädie der Holzgewächse: Handbuch und Atlas der Dendrologie; Roloff, A., Weisberger, H., Lang, U.M., Stimm, B., Schütt., P., Eds.; Wiley-VCH: Weinheim, Germany, 2011; p. 11. [Google Scholar]
  35. Simpson, M.G. Plant collecting and Documentation. In Plant Systematics; Simpson, M.G., Ed.; Academic Press: Cambridge, MA, USA, 2010; pp. 627–635. [Google Scholar]
  36. WinFolia TM; Version PRO 2005b; Regent Instruments Inc.: Quebec, QC, Canada, 2001.
  37. Sokal, R.R.; Rohlf, F.J. Biometry, 4th ed.; W.H. Freeman and Company: New York, NY, USA, 2012; p. 937. [Google Scholar]
  38. StatSoft, Inc. STATISTICA (Data Analysis Software System), Version 13; StatSoft, Inc.: Tulsa, OK, USA, 2018. [Google Scholar]
  39. Rohlf, F.J. NTSYS-pc: Numerical Taxonomy and Multivariate Analysis System, Version 2.2; Applied Biostatistics Inc.: New York, NY, USA, 2009; p. 44. [Google Scholar]
  40. Walters, S.M. Aceraceae. In Flora Europaea; Tutin, T.G., Heywood, V.H., Eds.; Cambridge University Press: Cambridge, UK, 1968; pp. 237–239. [Google Scholar]
  41. Idžojtić, M. Dendrologija—List [Dendrology—Leaf]; Fakultet Šumarstva i Drvne Tehnologije: Zagreb, Croatia, 2009; pp. 96–97. [Google Scholar]
  42. Kidner, C.A.; Umbreen, S. Why is leaf shape so variable? Int. J. Plant. Dev. Biol. 2010, 4, 64–75. [Google Scholar]
  43. Filartiga, A.L.; Klimeš, A.; Altman, J.; Nobis, M.P.; Crivellaro, A.; Schweingruber, F.; Doležal, J. Comparative anatomy of leaf petioles in temperate trees and shrubs: The role of plant size, environment and phylogeny. Ann. Bot. 2022, 129, 567–582. [Google Scholar] [CrossRef] [PubMed]
  44. Li, Y.; Kang, X.; Zhou, J.; Zhao, Z.; Zhang, S.; Bu, H.; Qi, W. Geographic Variation in the Petiole–Lamina Relationship of 325 Eastern Qinghai–Tibetan Woody Species: Analysis in Three Dimensions. Front. Plant Sci. 2021, 12, 748125. [Google Scholar] [CrossRef] [PubMed]
  45. Weijschedé, J.; Antonise, K.; de Caluwe, H.; de Kroon, H.; Huber, H. Effects of cell number and cell size on petiole length variation in a stoloniferous herb. Am. J. Bot. 2008, 95, 41–49. [Google Scholar] [CrossRef] [PubMed]
  46. Idžojtić, M. Dendrology: Cones, Flowers, Fruits and Seeds; Elsevier-Academic Press: London, UK; San Diego, CA, USA; Cambridge, UK; Oxford, UK, 2019; p. 800. [Google Scholar]
  47. Ren, J.; Ji, X.; Wang, C.; Hu, J.; Nervo, G.; Li, J. Variation and genetic parameters of leaf morphological traits of eight families from Populus simonii × P. nigra. Forests 2020, 11, 1319. [Google Scholar] [CrossRef]
  48. Krauze-Michalska, E.; Boratyńska, K. European geography of Alnus incana leaf variation. Plant Biosyst. 2013, 147, 601–610. [Google Scholar] [CrossRef]
  49. Koyama, K.; Smith, D.D. Scaling the leaf length-times-width equation to predict total leaf area of shoots. Ann. Bot. 2022, 130, 215–230. [Google Scholar] [CrossRef]
  50. Fishler, M.; Goldschmidt, E.E.; Monselise, S.P. Leaf Area and Fruit Size on Girdled Grapefruit Branches. J. Am. Soc. Hortic. Sci. 1983, 108, 218–221. [Google Scholar] [CrossRef]
  51. Wang, C.; Gong, H.; Feng, M.; Tian, C. Phenotypic variation in leaf, fruit and seed traits in natural populations of Eucommia ulmoides, a relict Chinese endemic tree. Forests 2023, 14, 462. [Google Scholar] [CrossRef]
  52. Li, H.; Wang, R.; Tian, Z.; Xu, J.; Sun, W.; Duan, R.; Fu, H.; Li, Y.; Zhang, Y.; Dong, L. Phenotypic variation and diversity in fruit, leaf, fatty acid, and their relationships to geoclimatic factors in seven natural populations of Malania oleifera Chun et S.K. Lee. Forests 2022, 13, 1733. [Google Scholar] [CrossRef]
  53. Chitwood, D.H.; Sinha, N.R. Evolutionary and Environmental Forces Sculpting Leaf Development. Curr. Biol. 2016, 26, R297–R306. [Google Scholar] [CrossRef] [PubMed]
  54. Wright, I.J.; Reich, P.B.; Westboy, M.; Ackerly, D.; Baruch, Z.; Bongers, F.; Cavender-Bares, J.; Chapin, T.; Cornelissen, J.H.C.; Diemer, M.; et al. The world-wide leaf economics spectrum. Nature 2004, 428, 821–827. [Google Scholar] [CrossRef] [PubMed]
  55. Shipley, B.; Lechowicz, M.J.; Wright, I.; Reich, P.B. Fundamental Trade-offs Generating the Worldwide Leaf Economics Spectrum. Ecology 2006, 87, 535–541. [Google Scholar] [CrossRef] [PubMed]
  56. Cescatti, A.; Zorer, R. Structural acclimation and radiation regime of silver fir (Abies alba Mill.) shoots along a light gradient. Plant Cell Environ. 2003, 26, 429–442. [Google Scholar] [CrossRef]
  57. Robakowski, P.; Wyka, T.; Samardakiewicz, S.; Kierzkowski, D. Growth, photosynthesis, and needle structure of silver fir (Abies alba Mill.) seedlings under different canopies. For. Ecol. Manag. 2004, 201, 221–227. [Google Scholar] [CrossRef]
  58. Petriţan, A.M.; von Lüpke, B.; Petriţan, I.C. Influence of light availability on growth, leaf morphology and plant architecture of beech (Fagus sylvatica L.), maple (Acer pseudoplatanus L.) and ash (Fraxinus excelsior L.) saplings. Eur. J. For. Res. 2009, 128, 61–74. [Google Scholar] [CrossRef]
  59. Ahmad, K.S.; Hameed, M.; Fatima, S.; Ashraf, M.; Ahmad, F.; Naseer, M.; Akhtar, N. Morpho-anatomical and physiological adaptations to high altitude in some Aveneae grasses from Neelum Valley, Western Himalayan Kashmir. Acta Physiol. Plant. 2016, 38, 93. [Google Scholar] [CrossRef]
  60. Liu, W.; Zheng, L.; Qi, D. Variation in leaf traits at different altitudes reflects the adaptive strategy of plants to environmental changes. Ecol. Evol. 2020, 10, 8166–8175. [Google Scholar] [CrossRef]
  61. Adamidis, G.C.; Varsamis, G.; Tsiripidis, I.; Dimitrakopoulos, P.G.; Papageorgiou, A.C. Patterns of leaf morphological traits of beech (Fagus sylvatica L.) along an altitudinal gradient. Forests 2021, 12, 1297. [Google Scholar] [CrossRef]
  62. Codarin, S.; Galopin, G.; Chasseriaux, G. Effect of air humidity on the growth and morphology of Hydrangea macrophylla L. Sci. Hortic. 2006, 108, 303–309. [Google Scholar] [CrossRef]
  63. Hovenden, M.J.; Schoor, J.K.V.; Osanai, Y. Relative humidity has dramatic impacts on leaf morphology but little effect on stomatal index or density in Nothofagus cunninghamii (Nothofagaceae). Aust. J. Bot. 2012, 60, 700–706. [Google Scholar] [CrossRef]
  64. Goba, K.A.E.; Kouonon, L.C.; Koffi, K.G.; Bony, B.S.; Diomandé, I.; Sylvère Sié, R. Morphological diversity within Pterocarpus erinaceus Poir. (Fabaceae), an overexploited species in the savannahs of Côte d’Ivoire. Am. J. Plant Sci. 2019, 10, 95460. [Google Scholar] [CrossRef]
  65. Sun, W.; Yuan, X.; Liu, Z.-J.; Lan, S.; Tsai, W.-C.; Zou, S.-Q. Multivariate analysis reveals phenotypic diversity of Euscaphis japonica population. PLoS ONE 2019, 14, e0219046. [Google Scholar] [CrossRef] [PubMed]
  66. Donskih, V.G.; Aniskina, T.S.; Kryuchkova, V.A. Variability and correlations of traits of the flower, leaf and fruit of Malus zumi. IOP Conf. Ser. Earth Environ. Sci. 2022, 1010, 012147. [Google Scholar] [CrossRef]
  67. Alcantra-Ayala, O.; Oyama, K.; Rios-Munoz, C.A.; Rivas, G.; Ramirez-Barahona, S.; Luna-Vega, I. Morphological variation of leaf traits in the Ternstroemia lineata species complex (Ericales: Penthyphylaceae) in response to geographic and climatic variation. Peer J. 2020, 8, e8307. [Google Scholar] [CrossRef] [PubMed]
  68. Lechowicz, M.J. Why do temperate deciduous trees leaf out at different times? Adaptation and ecology of forest communities. Am. Nat. 1984, 124, 821–842. [Google Scholar] [CrossRef]
  69. Pluess, A.R.; Schütz, W.; Stöcklin, J. Seed weight increases with altitude in the Swiss Alps between related species but not among populations of individual species. Oecologia 2005, 144, 55–61. [Google Scholar] [CrossRef]
  70. Li, F.; Bao, W. Elevational trends in leaf size of Campylotropis polyantha in the arid Minjiang River valley, SW China. J. Arid. Environ. 2014, 108, 1–9. [Google Scholar] [CrossRef]
  71. Xu, Y.; Woeste, K.; Cai, N.; Kang, X.; Li, G.; Chen, S.; Duan, A. Variation in needle and cone traits in natural populations of Pinus yunnanensis. J. For. Res. 2016, 27, 41–49. [Google Scholar] [CrossRef]
  72. Wang, L.; Li, J.; Zhao, J.; He, C. Evolutionary developmental genetics of fruit morphological variation within the Solanaceae. Front. Plant Sci. 2015, 6, 248. [Google Scholar] [CrossRef] [PubMed]
  73. Zhang, T.; Hong, Y.; Zhang, X.; Yuan, X.; Chen, S. Relationship between key environmental factors and the architecture of fruit shape and size in near-isogenic lines of cucumber (Cucumus sativus L.). Int. J. Mol. Sci. 2022, 23, 14033. [Google Scholar] [CrossRef] [PubMed]
  74. Bolvanský, M.; Užík, M. Morphometric variation and differentiation of European chestnut (Castanea sativa) in Slovakia. Biologia 2005, 60, 423–429. [Google Scholar]
  75. Idžojtić, M.; Zebec, M.; Poljak, I.; Medak, J. Variation of sweet chestnut (Castanea sativa Mill.) populations in Croatia according to the morphology of fruits. Sauteria 2009, 18, 232–333. [Google Scholar]
  76. Poljak, I.; Idžojtić, M.; Zebec, M.; Perković, N. The variability of European sweet chestnut (Castanea sativa Mill.) in the region of northwest Croatia according to morphology of fruits. Šumar. List 2012, 136, 479–489. [Google Scholar]
  77. Poljak, I.; Vahčić, N.; Liber, Z.; Šatović, Z.; Idžojtić, M. Morphological and chemical variation of wild sweet chestnut (Castanea sativa Mill.) populations. Forests 2022, 13, 55. [Google Scholar] [CrossRef]
  78. Zebec, M.; Idžojtić, M.; Mihaldinec, I.; Trogrlić, N. Morfološka varijabilnost brijesta veza (Ulmus laevis Pall.) u Hrvatskoj. In Proceedings of the 5th Croatian Botanical Symposium, Primošten, Croatia, 22–25 September 2016; Rešetnik, I., Ljubešić, Z., Eds.; Hrvatsko Botaničko Društvo: Zagreb, Croatia, 2016; pp. 23–24. [Google Scholar]
  79. Poljak, I.; Vahčić, N.; Liber, Z.; Tumpa, K.; Pintar, V.; Zegnal, I.; Vidaković, A.; Valković, B.; Kajba, D.; Idžojtić, M. Morphological and chemical diversity and antioxidant capacity of the service tree (Sorbus domestica L.) fruits from two eco-geographical regions. Plants 2021, 10, 1691. [Google Scholar] [CrossRef]
  80. Ramsey, M.W.; Cairns, S.C.; Vaughton, G.V. Geographic variation in morphological and reproductive characters of coastal and tableland populations of Blandfordia grandiflora (Liliaceae). Plant Syst. Evol. 1994, 192, 215–230. [Google Scholar] [CrossRef]
  81. Valladares, F.; Matesanz, S.; Guilhaumon, F.; Araújo, M.B.; Balaguer, L.; Benito-Garzón, M.; Cornwell, W.; Gianoli, E.; van Kleunen, M.; Naya, D.E.; et al. The effects of phenotypic plasticity and local adaptation on forecasts of species range shifts under climate change. Ecol. Lett. 2014, 17, 1351–1364. [Google Scholar] [CrossRef]
  82. Boratyński, A.; Jasińska, A.K.; Marcysiak, K.; Mazur, M.; Romo, A.M.; Boratyńska, K.; Sobierajska, K.; Iszkuło, G. Morphological differentiation supports the genetic pattern of the geographic structure of Juniperus thurifera (Cupressaceae). Plant Syst. Evol. 2013, 299, 773–784. [Google Scholar] [CrossRef]
  83. Jasińska, A.K.; Boratyńska, K.; Dering, M.; Sobierajska, K.I.; Ok, T.; Romo, A.; Boratyński, A. Distance between south-European and south-west Asiatic refugial areas involved morphological differentiation: Pinus sylvestris case study. Plant Syst. Evol. 2014, 300, 1487–1502. [Google Scholar] [CrossRef]
  84. Parisod, C.; Joost, S. Divergent selection in trailing- versus leading-edge populations of Biscutella laevigata. Ann. Bot. 2010, 105, 655–660. [Google Scholar] [CrossRef] [PubMed]
Plants 13 00320 g001
Figure 1. Locations of the eight sampled Tatar maple (Acer tataricum L. subsp. tataricum) populations. Populations: P1—Odransko polje; P2—Lipovljani; P3—Veliki Grđevac; P4—Mali Grđevac; P5—Grubišno Polje; P6—Virovitica; P7—Požega; P8—Županja. The blue area represents the natural habitat of Tatar maple according to Bartha [33].
Figure 1. Locations of the eight sampled Tatar maple (Acer tataricum L. subsp. tataricum) populations. Populations: P1—Odransko polje; P2—Lipovljani; P3—Veliki Grđevac; P4—Mali Grđevac; P5—Grubišno Polje; P6—Virovitica; P7—Požega; P8—Županja. The blue area represents the natural habitat of Tatar maple according to Bartha [33].
Plants 13 00320 g001
Plants 13 00320 g002
Figure 2. Measured leaf and fruit traits: LA—leaf area; LL—leaf length; MLW—maximum leaf width; PMLW—leaf length, measured from the leaf base to the point of maximum leaf width; LW1—leaf blade width at 50% of leaf blade length; LW2—leaf blade width at 90% of leaf blade length; PL—petiole length; LA1—angle closed by the main leaf vein (the centre of the leaf blade) and the line connecting the leaf blade base to a set point on the leaf margin at 10% of the total leaf blade length; and LA2—angle closed by the main leaf vein (the centre of the leaf blade) and the line connecting the leaf blade base to a set point on the leaf margin at 25% of the total leaf blade length; MA—mericarp area; ML—mericarp length; MMW—maximum mericarp width; PMMW—mericarp length, measured from the mericarp base to the point of maximum mericarp width; NW90—mericarp width at 90% of the fruit length; NL—nut length; NW—nut width; WA—wing angle.
Figure 2. Measured leaf and fruit traits: LA—leaf area; LL—leaf length; MLW—maximum leaf width; PMLW—leaf length, measured from the leaf base to the point of maximum leaf width; LW1—leaf blade width at 50% of leaf blade length; LW2—leaf blade width at 90% of leaf blade length; PL—petiole length; LA1—angle closed by the main leaf vein (the centre of the leaf blade) and the line connecting the leaf blade base to a set point on the leaf margin at 10% of the total leaf blade length; and LA2—angle closed by the main leaf vein (the centre of the leaf blade) and the line connecting the leaf blade base to a set point on the leaf margin at 25% of the total leaf blade length; MA—mericarp area; ML—mericarp length; MMW—maximum mericarp width; PMMW—mericarp length, measured from the mericarp base to the point of maximum mericarp width; NW90—mericarp width at 90% of the fruit length; NL—nut length; NW—nut width; WA—wing angle.
Plants 13 00320 g002
Plants 13 00320 g003
Figure 3. Isolation by distance (IBD) in Tatar maple populations. Scatter plot of simple Mantel test showing the relationship between geographic and leaf phenotypic distances (r = 0.585, p = 0.0092, R2 = 0.3421).
Figure 3. Isolation by distance (IBD) in Tatar maple populations. Scatter plot of simple Mantel test showing the relationship between geographic and leaf phenotypic distances (r = 0.585, p = 0.0092, R2 = 0.3421).
Plants 13 00320 g003
Table 1. Results of the descriptive statistical analysis for studied populations and leaf morphometric traits. Leaf morphometric traits analysed: LA—leaf area; LL—leaf length; MLW—maximum leaf width; PMLW—leaf length, measured from the leaf base to the point of maximum leaf width; LW1—leaf blade width at 50% of the leaf blade length; LW2—leaf blade width at 90% of the leaf blade length; PL—petiole length; FC—form coefficient; LA1—angle closed by the main leaf vein (the centre of the leaf blade) and the line connecting the leaf blade base to a set point on the leaf margin at 10% of the total leaf blade length; and LA2—angle closed by the main leaf vein (the centre of the leaf blade) and the line connecting the leaf blade base to a set point on the leaf margin at 25% of the total leaf blade length. Descriptive parameters: M—arithmetic mean and CV—coefficient of variation (%). Populations: P1—Odransko polje; P2—Lipovljani; P3—Veliki Grđevac; P4—Mali Grđevac; P5—Grubišno Polje; P6—Virovitica; P7—Požega; P8—Županja.
Table 1. Results of the descriptive statistical analysis for studied populations and leaf morphometric traits. Leaf morphometric traits analysed: LA—leaf area; LL—leaf length; MLW—maximum leaf width; PMLW—leaf length, measured from the leaf base to the point of maximum leaf width; LW1—leaf blade width at 50% of the leaf blade length; LW2—leaf blade width at 90% of the leaf blade length; PL—petiole length; FC—form coefficient; LA1—angle closed by the main leaf vein (the centre of the leaf blade) and the line connecting the leaf blade base to a set point on the leaf margin at 10% of the total leaf blade length; and LA2—angle closed by the main leaf vein (the centre of the leaf blade) and the line connecting the leaf blade base to a set point on the leaf margin at 25% of the total leaf blade length. Descriptive parameters: M—arithmetic mean and CV—coefficient of variation (%). Populations: P1—Odransko polje; P2—Lipovljani; P3—Veliki Grđevac; P4—Mali Grđevac; P5—Grubišno Polje; P6—Virovitica; P7—Požega; P8—Županja.
TraitDescriptive ParametersPopulationTotal
P1P2P3P4P5P6P7P8
LA (cm2)M20.5121.7618.2218.2017.3119.7619.6414.0418.38
CV34.9626.6640.8533.8238.3936.0829.3725.5036.45
LL (cm)M6.706.756.306.266.086.366.495.766.29
CV18.5914.0020.2716.0618.2817.8514.2013.7617.64
MLW (cm)M4.284.553.964.073.954.284.293.564.08
CV20.7215.7323.1819.7121.3620.1017.4114.6820.70
PMLW (cm)M2.332.322.102.051.932.082.151.872.08
CV26.7220.7829.4225.8028.2124.7524.2124.3026.75
LW1 (cm)M3.734.003.503.583.403.753.692.983.54
CV19.8417.1722.9820.5122.6020.7517.6515.1821.56
LW2 (cm)M0.790.890.830.810.810.820.800.600.78
CV31.0733.3835.2036.9140.9638.7436.9826.6137.40
PL (cm)M4.635.304.053.994.284.564.693.634.30
CV28.7522.3133.6030.8728.4429.9132.8525.9431.47
FCM0.650.620.600.620.670.690.630.550.63
CV13.5214.9515.3415.1812.2613.1917.3715.6116.29
LA1 (°)M67.4168.7367.2167.8568.1069.1068.1467.0967.91
CV3.863.724.424.834.134.074.294.734.44
LA2 (°)M49.5851.0548.9350.0050.3151.3250.7048.8350.00
CV5.105.766.806.666.705.586.476.476.50
Table 2. Results of the descriptive statistical analysis for studied populations and fruit morphometric traits. Fruit morphometric traits analysed: MA—mericarp area; ML—mericarp length; MMW—maximum mericarp width; PMMW—mericarp length, measured from the mericarp base to the point of maximum mericarp width; NW90—mericarp width at 90% of fruit length; NL—nut length; NW—nut width; WA—wing angle. Descriptive parameters: M—arithmetic mean and CV—coefficient of variation (%). Populations: P1—Odransko polje; P2—Lipovljani; P3—Veliki Grđevac; P4—Mali Grđevac; P5—Grubišno Polje; P6—Virovitica; P7—Požega; P8—Županja.
Table 2. Results of the descriptive statistical analysis for studied populations and fruit morphometric traits. Fruit morphometric traits analysed: MA—mericarp area; ML—mericarp length; MMW—maximum mericarp width; PMMW—mericarp length, measured from the mericarp base to the point of maximum mericarp width; NW90—mericarp width at 90% of fruit length; NL—nut length; NW—nut width; WA—wing angle. Descriptive parameters: M—arithmetic mean and CV—coefficient of variation (%). Populations: P1—Odransko polje; P2—Lipovljani; P3—Veliki Grđevac; P4—Mali Grđevac; P5—Grubišno Polje; P6—Virovitica; P7—Požega; P8—Županja.
TraitDescriptive ParametersPopulationTotal
P1P2P3P4P5P6P7P8
MA (cm2)M2.962.302.592.452.582.802.931.872.53
CV16.4721.1422.3816.9622.6719.9627.2426.4325.96
ML (cm)M3.483.153.223.183.213.353.182.753.18
CV7.559.2910.337.2410.389.9414.0813.9212.38
MMW (cm)M1.261.061.191.161.171.241.350.971.17
CV12.4117.4718.2913.5015.3013.1617.0115.6618.18
PMMW (cm)M2.432.262.302.302.292.442.261.982.27
CV9.8511.9513.797.8311.139.6114.9114.0313.33
NW90 (cm)M0.900.780.860.830.840.931.000.710.85
CV16.8417.3619.2814.2818.3514.6416.8716.3819.39
NL (cm)M1.301.251.311.241.191.281.261.031.23
CV9.898.8010.438.548.338.1512.8713.5912.67
NW (cm)M0.660.640.660.650.570.670.620.590.63
CV9.5613.0414.3115.3411.0513.6615.6213.0614.54
WAM59.6771.5670.7566.9955.0273.3565.9764.1366.35
CV15.1315.9420.9019.8920.3116.2317.7616.7519.88
Table 3. Hierarchical analysis of variance for leaf morphometric traits. Leaf morphometric traits analysed: LA—leaf area; LL—leaf length; MLW—maximum leaf width; PMLW—leaf length, measured from the leaf base to the point of maximum leaf width; LW1—leaf blade width at 50% of the leaf blade length; LW2—leaf blade width at 90% of the leaf blade length; PL—petiole length; FC—form coefficient; LA1—angle closed by the main leaf vein (the centre of the leaf blade) and the line connecting the leaf blade base to a set point on the leaf margin at 10% of the total leaf blade length; and LA2—angle closed by the main leaf vein (the centre of the leaf blade) and the line connecting the leaf blade base to a set point on the leaf margin at 25% of the total leaf blade length. n.s. not significant at p > 0.05; * significant at 0.01 < p < 0.05; ** significant at 0.001 < p < 0.01; *** significant at p < 0.001.
Table 3. Hierarchical analysis of variance for leaf morphometric traits. Leaf morphometric traits analysed: LA—leaf area; LL—leaf length; MLW—maximum leaf width; PMLW—leaf length, measured from the leaf base to the point of maximum leaf width; LW1—leaf blade width at 50% of the leaf blade length; LW2—leaf blade width at 90% of the leaf blade length; PL—petiole length; FC—form coefficient; LA1—angle closed by the main leaf vein (the centre of the leaf blade) and the line connecting the leaf blade base to a set point on the leaf margin at 10% of the total leaf blade length; and LA2—angle closed by the main leaf vein (the centre of the leaf blade) and the line connecting the leaf blade base to a set point on the leaf margin at 25% of the total leaf blade length. n.s. not significant at p > 0.05; * significant at 0.01 < p < 0.05; ** significant at 0.001 < p < 0.01; *** significant at p < 0.001.
TraitVariance Componentdf% VariationFp
LAAmong populations711.027.68***
Within populations8917.405.87***
Error 71.58
LLAmong populations76.695.03***
Within populations8916.495.29***
Error 76.82
MLWAmong populations711.348.12***
Within populations8916.515.58***
Error 72.15
PMLWAmong populations76.104.40***
Within populations8917.225.49***
Error 76.68
LW1Among populations713.397.97***
Within populations8921.697.69***
Error 64.92
LW2Among populations77.154.11***
Within populations8926.819.15***
Error 66.03
PLAmong populations713.377.57***
Within populations8920.757.30***
Error 65.88
FCAmong populations716.286.75***
Within populations8936.6416.60***
Error 47.08
LA1Among populations74.582.60**
Within populations8932.6111.40***
Error 62.81
LA2Among populations76.523.43**
Within populations8931.1911.04***
Error 62.29
Table 4. Hierarchical analysis of variance for fruit morphometric traits. Fruit morphometric traits analysed: MA—mericarp area; ML—mericarp length; MMW—maximum mericarp width; PMMW—mericarp length, measured from the mericarp base to the point of maximum mericarp width; NW90—mericarp width at 90% of fruit length; NL—nut length; NW—nut width; WA—wing angle. n.s. not significant at p > 0.05; * significant at 0.01 < p < 0.05; ** significant at 0.001 < p < 0.01; *** significant at p < 0.001.
Table 4. Hierarchical analysis of variance for fruit morphometric traits. Fruit morphometric traits analysed: MA—mericarp area; ML—mericarp length; MMW—maximum mericarp width; PMMW—mericarp length, measured from the mericarp base to the point of maximum mericarp width; NW90—mericarp width at 90% of fruit length; NL—nut length; NW—nut width; WA—wing angle. n.s. not significant at p > 0.05; * significant at 0.01 < p < 0.05; ** significant at 0.001 < p < 0.01; *** significant at p < 0.001.
TraitVariance Componentdf% VariationFp
MAAmong populations725.877.25***
Within populations8951.5946.83***
Error 22.54
MLAmong populations724.427.32***
Within populations8949.7239.52***
Error 25.86
MMWAmong populations725.876.90***
Within populations8952.6450.00***
Error 21.49
PMMWAmong populations719.687.23***
Within populations8941.1622.07***
Error 39.16
NW90Among populations723.657.73***
Within populations8941.8325.24***
Error 34.52
NLAmong populations731.509.86***
Within populations8947.4746.21***
Error 21.03
NWAmong populations711.343.40**
Within populations8954.7133.25***
Error 33.95
WAAmong populations715.743.84**
Within populations8962.1757.24***
Error 22.09
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Poljak, I.; Vidaković, A.; Benić, L.; Tumpa, K.; Idžojtić, M.; Šatović, Z. Patterns of Leaf and Fruit Morphological Variation in Marginal Populations of Acer tataricum L. subsp. tataricum. Plants 2024, 13, 320. https://doi.org/10.3390/plants13020320

AMA Style

Poljak I, Vidaković A, Benić L, Tumpa K, Idžojtić M, Šatović Z. Patterns of Leaf and Fruit Morphological Variation in Marginal Populations of Acer tataricum L. subsp. tataricum. Plants. 2024; 13(2):320. https://doi.org/10.3390/plants13020320

Chicago/Turabian Style

Poljak, Igor, Antonio Vidaković, Luka Benić, Katarina Tumpa, Marilena Idžojtić, and Zlatko Šatović. 2024. "Patterns of Leaf and Fruit Morphological Variation in Marginal Populations of Acer tataricum L. subsp. tataricum" Plants 13, no. 2: 320. https://doi.org/10.3390/plants13020320

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop