Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Phenotypic Effects of Salt and Heat Stress over Three Generations in Arabidopsis thaliana

Abstract

Current and predicted environmental change will force many organisms to adapt to novel conditions, especially sessile organisms such as plants. It is therefore important to better understand how plants react to environmental stress and to what extent genotypes differ in such responses. It has been proposed that adaptation to novel conditions could be facilitated by heritable epigenetic changes induced by environmental stress, independent of genetic variation. Here we assessed phenotypic effects of heat and salt stress within and across three generations using four highly inbred Arabidopsis thaliana genotypes (Col, Cvi, Ler and Sha). Salt stress generally decreased fitness, but genotypes were differently affected, suggesting that susceptibility of A. thaliana to salt stress varies among genotypes. Heat stress at an early rosette stage had less detrimental effects but accelerated flowering in three out of four accessions. Additionally, we found three different modes of transgenerational effects on phenotypes, all harboring the potential of being adaptive: heat stress in previous generations induced faster rosette growth in Sha, both under heat and control conditions, resembling a tracking response, while in Cvi, the phenotypic variance of several traits increased, resembling diversified bet-hedging. Salt stress experienced in earlier generations altered plant architecture of Sha under salt but not control conditions, similar to transgenerational phenotypic plasticity. However, transgenerational phenotypic effects depended on the type of stress as well as on genotype, suggesting that such effects may not be a general response leading to adaptation to novel environmental conditions in A. thaliana.

Introduction

In the face of climate change, many organisms may be forced to adapt to novel, potentially challenging environmental conditions that may often exceed their typical range of reaction [1]. Alternatively, organisms may have to migrate and track suitable habitats to escape extinction. Changing environmental conditions, such as increasing temperature, may thus have important consequences for species diversity in natural ecosystems [2], but also for crop production, since increasing temperatures can significantly decrease yield in many crop species [3]. Another major environmental factor affecting agriculture and natural vegetation is soil salinity, to date a problem in more than 100 countries worldwide [4]. To assess the potential of species to cope with such environmental challenges it is on the one hand important to know the species’ norms of reaction, i.e. the phenotypic response within one generation, and on the other hand the species’ ability to adapt over a number of generations to changing environmental conditions.

Molecular and phenotypic effects of heat or salt stress are well studied in Arabidopsis thaliana. Heat stress is known to affect gene expression [5,6] and many studies have assessed the roles of heat shock proteins in stress response [7,8]. Generally, heat stress is expected to have negative fitness consequences [9-11], and the timing of heat stress may be most important: heat stress during seed production can have highly detrimental fitness effects [9,10], whereas similar conditions at vegetative plant stages may have much milder fitness consequences [12]. Similarly, knowledge about molecular changes induced by salt stress is increasing [13,14] and a number of studies have examined the phenotypic consequences of saline conditions in A. thaliana [15-17]. Compared to other species, A. thaliana is highly sensitive to saline conditions [18], although remarkable differences between A. thaliana accessions have been reported [16,17].

While the consequences of environmental stress on phenotypes and molecular processes within generations are relatively well explored in A. thaliana, much less is known about transgenerational effects. A number of studies have suggested that later generations can be influenced by environmental conditions experienced by preceding generations [9,19-21]. However, the underlying mechanisms often remain elusive.

If genetic variability can be excluded, epigenetic inheritance is one candidate cause for environmentally induced transgenerational effects. Epigenetic inheritance encompasses anything heritable apart from DNA, although usually, three main mechanisms are distinguished: DNA methylation, histone modifications and inheritance of small RNA molecules [22-24]. Importantly, epigenetic patterns are not only heritable, but can also directly be influenced by the environment. For example, in apomictic dandelion, Verhoeven et al. [25] found altered DNA methylation upon environmental stress exposure, which was mostly heritable to the next generation. However, apart from epigenetic effects, environmental stress may also destabilize the genome and lead to an increase of genomic mutations [26], e.g. through mobilization of transposable elements in response to heat stress [27] or an increase of homologous recombination upon salt stress [28]. With phenotypic data alone it is difficult to disentangle genetic and epigenetic transgenerational effects of environmental stress. However, a relatively high number of genetically identical, independent replications and correlations between parent and offspring phenotypes may reveal whether observed transgenerational phenotypic effects are genetically or epigenetically inherited.

If genetic diversity can be excluded, heritable epigenetic variability could significantly contribute to phenotypic diversity and thus to the adaptive potential of a species [29]. Furthermore, epigenetic modifications can evolve much faster than genetic mutations and may be especially relevant when rapid adaptation is required [30,31]. Additionally, resetting of epigenetic marks would allow an organism to react to variable conditions [26,32]. The resulting transgenerational phenotypic effects of environmental stress could in principle fall into three categories: tracking, phenotypic plasticity and bet-hedging [33]. Tracking is typically understood as a change in mean of a phenotypic trait due to altered natural selection in a novel environment, resulting in a shift of allele frequencies within a population [33]. In the absence of genetic variability and selection on existing genetic variation, a similar effect could be observed if environmental stress induces heritable epigenetic alterations that lead to a shift in mean phenotype over generations. Phenotypic plasticity could similarly result from transgenerational epigenetic inheritance triggered by environmental stress experienced in preceding generations, but one would then expect a plastic interaction with the current environment, i.e. a phenotypic shift due to ancestral treatment could be observed in one present environment, but not in another. Phenotypic plasticity is advantageous when environments are variable, but can be predicted reliably within an organism’s lifetime [34,35]. However, if environments cannot be predicted adequately, e.g. due to rapid fluctuations or missing detectability from organisms, diversified bet-hedging can be adaptive [33,36]. Diversified bet-hedging refers to a strategy where one genotype produces phenotypically highly variable offspring, independent of environmental conditions. This way the risk is spread among offspring and a fraction of the offspring may express an appropriate phenotype, decreasing the variance in offspring fitness over different environments [33,37].

To what extent plants can express such transgenerational phenomena and what environmental conditions are necessary to trigger them, is presently unclear. A major problem is that different studies often apply environmental stresses differently, which may prevent meaningful comparisons among studies. For example, heat shock response has been studied at temperatures ranging from 38 °C to over 50 °C in A. thaliana [38], and it is not surprising that results of such studies are often incongruent. Additionally, the applied environmental stresses are often not mimicking natural conditions, such as in the case of sudden shifts in temperature to achieve maximal heat shock [39,40]. Not surprisingly, results of such studies can differ significantly from studies applying more realistic scenarios, such as for example gradual temperature increases [6]. This raises the question as to whether results from unrealistic experimental designs may uncover adaptive responses or instead report artifacts induced by unnatural stress conditions.

