Research Article
Research Article
Light intensity affects leaf morphology in a wild population of Adenostyles alliariae (Asteraceae)
expand article infoBailey Francis, Robert Tucker Gilman
‡ University of Manchester, Manchester, United Kingdom
Open Access


Low light conditions can impose environmental stress on plants, and plants often respond adaptively by increasing their leaf area. Light stress on plants can also result in developmental instability, which can manifest as increased fluctuating asymmetry in leaves or other organs. The relationship between light conditions and fluctuating asymmetry has been documented in experimental populations, but has been less frequently observed in the wild. Here, we studied how leaf surface area and fluctuating asymmetry correlate with light intensity in a wild population of Adenostyles alliariae (Asteraceae). We found strong evidence that leaf surface area increases and weak evidence that fluctuating asymmetry increases as light intensity decreases. Our results help to elucidate the relationship between light stress and developmental instability under naturally occurring conditions.


fluctuating asymmetry, light intensity, developmental instability, phenotypic plasticity, Adenostyles alliariae


Plant growth and survival are strongly influenced by abiotic factors such as light intensity, nutrient availability, and temperature (Agarwal et al. 2006). These abiotic factors can produce conditions that are unfavourable for normal plant development (Jan et al. 2012). Plants living in shaded or poor nutrient habitats or in areas outside their optimal temperature ranges can display increased levels of stress compared to plants living in more favourable environments (Sridevi et al. 1999; Lynch and Brown 2001; Yang et al. 2019). To survive in these stressful environments, plants have evolved the ability to respond to abiotic stressors with physiological and morphological alterations (Alpert and Simms 2002), resulting in intraspecific phenotypic variation (Miner et al. 2005; Rozendaal et al. 2006).

Plants require light for their development and metabolism, but they sometimes occur in areas where light availability is below the optimum level (i.e., light stress, Niinemets 2010; Cramer et al. 2011; Zhu et al. 2017). In these conditions, plants often respond by increasing their leaf area to intercept more sunlight (Dwyer et al. 2014; Liu et al. 2016), although this response can vary even among closely related species (Sultan 2003). As a result, plants growing in different natural light conditions are expected to have different leaf sizes (Marques et al. 1999).

Another consequence of light stress for plants is a phenomenon known as developmental instability (DI). DI occurs when an organism is unable to achieve its developmentally programmed phenotype, and can result from unfavourable environmental conditions (Markow 1995). One common manifestation of DI is fluctuating asymmetry (FA), or small deviations from perfect symmetry in otherwise symmetrical structures (Moller and Shykoff 1999; Cornelissen and Stiling 2010). These deviations are caused when organisms are exposed to severe stress at certain critical development stages, resulting in abnormal phenotypes and increased phenotypic variability among individuals (Waddington 1956; Graham et al. 1993). FA can reflect disturbances in developmental homeostasis at the molecular, chromosomal or epigenetic levels (Parsons 1990). Because FA is easily measured in the field, it is often used as a surrogate for DI in ecological studies (Palmer and Strobeck 1986; Graham et al. 2010). In studies of DI resulting from light stress in plants, researchers have often focused on FA in the leaves (Midgley et al. 1998; Roy and Stanton 1999; Venâncio et al. 2016).

If leaf FA is greater at low light intensity, this might be because low light intensity causes developmental stress. Alternatively, the relationship might be mediated by leaf growth. If leaves in low light conditions are larger than leaves in high light conditions, then they must grow faster or for longer, and faster or longer growth might result in DI (Moller and Eriksson 1994; Cuevas-Reyes et al. 2011). These potential drivers of FA have not been fully disentangled.

In the Carnic Alps, the herbaceous perennial Adenostyles alliariae (Asteraceae) is a useful model for investigating DI in response to light conditions. The species is common and grows in both coniferous forests and alpine meadows, and therefore occupies habitats with a wide variety of light conditions. Phenotypic variation and DI have not yet been examined in this species. Additionally, light stress is one of the most uncharacterized and least studied abiotic stresses that plants encounter (Yang et al. 2019). In this study we investigated both phenotypic variation and DI of A. alliariae leaves in areas of differing sunlight exposure. We measured phenotypic variation using leaf area, and used leaf FA as a surrogate for DI (Santos et al. 2013). We hypothesised that leaves growing in low light conditions would have increased area and increased FA, which may indicate phenotypic compensation for low light conditions and increased levels of stress in A. alliariae, respectively.


Site selection

