Next Article in Journal
Response Surface Optimization of Extraction Conditions and In Vitro Antioxidant and Antidiabetic Evaluation of an Under-Valued Medicinal Weed, Mimosa pudica
Previous Article in Journal
A Systemic View of Carbohydrate Metabolism in Rice to Facilitate Productivity
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Morphological and Chemical Diversity and Antioxidant Capacity of the Service Tree (Sorbus domestica L.) Fruits from Two Eco-Geographical Regions

1
Faculty of Forestry and Wood Technology, University of Zagreb, Svetošimunska Cesta 23, 10000 Zagreb, Croatia
2
Faculty of Food Technology and Biotechnology, University of Zagreb, Pierottijeva 6, 10000 Zagreb, Croatia
3
Department of Biology, Faculty of Science, University of Zagreb, Marulićev Trg 9a, 10000 Zagreb, Croatia
4
Centre of Excellence for Biodiversity and Molecular Plant Breeding, Svetošimunska Cesta 25, 10000 Zagreb, Croatia
5
Ministry of Economy and Sustainable Development, Institute for Environment and Nature, Nature Sector, Radnička Cesta 80, 10000 Zagreb, Croatia
6
Croatian Forest Research Institute, Cvjetno Naselje 41, 10450 Jastrebarsko, Croatia
*
Author to whom correspondence should be addressed.
Plants 2021, 10(8), 1691; https://doi.org/10.3390/plants10081691
Submission received: 21 July 2021 / Revised: 12 August 2021 / Accepted: 14 August 2021 / Published: 17 August 2021
(This article belongs to the Section Plant Ecology)

Abstract

:
Service tree, Sorbus domestica L., is a rare and neglected wild fruit tree species of southern and central Europe. Being distributed in different eco-geographical regions, with fragmented and low-density populations, S. domestica represents an interesting model case for investigating patterns of within- and between-population diversity at geographical and environmental scales. This study aimed to analyze the proximate composition, antioxidant activity, and morphometric fruit characteristics. We examined the diversity and population divergences of 49 S. domestica individuals originating from seven populations across continental and Mediterranean eco-geographical regions. In addition, tests of isolation by distance and environment were performed to detect the magnitude of divergence explained by geographic and environmental variables. Significant differences between the studied populations were found in almost all of the studied morphometric and chemical fruit characteristics. The studied service tree populations were characterized by high phenotypic variation despite the low number of trees per population. Model-based population structure analysis using morphometric and chemical fruit characteristics revealed three groups of service tree populations. We concluded that non-effective pollen and seed dispersal along with genetic drift and specific environmental factors resulted in a distinct phenotype with a specific chemical composition in the isolated island population. In addition, a pattern of isolation by the environment was revealed. We infer that morphological and chemical differences between the studied populations in the true service tree from different eco-geographical regions were mediated by adaptation to the specific environmental conditions.

1. Introduction

Forest fruit trees are species that play an important role in maintaining biodiversity, by enriching the gene pool of the forest ecosystem and with other trees and shrubs, ensuring the hardiness and health of the whole forest, as well as increasing the soil quality [1]. Due to their adaptability and resilience, the promotion of forest fruit tree species can be considered as one of the most promising efforts of climate change mitigation, and planting them is an efficient way of restoring degraded forests [2]. In addition, forest fruits represent an additional food source for wild animals, hence, being rich in antioxidants, minerals, and vitamins, they are recommended for human consumption as well [3]. As such, they have a large potential in the sustainable development of rural areas [4] and their usage and popularization are in line with sustainable resource management [5].
One of those forest fruit trees, almost forgotten in the previous decades, is the service tree. The service tree (Sorbus domestica L.) is an insect-pollinated, thermophile, and medium-sized deciduous fruit tree species of the Rosaceae family. It is found in southern and central Europe, with the center of distribution in the Balkan Peninsula, Italy, and southern France, sporadically appearing in North Africa and the Caucasus as well [6]. The extent of its original range is still unclear since it has been planted and sub-spontaneously spread since Roman times [7,8,9]. It is a light-demanding species, tolerant to different types of soil and occurring on a wide variety of sites [10,11]. Its distribution is scattered in two different eco-geographical regions, where it forms small and fragmented populations [12]. It can be found in a mild temperate climate on warm and south-facing slopes, and in the Mediterranean climate on dry to extremely dry, nutrient-poor sites. Although vigorous and relatively fast-growing, it is also a poor competitor, making it a rare species with low population density, with fruit-bearing trees being even scarcer.
Service tree truly does provide a number of services-high-quality wood, edible fruits, and health benefits, to name a few [13,14]. In addition, it is also a valuable species in ecosystems and landscapes, and it can be used as an ornamental tree in urban areas. Nevertheless, it is usually known as a fruit tree species. Its pear-shaped (f. domestica = f. pyriformis (Hayne) Rehder) or apple-shaped (f. pomifera (Hayne) Rehder, = f. maliformis (Kirchn. et J.Eichler) Hegi) fruits can be used for various culinary and medicinal purposes [15]. The fruits are pomes olive-brown or yellowish in color, flushed with red on the side facing the sun, often glaucous, with numerous large lenticels, maturing in September to October [16]. The weight of one pome is 15–32 g and its diameter is 1.2–3 (4.5) cm. Fruit flesh (enlarged fleshy hypanthium) is astringent and hard at maturity, edible after bletting by frost or in the state of incipient decay, with numerous sclereids. The fruit usually contains 1–6 (10) broadly obovoid, flattened, and reddish-brown seeds. The seeds are dispersed by birds and mammals [17]. The fruits can be eaten fresh or used to produce marmalades, juices, fruit wine, and to conserve apple cider. Although service tree fruits are nutrient-rich and suitable for processing, they are currently undervalued by the consumers and will need to gain in popularity to compete with the wide range of attractive and tasty fruits currently on the market.
Edible fruits characterizing this species make it an ideal model for the determination and further analysis of morphological and chemical diversity of populations. To date, only several studies of the morphological characteristics [18,19,20,21,22] and/or chemical composition [23,24,25,26,27] of the fruits have been carried out on the service tree. The fruit’s chemical composition was studied by different analyses including total phenols, flavonoids and tannins, macro- and micro-elements, and determination of antioxidant activity. However, to the best of our knowledge, until now no data have been published on the proximate composition in this species. The morphological variability of S. domestica fruit across a larger geographical range has not been studied to date either. In addition, several studies have been performed on S. domestica, including leaf morphometric analysis [28,29] and leaf antioxidant activity [30].
In this study, we combined morphological and chemical fruit traits to examine the diversity and divergence of seven service tree populations distributed across the Croatian continental and Mediterranean eco-geographical regions. In total, we included fruit samples from 49 trees. The small number of sampled trees was a direct result of the small density of populations [6,8,9,29], and the very low number of fruit-bearing trees. Our main objectives were to evaluate morphological diversity, fruit proximate composition, total phenols, antioxidants activity, and population structure, and to test the correlations between geographic, environmental, morphological, and chemical variation in S. domestica. Chemical differentiation was examined using ten chemical characteristics, and morphological variation was assessed for ten fruit traits. Finally, a variety of multivariate analyses and Mantel tests were used to examine the roles of geographic and environmental isolation in determining the morphological and chemical diversity in wild fruit tree species.

2. Results

2.1. Climate Differences among Sampling Sites

Nineteen bioclimatic variables were used to describe the environmental differences between seven studied populations (Table S1). The first principal component (PC1) explained 65.5% of the total variation and clearly separated the continental (P01 and P02) and Mediterranean (P03-P07) populations (Figure S1 and Table S2). The continental region is characterized by cold winters, with a more even distribution of rainfall throughout the year. The annual mean temperature is lower than that in the Mediterranean region, as well as the mean temperature of the coldest month. In contrast, the Mediterranean region features hot to very hot summers and mild winters. A distinctive feature of the Mediterranean region is the seasonality of rainfall distribution, with the majority of precipitation in the coldest quarter of the year.

2.2. Fruit Morphometric Characteristics

Based on all parameters of fruit size, the continental populations were shown to have larger fruits than the Mediterranean populations (Table 1, Figure S2). Fruit mass had a mean value of 7.32 g, and it was the second most variable trait, with a mean coefficient of variation of 36.04%. The least variable traits were fruit (FL) and seed (SL) length and fruit (MFW) and seed (SW) width.
Fruit width and fruit length ratio (FW/FL) values for the continental populations and Mediterranean population Istria were approximately 1.0, and for other south Mediterranean populations, 1.1. Although fruits of slightly oblong shape, e.g., pear-like circular or oval, characterized southern Mediterranean populations, and apple-like shape continental and north Mediterranean population Istria, both types of fruits were found in all studied populations. Fruit shape was a far less variable characteristic, with a mean CV of 10.17%.
Populations from the continental region had wider seeds than the populations from the Mediterranean eco-geographical region. The mean seed number (NS) varied from 1.25 in population Novi Vinodolski to 3.66 in population Brač, with an overall mean value of 1.88. Seed number was by far the most variable trait, with a CV of 60.06%. In general, Mediterranean populations were characterized by a larger number of smaller and more variable seeds.
All fruit and seed parameters were significantly different between individual trees and populations, except fruit width and fruit length ratio (FW/FL) and seed length (SL) (Table 2). The majority of tested parameters were characterized by having most of the variability linked to intra-population variability.

2.3. Proximate Analysis and Acidity

The analyzed service tree fruits were characterized by high water content and sugar levels and very low fat and acidity values (Table 3). Differences between the studied populations were confirmed for all studied proximate constitutes, except for acidity. The mean water content value was 65.14 g per 100 g dm. Values were relatively uniform, and the average coefficient of variation was low (CV = 6.70%). In general, Mediterranean populations had somewhat lower water content values than continental populations. Sugar content had a mean value of 49.36 g per 100 g dm and a coefficient of variation of 9.33%. The lowest values were found in population Brač (42.88 g per 100 g dm), and the highest in population Novi Vinodolski (55.60 g per 100 g dm). The mean protein value was 14.69 g per 100 g dm. Protein content was characterized by relatively high variability, with a mean coefficient of variation of 20.27%. Both ash and fat content were very low, with mean values of 2.12 and 0.72 g per 100 g dm, respectively. Both of these parameters were very variable, with mean coefficients of variation for fat of 19.52%, and ash content being the most variable parameter in general, with a coefficient of variation of 34.23%. Mediterranean populations were characterized by significantly higher fat content than the continental populations. Acidity values ranged from 0.64 to 0.74, with a mean value of 0.68. The coefficient of variations showed intermediate variability of data, with a mean value of 14.90% and a range of 5.24–28.20%. The cellulose content varied between the studied populations without a clear pattern. The average cellulose content was 6.24 g per 100 g dm, with a mean coefficient of variation of 27.95%. The population Brač was characterized by the highest cellulose content, 9.06 g per 100 g dm.