The goal of the present study was to assess the phenotypic responses of A. thaliana to realistic heat and salt stress conditions. To get an impression of how consistent plants responded to stresses within and across generations, we tested four different accessions, each with a relatively high number of replications over three generations of stress or control treatment. Specifically, we addressed the following questions: How does A. thaliana react to (realistic) environmental stress conditions and are there differences between genotypes? Does ancestral stress treatment influence phenotypes in stress or control treatments over three generations, and if yes, what are the modes of response (tracking, phenotypic plasticity, bet-hedging)? If three generations of stress lead to transgenerational phenotypic changes, are these changes adaptive?

Materials and Methods

Plant material and growth conditions

Four widely studied Arabidopsis thaliana accessions were used: Columbia (Col-0), Landsberg erecta (Ler-0), Cape Verde islands (Cvi-0) and Shahdara (Sha-0; hereafter referred to as Col, Ler, Cvi and Sha, respectively). One plant per accession was grown in generation 0 (G0) to minimize genetic diversity and all plants grown in generation 1 (G1) were descendants of these G0 plants.

Plants were grown in individual 7 x 7 x 8 cm pots filled with Biouniversalerde (Oekohum GmbH, Herrenhof, Switzerland), an all-purpose soil without peat. Pots were randomly arranged on 28-pot-trays (G0, G1) or 24-pot-trays (generations 2 and 3 (G2, G3)) that were randomized three times per week to avoid position effects. Randomization of trays was stopped upon maturation of siliques. Approximately five seeds were sown per pot to ensure successful germination. Throughout the experiment, plants were grown in climate chambers (Kälte 3000, Landquart, Switzerland) to equalize growth conditions over generations. Seeds were stratified at 4 °C in the dark for five days to break seed dormancy. After stratification, all plants were moved to climate chambers, and this day was counted as day 0 of the experiment for each generation. About one week after germination, seedlings were thinned to one plant per pot. Upon flowering, individual pots were packed in Arabisifter floral sleeves (Lehle Seeds, Round Rock, Texas, USA) to avoid cross-pollination between neighboring plants and to harvest seeds. Plants were watered once per week with tap water containing Solbac (Andermatt Biocontrol AG, Grossdietwil, Switzerland), diluted according to manufacturers instructions. Watering was stopped two weeks before seed harvest.

In G0, all plants were grown under control conditions (10 kLux light for 16 h, dark for 8 h; 22 °C/18 °C day/night temperatures; 50 %/60 % relative humidity at day/night). In G1 – G3, plants were either grown in control conditions or in one of two stress conditions.

Heat stress

To find the best setting for heat treatments, a pre-experiment was conducted. Col plants were heat treated at 32 °C, 36 °C or 40 °C for 1, 2 or 3 days, either 7 or 14 days after germination. Plants that were heat treated for more than 1 day were allowed to recover for 2 - 3 days between heat treatments. Overall, 18 different heat treatments with eight replications each were tested (for details see Table S1). As even the most extreme heat treatment (three days at 40 °C) had no major impact on plants fitness (data not shown) we decided to expose plants to heat stress at 40 °C for three successive days, starting at day 12 of the experiment. For this heat treatment, temperatures were gradually increased during the day over 7 h to reach 40 °C. This temperature was kept for 2 h followed by a gradual decrease over 7 h back to 18 °C. Light and relative humidity were identical to control conditions, as were night conditions. After three days of heat stress, plants were continued to be grown under control conditions.

Salt stress

A pre-experiment was conducted to identify suitable salt concentrations to apply salt stress. Eight replications per treatment of genotype Col were watered for the first four weeks with NaCl-solution with concentrations of 25 mM, 50 mM, 100 mM or 150 mM. 25 mM NaCl had no major effect on plant fitness, while none of the plants watered with 100 mM or 150 mM NaCl-solution survived to set seeds. Plants treated with 50 mM NaCl-solution showed somewhat reduced fitness, but still survived and produced seeds (data not shown). Therefore, we decided to water plants in the salt treatment with 50 mM NaCl for the first four weeks of each generation, followed by normal watering thereafter. In G2 of the full experiment, only one plant of genotype Cvi survived this salt treatment, therefore this genotype had to be excluded from the salt experiment.

Experimental design

In G1, 25 replications per accession were grown either under heat, salt or control conditions, and propagated by single seed descent to the next generation, G2, where plants were grown under the same conditions. In most accession × treatment combinations, not all 25 replicated plants of G1 survived to produce seeds. Missing lines were replaced in generation 2 (G2) by taking seeds of randomly chosen G1 lines such that again 25 plants per accession and treatment were grown in G2. The same procedure was applied for G3. Care was taken to maintain as many independent replications as possible in G3 (see Table S2). Due to space restrictions, it was not possible to grow salt and heat treated plants at the same time after G1. Consequently, control plants were also grown separately (together with each stress treatment) as from G2. Offspring of G2 plants were grown both under stress and under control conditions. Therefore, in G3, plants with four different histories were grown (see Figure 1): stress treatment in G1 – G3 (hereafter abbreviated SSS for salt stress, HHH for heat stress), control treatment in G1 and G2 followed by stress treatment in G3 (CCS, CCH), stress treatment in G1 and G2 and followed by control treatment in G3 (SSC, HHC) and control treatment in G1 – G3 (CCC).

thumbnail
Figure 1. Setup of the experiment.

G0 = generation 0, G1 = generation 1 etc. In blue, SSS refers to three generations of salt treatment, CSS indicates that the first generation was grown under control conditions followed by two generations in stress treatment etc., while in red HHH refers to three generations of heat treatment, etc. In G0, only one individual was grown, whereas in G1 to G3, 25 individuals per treatment were grown. Plants were propagated by single seed descent, and in case of mortality, missing plants were restocked by duplicating other lines of the same treatment.

https://doi.org/10.1371/journal.pone.0080819.g001

To estimate growth, phenological parameters, fitness and plant architecture, a large number of phenotypic traits, including flowering time, number of rosette leaves, and total number of siliques, were measured in G3 (Table S3), with a subset of traits also measured in G1 and G2. Some of these phenotypic traits best describe variation in life history traits among genotypes [41].