Fieldwork was conducted at the Baita Torino field station in the Passo del Pura, Ampezzo, Friuli-Venezia Giulia, Italy (46°25.5433'N, 12°44.5167'E (DDM), 1400 m asl) between 11:00 and 13:00 on 9, 11 and 12 July 2019. The area is part of the Carnic Alps, which are characterised by high plant biodiversity and endemism (Pignatti and Pignatti 2014). We selected four study sites in the vicinity of the Baita Torino. Two sites (i.e., the shaded sites) were covered by coniferous forest (primarily Abies alba and Picea abies) and two sites (i.e., the open sites) were unforested. Sites were selected for their coverage properties and for the presence of A. alliariae. Thus, site selection was not random.

Data collection

At each site we randomly selected 25 individual A. alliariae plants with no signs of herbivory for study (i.e., 100 plants in total). We identified the most basal leaf on each plant and we measured the light intensity at its apex and base, and on the right and left side of the leaf at its widest points using the Google Science Journal application (version 3.2) on an iPhone 7. We averaged these measurements to estimate the light intensity reaching the leaf. We then collected the leaf and photographed it on a white background with a 25 mm scale bar. We used ImageJ (version 1.51) to measure the adaxial surface area of each leaf. We also measured the distance from the midrib to the widest points on the left and on the right sides of the leaf. We called these distances Ls and Rs, respectively (Figure 1). We estimated the directed asymmetry, d, of each leaf using d = log(Ls) – log(Rs) (Waddington 1960; Venâncio et al. 2016). Thus, d > 0 if Ls > Rs and d < 0 if Ls < Rs. This measure is unitless and depends on the shape but not on the size of the leaf. We estimated the undirected asymmetry of each leaf as |d| (Waddington 1960; Venâncio et al. 2016). The undirected asymmetry captures the magnitude but not the direction of asymmetry in each leaf.

Figure 1. 

The most basal leaf from an A. alliariae plant sampled in this study. Ls and Rs represent the measures for the midrib to the widest points on left and right of the leaf, respectively.


Asymmetry in a population can be fluctuating or directional. Directional asymmetry occurs if one side of an organ or organism is consistently different from the other in all members of a population, and need not indicate DI (Graham et al. 1993). For example, in Tengmalm’s owl (Aegolius funereus), the right ear aperture is positioned higher on the head than the left, and this is believed to enhance directional hearing (Norberg 1978). In our study, asymmetry would be directional if the left side of the leaf were consistently larger than the right side, or vice versa. To test for directional asymmetry, we applied a one-sample t-test to ask whether d was significantly different from zero.

In the absence of directional asymmetry, we can interpret the undirected asymmetry of leaves in our study as FA. We regressed the logged leaf area and FA on the logged light intensity reaching each leaf using mixed effects models implemented in the R package lme4 (Bates et al. 2015). The models included random effects of site (to account for potential differences in site qualities such as soil chemistry) and day (to account for systematic differences in light measurements due to weather on different days). The residuals of the fitted models were not normally distributed, so we Box-Cox transformed logged leaf area and FA (Box and Cox 1964; Venâncio et al. 2016) using the powerTransform function in the R package car (Fox and Weisberg 2019). Box-Cox transformation requires that data to be transformed be strictly positive, and our FA measurements included two zeros. Therefore, we Winsorised the zero values to the smallest non-zero value in the data set (i.e., the smallest FA we know we could detect, Tukey 1962) prior to Box-Cox transformation. Then, we refit the mixed effects models to the transformed data. We calculated p-values for the effects of logged light intensity on logged leaf area and FA using the Satterthwaite approximation (Satterthwaite 1941; Luke 2017) implemented in the R package lmerTest (Kuznetsova et al. 2017).

If the relationship between FA and light intensity is mediated by leaf growth, then under the same light conditions we would expect leaves that grow faster or for longer (i.e., larger leaves) to show more FA. Thus, we would expect to see a relationship between leaf area and FA after controlling for light intensity. To test this, we regressed Box-Cox transformed FA on logged light intensity and logged leaf area in a mixed model with random effects of site and day, and estimated p-values using the Satterthwaite approximation.

Data are available in Suppl. material 1: File S1: aa_data.csv, and the R code used to analyze the data is available in Suppl. material 1: File S2: Francis_and_Gilman_analysis.R. Both files are also archived in the TRY Plant Trait Database (Kattge et al. 2011) at


The light intensity measured at individual leaves varied within sites (Figure 2), but was higher overall at open sites (8.26 ± 0.22 log(lx)) than at shaded sites (6.75 ± 0.25 log(lx)) (Monte Carlo randomisation, p < 0.0001). The mean logged leaf area was 9.53 ± 0.07 log(mm2) in open sites and 9.79 ± 0.07 log(mm2) in shaded sites. The mean undirected asymmetry was 0.069 ± 0.005 in open sites and 0.080 ± 0.008 in shaded sites.