2.4. Total Phenolics and Antioxidant Capacity

Total phenolic content mean values ranged from 3.50 mg GAE/g dm to 12.10 mg GAE/g dm (Table 3). Continental populations have the highest values, 11.11 mg GAE/g dm (Psunj) and 12.10 mg GAE/g dm (Tounj), whereas island population Brač and population Split have the lowest, 3.50 mg GAE/g dm and 5.62 mg GAE/g dm, respectively. Populations Psunj and Brač have the highest coefficient of variation, 63.19%, and 57.44%, while population Tounj has the lowest variability of 16.80%, followed by population Istria, 38.23%.
Measured DPPH values ranged from 4.40 to 17.76%. Populations Istria (17.76%) and Split (13.76%) have the highest values, and Brač (11.70%) and Tounj (4.40%) have the lowest. Measured DPPH values were less variable than those of total phenolic content, with a coefficient of variation of 44.00%. Mean FRAP values ranged from 2.41 to 2.91 mmol Fe2+. The highest value was calculated for population Tounj, and the lowest for population Psunj.

2.5. Population Structure, Isolation by Distance, and Environment

The structure of the seven service tree populations was inferred by the K-means clustering method. The most probable division was detected at K = 3, and the estimated population structure is shown in Figure 1A and Figure 2A. If the proportion of a certain population was equal to or higher than 0.70, it was assumed that the population belonged to one cluster, and if it was lower than 0.70, it was assumed that the population had a mixed origin. In both analyses, morphometric and chemical, we revealed the same geographical pattern. The samples from the populations Psunj and Tounj grouped into cluster A, the samples from island population Brač grouped into cluster B, and the samples from the coastal Mediterranean populations into cluster C. Only one population, P04 (Novi Vinodolski), was of mixed origin, with the dominant proportion from cluster C. The results obtained with the K-means clustering method were congruent with the results from the hierarchical clustering method (Figure 1B and Figure 2B) and Barrier software (Figure 3A,B). The existence of barriers was revealed between the continental and Mediterranean populations, as well as between the island population Brač and all other studied populations.
The first two components from a PC analysis of the phenotypic traits explained 59.2% and 20.6% of the total variation, respectively (Table S3 and Figure 4). PC1 was highly negatively correlated with the fruit size, and moderately positively with the seed number. On the other hand, PC2 was highly negatively correlated with the fruit shape, and moderately positively correlated with the seed size. Continental populations were associated with larger fruits with a smaller number of seeds, and Mediterranean populations with smaller fruits with a higher number of seeds.
PC analysis of all individuals showed that 50.2% of observed chemical variability was explained by the first two principal components (Table S4 and Figure 5). The first principal component (PC1) was strongly negatively correlated with the water and sugar content, and strongly positively with the ash and cellulose content. Strong negative correlations were observed between the second principal component (PC2) and the fat content and DPPH. Clear differences were revealed between the samples from the island Brač and other studied populations.
Discriminant analysis was performed to determine which of the morphological and chemical traits were the most useful for maximum discrimination between the three groups of service tree populations established by the K-means clustering method. Four out of nine morphological traits (number of seeds, fruit length, position of maximum fruit width, fruit width at 15% of fruit length), and four out of ten chemical characteristics (water, crude fat, cellulose, DPPH), were determined by stepwise discriminant analysis to be the best differentiating variables between the studied groups of service tree populations (Tables S5 and S6). The discriminant function based on morphometric and chemical traits showed a classification success of 87.75% and 97.96%, respectively (Figure 1C and Figure 2C). Overall, these results confirmed the usefulness of morphometric and chemical fruit traits in discrimination of the studied groups of service tree populations.
Neither morphological nor chemical distances were related to geographical distances (Figure 6A,C). On the other hand, there was an overall relationship between pairwise phenotypic and environmental distances matrices, suggesting a pattern of isolation by environment (Figure 6B). Although we did not find significant correlations between pairwise chemical and environmental distances matrices (Figure 6D), several chemical traits were closely related to the environmental variables, i.e., Mediterranean populations were characterized by lower water and total phenolic contents, and higher fat contents in comparison with the continental populations. In addition, the Mantel test identified significant correlations between the morphological and chemical distance matrices (Figure 6E).

3. Discussion

3.1. Fruit Morphometric Characteristics

Our results clearly demonstrate a significant phenotypic diversity of the S. domestica fruits. The analyzed fruit characters had diversity levels similar to those reported for other tree species with similar life characteristics, especially regarding the non-continuous population structure, such as Sorbus torminalis (L.) Crantz [31,32], Pyrus pyraster (L.) Burgsd. [33,34], Malus sylvestris (L.) Mill. [35], Prunus avium (L.) L. [36,37] and various Crataegus L. taxa [38,39].
In general, populations from the continental eco-geographical region had larger fruits than populations from the Mediterranean region. The values of particular morphological traits found in our study did not differ considerably from the values reported in botany and dendrology textbooks [16,40,41,42,43] or morphometric studies [18,19,20,21,22]. However, the range reported by other authors was significantly wider than the range obtained in this study, i.e., the maximum values for fruit size, noted by other authors, were somewhat higher than those obtained in this research. In previously published papers, analyzed fruits were collected in wild populations and orchards. Our samples, however, were collected only from wild populations, which most likely meant less human interference and sub-optimal growing conditions. Interestingly, the levels of diversity in cultivated populations were similar to those in our study for the wild populations. High diversity in cultivated populations was probably preserved due to the fact that true service tree has never been subjected to selective breeding intensively, even though it has often been grown as a fruit tree [28]. According to Bignami [44], it was only propagated generatively by seeds for centuries and therefore preserved wide and underexploited variation in orchards as well as in natural populations.
The results of the analysis of variance (ANOVA) were in line with the expectations of high morphological variation within populations and low differentiation between populations, as observed in rowans [28,29,31] and other insect-pollinated and animal-dispersed tree species [45,46,47,48]. The relatively high level of among-tree variation within the populations is probably a result of both phenotypic plasticity to specific micro-environmental conditions experienced by each tree, and genetic differentiation among individual trees [28,29,49]. The within-population variation was of similar magnitude in all studied populations, except for the island population Brač. Even though isolated island populations could be expected to have significantly lower variability compared to the continental populations [50,51,52], island population Brač demonstrated the highest variability levels. This is more likely linked to the microsite features of its habitat than to genetic variability. Habitat conditions in which island samples were collected varied greatly, ranging from deeper soils to shallow, skeletal soils, stretched across varying altitudes. This indicates that phenotypic plasticity might be of greater importance with regard to microhabitat variation within service tree populations [53]. In addition, the island population Brač was morphologically the most distinguished population out of those analyzed. This population was characterized by having the smallest fruit with the highest number of seeds per fruit.

3.2. Proximate Composition

To the best of our knowledge, this is the first report on the proximate chemical composition of service tree fruits. Service tree fruits are sweet due to their high sugar content, are rich in proteins, and have a high content of cellulose, which acts as dietary fiber and as such can enrich a diet. We revealed that the parameters with the highest mean values were the least variable, e.g., water and sugar, and the most variable parameters were those with the lowest values-ash, cellulose, fat, and acidity. In comparison with the closely related rowan (Sorbus aucuparia L.), service tree fruits are characterized by lower water content, significantly higher levels of carbohydrates and cellulose, and similar, very low, fat content [54,55,56]. Results of the chemical composition analysis in this paper demonstrate great potential for further research on the specifics and benefits of service tree fruits and indicates potential future uses of fruits in both fresh and processed food supplements and products, like low-fat, protein-rich food.
This research has found significant differences between the populations in the analyzed chemical composition of service tree fruits. In general, the continental populations have somewhat higher water content and significantly lower fat values from cuticular wax. The differences between continental and Mediterranean populations are readily explicable by the specific ecological conditions in which populations grow. The Mediterranean populations are exposed to drought and extremely high temperatures, meaning their fruit need protection from the thicker cuticular wax, which acts as a shield against desiccation. It is reported that this protective layer in many plant species is a key evolutionary innovation in plants [57], and it has been implicated in protection mechanisms against different environmental stressors, such as excessive ultraviolet radiation, high temperatures, and salinity, low temperature during the vegetation season, etc. The current relationships between fruit chemical composition and environmental factors in S. domestica may be the result of long-term adaptation to different habitats.
Similar to morphological variability findings, intrapopulation variability was similar among all analyzed populations, with population Brač being the most variable. In addition, Brač was also the most specific population in terms of chemical composition. This island population was characterized by the lowest water content and highest ash and cellulose contents.

3.3. Total Phenolic Content and Antioxidant Activity

Several studies have demonstrated the importance of service tree fruits as a good source of different bioactive phytochemicals [26,41,58,59]. It was reported that those constitutes are mostly affected by geographic origin and genotype [26,60], climatic environment [61,62,63], and stage of maturity [64,65,66]. In general, it is well known that phenolic compounds represent a chemical interface between plants and the environment [67] and are under the influence of various environmental factors, including soil composition, temperature, and rainfall [68].
Significant differences in total phenols were found between the Mediterranean and continental populations, with higher values observed in the continental populations. This was somewhat opposite to our expectations since higher total phenolic content is expected in the sites where water deficit, along with high temperatures and increased UV radiation, causes stress in plants [69,70,71]. When subjected to environmental stress, plants respond by increasing the production of secondary metabolites, including phenolic compounds. The effect of moderate drought stress on the increased production of these compounds has been reported for several fruit species, namely Punica granatum L. [72,73,74], Olea europaea L. [75,76], and Prunus dulcis (Mill.) D.A. Webb [77]. Increased production of antioxidants and secondary metabolites, including phenols, neutralizes reactive oxygen and radicals [78]. However, prolonged drought leads to a reduction in the overall content of the metabolites due to greater growth reductions [71]. This could explain low phenolic values in the southernmost populations, island Brač, Split, and Konavle, since these populations grow in the Mediterranean climate, with prolonged drought and intense heat during the summer months. On the other hand, S. domestica is a tree species that prefers the Mediterranean climate. Therefore, in the continental eco-geographical region, where service tree grows in sub-optimal sites, adverse conditions that are more likely to become extreme could result in a higher total phenolic content in the fruits from continental populations. These adverse conditions can also be notable in pedoclimate, i.e., soil type. The continental populations grow on shallow, extremely dry, and basic soils, with carbonate, mostly marl-type parent material. The Mediterranean populations, although also facing shallow and dry conditions, grow on terra rossa of flysch or dolomite parent material with clayey texture, enabling better water retention and distribution, despite the increased amount of soil skeleton.
Measured DPPH values did not follow the same trend as total phenolic content. For example, population Tounj was characterized by the highest total phenolic content values but lowest DPPH, whereas the values for Psunj were very similar. DPPH free radical scavenging is a method for the evaluation of the antioxidant activity of all compounds in a sample, e.g., plant tissue. Natural antioxidants in plants include polyphenols (phenolic acids, flavonoids, anthocyanins, lignans, and stilbenes), carotenoids (xanthophylls and carotenes), and vitamins (vitamin E and C) [79]. Research on S. domestica has largely been focused on polyphenols, thus providing only a partial picture of all the components with antioxidative abilities and their relative contents. Apart from polyphenols, the biochemical composition of the fruits of S. domestica is characterized by 4 to 22% of carotenoids, 22.3 to 98.3% of ascorbic acid, as well as 1.8% of vitamin C [80]. It is therefore possible to assume phenols compromise only a small proportion of antioxidative compounds in S. domestica fruits, which would explain the variable data and relations between total phenolic content and DPPH assays. Furthermore, FRAP has not been shown to have a significant correlation with geographical or environmental distribution and is overall not significant as a parameter of distinction.

