- Research Article
- Open Access
High altitude population of Arabidopsis thaliana is more plastic and adaptive under common garden than controlled condition
BMC Ecology volume 17, Article number: 39 (2017)
Population differentiation and their adaptation to a particular environment depend on their ability to respond to a new environment. This, in turn is governed to an extent, by the degree of phenotypic plasticity exhibited by the populations. The populations of same species inhabiting different climatic conditions may differ in their phenotypic plasticity. Himalayan populations of Arabidopsis thaliana originating from a steep altitude are exposed to different climatic conditions ranging from sub-tropical to temperate. Thus they might have experienced different selection pressures during evolution and may respond differently under common environmental condition.
Phenotypic plasticity and differentiation of natural populations of A. thaliana grown under common garden and controlled conditions were determined. A total of seventeen morphological traits, their plasticity, association between traits and environment were performed using 45 accessions from three populations. Plants from different altitudes differed in phenotypes, their selection and fitness under two conditions. Under both the conditions lower altitude population was characterized by higher leaf count and larger silique than higher and middle altitude population. Flowering time of high altitude population increased while that of low and medium altitude decreased under controlled condition compared to open field. An increase in seed weight and germination was observed for all the population under open field than controlled. Rosette area was under divergent selection in both the condition. The heritability of lower altitude population was the highest under both the conditions, where as it was the least for higher altitude population further indicating that the high altitude populations are more responsive towards phenotypic changes under new environmental conditions. Ninety-nine percent of variability in traits and their plasticity co-varied with the altitude of their origin. The population of high altitude was more plastic and differentiated as compared to the lower altitude one.
Arabidopsis thaliana population native to different altitudes of the west Himalaya responds differently when grown under common environments. The success of high altitude population is more in common garden than the controlled conditions. The significant variability in phenotype and its association with altitude of origin predicts for non-random genetic differentiation among the populations.
Plants being sessile organisms are frequently exposed to heterogeneous environmental conditions. The variation in environmental condition influences the evolution of traits and their differentiation. It can favour either local adaptation (adaptive differentiation/specialization) or evolution of plastic genotypes that respond through phenotypic plasticity (generalization). Both adaptive differentiation and phenotypic plasticity acts as major factor in generation of inter-individual variation and plant diversity. The amount of variation both within and among populations influences the response towards the changing environment. There has been a long standing interest to understand the role of plasticity on the performance of individuals, populations and species in new environment . Variations in the environmental conditions can lead to both phenotypic and genotypic variations among the populations over a long period of time. Natural selection acting on phenotypes and their plasticity can cause evolution of populations leading to local adaptation.
Arabidopsis thaliana (L.) Heynh., has recently been studied as a model plant species for the ecological and eco-genomics studies [2, 3]. This is mainly because of its ease of growth and ubiquitous presence as well as publicly available large information resources. Being a model organism, it is important to quantify the diversity of this plant contributing to adaptation throughout its geographical distribution. There are reports of environmental influence on phenotypic traits like number of fruits, germination, length and width of fruits, flowering time, flooding response, etc. using field grown Arabidopsis plants [4,5,6,7]. Similarly, a few studies investigated the fitness of Arabidopsis mutants grown under field conditions and commented on morphological changes that are likely to result from varied environmental influence [3, 4, 8]. There are also reports on phenotypic, physiological and gene expression profiling of plants grown either in simulated or controlled condition [9,10,11,12]. Phylogeography of A. thaliana in western Eurasia [13, 14] and its rapid expansion throughout the world has also been studied . Though, a lot of studies have been conducted throughout the world including other European populations [6, 11, 13, 14, 16,17,18,19], no such detailed studies have been conducted on highly diverse populations from West Himalaya. This incomplete sampling could limit our knowledge on extent of variation known in A. thaliana.
The West Himalayan populations are inhabited along wide altitude ranging from ~ 700 to ~ 3400 m amsl, which is also the altitudinal maxima of A. thaliana distribution reported so far. In contrast to predominantly temperate climatic condition in other European and American countries, A. thaliana growing along Himalayan altitude are exposed to sub-tropical to temperate climatic conditions . Thus historically, these populations might have experienced different evolutionary pressure than the other world populations. However, these populations were never described earlier in terms of their phenotypic response towards different environmental conditions. These populations provide an opportunity for studying the effect of different environmental factors and their combinations on phenotypic diversity and differentiation. The strong variation in environmental factors along the steep altitude [17, 21] can impose strong selection on traits leading to their adaptation . The highest altitude population (~ 3400 m amsl) which is exposed to extreme weather conditions might have evolved differently as compared to the more benign lower altitude population. Thus the fitness of this population may be higher under a common but heterogeneous environment as compared to lower altitude population. These populations are relatively old and harbour significant amount of genetic variations between the populations [23, 24]. Interestingly, these populations were also found to be genetically quite different from the rest of the world populations .
Here, we measured the phenotypic differentiation and plasticity of natural populations of A. thaliana from West Himalaya. A few studies had shown and argued that in order to understand population differentiation a relative analysis of the field and controlled grown plants is essential [5, 25, 26]. A significant variation in fitness has also been observed in mutants grown in the two conditions . We compared the level of plasticity of these populations exhibited due to interaction between genotypes and environment. Specifically, we asked the following questions: Do the populations from different altitudes differ in the expression of traits in the common garden and controlled conditions (i.e. plasticity)? Are these differentiations population specific in either conditions, and related to altitude? How the populations differ in the selection of the traits in the two conditions? Under which of the two novel conditions plants show more fitness? The variations in phenotype could be either due to adaptive differentiation or random by genetic drift. An association of phenotype under common growth condition to its native altitude and climate predicts for an adaptive genetic differentiation of these populations reducing chance of being a random processes . These differences are also expected to indicate the abiotic stress experienced by the populations in the field conditions.
Trait variation among populations in either condition
The three populations grown in common garden (CG) and controlled condition (GH) varied significantly for all the vegetative and reproductive traits, except a few (Table 1). Besides being significantly variable, some of the traits also followed a trend of either increase or decrease with altitude. For example, in CG the lower altitude population (Deh) was characterized by lower leaf area as compared to middle (Mun) and higher (Chit) altitude population. Whereas, in the GH Deh was characterized by wider leaves, longer petioles and larger leaf area as compared to other two populations. Leaf counts and silique length were higher in Deh as compared to Mun and Chit in both the conditions. Chit flowered the most late, followed by Mun and the Deh in both the conditions. The post hoc Bonferroni pair wise comparison between the populations indicated that the lower (Deh) and the highest altitude populations (Chit) were most differentiated under both the conditions. However, the middle altitude population (Mun) showed a variable trend, being more similar to the Deh and Chit under GH and CG condition, respectively (Additional file 1: Table S1).
Trait variations between the two conditions
The differences in phenotypic trait expression between the two conditions varied at different level (Fig. 1; Additional file 1: Table S2). All the vegetative traits showed a significant decrease in GH as compared to CG for all the populations, except leaf length:width ratio and rosette major:minor axis ratio. The leaf length:width ratio (leaf shape) did not varied significantly for Deh and Mun in the two conditions. However, the leaf length:width ratio increased for Chit in GH. Thus the leaf was more rounded in Chit under CG as compared to GH condition. An increase in rosette major:minor axis ratio (rosette shape) was observed for all the populations in GH as compared to CG. Flowering time, expressed as days to bolting increased in Chit while that of Deh and Mun decreased in GH as compared to CG.
Further, among the reproductive traits, there was a decrease in plant height, inflorescence length and silique length in CG as compared to GH in all the populations. The length of pedicel of silique increased in CG as compared to GH in all the populations. A decrease in distance between two siliques was observed in lower altitude population (Deh) while it increased for Mun and Chit in CG as compared to GH. The total number of flowers produced decreased for Deh and Chit while it increased for Mun in CG. The number of fruits produced also differed significantly among the populations grown in the two conditions. Chit and Mun produced more number of fruits in CG as compared to GH where as Deh produced less fruits in CG as compared to GH. An increase in seed weight was observed for all the population in CG [population (P), P = 0.261; environment (E), P = 0.028; P × E, P = 0.049]. Germination percentage was significantly higher in CG produced seeds than GH [population (P), P = 0.233; environment (E), P < 0.0001; P × E, P = 0.183] in all populations. Overall, in either condition the populations had sufficient amount of phenotypic variations to distinguish the three populations as indicated by the discriminant analysis (Fig. 2a).
Correlation between traits within and between two conditions
The expression of traits was generally similar in the two conditions except for leaf length and rosette area (Table 2). The leaf area was found to be negatively correlated between the two conditions (r = − 0.335, P = 0.025). However the correlation between the phenotypic traits within the two conditions differed both in the strength and direction (Additional file 1: Table S3). Among the vegetative traits, leaf number was negatively correlated with leaf area in CG, while it was positively correlated with leaf width, leaf area, petiole length, rosette area and negatively with leaf shape. Leaf shape showed a negative correlation with leaf area and positive with leaf length and width in CG, while in GH it was positively correlated with leaf length and negatively with leaf number, leaf width and rosette shape. Rosette shape was not correlated with any of the vegetative traits in CG, while it was negatively correlated with leaf length and leaf shape in the GH. Similar variation in both strength and direction of correlation was also observed for the reproductive traits (Additional file 1: Table S3). Most of the reproductive traits were positively correlated in the CG and negatively correlated in GH. The fruits produced in CG were positively correlated with most of the traits except (leaf number; r = − 0.375, P = 0.011) while in GH, it was correlated only to the number of flowers (r = 0.812, P < 0.001).
Correlation of both vegetative and reproductive traits with flowering time within the two conditions also differed both in strength and direction (Fig. 3). In the CG, lower leaf area, rounded rosette, elongated leaf larger siliques with smaller pedicels, lesser distance between two siliques, lesser plant height and inflorescence length were correlated with early flowering (Fig. 3a; Additional file 1: Table S3). While in the GH, flowering time was negatively correlated with leaf number, area and width, petiole length, silique length and fertility percentage and positively with leaf length, leaf shape, plant height, inflorescence length and distance between two siliques (Fig. 3b; Additional file 1: Table S3).
Phenotypic plasticity of the populations
The P × E interaction term was significant for all the traits indicating for the significant variation in the phenotypic plasticity. In addition to using P × E interaction as a measurement of phenotypic plasticity, we also quantified the amount of plasticity. Significance of variation in the amount of phenotypic plasticity of each trait and population was tested by ANOVA (Additional file 1: Table S4). Interestingly, we found that the variation in phenotypic plasticity between the extreme populations (Chit and Deh) were significant for all the traits. This differentiation decreased as the difference in altitude between the populations decreased. This is because there was no significant variation in the amount of plasticity for leaf shape, number of fruits, number of flowers and flowering time between Deh (lower altitude) and Mun (middle altitude). Although there was variations in slope of reaction norms for the number of fruits and flowers between Deh and Mun but there was no significant variation in their strength i.e., amount of plasticity. The plasticity of all the traits between Mun and Chit (middle and higher altitude) was significant, except plant height (Additional file 1: Table S5). Further most of the reproductive traits of Chit were less plastic as compared to Deh and Mun. On the other hand Mun exhibited greater plasticity for most of the reproductive traits as compared to Deh and Chit (Additional file 1: Table S6). A scatter plot representing the measure of plasticity as residual variation from the best fit line also indicated the same (Additional file 2: Figure S1). The overall plasticity of Chit [0.291 ± 0.038 (mean ± SD)] was significantly higher than Mun (0.262 ± 0.012) followed by Deh (0.181 ± 0.014). The overall phenotypic plasticity was positively correlated with the flowering time in both the conditions (Fig. 3). A clear differentiation of the populations on the basis of the expression of the plasticity of the traits was observed (Fig. 2b).
Reproductive fitness and path analysis for fitness in the two growth conditions
The expression of the reproductive component of fitness differences between the populations depended on the growth conditions as indicated by significant ‘P × E’ interaction term (F2,84 = 9.301; P < 0.0001). Pair- wise comparison of the same population grown under two conditions revealed that the expression of reproductive component of fitness in Deh decreased significantly whereas it increased in Mun and Chit in CG (though not significant). Further, no significant differences in reproductive component of fitness between the populations was found in both the conditions after adjusting for phenotypes and their plasticity (CG − F2,10 = 0.802, P = 0.475; GH − F2,10 = 1.569, P = 0.255).
A partial least square path analysis was performed to determine the traits those affect the reproductive component of fitness. The traits contributing to the fitness varied both in strength and direction in the two conditions (Fig. 4; Additional file 1: Table S7). In the CG, the vegetative traits showed a direct positive effect on flowering time and reproductive traits. In contrast, vegetative traits were negatively affecting both these traits in the GH. Reproductive fitness was affected positively by vegetative traits in the CG and negatively in GH. On the other hand, the flowering time negatively affected fitness in the CG. The reproductive traits contributed positively to the fitness in the two conditions but with different strengths. The loading of the individual traits on the vegetative, flowering and reproductive components also varied (Additional file 1: Table S8). Among vegetative traits, petiole length and leaf area were significantly loaded on vegetative component in the CG. In contrast, leaf number, leaf width and leaf area was significantly loaded in the GH. In reproductive component, plant height, length of inflorescence, distance between two siliques, and number of flowers were loaded in CG. While in the GH, plant height, length of inflorescence, size of silique, distance between two siliques and fitness percentage was significantly loaded on reproductive component.
Broad sense heritability (H2) and Qst–Fst comparison
The values of heritability were generally higher for the GH grown plants than the CG [GH = 0.151025 ± 0.040659 (mean ± SE); CG = 0.033834 ± 0.012189 (mean ± SE)]. The mean heritability was generally lower for vegetative traits than reproductive traits in both the conditions (Table 3). The lowest mean heritability was shown by petiole length (mean ± SE = 0.000779 ± 0.000779) under CG condition. In contrast, the highest mean heritability was observed for flowering time under GH condition (0.645 ± 0.204093). At population level, Deh showed the highest mean heritability (0.052575 ± 0.020291) followed by Mun (0.034178 ± 0.02102) and Chit (0.014749 ± 0.008268) under CG condition. In the GH condition Mun (0.071429 ± 0.02444) showed the lowest and Deh (0.213917 ± 0.073223) showed the highest mean heritability.
Previously, we reported the genetic differentiation as FST for these populations . The estimated FST for the three populations was 0.543 (95% CI = 0.44–0.63). Overall the total variance of traits explained as genetic component (QST) was higher in the CG than GH [(mean ± SE) CG = 0.6411 ± 0.04; GH = 0.447 ± 0.075]. The QST estimated ranged from 0 to 1 for GH grown plants and from 0.37 to 1 for CG grown plants representing the level of differential quantitative differentiation (Fig. 5). QST was found to be significantly higher than FST for leaf length, rosette area, rosette shape and fertility percentage under the CG condition. Under GH condition, QST of only rosette area was higher than FST representing complete differentiation and divergent selection. Few traits like rosette shape, pedicel of silique, number of flowers and fruits did not varied at population level in the GH condition. For the rest of the traits, the confidence interval of QST and FST overlapped, suggesting for neutral differentiation.
Effect of altitude of origin of populations in the two conditions
The phenotypic variation observed in both the growth conditions was significantly associated with the altitude of origin of the populations (CG − Wilk’s lambda = 0.0001; F = 139.9; P = < 0.0001; GH − Wilk’s lambda = 0.002; F = 34.9; P = < 0.0001). 99% of the total multivariate trait variation co-varied with the altitude as described by the canonical correlation analysis for both the sites (CG-R 2c = 0.99, Wilk’s statistic = 0.005; GH-R 2c = 0.99, Wilk’s statistic = 0.020). The variation in the expression of the phenotypic plasticity by the three populations was also significantly associated with the altitude of origin (Wilk’s lambda = 0.0001; F = 712.7; P = < 0.0001; R 2c = 0.99, Wilk’s statistic = 0.002).
Phenotypic trait expression is the effect of interaction between genotype and environment. Populations that grow under different climatic conditions are good source for studying these interactions and their effect on phenotype. In this study, we determined the phenotypic differentiation of natural populations of A. thaliana originating from different altitudes. These Indian populations are genetically diverse varying significantly from other world populations and are characterised by strong population structure . The populations showed a clear differentiation on the basis of both absolute phenotype and their plasticity as demonstrated by the discriminant analysis. The existing genetic variations in these populations might be responsible for observed phenotypic variation of traits. The phenotypic variation can be achieved either in the form of adaptive differentiation (for traits that show similar trends in the phenotypic expression in the two conditions) or random phenotypic variation (for traits that show different trends in phenotype in the two conditions). For example, in Deh higher number of leaves and longer siliques were observed in both the conditions. In our earlier study, involving the same but native populations also showed that Deh had higher leaf numbers and longer siliques . This indicates that these traits are fixed. On the other hand, leaf and rosette shape were variable (varying non-significantly in CG and significantly in GH) when compared between CG and GH. Again, these two traits were not variable in their native sites . Differences in the level of significance among the populations and the environmental conditions indicate that these characters are merely random in response and lack any adaptive differentiation. These populations, thus show both an adaptive differentiation and random response depending on environmental condition.
Among the vegetative traits, leaf size and numbers are important adaptive traits. These traits in turns are correlated with other associated rosette and leaf traits. Larger and greater number of leaves was observed in the CG than the GH. Size variation have a direct effect on whole-plant growth rate, mainly through changes in the conductance of the boundary layer, affecting exchange of heat, uptake of carbon dioxide and transpiration . This is mainly due to the fact that larger leaves having thicker boundary layer of slower convective heat loss tends to be hotter than ambient air temperature . In CG there was significant difference in day and night temperature. The large leaf in CG might help to maintain hotter temperature than the ambient air temperature prevalent during night . However, this might be a limitation during water deficit conditions, but we provided sufficient water throughout the growth period. Thus, though edaphic factors and day temperature in field (~ 22–23 °C) were similar with that of GH, there was significant difference in the light intensity and day/night temperature difference between the two conditions, responsible for above observation.
Besides vegetative growth, flowering time is an important indicator of plant fitness. We observed significant differences in the flowering time between the two conditions as well as among the populations in a particular condition. Flowering time of Deh and Mun increased in the CG condition. This was in accordance to the earlier study that showed an increase in flowering time in the field conditions . In contrast, we observed a decrease in flowering time in the field conditions as compared to GH for higher altitude population (Chit). This may be due to the fact that while delayed bolting under benign conditions can lead to higher fecundity , early flowering is advantageous in habitats where mortality is likely to be early [31,32,33,34]. The winter (Dec–Jan) and spring (Feb–March) seasons in Lucknow, India is followed by hot and dry summer, where temperature rises beyond 40 °C which is unfavourable for growth of A. thaliana as against more suitable benign condition in GH. Thus, the CG grown Chit plants completed their life cycle earlier before the onset of unfavourable condition. On the other hand, the climatic factors in field condition of earlier study  was more favourable where the maximum temperature recorded during experiments was 28 °C. Further, this trend of flowering time indicates population specific response to climatic condition in context of flowering time. The high altitude population (Chit) might flower early, sensing the unfavourable condition in the CG whereas in GH it has favourable condition to prolong the life cycle. These findings indicate population specific dependency in flowering time under different conditions. This variation in flowering time among the populations is commonly reported in A. thaliana [16, 35, 36].
The reproductive traits (measured at the end of life cycle), represents the overall productivity, reproductive allocation and development of the plants. Among other reproductive traits, plant height and length of inflorescence are important traits. These traits were also correlated either positively or negatively with the most of the other traits. A number of studies describing shade avoidance syndrome have also shown an increase in plant height under shade conditions (e.g. [25, 37, 38]). It has also been suggested that low light intensity favours the internodal elongation and hence an increase in plant height . The increase in plant height observed in GH grown plants might be due to the low light in GH as compared to the CG. Seed sizes and germination ability, a measure of maternal investment also determine the plants adaptability under a particular environment [40, 41]. Though there was no significant differences in these traits between any pair of populations in either conditions but when compared between the two environments, seed size (weight) increased in all the populations in CG. More particularly, the seed size differences were more in higher altitude populations (higher plasticity). This may be due to the higher adaptability of high altitude plants in more heterogeneous climatic conditions that is prevalent at high altitude as compared to low altitude plants. High altitude plants may produce more of these reproductive traits for better performance under such heterogeneous conditions.
The observed significance of variation in phenotypic plasticity of the populations implies an underlying genetic polymorphism commonly found in A. thaliana [42, 43]. However, the greater overall plasticity observed in the highest altitude plant (Chit) indicates that being exposed to more heterogeneous and stressful environment might have selected individuals having greater plasticity. The greater plasticity in the high altitude population in turn helps them to respond to more stressful and heterogeneous condition of CG than GH. The Deh population experiencing less heterogeneous and stressful conditions suffered a comparative loss of reproductive fitness in the CG. This was in accordance with the earlier study where gene ontology term enrichment of SNP containing genes for the lowest altitude population (Deh) did not retrieved any abiotic stress related term, while high altitude population (Mun and Chit) showed enrichment of GO terms related to high altitude stresses . Plasticity of the traits increases the ecological tolerance [44, 45] which results in higher fitness. Thus the differences in the reproductive fitness observed between the populations were nullified when all the measured traits and their plasticity were included. Selection of more plastic genotypes for environmental heterogeneity has also been reported in other studies [46, 47].
Further, partial least square method of path analysis showed that the different traits were differentially regulating reproductive fitness in the two conditions. The flowering time has been shown to significantly affect fitness of plants [4, 16, 35, 48]. But this association was observed in the CG only. The adverse effect of increase in flowering time on the fitness was in accordance with earlier studies [48,49,50]. The vegetative traits were shown to affect reproductive fitness through its correlated effects on flowering time and reproductive traits.
The broad sense heritability of the traits differed between the two growth conditions. This observation for heritability of traits in two conditions suggests that heritability was a function of both genotype and environment. Overall the total variance of traits explained as genetic component (QST) was higher in the CG than GH. This also suggests that the expression of genetically driven phenotypes depends on the environmental condition. The QST–FST comparison of these traits helps in predicting the role of neutral or selective divergence in shaping phenotypic variation. Although the method of estimation has been criticized for its biasness due to its inherent theoretical caveats, but still it is commonly used to identify the adaptive traits [51, 52]. Based on this approach we found that the different traits were under divergent selection in the two conditions. The only trait that was under selection in both the condition was rosette area, predicting it to be relevant under two conditions. Interestingly, the confidence interval of QST and FST of flowering time overlapped suggesting its neutral divergence in these populations. In other studies this trait has been shown to be under local divergent selection both in the field and controlled conditions [16, 53]. The significance of variation of this trait in the two conditions and neutrality of its selection suggests for genetic drift effect. However, one needs to be careful in interpreting the data as accession used for estimating QST and FST were different but from the same populations.
Further, there was a significant multivariate association of all the phenotypic traits and their plasticity with the altitude of their native sites. This association of native altitude with the traits that were measured in the uncorrelated environmental conditions (CG and GH) predicts for an altitudinal divergence commonly found by other studies [17, 35, 37, 54, 55]. The similar association of altitude and climate with phenotypic and genetic data were observed earlier using same set of populations [23, 24, 27]. However, our analysis is limited due to use of low number of populations and thus caution may be taken to interpret this particular result.
West Himalayan populations of A. thaliana showed significantly variable morphological response when grown in two different conditions. The traits that could be important for fitness under field condition may not play significant role under controlled conditions. These populations also showed differential genotype–environment interaction for all the traits studied. Further, the high altitude population was more plastic and had higher reproductive fitness in CG as compared to the low altitude populations. The high altitude populations of Himalaya are exposed to harsh environmental condition and thus might have evolved to fit well in the heterogeneous condition (like CG), hence showing maximum plasticity and reproductive fitness as compared to the lower altitude populations. However, experimental validation needs to be performed using reciprocal transplant experiment to prove the observed phenomenon. We also emphasise that caution may be taken to interpret these results due to small population size. However, more efforts are being made to phenotype more populations from this region.
Study population and experimental design
The details of the populations and climatic condition of the native sites have been reported earlier . Seeds were collected from 15 accessions from each of the three sites [viz., Deh (700 m above msl); Mun (1800 m above msl); and Chit (3400 m above msl)] maintaining a minimum distance of 2 m between the two accessions. Seeds were grown for one generation under controlled conditions in order to eliminate effect due to level of seed maturity in field conditions. Stratified seeds were grown under two environmental conditions: (1) controlled conditions (GH) with 16 h/8 h light/dark; temperature 20–22 °C, and light intensity—120 PPFD; (2) common garden (CG) at Council of Scientific and Industrial Research-National Botanical Research Institute (CSIR-NBRI), Lucknow, India (details in Additional file 1: Table S9, Additional file 2: Figure S2). To check the best environmental condition for germination of seeds and growth of the plants in CG, batches of seed were sown at weekly interval starting from mid November to mid December in the previous year. The appropriate environmental condition for the growth of A. thaliana in CG was found to be from the last week of November to March. Subsequently, experiment was set up during last week of November, 2013 in randomized-block factorial design consisting of two factors: population (P) and growth condition (E) viz. CG and GH. 180 pots (4 replicate × 15 accessions × 3 populations) filled with Soilrite Mix® (Keltech Energies Ltd. Perlite Division, India.) were planted with five seeds per pot for each environmental condition. After germination, only one healthy seedling per pot was maintained. The plants were watered regularly. The accessions were randomised within the trays per block to avoid positional effects. The CG, GH and native sites of the populations is shown in Fig. 6.
Phenotypic trait and plasticity measurement
A total of 17 vegetative, flowering and reproductive traits were measured during the experiment. Vegetative traits were measured after 5 weeks of germination while reproductive traits were measured after completion of life cycle (Additional file 1: Tables S10, S11). Leaf area was estimated as product of length and width, while leaf shape was represented as a ratio of length to width. Rosette area was measured as product of major and minor axis and rosette shape as a ratio of major to minor axis. We considered both the leaf and rosette shape as rounded when the ratio is close to one and elongated if the ratio was more than one. The plants were censused in every two days for the estimation of flowering time until all plants have bolted. Flowering time was estimated as number of days from germination to bolting. After the completion of life cycle, percentage of total number of flowers produced that converted into fruits was used to estimate fertility percentage. Being highly correlated with seed production, total reproductive fitness was measured as the number of fruits produced [16, 56]. The seed weight was estimated by measuring the weight of 100 seeds. This measurement was taken randomly for six accessions from each population and for each condition (6 accessions × 3 population × 2 environment). For each condition, equal number of seeds from each of the 15 accessions per population was pooled in a replicate of four (3 population × 2 environment × 4 replicate) for the estimation of germination percentage after 5 days.
Amount of phenotypic plasticity was calculated for each trait as the change in mean trait value in one environment to the other . Least square means of each accession was used for calculating phenotypic plasticity. Overall plasticity of each accession was estimated as the average of plasticity of each trait. A bootstrapping procedure, with 1000 iterations was used to obtain confidence interval and error estimate for the plasticity, using bootMer function of lme4 package in R . The error estimate for overall plasticity was calculated as average of the errors obtained for all traits. A scatter plot of raw trait values between the two conditions was made to visualize the phenotypic plasticity as residual deviation from the best fit and expected line (observed–expected).
Statistical data analysis
All the statistical analysis (except mixed model ANOVA) was performed with statistical software IBM-SPSS 22 (IBM Corp., released 2013, IBM SPSS Statistics for windows, Version 22.0. IBM Corp, Armonk, NY). Leaf shape was square root transformed to meet the assumption of normal distribution. All the analysis was performed on least square mean of replicates estimated for each accession per population. In our previous study, we have shown that these populations harbour considerable amount of genetic variation within population . Thus, each accession was considered as replicate within population for statistical analysis. ANOVA was performed to find the trait differentiations between the populations in the two conditions separately. Post hoc bonferroni test was used to compare the population pair wise.
A mixed model ANOVA was conducted to identify the differences in the expression of phenotypic traits with the interaction of the environment by nlme package of R . ‘Population (P)’, ‘environment (E)’ and ‘population × environment (P × E)’ interaction were used as fixed factors and accessions as random factor. Replicates were considered as blocks. Discriminant analysis was performed to aid the visualization of population differentiation. The significance of variation in the amount of plasticity was also analysed using one way ANOVA and visualized via discriminant analysis.
Correlation among the populations for different traits within and between the two conditions was analysed using Pearson’s correlation. A correlation of all the phenotypic traits with the flowering time was also calculated to predict whether the difference in the flowering time was related to the phenotypic traits and their plasticity.
To test the differences in performance of the three populations in the two sites, ANOVA, for relative and absolute reproductive fitness was conducted with population (P) and growth condition (E) as fixed effects and accessions as random effect. A significant ‘P × E’ term for relative reproductive fitness was used as indication of fitness differences present between the populations that depends on the environment. To find out the extent to which environment dependent selection on plastic characters explains these fitness differences between the populations, ANCOVA was conducted on relative reproductive fitness as dependent variable and all the measured traits and their plasticity as covariate. Further in order to study potential interrelationship and effect of vegetative, flowering and reproductive traits on the reproductive fitness of the plants in the two conditions, partial least square regression path analysis was performed. We utilized partial least square method of path analysis. This method is different from conventional covariance-based method as it do not impose any assumption on data distribution . The traits were divided into vegetative, flowering and reproductive variables to predict fitness as number of fruits produced. All these traits are supposed to influence fitness of plants. The hypothesized model for path analysis is shows in (Additional file 2: Figure S3). Loading of traits greater than 0.7 on their respective variables were considered to be significant. Confidence intervals of the path coefficients were calculated using 1000 bootstraps. This analysis was performed in plspm package of R .
Broad sense heritability was estimated for each population and each trait as H2 = VG/(VG + VE), where VG is among accession variance and VE is residual variance [16, 53]. The amount of quantitative genetic variation, QST  was estimated for all the traits. QST was estimated assuming complete selfing as VB/(VB + VW), where VB is among population variation and VW is within population variation [16, 53, 63]. 95% confidence intervals were calculated for both QST and H2 using restricted maximum likelihood (REML) variance components . FST used as a measure of genetic differentiation was calculated from previously published microsatellite data of these populations . The comparison of QST and FST is commonly used method to differentiate selection from neutral divergence for the traits. QST > FST suggests for divergent selection, while QST < FST predicts for uniform selection. The traits that show QST = FST are under neutral divergence [65, 66]. Though, this method (QST–FST comparison) of identifying selection has been predicted to be biased but is still used to identify traits under selection . The QST was considered to differ significantly from FST only when their confidence intervals did not overlap.
The significance of impact of altitude of native sites of the populations on the expression of traits and their plasticity was assessed using MANOVA. Canonical correlation was also performed to account for percentage of variance in the phenotypic traits explained (using squared canonical correlation coefficient) by altitude of the site of origin. These analyses were performed individually for CG, GH and phenotypic plasticity.
common garden set up CSIR-National Botanical Research Institute, Lucknow
controlled environmental conditions used for plant growth
growth condition (environment factor)
- P × E:
population and growth condition interaction term (factor)
- H2 :
broad sense heritability
- QST :
quantitative population differentiation
- FST :
Forsman A. Rethinking phenotypic plasticity and its consequences for individuals, populations and species. Heredity. 2015;115(4):276–84.
Mitchell-Olds T. Arabidopsis thaliana and its wild relatives: a model system for ecology and evolution. Trends Ecol Evol. 2001;16(12):693–700.
Weigel D. Natural variation in Arabidopsis: from molecular genetics to ecological genomics. Plant Physiol. 2012;158(1):2–22.
Malmberg RL, Held S, Waits A, Mauricio R. Epistasis for fitness-related quantitative traits in Arabidopsis thaliana grown in the field and in the greenhouse. Genetics. 2005;171(4):2013–27.
Mishra Y, Jankanpaa HJ, Kiss AZ, Funk C, Schroder WP, Jansson S. Arabidopsis plants grown in the field and climate chambers significantly differ in leaf morphology and photosystem components. BMC Plant Biol. 2012;12:6.
Moharekar SL, Moharekar S, Kobayashi T, Ishii H, Sumida A, Hara T. Phenotypic plasticity and ecotypic variations in growth and flowering time of Arabidopsis thaliana (L.) under different light and temperature conditions. Indian J Exp Biol. 2014;52(4):344–51.
Pigliucci M, Kolodynska A. Phenotypic plasticity and integration in response to flooded conditions in natural accessions of Arabidopsis thaliana (L.) Heynh (Brassicaceae). Ann Bot. 2002;90(2):199–207.
Pigliucci M. Ecology and evolutionary biology of Arabidopsis. In: The Arabidopsis book. vol. 1. Rockville: American Society of Plant Biologists; 2002.
Gong Q, Li P, Ma S, Indu Rupassara S, Bohnert HJ. Salinity stress adaptation competence in the extremophile Thellungiella halophila in comparison with its relative Arabidopsis thaliana. Plant J. 2005;44(5):826–39.
Kazachkova Y, Batushansky A, Cisneros A, Tel-Zur N, Fait A, Barak S. Growth platform-dependent and-independent phenotypic and metabolic responses of Arabidopsis and its halophytic relative, Eutrema salsugineum, to salt stress. Plant Physiol. 2013;162(3):1583–98.
Li Y, Cheng R, Spokas KA, Palmer AA, Borevitz JO. Genetic variation for life history sensitivity to seasonal warming in Arabidopsis thaliana. Genetics. 2014;196(2):569–77.
Provart NJ, Gil P, Chen W, Han B, Chang H-S, Wang X, Zhu T. Gene expression phenotypes of Arabidopsis associated with sensitivity to low temperatures. Plant Physiol. 2003;132(2):893–906.
François O, Blum MG, Jakobsson M, Rosenberg NA. Demographic history of European populations of Arabidopsis thaliana. PLoS Genet. 2008;4(5):e1000075.
Sharbel TF, Haubold B, Mitchell-Olds T. Genetic isolation by distance in Arabidopsis thaliana: biogeography and postglacial colonization of Europe. Mol Ecol. 2000;9(12):2109–18.
Provan J, Campanella JJ. Patterns of cytoplasmic variation in Arabidopsis thaliana (Brassicaceae) revealed by polymorphic chloroplast microsatellites. Syst Bot. 2003;28(3):578–83.
Méndez-Vigo B, Gomaa NH, Alonso-Blanco C, Xavier Pico F. Among-and within-population variation in flowering time of Iberian Arabidopsis thaliana estimated in field and glasshouse conditions. New Phytol. 2013;197(4):1332–43.
Montesinos-Navarro A, Wig J, Xavier Pico F, Tonsor SJ. Arabidopsis thaliana populations show clinal variation in a climatic gradient associated with altitude. New Phytol. 2011;189(1):282–94.
Rutter MT, Fenster CB. Testing for adaptation to climate in Arabidopsis thaliana: a calibrated common garden approach. Ann Bot. 2007;99(3):529–36.
Wolfe MD, Tonsor SJ. Adaptation to spring heat and drought in northeastern Spanish Arabidopsis thaliana. New Phytol. 2014;201(1):323–34.
Rawat D. Elevational reduction of plant species diversity in high altitudes of Garhwal Himalaya, India. Curr Sci. 2011;100(6):833–6.
Körner C. The use of ‘altitude’ in ecological research. Trends Ecol Evol. 2007;22(11):569–74.
Endler JA. Natural selection in the wild. Princeton: Princeton University Press; 1986.
Tyagi A, Singh S, Mishra P, Singh A, Tripathi AM, Jena SN, Roy S. Genetic diversity and population structure of Arabidopsis thaliana along an altitudinal gradient. AoB Plants. 2016;8:plv145.
Tyagi A, Yadav A, Tripathi AM, Roy S. High light intensity plays a major role in emergence of population level variation in Arabidopsis thaliana along an altitudinal gradient. Sci Rep. 2016;6:26160.
Donohue K, Messiqua D, Pyle EH, Heschel MS, Schmitt J. Evidence of adaptive divergence in plasticity: density-and site-dependent selection on shade-avoidance responses in Impatiens capensis. Evolution. 2000;54(6):1956–68.
Kulheim C, Agren J, Jansson S. Rapid regulation of light harvesting and plant fitness in the field. Science. 2002;297(5578):91–3.
Singh A, Tyagi A, Tripathi AM, Gokhale SM, Singh N, Roy S. Morphological trait variations in the west Himalayan (India) populations of Arabidopsis thaliana along altitudinal gradients. Curr Sci. 2015;108(12):2213.
Givnish TJ. Comparative studies of leaf form: assessing the relative roles of selective pressures and phylogenetic constraints. New Phytol. 1987;106(s1):131–60.
Malhado A, Malhi Y, Whittaker R, Ladle R, Ter Steege H, Phillips O, Butt N, Aragão L, Quesada C, Araujo-Murakami A. Spatial trends in leaf size of Amazonian rainforest trees. Biogeosciences. 2009;6(8):1563–76.
Tonsor SJ, Scheiner SM. Plastic trait integration across a CO2 gradient in Arabidopsis thaliana. Am Nat. 2007;169(5):E119–40.
Alonso-Blanco C, Méndez-Vigo B. Genetic architecture of naturally occurring quantitative traits in plants: an updated synthesis. Curr Opin Plant Biol. 2014;18:37–43.
Biere A. Genotypic and plastic variation in plant size: effects on fecundity and allocation patterns in Lychnis flos-cuculi along a gradient of natural soil fertility. J Ecol. 1995;83(4):629–42.
Fox GA. Drought and the evolution of flowering time in desert annuals. Am J Bot. 1990:1508–18.
Lacey EP. Onset of reproduction in plants: size-versus age-dependency. Trends Ecol Evol. 1986;1(3):72–5.
Méndez-Vigo B, Picó FX, Ramiro M, Martínez-Zapater JM, Alonso-Blanco C. Altitudinal and climatic adaptation is mediated by flowering traits and FRI, FLC, and PHYC genes in Arabidopsis. Plant Physiol. 2011;157(4):1942–55.
Stinchcombe JR, Weinig C, Ungerer M, Olsen KM, Mays C, Halldorsdottir SS, Purugganan MD, Schmitt J. A latitudinal cline in flowering time in Arabidopsis thaliana modulated by the flowering time gene FRIGIDA. Proc Natl Acad Sci USA. 2004;101(13):4712–7.
Botto JF. Plasticity to simulated shade is associated with altitude in structured populations of Arabidopsis thaliana. Plant Cell Environ. 2015;38(7):1321–32.
Schmitt J, Stinchcombe JR, Heschel MS, Huber H. The adaptive evolution of plasticity: phytochrome-mediated shade avoidance responses. Integr Comp Biol. 2003;43(3):459–69.
Callahan HS, Pigliucci M. Shade-induced plasticity and its ecological significance in wild populations of Arabidopsis thaliana. Ecology. 2002;83(7):1965–80.
Mousseau TA, Fox CW. The adaptive significance of maternal effects. Trends Ecol Evol. 1998;13(10):403–7.
Wulff RD, Cáceres A, Schmitt J. Seed and seedling responses to maternal and offspring environments in Plantago lanceolata. Funct Ecol. 1994;8(6):763–9.
Pigliucci M. Phenotypic plasticity: beyond nature and nurture. Baltimore: JHU Press; 2001.
West-Eberhard MJ. Developmental plasticity and evolution. Oxford: Oxford University Press; 2003.
Andersson S, Widén B. Reaction norm variation in a rare plant, Senecio integrifolius (Asteraceae). Heredity. 1994;73(6):598–607.
Bradshaw AD. Evolutionary significance of phenotypic plasticity in plants. Adv Genet. 1965;13(1):115–55.
Lázaro-Nogal A, Matesanz S, Godoy A, Pérez-Trautman F, Gianoli E, Valladares F. Environmental heterogeneity leads to higher plasticity in dry-edge populations of a semi-arid Chilean shrub: insights into climate change responses. J Ecol. 2015;103(2):338–50.
Scheiner SM. The genetics of phenotypic plasticity. XII. Temporal and spatial heterogeneity. Ecol Evol. 2013;3(13):4596–609.
Springate DA, Kover PX. Plant responses to elevated temperatures: a field study on phenological sensitivity and fitness responses to simulated climate warming. Glob Change Biol. 2014;20(2):456–65.
Amasino R. Seasonal and developmental timing of flowering. Plant J. 2010;61(6):1001–13.
Kover PX, Rowntree J, Scarcelli N, Savriama Y, Eldridge T, Schaal B. Pleiotropic effects of environment-specific adaptation in Arabidopsis thaliana. New Phytol. 2009;183(3):816–25.
Leinonen T, McCairns RS, O’hara RB, Merilä J. QST–FST comparisons: evolutionary and ecological insights from genomic heterogeneity. Nat Rev Genet. 2013;14(3):179.
O’Hara RB, Merilä J. Bias and precision in Q ST estimates: problems and some solutions. Genetics. 2005;171(3):1331–9.
Corre V. Variation at two flowering time genes within and among populations of Arabidopsis thaliana: comparison with markers and traits. Mol Ecol. 2005;14(13):4181–92.
Ågren J, Schemske DW. Reciprocal transplants demonstrate strong adaptive differentiation of the model organism Arabidopsis thaliana in its native range. New Phytol. 2012;194(4):1112–22.
Molina-Montenegro MA, Naya DE. Latitudinal patterns in phenotypic plasticity and fitness-related traits: assessing the climatic variability hypothesis (CVH) with an invasive plant species. PLoS One. 2012;7(10):e47620.
Pigliucci M, Schmitt J. Genes affecting phenotypic plasticity in Arabidopsis: pleiotropic effects and reproductive fitness of photomorphogenic mutants. J Evol Biol. 1999;12(3):551–62.
Valladares F, Sanchez-Gomez D, Zavala MA. Quantitative estimation of phenotypic plasticity: bridging the gap between the evolutionary concept and its ecological applications. J Ecol. 2006;94(6):1103–16.
Bates D, Maechler M, Bolker B, Walker S. lme4: Linear mixed-effects models using Eigen and S4. R Package Version 1.7; 2014. p. 1–23.
Pinheiro J, Bates D, DebRoy S, Sarkar D, R Core Team. nlme: linear and nonlinear mixed effects models. R package version 3.1-117. 2014. Available at http://CRAN.r-project.org/package=nlme. Accessed 15 Oct 2015.
Sanchez G. PLS path modeling with R. Berkeley: Trowchez Editions; 2013.
Sanchez G, Trinchera L, Russolillo G. plspm: tools for partial least squares path modeling (PLS-PM). R Package Vers. 2013;04:1.
Spitze K. Population structure in Daphnia obtusa: quantitative genetic and allozymic variation. Genetics. 1993;135(2):367–74.
Bonnin I, Prosperi J-M, Olivieri I. Genetic markers and quantitative genetic variation in Medicago truncatula (Leguminosae): a comparative analysis of population structure. Genetics. 1996;143(4):1795–805.
Lynch M, Walsh B. Genetics and analysis of quantitative traits, vol. 1: Sunderland: Sinauer; 1998.
Merilä J, Crnokrak P. Comparison of genetic differentiation at marker loci and quantitative traits. J Evol Biol. 2001;14(6):892–903.
Whitlock MC. Evolutionary inference from QST. Mol Ecol. 2008;17(8):1885–96.
SR conceived and designed the experiments; AS conducted the experiments; AS and SR analysed the results; SR and AS wrote the manuscript. Both authors read and approved the final manuscript.
The authors are thankful to the Director, CSIR-NBRI, Lucknow for providing facilities and support. The authors also acknowledge Council for Scientific and Industrial Research, New Delhi for fellowship provided to AS.
The authors declare that they have no competing interests.
Availability of data and materials
All data generated and analysed during the study are included in this published article and its additional files.
Consent for publication
Ethics approval and consent to participate
The project was funded by Council of Scientific and Industrial Research, New Delhi under project: SIMPLE (BSC-109). The funding body played no role in design of the study and collection, analysis, and interpretation of the data or writing of the manuscript.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Additional file 1: Table S1. Post-hoc (Bonferroni) test (P values) for pair wise comparison of traits between the three populations of Arabidopsis thaliana in common garden (CG) and controlled condition (GH). Significant figures are in bold. Table S2. Analysis of significance of effect of population and growth environment (CG vs. GH) on vegetative, flowering and reproductive traits for three populations of Arabidopsis thaliana using mixed model ANOVA. Table S3. Correlation between morphological traits measured within the two conditions. The upper triangle shows Pearson correlation (r) between traits within the common garden (CG) and lower triangle for within the controlled condition (GH). Significant values are shown in bold (*P < 0.05; **P < 0.01; ***P < 0.001). Table S4. Analysis of variation for amount of phenotypic plasticity. The significance of variation in phenotypic plasticity among the three populations was tested using ANOVA. All traits were significant at P value < 0.05. Table S5. Post-hoc (Bonferroni) test (P values) for pair wise comparison of plasticity for all the traits between the three populations of Arabidopsis thaliana. Table S6. Phenotypic plasticity of the three populations of Arabidopsis thaliana. The confidence interval of plasticity was calculated using 1000 iterations. Table S7. Coefficients for direct, indirect and bootstrap estimated total effect for the latent variables used for path analysis in the two condition. Total effect is sum of direct and indirect effect. Numbers in brackets represent 95% confidence interval for total effect. Table S8. Loading and weight of all the traits on their latent variable (blocks) estimated using partial least square path analysis. 95% confidence intervals estimated using 1000 bootstrap are shown in brackets. Table S9. Average monthly weather conditions of Common garden site at Lucknow, India during the growth period. Altitude: 113 m above mean sea level; Geographical coordinates: 260 55′ N, 800 59′ E. Table S10. Ls mean of phenotypic traits for all accessions with 95% confidence interval. Table S11. Raw data of all the phenotypic traits measured for all accessions under the two growth conditions.
Additional file 2: Figure S1. Scatter plot of accession means in CG vs GH. The observed line is shown as dotted red line and best fit line by solid black colour. Figure S2. Log climatic conditions of Common garden site at Lucknow, India during the growth period. (a) Daily Average Temperature; (b) Daily Maximum Temperature; (C) Daily Minimum Temperature; (D) Daily Light Intensity. Altitude: 113 m above mean sea level; Geographical coordinates: 260 55′ N, 800 59′ E. The weather data of CG was obtained from Amausi weather station, Lucknow situated at around 16 km from CG Freely available from- http://en.tutiempo.net/climate/ws-423690.html, Accessed 5th October 2015. Figure S3. Model used for path analysis. Oval circles represent the latent variables with their measured variables in the square boxes.
About this article
Cite this article
Singh, A., Roy, S. High altitude population of Arabidopsis thaliana is more plastic and adaptive under common garden than controlled condition. BMC Ecol 17, 39 (2017). https://doi.org/10.1186/s12898-017-0149-5
- West Himalaya
- Arabidopsis thaliana
- Phenotypic differentiation
- Phenotypic plasticity
- Common garden
- Controlled condition