Baranov S., Zykov I., Fyodorova L. INFLUENCE OF CLIMATIC FACTOR ON PLASTICITY AND ASYMMETRY OF LAMINAS OF GENUS TILIA L. SPECIES // Principy èkologii. 2019. № 1. P. 4‒18. DOI: 10.15393/

Issue № 1

Original research



   Sergey Gennadevich
Ph.D., A.G. and N.G.Stoletov Vladimir State University, Russia 600000, Vladimir Gorky st., 87,
   Igor Evgenyevich
State humanitarian technological University, 142611 Orekhovo-Zuyevo Zelenaya st., 22,
   Lyubov Valeryevna
Federal State Autonomous Educational Institution of Higher Education I.M. Sechenov First Moscow State Medical University of the Ministry of Health of the Russian Federation (Sechenov University), 119991, Moscow, Trubetskay st., 8-2,
Tilia cordata
Tilia platyphyllos
fluctuating asymmetry
phenotypic plasticity
stability of development
fluctuation variability
Summary: The properties of bilaterally symmetric laminas of two types: small-leaved limes (Tiliacordata Mill.) and linden broadleaf ones (Tilia platyphyllos Scop.) were compared in the contaminated zone and in the control one (Moscow, 2016–2018). The fluctuating asymmetry (FA) depended on the year and gathering place of laminas, while the plastic variability depended on the year of gathering. The value of FA was higher in T. platyphyllos than in T. cordata. In 2017, the FA of Tilia platyphyllos in the park area (control) was lower than on a traffic highway. The reasons for the increase of the FA were: a) the cold summer of 2017, which inhibited the laminas growth and increased the plastic variability of metric traits; b) motor vehicle pollution. Thus, the climatic conditions of the year contributed to both fluctuation and plastic variability that did not mean a change in the stability of development.

© Petrozavodsk State University

Reviewer: I. Korotchenko
Received on: 01 September 2018
Published on: 28 March 2019


The index of fluctuating asymmetry (FA) is a popular indicator, which is characterized by a slight deviation from zero of the difference between the values of the right and left parts of the homologous bilaterally symmetric trait. The dimensional or countable bilaterally symmetric traits with a wide range of response to stress factors are used to determine the value of FA and developmental stability / instability (Palmer, Strobeck, 2003; Baranov et al., 2015). They include the most genotypically or phenotypically variable traits. In representatives of the genus Tília the dimensional, bilaterally symmetrical traits of leaf plates have high variability, which is associated with species-specific features. It is known that developmental stability (DS) is the ability of an organism to compensate and reduce deviation from normal ontogenetic development along a certain channeling pathway (Debat, David, 2001; Lens et al., 2002; Klingenberg, 2016).

Fluctuating asymmetry is of interest both in ecological terms and in the aspect of evolution and taxonomy of plant organisms. Representatives of the Birch and Linden families are widely used to determine the level of stability of development (Zorina, Korosov, 2015; Zykov et al., 2015), however, the causes of fluctuating changes in the morphological structures of bilateral organs remain poorly understood. The use of tree species as an indicator of development stability in points faces difficulties related to the genotypic purity of the used species, population heterogeneity, the influence of factors, including lighting, water-salt composition and soil regime, terrain, other physical-chemical and other features of the environment.

Methodologically, the determining the level of developmental stability raises some questions. Methods of measurement and preparation of herbarium material occupy an important part and imply a "blind" unbiased approach (Kozlov, Niemelä, 2003; Kozlov, Zvereva, 2015). Traditional methods, such as normalizing difference, are supposed to work with the use of individual traits with a normal distribution of the difference between the right and left values. This approach still finds application, although the evaluation using one or two or many traits does not often provide a representative recording of the signal in the form of an increased value of FA. The rapidly developing method of geometric morphometry allows giving a more harmonious integrative assessment, bringing researchers closer to the true causes for the increase or decrease in FA. The development of epigenetics has led to a number of new studies. For example, works carried out on several generations indicate genotypic conditionality of the processes of decrease-increase in FA (Clarke, 1993; Freeman et al., 1993).