3.4. Population Structure, Isolation by Distance and Environment

Our results provide new insights into the morphological and chemical diversity and structure of service tree populations from two eco-geographical regions. In both cases, the studied populations were dived into three groups. The first group encompasses populations from the continental eco-geographical region, the second group samples from the island population Brač, and the third group populations from the coastal Mediterranean area. In both cases, island population Brač was specific, compared to other populations included in this research. Overall, these findings are consistent with the well-known fact that island populations are distinct because they have a unique array of traits/genes when compared with the mainland populations [81,82]. Scattered distribution and geographical isolation between those populations have probably resulted in a limited long-distance gene flow. Under such conditions, genetic drift may play a very important role in shaping the structure of the populations [28,83,84,85]. This microevolutionary process is especially pronounced in such low-density populations. Similar results were reported for leaf morphometric analysis of service tree populations [29]. Island populations were highly divergent in leaflets shape in comparison with adjacent mainland populations. However, it is reported that in temperate, insect-pollinated, and animal-dispersed tree species such as S. domestica [17,86,87] and S. torminalis [88,89,90], fragmented subpopulations are functionally connected by gene flow through both pollen and seed. This can explain connectivity between other coastal populations since they were in both cases classified in the same group. Also, we cannot exclude the possibility of human-mediated gene-flow as this species has been cultivated since Roman times [8,9]. Nevertheless, the clear pattern of morphological and chemical diversity indicates a natural origin of the studied populations.
The ‘mixed origin’ of the Novi Vinodolski population (P04) is of particular interest because our data indicate a certain gene flow between this population and the populations in the continental eco-geographical region. The Northern Adriatic region, stretching from Novi Vinodolski to Senj, has been linked to Central Croatia throughout history by trade and the exchange of goods. It is, therefore, safe to assume human-mediated transplanting and exchange of plant material happened between the continental population of Tounj and the Mediterranean population Novi Vinodolski, which would explain the mixed origin of the Novi Vinodolski population.
The tests for isolation-by-distance were not significant for the studied populations, and the morphological structure between the Mediterranean and continental populations was largely explained by the environmental conditions and fits the isolation-by-environment pattern. In addition, we revealed that water, crude fat, and total phenolic contents of the fruits were associated with environmental variables. Overall, our results suggest that populations from the continental and Mediterranean eco-geographical regions, despite their admixture in some cases, are isolated by the environment. It is well known that environmental heterogeneity can have a fundamental impact on phenotypic diversity [91,92,93]. In particular, different environments can result in spatially variable selection pressures, thereby contributing to phenotypic divergence among populations via phenotypic plasticity or local adaptation [92,94]. Lower annual mean temperatures accompanied by high precipitation during the warmest quarter in the continental eco-geographical region resulted in larger fruits with higher water content. On the other hand, higher annual temperatures, especially during the warmest quarter, resulted in smaller fruits with lower water content and higher fat content from cuticular wax. Overall, our results indicate that multiple morpho- and chemo-types, each having higher fitness in its native habitat than the others, are probably the result of local adaptation [95].

4. Materials and Methods

4.1. Plant Material and Study Area

The study encompassed two populations from the continental eco-geographical region (P01—Psunj, P02—Tounj), and five populations from the Mediterranean eco-geographical region (P03—Istria, P04—Novi Vinodolski; P05—Split; P06—Brač; P07—Konavle). In each population, samples were collected from seven trees, which were then used in morphometric and proximate analyses as well as studies of total phenolics and antioxidant activity. The small number of sampled trees is a direct result of the small density of populations [6,8,9,29], i.e., a small number of mature trees spread across large areas in the habitat. This fragmentation made collecting challenging and caused a smaller overall sample size.
To describe the environmental differences of the studied populations, for the principal component (PC) analysis and the calculation of environmental distances, 19 bioclimatic variables were obtained for each collection site.

4.2. Morphometric Analysis

In order to conduct the morphological analysis, 50 fruits were collected from each tree. Firstly, fruit mass (m) was determined with 0.1 g precision weighing. After weighing, fruits were cut longitudinally, and the dimensions were measured automatically using the WinFolia software [96]. The accuracy of measurements was 0.1 mm, and the following morphological characteristics were measured: fruit length (FL); fruit length, measured from the fruit base to the point of maximum fruit width (PMFW); maximum fruit width (MFW); and fruit width at 15% (FW1) and 85% (FW2) of fruit length. After the fruits were measured, seeds were extracted, the pulp removed, and the seed dimensions, i.e., length (SL) and width (SW), were measured automatically using the WinSEEDLE software [97]. In addition, the number of seeds per fruit (NS) was counted manually. Finally, ten morphological characters were examined in order to assess the variation within and between populations.

4.3. Proximate Analysis

Water, crude protein, crude fat, ash, cellulose, and sugar contents were determined according to the procedures established by the Association of Official Analytical Chemists (AOAC). All analyses were performed in duplicate. Water content in the samples was determined by a physical, indirect method, in which a sample of known mass was dried in an air dryer (Instrumentaria, Zagreb, Croatia) at 105 °C until a constant weight was achieved (4 h) [98]. Total mineral content was measured as ash content, which represented inorganic fruit compounds left after the organic matter was combusted [99]. The 4 g fruit samples, previously carbonized on open gas burners, were ignited in a muffle furnace (Nabertherm GmbH, Lilienthal, Germany) at 580 °C, until a constant weight was achieved, i.e., until uniform, light-grey ash with no black admixtures formed. The Kjeldahl method was employed to determine the total nitrogen content, in combination with a copper catalyst using the block digestion system Foss Tecator 6—1007 Digestor (Foss Tecator, Höganäs, Sweden) and the Foss Kjeltec™ 8100 Auto Distillation unit (Foss Tecator, Höganäs, Sweden). Crude protein content was obtained by multiplying total nitrogen by a conversion factor of 6.25 [100]. Total crude fat extraction was performed by the Soxhlet apparatus (Inkolab d.o.o., Zagreb, Croatia); medical-grade petroleum was used for extraction during 5 h [101]. Crude cellulose was determined by the method of Kürschner and Hanak [102]. Sugar content in the fruits was determined according to the AOAC 925.35 [103] method. Reducing sugars, glucose and fructose were determined using Fehling solution A, which was reduced to copper (I) oxide by the sugars, under specific conditions. Non-reducing sugar, e.g., sucrose, was firstly hydrolyzed down to reducing sugars, whose content was determined using Fehling solutions and used to calculate the total sugar content (total inverted sugars).

4.4. Determination of Total Phenolics and Antioxidant Capacity

4.4.1. Fruit Extraction Procedure

Extraction was conducted according to Benvenuti et al. [104]. After the seeds were removed, fruits were finely chopped and 6 g (±0,01) were placed in a 250 mL Erlenmeyer flask. 20 mL of methanol (MeOH)/2% hydrochloric acid (HCl) (95:5) extraction solvent was added. The flask was sealed, wrapped in aluminum foil, and left for 1 h on the shaker. After 1 h on the shaker, extracts were filtered in a Büchner flask, using vacuum, and transferred into 50 mL volumetric flasks. The remaining extract in the Büchner flask was rinsed with a few mL of extraction solvent. Residue from the Büchner funnel was transferred into the same 250 mL Erlenmeyer flask and the procedure was repeated, resulting in a double volume of extract in the 50 mL volumetric flask. Flask was filled with extraction solvent up to the 50 mL mark and left refrigerated, at +4 °C, until further analysis. These extracts were used in determining the total phenolic content, the antioxidant capacity assay by DPPH method, as well as the reducing power assay using the FRAP method.

4.4.2. Estimation of Total Phenolic Content Using the Folin-Ciocalteu Method

Total phenolic content in the samples was determined according to the Folin-Ciocalteau method [105]. This colorimetric method is based on the reduction of the phosphotungstic-phosphomolybdenum acid complex, creating a blue chromophore. Maximum absorption of the chromophores depends on the alkaline solution and the concentration of phenolic compounds, with a more intense blue color indicating higher absorption and concentration of antioxidants [106].
In a tube, 40 µL of diluted extract sample, 3.16 mL of distilled water, and 0.2 mL FC reagent were added and vortexed. After 3 min, 0.6 mL of aqueous sodium carbonate (20%) were added, the sample was vortexed again and placed in a water bath at 40 °C. After 30 min, the absorbance was read, at 765 nm. Total phenol concentration was calculated from the calibration curve, where gallic acid was the standard. Results were expressed as mg of gallic acid per g of dried fruit matter (seeds excluded) [24].

4.4.3. Evaluation of Antioxidant Activity Using the DPPH Method

DPPH method is based on the reduction of stable radical 2,2-Diphenyil-picrylhydrazyl (DPPH*) in the presence of antioxidants. Due to its unpaired electron, DPPH* strongly absorbs in the spectrum of 517 nm [107]. Reduction causes stable, purple-colored DPPH* to change into yellow-colored DPPH-H [108,109,110]. Phenol compounds act as hydrogen donors, thus making this method suitable as an antioxidant assay.
A volume of 2 mL of diluted extract, 2 mL of methanol, and 1 mL of 0.5 mM DPPH methanolic solution were placed into a tube and shaken on a shaker. The blank was methanol without DPPH. Closed tubes were then kept in the dark, at room temperature, for 20 min. Absorbance was measured at a wavelength of 517 nm, against the blank of DPPH-free methanol. The DPPH* concentration in the reaction medium was calculated from the calibration curve for Trolox. Results were expressed as a percentage (%) of residual DPPH.

4.4.4. Determination of FRAP