Data analysis

All statistical analyses were conducted using R [42]. As heat and salt treatments had very different effects on phenotypic traits, data were analyzed separately for heat and salt treatments. Means and standard deviations of all phenotypic traits computed separately for each genotype and treatment combination (HHH, CCH, etc.) are presented in Table S4 and S5.

Selection of phenotypic traits – Many of the measured traits strongly correlated with each other, e.g. diameter and rosette leaves after two and after three weeks (Figure 2). We were primarily interested in how many independent, i.e. non-correlated traits were affected by the applied stress treatment, as they may also be independent from each other on a molecular level. To select non-correlated traits, pairwise Pearson correlations were computed and plotted to a heatmap, where traits were reordered according to their similarities determined with a dendrogram. From the resulting clusters of traits, five traits that were both responsive to the applied stress and not correlated with each other were selected both for salt and heat stress (Figure 2).

thumbnail
Figure 2. Correlation analysis to select a subset of independent phenotypic traits.

The heat map shows pairwise Pearson correlations between phenotypic traits (darker colors denote stronger correlations), with data from the heat stress experiments above the diagonal in red and data from the salt stress experiment below the diagonal in blue. Phenotypic traits written in bold letters were selected for statistical analyses.

https://doi.org/10.1371/journal.pone.0080819.g002

Effect of genotypes- initially we performed a linear mixed effect model with genotype, previous treatments (G1 and G2) and G3 treatments as well as their pairwise interactions as fixed and tray as random factors. As genotype had a dominant effect on all measured traits (Tables S6 and S7), principal component analyses with the full set of phenotypic traits were performed to further assess the effect of genotypes. Both for the heat and the salt data the genotypes clearly clustered separately (Figure 3 and Figure S1), thus all further analyses were performed for individual genotypes.

thumbnail
Figure 3. Principal component analysis of data from the heat experiment.

Colors denote different genotypes, indicating their distinct phenotypes.

https://doi.org/10.1371/journal.pone.0080819.g003

Effect of stress treatments – Linear mixed effect models with previous treatment (G1 and G2) and G3 treatment as well as their interactions as fixed and tray as random factor were computed separately for each genotype. P-values were adjusted for multiple testing separately for each genotype following Benjamini and Hochberg [43]. Where significant effects of previous (G1 and G2) treatments or interactions between previous and G3 treatments were observed, effects of previous treatments were calculated separately for each G3 treatment. To assess the effect of previous treatments on the variance of phenotypic traits, F tests were performed separately for each G3 treatment. As trays had a significant effect on several traits, we first estimated tray effects using a linear model and used the residuals of these tests for all further calculations. P-values were adjusted for multiple testing following Benjamini and Hochberg [43].To test whether increased variance could be explained by inheritance of extreme phenotypes in G2, a Spearman correlation for a subset of traits was calculated between G3 plants and their progenitors grown in G2. To assess whether increased variance may allow for different life history strategies (bet-hedging), five Cvi plants with the highest fitness (number of siliques) in treatments HHH, HHC, CCH and CCC were identified. Plants with the same ancestral treatment were compared with a linear model to test whether traits differed between G3 treatments (e.g. HHH vs. HHC). Similarly, the five plants with the lowest number of siliques were compared, as well as all Cvi plants independent of silique numbers.

Results

Phenotypic effects of environmental stress in generation 3

Salt treatment – Exposure to salt stress generally reduced fitness in all genotypes when compared to plants grown under control conditions (Table 1). Plants were shorter and had reduced numbers of branches in all genotypes. Additionally, in genotypes Col and Ler, the number of rosette leaves at FFD, a trait that indicates at what physiological age the transition from vegetative to reproductive stage takes place, was increased under salt conditions, indicating a delay of flowering when compared to control plants (Figure 4A, Table 2). Overall, the response of genotype Sha to salt treatment appeared weaker than in the other genotypes, with no significant delay in flowering and a relatively small decrease in total branches (Figure 4A, Table 2).

G3 treatmentG1G2 treatmentG3 x G1G2
GenotypePhenotypic traitFdFPbFdFPbFdFPb
ColRosette diameter 3 weeks5.4761,80.0950.0281,730.8681.6781,730.266
Rosette leaves FFD9.2661,80.032*3.9071,730.0690.4371,730.511
Height45.2671,8<0.001***2.1591,720.1951.2781,720.262
Siliques per branch2.1351,80.3071.4631,720.3070.0421,720.838
Total branches18.4251,80.005**0.0741,730.9670.0021,730.967
LerRosette diameter 3 weeks6.5651,80.0670.0021,820.9650.4481,820.673
Rosette leaves FFD29.0301,80.001**0.3031,820.5831.8631,820.235
Height13.0211,80.014*0.0011,810.9711.9681,810.219
Siliques per branch0.3121,80.8970.0561,780.8970.0171,780.897
Total branches20.9871,80.004**0.3211,820.5720.3951,820.572
ShaRosette diameter 3 weeks2.9931,80.2440.1851,840.6680.2091,840.668
Rosette leaves FFD1.5861,80.3240.1971,840.6591.4251,840.324
Height47.1871,8<0.001***4.3151,840.0540.3861,840.536
Siliques per branch0.4801,80.5082.5641,840.1516.4741,840.026*
Total branches14.0241,80.011*0.5201,840.5040.4511,840.504

Table 1. Effects of G3 and G1G2 salt treatments and interactions on phenotypic traits were calculated for three genotypes using linear mixed models with trays as random factor.

bP-values were adjusted for multiple testing according to Benjamini and Hochberg [43], separately for each genotype.
*** P-value <0.001, **P-value < 0.01, *P-value < 0.05, ⋅ P-value < 0.1
CSV
Download CSV
thumbnail
Figure 4. Effect of stress treatments on rosette leaves at FFD.

A: Col, Ler and Sha grown in G3 under salt (blue) and control (white) conditions. B: Col, Ler, Sha and Cvi grown in G3 under heat (red) and control (white) conditions. ** P-value < 0.01; * P-value < 0.05, · P value < 0.1, ns: P-value not significant. P-values were adjusted for multiple testing according to Benjamini and Hochberg [43].

https://doi.org/10.1371/journal.pone.0080819.g004