It remains to be recognized that "the health of the environment» is a term which hardly pave the way in bioindication. This opinion is derived from the results of numerous studies in the field of dynamics of the value of FA, modeling experiments in vitro, comparative studies on FA in different species, populations, works with negative results (Zverev et al., 2018), finally, works with justification of paradoxical effects of a dose of stress on FA and developmental stability (Erofeeva, 2015). Thus, studies of FA and DS depending on the gradient of pollution by technogenic toxins show contradictory results: both reverse (Kaligarič et al., 2008, Gostin, 2015) and direct relationship (Belyaeva, 2013; Sobchak et al., 2013; Koroteeva et al., 2015; Ivanov et al., 2015), as well as uncertain or conditionally positive results (Leung, 1996; Lung, 2000). At the same time, the study of FA of leaf plates in intact ecosystems showed a serious impact of terrain height and other abiotic (Zvereva et al., 1997; Wuytack et al., 2011; Baranov, 2014a; Baranov, 2014b) and biotic factors (Zvereva, 1997; Kozlov et al., 2017).

The influence of climatic factors was studied in detail on the example of oak (Hódar, 2002), pine (Kozlov, Niemelä, 2003), white birch (Valkama, Kozlov, 2001) and some herbaceous plants (Alados et al., 2001). The combined effect of many factors, such as humidity and low temperature remains unclear, so these factors against the background of anthropogenic stress are of interest. The logical continuation of the research is background monitoring, studying the stages of morphogenesis of leaf plates (Anastasia, 2014) and work in the field of evolutionary development along a certain, channeling path. The latter term is associated with the property of plasticity of development (plasticity of development; Graham, 2010). In Russian literature, the familiar terms continue to be used – phenotypic or ecological plasticity. It is known that the ecological plasticity of plants is caused by the buffer capacity of morphological structures, which allows them to actively adapt to environmental conditions. The relationship between plasticity and developmental stability continues to be debated. At the same time, the question of fluctuating asymmetry, as an adaptive signal manifested in the trajectory of population development, is debated (Clarke, 1993; Debat, David, 2001; Graham et al., 2010). The modular concept of ecological plasticity is detailed in many sources, for example, De Kroon (De Kroon et al., 2005). According to this concept, the norm of reaction, as the main ecological characteristic of the organism, covers the sum of the plasticity modules of the root, leaf and stem. Evolution and adaptive properties of plastic variability (PV) under the conditions of environmental perturbation are also described in many review papers (Lande, 2009; Gilbert et al., 2015; Abley et al., 2016), including paleontological studies of the size and shape of leaf plates on the example of maple subspecies in the temperature gradient (Royer et al., 2009). A well-established opinion is the correlation between developmental stability, phenotypic plasticity and adaptability at the genetic level (Tonsor et al., 2013). Therefore, plastic variability is of interest as a companion of fluctuation variability both under anthropogenic stress and under severe climatic conditions. Despite the high variation in the characteristics of linden leaves and the frequent presence of directional asymmetry that prevents the determination of FA compared to species of genus Betula, representatives of the genus Linden are more convenient due to a more stable chromosomal karyotype. It is known that close species differ in FA value in response to carbon dioxide content in the environment and emissions of industrial enterprises that confirms the genotypic basis of development stability (Hochwender, Fritz, 1999; Zorina, Korosov, 2015; Koroteeva et al., 2015).

The aim of this work was to determine the stability of development and phenotypic plastic variability in two close species in two urbanized stations with different degrees of anthropogenic load. The task was to study the influence of climate-temperature factors on the stability of development and plastic variability of leaf plates of close species of the genus Linden. The hypothesis being tested was as follows: fluctuation variability is associated with plastic variability, but in different species this relationship is manifested differently, which is explained by the genotypic features of the species.


Collecting places and preparation of herbarium material

It is known that the big-leaf linden (Tilia platyphyllos Scop.) is a more eurybiontic species in relation to soil and climatic conditions. Both species are equally demanding to light, which is convenient for studying genotypic and phenotypic properties in the context of asymmetry and plastic variability of leaf plates. Populations of both species grow in the Eastern part of Moscow, in Izmailovsky Park and on Izmailovsky Boulevard. Izmailovo natural and historical Park (55°47'N; 37°47'E) is a biocenosis with more than 500 species of plants. Based on this, individuals of trees from different parts of the park, i.e. from different phytocenoses, were selected.