Ferric reducing antioxidant power (FRAP) assay is a method in which yellow-colored Fe3+-2,4,6-tripyridyl-s-triazine (TPTZ) complex is reduced to the blue-colored ferrous form, in a low pH-value medium [109]. The FRAP reagent was prepared by mixing 50 mL of acetate buffer with 5 mL of TPTZ and 5 mL of FeCl3 (10:1:1). For sample analysis, 240 µL of distilled water was mixed with 80 µL of appropriately diluted sample and 80 µL of FRAP reagent. Blank was prepared in the same manner, substituting sample with extraction solution (methanol/HCl 2%). The mixture was left to stand for 5 min in a water bath, at 37 °C. Absorbance was measured at 595 nm. A calibration curve for FeSO4·7H2O was used to calculate FRAP values of the measured absorbances. Results were expressed as FRAP values, e.g., mmol Fe2+.

4.5. Statistical Analysis

Descriptive statistical parameters, arithmetic means, and coefficient of variation were calculated for the particular fruit morphometric characters and chemical composition for the studied populations. To assess the possibility of conducting multivariate statistical analyses and parametric tests, the symmetry, unimodality, and homoscedasticity of data were verified [111]. 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 studied populations and within populations. In addition, one-way ANOVA was used to test the differences in chemical composition between the studied populations. For each analysis, statistically significant differences between means were identified using the Fisher’s LSD multiple comparison test, at p ≤ 0.05. Descriptive statistics and ANOVA analyses were carried out using the program package STATISTICA [112].
To evaluate the correlation between multicharacter differences among populations, a Mantel test [113] was performed on the matrices of Euclidean distances. To assess isolation-by-distance and isolation-by-environment, morphological and chemical distance matrices were compared to the geographic and environmental distance matrices using the simple Mantel test [114,115]. To calculate the environmental distance matrix, climate data were obtained from the WorldClim 2 database with a spatial resolution close to a square km [116,117]. In addition, a simple Mantel test was performed between morphological and chemical distance matrices. The significance level was assessed after 10,000 permutations, and the Mantel test was performed with the R package “Vegan” [118].
To identify the divergence and structure of the studied populations, multivariate statistical methods were used [119]. Agglomerative hierarchical clustering algorithms were used to construct a tree diagram. Pairwise Euclidean distances were calculated, and cluster analysis was performed using the unweighted pair-group method with arithmetic mean (UPGMA). The K-means method was applied to detect structure and define the number of K-groups that best explained the morphological and chemical variation of populations (e.g., [49,120,121,122]). The eco-geographical structure of the studied populations was further analyzed using the Monmonier’ maximum difference algorithm, implemented in Barrier 2.2 software [123]. In addition, principal component (PC) analysis was used to calculate principal components across all individuals and all studied morphometric and chemical traits. The biplots were constructed by two principal components showing analyzed individuals and traits. To calculate the discriminatory power of characters among three groups of service tree populations established by the K-means clustering method, a canonical discriminant analysis was performed. The proportion of individuals correctly classified into the above-mentioned groups was determined using classificatory discriminant analyses [118,124]. 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 statistical program R v.3.2.2 [125].

5. Conclusions

Our study showed high levels of phenotypic and chemical diversity for service tree populations despite its scattered nature and low-density populations. We revealed the existence of three morphologically and chemically distinct and well-defined groups of service tree populations. We suggest that those differences are the result of both neutral and adaptive microevolution processes. With regard to the island population Brač, non-effective pollen and seed dispersal, along with genetic drift and specific environmental factors, resulted in a distinct phenotype with specific chemical characteristics. On the other hand, environmental heterogeneity, i.e., differences between the Mediterranean and continental eco-geographical regions, resulted in an adaptive phenotypic and chemical evolution. Overall, the patterns found in this study confirmed that multiple evolutionary processes influence the morphological and chemical diversity and structure of the populations.

Supplementary Materials

The following are available online at www.mdpi.com/article/10.3390/plants10081691/s1, Table S1: Environmental variables for seven Sorbus domestica populations. Acronyms of populations: P01—Psunj; P02—Tounj; P03—Istria; P04—Novi Vinodolski; P05—Split; P06—Brač; P07—Konavle. Environmental 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, Figure S1: Biplot of the principal component analysis based on 19 bioclimatic variables in seven studied Sorbus domestica populations. Acronyms of populations: P01 (Psunj), P02 (Tounj), P03 (Istria), P04 (Novi Vinodolski), P05 (Split), P06 (Brač), and P07 (Konavle). Acronyms for environmental variables as in Table S1, Table S2: Pearson’s correlation coefficients between 19 bioclimatic variables and scores of the first three principal components. Acronyms for environmental variables as in Table S1, Figure S2: Sorbus domestica fruit samples from seven studied populations, Table S3: Pearson’s correlation coefficients between ten morphological traits and scores of the first three principal components, Table S4: Pearson’s correlation coefficients between ten chemical traits and scores of the first three principal components, Table S5: Results of the stepwise discriminant analyses for morphometric traits. p(λ), significance of Wilks’λ: *** significant at p < 0.001, ** significant at 0.001 < p < 0.01, * significant at 0.01 < p < 0.05, ns non-significant values (p > 0.05), Table S6: Results of the stepwise discriminant analyses for chemical traits. p(λ), significance of Wilks’λ: *** significant at p < 0.001, ** significant at 0.001 < p < 0.01, * significant at 0.01 < p < 0.05, ns non-significant values (p > 0.05).

Author Contributions

Conceptualization: I.P., N.V. and M.I.; Methodology, I.P. and N.V.; Software, I.P.; Validation, Z.L. and D.K.; Morphometric Analysis, I.Z. and A.V.; Proximate Analysis, I.P. and N.V.; Total Phenolics and Antioxidant Analysis, N.V., V.P. and B.V.; Statistical Analysis, I.P.; Resources, I.P., N.V. and M.I.; Data Curation, I.P., N.V. and Z.L.; Writing—Original Draft Preparation, I.P., N.V. and K.T.; Writing—Review & Editing, Z.L., I.Z., A.V., V.P., D.K. and M.I.; Visualization, I.P. 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.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

We would like to thank Ivan Ljubić, for his help in collecting samples in the area of Split, and Ivan Andrić, for his help in preparing the seed samples for morphometric analysis. Furthermore, we would like to thank the technicians, Valentina Hohnjec and Renata Petrović, from the Laboratory for Quality Control in the Food Industry of the Department of Food Quality Control, Faculty of Food Technology and Biotechnology of the University of Zagreb, for their assistance in chemical analyses.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