G3 treatmentG1G2 treatmentG3 x G1G2
GenotypePhenotypic traitFdFPbFdFPbFdFPb
ColRosette leaves 3 weeks0.4461,80.6982.3631,650.2580.0001,651.000
Rosette leaves FFD13.5461,80.012*0.6731,650.5530.0501,650.824
Height0.1171,80.7410.8311,630.4871.2751,630.487
Total siliques0.0031,80.9580.5491,650.6151.4041,650.481
Mean length siliques6.4131,80.0700.1171,640.8270.0481,640.827
LerRosette leaves 3 weeks0.4451,80.7950.0681,520.7950.2061,520.795
Rosette leaves FFD5.9001,80.0830.6021,550.5880.0301,550.863
Height0.3551,80.5680.6521,550.5641.1921,550.559
Total siliques0.3951,80.5473.6111,550.1250.5781,550.547
Mean length siliques0.5591,80.6350.0041,550.9492.9421,550.184
ShaRosette leaves 3 weeks4.7701,80.08110.6301,810.003**0.1551,810.695
Rosette leaves FFD14.9221,80.010**1.9951,820.2151.3141,820.255
Height0.3661,80.5621.7361,820.3831.1351,820.386
Total siliques0.1111,80.9770.0371,820.9770.0011,820.977
Mean length siliques0.1711,80.6900.2281,820.6900.1691,820.690
CviRosette leaves 3 weeks1.6571,60.4910.3231,420.7640.0791,420.780
Rosette leaves FFD1.1911,60.3870.7631,420.3870.7751,430.387
Height0.6311,60.7920.0081,420.9280.2881,430.792
Total siliques0.2771,60.6180.4941,420.6184.1871,430.094
Mean length siliques0.0191,60.8960.6091,420.8790.0881,420.896

Table 2. Effects of G3 and G1G2 heat treatments including interactions on phenotypic traits separately analyzed for each genotype using linear mixed models with trays as random factor.

bP-values were adjusted for multiple testing according to Benjamini and Hochberg [43] separately for each genotype.
** P-value < 0.01, *P-value < 0.05, ⋅ P-value < 0.1
CSV
Download CSV

Heat treatment – After three days of heat treatment at an early rosette stage genotypes Col, Sha and Ler (trend) had fewer rosette leaves at first flowering day (FFD) when compared to control treatments, while no difference was found in genotype Cv (Table 2, Figure 4B).

Phenotypic effects of environmental stress across generations

Salt treatment –, We found interactions between previous treatments and G3 treatments for the number of siliques per branch in genotype Sha (Table 1). Specifically, under G3 salt conditions offspring of salt treated plants (SSS) produced more siliques per branch than plants who’s progenitors were grown under control conditions (CCS, Fdf = 7.1911,41, P = 0.011), while no differences were found under G3 control conditions between SSC and CCC (Fdf = 0.5871,43, P = 0.448; Figure 5B). Furthermore, SSS plants grew taller than CCS plants (Fdf = 4.2441,41, P = 0.046), while SSC and CCC plants grew similarly high (Fdf = 1.071,43, P = 0.307), indicating that the plants architecture was altered depending both on previous and G3 treatments.

thumbnail
Figure 5. Transgenerational effects of two stress generations displaying three different modes of adaptation to novel environments.

A: Tracking: in genotype Sha two generations of ancestral heat treatment (orange) led to more rosette leaves after three weeks in comparison to ancestral control treatment (white). B: Transgenerational phenotypic plasticity: in genotype Sha, two generations of ancestral salt treatment (blue) led to increased number of siliques per branch when compared to ancestral control conditions (white) under G3 salt conditions, but not G3 control conditions. C: Diversified bet-hedging: in genotype Cvi ancestral heat treatment (orange) increased the variance of total siliques in comparison to ancestral control treatment (white).

https://doi.org/10.1371/journal.pone.0080819.g005

Previous salt treatments did not affect phenotypic traits measured in genotypes Col and Ler, nor did we find interactions between previous and G3 treatments for these genotypes.

While ancestral treatment had an effect on the variance of a few traits (Table S8), no genotype showed increased variance in more than one trait per G3 treatment, suggesting that ancestral salt stress overall had no major impact on phenotypic variance. Nevertheless, when we did find an effect, variance was usually increased in plants growing in a novel environment compared to the ancestral treatment, i.e. SSC in comparison to CCC; and CCS in comparison to SSS.

Heat treatment – In genotype Sha, heat treatment in G1 and G2 increased the number of rosette leaves at day 21 in G3, both in heat and control treatment (Table 2, Figure 5A), while none of the other genotypes showed a similar effect (Table 2).

Two generations of ancestral heat treatment increased the variance of a number of traits for Cvi, both under G3 control and heat conditions when compared to plants with ancestral control treatment. Under both treatments these effects were moderate (p-values not significant after correction for multiple testing; Table 3; Figure 5C). None of the other genotypes showed increased variance for more than one trait (Table S9). Where increased variance was observed we tested whether G3 values correlated with values of their progenitors grown in G2. However, no such correlations could be observed, indicating that increased variance in G3 was independent of phenotypes in G2 (data not shown).

G3 treatmentPhenotypic traitFdFPPb
HeatRosette leaves d210.27611,100.046*0.114
Rosette leaves FFD1.37011,100.6270.697
Height0.88511,100.8390.839
Total siliques0.25211,100.033*0.110
Mean silique length0.21711,100.019*0.095
ControlRosette leaves d210.40016,120.089·0.178
Rosette leaves FFD0.68916,120.4790.599
Height0.56216,120.2800.400
Total siliques0.23216,120.008 **0.078·
Mean silique length0.55116,120.2640.400

Table 3. Effect of two generations of heat treatment (G1 and G2) vs two generations of control treatment on variances of traits measured under G3 heat and control conditions in genotype Cvi.

bP values were adjusted for multiple testing according to Benjamini and Hochberg [43].
** P value < 0.01; * P value < 0.05; ·:P value < 0.1
CSV
Download CSV

We selected five Cvi plants with the highest fitness (number of siliques), each for HHH, HHC, CCH and CCC, to test whether increased variance due to ancestral heat stress (HHH and HHC) would allow different life history strategies to appear and be successful in different G3 environments. Of these five plants, those with ancestral heat treatment differed significantly in number of rosette leaves at FFD when grown in G3 heat treatment (HHH, few rosette leaves at FFD) compared to G3 control treatment (HHC, many rosette leaves at FFD; Figure 6B, left) indicating that through increased variance contrasting life history strategies could lead to highest fitness in different environments. No such effects were observed in plants with ancestral control treatment (Figure 6B, right), when the five plants with the lowest fitness were compared with each other (Figure 6C) or when the whole data set was used (Figure 6A).