Izmailovsky Boulevard (55°48 '13.8"N; 38°58'23.8" E) is the center of a busy highway, the area of the station on the Boulevard was about 5 km2. The characteristic features of this area were the uniformity of lighting, aeration, lack of competition and a high level of atmosphere pollution by carbon oxides, nitrogen, sulfur and hydrocarbons.

The leaf plates were collected at the generative stage of trees development in 2016-2018. In each station, leaf plates 3-4 cm wide were uniformly collected from the lower parts of the crowns of ten trees of the same age, 10 samples from each of the 10 trees. For measurement we used the traits originally developed by the laboratory of V.M.Zakharov (Zakharov, Chubinishvili, 2001). As previous studies have shown (Baranov et al., 2015; Zykov et al., 2015), the angular trait was clearly inconvenient to measure due to the high degree of curvature of the first lateral vein. It was replaced by an alternative to it – linear: the distance between the base of the second vein of the 1st order and the base of the first vein of the 2nd order on the first vein of the 1st order (Fig. 1).


Fig. 1.  The traits used to determine the fluctuating asymmetry of species of the genera Tilia L.: 1) width in the middle of a leaf; 2) distance between the bases of the first vein of the 1st order and the second vein of the 2nd order; 3) distance between the bases of the second and third veins of the 2nd order; 4) distance between the bases of the first and second veins of the 1st order; 5) distance between the base of the second vein of the 1st order and the base of the first vein of the 2nd order on the first vein of the 1st order.


The leaves were dried under a paper press, and after 2-3 days were measured with a measuring ruler with an accuracy of 0.5 mm (GOST 427-75), the data were stored in Excel tables. The triple measurement showed an extremely low error, less than 1 % of the FA value, which is considered acceptable when testing the asymmetry (Palmer, Strobeck, 2003). Taking into account the data of Velicković (Velicković, 2010) on the greatest sensitivity of the trait of leaf length and previous results, in which another trait was the most suitable (No. 4), we decided to use the whole set of traits, which seems to be an advantage of the integrative approach.


Preparatory screening of L and R values (values of the left and right traits) plays a significant role in testing FA of leaf plates, no less important than the preparation of herbarium material and improving the accuracy of measurement. First, the data were inspected for the presence of clearly asymmetric features; for example, if L and R differed more than 8-10 times, then such (unit) data were rejected. Then the normality of the distribution (L – R) was checked by the Kolmogorov – Smirnov test with the Lilliefors correction, since this difference is the starting point of the FA testing. The presence of antisymmetry was checked by the kurtosis of values (L-R), since the value of the kurtosis γ < -2 indicate a deviation from the normal distribution and the presence of antisymmetry as a possible form of bilateral asymmetry, which affects the value of FA. The directional asymmetry (DA) was determined by a t-test with the null hypothesis H0: (L – R) = 0.

Taking into account the possible relationship between the trait value (L + R)/2 and the fluctuating asymmetry, correlation analysis of pair values by Spearman was carried out (Palmer, Strobeck, 2003). Testing of statistical significance of the factors influencing fluctuation and plastic variability was carried out by factorial dispersion analysis. Plastic variability was determined using the formula PL = 1-x/X, where PL was the value of plastic variability, and x and X corresponded to the minimum and maximum values of the trait (L + R)/2 of the leaf plate from each tree (Bruschi, 2003). The average values of PL for each tree were used, i.e. the individual played the role of an experimental unit (n = 10). Most of the statistical analyses were carried out in STATISTICA10 (StatSoft Ink) using the level of statistical significance α = 0.05 %.


Weather conditions

 The 2017 vegetative season was colder and wetter. The average temperature in May-August 2017 was 16-18 % lower, and humidity 2-3 % higher compared to 2016 and 2018. According to the Russian research Institute of hydrometeorological information, May 2017 was the coldest since the beginning of the XX century, and June – the coolest in the last 70 years (table. 1).

Table 1. Climatic data, 2016–2018 

  2016 2017 2018 2016 2017 2018


May 60.30 60.21 57.96 15.03 10.91 16.19
June 62.36 69.20 58.87 18.19 14.45 16.32
July 69.54 73.51 72.63 20.96 17.94 20.47

72.92 69.82 63.51 19.48 18.84 19.64
Average 66.28 68.18 63.24 18.41 15.53 18.15

 Note. T – average air temperature Celsius at an altitude of 2 m above ground level; H – relative humidity in% at an altitude of 2 m above ground level.

The average air temperature in May was below the climatic norm by more than three degrees, temperature drops below zero and the appearance of snow in May and June were noted, which affects the phases of meristematic activity, with following budding, leafing and the beginning of Linden flowering in Russia (Kishchenko, 2015; Erofeeva, 2011; Ufimtseva, Terekhina, 2017).


Verification of preliminary data