References

  1. Drvodelić, D.; Oršanić, M.; Grahovac-Tremski, M. Rasadnička proizvodnja šumskih voćkarica u rasadnicima Hrvatskih šuma doo. Šumar. List 2020, 11–12, 597–606. [Google Scholar] [CrossRef]
  2. Sardeshpande, M.; Shackleton, C. Wild edible fruits: A systematic review of an under-researched multifunctional NTFP (Non-Timber Forest Products). Forests 2019, 10, 467. [Google Scholar] [CrossRef] [Green Version]
  3. Bvenura, C.; Sivakumar, D. The role of wild fruits and vegetables in delivering a balanced and healthy diet. Food Res. Int. 2017, 99, 15–30. [Google Scholar] [CrossRef]
  4. Stryamets, N. Non-wood forest products for livelihoods. Bosque 2012, 33, 329–332. [Google Scholar]
  5. Sheppard, J.P.; Chamberlain, J.; Agundez, D.; Bhattacharya, P.; Chirwa, P.W.; Gontcharov, A.; Sagona, W.C.J.; Shen, H.; Tadesse, W.; Mutke, S. Sustainable forest management beyond the timber-oriented Status Quo: Transitioning to coproduction of timber and non-wood forest products—A global perspective. Curr. For. Rep. 2020, 6, 26–40. [Google Scholar] [CrossRef] [Green Version]
  6. Rotach, P. EUFORGEN Technical Guidelines for Genetic Conservation and Use for Service Tree (Sorbus domestica); International Plant Genetic Resources Institute: Rome, Italy, 2003; p. 6. [Google Scholar]
  7. Brütsch, U.; Rotach, P. Der Speierling (Sorbus domestica L.) in der Schweiz: Verbreitung, ökologie, standsortsanspruche, konkurrenzkraft und waldbauliche eignung. Schweiz. Z. Forstwes. 1993, 144, 967–991. [Google Scholar]
  8. George, J.-P.; Konrad, H.; Collin, E.; Thevenet, J.; Ballian, D.; Idžojtić, M.; Kamm, U.; Zhelev, P.; Geburek, T. High molecular diversity in the true service tree (Sorbus domestica) despite rareness: Data from Europe with special reference to the Austrian occurrence. Ann. Bot. 2015, 115, 1105–1115. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  9. George, J.-P.; Woodman, J.; Hampton, M.; Konrad, H.; Geburek, T. True service-tree (Sorbus domestica, Rosaceae) in the British Isles: Rare but diverse. New J. Bot. 2016, 6, 21–30. [Google Scholar] [CrossRef]
  10. Matić, S.; Vukelić, J. Speierling und Elsbeere in den Waldern Kroatiens. Corminaria 2001, 16, 31–33. [Google Scholar]
  11. Paganová, V. Ecology and distribution of Service tree Sorbus domestica L. in Slovakia. Ekológia 2008, 27, 152–167. [Google Scholar]
  12. Poljak, I.; Idžojtić, M.; Zebec, M.; Šapić, I. Diversity and distribution of Sorbus spp. in Croatia. In Book of Abstracts of the 34th International Symposium—Eastern Alpine and Dinaric Society for Vegetation Ecology, Proceedings of International symposium—Eastern Alpine and Dinaric Society for Vegetation Ecology, Camerino, Italy, 24–28 May 2011; Carimini, L., Ed.; University of Camerino: Camerino, Italy, 2011. [Google Scholar]
  13. Enescu, C.M.; de Rigo, D.; Houston Durrant, T.; Caudullo, G. Sorbus domestica in Europe: Distribution, habitat, usage and threats. In European Atlas of Forest Tree Species; San-Miguel-Ayanz, J., de Rigo, D., Caudullo, G., Houston Durrant, T., Mauri, A., Eds.; Publication Office of the European Union: Luxembourg, 2016; p. e019db5+. [Google Scholar]
  14. Jeran, N.; Židovec, V.; Dujmović Purgar, D.; Vokurka, A.; Mešić, A.; Karlović, K.; Duralija, B. Utilization value of Sorbus s.l. in Croatia and its potential enhancement. Acta Hortic. 2021, 1308, 81–88. [Google Scholar] [CrossRef]
  15. Pietzarka, U.; Heidecke, C.; Roloff, A. Sorbus domestica L. In Enzyklopädie der Holzgewächse; Roloff, A., Weisgeber, H., Lang, U.M., Stimm, B., Eds.; Wiley-VCH: Weinheim, Germany, 2010; pp. 1–16. [Google Scholar]
  16. Idžojtić, M. Dendrology: Cones, Flowers, Fruits and Seeds; Elsevier—Academic Press: London, UK; San Diego, CA, USA; Cambridge, MA, USA; London, 2019; 800p. [Google Scholar]
  17. Kamm, U.; Rotach, P.; Gugerli, F.; Siroky, M.; Edwards, P.; Holderegger, R. Frequent long-distance gene flow in a rare temperate forest tree (Sorbus domestica) at the landscape scale. Heredity 2009, 103, 476–482. [Google Scholar] [CrossRef] [Green Version]
  18. Miko, M.; Gažo, J. Morphological diversity of Sorbus domestica at the level of fruits and leaves in the selected localities of Slovakia. Biologia 2003, 58, 35–39. [Google Scholar]
  19. Miko, M.; Gažo, J. Morphological and biological characteristics of fruits and seed of the service tree (Sorbus domestica L.). J. Fruit Ornam. Plant Res. 2004, 12, 139–146. [Google Scholar]
  20. Mikić, T.; Orlović, S.; Marković, M.; Kovačević, B.; Pilipović, A. Variability in service tree (Sorbus domestica L.) populations in Serbia. Lesn. Čas.–For. J. 2008, 54, 61–67. [Google Scholar]
  21. Miletić, R.; Paunović, S.M. Research into service tree (Sorbus domestica L.) population in Eastern Serbia. Genetika 2012, 44, 483–490. [Google Scholar] [CrossRef]
  22. Drvodelić, D.; Oršanić, M.; Vuković, M.; Jatoi, M.A.; Jemrić, T. Correlation of fruit size with morphophysiological properties and germination rate of the seeds of Service tree (Sorbus domestica L.). SEEFOR 2018, 9, 47–54. [Google Scholar] [CrossRef] [Green Version]
  23. Ölschläger, C.; Milde, J.; Schempp, H.; Treuter, D. Polyphenols and antioxidant capacity of Sorbus domestica L. fruits. J. Appl. Bot. Food. Qual. 2004, 78, 112–116. [Google Scholar]
  24. Termentzi, A.; Kefalas, P.; Kokkalou, E. Antioxidant activities of various extracts and fractions of Sorbus domestica fruits at different maturity stages. Food Chem. 2006, 98, 599–608. [Google Scholar] [CrossRef]
  25. Termentzi, A.; Kefalas, P.; Kokkalou, E. LC–DAD–MS (ESI+) analysis of the phenolic content of Sorbus domestica fruits in relation to their maturity stage. Food Chem. 2008, 106, 1234–1245. [Google Scholar] [CrossRef]
  26. Piagnani, M.C.; Debellini, C.; LoScalzo, R. Ohyllometry and carpometry, chemical and functional characterization of fruits of Sorbus domestica L. (service tree) selections. J. Berry Res. 2012, 2, 7–22. [Google Scholar] [CrossRef] [Green Version]
  27. Majić, B.; Šola, I.; Likić, S.; Juranović Cindrić, I.; Rusak, G. Characterisation of Sorbus domestica L. bark, fruit and seeds: Nutrient composition and antioxidant activity. Food. Technol. Biotechnol. 2015, 53, 463–471. [Google Scholar] [CrossRef] [PubMed]
  28. Brus, R.; Ballian, D.; Bogunić, F.; Bobinac, M.; Idžojtić, M. Leaflet morphometric variation of service tree (Sorbus domestica L.) in the Balkan Peninsula. Plant Biosyst. 2011, 145, 278–285. [Google Scholar] [CrossRef]
  29. Poljak, I.; Kajba, D.; Ljubić, I.; Idžojtić, M. Morphological variability of leaves of Sorbus domestica L.in Croatia. Acta Soc. Bot. Pol. 2015, 84, 249–259. [Google Scholar] [CrossRef] [Green Version]
  30. Rutkowska, M.; Olszewska, M.A.; Kolodziejczyk-Czepas, J.; Nowak, P.; Owczarek, A. Sorbus domestica leaf extracts and their activity markers: Antioxidant potential and synergy effects in scavenging assays of multiple oxidants. Molecules 2019, 24, 2289. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  31. Bednorz, L. Morphological variability of leaves of Sorbus torminalis (L.) Crantz in Poland. Acta Soc. Bot. Pol. 2006, 75, 233–243. [Google Scholar] [CrossRef] [Green Version]
  32. Bednorz, L. Morphological variability of fruits and seeds of Sorbus torminalis in Poland. Dendrobiology 2007, 57, 3–14. [Google Scholar]
  33. Paganová, V. Taxonomic reliability of leaf and fruit morphological characteristics of the Pyrus L. taxa in Slovakia. Hort. Sci. 2003, 30, 98–107. [Google Scholar] [CrossRef] [Green Version]
  34. Paganová, V. The occurrence and morphological characteristics of the wild pear lower taxa in Slovakia. Hort. Sci. 2009, 36, 1–13. [Google Scholar] [CrossRef] [Green Version]
  35. Reim, S.; Proft, A.; Heinz, S.; Höfer, M. Diversity of the European indigenous wild apple Malus sylvestris (L.) Mill. in the East Ore Mountains (Osterzgebirge), Germany: I. Morphological characterization. Genet. Resour. Crop Evol. 2012, 59, 1101–1114. [Google Scholar] [CrossRef]
  36. Rakonjac, V.; Mratinić, E.; Jovković, R.; Fotirić Akšić, M. Analysis of morphological variability in wild cherry (Prunus avium L.) genetic resources from Central Serbia. J. Agric. Sci. Tech. 2014, 16, 151–162. [Google Scholar]
  37. Khadivi, A.; Mohammadi, M.; Asgari, K. Morphological and pomological characterizations of sweet cherry (Prunus avium L.), sour cherry (Prunus cerasus L.) and duke cherry (Prunus × gondouinii Rehd.) to choose the promising selections. Sci. Hort. 2019, 257, 108719. [Google Scholar] [CrossRef]
  38. Albarouki, E.; Peterson, A. Molecular and morphological characterization of Crataegus L. species (Rosaceae) in southern Syria. Bot. J. Linn. Soc. 2007, 153, 255–263. [Google Scholar] [CrossRef]
  39. Yanar, M.; Ercisli, S.; Yilmaz, K.U.; Sahiner, H.; Taskin, T.; Zengin, Y.; Akgul, I.; Celik, F. Morphological and chemical diversity among hawthorn (Crataegus spp.) genotypes from Turkey. SRE 2011, 6, 35–38. [Google Scholar]
  40. Krüssmann, G. Handbuch der Laubgehölze; Paul Parkey in Berlin und Hamburg, Verlag für Landwirtschaft, Veterinärmedizin, Gartenbau und Forstwesen: Berlin, Germany, 1962; Volume 2. [Google Scholar]
  41. Hegi, G. Illustrierte Flora von Mitteleuropa; Blackwell Wissenschaft-Verlag: Berlin, Germany, 1995; pp. 353–356. [Google Scholar]
  42. Barengo, N.; Rudow, A.; Schwab, P. Förderung seltener Baumarten auf der Schweizer Alpennordseite; Buwal/Eidgenössische Forstdirektion and ETHZ/Professur Waldbau: Bern, Switzerland, 2001. [Google Scholar]
  43. Brus, R. Dendrologija za Gozdarje (Dendrology for Foresters), 2nd ed.; Biotechnical Faculty, University of Ljubljana: Ljubljana, Slovenia, 2008. [Google Scholar]
  44. Bignami, C. Service tree (Sorbus domestica L.). Description and use of service tree. Inf. Agrar. 2000, 56, 55–58. [Google Scholar]
  45. Ferrazzini, D.; Monteleone, I.; Belletti, P. Small-scale genetic diversity in oneseed hawthorn (Crataegus monogyna Jacq.). Eur. J. For. Res. 2008, 127, 407–414. [Google Scholar] [CrossRef]
  46. Rahmani, M.-S.; Shabanian, N.; Khadivi-Khub, A.; Woeste, K.E.; Badakhshan, H.; Alikhani, L. Population structure and genotypic variation of Crataegus pontica inferred by molecular markers. Gene 2015, 572, 123–129. [Google Scholar] [CrossRef]
  47. Popović, V.; Kerkez, I. Varijabilnost populacija divlje trešnje (Prunus avium L.) u Srbiji prema morfološkim svojstvima listova. Šumar. List 2016, 140, 347–355. [Google Scholar] [CrossRef] [Green Version]
  48. Popović, V.; Lučić, A.; Kerkez Janković, I.; Rakonjac, L.; Bogdan, S. Varijabilnost svojstava plodova kod provenijencija divlje trešnje (Prunus avium L.) u Srbiji. Šumar. List 2020, 144, 585–596. [Google Scholar] [CrossRef]
  49. Poljak, I.; Idžojtić, M.; Šapić, I.; Korijan, P.; Vukelić, J. 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. Šumar. List 2018, 142, 19–32. [Google Scholar]
  50. Frankham, R. Do island populations have less genetic variation than mainland populations? Heredity 1997, 78, 311–327. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  51. García-Verdugo, C.; Sajeva, M.; La Mantia, T.; Harrouni, C.; Msanda, F.; Caujapé-Castells, J. Do island plant populations really have lower genetic variation than mainland populations? Effects of selection and distribution range on genetic diversity estimates. Mol. Ecol. 2015, 24, 726–741. [Google Scholar] [CrossRef] [PubMed]
  52. Morente-López, J.; García, C.; Lara-Romero, C.; García-Fernández, A.; Draper, D.; María Iriondo, J. Geography and environment shape landscape genetics of Mediterranean Alpine species Silene ciliata Poiret. (Caryophyllaceae). Front. Plant. Sci. 2018, 9, 1698. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  53. Dowle, E.J.; Morgan-Richards, M.; Brescia, F.; Trewick, S.A. Correlation between shell phenotype and local environment suggests a role for natural selection in the evolution of Placostylus snails. Mol. Ecol. 2015, 24, 4205–4221. [Google Scholar] [CrossRef]
  54. Souci, S.W.; Fachmann, W.; Kraut, H. Die Zusammensetzung der Lebensmittel, Nährwert-Tabellen; Wissenschaftliche Verlagsgesellschaft MBH: Stuttgart, Germany, 1979. [Google Scholar]
  55. Raspé, O.; Findlay, C.; Jacquemart, A.-L. Sorbus aucuparia L. J. Ecol. 2000, 88, 910–930. [Google Scholar]
  56. Drvodelić, D.; Jemrić, T.; Oršanić, M. Jarebika (Sorbus aucuparia L.): Važnost, Uzgoj i Uporaba; University of Zagreb, Faculty of Forestry: Zagreb, Croatia, 2019. [Google Scholar]
  57. Xue, D.; Zhang, X.; Lu, X.; Chen, G.; Chen, Z. Molecular and evolutionary mechanisms of cuticular wax for plant drought tolerance. Front. Plant Sci. 2017, 8, 621. [Google Scholar] [CrossRef] [PubMed]
  58. Von Kausch-Blecken Schmeling, W. Der Speierling, 2nd ed.; Druckhaus Gottingen im Gottinger Tageblatt GmbH & Co.: Göttinger, Germany, 2000. [Google Scholar]
  59. Sarv, V.; Venskutonis, P.R.; Bhat, R. The Sorbus spp.—Underutilised plants for foods and nutraceuticals: Review on polyphenolic phytochemicals and antioxidant potential. Antioxidants 2020, 9, 813. [Google Scholar] [CrossRef]
  60. Nielsen, E.; Temporiti, M.E.E.; Cella, R. Improvement of phytochemical production by plant cells and organ culture and by genetic engineering. Plant Cell. Rep. 2019, 38, 1199–1215. [Google Scholar] [CrossRef]
  61. Oh, M.-M.; Carey, E.E.; Rajashekar, C.B. Environmental stresses induce health-promoting phytochemicals in lettuce. Plant Physiol. Biochem. 2009, 47, 578–583. [Google Scholar] [CrossRef]
  62. Bian, Z.H.; Yang, Q.C.; Liu, W.K. Effects of light quality on the accumulation of phytochemicals in vegetables produced in controlled environments: A review. J. Sci. Food Agric. 2015, 95, 869–877. [Google Scholar] [CrossRef]
  63. Liebelt, D.J.; Jordan, J.T.; Doherty, C.J. Only a matter of time: The impact of daily and seasonal rhythms on phytochemicals. Phytochem. Rev. 2019, 18, 1409–1433. [Google Scholar] [CrossRef]
  64. Howard, L.R.; Talcott, S.T.; Brenes, C.H.; Villalon, B. Changes in phytochemical and antioxidant activity of selected Pepper cultivars (Capsicum species) as influenced by maturity. J. Agric. Food. Chem. 2000, 48, 1713–1720. [Google Scholar] [CrossRef] [PubMed]
  65. Seguin, P.; Zheng, W.; Souleimanov, A. Alfalfa phytoestrogen content: Impact of plant maturity and herbage components. J. Agron. Crop. Sci. 2004, 190, 211–217. [Google Scholar] [CrossRef]
  66. Figueroa-Pérez, M.G.; Pérez-Ramírez, I.F.; Paredes-López, O.; Mondragón-Jacobo, C.; Reynoso-Camacho, R. Phytochemical composition and in vitro analysis of Nopal (O. Ficus-Indica) cladodes at different stages of maturity. Int. J. Food Prop. 2018, 21, 1728–1742. [Google Scholar] [CrossRef] [Green Version]
  67. Gobbo-Neto, L.; Lopes, N.P. Medicinal plants: Factors of influence on the content of secondary metabolites. Quim. Nova 2007, 30, 374–381. [Google Scholar] [CrossRef]
  68. Borges, L.L.; Alves, S.F.; Sampaio, B.L.; Conceição, E.C.; Bara, M.T.F.; Paula, J.R. Environmental factors affecting the concentration of phenolic compounds in Myrcia tomentosa leaves. Rev. Bras. Pharmacogn. 2013, 23, 230–238. [Google Scholar] [CrossRef] [Green Version]
  69. Rivero, R.M.; Ruiz, J.M.; Romero, L. Oxidative metabolism in tomato plants subjected to heat stress. J. Hortic. Sci. Biotech. 2004, 79, 560–564. [Google Scholar] [CrossRef]
  70. Sanchez-Rodriguez, E.; Rubio-Wilhelmi, M.; Cervilla, L.M.; Blasco, B.; Rios, J.J.; Rosales, M.A. Genotypic differences in some physiological parameters symptomatic for oxidative stress under moderate drought in tomato plants. Plant Sci. 2010, 178, 30–40. [Google Scholar] [CrossRef]
  71. Toscano, S.; Trivellini, A.; Cocetta, G.; Bulgari, R.; Francini, A.; Romano, D.; Ferrante, D. Effect of preharvest abiotic stresses on the accumulation of bioactive compounds in horticultural produce. Front. Plant Sci. 2019, 10, 1212. [Google Scholar] [CrossRef] [Green Version]
  72. Hadi Rad, M.; Reza Asghari, M.; Hassan Asareh, M. Morphological, chemical and biochemical response to drought stress in pomegranate (Punica granatum L.) fruit. IJBPAS 2015, 4, 64–78. [Google Scholar]
  73. Hassani Moghaddam, E.; Esna-Ashari, M.; Shaaban, M. An investigation of the secondary metabolites and antioxidant capacity of some commercial Iranian pomegranate (Punica granatum L.) cultivars under drought stress. Herb. Med. J. 2018, 3, 14–25. [Google Scholar]
  74. Faraji, S.; Hadadinejad, M.; Abdosii, V.; Basaki, T.; Karami, S. Screening pomegranate (Punica granatum L.) genotypes for drought tolerance using physiological and phytochemical characteristics. Fruits 2020, 75, 130–140. [Google Scholar] [CrossRef]
  75. Brito, C.; Dinis, L.-T.; Moutinho-Pereira, J.; Correia, C.M. Drought stress effects and Olive tree acclimation under a changing climate. Plants 2019, 8, 232. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  76. Jiménez-Herrera, R.; Pacheco-López, B.; Peragón, J. Water stress, irrigation and concentrations of pentacyclic triterpenes and phenols in Olea europaea L. cv. Picual Olive trees. Antioxidants 2019, 8, 294. [Google Scholar] [CrossRef] [Green Version]
  77. Lipan, L.; Cano-Lamadrid, M.-C.; Hernández, F.; Sendra, E.; Corell, M.; Vázquez-Araújo, L.; Moriana, A.; Carbonell-Barrachina, A.A. Long-term correlation between water deficit and quality markers in hydroSOStainable almonds. Agronomy 2020, 20, 1470. [Google Scholar] [CrossRef]
  78. Weidner, S.; Karolak, M.; Kramarać, M.; Kosinska, A.; Amarowicz, R. Phenolic compounds and properties of antioxidants in grapevine roots (Vitis vinifera L.) under drought stress followed by recovery. Acta Soc. Bot. Pol. 2009, 78, 97–103. [Google Scholar] [CrossRef] [Green Version]
  79. Dong-Ping, X.; Meng, X.; Zhou, T.; Zhou, Y.; Zheng, J.; Zhang, J.; Li, H. Natural antioxidants in foods and medicinal extractions, assessment and resources. Int. J. Mol. Sci. 2017, 18, 96. [Google Scholar]
  80. Drvodelić, D.; Jemrić, T.; Oršanić, M. Oskoruša: Važnost, Uporaba i Uzgoj; University of Zagreb, Faculty of Forestry: Zagreb, Croatia, 2015. [Google Scholar]
  81. Maki, M.; Yamashiro, T.; Matsumura, S. High levels of genetic diversity in island populations of the island endemic Suzukia luchuensis (Labiatae). Heredity 2003, 91, 300–305. [Google Scholar] [CrossRef]
  82. Hufford, K.M.; Mazer, S.J.; Hodges, S.A. Genetic variation among mainland and island populations of a native perennial grass used in restoration. AoB PLANTS 2014, 6, plt055. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  83. Wright, S. Isolation by distance under diverse systems of mating. Genetics 1946, 31, 39–59. [Google Scholar] [CrossRef]
  84. Savolainen, O.; Kuittinen, H. Small Population Processes. In Forest Conservation Genetics—Principles and Practice; Young, A., Boshier, D., Boyle, T., Eds.; CABI Publishing: Wallingford, UK, 2000; pp. 91–100. [Google Scholar]
  85. Arnaud-Haond, S.; Teixeira, S.; Massa, S.I.; Billot, C.; Saenger, P.; Coupland, G.; Duarte, C.M.; Serrão, A. Genetic structure at range edge: Low diversity and high inbreeding in Southeast Asian mangorve (Avicennia marina) populations. Mol. Ecol. 2006, 15, 3515–3525. [Google Scholar] [CrossRef]
  86. Kamm, U.; Gugerli, F.; Rotach, P.; Edwards, P.; Holderegger, R. Open areas in a landscape enhance pollen-mediated gene flow of a tree species: Evidence from northern Switzerland. Landsc. Ecol. 2010, 25, 903–911. [Google Scholar] [CrossRef] [Green Version]
  87. Kamm, U.; Gugerli, F.; Rotach, P.; Edwards, P.; Holderegger, R. Genetic consequences of local tree densities—Implications for the management of naturally rare, insect pollinated species in temperate forests. Forest Ecol. Manag. 2011, 262, 1047–1053. [Google Scholar] [CrossRef]
  88. Oddou-Muratorio, S.; Petit, R.J.; Le Guerroue, B.; Guesnet, D.; Demesure, B. Pollen- versus seed-mediated gene flow in a scattered forest tree species. Evolution 2001, 55, 1123–1135. [Google Scholar] [CrossRef]
  89. Oddou-Muratorio, S.; Klein, E.K.; Austerlitz, F. Pollen-flow in the wild service tree, Sorbus torminalis (L.) Crantz. II. Pollen dispersal and heterogeneity in mating success inferred from parent-offspring analysis. Mol. Ecol. 2005, 14, 4441–4452. [Google Scholar] [CrossRef] [PubMed]
  90. Oddou-Muratorio, S.; Klein, E.K.; Demesure-Musch, B.; Austerlitz, F. Real-time patterns of pollen flow in the wild service tree, Sorbus torminalis (Rosaceae). III. Mating patterns and the ecological maternal neighbourhood. Am. J. Bot. 2006, 93, 1650–1659. [Google Scholar] [CrossRef] [PubMed]
  91. Xu, F.; Guo, W.; Xu, W.; Wei, V.; 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]
  92. Vitasse, Y.; Lenz, A.; Kollas, C.; Randin, C.F.; Hoch, G.; Korner, C. Genetic vs. non-genetic responses of leaf morphology and growth to elevation in temperate tree species. Funct. Ecol. 2013, 28, 243–252. [Google Scholar] [CrossRef] [Green Version]
  93. Sander, N.L.; da Silva, C.J.; Duarte, A.V.M.; Zago, B.W.; Galbiati, C.; Viana, I.G.; de Arruda, J.C.; Dardengo, J.E.; Poletine, J.P.; Siqueira Leite, M.H.; et al. The influence of environmental features on the morphometric variation in Mauritia flexuosa L.f. fruits and seeds. Plants 2020, 9, 1304. [Google Scholar] [CrossRef]
  94. Dubuc-Messier, G.; Réale, D.; Perret, P.; Charmantier, A. Environmental heterogeneity and population differences in blue tits personality traits. Behav. Ecol. 2017, 28, 448–459. [Google Scholar] [CrossRef] [Green Version]
  95. DeWoody, J.; Trewin, H.; Taylor, G. Genetic and morphological differentiation in Populus nigra L.: Isolation by colonization or isolation by adaptation? Mol. Ecol. 2015, 24, 2641–2655. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  96. WinFolia, T.M. WinFOLIA Software Program; Version PRO 2005b; Regent Instruments Inc.: Quebec, QC, Canada, 2001. [Google Scholar]
  97. WinSeedle, T.M. For Seed and Conifer Needles Analysis; Version 2011a; Regent Instruments Inc.: Quebec, QC, Canada, 2010. [Google Scholar]
  98. AOAC International. Fruits and Fruit Products—Moisture in Fruits and Fruit Products; AOAC Official method 925.03; AOAC International: Gaithersburg, MD, USA, 2000. [Google Scholar]
  99. AOAC International. Fruits and Fruit Products—Ash Content in Fruits and Fruit Products; AOAC Official method 923.03; AOAC International: Gaithersburg, MD, USA, 2000. [Google Scholar]
  100. AOAC International. Fruits and Fruit Products—Crude Proteins in Fruits and Fruit Products; AOAC Official method 992.15; AOAC International: Gaithersburg, MD, USA, 2000. [Google Scholar]
  101. AOAC International. Fruits and Fruit Products—Total Fat in Fruits and Fruit Products; AOAC Official method 989.05; AOAC International: Gaithersburg, MD, USA, 2000. [Google Scholar]
  102. AOAC International. Fruits and Fruit Products—Cellulose in Fruits and Fruit Products; AOAC Official method 973.18; AOAC International: Gaithersburg, MD, USA, 2000. [Google Scholar]
  103. AOAC International. Fruits and Fruit Products—Sucrose in Fruits and Fruit Products; AOAC Official method 925.35; AOAC International: Gaithersburg, MD, USA, 2000. [Google Scholar]
  104. Benvenuti, S.; Pellati, F.; Melegari, M.; Bertelli, D. Polyphenols, anthocyanins, ascorbic acid and radical scavenging activity of Rubus, Ribes and Aronia. JFS Food Chem. Toxycol. 2004, 69, 164–169. [Google Scholar] [CrossRef]
  105. Waterman, P.G.; Mole, S. Analysis of Phenolic Plant Metabolites; Methods in Ecology; Blackwell Scientific Publications: Hoboken, NJ, USA; Oxford, UK, 1994. [Google Scholar]
  106. Stratil, P.; Klejdus, B.; Kuban, V. Determination of total content of phenolic compounds and their antioxidant activity in vegetables: Evaluation of spectrophotometric methods. J. Agric. Food Chem. 2006, 54, 607–616. [Google Scholar] [CrossRef] [PubMed]
  107. Kazazić, S.P. Antioksidacijska i antiradiklaska aktivnost flavonoida. Arh. Hig. Rada Toksikol. 2004, 55, 279–290. [Google Scholar] [PubMed]
  108. Brand-Williams, W.; Cuvelier, M.E.; Berset, C. Use of a free radical method to evaluate antioxidant activity. LWT 1995, 28, 25–30. [Google Scholar] [CrossRef]
  109. Tolić, M.T.; Jurčević, I.L.; Panjkota Krbavčić, I.; Marković, K.; Vahčić, N. Antioxidant properties of Chokeberry products. Food Technol. Biotechnol. 2015, 53, 171–179. [Google Scholar] [PubMed]
  110. Mrkonjić, Z. Fitohemijska Karakterizacija i Biohemijska Ispitivanja Plodova Vrsta Roda Sorbus L. 1753 (Rosaceae, Maloideae) Kao Izvora Prirodnih Nutraceutika. Ph.D. Thesis, Faculty of Science, University of Novi Sad, Novi Sad, Serbia, 2017. [Google Scholar]
  111. Sokal, R.R.; Rohlf, F.J. Biometry; Freeman and Co.: San Francisco, CA, USA, 1989. [Google Scholar]
  112. StatSoft, Inc. STATISTICA (Data Analysis Software System); Version 8.0; StatSoft, Inc.: Tulsa, OK, USA, 2001. [Google Scholar]
  113. Mantel, N. The detection of disease clustering and a generalized regression approach. Cancer Res. 1967, 27, 209–220. [Google Scholar]
  114. Smouse, P.E.; Long, J.C.; Sokal, R. Multiple regression and correlation extensions of the Mantel test of matrix correspondence. Syst. Zool. 1986, 35, 627–632. [Google Scholar] [CrossRef]
  115. Manly, B.F.J. Randomization, Bootstrap and Monte Carlo Methods in Biology, 3rd ed.; Chapman & Hall/CRC, Taylor & Francis Group: Boca Raton, FL, USA, 2007; p. 480. [Google Scholar]
  116. Hijmans, R.J.; Cameron, S.E.; Parra, J.L.; Jones, P.G.; Jarvis, A. Very high resolution interpolated climate surfaces for global land areas. Int. J. Climatol. 2005, 25, 1965–1978. [Google Scholar] [CrossRef]
  117. Fick, S.E.; Hijmans, R.J. WorldClim 2: New 1km spatial resolution climate surfaces for global land areas. Int. J. Climatol. 2017, 37, 4302–4315. [Google Scholar] [CrossRef]
  118. Oksanen, J.; Blanchet, F.G.; Friendly, M.; Kindt, R.; Legendre, P.; McGlinn, D.; Minchin, P.R.; O’Hara, R.B.; Simpson, G.L.; Solymos, P.; et al. Package ‘Vegan’: Community Ecology Package. Version 2.4-2. 2017. Available online: https://github.com/vegandevs/vegan (accessed on 24 June 2021).
  119. McGarigal, K.; Cushman, S.; Stafford, S. Multivariate Statistics for Wildlife and Ecology Research; Springer Verlag: New York, NY, USA, 2000; 283p. [Google Scholar]
  120. Douaihy, B.; Sobierajska, K.; Jasińska, A.K.; Boratyńska, K.; Ok, T.; Romo, A.; Machon, N.; Didukh, Y.; Dagher-Kharrat, M.B.; Boratyński, A. Morphological versus molecular markers to describe variability in Juniperus excelsa subsp. excelsa (Cupressaceae). AoB Plants 2012, 2012, pls003. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  121. 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] [Green Version]
  122. Sobierajska, K.; Boratyńska, K.; Jasińska, A.; Dering, M.; Ok, T.; Douaihy, B.; Dagher-Kharrat, M.B.; Romo, A.; Boratyński, A. Effect of the Aegean Sea barrier between Europe and Asia on differentiation in Juniperus drupacea (Cupressaceae). Bot. J. Linn. Soc. 2016, 180, 365–385. [Google Scholar] [CrossRef] [Green Version]
  123. Manni, F.; Guerard, E.; Heyer, E. Geographic patterns of (genetic, morphologic, linguistic) variation: How barriers can be detected by ‘Monmonier’s algorithm’. Hum. Biol. 2004, 76, 173–190. [Google Scholar] [CrossRef] [PubMed]
  124. Koutecky, P. MorphoTools: A set of R functions for morphometric analysis. Plant Syst. Evol. 2015, 301, 1115–1121. [Google Scholar] [CrossRef]
  125. R Core Team. R: A Language and Environment for Statistics in Biological Research, 4th ed.; W.H. Freeman and Co.: New York, NY, USA, 2012; p. 973. [Google Scholar]
