Some 49 species of woody plants were identified among the tree plantations in sanitary protection zones of industrial enterprises in Zaporizhzhya (southeastern Ukraine); 10 species of tree plants are present in all the sanitary protection zones. The number of plant species in protective plantings varies considerably from 11 ("Zaporizhskloflus") to 30 species ("Zaporizhstal"). The similarity of the species composition in plantations was assessed using the Jaccard and Sørensen indices. We calculated the species richness by Margaleff and Shannon indices, the dominance indices were calculated according to Simpson and Berger-Parker. The vast majority of forest belts in the sanitary protection zones of Zaporizhzhya industrial enterprises require the species enrichment, taking into account the tolerance of species to various pollutants.
In nature, the creation of infra‒specific difference can occur through a wide range of mechanisms such as local adaptation, phenotypic plasticity, parental conditions and artificial selection . West-Eberhard  has believed that when it is resulted from plasticity, features can differ rapidly within generations and vary drastically across populations in dissimilar habitats. Recently, studies have broadened in the identification of variation to incorporate the considerable phenotypic difference within and among populations of the same species .
Investigations have showed that difference of phenotypic feature within species can be as extreme as the characteristic difference across species . Moreover, infra‒specific difference may influence structure of community and function of ecosystem as much as difference among species [6,14].
Salvia limbata is an aromatic herb of Lamiaceae that grows naturally in various regions of Iran and some neighboring countries [7,16]. Different flavone compounds such as ladanein, salvigenin, luteolin 7‒methyl ether, cirsiliol, eupatorin, luteolin 7‒O‒glucoside and rosmarinic acid were isolated from extracts of this plant . Moreover, Saeidnia et al.  have suggested that S. limbata can accumulate tryptophan, sterols and their glucosides, therefore, consumption of this plant as an herbal tea or other preparations might be useful for dietary deficiency of this amino acid. In mice, S. limbata has beneficial effects to decrease dependence sign produced by morphine and increase in pain threshold compared with the control .
Although, comparative infraspecific anatomical  and phytochemical studies [10,13] have been performed in this species, as far we could search, no infra‒specific morphological investigation is available for its populations.
In current study, we studied morphological characteristics of inflorescence in eight populations of S. limbata, in order to evaluate pattern of morphological variations among the populations and also find possible phenotype(s) among the studied populations.
Material and methods
In current study, eight natural populations of S. limbata were harvested from different regions of the country (Table 1), and were identified according to descriptions provided in Flora of Iran . Morphological differences were carefully evaluated for reproductive organs from collected samples of the studied populations. The morphological characteristics of flowers were studied under dissecting stereo microscopes. For morphometric analysis, 80 plant specimens from 8 populations were used. Voucher specimens were deposited at Herbarium of Arak University.
|1||Markazi province, arak, Sardasht, 1850 m.|
|2||Markazi province, Saveh, Sangak, 1940 m.|
|3||Khorasan Razavi, Semnan to Ghochan, 1730 m.|
|4||Markazi province, Zarandiyeh, Vidar, 1750 m.|
|5||Semnan province, Ahovan, 1720 m.|
|6||Markazi province, Zarandiyeh, 1430 m.|
|7||Semnan province, Ghochan to Semnan|
|8||Tehran province, Gajroud, 1600 m.|
Measurements of morphological traits were performed on each plant, including its flowering stem with flowers in bloom. In total, 21 qualitative features were measured on each specimen. Each character was measured four times per each plant sample, and their average determined. The studied characteristics were: calyx width, length and calyx length / width ratio, petal length, width and petal length / wide ratio, calyx short teeth length, calyx short teeth width, calyx short teeth length / width ratio, style length, stigma length, calyx long teeth length, calyx long teeth width, calyx long teeth length / width ratio, anther length, anther width, anther length /width ratio, pedicle length, long filament length, short filament length, long /short filament length ratio.
Morphometric data was subjected to one‒way analysis of variance (ANOVA) to determine if significant difference existed among populations for each feature measured. Mean and standard deviations of characteristics were tacalculated. These analyses were performed using SPSS ver. 17. Cluster analysis was carried out based on quantitative features using UPGMA, PCO, PCA and C.A‒Joined plots clustering in Multivariate Statistical Package (MVSP) program .
The mean and standard deviation of the studied morphological variables were presented in Table 2.
|Population||Anther length||Antler width||Stigma length||Style length||Short Filament length||Long filament length||Petal length||Petal width||calyx Short teeth length||calyx Short teeth width||calyx long teeth length||calyx long teeth width||Calyx length||Calyx width|
Morphological traits differed among the studied populations. Biggest and smallest anther lengths were recorded in population no. 5 and 1, respectively. While, reverse pattern were recorded for long filament length. Therefore, largest and shortest long filament lengths were reported for populations no. 1 and 5, respectively. Maximum and minimum lengths of pedicle were belonged to population’s no. 3 and 2, respectively. Largest stigma was observed in population no. 1, but population no. 4 had smallest stigma. Longest style was found in population no. 4, and shortest in population no. 5. Maximum and minimum petal lengths were registered in populations no. 5 and 2, respectively. Population no. 1 had longest calyx and widest calyx short teeth; however, shortest calyx and narrowest calyx short teeth were belonged to population no. 7. Largest and smallest short filament length was found in populations no. 6 and 5, respectively (Fig. 1).
Furthermore, the ANOVA test confirmed significant variations (p ≤ 0.01) for all studied variables, except for anther width (Table 3).
|Characteristics||Sum of Squares||df||Mean Square||F||Sig.|
|Anther length||Between Groups||9.573||9||1.064||5.741||0.000|
|Anther width||Between Groups||1.010||9||0.112||1.836||0.072|
|Anther length/width ratio||Between Groups||2.387||9||0.265||3.696||0.001|
|Pedicle length||Between Groups||171.202||9||19.022||22.698||0.000|
|Stigma length||Between Groups||10.962||9||1.218||5.420||0.000|
|Style length||Between Groups||277.690||9||30.854||2.654||0.009|
|Short Filament length||Between Groups||6.772||9||0.752||4.463||0.000|
|Long filament length||Between Groups||11.089||9||1.232||6.199||0.000|
|Petal length||Between Groups||171.200||9||19.022||3.524||0.001|
|Petal width||Between Groups||18.840||9||2.093||6.587||0.000|
|Petal length/width ratio||Between Groups||4.134||9||0.459||2.408||0.017|
|Calyx Short teeth length||Between Groups||22.404||9||2.489||13.065||0.000|
|Calyx Short teeth width||Between Groups||4.532||9||0.504||6.069||0.000|
|Calyx Short teeth length/width ratio||Between Groups||2290.906||9||254.545||11.111||0.000|
|Calyx long teeth length||Between Groups||41.789||9||4.643||8.759||0.000|
|Calyx long teeth width||Between Groups||11.206||9||1.245||9.699||0.000|
|Calyx long teeth length/width ratio||Between Groups||31225.271||9||3469.475||11.935||0.000|
|Calyx length||Between Groups||84.810||9||9.423||6.812||0.000|
|Calyx width||Between Groups||106.462||9||11.829||10.622||0.000|
|calyx length/width ratio||Between Groups||7.300||9||0.811||10.565||0.000|
Fig. 4. PCA plot of the studied populations based on morphological characteristics.
CA. joined plot (Fig. 5) revealed each of these groups were characterized by special features that were useful in identification of them. For example, calyx length, calyx short teeth length and calyx long teeth length were prominent variables for populations no.1. Moreover, populations no. 2 and 4 were characterized by style length, calyx short teeth length/width ratio and calyx long teeth length/width ratio. Populations no. 5 was identified according to calyx length/width ratio.
In this research, we evaluated morphological variables of inflorescence in the selected populations of S. limbata. Because, these characteristics are more stable than vegetative characteristics , in addition, Jamzad  used inflorescence characteristics along with palynological as well as ITS data in her valuable work on Iranian Lamiaceae species, Flora of Iran.
Morphological features highly varied among the studied populations. Furthermore, previous studies have showed phytochemical and anatomical variations among various populations of S. limbata. For example, Mirza et al.  identified germacrene D (25.7%), linalool (17.5%) and linalyl acetate (16.1%) as the major components of essential oil in Iranian population of S. limbata, while in Turkey samples ᾳ‒pinene (11.2–24.3%), β‒pinene (10.0–20.9%), and sabinene (14.6–17.4%) were the major compounds . As observed, the type and amount of major essential oil compounds were completely different between these populations. Moreover, Nejadhabibvash et al.  studied effect of ecological gradients on S. limbata anatomical traits. Several quantitative anatomical characteristics related to cuticle, collenchymas, pith, phloem, xylem, stomata frequency, epidermal number and stomata length/ width ratio differed among populations.
Fig. 5. CA‒Joined plot of the studied populations with morphological characteristics. Abbreviations: blue symbols are population’s number according to tables 1, calwi: calyx width, callen: calyx length, stle: style length, calshtele: calyx short teeth length, callontele: calyx long teeth length, callenwidrat: calyx length/width ratio, callotelewira: calyx long teeth length /width ratio, calshtelewira: calyx short teeth length /width ratio, petlenwira: petal length /width ratio.
According to UPGMA tree, PCA and PCO plots of morphological characteristics, we found four primitive phenotypes among these populations: phenotype no. 1 (population no. 1), phenotype no. 2 (population no. 5), phenotype no. 3 (populations no. 2 and 4), and phenotype no. 4 (including the rest populations). Each phenotype was characterized by special morphological characteristics.
S. limbata has large distribution range in Iran and grows in environmental heterogeneity habitats. Its populations achieve adaptations for living under different ecological conditions, therefore phenotype plasticity occurs in order to adaptation with ecological conditions.
Wherever the population’s habitats are the same, the morphological characteristics of populations will be similar, regardless of geographical distance of populations. Studies of population on the pattern of variation in several plant species have revealed the existence of localized populations each adapted to the particular environmental conditions of their habitat .
Phenotype no.3 had two populations. The habitats of these populations were close together with similar ecological conditions; therefore, morphological similarity of them was very probable. Moreover, low geographical distance facilitates gene flow between populations and leads to lower differentiation. Hamrick and Godt  have suggested that higher inter‒populational gene flow leads to lower differentiation among populations.
Phenotype no. 4 consisted of four populations. The distance between populations no. 3 and 6 was more than 700 km. In addition, it was ca 350 km between populations no. 7 and 8. Contrary with very large geographical distance these populations were clustered as a group. Moreover, in comparison to phenotype no. 3, phenotype no. 2 was less distant with phenotype no. 4, while according to UPGMA tree, PCO and PCA plots phenotype 3 was closer to phenotype no. 4 rather than phenotype no. 2.
These authors propose two main possible reasons for this similarity, despite the existence of a long geographical distance. The genetic similarity is the first possible case, Sheidai et al.  have suggested that the morphological similarity among populations may be due to genetic likeness, although, it is more important to know that phenotype and genotype are not necessarily tightly linked. For example, Liston  and Whitkus  have revealed high levels of genetic variations among morphologically similar populations of Astragalus sect. Leptocarpi and Carex pachystach, respectively.
What, unlike the long distance, has led to the similarity of populations is similarity of ecological conditions of habitats. This is well suited to the idea of ecological mosaic. Jain and Bradshaw  have suggested that if a plant species encounters a mosaic of habitats while expanding its range, the selective influences of each habitat act in turn upon the whole available morphologic difference.
Talebi et al.  studied morphological and ecological features of several populations of Linum album in Iran. They found that there is a strong relationship between morphology of populations and their ecological conditions. Wherever the ecological conditions were similar, the morphological characteristics of the plants were alike.
So, it seems that in the case of our studied populations, ecological factors have a very strong effect on population’s similarity or difference. This condition creates morphological polymorphism among populations. Of course, this situation gives the species that has the genetic background to occupy a variety of environments and emergence of morphological polymorphism.
We evaluated inflorescence morphological characteristics among eight populations of S. limbata. Morphological traits highly differed among the populations and ANOVA test revealed significant difference among the populations. According to UPGMA tree, PCA and PCO plots, we found four primitive phenotypes, which were characterized by special morphological features. Two phenotypes were monotypic, while one phenotype consisted of two, and another one has four populations. The distribution domain of a species is not composed exactly identical ecological environments, but also a set of different environment with different ecological conditions. So species must adapt its phenotype structure in accordance with those conditions in order to be able to live under those conditions.
- Albert CH, Thuiller W, Yoccoz NG, Douzet R, Aubert S, Lavorel S. A multi‒trait approach reveals the structure and the relative importance of intra‒ vs. interspecific variability in plant traits Functional Ecology. 2010; 24:1192-1201. DOI
- Alemy SH, Karami M, Hossini E, Ebrahimzadeh MA, Shahbi-Majd N. Antinociceptive activity and effect of methanol extract of Salvia limbata on withdrawal syndrome in mice. European Review for Medical and Pharmacological Sciences. 2012; 16:38-42.
- Bradshaw AD. Genostasis and the limits to evolution. Philosophical Transactions of the Royal Society of London Series Biological Sciences. 1991; 333:289-305.
- Gohari AR, Saeidnia S, Malmir M, Hadjiakhoondi A, Ajani Y. Flavones and rosmarinic acid from Salvia limbata. Natural Product Research. 2010; 24(20):1902-1906. DOI
- Govaert L, Pantel JH, De Meester L. Eco‒evolutionary partitioning metrics: assessing the importance of ecological and evolutionary contributions to population and community change. Ecology Letters. 2016; 19:839-853. DOI
- Guner A, Ozhatay N, Ekim T, Baser KHC. Flora of Turkey and the East Aegean Islands (11). Edinburgh: Edinburgh University Press; 2000.
- Jain SK, Bradshaw AD. Evolutionary divergence among adjacent plant populations. I The evidence and its theoretical analysis Heredity. 1996; 21:407-441. DOI
- Jamzad Z, Lamiaceae Flora of Iran. Research Institute of Forest and Rangelands. Iran: Iran. 2012.
- Kurkcuoglu M, Demirci B, Baser KHC, Dirmenci T, Tumen G, Ozgen U. The essential oil of Salvia limbata C. A Meyer Growing in Turkey Journal of Essential Oil Research. 2005; 17:192-195. DOI
- Liston A. Isozyme systematics of Astragalus sect. Leptocarpi subsect Californici (Fabaceae) Systematic Botany. 1992; 17:367-379. DOI
- Nejadhabibvash F, Rezaei-Chiyaneh E, Pirzad A.
- Mirza M, Mozaffarian V, Nik ZB. Composition of the Essential Oil of Salvia limbata C. A Mey Journal of Essential Oil Research. 2005; 17(1):10-11. DOI
- Palkovacs EP, Post DM. Experimental evidence that phenotypic divergence in predators drives community divergence in prey. Ecology. 2009; 90:300-305. DOI
- Podani J.
- Rechinger KH. Labiatae in Flora Iranica. Astria: Akademische Druck-u Verlagsanstalt; 1982.
- Saeidnia S, Gohari AR, Malmir M, Moradi-Afrapoli F, Ajani Y. Tryptophan and Sterols from Salvia limbata. Journal of Medicinal Plants. 2011; 10(37):41-47.
- Sheidai M, Afshar F, Keshavarzi M, Talebi SM, Noormohammadi Z, Shafaf T. Genetic diversity and genome size variability in Linum austriacum (Linaceae) populations. Biochemical Systematics and Ecology. 2014; 57:20-26. DOI
- Stace CA. Plant Taxonomy and Biosystematics. Edward Arnold: London; 1989.
- Talebi SM, Atri M, Sheidai M, Sharifnia F, Noormohammadi Z. Infraspecific variations in Linum album based on the determination of special stations approach in Iran. Phytologia Balcanica. 2014; 20(1):9-22.
- Talebi SM, Rezakhanlou A, Matsyura A. Do we have infraspecific taxa of Salvia multicaulis Vahl (Lamiaceae) in Iran? Ukrainian Journal of Ecology. 2017; 7(4):432-439. DOI
- Violle C, Enquist BJ, McGill BJ, Jiang L, Albert CH, Hulshof C, Jung V, Messier J. The return of the variance: intraspecific variability in community ecology. Trends in Ecology & Evolution. 2012; 27(4):244-252. DOI
- West-Eberhard MJ. Phenotypic plasticity and the origins of diversity. Annual Review of Ecology and Systematics. 1989; 20:249-278. DOI
- Whitkus R. Allozyme variation within the Carex pachystachya complex (Cyperaceae). Systematic Botany. 1992; 17:16-24. DOI