Morpho-Physiological Characterization of Winter Wheat “Buster” Population during the Vegetative Stage under Heat Stress ()
1. Introduction
The global temperature increase is projected to range, depending on the location, between 1.5˚C and 11˚C by the year 2100 [1]. The high temperature at the beginning of spring season, coinciding with anthesis and grain-filling stages of the wheat crop, substantially reduces grain number and size [2]. This leads to a decrease in overall wheat productivity. Several studies have screened genotypes for heat tolerance [3] [4] [5]. However, most of these studies have considered either the whole plant life cycle or the post-anthesis period of crop growth. Hence, screening plants for the heat tolerance traits at early plant growth stages can help shorten the time period of the selection process, thus shortening breeding cycles and increasing genetic gain at a reduced time and cost. In this study, we screen 100 genotypes from a double haploid population for heat tolerance at the vegetative growth stage based on morpho-physiological traits.
Morpho-physiological improvements are one of the reasons behind increased productivity in winter wheat [6]. Selection of genotypes for higher yield based on their morphological characteristics including plant height in wheat [7], tiller number in wheat [8] and leaf area [9] in several horticultural crops such as cherries, grapes, and strawberries have been a successful approach for crops improvement. Similarly, morphological attributes are also taken into account while developing a tolerant variety for abiotic stresses such as heat and drought in wheat [10] [11] [12]. The number of effective tillers (fertile tillers) is an important yield attribute in wheat [13], and a source of difference between wheat genotypes [14]. Likewise, short plant height is an ideotype for wheat and one of the main reasons for the increase in wheat yields in the last five decades [6] [15]. Another important trait used for genotypic selection is the leaf area as the per unit area leaf traits may not represent the actual differences and produce misleading results [9] [16]. A negative correlation was recorded between leaf area and photosynthesis per unit leaf area as indicated by correlation analysis [16] [17]. Balota et al. [11] found that the drought-tolerant wheat varieties have significantly smaller leaves under both irrigated and drought conditions as compared to drought susceptible varieties. We take into account plant height, number of tillers, leaf area, and leaf number in order to make accurate inferences on the population studied.
In addition to the morphological traits, plant physiological processes have a deterministic effect on crop yield. Photosynthesis is one of the major factors influencing crop growth, biomass, and yield [18] [19] [20]. Plants are able to survive the climate extremes because of plasticity and resiliency of photosynthesis [21]. Therefore, understanding the response of photosynthesis to the changing environment is necessary to correctly assess the changes in plant productivity [22]. Xue et al. [23] found a positive correlation between leaf photosynthetic rates and grain yield in a few studies but no relation was observed in several studies identified in the review. According to Long et al. (2006) [24], leaf photosynthetic rates correlated poorly with yield in the past, but several recent studies are showing an increase in yield with an increase in photosynthetic rates.
In addition to exploring differences in photosynthetic rates between the genotypes, our study also pays attention to the stomatal conductance since it affects all gas exchange processes. Stomatal conductance is the rate of CO2 moving in and water vapor moving out of the stomatal apertures in leaf. The rates of diffusion of CO2 into leaf for photosynthesis and water vapor out of the leaf for transpiration are controlled by the stomatal aperture openings [25]. Variation in stomatal conductance among genotypes can be utilized in the selection of wheat lines for improved adaptation in a wide range of growing conditions [3]. Higher stomatal conductance is associated with an increase in grain yield in wheat [26]. The three processes that light can undergo in a leaf after the chlorophyll molecules receive light are photosynthesis, heat dissipation, and chlorophyll fluorescence. This chlorophyll fluorescence is the process of dissipating excess light as re-emission by chlorophyll A after fulfilling the photosynthetic demands [27]. These three processes always counterbalance each other’s efficiency increasing one of them while the others decrease [28]. Therefore, chlorophyll fluorescence ultimately reflects the photosynthetic activities of a plant in a complex manner [29]. Chlorophyll fluorescence measurement is one of the well-established techniques to evaluate the integrity of photosynthetic apparatus for stress detection in plants [30] [31].
During photosynthesis, chlorophyll absorbs photon for CO2 fixation [32]. If excess photons are absorbed by chlorophyll (more than a leaf can use for fixing CO2), then ROS are formed which cause photo-oxidative damage to the leaves [33] [34]. There are antioxidant compounds present in the leaves that scavenge the ROS and protect the photosynthetic apparatus (Photosystem-I + photosystem-II) [35] [36]. A category of such antioxidant compounds is phenolics. Phenolic compounds scavenge the ROS produced during light reactions in photosynthesis under moderate and high irradiance [32]. The concentration of phenolic compounds correlates positively to antioxidant activities [37]. [38] showed that phenolic compounds are reliable indicators for differences in genotypes in Triticale spp., especially in water deficit conditions where resistant genotypes had higher phenolic content compared to susceptible genotypes. Likewise, carotenoids are one of the indispensable components of photosynthetic mechanism in plants and several studies have demonstrated their importance [39] [40] [41]. Carotenoids play a major role in photosynthesis by harvesting light to extend the spectral range and protecting chlorophyll from photo-oxidative damage [39] [42] [43]. Most of the carotenoids are present in the thylakoid membrane of leaves, which is the site for light reactions of photosynthesis. Carotenoids improve electron transfer and light-harvesting efficiency of plants to stabilize the photosynthetic apparatus and protect it from photo-destruction [41].
This study aims to identify lines with desirable morpho-physiological traits for heat stress tolerance from a double haploid population called “Buster”. We hypothesize that double haploid lines of “Buster” population will enable us to amplify the responses to heat stress.
2. Materials and Methods
This study was conducted at Oklahoma State University in Stillwater, OK, USA. Two experiments were conducted, one in a greenhouse (Feb.-May, 2015) and in growth chambers. A subset of a hundred double haploid (DH) lines from a total of 256 Buster DH lines were used for both the experiments.
Buster Population
This study utilized the plant materials from a double haploid (DH) population developed by the OSU WIT. This population resulted from 32 F1’s obtained by crossing two popular wheat varieties “Duster” and “Billings”. From an ancestral perspective, “Duster” and “Billings” probably account for the largest segment of the elite germplasm currently flowing through the OSU WIT variety development program. These two parent lines demonstrate high yield potential with impressive disease resistance and end-use quality performance. However, they reach their yield in different and complementary ways with “Duster” having high kernel number and drought resistance, while “Billings” has large kernel size and is susceptible to drought. In addition, “Duster” and “Billings” show wide pattern differences in reproductive development, yet all known genes for reproductive development were identical between them. A population developed combing these varieties would have extremely high potential value to variety development. The OSU WIT envisions that a DH population would lead to trait discoveries, marker discoveries, knowledge of inheritance, and reduce the breeding time that would have far-reaching impact in further manipulating the pipeline (B. Carver, personal communication).
To this effect, 36 F1 seeds from the single cross Duster/Billings (OK10x994) were provided to Heartland Plant Innovations (HPI, Manhattan, KS) on 10/26/10, with the expectation to produce 300 haploids (DHs). Colchicine treatment was used to develop the DHs. At HPI, the D0 and D1 plant generations were reared and D2 seed was provided to WIT at OSU in 2012. A total of 278 DHs were on sufficient supply to plant back in non-replicated single-row observation plots in 2012-2013 at Stillwater. About 271 DHs were then advanced for further evaluation in 2013-2014. The 271 lines were arbitrarily assigned to 6 sets of 42 lines each, plus one overflow set of 19 DHs, to evaluate in replicated field plots in 2014, 2015, and 2016 at Stillwater. Sets were created to reduce block size in the field, and the two parents were included in each set as a common check. Seed yield and seed of 256 DH lines including the parental varieties were available from the 2013-2014 season. The Billings/Duster DH population will be referred to as Buster population (B. Carver, personal communication).
The 100 Buster lines from the total of 256 were selected from a yield trial conducted during the 2013-2014 growing season in Stillwater, OK. This was an extreme drought year with a 270% yield difference between the low and high yielding lines. Each of the 6 sub-group (described earlier) was divided into high, average and low yield based on the mean yield ±1 standard deviation. From each yield group, five lines were selected resulting in 15 lines for each subgroup. A few additional lines with extreme yield values along with parents were selected to create the set of 100 Buster lines for this research.
Greenhouse conditions
The greenhouse study was conducted without imposing any artificial stress. Five seeds of each line were sown in PVC pots 50 cm deep and 15 cm in diameter. There were two replications i.e. two pots per genotype. Pure sand was used as a rooting medium to obtain optimum control of water and nutrient supply to roots. An automatic drip irrigation system was used to supply 0.3 L of Hoagland’s nutrient solution to the plants each time, four times a day at 8:00 AM, 12:00 PM, 4:00 PM and 8:00 PM. In this study, data on leaf morphology (length, width, and area), plant developmental changes (plant height, tiller number, and leaf number) and pigment concentrations (chlorophyll A, chlorophyll B, carotenoids and phenolic content) were collected.
Photosynthetic pigments
Chlorophyll A, chlorophyll B, and carotenoids were extracted using dimethyl sulfoxide (DMSO). Five leaf discs, 1 cm2 each, were punched from five randomly selected leaves from each pot. The leaf discs were stored in 5 ml of DMSO for 24 hours in the dark. The concentrations of the pigments were calculated from absorbance values obtained with a spectrophotometer (Genesys 10 Bio Spectrophotometer, Thermo Scientific) at 664 nm, 648 nm and 470 nm for chlorophyll A, chlorophyll B, and carotenoids respectively using equations by [44]:
,
,
where,
A = absorbance at respective wavelengths;
c = pigment concentration (µg/mL of extract).
To determine phenolic compounds concentration, another set of five-leaf discs, 1 cm2 each was punched from the same leaves. The leaf discs were placed in the extractant solution for 24 hours at room temperature. The solution used for the extraction of phenolic compounds was composed of methanol, water, and hydrochloric acid in the ratio of 79:20:1.
Absorbance values were obtained at 330 nm for phenolic compounds and the concentration was calculated as given by [45]:
where, C = concentration of phenolic compounds (µg/mL of extract), A = absorbance at 330 nm.
Leaf area
Three leaves were randomly selected per pot for leaf morphological data measurements after 75 days after sowing (DAS). LI-3000 (Licor Inc., NE, USA) portable leaf area meter was used for measuring leaf area, leaf length, maximum width and the average width of each selected leaf. Leaves were carefully selected from the same position on three different plants to avoid differences in physiological age.
Growth attributes
Two plants per pot were randomly selected and marked during their seedling stage. Data on tiller number, leaf number, and plant height were collected from the marked plants on a weekly basis for five weeks starting when leaf nodes were visible in the sampled plants. Plant heights were recorded from the base of the plant to the uppermost collar on the main stem. The rate of increase in plant height per day was calculated. Tiller number was counted for each of the two plants and leaf numbers were counted in the main stem of the same plants.
Growth Chamber experiment
The growth chamber experiment was conducted in a controlled environment research laboratory (CERL) at Oklahoma State University. Four identical growth chambers each with 50 PVC pots (15 cm in diameter and 35 cm in depth) were used for the study. The set of 100 Buster lines were thus split between two growth chambers. Fifty Buster lines were planted in each chamber, four seeds of one Buster line in each pot. An automated irrigation system was used to provide 0.3 L of Hoagland’s nutrient solution to the plants each time, three times a day (8:00 AM, 1:00 PM and 6:00 PM), after germination. Pure sand was used as the growth medium. Plants in all chambers were grown at temperatures (22/16˚C day/night) up to 65 days after sowing (DAS) and continued to grow at the same temperature in two chambers which were designated as controls. The temperature was raised to 32/26˚C (day/night) in two of the chambers after 65 DAS to impose heat stress on one set of Buster lines. A gradual increase in temperature from night to day and vice-versa was achieved through temperature ramping. The photoperiod was adjusted to 14 hours light period and 10 hours dark period. Thus, one set of the 100 lines was under heat stress treatment, and the other set was grown under control conditions. Data on plant physiological parameters were collected from this experiment.
Gas exchange parameters and fluorescence
Parameters such as Pn, gs, E, Ci, ETR and Fv'/Fm' were determined on leaves between 9 AM to 1 PM using an infrared gas analyzer (IRGA) in an open photosynthesis system, LI-6400 XT (Licor Inc., NE, USA). The two youngest fully opened leaves from adjacent plants (similar leaf position) were used for the measurements in order to cover the 2 cm2 area of the leaf cuvette. The leaves were artificially irradiated at 1200 µmol·m−2·s−1 with a blue-red LED radiation source attached to the sensor head set. The temperature in the leaf cuvette was set in accordance with the daytime temperature of the treatment chambers. The leaf chamber reference CO2 was set to 400 µL·L−1.
The efficiency of energy harvesting by photosystem II (PSII) was calculated by built-in algorithms in the LI-6400XT system using the equation:
where,
Fo' = minimal fluorescence of a momentarily darkened leaf;
Fm' = maximum fluorescence during a saturating flash light;
Fv' = variable fluorescence during a saturating flash light.
Instantaneous water use efficiency (IWUE) was calculated as the ratio of net photosynthesis (Pn) to transpiration (T).
The gas exchange parameters and fluorescence were measured three times: first before starting the heat stress treatment, second after three days of the treatment, and third a week after the treatment. The measurements taken after the introduction of heat stress are expressed as an average and compared to the average before heat stress treatment.
Statistical analysis
Data collected was analyzed using SAS Version 9.4 (SAS Institute, Cary, NC). The Analysis of Variance (ANOVA) was performed using PROC GLM to see if the Buster lines are statistically significant at p = 0.05 probability level for the recorded parameters. PROC CORR was used to obtain correlation coefficients between the different parameters. Graphs were constructed using Sigma Plot.
Principal component analysis (PCA) was conducted using PROC PRINCOMP on the gas-exchange parameters data. The PCA was performed on the differences between values of the parameters in control and treatment conditions to identify the variables that were mainly causing the differences. A biplot, graphical representation of eigenvectors, also known as loadings, of the first two PC scores, was constructed using PROC PRINQUAL.
3. Results
Greenhouse study
Growth attributes (plant height, tiller number, and leaf number)
The average plant height ranged from 8.4 cm in line DH97 to 13.3 cm in DH14 and DH29. Significant differences (P < 0.05) were observed between Buster lines for the rate of increase in plant height and tiller number. The highest rate of increase in plant height was found in a parental line “Billings”, and the lowest rate was observed in a Buster line “DH231” followed by another parental line “Duster”. The number of tillers and rate of increase in tiller number were highest in Buster line “DH136” and lowest in “DH224”.
The number of leaves on the main stem was not significantly different among the Buster lines. The final plant height is weakly negatively correlated to the rate of increase in leaf number (Table 1).
Photosynthetic pigments
Differences in concentrations were observed for different pigments among the Buster lines but were not statistically significant (P > 0.05). The p-values for chlorophyll A, chlorophyll B, carotenoids and phenolic contents are 0.75, 0.95, 0.11 and 0.29, respectively. Chlorophyll A:B ratio ranged from 1.8 in DH55 to 4.3 in DH54. The photosynthetic pigments positively correlated to each other at 0.01 levels of significance (Table 1). Among those, chlorophyll A and B, and carotenoids showed a stronger correlation with each other.
To depict the observed differences, the Buster lines are grouped for individual pigment components chlorophyll A, chlorophyll B, carotenoids and phenolic content based on mean ± 1 and 2 standard deviations (Figure 1).
Leaf morphological attributes
Significant (P < 0.05) differences between Buster lines were observed in leaf area, leaf length and leaf width. The leaf area ranged from 20.94 cm2 in parental line “Duster” to 38.56 cm2 in Buster line “DH73” with an average of 28.85 cm2 and s.d. of 3.59 cm2 in the population. The Buster lines are grouped based on leaf area mean ± 1 s.d. and mean ± 2 s.d. (Figure 2). Leaf area was significantly positively correlated to the height attributes whereas it negatively correlated to leaf number attributes (Table 1).
Figure 2. Number of Buster lines for each group of leaf area based on mean ± s.d.
Growth Chamber study
Gas exchange parameters and leaf fluorescence
Before heat stress, the Buster lines were not statistically different for gas exchange and fluorescence parameters (data not shown). Whereas after heat stress, significant differences (P < 0.05) were observed between the genotypes in interaction with temperature for Pn, gs, T, Fv'/Fm', Ci and IWUE with p-values of 0.0041, 0.0095, 0.0011, 0.0166, 0.0002, and less than 0.0001 respectively (Table 2).
The Pn, gs, T and Ci increased in the Buster lines in response to increased temperature. The parameters Fv'/Fm' and ETR did not show consistent responses to the increase in temperature for the studied Buster lines. In response to the increase in temperature, the values for Fv'/Fm' decreased in 58 lines and increased in 37 lines, whereas ETR increased in 78 lines and decreased in 17 lines.
Upon studying the relationship between the recorded parameters, the strongest correlation is observed between gs and T with a correlation coefficient of 0.97. The variables Pn, gs, E and Ci all were strongly positively correlated (correlation coefficients greater than 0.8) with each other (Table 3). The IWUE was negatively correlated to Pn, gs, T, and Ci. The decrease in IWUE in response to heat stress ranged from 1.28 µmol CO2/µmol H2O (DH263) to 8.45 µmol CO2/µmol H2O (DH102). The decrease was 2.6 and 3.4 µmol CO2/µmol H2O for the parental lines “Duster” and “Billings” respectively. The electron transport rate was not correlated to any of the parameters and Fv'/Fm' had a weak positive correlation with Pn and Ci.
From the results of PCA, it was observed that more than 80% of the variability was explained by the first two PC scores (Table 4). The biplot shows that gs and IWUE explain most of the variability in the first and second axis, respectively (Figure 3). The Buster lines in the upper right (first) quadrant have a comparatively small increase in gs, T and Pn and the Buster lines in lower (third and fourth) quadrants are relatively less affected by heat stress as indicated by less difference in IWUE between control and treatment conditions.
Table 2. P-values for photosynthesis, stomatal conductance, transpiration, ETR, fluorescence, instantaneous WUE and intercellular CO2 showing significant differences for main factors (genotype and heat stress treatment) and their interaction.
* Significant at α = 0.05, ** Significant at α = 0.01, *** Significant at α = 0.001, NS not significant.
Table 3. Matrix of Pearson correlation coefficients showing correlation between photosynthesis (Pn), stomatal conductance (gs), transpiration (E), intercellular CO2 (Ci), ETR, fluorescence (Fv'/Fm') and instantaneous water use efficiency (IWUE).
*Significant at α = 0.05 **Significant at α = 0.01 ***Significant at α = 0.001, NSnot significant.
Table 4. Eigenvectors (loadings) of the principal components and proportional and cumulative variance explained by the principal components.
Figure 3. Biplot of the eigenvectors of first two principal component scores. The genotype numbers correspond to genotype numbers in Table 4.
4. Discussion
Greenhouse study
Growth attributes (plant height, tiller number, and leaf number)
The wide range in average plant height within the population indicates that the population is well segregated for plant height. “Duster” and “Billings” are both categorized as intermediate semi-dwarfs but differ in plant heights; “Duster” (71 cm) [46] and “Billings” (73 cm) [47]. Although these plant heights represent the heights from ground level to spike tip, these differences were visible during the vegetative plant growth stages in our study.
A higher number of tillers contribute towards higher harvest index in normal conditions but a reduced number of tillers are desirable under water deficit [8]. The differences in tiller numbers, in this case, can be explained by probable segregation of genes in the population for tiller number because the parental line “Duster” is known to have a high tiller number while “Billings” lacks this attribute. Duster and Billings showed wide pattern differences in reproductive development, yet all known genes for reproductive development were identical between them (B. Carver, personal communication). The number of tillers observed is an indicator of the tillering capacity of a genotype because the plants were putting new tillers for a long time since there was no vernalization imposed for the plants to start the reproductive phase.
The parental lines “Duster” and “Billings” do not have reported differences for the leaf number or rates of increase in leaf number. This might be a result of a limited time period in which the data was collected i.e. five weeks. Unlike tiller number and plant height, which gain measurable increments in a short time period, it requires a long time for a leaf to be fully developed and thus the data collection duration might not have been sufficient enough to reflect the differences in leaf number among the genotypes if any.
Photosynthetic pigments
A study done by [48] found no significant differences for chlorophyll content among genotypes when the first spikelet of the inflorescence was visible, close to the vegetative stage, but differences were recorded for later growth stages [48]. Studies have shown differences in chlorophyll content in different wheat cultivars and under different stress conditions [49] [50]. There is very limited information on carotenoids analysis on leaves, especially in the context of wheat where studies are concentrated towards grain carotenoids content in durum wheat. Studies show differences in carotenoids content in leaves for different genotypes and stress combinations in crops like soybeans [51] and tomatoes [52]. Similarly, phenolic content in leaves and their role in scavenging ROS have been studied under various abiotic stress conditions like drought and salinity but not yet studied for non-stressed conditions. Differences in phenolic compounds concentration have been found in response to stress at various growth stages where differences were mostly expressed at reproductive stages of plant growth [53]. The heterogeneity of the Buster lines is deemed to be not enough to reflect differences in the photosynthetic pigments to a significant effect.
Leaf morphological attributes
Most of the Buster lines demonstrated greater leaf area than the parental lines because of the segregation of genes in the DH population. In a study by [54], significant differences were observed in leaf areas and a weak negative correlation between leaf area and water use efficiency was reported.
Correlation between parameters
No such correlations have been studied specifically, but studies have been done with the plant height and leaf area as the selection criteria. Our results showing the plant height and leaf area positively correlated to each other provides an opportunity to select Buster lines with shorter plant height and lower leaf area at the same time, which are the desirable characteristics of wheat in drought-prone areas.
Growth Chamber study
Gas exchange parameters and leaf fluorescence
In the past, genotypic differences for gas exchange parameters among different wheat varieties are commonly observed [23] [55] [56]. Our results indicate that the differences in Buster lines are prominently expressed under stress than in normal conditions.
Since all gas exchange parameters are directly related to stomatal opening, their increase under heat stress is directly related to increased stomatal apertures. The stomatal conductance is not limited by the high temperature unless water stress is associated with it [57]. The higher stomatal conductance under high temperatures also explains the transpirational cooling mechanism of plants in response to high temperatures. Higher transpiration rates in higher temperatures allow more water vapor to exit the leaves ultimately having a cooling effect. On the other hand, increased stomatal openings allow more CO2 to enter the leaves, which increases photosynthesis. Furthermore, an increase in enzymatic activity of Rubisco with an increase in temperature also influences photosynthesis in high temperatures and sufficient water conditions [22]. A review done by [58] suggests that the yields of cotton and wheat are directly correlated to the stomatal conductance under supra optimal temperatures when there’s no influence of other stresses like drought and vapor pressure deficit. They also concluded that an increase in stomatal conductance is an avoiding type of resistance in response to the high temperatures but the water use efficiency is decreased with the increase in temperature because of wasteful water use, which is in accordance with the results of this experiment. However, in the literature, the gas exchange parameters are found to have no correlation with the yield parameters [59]. This study gives us an idea of the potential performance of Buster lines but cannot conclude on the plant responses to naturally occurring heat stress that is associated with water stress most of the time.
The correlation of gs with Pn and T is expected because the rate of flow of CO2 and water vapor in and out of the leaves is controlled by the stomatal aperture. However, the negative correlation of IWUE with the gas exchange parameters suggests that water is not being efficiently used and photosynthesis is increased at a very high cost of water.
A study done by [23] under drought reported no correlation between the gas exchange parameters and IWUE.
5. Conclusion
The Buster lines were significantly different in the morphological traits such as plant height, tiller number, and leaf area under heat stress. The tiller number recorded in this study is an indicator of the potential tillering capacity of the genotype. The genotypic differences based on gas exchange parameters were discernible only under heat stress. The variables gs and IWUE explained most of the differences between the treatments. The IWUE decreased in response to heat stress in all Buster lines at different rates, whereas the values for gas exchange parameters increased under heat stress compared to controlled conditions. Thus, the Buster lines can potentially be selected for further breeding programs based on their morphological differences and physiological response to stress.
Acknowledgements
Authors would like to acknowledge support from Oklahoma Agricultural Experiment Station, Hatch project no. 2969 and Oklahoma Wheat Commission for this study.
Abbreviations
DH, Double Haploid; Pn, net photosynthesis; gs, stomatal conductance; Ci, internal CO2 concentration; T, transpiration; IWUE, instantaneous water use efficiency; ETR, electron transport rate; Fv'/Fm', Photosystem II efficiency; IRGA, infrared gas analyzer; CERL, controlled environment research laboratory; µL·L−1, microliter per liter; µmol·m−2·s−1, micro mole per meter squared per second; ROS, reactive oxygen species; OSUWIT, Oklahoma State University Wheat Improvement Team; DMSO, Dimethyl Sulfoxide.