thumbnail
Figure 6. Transgenerational effects of heat on rosette leaves at FFD in genotype Cvi indicate diversified bet-hedging.

A: Entire data-set. No differences in number of rosette leaves at FFD were found between G3 treatments when plants experienced the same ancestral treatment. B: For each G1G2 × G3-treatment combination the five plants with the highest silique number were chosen. In plants with ancestral heat treatment, number of rosette leaves differed significantly between G3 heat (HHH) and G3 control (HHC) treatment. When plants were control treated in G1 and G2, no differences were found between G3 treatments. C: For each G1G2 × G3-treatment combination the five plants with the lowest silique number were chosen. No differences in numbers of rosette leaves were found between G3 plants when plants experienced the same ancestral treatments.

https://doi.org/10.1371/journal.pone.0080819.g006

Discussion

Phenotypic responses to environmental stress within generations

Salt –Salt stress clearly reduced plant fitness, but genotypes differed substantially in the strength of this response. Genotypes Col and Ler responded strongly in most observed traits, whereas in genotype Sha fewer traits were responsive, and overall fitness seemed less affected, similar to previous studies [16,17]. Comparable to other studies [15,44] we also observed delayed flowering under saline conditions for Col and Ler. However, in Sha flowering was not delayed, making this genotype an interesting candidate to further study responses to saline conditions.

Heat – Overall the phenotypic response to the applied heat treatment was weak, without decreasing overall fitness of heat treated plants. However, in three of the four genotypes, flowering was accelerated, similar to findings of Balasubramanian et al. [45], who grew about 100 A. thaliana accessions in slightly elevated temperatures (27 °C) and found induction of flowering due to elevated temperatures, although this effect varied extensively among accessions. Amongst the accession with moderate to strong induction were Col, Ler and Sha, corroborating our results. They further found that accessions with non-functional FRIGIDA (FRI) or FLOWERING LOCUS C (FLC) alleles showed induction of flowering, whereas accessions with functional FLC were insensitive to elevated temperatures. One exception was genotype Cvi, which has a non-functional FRI but nevertheless expresses FLC and as a result was insensitive to elevated temperatures. Because we also observed that Cvi was unresponsive to high temperatures, we assume that the same mechanism as described by Balasubramanian et al. [45] and Balasubramanian and Weigel [46] also induced flowering of Cvi in our experiments. However, there are important differences in our experimental setup compared to these earlier experiments. Balasubramanian et al. [45] grew their plants throughout the growth phase of the plant at only slightly elevated temperatures (27 °C) and short day conditions, whereas our plants grew mainly under control conditions (long day) and we applied heat stress only during three days at an early rosette stage by gradually increasing temperature to 40 °C to mimic hot days in natural environments. Such extreme temperatures are generally regarded as heat shock conditions for A. thaliana, although a gradual increase of temperature may affect plants differently than sudden exposure to high temperatures [6]. While it was proposed that heat shock may not induce flowering [45], we find that a relatively short exposure to high temperatures may suffice to effectively induce flowering in A. thaliana.

Phenotypic effects of environmental stress across generations

Organisms can respond to changing environments in three different ways: a) tracking, i.e. a gradual change of mean trait values in response to altered natural selection in novel environments, b) phenotypic plasticity, i.e. different means of traits in different environments, or c) diversified bet-hedging, i.e. an increase of variance of a trait independent of the environment [33]. In our study we detected evidence for all three modes over the short time of three generations: ancestral heat conditions accelerated growth in genotype Sha, independent of G3 conditions, resembling a tracking response, while interaction between ancestral and G3 salt and control conditions led to increased number of siliques per branch in genotype Sha, comparable to transgenerational phenotypic plasticity. Furthermore, ancestral heat treatment increased the variance of a number of phenotypic traits in genotype Cvi, allowing for contrasting life history strategies to be successful under different G3 conditions, similar to diversified bet-hedging: while highest fitness under G3 heat conditions was reached by the fraction of offspring of heat treated plants that flowered earliest (i.e. with very few rosette leaves at FFD), the opposite fraction, i.e. the one with most rosette leaves at FFD, was fittest under G3 control conditions. Such different life history strategies were not found in offspring of control plants. However, it is important to emphasize that all these observations depended both on the type of stress applied as well as on the studied genotype, e.g. neither Ler nor Col exhibited any transgenerational phenotypic responses, suggesting that our findings of short-time phenotypic alterations may not be a general response providing adaptation to novel environments in A. thaliana.

Other studies addressing similar questions have found different traits to be responsive to transgenerational stress, however, we think that most of this can be assigned to different experimental designs. E.g. Whittle et al. [9] found higher fitness upon heat exposure due to ancestral heat stress, however, they used a design where plants were exposed to heat during seed development, which might affect plants very differently than the heat stress we applied during early vegetative growth. Similarly, Suter and Widmer [21] and Boyko et al. [47] found improved growth under salt conditions due to ancestral salt stress, suggesting an acquired tolerance to salt. Boyko et al. [47] used a very different approach by stressing plants in petri dishes, whereas Suter and Widmer [21] used more generations of stress, potentially reinforcing effects not observable after only three generations of stress.

By using a relatively large number of genetically virtually homozygous replications, we could exclude genetic diversification between stressed and control lines as important source of phenotypic variation. This suggests that a different mode of inheritance underlies the observed phenotypic changes in our study, potentially involving epigenetic processes. The finding of interaction between ancestral and G3 treatment could further corroborate the hypothesis of epigenetic inheritance: ancestral treatments could alter epigenetic patterns, which then are reversible upon exposure to present treatments (SSS differs from CCS, but SSC not from CCC in number of siliques per branch). Like in other studies suggesting similar epigenetic involvement in transgenerational phenotypic inheritance [9,21], further studies involving molecular analyses would be required to identify the molecular mechanisms of inheritance.