Figure 2. 

Light intensity measured at individual leaves in open and shaded sites.

Adaxial surface area increased with decreasing light intensity (p < 0.0001, b = 0.27 log(mm2) log(lx)-1 calculated at the mean-centered logged light intensity of 7.50 log(lx); Figure 3A). We found no evidence that asymmetry in the study population was directional rather than fluctuating (p = 0.2032) and thus we treated undirected asymmetry as FA. We found weak evidence that FA increased with decreasing light intensity (p = 0.0730, b = 4.8 x 10-3 log(lx)1 calculated at the mean-centered logged light intensity; Figure 3B). If we forced site type (i.e., open or shaded) into the model, we found no relationship between site type and either logged leaf area (p = 0.70) or FA (p = 0.99). Thus, leaf morphology is better explained by light intensity at the individual leaf than by light intensity at the study site. We found no relationship between FA and logged leaf area after controlling for the effect of light intensity (p = 0.532). Therefore, we have no evidence that the effect of light intensity on FA was mediated by leaf growth.

Figure 3. 

Effects of light intensity on adaxial surface area A and FA B in A. Alliariae leaves. Each point represents the most basal leaf of a plant growing in a shaded (filled circles) or open (open circles) site. The solid line indicates a significant relationship (p < 0.0001) and the dashed line indicates a marginally significant relationship (p = 0.0732). Lines of best fit are back-transformed from models fitted to Box-Cox transformed data.


Plants can respond to environmental stresses by modifying their physiology and morphology, and this is reflected in their phenotypic variation (Miner et al. 2005; Rozendaal et al. 2006). Such modifications are common in plants growing in conditions of poor light availability (Markesteijn et al. 2007). In this study, we found that A. alliariae growing at lower light intensity produce larger leaves. A plastic response of leaf size to light intensity is likely to be adaptive for plants that experience variable light conditions between generations or across small spatial scales (Markesteijn et al. 2007).

We found weak evidence that leaf FA increases with decreasing light intensity in A. alliariae. The relationship between light intensity and FA has found mixed support in the literature. Several studies have reported increased leaf FA in plants growing at low light intensity under manipulated (Sinapis arvensis Roy & Stanton, 1999; Quercus pyrenaica Puerta-Pinero et al., 2008; Silene vulgaris Sandner & Matthies, 2017) or natural conditions (Quercus alba Kusi, 2013; Miconia fallax Alves-Silva, 2012). In contrast, Wuytack and colleagues (2011) found no effect of light intensity on leaf FA in Salix alba and Venâncio and colleagues (2016) found no effect of light intensity on leaf FA in Bauhinia brevipes, albeit in a study of only 26 plants versus the 100 plants in our study. Other studies have found that FA increases with increasing light intensity. For example, Midgley and colleagues (1998) found that leaf FA in Dimorphotheca sinuate increased with UV-B exposure, although their UVB exposure treatments simulated depleted ozone conditions and were therefore greater than plants regularly encounter. Tucic and Miljkovic (2010) found increased FA in Iris pumila growing higher light conditions in nature, but in flowers rather than leaves. To our knowledge, our study provides the first evidence of increasing FA with decreasing light intensity in a naturally occurring population of an herbaceous plant. Given the variability of reported results, systematic work to identify the ecological traits or genotypes that mediate the relationship between leaf FA and light intensity might reward effort.

Some authors have suggested that the relationship between low light intensity and increased leaf FA might be mediated by leaf growth (Moller and Eriksson 1994; Cuevas-Reyes et al. 2011). That is, leaves grow to greater size under low light conditions, and the demands of longer or faster growth might result in stress and thus developmental instability. However, if that were the case, we would expect larger leaves to show more FA under similar light conditions, and we did not see this pattern in our study. Many environmental conditions other than light can also affect leaf FA, and some of these conditions may have varied among sites in our study. However, we controlled for site-specific environmental variation by including site as a random effect in our analyses. Therefore, we believe that differences in light intensity are likely to explain the differences we observed in leaf FA among A. alliariae plants.

Our study examined only one leaf per plant. Therefore, we cannot determine whether light intensity at each leaf affects the morphology of that leaf, or whether the light intensity experienced by each plant affects the morphology of all leaves on the plant and the light intensity at the most basal leaf is correlated with the light intensity experienced by the plant as a whole. These two possibilities are not mutually exclusive, and the answer need not be the same for leaf area and FA. Studies that assess the light intensity at and morphology of multiple leaves per plant could answer this question. To our knowledge, this analysis has rarely been attempted for FA (but see Kusi 2013).