The samples (L – R) deviated from the normal distribution, so to determine the fluctuating asymmetry the formula |Log L – Log R|was used. After normalization, only two samples did not satisfy the normal distribution conditions and showed directional asymmetry (t-test; p = 0.01). Considering the slight deviation from the criterion t, it was decided to include these data in the study. In 90% of samples (L – R) the kurtosis value γ was above zero, in the rest - the kurtosis was less than zero, but not less than -0.2. The critical value γ indicating antisymmetry was -0.68, when the number of observations n = 100, therefore, no antisymmetry was detected in the samples (L – R). Paired correlation coefficients, indicating the strength of the relationship between the magnitude of the trait and FA, showed a weak correlation dependency in Tilia cordata (r = -0.11 ÷ 0.1; p < 0.05). The weak correlation between the FA values of the traits allowed us to conclude that the traits were independent and suitable for determining development stability. The influence of various factors on the value of FA is shown in table. 2.


Table 2. Effect of species, year, gathering place and combined effect of factors on FA value 

Sources of variation






Species 0.161 1 0.161 149.77 0.000
Year of collection

0.016 2 0.008 7.62 0.001
Species x year of collection

0.012 2 0.006 5.79 0.003
Collecting place

0.003 1 0.003 2.74 0.098
Species x collecting place

0.005 1 0.005 4.31 0.038
Error 1.168 943 0.001

 Note. Conventional signs: SS – sum of squares; df – degree of freedom; MS – mean square; F – F-test; p – level of statistical significance


In fact, the value of FA depended on the species, the year of collection of leaf plates, but did not depend on the location of populations.

Fluctuation variability

As shown in Fig. 2, on the Boulevard, the FA value of only one species – cut-leaf linden– was significantly higher compared to the control.


Fig. 2. Dependence of FA on collecting place, F (1;945) = 4.31; p = 0.04. Vertical segments represent a 95 % confidence interval.


Thus, the cut-leaf linden had fluctuation variability explained by the action of 2 factors: the year of collection and collecting place. Observations of 2017 revealed critical FA values only for cut-leaf linden (Fig. 3).


Fig. 3. Dependence of FA on the year of collection. In the field of the diagram – average annual tempera­ture and humidity, average value for May – August. F (2;943) = 5.79; p = 0.003. Vertical segments represent a 95 % confidence interval



Plastic variability

Significant factors affecting the plastic variability of leaf plates were species and year of collection of leaf plates (table. 3).

    Table 3. Effect of some factors on phenotypic plasticity

Sources of variation






Species 0.037 1 0.037 12.351 0.001
Year of collection

0.053 2 0.027 8.780 0.000
Collecting place

0.001 1 0.001 0.403 0.527
Species x year of collection

0.039 2 0.019 6.412 0.003
Species x collecting place

0.016 1 0.016 4.501 0.037
Error 0.58 178 0.007    

    Note. Conventional signs: SS – sum of squares; df – degree of freedom; MS – mean square; F – F-test; p – level of statistical significance.


Plastic variability was not the same in the two species (df = 62; p = 0.0002) and depended on year of collection (df = 62;    p = 0.006). In the big-leaf linden indicators PV, as well as FA, were higher than in small-leaved linden (p < 0.05; Fig. 4).


Fig. 4.  Plastic variability, F (2; 88) = 6.41; p = 0.003. Vertical segments mean a 95 % confidence interval


2017 led to an increase in FA only in big-leaf linden, this year the leaf plate was significantly smaller: 44.34 ± 0.83 – in 2017; 48.74 ± 1.04 - in 2016 and 46.85 ± 0.62-in 2018 (p < 0.001). There was no statistically significant correlation between PV and the mean value of the trait. Characteristics of the leaf plate of small-leaved linden remained invariant during the entire observation period.


Verification of results

 Since the leaf plates were unequal in size (about 30 % of samples), the verification task was to determine whether the standard error ME in determining the FA significantly changed depending on the value of the trait. If the FA–trait relationship was stronger than the error–trait relationship, this gave the right to positively verify the statistical significance of the FA value (Palmer, Strobeck, 2003). A regression analysis was carried out to estimate the Pearson correlation coefficient, since the samples containing the average values per tree had normal distribution. The ratio of the difference between correlation coefficients to the measurement error of FA was tested using a t-test based on the assumption of the normality of the distribution of such a ratio:

ts = (rFA – rME)/ SEFA,

where ts is statistical test of differences of correlation coefficients; rFA– value of r correlation between FA and (R+ L)/2; rMEr correlation between FA and ME; SEFA – standard error of the FA.

 For both species, the statistical value of the ts criterion exceeded the tabular criterion value ts = 2.02; n = 40 (small-leaved linden – 62.4, big-leaf linden– 43.4). Thus, the effect of measurement error on the value of FA at the level of individual variability was insignificant. The generalized linear model (GLM) showed no difference in the obtained dependence in the samples grouped by year of collection or by place of collection (p > 0.05), i.e. in different years and in different populations, the measurement error did not affect the FA equally.The authors consider the described verification of the statistical significance of the standard error to be fundamental in the case of working with samples of plates of different sizes, which is often found in such studies.