Adaptation to environmental stress is often divided into tolerance and avoidance strategies. The former requires a plant to endure certain environmental conditions, whereas in the latter case a plant needs to sense adverse conditions and find a way to circumvent them. While the adaptive value of the transgenerational phenotypic changes observed in our climate chamber experiments is difficult to assess due to artificial conditions, we nevertheless propose that the observed effects of transgenerational phenotypic plasticity may be similar to avoidance strategies. SSS plants had more siliques per branch and grew taller than CCS plants, suggesting that the plants architecture underwent significant changes due to ancestral salt treatment. Similar phenotypic characteristics have been found to increase seed dispersal [48], which is one way how plants can affect their environment: if a plant grows under harsh conditions, dispersing seeds as far as possible may enhance chances of offspring to experience a more favorable environment [49]. Therefore, this altered architecture may help offspring of SSS plants to avoid saline conditions if the environment is sufficiently patchy and may therefore ultimately increase fitness.

The adaptive value of the tracking response is even more difficult to assess, as normally selection would be regarded as the driver of tracking [33], and no selection was consciously applied in this study. Nonetheless, we hypothesize that increased rosette growth due to ancestral heat treatment might indicate an acquired tolerance to heat stress, e.g. through stabilization of genes involved in the functioning of photosynthesis [50], known to be highly sensitive to heat stress [51]. This would allow the plant to grow faster even under elevated temperatures. Further studies taking e.g. the stability of Rubisco into account could clarify the underlying mechanisms.

The diversified bet-hedging resulting from ancestral heat treatment observed for genotype Cvi could be adaptive if different environmental conditions require contrasting phenotypes to reach maximal fitness, and if environments cannot be predicted reliably. This would ensure that at least a fraction of the offspring expresses the right phenotype in the encountered environment [33,36]. As the increased variance could not be explained by inheritance of extreme phenotypes in G2, we hypothesize that ancestral heat treatment may have destabilized epigenetic patterns. Increased epigenetic variation has recently been associated with a wider range of phenotypes [52], matching our observations. While bet-hedging is an evolutionary strategy often found under natural conditions in plants, e.g. in the timing of seed germination [37], de novo evolution of bet-hedging has, to our knowledge, so far mainly been studied in bacteria [36,53]. Further studies will have to establish under what conditions evolution or induction of bet-hedging can occur in plants and what molecular mechanisms are involved.

In conclusion, we found that salt stress decreases fitness and delays flowering, while the exposure to heat early in the life cycle can accelerate flowering, although genotypes differ in the extent of their responses to stress. We also detected three different modes by which A. thaliana can respond to environmental change over the short time scale of three generations, although these were highly dependent on the type of stress as well as on the genotypes. We suggest that even in the absence of genetic variability, A. thaliana still maintains a surprisingly high potential to react to environmental stress and that these effects may at least partly be inherited to subsequent generations, potentially through epigenetic mechanisms.

Supporting Information

Figure S1.

Principal component analysis of the salt experiment, colored according to genotype.

https://doi.org/10.1371/journal.pone.0080819.s001

(PDF)

Table S1.

Setup of the pre-experiment to determine appropriate heat stress conditions.

https://doi.org/10.1371/journal.pone.0080819.s002

(DOCX)

Table S2.

Numbers of independent replicated lines that survived from G1 to G3 and the total number of plants that survived in G3.

https://doi.org/10.1371/journal.pone.0080819.s003

(DOCX)

Table S3.

Phenotypic traits measured in generation 3 (G3).

https://doi.org/10.1371/journal.pone.0080819.s004

(DOCX)

Table S4.

Mean ± standard deviation of phenotypic traits measured for four genotypes under heat and control conditions in G3.

https://doi.org/10.1371/journal.pone.0080819.s005

(DOCX)

Table S5.

Mean ± standard deviation of phenotypic traits measured for three genotypes under salt and control conditions in G3.

https://doi.org/10.1371/journal.pone.0080819.s006

(DOCX)

Table S6.

Effect of genotype, G3 heat treatment and G1G2 heat treatment and their pairwise interactions on phenotypic traits.

https://doi.org/10.1371/journal.pone.0080819.s007

(DOCX)

Table S7.

Effect of genotype, G3 salt treatment and G1G2 salt treatment and their pairwise interactions on phenotypic.

https://doi.org/10.1371/journal.pone.0080819.s008

(DOCX)

Table S8.

Effect of two generations of salt treatment (G1 and G2) vs. two generations of control treatment on variances of traits measured under G3 salt and control conditions.

https://doi.org/10.1371/journal.pone.0080819.s009

(DOCX)

Table S9.

Effect of two generations of heat treatment (G1 and G2) vs. two generations of control treatment on variances of traits measured under G3 heat and control conditions.

https://doi.org/10.1371/journal.pone.0080819.s010

(DOCX)

Acknowledgments

We thank Maja Frei for handling the plants and Sabine Güsewell for statistical advice. Many people were involved in harvesting and measuring the plants: many thanks to Franziska Berger, Bryan Berger, Ana Florez, Maja Frei, Sonja Hassold, Céline Janz, Margot Paris, Sibil Rentsch, Adriano Spiccia and Niklaus Zemp.

Author Contributions

Conceived and designed the experiments: LS AW. Performed the experiments: LS. Analyzed the data: LS. Contributed reagents/materials/analysis tools: LS AW. Wrote the manuscript: LS AW.