Figure 1. Results of the multivariate statistical methods for fruit and seed morphometric analysis and locations of the seven sampled Sorbus domestica populations. (A) Geographical distribution of three groups of populations detected from K-means clustering (the proportions of the ancestry of each population in each of the defined clusters are color-coded: cluster A– green, cluster B–yellow, cluster C–blue); (B) Hierarchical tree dendrogram; (C) Barplot with posterior probabilities of classification of each individual into each group from the results of the classification discriminant analysis. Acronyms of populations: P01—Psunj; P02—Tounj; P03—Istria; P04—Novi Vinodolski; P05—Split; P06—Brač; P07—Konavle.
Figure 1. Results of the multivariate statistical methods for fruit and seed morphometric analysis and locations of the seven sampled Sorbus domestica populations. (A) Geographical distribution of three groups of populations detected from K-means clustering (the proportions of the ancestry of each population in each of the defined clusters are color-coded: cluster A– green, cluster B–yellow, cluster C–blue); (B) Hierarchical tree dendrogram; (C) Barplot with posterior probabilities of classification of each individual into each group from the results of the classification discriminant analysis. Acronyms of populations: P01—Psunj; P02—Tounj; P03—Istria; P04—Novi Vinodolski; P05—Split; P06—Brač; P07—Konavle.
Plants 10 01691 g001
Figure 2. Results of the multivariate statistical methods for fruit chemical analysis and locations of the seven sampled Sorbus domestica populations. (A) Geographical distribution of three groups of populations detected from K-means clustering (the proportions of the ancestry of each population in each of the defined clusters are color-coded: cluster A– green, cluster B–yellow, cluster C–blue); (B) Hierarchical tree dendrogram; (C) Barplot with posterior probabilities of classification of each individual into each group from the results of the classification discriminant analysis. Acronyms of populations: P01—Psunj; P02—Tounj; P03—Istria; P04—Novi Vinodolski; P05—Split; P06—Brač; P07—Konavle.
Figure 2. Results of the multivariate statistical methods for fruit chemical analysis and locations of the seven sampled Sorbus domestica populations. (A) Geographical distribution of three groups of populations detected from K-means clustering (the proportions of the ancestry of each population in each of the defined clusters are color-coded: cluster A– green, cluster B–yellow, cluster C–blue); (B) Hierarchical tree dendrogram; (C) Barplot with posterior probabilities of classification of each individual into each group from the results of the classification discriminant analysis. Acronyms of populations: P01—Psunj; P02—Tounj; P03—Istria; P04—Novi Vinodolski; P05—Split; P06—Brač; P07—Konavle.
Plants 10 01691 g002
Figure 3. The boundaries obtained with Monmoniers’s maximum difference algorithm of Euclidean distances among populations revealed by program BARRIER. (A) The first three boundaries (a–c) for morphometric variation; and (B) the first two boundaries (a and b) for chemical variation. The thickness of barriers is related to their importance. Acronyms of populations: P01—Psunj; P02—Tounj; P03—Istria; P04—Novi Vinodolski; P05—Split; P06—Brač; P07—Konavle.
Figure 3. The boundaries obtained with Monmoniers’s maximum difference algorithm of Euclidean distances among populations revealed by program BARRIER. (A) The first three boundaries (a–c) for morphometric variation; and (B) the first two boundaries (a and b) for chemical variation. The thickness of barriers is related to their importance. Acronyms of populations: P01—Psunj; P02—Tounj; P03—Istria; P04—Novi Vinodolski; P05—Split; P06—Brač; P07—Konavle.
Plants 10 01691 g003
Figure 4. Biplot of the principal component analysis based on ten fruit and seed morphometric characteristics in studied Sorbus domestica populations. Acronyms of populations: P01 (Psunj), P02 (Tounj), P03 (Istria), P04 (Novi Vinodolski), P05 (Split), P06 (Brač), and P07 (Konavle). Each point represents a single tree, with a color corresponding to cluster from K-means clustering method and origin of population: continental populations—cluster A (green); island population Brač—cluster B (yellow); and coastal Mediterranean populations—cluster C (blue).
Figure 4. Biplot of the principal component analysis based on ten fruit and seed morphometric characteristics in studied Sorbus domestica populations. Acronyms of populations: P01 (Psunj), P02 (Tounj), P03 (Istria), P04 (Novi Vinodolski), P05 (Split), P06 (Brač), and P07 (Konavle). Each point represents a single tree, with a color corresponding to cluster from K-means clustering method and origin of population: continental populations—cluster A (green); island population Brač—cluster B (yellow); and coastal Mediterranean populations—cluster C (blue).
Plants 10 01691 g004
Figure 5. Biplot of the principal component analysis based on ten chemical traits in studied Sorbus domestica populations. Acronyms of populations: P01 (Psunj), P02 (Tounj), P03 (Istria), P04 (Novi Vinodolski), P05 (Split), P06 (Brač), and P07 (Konavle). Each point represents a single tree, with a color corresponding to cluster from K-means clustering method and origin of population: continental populations—cluster A (green); island population Brač—cluster B (yellow); and coastal Mediterranean populations—cluster C (blue).
Figure 5. Biplot of the principal component analysis based on ten chemical traits in studied Sorbus domestica populations. Acronyms of populations: P01 (Psunj), P02 (Tounj), P03 (Istria), P04 (Novi Vinodolski), P05 (Split), P06 (Brač), and P07 (Konavle). Each point represents a single tree, with a color corresponding to cluster from K-means clustering method and origin of population: continental populations—cluster A (green); island population Brač—cluster B (yellow); and coastal Mediterranean populations—cluster C (blue).
Plants 10 01691 g005
Figure 6. Isolation-by-distance and Isolation-by-environmental distance. Scatter plots of simple Mantel tests showing the relationships between (A) geographic and morphological distances (r = –0.0729, p = 0.58), (B) environmental and morphological distances (r = 0.3814, p = 0.05), (C) geographic and chemical distances (r = –0.0848, p = 0.47), (D) environmental and chemical distances (r = 0.0706, p = 0.37), and (E) morphological and chemical distances (r = 0.6129, p = 0.02).
Figure 6. Isolation-by-distance and Isolation-by-environmental distance. Scatter plots of simple Mantel tests showing the relationships between (A) geographic and morphological distances (r = –0.0729, p = 0.58), (B) environmental and morphological distances (r = 0.3814, p = 0.05), (C) geographic and chemical distances (r = –0.0848, p = 0.47), (D) environmental and chemical distances (r = 0.0706, p = 0.37), and (E) morphological and chemical distances (r = 0.6129, p = 0.02).
Plants 10 01691 g006
Table 1. Arithmetic means (M) and coefficient of variations (CV/%) for fruit and seed morphometric characteristics. Different letters within one column denote statistically significant differences (p < 0.05) by Fisher’s LSD test. Acronyms of populations: P01—Psunj; P02—Tounj; P03—Istria; P04—Novi Vinodolski; P05—Split; P06—Brač; P07—Konavle. Fruit and seed morphometric characteristics: m (g)—fruit mass; FL (cm)—fruit length; MFW (cm)—maximum fruit width; PMFW (cm)—fruit length, measured from the fruit base to the point of maximum fruit width; FW1 (cm)—fruit width at 15% of fruit length; FW2 (cm)—fruit width at 85% of fruit length; FW/FL—maximum fruit width/fruit length; SL (mm)—seed length; SW (mm)—maximum seed width; NS—number of seeds per fruit.
Table 1. Arithmetic means (M) and coefficient of variations (CV/%) for fruit and seed morphometric characteristics. Different letters within one column denote statistically significant differences (p < 0.05) by Fisher’s LSD test. Acronyms of populations: P01—Psunj; P02—Tounj; P03—Istria; P04—Novi Vinodolski; P05—Split; P06—Brač; P07—Konavle. Fruit and seed morphometric characteristics: m (g)—fruit mass; FL (cm)—fruit length; MFW (cm)—maximum fruit width; PMFW (cm)—fruit length, measured from the fruit base to the point of maximum fruit width; FW1 (cm)—fruit width at 15% of fruit length; FW2 (cm)—fruit width at 85% of fruit length; FW/FL—maximum fruit width/fruit length; SL (mm)—seed length; SW (mm)—maximum seed width; NS—number of seeds per fruit.
PopulationmFLMFWPMFWFW1FW2FW/FLSLSWNS
MCVMCVMCVMCVMCVMCVMCVMCVMCVMCV
P019.53f27.342.37f10.242.40f11.150.99f17.721.91c12.671.56c16.121.02a9.847.26e9.065.51e10.831.35ab51.85
P028,91e32.222.26e16.032.31a12.890.92e19.301.87d15.671.53c14.221.03a10.756.68c13.665.23d14.331.48c54.40
P037.01a26.032.09a7.362.04c8.920.85a15.321.60a13.151.28a14.710.98e10.416.72cd11.054.85a13.721.46bc37.14
P047.49d31.452.15d11.562.31a12.020.85a21.121.92c13.811.46b17.331.08d6.976.39a9.425.07b10.651.25a47.59
P056.86a30.952.06a11.442.21e12.560.88d18.421.76b15.551.46b15.201.07cd8.536.80d9.545.12b9.591.97d46.18
P065.09b36.121.88b14.961.98b12.430.75b20.901.61a13.461.25a16.141.06bc9.186.48ab10.194.42c14.563.66e25.49
P076.34c28.541.99c11.842.09d10.920.78c19.471.75b12.671.32d18.411.06b11.836.53b13.554.87a14.092.00d56.25
Total7.3236.042.3710.242.1913.480.8620.781.7715.621.4118.071.0410.176.7011.715.0114.081.8860.06
Table 2. Results of the hierarchical analysis of variance for fruit and seed morphometric characteristics. Acronyms for fruit and seed morphometric characteristics as in Table 1.
Table 2. Results of the hierarchical analysis of variance for fruit and seed morphometric characteristics. Acronyms for fruit and seed morphometric characteristics as in Table 1.
VariableComponents of the VariancedfFPercent of Variabilityp-Value
mPopulation65.2725.440.000401
Tree (Population)4262.2241.040.000000
Error 33.52
FLPopulation64.7522.680.000898
Tree (Population)4259.4741.680.000000
Error 35.64
MFWPopulation65.2422.490.000422
Tree (Population)4245.1336.340.000000
Error 41.17
PMFWPopulation66.3416.570.000084
Tree (Population)4217.2720.490.000000
Error 62.94
FW1Population64.4817.660.001365
Tree (Population)4237.1934.570.000000
Error 47.77
FW2Population64.5018.360.001325
Tree (Population)4240.0835.820.000000
Error 45.82
FW/FLPopulation61.322.480.270911
Tree (Population)4262.5153.800.000000
Error 43.73
SLPopulation61.524.480.196602
Tree (Population)4285.9360.120.000000
Error 35.40
SWPopulation65.2018.650.000445
Tree (Population)4230.3530.090.000000
Error 51.26
NSPopulation613.2747.020.000000
Tree (Population)4250.2426.290.000000
Error 26.69
Table 3. Arithmetic means (M) and coefficient of variations (CV/%) for water (g/100 g), crude protein (g/100 g), sugar (g/100 g), ash (g/100 g), crude fat (g/100 g), cellulose (g/100 g), acidity (%), total phenols (mg GAE/g), and antioxidant activity (DPPH %; FRAP mmol Fe2+). All mass fractions were determined on a dry mass basis. Different letters within one column denote statistically significant differences (p < 0.05) by ANOVA and Fisher’s LSD test. Acronyms of populations: P01—Psunj; P02—Tounj; P03—Istria; P04—Novi Vinodolski; P05—Split; P06—Brač; P07—Konavle.
Table 3. Arithmetic means (M) and coefficient of variations (CV/%) for water (g/100 g), crude protein (g/100 g), sugar (g/100 g), ash (g/100 g), crude fat (g/100 g), cellulose (g/100 g), acidity (%), total phenols (mg GAE/g), and antioxidant activity (DPPH %; FRAP mmol Fe2+). All mass fractions were determined on a dry mass basis. Different letters within one column denote statistically significant differences (p < 0.05) by ANOVA and Fisher’s LSD test. Acronyms of populations: P01—Psunj; P02—Tounj; P03—Istria; P04—Novi Vinodolski; P05—Split; P06—Brač; P07—Konavle.
PopulationWaterCrude ProteinSugarAshCrude FatCelluloseAcidityTotal PhenolsDPPHFRAP
MCVMCVMCVMCVMCVMCVMCVMCVMCVMCV
P0169.50d2.7815.32ac14.4749.94a4.811.90ab16.390.59d16.015.93ab10.680.74a13.1411.11c63.1911.99a10.502.41ac12.58
P0267.81ad2.1812.49b20.2949.44a8.492.01ab44.110.60bd9.226.01ab20.120.64a19.7912.10c16.804.40b54.742,91bc1.71
P0366.52ac1.5916.85a18.6150.70a6.952.19a14.230.80ac12.985.54ab6.420.70a8.259.44bc38.2317.76c19.422.49ab4.20
P0467.05a2.0912.50b8.4155.60c4.991.44b20.050.85c19.535.33ab8.580.68a10.046.57ab50.3313.61a47.362.69abc5.94
P0564.24bc3.0116.67a14.7349.47a3.902.29ac8.090.78ac7.475.13a11.200.67a5.245.62a42.1713.76a25.882.76bc12.32
P0656.80e6.2815.63ac19.5642.88b7.662.91c33.900.72ab15.629.06c32.340.64a28.203.50a57.4411.70a12.132.50ab12.77
P0764.05b3.9913.41bc21.9647.50a6.832.13ab41.780.71ab22.046.65b10.830.67a11.266.82ab53.6512.57a39.972.50ab23.37
Total65.146.7014.6920.2749.369.332.1234.230.7219.526.2427.950.6814.907.8857.4811.8144.002.6113.34
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

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. https://doi.org/10.3390/plants10081691

AMA Style

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(8):1691. https://doi.org/10.3390/plants10081691

Chicago/Turabian Style

Poljak, Igor, Nada Vahčić, Zlatko Liber, Katarina Tumpa, Valentino Pintar, Ivana Zegnal, Antonio Vidaković, Bernarda Valković, Davorin Kajba, and Marilena Idžojtić. 2021. "Morphological and Chemical Diversity and Antioxidant Capacity of the Service Tree (Sorbus domestica L.) Fruits from Two Eco-Geographical Regions" Plants 10, no. 8: 1691. https://doi.org/10.3390/plants10081691

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