The choice of individual traits to determine the stability of development in contrast to the integrative path is questionable and can be justified only if there is an accurate assessment of the gradient of factors on the value of FA of the trait (Baranov, 2014b; Baranov et al., 2015, Zykov et al., 2015). In our case, normalization with the use of logarithmation allowed to operate with normally distributed values in all 5 traits, which can be considered, rather, as luck. At the same time, an intrigue to use separate traits (one or two) the most sensitive to pollutants remains (Baranov, 2014a; Koroteeva et al., 2015). It is known that the plasticity of plants is caused by the buffer capacity of morphological structures, which allows them to adapt actively to environmental conditions (Debat, David, 2001; Klingenberg, 2003; Tonsor et al., 2013; Venâncio et al., 2016). Consequently, stressful climatic effects and chemical pollution on the Boulevard affected the PV, adapting the plant. It is known that on the Kola Peninsula in the cold summer there was an increase in the FA of the leaves of the white birch, while the temperature was a more significant factor than pollution with sulfur dioxide and nickel (Valkama, Kozlov, 2001). As in the case of birch, the climatic factor played a more important role than the action of chemical pollutants.

Big-leaf linden is a more thermophilic species with more pronounced stenobionticy to climatic factors and physical and to chemical properties of the soil. We see the reason for the increase in FA of leaf plates in the decrease in the growth of plates and the decrease the photosynthesis activity during the cold summer. This is confirmed by the indicator values of thermophilic big-leaf linden, which are shifted to the right positive side (8-13) compared to small-leaf linden (6-12 on the Tsyganov scale).

The authors attribute the adaptation of big-leaf linden to the conditions of the Central Russian plain to one of the phases of the reaction to stress with an increase in the FA of leaf plates. Consequently, the hypothesis of the relationship of 2 types of variability was confirmed. We consider the specific genotypic feature of this introducer in the form of increased plastic and fluctuation variability as a "fee" for adaptation in the condition of climatic stress factors.


The phenotypic plasticity influences fluctuation variability (Houle, 2000; Sultan, 2003; Scheine, 2004; Lajus, Alekseev, 2003; Tucić et al., 2018). Stability of development is a category that explains the processes caused by the peculiarities of the genotype and the environment. PV is an optional indicator correlated with FA, but a frequent manifestation accompanying FA. In this case, the increased plastic variability was a phenotypic deviation caused by low air temperatures, while the increased fluctuating asymmetry was due to the small size of leaf plates. The growing plate is known to develop intermittently in the right, then in the left direction and has increased asymmetry (Freeman et al., 1993). It is known that the oscillation of the variance of the difference (L-R) within the same value of the trait should be the attribute of destabilization of development (Palmer, Strobeck, 2003). In our case, the plasticity of the development of leaf plates was increased, however, FA, although increased, was not evidence of a decrease in the stability of development as a deviation in the homeostasis of development at the biochemical level (Klingenberg, 2003; Graham et al., 2010; Klingenberg, 2016).

Increased plastic variability in the leaf plates of the big-leaved linden as a species-specific property, without imbalance in the homeostasis of development, should be taken into account in studies related to long-term monitoring of the stability of the development of woody plants.


Takahashi K. H. Multiple modes of canalization: links between genetic, environmental canalizations and developmental stability, and their trait-specificity, Seminars Cell & developmental biology. 2018. May 22. DOI: 10.1016/j.semcdb.2018.05.018.

Abley K., Locke J. C., Leyser H. O. Developmental mechanisms underlying variable, invariant and plastic phenotypes, Annals of botany. 2016. Apr. 1. Vol. 117 (5). P. 733–748.

Alados C. L., Navarro T., Escós J., Cabezudo B., Emlen J. M. Translational and fluctuating asymmetry as tools to detect stress in stress-adapted and nonadapted plants, International Journal of Plant Sciences. 2001. May. Vol. 162 (3). P. 607–616.

Anastasia Z. The formation of fluctuating asymmetry during individual development of Vetula pendula, Principy Èkologii. 2014. Vol. 4. Jan. 1. P. 31–52.

Baranov S. G. Zykov I. E. Fedorova L. V. Investigation of intraspecific variability of little-leaved linden (Tilia cordata Mill.) on the basis of bilateral asimmetry of laminas, Vestnik Tomskogo gosudarstvennogo universiteta. Biologiya. 2015. No. 2 (30) C. 134–145.

