|Home | About | Journals | Submit | Contact Us | Français|
Background Discussions of phenotypic robustness often consider scenarios where invariant phenotypes are optimal and assume that developmental mechanisms have evolved to buffer the phenotypes of specific traits against stochastic and environmental perturbations. However, plastic plant phenotypes that vary between environments or variable phenotypes that vary stochastically within an environment may also be advantageous in some scenarios.
Scope Here the conditions under which invariant, plastic and variable phenotypes of specific traits may confer a selective advantage in plants are examined. Drawing on work from microbes and multicellular organisms, the mechanisms that may give rise to each type of phenotype are discussed.
Conclusion In contrast to the view of robustness as being the ability of a genotype to produce a single, invariant phenotype, changes in a phenotype in response to the environment, or phenotypic variability within an environment, may also be delivered consistently (i.e. robustly). Thus, for some plant traits, mechanisms have probably evolved to produce plasticity or variability in a reliable manner.
Robustness is often defined as the ability of an organism to produce a constant phenotype in the face of environmental, genetic and stochastic perturbations (Debat and David, 2001; Lempe et al., 2013). Waddington coined a similar term, canalization, to describe how members of a population may have the same phenotype despite environmental and genetic variation (Waddington, 1942). Plastic traits, which vary with the environment, are sometimes described as uncanalized and may be considered to result from a lack of robustness (Wagner et al., 1997; Debat and David, 2001). Likewise, traits that are highly variable within an environment are often considered to be a result of poorly constrained development (Debat and David, 2001). However, as sessile organisms, sensitivity to the environment is an essential feature of plant development, and a number of well-characterized plastic responses allow plants to tune their phenotype reliably to prevailing environmental conditions. For example, in response to shading from neighbours, some species undergo increased stem elongation in a well-characterized shade avoidance response (Ballaré et al., 1987). Additionally, there are scenarios where increased phenotypic variability within a single environment may be advantageous. For instance, in environments with unpredictable weather patterns, there is evidence that increased variability in the timing of germination is adaptive in some species (Venable, 2007; Simons, 2009). Therefore, mechanisms that reproducibly (i.e. robustly) deliver increased phenotypic variance between or within environments may be under selection. Here we argue that it is useful to categorize phenotypes based on their variability within and between environments, rather than based on the narrower concept of robustness.
We use three terms to describe different types of phenotype: invariant, plastic and variable (Table 1). Throughout the text we use the term ‘phenotype’ to refer to characteristics of specific traits (e.g. stem length, number of branches and flower colour) rather than to describe the composite of a plant’s characteristics.
For a given genotype, an invariant phenotype is relatively constant and insensitive to differences in the environment (Table 1; Fig. 1). On the other hand, a plastic phenotype changes reproducibly in response to a change in the environment (Fig. 1). A variable phenotype varies substantially even between genetically identical individuals within the same environment (Fig. 1). This may be due to properties of regulatory networks underlying development of the trait which amplify stochasticity (see text below and Fig. 2). Within a given environment, plastic phenotypes may be more or less variable (Fig. 1).
For each of these phenotypic categories, we discuss the ecological scenarios under which it might be advantageous. In addition, drawing from a large body of work on unicellular organisms and from work on variability and plasticity in multicellular systems, we review the mechanistic principles that could underlie the development of each type of phenotype.
In addition to developmental mechanisms of plasticity and variability, the ability of a given genotype to produce multiple phenotypes may be attributed to heritable epigenetic changes which can be influenced by the environment and cause transgenerational changes in gene expression (Richards, 2006; Bossdorf et al., 2008). By mediating maternal effects on offspring gene expression, such transgenerational epigenetic mechanisms are likely to cause differences between families of individuals. An in-depth discussion of specifically epigenetic mechanisms underlying phenotypic variation is beyond the scope of this review. Here we focus on mechanisms that could generate phenotypic variability even between genetically identical individuals descended from the same parents.
Although plant development is remarkably plastic, some phenotypes are relatively constant, even between different environments. Early on in arabidopsis development, embryogenesis consistently generates mature embryos with the same basic body plan including correctly positioned shoot and root apices and two cotyledons (Park and Harada, 2008). This sets the foundation for growth and patterning at later stages of development which are often more environmentally sensitive. Flower patterning and morphology are also often relatively invariant. For example, the numbers of each type of floral organ (sepals, petals, stamens and carpels) are remarkably constant between individuals (Baker et al., 2005) and environments. The relative invariance of flower development may help to ensure that the different organs of the flower remain co-ordinated in shape and size throughout flower development, allowing protection of the reproductive organs and successful pollination. Invariance in flower morphology is also likely to be important from an ecological perspective. For example, within species that rely on a specific animal pollinator, individuals consistently generate specific flower morphologies that are recognized by specific pollinators, and, in some cases, restrict entry to the desired pollinator (Whitney and Glover, 2007). In this scenario, invariance in flower morphology and patterning is likely to be important for successful attraction, entry and pollination by insects. Invariant phenotypes may also be selected for in extreme and relatively invariant environments where a consistent and environmentally insensitive phenotype is likely to be favourable. For example, in contrast to shade-avoiding species that exhibit plasticity in stem length in response to shading, shade-tolerant species have relatively invariant stem lengths, regardless of light quality. A short, invariant stem length is hypothesized to reduce carbon wastage due to futile stem elongation and maintenance in scenarios where stem elongation will not enable light capture (Givnish, 1988; Valladares and Niinemets, 2008). Such invariant phenotypes are likely to have evolved under stabilizing selection that removes genotypes giving rise to deviations away from an optimum phenotype (Gibson and Wagner, 2000).
The development of invariant phenotypes is probably controlled by mechanisms that limit the effects of environmental variation and stochasticity on the final phenotype. As most of plant development is post-embryonic, developing organs are often exposed to environmental fluctuations much larger than those that normally occur during the development of many animal species (Lempe et al., 2013). Additionally, the molecular processes of transcription and translation are inherently noisy and can generate random differences in protein levels between cells (Raj and van Oudenaarden, 2008). Stochastic differences between cells can also be generated by spontaneous epi-mutations that alter the epigenetic regulation of gene expression (Blewitt et al., 2004). Noise may also arise independently of differences in gene expression during the allocation of cellular components to daughter cells during cell division. If components such as proteins and organelles are present at low copy numbers per cell, their stochastic allocation to daughters will lead to variability in a cell’s offspring (Berg, 1978; Oates, 2011). When all these processes are considered, the development of an invariant trait is clearly non-trivial. Such processes also have the potential to disrupt the generation of consistent phenotypic changes in response to specific environmental cues during the development of plastic phenotypes.
The extent of variability in development is under genetic control, and, in arabidopsis and maize, quantitative trait loci (QTL) controlling the extent of variability within individual environments have been mapped (Hall et al., 2007; Ordas et al., 2008; Jimenez-Gomez et al., 2011; Shen et al., 2012; Joseph et al., 2015). These studies have shown that some loci influence phenotypic variance without affecting the mean phenotype, and that distinct loci influence the variance of different phenotypes (Joseph et al., 2015). This suggests that specific molecular mechanisms exist to regulate the variability of particular traits.
Developmental studies addressing how phenotypes can be buffered against environmental and stochastic influences are beginning to emerge, and a large body of work on microbes has elucidated mechanisms that can promote cellular homogeneity by limiting the effects of noise in gene expression. Below, we discuss sources of noise in gene expression and how the architecture of gene regulatory networks, and the regulatory environment in which they operate, may buffer against stochasticity and environmental variability. Such mechanisms are likely to be important in the development of both invariant and plastic phenotypes. In some cases, the mechanisms discussed might also maintain invariant phenotypes in the face of genetic variation, although this is not covered here (Stearns et al., 1995; Rutherford and Lindquist, 1998; Gibson and Wagner, 2000; Meiklejohn and Hartl, 2002; Queitsch et al., 2002; Fraser and Schadt, 2010; Lempe et al., 2013).
Changes in gene expression that occur in response to spatial and temporal cues form the core of most developmental mechanisms. However, the small numbers of molecules that are involved in these processes make them inherently noisy (Li and Xie, 2011). Transcription initiation occurs stochastically due to a low copy number of genes per cell, and the random arrival of transcription factors at promoters (Blake et al., 2003; Raser and O’Shea, 2004; Newman et al., 2006). Following initiation of transcription, a burst of mRNA production may occur (Golding et al., 2005; Blake et al., 2006; Chubb et al., 2006; Raj et al., 2006; Raj and van Oudenaarden, 2008). Each mRNA molecule can then be translated multiple times, resulting in a burst of protein production following the stochastic initiation of transcription (Ozbudak et al., 2002; Fraser et al., 2004; Yu et al., 2006).
Stochasticity in protein production not only means that an individual cell will vary over time, but also that at a given time in development, cells within a population will be heterogeneous in their protein levels (Lander, 2011). This could disrupt developmental processes when groups of cells are required to respond to upstream signals. Additional variability can arise since, when a protein regulates the expression of downstream genes, fluctuations in its levels can propagate through the gene-regulatory network (Thattai and van Oudenaarden, 2001; Becskei et al., 2005; Hooshangi et al., 2005). Stochasticity may even be amplified in gene expression cascades where multiple noisy components are placed in series (Lander, 2011). Noise can also disrupt the cellular memory of previously received cues (Acar et al., 2005). Thus, stochastic fluctuations in gene expression have the potential to disrupt developmental signalling and increase phenotypic variability between individuals.
These stochastic fluctuations in protein levels can be minimized if rates of transcription are relatively high and the rates of translation per mRNA are low (e.g. due to elevated mRNA degradation) since this prevents bursts of proteins being generated from rarely produced mRNAs (Ozbudak et al., 2002; Fraser et al., 2004). Consistent with the relative rates of transcription and translation being regulated to limit stochastic effects, in yeast, essential genes tend to have higher rates of transcription and lower rates of translation than non-essential genes which produce proteins at similar rates (Fraser et al., 2004). Protein turnover rates may also be regulated to reduce stochastic fluctuations. If a protein decays relatively slowly, its total level in the cell will be governed by the integral of its synthesis over time, allowing the effects of stochastic and environmental fluctuations on production rates to be averaged out (Raj and van Oudenaarden, 2008).
The interactions between genes in a cell, or between gene-regulatory networks in different cells, may be structured to minimize variation in development despite stochastic or environmentally induced fluctuations in protein production.
Negative autoregulation is a common motif in genetic networks and involves a protein feeding back to downregulate its own expression, either directly or through the production of a second regulatory protein (Fig. 2A) (Alon, 2007). This feedback may help to restrict fluctuations in a protein’s levels and reduce noise propagation through a network: if the protein’s levels become elevated, its production will decrease; but, if its levels drop, its production rate will increase (Becskei and Serrano, 2000; Thattai and van Oudenaarden, 2001; Nevozhay et al., 2009). However, due to the noise inherent in gene expression, the ability to suppress noise through negative feedback involving regulatory proteins is likely to be limited to cases where large numbers of regulatory molecules are involved (Lestas et al., 2010). This is because, in a system where protein X promotes the production of a regulatory protein (Y) which in turn downregulates X (Fig. 2A), the production of Y in response to X will be a stochastic event. Unless large numbers of Y molecules are made in response to X (minimizing stochastic effects), the total levels of Y are unlikely accurately to reflect the amount of X, leading to imperfect control of X by Y (Lestas et al., 2010). Also, if there is a time delay in the production of Y and its effect on X's expression, this could lead to oscillations in the levels of both proteins (Novák and Tyson, 2008).
Many flowering plants have large genomes as a result of both small- and large-scale genome duplications, with two whole-genome duplications believed to have occurred during the evolution of angiosperms (Cui et al., 2006; Vanneste et al., 2014). As a consequence, many genes have multiple paralogues, which, in some cases, have at least partly overlapping functions (Veitia, 2005). Such redundancy can make networks resistant to mutation. For example, a number of families of transcription factors that each have multiple redundant members contribute to the development of distinct upper (adaxial) and lower (abaxial) surfaces of leaves and other plant lateral organs (McConnell and Barton, 1998; Eshed et al., 2001; McConnell et al., 2001; Emery et al., 2003; Yamaguchi et al., 2012). Adaxial and abaxial identity determinants mutually inhibit each other to specify precisely an adaxial–abaxial boundary along the mid-plane of the leaf. Loss-of-function mutations in single components rarely cause a complete loss of patterning, indicating that the network is robust to changes in the levels of its components (Eshed et al., 2001; Emery et al., 2003; Husbands et al., 2009). The presence of redundant copies of the transcription factors involved could also reduce variability due to stochastic effects. If multiple transcription factors perform a common function, stochastic bursts in the production of each individual protein would be buffered at the level of downstream gene expression. This is because this output would be proportional to the total level of all the transcription factors at a given time, allowing the effects of stochastic bursts in the production of each one to be averaged out across the different proteins.
The development and correct behaviour of multicellular structures requires that the gene regulatory networks operating in individual cells can be co-ordinated and regulated to achieve a reliable tissue-scale output. This can be achieved through intercellular signalling which links the behaviour of the gene-regulatory networks operating in each cell to its environment within the tissue.
Intercellular signalling may allow the dynamics of intracellular processes to be synchronized across a tissue. This occurs in the mammalian suprachiasmatic nucleus (part of the hypothalamus) where neurons show circadian oscillations in firing rates. If these cells are dissociated, their circadian rhythms become desynchronized (Welsh et al., 1995). Somewhat similarly, in arabidopsis leaves, the presence of spatio-temporal waves of circadian clock gene expression suggests that intercellular coupling helps to co-ordinate the circadian clocks of individual cells (Wenden et al., 2012). A similar phenomenon of intracellular coupling was observed in the context of somite development in the vertebrate mesoderm. The segmentation of the mesoderm into somites is controlled by oscillations in gene expression that provide temporal cues for somite formation, referred to as the segmentation clock (Palmeirim et al., 1997). If mesodermal cells are dissociated from each other, the oscillations within individual cells become unstable, with highly variable amplitudes (Masamizu et al., 2006). These examples suggest that signalling between cells synchronizes and stabilizes their oscillations, which would otherwise be noisy.
In contrast to the generation of tissue-scale oscillations, the maintenance of a relatively constant average level of gene expression within a tissue may sometimes be best achieved by regulatory networks that prevent synchronization of groups of cells. This has been proposed to be the case for nuclear factor-κB (NF-κB) signalling, which plays an important role in inflammation (Paszek et al., 2010). Upon perception of cytokines (paracrine signalling factors), the NF-κB transcription factor oscillates between being located in the cytoplasm and in the nucleus (Nelson et al., 2004) where it influences the expression of downstream genes, including those encoding cytokines (Tian et al., 2005; Ashall et al., 2009; Paszek et al., 2010). If neighbouring cells were co-ordinated in their perception and production of cytokines, this system would provide the potential for positive feedback (between the production of cytokines and NF-κB expression), leading to a strong inflammatory response. However, through computational modelling, it was shown that the signalling network involved appears to be optimal for giving rise to NF-κB oscillations that are unco-ordinated between cells (Paszek et al., 2010). It has been proposed that by desynchronizing the oscillations of individual cells, the gene-regulatory network restricts the extent of positive feedback and prevents large fluctuations in amounts of cytokines at the tissue level (Paszek et al., 2010).
Intercellular signalling via morphogen gradients plays an important role in large-scale developmental patterning by generating reproducible differences in gene expression across a tissue (Lander, 2011). The classical model of a morphogen gradient involves a graded signal inducing separate domains of gene expression in a field of cells that respond to threshold concentrations of the signal (Wolpert, 1969). To generate boundaries of gene expression that are reproducible between individuals, neighbouring cells must consistently be able to perceive differences in morphogen concentration and induce different gene expression states (Gregor et al., 2007). Therefore, stochastic fluctuations in morphogen concentrations and in the ability of cells to perceive the morphogen are likely to impact on the precision of this process. However, a number of different strategies appear to be employed in plant and animal development to generate invariant outcomes from graded intercellular signals.
Early in Drosophila embryogenesis, a maternally controlled anterior–posterior gradient in Bicoid (Bcd) protein activates zygotic expression of hunchback (hb) in a concentration-dependent manner (Driever and Nusslein-Volhard, 1988a, b; Struhl et al., 1989). hb expression is activated in the anterior half of the embryo, where the Bcd concentration is highest. The relationship between Bcd and Hb protein concentrations is remarkably precise despite a relatively high degree of stochasticity in the transcription of hb (Gregor et al., 2007; Little et al., 2013). This has been proposed to be due to spatial and temporal averaging of Bcd concentrations during the control of hb expression (Gregor et al., 2007; Little et al., 2013). hb transcripts have a relatively long lifetime and therefore Hb protein levels are a read-out of the cumulative production of hb mRNAs over time. As described above, time averaging of hb transcription in this way has the potential to reduce the effects of stochastic transcription initiation and noise in Bcd concentration (Little et al., 2013). Also, as the embryo is a syncitium at this stage of development, averaging may occur through short-range diffusion of Hb between neighbouring nuclei, which can smooth out fluctuations in Bcd and Hb levels and yield more precise boundaries (Erdmann et al., 2009). Boundaries specified by the Hb morphogen gradient appear then to be further refined by cross-regulation between gap genes induced in neighbouring tissue domains (Jaeger, 2011; Gursky et al., 2012).
Interestingly, plants appear to employ a unique strategy to overcome the disruptive effects of noise on the sensing of gradients: they maximize the steepness of the gradient over the whole region of the tissue in which patterning occurs. This allows them to solve the problem of generating large-scale positional information in the face of stochastic perturbations, since the difference in signal between adjacent cells will be larger than the noise. This appears to be the case in the arabidopsis root, where an auxin concentration gradient, with its maximum at the root tip, influences specification of the boundaries between different zones of cellular behaviour (Aida et al., 2004; Blilou et al., 2005; Galinha et al., 2007; Grieneisen et al., 2007; Petersson et al., 2009). Rather than being generated through a localized source combined with global decay (as is the case for many morphogen gradients in animals), the auxin gradient is generated through a reflux loop of auxin transport created by the pattern of polarly localized PIN proteins (auxin efflux carriers) in the root (Blilou et al., 2005; Grieneisen et al., 2007). In the internal tissue of the root, auxin is transported towards its tip, while, in outer tissue, transport is directed back up towards the shoot. In intermediate cell layers, PIN proteins are localized inwards, returning auxin back to inner tissue from the outer cell layers (Blilou et al., 2005). Taking the known biophysical parameters involved in auxin transport into account, the shape and steepness of the auxin gradient generated through the reflux loop mechanism was compared with other possible mechanisms, such as source decay or unidirectional transport towards the root tip (Grieneisen et al., 2012). This revealed that the auxin transport reflux loop is probably the most effective mechanism for generating the observed steep auxin gradient over the large spatial range of the root.
Another approach to achieve invariant tissue patterning is to use two or more signalling systems in parallel to provide positional information. Indeed, recent work on phyllotaxis revealed that plants achieve relatively reproducible spacing of successive organs through the action of two interacting signalling systems, one involving polarized transport of auxin, and the other involving diffusion of a cytokinin signalling inhibitor (Besnard et al., 2014). The auxin efflux carrier, PIN1, forms polarity convergence points with high intracellular auxin that specify the position of each new organ in the shoot apical meristem (Reinhardt et al., 2000, 2003). The polarized transport of auxin towards existing primordia is predicted to cause depletion of auxin in surrounding tissue, meaning that new primordia can only be initiated at a distance from existing ones where the auxin concentration becomes sufficiently high (Jönsson et al., 2006; Smith et al., 2006; Stoma et al., 2008). However, in the presence of variability in auxin concentrations and primordium size and upon increases in meristem size, the positional information generated by this mechanism may be ambiguous, causing the next two sites for organ initiation to be indistinguishable (Mirabet et al., 2012; Landrein et al., 2014). This can cause pairs of organs to initiate simultaneously or the order of organ initiation to be swapped. Mutations in AHP6, a cytokinin signalling inhibitor, cause an increased frequency of these aberrations, suggesting that this second signal contributes to the invariance of phyllotaxis (Besnard et al., 2014). Indeed, accumulation of auxin in new primordia activates AHP6 expression, and intercellular movement of AHP6 generates a graded distribution around each young primordium. This additional layer of information can be used to distinguish between the next two sites for primordium initiation, decreasing the variability of phyllotactic patterning.
In the context of growth, intercellular signalling may ensure that a precise final state is reached despite stochasticity throughout development. For example, a signal that terminates growth may be increasingly produced in a tissue as it approaches its final state, ensuring an accurate final outcome. Such a mechanism functions in the mouse olfactory epithelium, where stem cells divide to produce progenitors that pass through distinct lineage stages before differentiating into mature neurons (Lander et al., 2009; Lander, 2011). Signalling molecules are produced by differentiated cells that feed-back to inhibit the self-renewal of progenitor cells, terminating proliferation at the correct time. A model of such a system was shown to produce a reproducible final outcome even in the face of variation in cell cycle speeds and initial stem cell number (Lander et al., 2009). A role for non-cell-autonomous signals in regulating the final size of determinate plant organs is indicated by the phenomenon of compensation, whereby mutants with reduced cell division still reach a similar final organ size by adjusting the rates or duration of organ expansion (Ferjani et al., 2007; Powell and Lenhard, 2012). The mechanism underlying compensation could contribute to the consistent generation of similar organ sizes in the face of stochastic variations in individual cell behaviour.
A role for tissue stresses in regulating growth heterogeneity has also been proposed. Heterogeneous growth of neighbouring regions of tissues is predicted to create local stresses. If these were sensed and used to influence cellular behaviour, this could in turn influence growth variability (Shraiman, 2005). Such feedback appears to occur in plants where cells orient their microtubules parallel to the local principle direction of stress (Hamant et al., 2008). As a consequence, cell walls become reinforced along the principle direction of stress as microtubules guide the synthesis of cellulose microfibrils (Lloyd and Chan, 2008). Thus, growth in this direction is expected to be restricted (Hamant et al., 2008). Via computational modelling, it was shown that when cell-autonomous growth rates (the rate at which a cell would grow if removed from the tissue) are variable, such mechanical feedback can reduce overall growth rate variability in the tissue (Uyttewaal et al., 2012). However, if the strength of feedback is increased, it can effectively overshoot and amplify differences in growth rate between regions. Surprisingly, this latter scenario seems to be the case in the shoot apical meristem, where an apparent reduction in mechanical feedback (in a mutant with a reduced ability to align its microtubules with the principle direction of stress) causes reduced growth rate variability (Uyttewaal et al., 2012). In the meristem, amplified growth variability due to mechanical feedback was proposed to promote organ emergence (Uyttewaal et al., 2012). However, little is known about the effects of mechanical feedback on growth heterogeneity in other tissues.
The examples discussed above illustrate how changes in the properties of specific gene-regulatory networks can influence the variability of their output on a multicellular scale. However, there are other more general processes that may reduce variability across a large number of different gene-regulatory networks in parallel.
Chaperone proteins bind to a number of target proteins and influence their stability. The HSP90 chaperone binds to developmental regulators and other signalling proteins when they are in unstable states. When bound, the chaperone prevents its target protein from misfolding until a correct signalling event occurs and triggers a conformational change (Zuehlke and Johnson, 2010). Disruption of HSP90 function in Drosophila or arabidopsis leads to increased variability in many phenotypes, suggesting that its role in protein stability contributes to invariance in a number of developmental processes (Rutherford and Lindquist, 1998; Queitsch et al., 2002).
Recent work on yeast suggests that, within highly connected networks, there are many other genes that play a similar role to HSP90 in reducing variability. When thousands of haploid yeast single-gene knockout strains were screened for their phenotypic variance, about 5 % of genes were required for achievement of a wild-type level of phenotypic variability, causing increased variability when knocked out (Levy and Siegal, 2008). A large proportion of these genes was involved in core cellular processes such as chromosome organization, mRNA processing and protein modification. Also, their products tended to be highly connected within protein–protein interaction networks, suggesting that they may have general functions with many client proteins, as is the case with HSP90.
In plant tissues, neighbouring cells are connected via cell walls and therefore cannot slide relative to each other. This connectedness of tissues constrains the growth of each local region so that each cell or local region of the tissue may not grow at the same rate as it would if it were in isolation (Coen et al., 2004). The constraints imposed by such mechanical connectivity can restrict the overall shape changes occurring in a tissue as a result of locally altered growth (Boudon et al., 2015) and therefore could result in locally variable growth being smoothed out on the tissue scale, increasing homogeneity. In contrast to the mechanism of mechanical feedback discussed above, this could occur without mechanical forces being sensed by cells to alter directly specific aspects of their behaviour. This mechanism would therefore act as a more general buffer for growth variability.
An important factor influencing whether plant phenotypes can be buffered against environmental variability appears to be the effect of stress on growth rate (Achard et al., 2006). The perception of multiple types of stress can lead to a generally reduced growth rate throughout the plant via increased levels of DELLA proteins, which repress the activity of growth-promoting transcription factors (Achard et al., 2006; Harberd et al., 2009). This DELLA-dependent mechanism of stress-induced growth reduction was shown to increase survival in stressful environments (Achard et al., 2006). There are several ways in which slow growth might confer stress resistance. One possibility is that by slowing down growth in adverse environmental conditions, resource consumption by growth is temporarily reduced, allowing resources to be diverted towards homoeostatic mechanisms that counteract the effect of stress (Harberd et al., 2009). Alternatively, slow growth may be inherently less sensitive to stress and stochasticity because of the buffering provided by integrating regulatory processes over longer periods of time. Slow growth could therefore contribute to the development of invariant phenotypes by reducing the phenotypic effects of potentially detrimental environmental fluctuations.
In summary, there are many mechanisms that might contribute to the development of invariant phenotypes by generating reproducible developmental outputs despite stochastic effects. Since environmental changes probably cause changes in gene expression, developmental mechanisms that confer robustness to stochastic changes in gene expression may also mitigate the effects of environmentally induced changes (Meiklejohn and Hartl, 2002). A number of these mechanisms relate to properties of specific gene-regulatory networks controlling developmental processes. Such networks may include negative feedback to reduce stochastic bursts in protein levels. The presence of paralogous genes with similar functions may also reduce stochasticity at the level of particular protein functions. At the tissue level, the gene-regulatory networks of individual cells may be coupled via intercellular signalling or maintained in a decoupled state to reduce variability across the tissue. During patterning, noise in developmental signals and their perception may be compensated for by employment of spatio-temporal averaging. The generation of steep signal gradients and the use of parallel patterning systems may also help to provide unambiguous positional information. Additionally, the regulatory context may act to reduce variability by providing a more generalized buffering of a number of different developmental processes. Chaperone proteins, physical constraints of growing tissues and the regulated repression of growth to aid the maintenance of homeostasis may each restrict the phenotypic variability arising from a number of different developmental pathways. These aspects of the regulatory environment might be seen as imposing a 30mph speed limit on a number of roads in parallel to prevent accidents.
The development of many plant phenotypes changes in response to the environment. For some traits, plasticity may be an adaptive response and involve signalling pathways that were selected to cause a specific change in development in response to the environment. However, not all plastic responses are adaptive, and some may be due to the inevitable consequences of resource limitation or stress on growth and physiology (Sultan, 2000; van Kleunen and Fischer, 2005). Many instances of phenotypic plasticity appear to be functionally appropriate and therefore examples of the adaptive scenario (Sultan, 2000). Here we illustrate types of molecular mechanisms that can underlie adaptive plasticity using two well-characterized examples in plants: shade avoidance and the control of flowering time by the perception of winter. We then review some empirical evidence for the adaptive significance of plastic phenotypes. Finally, we explore the relationship between plasticity and variability and discuss the benefits, limitations, and costs involved in the evolution of phenotypic plasticity, illustrating that under certain conditions, invariant or variable phenotypes may be favoured over plastic ones.
Plants often change their development in response to cues that are likely to vary spatially in their natural environment and impact upon fitness. One well-studied example is the avoidance of shading. When plants are present at high densities, competition for light can reduce their growth and fitness (Schmitt and Wulff, 1993; Schmitt et al., 1999). Plants are able to use the red to far-red light ratio to predict future shading, as nearby neighbours (which do not yet cause shading) absorb red light but reflect far-red light (Ballaré et al., 1987). In response, early flowering and a range of morphological changes can be triggered that increase chances of light capture, including increased elongation of stems and reduced branching (Schmitt and Wulff, 1993; Ruberti et al., 2012).
Plants employ a number of different photoreceptors to perceive their light environment and control development accordingly (Fig. 2B). The phytochrome B clade of photoreceptors is the most important in perception of the red to far-red light ratio, which is indicative of the degree of shading (Casal, 2013). Active phytochrome B physically interacts with the PHYTOCHROME INTERACTING FACTOR (PIF) family of transcription factors, causing either their degradation or inactivation. Increased far-red light promotes the conversion of phytochrome B to an inactive form, which allows PIF levels to increase rapidly and enables their binding to target promoters. PIFs control the expression of a number of genes required for shade avoidance, including cell wall-modifying enzymes and auxin biosynthetic enzymes. Indeed, many of the phenotypes associated with shade avoidance, including reduced branching, stem elongation and upwards bending of leaves, are characteristic of elevated auxin levels. Phytochrome B also feeds more directly into the pathways influencing flowering and branching (Jang et al., 2008; Finlayson et al., 2010). Thus, in the case of shade avoidance, phenotypic plasticity is caused by a specific mechanism of environmental perception that feeds into core developmental pathways.
As well as changing development in response to cues that may be spatially heterogeneous, plants use environmental cues to control the timing of key life history events such as flowering time and germination. Flowering time is influenced by a number of seasonal cues including daylength, temperature and prolonged periods of cold, as well as responding to indicators of unfavourable conditions such as salinity, shading and a lack of nutrient availability (de Jong and Leyser, 2012; Kim and Sung, 2014). Following their perception, these cues converge to influence the expression of three floral integrator genes. These integrators control flowering time by regulating expression of floral meristem identity genes, which when expressed result in a switch from vegetative to reproductive development (de Jong and Leyser, 2012; Kim and Sung, 2014).
One of the best characterized pathways influencing flowering time is the vernalization pathway, which promotes flowering in response to prolonged periods of cold (Ietswaart et al., 2012; Song et al., 2012; Kim and Sung, 2014). Arabidopsis accessions differ in their requirement for vernalization, and this largely depends on genetic variation at the FRIGIDA (FRI) locus. FRI promotes the expression of a floral repressor, FLOWERING LOCUS C (FLC), which becomes stably downregulated in the cold. Accessions from northern Europe tend to have functional alleles of FRI, and therefore require vernalization to flower. These winter annual accessions tend to over-winter in a vegetative state (when vernalization occurs) and then flower the following year in response to increasing daylength. For this to occur, FLC downregulation must be mitotically stable for several months, providing a memory of perceived winter.
Recent work has shown that FLC expression is modulated by multiple regulators to ensure stable and robust downregulation following winter (Ietswaart et al., 2012; Song et al., 2012). After around 2 weeks of cold treatment, FLC undergoes a decrease in expression levels. This appears to involve a number of overlapping pathways that recruit Polycomb repressive complex 2 (PRC2), a highly conserved chromatin-modifying complex, to the FLC locus (Heo and Sung, 2011; Song et al., 2012). Through the action of PRC2, continued cold causes a gradual build-up of repressive histone modifications at a particular region (termed the nucleation region) of the FLC locus. Upon return to warm conditions, repressive modifications can spread throughout the rest of the locus, and this spreading is required for the maintenance of FLC repression and memory of vernalization (Angel et al., 2011).
The stable silencing of FLC upon return to warm appears to be a probabilistic event in each cell: following a moderate duration of cold, only a small proportion of nuclei stably silence FLC, whereas a longer duration of cold allows a higher proportion of nuclei to establish stable silencing (Angel et al., 2011). The floral integrator FLOWERING-LOCUS T (FT ) is repressed by FLC and can act as a long-range signal. FT levels are therefore thought to provide an average read-out of FLC expression across a number of cells (Corbesier et al., 2007; Angel et al., 2011). Thus, the extent of flowering acceleration is proportional to the time spent in the cold and the proportion of nuclei that stably repress FLC expression. The slow and quantitative response of the FLC locus to cold means that several weeks of cold are needed for flowering to be significantly accelerated. This is likely to remove sensitivity to fluctuations in temperature occurring at short time scales (e.g. a one-off frost in autumn) but allows a robust plastic response upon perception of meaningful environmental cues.
As mentioned above, not all instances of plasticity will be adaptive, and, in fact, evidence for plasticity conferring a fitness advantage is relatively difficult to acquire. This is because it is necessary to show that the plastic phenotype confers an increased fitness across relevant environments relative to a non-plastic phenotype (Schmitt et al., 1999; Sultan, 2000). One way to do this is to compare the fitness of a wild-type plant, which shows a plastic response, with that of mutants or transgenic plants that lack the plasticity. This was successfully achieved for the shade avoidance response (Schmitt et al., 1995, 1999). Wild-type tobacco plants were sown at high and low density alongside transgenic tobacco plants overexpressing PHYTOCHROME A (PHYA) which are defective in the shade avoidance response. At low density, PHYA overexpressor and wild-type plants had comparable fitness, indicating that PHYA overexpression does not cause deleterious pleiotropic effects. However, at high density, wild-type plants outcompeted the PHYA overexpressors, indicating that shade avoidance at high density confers an advantage. In a complementary experiment, wild-type plants of Brassica rapa were grown alongside mutants deficient in phytochrome B, which have a constitutive shade avoidance response. These mutants were less fit than the wild type when grown at low densities. Thus, the plasticity of shade avoidance-related phenotypes in wild-type plants confers an advantage at both low and high densities.
In the case of plasticity in flowering time, the requirement for vernalization in winter annuals from northern Europe is likely to be adaptive. Evidence for this comes from a study of 360 arabidopsis accessions that were sown in a field site in Rhode Island (northern USA) in autumn or spring (Korves et al., 2007). When sown in autumn, accessions carrying a functional allele of FRI (conferring a vernalization requirement) had better survival over winter (due to delayed flowering) and therefore increased fitness compared with accessions with a non-functional FRI allele (which flowered too early and thus had reduced fitness). Therefore, for plants native to northern latitudes, which will experience cold winters, acceleration of flowering following a period of cold is likely to be adaptive. However, in this case, it is likely to be environmental sensitivity rather than plasticity per se which is adaptive. Any one plant germinating close to its parent plant is unlikely to be required to respond appropriately to diverse lengths and severities of winter. Thus the delay in flowering of vernalization-requiring accessions in the absence of winter may not be advantageous (Korves et al., 2007). Indeed, there is evidence that non-functional alleles of FRI have undergone multiple selective sweeps in central Europe, indicating that a lack of vernalization requirement (and thus a lack of plasticity) might be advantageous in the absence of very cold winters (Le Corre et al., 2002; Le Corre, 2005).
Plastic phenotypes have sometimes been described as uncanalized against environmental variation, implying that plasticity results as a consequence of a lack of mechanisms conferring reliable developmental outcomes (Debat and David, 2001; Meiklejohn and Hartl, 2002). However, the shade avoidance and vernalization examples show that specific mechanisms have evolved that generate consistent and pre-emptive responses to the environment. In some cases, plasticity may therefore be considered to be robust.
However, sometimes plasticity of a phenotype may be correlated with the variability of this phenotype within an environment. A positive correlation was found within a set of arabidopsis recombinant inbred lines, where lines that showed more plastic responses to nitrate availability were also more variable within an environment (Tonsor et al., 2013). Consistent with these findings, in yeast, genes with highly plastic expression often show highly variable expression within a constant environment (Landry et al., 2007; Choi and Kim, 2009; Lehner, 2010). Genes with high expression variability both within and between environments tend to have specific mechanisms of gene expression regulation including the presence of a nucleosome in a specific regulatory region of the promoter, higher rates of histone exchange and the presence of TATA box elements in promoters (Choi and Kim, 2009; Lehner, 2010). However, at the gene expression level, plasticity and variability are not always correlated as essential genes tend to lack these mechanisms of gene expression regulation, but can still show plasticity whilst maintaining low variability (Lehner, 2010). Additionally, studies on Drosophila have not found a clear relationship between phenotypic plasticity and variability (Scheiner et al., 1991). In a study on seed germination in populations of Lobelia inflata, it was found that plasticity of germination timing in response to temperature was negatively correlated with the variability in germination time (Simons, 2014).
There are several mechanistic explanations for why plasticity and variability might sometimes be positively correlated. It may be the case that the development of some phenotypes is poorly constrained due to a lack of mechanisms buffering against environmental or stochastic perturbations. These phenotypes may therefore be variable within an environment and change in response to environmental fluctuations, showing plastic responses that are not necessarily adaptive. In other cases, the evolution of adaptive plasticity might lead to increases in variability. This could occur if the evolution of plasticity involves changes in regulatory architecture that allow an increased range of phenotypes to be produced in response to the environment. As a consequence of these changes, the phenotype may be less buffered against stochastic effects, leading to increased variability. Another possibility is that the evolution of sensitivity to meaningful environmental cues allowing a plastic response may, as a side effect, lead to increased sensitivity to microenvironmental cues, giving rise to increased variability within an apparently homogeneous environment. Alternatively, in some cases, separable mechanisms may give rise to increased plasticity and variability, and both may be favoured by selection in particular environments (see the section on variable phenotypes). In the case of the negative correlation between plasticity and variability seen in L. inflata seed germination, it has been proposed that this may be due to a constraint on the total amount of phenotypic variability for a specific trait that can be produced through the combination of plasticity or variability (Simons, 2014). To test these different possibilities, it is necessary to improve our understanding of the developmental mechanisms involved in phenotypic plasticity and variability and how these mechanisms change during the evolution of increased plasticity.
If plants could instantly change their development to suit best the prevailing environment this would appear to be an ideal strategy. Indeed, theoretical work shows that phenotypic plasticity increases fitness in spatially and temporally varying environments where environmental cues are good predictors of the optimal strategies in the future (Van Tienderen, 1997; Tufto, 2000; Givnish, 2002). Plasticity may also allow successful colonization of new or changed environments by allowing the persistence of a genotype in the altered conditions (Crispo, 2008). At early stages of colonization of a new environment, an increase in plasticity may be selected for and the plastic phenotype may become fixed (Ghalambor et al., 2007; Lande, 2009). If the plastic phenotype is sub-optimal, it may then be refined by natural selection operating in the new habitat (Pigliucci et al., 2006; Crispo, 2008; Lande, 2009). Consistent with a role for plasticity in successful colonization of new habitats, more invasive species tend to occupy a greater range of habitats, and have relatively high phenotypic plasticity (Baker, 1974; Williams et al., 1995; Sultan, 2003). However, despite the advantages of plasticity, many factors can constrain its evolution and influence the likelihood that a plastic phenotype is favoured over invariant or variable phenotypes in a given environment (DeWitt et al., 1998; Alpert and Simms, 2002; Givnish, 2002).
An inherent constraint to the evolution of plasticity is that genes conferring phenotypic plasticity, which are usually expressed only in certain environmental conditions, are likely to be under relaxed selection compared with genes expressed in all environments. This is because only the copies of the gene that are expressed will experience selection, whereas, if the inductive environment is never encountered, the gene will not be expressed, and its effect will be neutral. This means that advantageous alleles conferring a plastic phenotype will take longer to be fixed than those involved in the generation of an adaptive invariant phenotype (Snell-Rood et al., 2010). Thus, in environments where spatial or temporal variability occurs rarely (causing plasticity genes to be rarely expressed), non-plastic specialists may be favoured over generalists with relatively high plasticity in a number of traits.
Plasticity of a trait may be considered to have a cost if it causes a reduced mean fitness in a given environment compared with an invariant phenotype with the same mean trait value (DeWitt et al., 1998). Costs of plasticity could be related to detrimental pleiotropic or epistatic effects when a developmental process comes under environmental control. Because of this, increased modularity of developmental processes has been proposed to enable evolution of plasticity (Snell-Rood et al., 2010). If the development of different traits is controlled by relatively independent genetic networks, then a developmental module could evolve to become environmentally sensitive without disrupting other developmental processes. It has also been proposed that energetic costs associated with machinery that allows perception, integration and response to environmental cues could provide a cost to plasticity (DeWitt et al., 1998; Givnish, 2002), although it has been argued that these costs would probably be minimized by natural selection (Schlichting and Pigliucci, 1998; Givnish, 2002).
Another potential cost of plasticity is that the change in phenotype produced in response to a specific environmental cue might be detrimental in the presence of other co-occurring abiotic and biotic variables (Valladares et al., 2007). For example, tree seedlings that exhibit a strong stem elongation response to shading are more susceptible to cold temperatures than seedlings with a less plastic stem elongation phenotype (Valladares et al., 2007).
A limitation to the benefits of plasticity arises if, due to developmental constraints, a plastic phenotype cannot achieve as extreme a trait value in a given environment as an invariant phenotype can (DeWitt et al., 1998; Alpert and Simms, 2002). In severe environmental conditions where extreme trait values confer a selective advantage, invariant phenotypes may therefore be more advantageous. This is consistent with the observation that generalist species that inhabit a wide variety of relatively moderate environment conditions frequently exhibit developmental plasticity, while specialists that inhabit more extreme environments tend to have invariant phenotypes (Wilson and Yoshimura, 1994; Sultan, 2000, 2003).
Another limit to the effectiveness of plasticity comes from how well temporal environmental variability can be matched with changes in development. The lag time between the perception of an environmental change and a phenotypic outcome may be large compared with the time scale on which the environment changes. In this case, by the time a plastic phenotype is produced, it may be too late to prevent detrimental effects on fitness or it may no longer be appropriate for the current environment (DeWitt et al., 1998; Alpert and Simms, 2002). One solution to this problem is to use cues to influence development that predict future environmental changes (Alpert and Simms, 2002). The two instances of plasticity focused on here provide examples of such a solution. In the case of shade avoidance, morphological changes can be triggered before shading occurs through changes in the red:far-red light ratio due to reflection from neighbours (Ballaré et al., 1987). Flowering time is also influenced by a number of cues (such as cold) that indirectly indicate future optimal conditions for flowering (de Jong and Leyser, 2012).
However, in environments that vary stochastically, predictable cues about the environment to be encountered might not exist and plasticity is unlikely to be advantageous (Leimar et al., 2006). Indeed, using stochastic simulations of populations of individuals in variable environments, Reed et al. (2010) showed that strong plasticity is favoured when environmental cues provide good predictors of future selective pressures, but increases the risk of extinction when cues are unreliable. Low predictability of the optimum strategy for the future environment may have especially important implications for irreversible developmental events duch as the transition to flowering or seed germination (Cohen, 1966; Simons and Johnston, 2003). In the absence of reliable cues, variable phenotypes, where an individual genotype produces a range of different phenotypes relatively independently of environmental cues, have been proposed to be favourable over phenotypic plasticity (Simons, 2011).
Variable phenotypes may be produced as bet-hedging strategies in unpredictable environments where fluctuating natural selection favours different phenotypes across generations (Cohen, 1966; Lewontin and Cohen, 1969; Simons, 2011). As discussed above, in such an environment, the advantages of phenotypic plasticity may be limited due to the lack of reliable cues that can be used to provide information about the optimal future phenotype. An invariant phenotype may also be disadvantageous. This is because, even if it was optimal for the mean environment across generations, an invariant phenotype could result in extinction if extreme variations in the environment occurred. In contrast, a variable phenotype, produced as a bet-hedging strategy (see Table 1), would allow continued survival of a fraction of the population, even though the increased phenotypic variability would probably reduce population fitness in the mean environment. Thus, in unpredictable environments, variance in a phenotype may confer a selective advantage by reducing variance in fitness from generation to generation. This reduction in fitness variance across generations can have a large positive effect on long-term population fitness, which is severely reduced by occasional drops in population size (Gillespie, 1974; Seger and Brockmann, 1987; Philippi and Seger, 1989; Simons, 2011).
As discussed previously, phenotypic variability is under genetic control, and molecular mechanisms exist that can modulate its levels (Queitsch et al., 2002; Hall et al., 2007; Ansel et al., 2008; Viney and Reece, 2013). Thus, it is possible that increased variability as a bet-hedging strategy could be under positive selection in unpredictable environments. Indeed, there is empirical evidence that this has occurred in a number of systems (Simons, 2011). One of the best studied examples is the phenomenon of bacterial persistence, whereby in stressful conditions (such as in the presence of antibiotics) a slow-growing fraction of cells in a clonal population can survive and therefore later regenerate the population. Even in the absence of stress, bacteria can switch stochastically into a slow-growing mode, which can persist if conditions become stressful (Balaban et al., 2004; Maisonneuve et al., 2013). The slow-growing phenotype decreases the mean fitness of the genotype in favourable conditions, but reduces the risk of extinction should stressful conditions arise (Martins and Locke, 2015).
In plants, seedling mortality is often very high, and ecological theory suggests that bet-hedging germination strategies are likely to be beneficial in temporally variable and unpredictable environments (e.g. deserts with sporadic rainfall patterns) (Cohen et al., 1966; Simons and Johnston, 1997). Thus, species in unpredictable environments are expected to have increased variances in germination time. This hypothesis has been tested in desert annuals, which spread germination between seasons by producing a fraction of seeds that is dormant and resists germination in the season of seed production. Consistent with bet-hedging theory, a study of ten desert annual species over 22 years revealed that species experiencing more unpredictable reproductive success upon germination had higher fractions of dormant seeds (Venable, 2007). The species with increased variability in germination time tended to be those most sensitive to variation in precipitation and thus with greater variability in survival.
Spreading offspring germination times within a season to cope with an unpredictable environment can also constitute a bet-hedging strategy. This phenomenon has been studied in L. inflata, an annual species from a temperate habitat in which most seeds are non-dormant and show relatively high variation in time to germination within a season (Simons and Johnston, 2006). Naturally inbred lines of L .inflata show high variance in germination time even when grown under controlled laboratory conditions (Simons and Johnston, 2006). Simons et al. (2009) estimated the extent to which selection on L. inflata germination time fluctuates between years by germinating seeds in the lab and transferring the seedlings to field sites at different times of the year, for 5 years. The survival of plants transferred at different times of the year was measured so that for each year the optimum timing of germination could be estimated. This revealed that the optimal germination time varied widely between the 5 years. Simulations to determine the optimum germination strategy in these conditions revealed that variability in germination time, rather than synchronous germination at the mean optimum time across years, maximized long-term fitness (Simons, 2009). The extent of variability predicted to be optimal was in fact similar to the observed variability, suggesting that the observed germination strategy is consistent with bet-hedging (Simons and Johnston, 2006; Simons, 2009).
Selection for variable germination rates may also minimize intersibling competition for resources. Selection for increased variance in germination rates is expected to be particularly strong for selfing species that produce large numbers of seeds simultaneously and have short seed dispersal distances (Cheplick, 1992; Hyatt and Evans, 1998; Kobayashi and Yamamura, 2000). In such a species, genetically identical offspring are likely to undergo competition, which could reduce fitness of the maternal genotype. Spreading seed germination over time would be likely to reduce this competition (Ellner, 1986; Cheplick, 1992). There is some experimental evidence that plants may be able to regulate the dormancy of their seeds (and therefore the extent to which germination is spread over several seasons) according to the likelihood that sibling competition will occur, with larger families and seeds with lower dispersal having increased variability in germination rate (Cheplick, 1996; Hyatt and Evans, 1998).
It should be noted that the timing of germination is sensitive to a number of environmental cues experienced by the mother plant and by the seed, and thus is likely to be both plastic and variable in a number of species (Fenner, 1991; Finch-Savage and Leubner-Metzger, 2006; Graeber et al., 2012; Footitt et al., 2014). Indeed, in the presence of low environmental predictability, but where some cues about the future are available, a plastic and variable phenotype is predicted to be optimal (Wong and Ackerly, 2005). In such a scenario, a bet-hedging strategy may evolve as an increase in the variance around reaction norms for a plastic phenotype (Simons and Johnston, 1997; Simons, 2011).
In addition to germination timing, other plant phenotypes may also have been selected to have high variability within an environment. This has been proposed to be the case for metabolic pathways that produce toxins to defend against herbivory. A prolonged period of exposure of an insect population to a particular toxin can encourage the evolution of toxin resistance. It may therefore be beneficial for plants to prevent prolonged exposure of herbivores to a particular toxin by stochastically varying toxin production (Shelton, 2004). Indeed, variability in the levels of defence compounds produced by arabidopsis plants in a common environment is genetically controlled and is greater than the variability in production of primary metabolites (Joseph et al., 2015). Variability in lateral root emergence has also been hypothesized to be advantageous (Forde, 2009). However, in many cases, phenotypes might be variable between individuals because the system is poorly constrained and mechanisms to buffer against noise in development have not been selected.
Currently, little is known about the mechanisms that underlie functional variability between genetically identical plants within an environment. One possibility is that what appears to be variability within a homogeneous environment is actually due to plasticity in response to unknown microenvironmental cues (Bradshaw, 1965). Indeed, in a study where almost genetically identical L. inflata seeds were grown under growth chamber conditions, the position of seeds within Petri dishes could account for around a third of the variance in the time to germination, indicating that microenvironmental plasticity can account for at least some of the variability (Simons and Johnston, 2006). However, it is also possible that interindividual variation exists due to a lack of noise-buffering mechanisms in development, as discussed in the first section. Molecular mechanisms may also exist that amplify stochastic fluctuations to generate increased variability within a homogeneous environment. Although little is known about whether such noise-harnessing mechanisms exist in plants, they have been described for microbes and during the development of some multicellular organisms (Raj and van Oudenaarden, 2008; Eldar and Elowitz, 2010; Martins and Locke, 2015).
In microbes, it has been shown that gene-regulatory networks can amplify stochastic fluctuations to generate sub-populations of cells with distinct gene expression states (Alon, 2007; Raj and van Oudenaarden, 2008; Eldar and Elowitz, 2010; Viney and Reece, 2013; Martins and Locke, 2015). Bistability (the existence of two stable states) is often generated by positive feedback acting upon gene expression, i.e. an increase in a protein’s levels causes an increase in the production rate of that protein (Fig. 2C). In this case, the presence of low protein levels will be stable in a fraction of cells. However, stochastic fluctuations may increase protein levels above a threshold, causing a transition to a new stable state of gene expression (Choi et al., 2008). Thus, two sub-populations of cells may be generated, one with high and one with low expression levels of the protein. More complex network architectures which combine a number of interlinked network motifs can give rise to a greater number of possible stable states, allowing increased heterogeneity in the cell population (Alon, 2007). As well as possibly playing a role in creating functional interindividual differences, stochastically induced gene expression states may be used to generate differences between equivalent cells during tissue patterning or cellular differentiation. In this scenario, the realizations of the patterns produced may be variable between individuals, but functionally equivalent (Meyer and Roeder, 2014).
These examples can produce alternative stable states of gene expression; however, genetic networks also exist that can cause unstable, temporary transitions into alternative cellular states. These networks are called excitable circuits and combine a fast-acting positive feedback loop with slower acting negative feedback (Süel et al., 2006; Eldar and Elowitz, 2010; Martins and Locke, 2015). Stochastic fluctuations may trigger the positive feedback, causing a transition to an alternative state. The activation of negative feedback may then bring the system back to the original state. Since this mechanism allows independent stochastic transitions between states in each cell, it provides a means for generating variable phenotypes within a population of genetically identical cells at a given time.
Mechanisms may also exist to create cell heterogeneity in response to environmental cues, thus generating variability in a plastic phenotype. For example, in bacteria, stressful conditions such as starvation cause a large fraction of cells to undergo sporulation, in which highly resistant, dormant spores are produced (Phillips and Strauch, 2002; Schultz et al., 2009). The rest of the population may remain in an active state, allowing them to continue to grow if environmental conditions improve. In Bacillus subtilis, sporulation is thought to be triggered when the levels of a phosphorylated protein, Spo0A, exceed a particular threshold (Chastanet et al., 2010). Phosphorylation of Spo0A involves a four-step phosphorelay mechanism that is activated upon nutrient deprivation (Phillips and Strauch, 2002). Upon induction of sporulation, the activity of phosphorylated Spo0A is heterogeneous amongst cells within a population, which is proposed to explain why only a fraction of cells undergo sporulation (Chastanet et al., 2010). This heterogeneity has been proposed to be due to noise in the transfer of phosphate along the phosphorelay pathway. The presence of four steps in the phosphorelay (rather than a lower number) has been hypothesized to increase noise accumulation along the signal transduction pathway, maximizing cellular heterogeneity, and ensuring that only a fraction of the population undergoes sporulation (Chastanet et al., 2010). Although unknown, one might speculate that in plants, such a mechanism of differential regulator accumulation could underlie the discussed variability in germination timing upon the perception of permissive environmental cues by a population of seeds.
Variability can also be generated through a deterministic mechanism of differential inheritance. Again, yeast provides an example, concerning its dependence on metal ions for survival. If a cell experiences a depletion of metal ions when it undergoes budding (which generates a smaller daughter cell), most of the vacuole containing essential ions is retained by the mother (Avraham et al., 2013). In this situation, the mother cell continues to divide, whereas the daughter’s cell cycle arrests in G1 phase. This has been proposed to prevent the dilution of ions by cell division (Avraham et al., 2013). In this case, differential inheritance occurs deterministically (the vacuole is always inherited by the older mother cell) and generates two populations of cells with different behaviours. Cell fate diversity is also frequently generated by asymmetric cell divisions in the development of multicellular organisms (Macara and Mili, 2008; De Smet and Beeckman, 2011). One hypothesis to explain variability in germination time is that differential inheritance by a mother’s offspring creates inter-seed differences in a deterministic manner.
In contrast to traditional views that consider plasticity and variable phenotypes as resulting from poorly buffered development, we emphasize that specific genetic mechanisms are likely to be involved in reproducibly (i.e. robustly) delivering invariant, plastic and variable phenotypes. Each of these types of phenotype may be favoured in different ecological scenarios. Invariant phenotypes are likely to have evolved in environments where a constant phenotypic optimum exists, whereas plastic phenotypes are expected to evolve in the presence of predictable changes in the environment. When the environment varies unpredictably, variable phenotypes probably reduce the risk of extinction and promote the long-term survival of the population.
Plants of the same genotype in a given environment are likely to show more variability or plasticity for some traits than others. For example, for a hypothetical plant with a specific insect pollinator, living in an environment where the timing of rainfall is highly unpredictable and neighbouring plants often shade each other, flower morphology may be predicted to be relatively invariant, while germination might be expected to be variable and stem length plastic.
Due to their ecological significance, we argue that developmental mechanisms have probably evolved to allow the robust formation of each of the three types of phenotype. We have illustrated how a number of mechanisms can buffer development against stochastic and environmental perturbations, allowing invariant phenotypes to be produced. There are also well-characterized mechanisms that alter development in response to environmental cues. Finally, the developmental mechanisms involved in producing variable plant phenotypes are the least well understood, although QTL mapping has shown that variance in a number of phenotypes in arabidopsis is under genetic control (Hall et al., 2007). An important future step is to extend this work towards gaining a mechanistic understanding of phenotypic variability in plants.
Another important question is the extent to which plastic, variable and invariable phenotypes are beneficial in the face of climate change. This has been discussed for plasticity (Jump and Penuelas, 2005; Ghalambor et al., 2007; Chevin et al., 2010; Nicotra et al., 2010) which, as mentioned previously, may facilitate survival in a new environment, allowing the phenotype then to be optimized by natural selection (Ghalambor et al., 2007). This is likely to be the case if the phenotypic change that is induced in the new environment changes in the direction of the optimum phenotype in that environment (Ghalambor et al., 2007). However, in a changing climate, the environmental cues that trigger a plastic response may no longer be good indicators of the future environment, and this could lead to the induction of inappropriate phenotypes, which could be detrimental for survival (Ghalambor et al., 2007; Nicotra et al., 2010). In the case of changes in flowering time in response to climate change in the last century, there is evidence for the former scenario. Plant species with flowering times that are plastic with respect to temperature have undergone less rapid declines in abundance than species in which flowering time is less plastic (Willis et al., 2008). This may be due to a better ability of such species to co-ordinate their flowering time with the availability of insect pollinators or appropriate growing conditions. Phenotypic variability could also enable survival in future climates in which weather patterns are expected to become more extreme and unpredictable (IPCC, 2012). More variable phenotypes, produced as bet-hedging strategies, are likely to allow at least some individuals of a population to survive in unpredictable and extreme conditions (Cohen, 1966; Lewontin and Cohen, 1969; Simons, 2011).
We thank Hugo Tavares for useful suggestions and critical reading of the manuscript. This work was supported by the Gatsby Foundation.