References

  1. 1. Visser ME (2008) Keeping up with a warming world; assessing the rate of adaptation to climate change. Proc R Soc Lond B 275: 649-659. doi:https://doi.org/10.1098/rspb.2007.0997. PubMed: 18211875.
  2. 2. Dawson TP, Jackson ST, House JI, Prentice IC, Mace GM (2011) Beyond predictions: biodiversity conservation in a changing climate. Science 332: 53-58. doi:https://doi.org/10.1126/science.1200303. PubMed: 21454781.
  3. 3. Hedhly A (2011) Sensitivity of flowering plant gametophytes to temperature fluctuations. Environ Exp Bot 74: 9-16. doi:https://doi.org/10.1016/j.envexpbot.2011.03.016.
  4. 4. Rengasamy P (2006) World salinization with emphasis on Australia. J Exp Bot 57: 1017-1023. doi:https://doi.org/10.1093/jxb/erj108. PubMed: 16510516.
  5. 5. Zeller G, Henz SR, Widmer CK, Sachsenberg T, Rätsch G et al. (2009) Stress-induced changes in the Arabidopsis thaliana transcriptome analyzed using whole-genome tiling arrays. Plant J 58: 1068-1082. doi:https://doi.org/10.1111/j.1365-313X.2009.03835.x. PubMed: 19222804.
  6. 6. Larkindale J, Vierling E (2008) Core genome responses involved in acclimation to high temperature. Plant Physiol 146: 748-761. PubMed: 18055584.
  7. 7. Vierling E (1991) The roles of heat shock proteins in plants. Annu Rev Plant Biol 42: 579-620. doi:https://doi.org/10.1146/annurev.pp.42.060191.003051.
  8. 8. Wang W, Vinocur B, Shoseyov O, Altman A (2004) Role of plant heat-shock proteins and molecular chaperones in the abiotic stress response. Trends Plant Sci 9: 244-252. doi:https://doi.org/10.1016/j.tplants.2004.03.006. PubMed: 15130550.
  9. 9. Whittle CA, Otto SP, Johnston MO, Krochko JE (2009) Adaptive epigenetic memory of ancestral temperature regime in Arabidopsis thaliana. Botany 87: 650-657. doi:https://doi.org/10.1139/B09-030.
  10. 10. Zinn KE, Tunc-Ozdemir M, Harper JF (2010) Temperature stress and plant sexual reproduction: uncovering the weakest links. J Exp Bot 61: 1959-1968. doi:https://doi.org/10.1093/jxb/erq053. PubMed: 20351019.
  11. 11. Tonsor SJ, Scott C, Boumaza I, Liss TR, Brodsky JL et al. (2008) Heat shock protein 101 effects in A. thaliana: genetic variation, fitness and pleiotropy in controlled temperature conditions. Mol Ecol 17: 1614-1626. doi:https://doi.org/10.1111/j.1365-294X.2008.03690.x. PubMed: 18321256.
  12. 12. Wassmann R, Jagadish SVK, Heuer S, Ismail A, Redona E et al. (2009) Climate change affecting rice production: the physiological and agronomic basis for possible adaptation srategies. Adv Agron 101: 59-122. doi:https://doi.org/10.1016/S0065-2113(08)00802-X.
  13. 13. Galpaz N, Reymond M (2010) Natural variation in Arabidopsis thaliana revealed a genetic network controlling germination under salt stress. PLOS ONE 5: e15198. doi:https://doi.org/10.1371/journal.pone.0015198. PubMed: 21188145.
  14. 14. Conde A, Manuela Chaves M, Gerós H (2011) Membrane transport, sensing and signaling in plant adaptation to environmental stress. Plant Cell Physiol 52: 1583-1602. doi:https://doi.org/10.1093/pcp/pcr107. PubMed: 21828102.
  15. 15. Kim S-G, Kim S-Y, Park C-M (2007) A membrane-associated NAC transcription factor regulates salt-responsive flowering via FLOWERING LOCUS T in Arabidopsis. Planta 226: 647-654. doi:https://doi.org/10.1007/s00425-007-0513-3. PubMed: 17410378.
  16. 16. Vallejo AJ, Yanovsky MJ, Botto JF (2010) Germination variation in Arabidopsis thaliana accessions under moderate osmotic and salt stresses. Ann Bot 106: 833-842. doi:https://doi.org/10.1093/aob/mcq179. PubMed: 20861274.
  17. 17. DeRose-Wilson L, Gaut BS (2011) Mapping salinity tolerance during Arabidopsis thaliana germination and seedling growth. PLOS ONE 6: e22832. doi:https://doi.org/10.1371/journal.pone.0022832. PubMed: 21857956.
  18. 18. Munns R, Tester M (2008) Mechanisms of salinity tolerance. Annu Rev Plant Biol 59: 651-681. doi:https://doi.org/10.1146/annurev.arplant.59.032607.092911. PubMed: 18444910.
  19. 19. Galloway LF, Etterson JR (2007) Transgenerational plasticity is adaptive in the wild. Science 318: 1134-1136. doi:https://doi.org/10.1126/science.1148766. PubMed: 18006745.
  20. 20. Molinier J, Ries G, Zipfel C, Hohn B (2006) Transgeneration memory of stress in plants. Nature 442: 1046-1049. doi:https://doi.org/10.1038/nature05022. PubMed: 16892047.
  21. 21. Suter L, Widmer A (2013) Environmental heat and salt stress induce transgenerational phenotypic changes in Arabidopsis thaliana. PLOS ONE 8: e60364. doi:https://doi.org/10.1371/journal.pone.0060364. PubMed: 23585834.
  22. 22. Goldberg AD, Allis CD, Bernstein E (2007) Epigenetics: a landscape takes shape. Cell 128: 635-638. doi:https://doi.org/10.1016/j.cell.2007.02.006. PubMed: 17320500.
  23. 23. Saze H (2008) Epigenetic memory transmission through mitosis and meiosis in plants. Semin Cell Dev Biol 19: 527-536. doi:https://doi.org/10.1016/j.semcdb.2008.07.017. PubMed: 18708156.
  24. 24. Lauria M, Rossi V (2011) Epigenetic control of gene regulation in plants. Biochim Biophys Acta 1809: 369-378. doi:https://doi.org/10.1016/j.bbagrm.2011.03.002. PubMed: 21414429.
  25. 25. Verhoeven KJF, Jansen JJ, van Dijk PJ, Biere A (2010) Stress-induced DNA methylation changes and their heritability in asexual dandelions. New Phytol 185: 1108-1118. doi:https://doi.org/10.1111/j.1469-8137.2009.03121.x. PubMed: 20003072.
  26. 26. Boyko A, Kovalchuk I (2011) Genome instability and epigenetic modification - heritable responses to environmental stress? Curr Opin Plant Biol 14: 260-266. doi:https://doi.org/10.1016/j.pbi.2011.03.003. PubMed: 21440490.
  27. 27. Ito H, Gaubert H, Bucher E, Mirouze M, Vaillant I et al. (2011) An siRNA pathway prevents transgenerational retrotransposition in plants subjected to stress. Nature 472: 115-119. doi:https://doi.org/10.1038/nature09861. PubMed: 21399627.
  28. 28. Boyko A, Blevins T, Yao Y, Golubov A, Bilichak A et al. (2010) Transgenerational adaptation of Arabidopsis to stress requires DNA methylation and the function of Dicer-like proteins. PLOS ONE 5: e9514. doi:https://doi.org/10.1371/journal.pone.0009514. PubMed: 20209086.
  29. 29. Yi C, Zhang S, Liu X, Bui HTN, Hong Y (2010) Does epigenetic polymorphism contribute to phenotypic variances in Jatropha curcas. p. L.? BMC Plant Biol 10: 259.
  30. 30. Bossdorf O, Richards CL, Pigliucci M (2008) Epigenetics for ecologists. Ecol Lett 11: 106-115. PubMed: 18021243.
  31. 31. Rando OJ, Verstrepen KJ (2007) Timescales of genetic and epigenetic inheritance. Cell 128: 655-668. doi:https://doi.org/10.1016/j.cell.2007.01.023. PubMed: 17320504.
  32. 32. Paszkowski J, Grossniklaus U (2011) Selected aspects of transgenerational epigenetic inheritance and resetting in plants. Curr Opin Plant Biol 14: 195-203. doi:https://doi.org/10.1016/j.pbi.2011.01.002. PubMed: 21333585.
  33. 33. Simons AM (2011) Modes of response to environmental change and the elusive empirical evidence for bet hedging. Proc R Soc Lond B 278: 1601-1609. doi:https://doi.org/10.1098/rspb.2011.0176. PubMed: 21411456.
  34. 34. Reed TE, Waples RS, Schindler DE, Hard JJ, Kinnison MT (2010) Phenotypic plasticity and population viability: the importance of environmental predictability. Proc R Soc Lond B 277: 3391-3400. doi:https://doi.org/10.1098/rspb.2010.0771. PubMed: 20554553.
  35. 35. Nicotra AB, Atkin OK, Bonser SP, Davidson AM, Finnegan EJ et al. (2010) Plant phenotypic plasticity in a changing climate. Trends Plant Sci 15: 684-692. doi:https://doi.org/10.1016/j.tplants.2010.09.008. PubMed: 20970368.
  36. 36. Beaumont HJE, Gallie J, Kost C, Ferguson GC, Rainey PB (2009) Experimental evolution of bet hedging. Nature 462: 90-93. doi:https://doi.org/10.1038/nature08504. PubMed: 19890329.
  37. 37. Childs DZ, Metcalf CJE, Rees M (2010) Evolutionary bet-hedging in the real world: empirical evidence and challenges revealed by plants. Proc R Soc Lond B 277: 3055-3064. doi:https://doi.org/10.1098/rspb.2010.0707. PubMed: 20573624.
  38. 38. Larkindale J, Mishkind M, Vierling E (2005) Plant responses to high temperature. Plant Abiotic Stress. Blackwell Publishing Ltd.. pp. 100-144.
  39. 39. Tittel-Elmer M, Bucher E, Broger L, Mathieu O, Paszkowski J et al. (2010) Stress-induced activation of heterochromatic transcription. PLOS Genet 6: e1001175. PubMed: 21060865.
  40. 40. Lang-Mladek C, Popova O, Kiok K, Berlinger M, Rakic B et al. (2010) Transgenerational inheritance and resetting of stress-induced loss of epigenetic gene silencing in Arabidopsis. Mol Plants 3: 594-602. doi:https://doi.org/10.1093/mp/ssq014. PubMed: 20410255.
  41. 41. Reboud X, Corre VL, Scarcelli N, Roux F, David J et al. (2004) Natural variation among accessions of Arabidopsis thaliana: beyond the flowering date, what morphological traits are relevant to study adaptation? Plants Adaptation Molecular Genet Ecol: 135–142.
  42. 42. R Development Core Team (2010) R: a language and environment for statistical computing. Vienna, Austria: R Foundation for Statistical Computing.
  43. 43. Benjamini Y, Hochberg Y (1995) Controlling the false discovery rate - a practical and powerful approach to multiple testing. J R Stat Soc B Stat Methodol 57: 289-300.
  44. 44. Achard P, Cheng H, De Grauwe L, Decat J, Schoutteten H et al. (2006) Integration of plant responses to environmentally activated phytohormonal signals. Science 311: 91-94. doi:https://doi.org/10.1126/science.1118642. PubMed: 16400150.
  45. 45. Balasubramanian S, Sureshkumar S, Lempe J, Weigel D (2006) Potent induction of Arabidopsis thaliana flowering by elevated growth temperature. PLOS Genet 2: e106. doi:https://doi.org/10.1371/journal.pgen.0020106. PubMed: 16839183.
  46. 46. Balasubramanian S, Weigel D (2006) Temperature induced flowering in Arabidopsis thaliana. Plant Signal Behav 1: 227-228. doi:https://doi.org/10.4161/psb.1.5.3452. PubMed: 19704664.
  47. 47. Boyko A, Golubov A, Bilichak A, Kovalchuk I (2010) Chlorine Ions but not Sodium Ions Alter Genome Stability of Arabidopsis thaliana. Plant Cell Physiol 51: 1066-1078. doi:https://doi.org/10.1093/pcp/pcq048. PubMed: 20385609.
  48. 48. Wender NJ, Polisetty CR, Donohue K (2005) Density-dependent processes influencing the evolutionary dynamics of dispersal: a functional analysis of seed dispersal in Arabidopsis thaliana (Brassicaceae). Am J Bot 92: 960-971. doi:https://doi.org/10.3732/ajb.92.6.960. PubMed: 21652479.
  49. 49. Donohue K (2005) Niche construction through phenological plasticity: life history dynamics and ecological consequences. New Phytol 166: 83-92. doi:https://doi.org/10.1111/j.1469-8137.2005.01357.x. PubMed: 15760353.
  50. 50. Kurek I, Chang TK, Bertain SM, Madrigal A, Liu L et al. (2007) Enhanced thermostability of Arabidopsis Rubisco activase improves photosynthesis and growth rates under moderate heat stress. Plant Cell 19: 3230-3241. doi:https://doi.org/10.1105/tpc.107.054171. PubMed: 17933901.
  51. 51. Salvucci ME, Crafts-Brandner SJ (2004) Inhibition of photosynthesis by heat stress: the activation state of Rubisco as a limiting factor in photosynthesis. Physiol Plant 120: 179-186. doi:https://doi.org/10.1111/j.0031-9317.2004.0173.x. PubMed: 15032851.
  52. 52. Johannes F, Porcher E, Teixeira FK, Saliba-Colombani V, Simon M et al. (2009) Assessing the impact of transgenerational epigenetic variation on complex traits. PLOS Genet 5: e1000530. PubMed: 19557164.
  53. 53. Rainey PB, Beaumont HJ, Ferguson GC, Gallie J, Kost C et al. (2011) The evolutionary emergence of stochastic phenotype switching in bacteria. Microb Cell Factories 10: S14. doi:https://doi.org/10.1186/1475-2859-10-14. PubMed: 21995592.