Baranov S. G. Littleleaf Linden Tilia cordata (Mill.): Only Some Bilateral Traits Indicate Chemical Pollution Induced by Chemical Plant, Advances in Biological Research. 2014a. Vol. 8 (4). P. 143–148.

Baranov S. G. Use of MorphoJ Soft Package for Phenotypic and Genotypic Variety Testing (English Oak Case Study), Russian Journal of Genetics: Applied Research. 2018. Vol. 8. No 1. P. 22–30.

Baranov S. G. Use of a geometric morphometric method to determine the developmental stability of Betula pendula Roth., Biology Bulletin. 2017. Sep. 1. Vol. 44 (5). P. 547–551.

Baranov S. G. Use of morphogeometric method for study fluctuating asymmetry in leaves Tilia cordata under industrial pollution, Adv. Environ. Biol. 2014b. Vol. 8 (7). P. 2391–2398.

Belyaeva Yu. V. Indicators of fluctuating asymmetry Betula pendula Roth. in conditions of anthropogenic impact, Izvestiya Samarskogo nauchnogo centra Rossiyskoy akademii nauk. 2013. No. 15. P. 3–7.

Bruschi, P., Grossoni P., Bussotti F. Within-and among-tree variation in leaf morphology of Quercus petraea (Matt.), Liebl. natural populations. Trees. 2003. Vol. 17 (2). P. 164–172.

Clarke G. M. Fluctuating asymmetry of invertebrate populations as a biological indicator of environmental quality, Environmental Pollution. 1993. Jan. 1. Vol. 82 (2). P. 207–211.

Cornelissen T., Stiling P., Drake B. Elevated CO2 decreases leaf fluctuating asymmetry and herbivory by leaf miners on two oak species, Global Change Biology. 2004. Jan. 1. Vol. 10 (1). P. 27–36. DOI: 10.1046/j.1529-8817.2003.00712.x.

De Kroon H., Huber H., Stuefer J. F., and Van Groenendael J. M. A modular concept of phenotypic plasticity in plants, New Phytol. 2005. Vol. 166. P. 73–82. DOI: 10.1111/j.1469-8137.2004.01310.x.

Debat V., David P. Mapping phenotypes: canalization, plasticity and developmental stability, Trends in Ecology & Evolution. 2001. Oct. 1. Vol. 16 (10). P. 555–561. DOI: org/10.1016/S0169-5347(01)02266-2.

Erofeeva E. A. Effect of motor vehicle pollution on the rate of exit from winter dormancy and termination of vegetation in little-leaved linden, Vestnik NNGU. 2011. No. 2-2. URL: (data obrascheniya: 25.08.2018).

Erofeeva E. A. Hormesis and paradoxical effects of Drooping Birch (Betula pendula Roth) parameters under motor traffic pollution, Dose-Response. 2015. Jun. 1. Vol. 13 (2). DOI: 10.1177/1559325815588508.

Freeman D. P., Graham J. H., Emlen J. M. Developmental stability in plants: symmetries, stress and epigenesist, Genetica. 1993. No 89. P. 97–119.

Gilbert S. F., Bosch T. C., Ledón-Rettig C. Eco-Evo-Devo: developmental symbiosis and developmental plasticity as evolutionary agents, Nature Reviews Genetics. 2015. Oct. Vol. 16 (10). P. 611.

Gostin I. Air pollution stress and plant response, Plant Responses to Air Pollution. 2016. P. 99–117.

Graham J. H., Raz S., Hel-Or H., Nevo E. Fluctuating Asymmetry: Methods, Theory, and Applications, Symmetry. 2010. Vol. 2. P. 466–540. DOI: 10.3390/sym2020466.

Hódar J. A. Leaf fluctuating asymmetry of Holm oak in response to drought under contrasting climatic conditions, Journal of arid environments. 2002. Oct. 1. Vol. 52 (2) P. 233–243.

Hochwender C. G., Fritz R. S. Fluctuating asymmetry in a Salix hybrid system: the importance of genetic versus environmental causes, Evolution. 1999. Apr. Vol. 53 (2). P. 408–416.

Houle D. A simple model of the relationship between asymmetry and developmental stability, Journal of Evolutionary Biology. 13.4.2000. P. 720–730. Instab.: Causes and Consequences. Polak M, editor. Oxford University Press, 2003. DOI: 10.1046/j.1420-9101.2000.00195.x.