A broad understanding of how light intensity affects developmental instability in plants remains to be achieved. By linking light intensity to fluctuating asymmetry in a naturally occurring population, albeit weakly, we believe our study offers a step toward that goal.


The authors thank Prof. Pier Luigi Nimis and James Dinsley for help with plant identification; Mariana Villalba de la Pena, Estevao Alvez-Silva and an anonymous reviewer for advice; and Natural Environment Research Council grant NE/L002469/1 for support.


  • Agarwal P, Agarwal P, Reddy M, Sopory S (2006) Role of DREB transcription factors in abiotic and biotic stress tolerance in plants. Plant Cell Reports 25(12): 1263–1274.
  • Alpert P, Simms EL (2002) The relative advantages of plasticity and fixity in different environments: when is it good for a plant to adjust? Evolutionary Ecology 16: 285–297.
  • Alves-Silva E (2012) The influence of Ditylenchus (Nematoda) galls and shade on the fluctuating asymmetry of Miconia fallax (Melastomataceae). Ecologia Austral 22: 53–61.
  • Cornelissen T, Stiling P (2010) Small Variations over Large Scales: Fluctuating Asymmetry over the Range of Two Oak Species. International Journal of Plant Sciences 171(3): 303–309.
  • Cramer G, Urano K, Delrot S, Pezzotti M, Shinozaki K (2011) Effects of abiotic stress on plants: a systems biology perspective. BMC Plant Biology 11(1): 1–163.
  • Cuevas-Reyes P, Oyama K, Gonzalez-Rodriguez A, Wilson-Fernandes G (2011) Contrasting herbivory patterns and leaf fluctuating asymmetry in Heliocarpus pallidus between different habitat types within a Mexican tropical dry forest. Journal of Tropical Ecology 27(4): 383–391.
  • Dwyer JM, Hobbs RJ, Mayfield MM (2014) Specific leaf area responses to environmental gradients through space and time. Ecology 95(2): 399–410.
  • Fox J, Weisberg S (2019) An {R} Companion to Applied Regression (3th edn.). Sage.
  • Kattge J, Díaz S, Lavorel S, Prentice IC, Leadley P, Bönisch G, Garnier E, Westoby M, Reich PB, Wright IJ, Cornelissen JHC, Violle C, Harrison SP, Van Bodegom PM, Reichstein M, Enquist BJ, Soudzilovskaia NA, Ackerly DD, Anand M, Atkin O, Bahn M, Baker TR, Baldocchi D, Bekker R, Blanco CC, Blonder B, Bond WJ, Bradstock R, Bunker DE, Casanoves F, Cavender‐Bares J, Chambers JQ, Chapin III FS, Chave J, Coomes D, Cornwell WK, Craine JM, Dobrin BH, Duarte L, Durka W, Elser J, Esser G, Estiarte M, Fagan WF, Fang J, Fernández‐Méndez F, Fidelis A, Finegan B, Flores O, Ford H, Frank D, Freschet GT, Fyllas NM, Gallagher RV, Green WA, Gutierrez AG, Hickler T, Higgins SI, Hodgson JG, Jalili A, Jansen S, Joly CA, Kerkhoff AJ, Kirkup D, Kitajima K, Kleyer M, Klotz S, Knops JMH, Kramer K, Kühn I, Kurokawa H, Laughlin D, Lee TD, Leishman M, Lens F, Lenz T, Lewis SL, Lloyd J, Llusià J, Louault F, Ma S, Mahecha MD, Manning P, Massad T, Medlyn BE, Messier J, Moles AT, Müller SC, Nadrowski K, Naeem S, Niinemets Ü, Nöllert S, Nüske A, Ogaya R, Oleksyn J, Onipchenko VG, Onoda Y, Ordoñez J, Overbeck G, Ozinga WA, Patiño S, Paula S, Pausas JG, Peñuelas J, Phillips OL, Pillar V, Poorter H, Poorter L, Poschlod P, Prinzing A, Proulx R, Rammig A, Reinsch S, Reu B, Sack L, Salgado‐Negret B, Sardans J, Shiodera S, Shipley B, Siefert A, Sosinski E, Soussana J-F, Swaine E, Swenson N, Thompson K, Thornton P, Waldram M, Weiher E, White M, White S, Wright SJ, Yguel B, Zaehle S, Zanne AE, Wirth C (2011) TRY – a global database of plant traits. Global Change Biology 17(9): 2905–2935.
  • Kusi J (2013) Variations in Phenotypic Plasticity and Fluctuating Asymmetry of Leaf Morphology of Three Quercus (Oak) Species in Response to Environmental Factors. MS Thesis, East Tennessee State University, USA.
  • Kuznetsova A, Brockhoff PB, Christensen RHB (2017) lmerTest package: tests in linear mixed effects models. Journal of Statistical Software 82(13): 1–26.
  • Liu Y, Dawson W, Prati D, Haeuser E, Feng Y, Kleunen MV (2016) Does greater specific leaf area plasticity help plants to maintain a high performance when shaded? Annals of Botany 118(7): 1329–1336.
  • Markesteijn L, Poorter L, Bongers F (2007) Light-dependent leaf trait variation in 43 tropical dry forest tree species. American Journal of Botany 94: 515–525.
  • Marques AR, Garcia QS, Fernandes GW (1999) Effects of sun and shade on leaf structure and sclerophylly of Sebastiania myrtilloides (Euphorbiaceae) from Serra do Cipó, Minas Gerais, Brazil. Boletim de Botânica da Universidade de São Paulo 18: 21–27.
  • Midgley GF, Wand SJE, Musil CF (1998) Repeated exposure to enhanced UV-B radiation in successive generations increases developmental instability (leaf fluctuating asymmetry) in a desert annual. Plant, Cell & Environment 21: 437–442.
  • Moller AP, Shykoff JA (1999) Morphological developmental instability in plants: patterns and causes. International Journal of Plant Sciences 160: 135–146.
  • Niinemets Ü (2010) A review of light interception in plant stands from leaf to canopy in different plant functional types and in species with varying shade tolerance. Ecological Research 25: 693–714.
  • Norberg RA (1978) Skull asymmetry, ear structure and function, and auditory localization in Tengmalm’s owl, Aegolius funereus (Linné). Philosophical Transactions of the Royal Society B: Biological Sciences 282(991): 325–410.
  • Puerta‐Piñero C, Gómez J, Hódar J (2008) Shade and Herbivory Induce Fluctuating Asymmetry in a Mediterranean Oak. International Journal of Plant Sciences 169(5): 631–635.
  • Rozendaal DMA, Hurtado VH, Poorter L (2006) Plasticity in leaf traits of 38 tropical tree species in response to light; relationships with light demand and adult stature. Functional Ecology 20: 207–216.
  • Sander TB, Matthies D (2017) Fluctuating asymmetry of leaves is a poor indicator of environmental stress and genetic stress by inbreeding in Silene vulgaris. Ecological Indicators 79: 247–253.
  • Santos JC, Alves-Silva E, Cornelissen TG, Fernandes GW (2013) The effect of fluctuating asymmetry and leaf nutrients on gall abundance and survivorship. Basic and Applied Ecology 14: 489–495.
  • Sridevi V, Satyanarayana NV, Madhava Rao KV (1999) Induction of heat shock proteins and acquisition of thermotolerance in germinating pigeonpea seeds. Biologia Plantarum 42: 589–597.
  • Venâncio H, Alves-Silva E, Santos J (2016) Leaf phenotypic variation and developmental instability in relation to different light regimes. Acta Botanica Brasilica 30(2): 296–303.
  • Waddington CH (1956) Principles of Embryology, George Allen & Unwin, London.
  • Yang B, Yu Z, Srivastav A, Datir S, Kumar V (2019) Light Stress Responses and Prospects for Engineering Light Stress Tolerance in Crop Plants. Journal of Plant Growth Regulation: 1–18.
  • Wuytack T, Wuyts K, Van Dongen S, Baeten L, Kardel F, Verheyen K, Samson R (2011) The effect of air pollution and other environmental stressors on leaf fluctuating asymmetry and specific leaf area of Salix alba L. Environmental Pollution 159(10): 2405–2411.
  • Zhu H, Li X, Zhai W (2017) Effects of low light on photosynthetic properties, antioxidant enzyme activity, and anthocyanin accumulation in purple pak-choi (Brassica campestris ssp. Chinensis Makino). PLoS ONE 12(6): e0179305.

Supplementary material

Supplementary material 1 

Supplementary files

Bailey Francis, Robert Tucker Gilma

Data type: morphometric data

This dataset is made available under the Open Database License ( The Open Database License (ODbL) is a license agreement intended to allow users to freely share, modify, and use this Dataset while maintaining this same freedom for others, provided that the original source and author(s) are credited.
Download file (5.79 kb)
login to comment