Ivanov V. P., Ivanov Y. V., Marchenko S. I., Kuznetsov V. V. Application of fluctuating asymmetry indexes of silver birch leaves for diagnostics of plant communities under technogenic pollution, Russian journal of plant physiology. 2015. May 1. Vol. 62 (3). P. 340–348.

Kaligarič M., Tognetti R., Janžekovič F., Raschi A. Leaf Fluctuating Asymmetry of Myrtus Communis L., Affected by Increases in Atmospheric CO2 Concentration: Evidence from a Natural CO2 Spring, Polish Journal of Environmental Studies. 2008. Jul. 1. Vol. 17 (4).

Kischenko I. T. Tilia L. Growth and development of introduced species of genus Tilia L. (tiliaceae) in conditions of Karelia, Arctic Evironmental Research. 2015. No. 4. URL: (data obrascheniya: 25.08.2018).

Klingenberg C. P. A developmental perspective on developmental instability: theory, models and mechanisms, Developmental instability: causes and consequences. 2003. P. 14–34.

Klingenberg C. P. Size, shape, and form: concepts of allometry in geometric morphometrics, Development genes and evolution. 2016. Jun. 1. Vol. 226 (3). P. 113–137. DOI: 10.1007/s00427-016-0539-2.

Klingenberg, C. P., Duttke S., Whelan S., Kim M. Developmental plasticity, morphological variation and evolvability: A multilevel analysis of morphometric integration in the shape of compound leaves, J. Evol. Biol. 2012. Vol. 25. P. 115–129.

Koroteeva E. V., Veselkin D. V., Kuyantseva N. B., Chashchina O. E. The size, but not the fluctuating asymmetry of the leaf, of silver birch changes under the gradient influence of emissions of the Karabash Copper Smelter Plant, Doklady Biological Sciences. 2015. Vol. 460. No 3. P. 364–367. DOI: 10.1134/S0012496615010056.

Kozlov M. V., Cornelissen T., Gavrikov D. E., Kunavin M. A., Lama A. D., Milligan J. R., Zverev V., Zvereva E. L. Reproducibility of fluctuating asymmetry measurements in plants: Sources of variation and implications for study design, Ecological indicators. 2017. Feb. 1. Vol. 73. P. 733–740/ DOI: 10.1016/j.ecolind.2016.10.033

Kozlov M. V., Niemelä P. Drought is more stressful for Northern populations of Scots pine than low summer temperatures, Silva Fennica. 2003. Vol. 37. P. 175–180.

Kozlov M. V., Zverev V., Zvereva E. L. Do defoliating insects distinguish between symmetric and asymmetric leaves within a plant?, Ecological Entomology. 2018. Vol. 43 (5). P. 656–664.

Kozlov M. V., Zvereva E. L. Confirmation bias in studies of fluctuating asymmetry, Ecological indicators. 2015. Oct. 1. Vol. 57. P. 293–297. DOI: 10.1016/j.ecolind.2015.05.014.

Lajus D. L., Alekseev V. R. Phenotypic variation and developmental instability of life-history traits: a theory and a case study on within-population variation of resting eggs formation in Daphnia, Journal of Limnology. 2004. Sep. 1. Vol. 63 (1). P. 37–44. DOI: 10.4081/jlimnol.2004.s1.37.

Lande R. Adaptation to an extraordinary environment by evolution of phenotypic plasticity and genetic assimilation, J. Evol. Biol. 2009. Vol. 22. P. 1435–1446. DOI: 10.1111/j.1420-9101.2009.01754.x.

Lens L. U., Van Dongen S., Kark S. et al. Fluctuating asymmetry as an indicator of fitness: can we bridge the gap between studies?, Biological Reviews. 2002. Feb. Vol.77 (1). P. 27–38. DOI: 10.1017S1464793101005796.

Leung B., Forbes M. R. Fluctuating asymmetry in relation to stress and fitness: effects of trait type as revealed by meta-analysis, Ecoscience. 1996. Jan 1. Vol. 3 (4). P. 400–413. DOI: 10.1080/11956860.1996.11682357.

Leung B., Forbes M. R., Houle D. Fluctuating asymmetry as a bioindicator of stress: comparing efficacy of analyses involving multiple traits, The American Naturalist. 2000. Jan. Vol. 155 (1). P. 101–115.

Palmer A. R., Strobeck C. H. Fluctuating Asymmetry Analyses Revisited, Developmental Instab.: Causes and Consequences. Polak M, editor. Oxford University Press, 2003. P. 279–319.

Royer D. L., Meyerson L. A., Robertson K. M., and Adams J. M. Phenotypic plasticity of leaf shape along a temperature gradient in Acer rubrum, PLoS ONE. 2009. Vol. 4(10). DOI: 10.1371/journal.pone.0007653.

Scheiner S. M. The genetics of phenotypic plasticity. XIII. Interactions with developmental instability, Ecology and evolution. 2014. Apr. 1. Vol. 4 (8). P. 1347–1360. DOI: 10.1002/ece3.1039.

Sobchak R. O. Afanas'eva T. G. Kopylov M. A. Estimation of the ecological state of recreational zones by the method of fluctuating asymmetry of leaves Betula pendula Roth., Vestnik Tomskogo gosudarstvennogo universiteta. 2013. No. 368. P. 195–199.

Sultan S. E. Phenotypic plasticity in plants: a case study in ecological development, Evolution & development. 2003. Jan. 1. Vol. 5 (1). P. 25–33. DOI: 10.1046/j.1525-142X.2003.03005.x.

Tonsor S. J., Elnaccash T. W., Scheiner S. M. Developmental instability is genetically correlated with phenotypic plasticity, constraining heritability, and fitness, Evolution. 2013. Oct. 1. Vol. 67 (10). P. 2923–2935. DOI: 10.1111/evo.12175.

Tucić B., Budečević S., Manitašević Jovanović S. et al. Phenotypic plasticity in response to environmental heterogeneity contributes to fluctuating asymmetry in plants: first empirical evidence, Journal of Evolutionary Biology. 2018. Vol. 31 (2). P. 197–210. DOI: 10.5061/dryad.8th5m.

Ufimceva M. D. Terehina N. V. Estimation of ecological state of the central region of Saint-Petersburg on the basis of ecophitoindication, Vestnik SPbGU. Nauki o Zemle. 2017. No. 2. URL: (data obrascheniya: 25.08.2018).

Valkama J., Kozlov M. V. Impact of climatic factors on the developmental stability of mountain birch growing in a contaminated area, Journal of Applied Ecology. 2001. Vol. 38. P. 665–673.

Veličković M. V. Reduced developmental stability in Tilia cordata leaves: effects of disturbed environment, Periodicum biologorum. 2010. Sep. 1. Vol. 112 (3). P. 273–281. P. 296–303.

Venâncio Henrique, Estevao Alves-Silva and Jean Carlos Santos. Leaf phenotypic variation and developmental instability in relation to different light regimes, Acta Botanica Brasilica. 2016. Vol.30 (2). P. 296–303. DOI: 10.1590/0102-33062016abb0081

Wuytack T., Wuyts K., Van Dongen S., Baeten L., Kardel F., Verheyen K., Samson R. The effect of air pollution and other environmental stressors on leaf fluctuating asymmetry and specific leaf area of Salix alba L., Environmental pollution. 2011. Oct. 1. Vol. 159 (10). P. 2405–

Zaharov V. M. Chubinishvili A. T. Environmental health monitoring in protected natural areas. M.: Centr ekologicheskoy politiki Rossii, 2001. 148 p.

Zorina A. A. Korosov A. V. Characteristic of the fluctuation asymmetry of leaf of two species of birches in Karelia, Trudy Karel'skogo nauchnogo centra Rossiyskoy akademii nauk. 2007. Vyp. 11. P. 28–36

Zverev V., Lama A. D., Kozlov M. V. Fluctuating asymmetry of birch leaves did not increase with pollution and drought stress in a controlled experiment, Ecological indicators. 2018. Jan. 31. Vol. 84. P. 283–289.DOI: 10.1016/j.ecolind.2017.08.058

Zvereva E. L., Kozlov M. V., Haukioja E. Stress responses of Salix borealis to pollution and defoliation, Journal of Applied Ecology. 1997. Vol. 34 (6) Dec. 1. P. 1387–1396. DOI: 10.2307/2405256.

Zvereva E. L., Kozlov M. V., Niemelä P., Haukioja E. Delayed induced resistance and increase in leaf fluctuating asymmetry as responses of Salix borealis to insect herbivory, Oecologia. 1997. Feb. 1. Vol. 109 (3). P. 368–373.

Zykov I. E. Fedorova L. V. Baranov S. G. Assessment of the biological value of the level of variability of the parameters of leaf blades of small-leaved linden (Tilia cordata Mill.) in the city of Orekhovo-Zuevo and Orekhovo-Zuevo region, Vestnik Moskovskogo gosudarstvennogo oblastnogo universiteta. Seriya: Estestvennye nauki. 2015. No. 1. P. 15–21.

Displays: 5056; Downloads: 914;