Jiemeng Xu1
Claudia Lowe1
Sergio G. Hernandez-Leon2
Susanne Dreisigacker3
Matthew P. Reynolds3
Elisa M. Valenzuela-Soto2
Matthew J. Paul1
Sigrid Heuer1,4*
- 1Plant Science Department, Rothamsted Research, Harpenden, United Kingdom
- 2Centro de Investigación en Alimentación y Desarrollo A.C., Carretera Gustavo Enrique Aztiazarán Rosas, Hermosillo, Mexico
- 3International Maize and Wheat Improvement Center (CIMMYT), Texcoco, Mexico
- 4Pre-Breeding Department, National Institute of Agricultual Botany (NIAB), Cambridge, United Kingdom
Rising temperatures due to climate change threaten agricultural crop productivity. As a cool-season crop, wheat is heat-sensitive, but often exposed to high temperatures during the cultivation period. In the current study, a bread wheat panel of spring wheat genotypes, including putatively heat-tolerant Australian and CIMMYT genotypes, was exposed to a 5-day mild (34°C/28°C, day/night) or extreme (37°C/27°C) heat stress during the sensitive pollen developmental stage. Worsening effects on anther morphology were observed, as heat stress increased from mild to extreme. Even under mild heat, a significant decrease in pollen viability and number of grains per spike from primary spike was observed compared with the control (21°C/15°C), with Sunstar and two CIMMYT breeding lines performing well. A heat-specific positive correlation between the two traits indicates the important role of pollen fertility for grain setting. Interestingly, both mild and extreme heat induced development of new tillers after the heat stress, providing an alternative sink for accumulated photosynthates and significantly contributing to the final yield. Measurements of flag leaf maximum potential quantum efficiency of photosystem II (Fv/Fm) showed an initial inhibition after the heat treatment, followed by a full recovery within a few days. Despite this, model fitting using chlorophyll soil plant analysis development (SPAD) measurements showed an earlier onset or faster senescence rate under heat stress. The data presented here provide interesting entry points for further research into pollen fertility, tillering dynamics, and leaf senescence under heat. The identified heat-tolerant wheat genotypes can be used to dissect the underlying mechanisms and breed climate-resilient wheat.
Introduction
Wheat is one of the most important crops for human consumption, grown on 220 million hectares with a total production of 760 million tons in 2020 (FAOSTAT). In scenarios of climate change, wheat plants are prone to be exposed to warmer and more variable temperatures (Trnka et al., 2014). Beyond a physiological threshold, high temperatures cause stress and impair plant growth and development. Both historical data and future predictions have revealed the negative effects of heat on wheat productivity at the global and regional scale (Liu et al., 2016; Zampieri et al., 2017; Pequeno et al., 2021). Therefore, it is crucial to identify and breed heat-adapted varieties to sustain wheat production and ensure food security.
In nature, the adverse effects of heat stress on plants can be variable depending on the intensity, duration, and developmental stage (Stone and Nicolas, 1995; Yeh et al., 2012). Most of the heat-related studies in wheat have been field-based and used late sowing to expose plants to high temperatures during the flowering and grain filling stages; however, short episodes of heat during earlier reproductive stages can also cause significant damage (Zampieri et al., 2017). Indeed, anther and pollen development are considered to be the stages most vulnerable to heat stress (Zinn et al., 2010; Rieu et al., 2017). Stage-specific treatments have found that wheat is particularly sensitive to heat around 8 days before anthesis, which coincides with the early meiosis to tetrad stage of pollen development (Saini and Aspinall, 1982; Prasad and Djanaguiraman, 2014). Because pollen development occurs during the booting stage, while spikes are still inside the developing pseudostem in wheat, the length of the auricles between the flag leaf and the penultimate leaf (referred as auricle interval length, AIL) has been used as a proxy for pollen development. AIL between 3 and 6 cm has been associated with this sensitive stage (Bokshi et al., 2021; Erena et al., 2021). Brief heat exposure during this sensitive period resulted in abnormal meiosis behavior (Omidi et al., 2014; Draeger and Moore, 2017) and a significant reduction in pollen fertility (Prasad and Djanaguiraman, 2014; Begcy et al., 2018; Browne et al., 2021). A few studies have examined the natural variation in pollen viability under heat stress and its association with yield, as booting usually occurs during the cooler time of the cropping season and it is difficult to apply precise stage-specific heat stress (Bheemanahalli et al., 2019; Bokshi et al., 2021). However, considering a warmer and increasingly erratic climate, this area warrants further investigation.
Spike number is one of the main components in determining wheat yield; it is highly variable and responsive to the environmental factors (Slafer et al., 2014). Interestingly, contrasting responses of spike number under heat stress have been reported. When exposed to continuous high temperatures during the terminal flowering and grain filling stages, spike formation and tillering were always reduced (Cai et al., 2016; Sharma et al., 2016; Dwivedi et al., 2017; Kumar et al., 2021). In contrast, after a short episode of heat stress during earlier developmental stages, spike numbers increased (Bányai et al., 2014; Chavan et al., 2019; Hütsch et al., 2019). Enhanced spike formation after early heat stress is surprising, but the underlying tillering dynamics and impact on final yield have not been discussed.
In addition to the effects on pollen fertility and spike number, heat-induced yield loss has also been ascribed to accelerated leaf senescence, shortening the duration of grain filling (Cossani and Reynolds, 2012; Pinto et al., 2016; Shirdelmoghanloo et al., 2016; Bergkamp et al., 2018; Sade et al., 2018). As indicators of senescence, chlorophyll soil plant analysis development (SPAD) (chlorophyll content index) (Richardson et al., 2002) and Fv/Fm (the maximum potential quantum yield of photosystem II) (Murchie and Lawson, 2013) have been widely used to evaluate this trait. Under terminal heat, SPAD and Fv/Fm were often reduced in leaf tissue during senescence and were closely related to yield-contributing traits, such as thousand grain weight (Talukder et al., 2014; Hassan et al., 2018; Mirosavljević et al., 2021; Touzy et al., 2022). Studies for the genetic analysis of these leaf senescence related traits are also available (Azam et al., 2015; Bhusal et al., 2018; Touzy et al., 2022). Nevertheless, time course measurements of SPAD and Fv/Fm, which enable model fitting and senescence parameter prediction, have rarely been captured in wheat under heat stress (Pinto et al., 2016; Šebela et al., 2020; Touzy et al., 2022), especially after brief heat during the early reproductive stage. In the present study, a wheat heat panel, including putatively heat-tolerant Australian and CIMMYT-nominated spring wheat genotypes, was exposed to a 5-day extreme (37°C/27°C) (Erena, 2018) or mild (34°C/28°C) heat stress during the early pollen developmental stage and analyzed for the effects on (i) pollen viability and seed set; (ii) spike formation and underlying tillering dynamics; (iii) leaf senescence measured with SPAD and Fv/Fm; and (iv) relationships among the reproductive, developmental, physiological, and yield-related traits.
Materials and Methods
Plant Material
In this study, three greenhouse experiments (i.e., Exp1, Exp2, and Exp3) were conducted during 2020–2021 at the controlled environment and glasshouse facilities at Rothamsted Research, Harpenden, United Kingdom (51.8094°N, 0.3561°W). For Exp1 and Exp2, the same set of 14 wheat genotypes was used, and 22 lines (seven overlapping with Exp1 and Exp2) were grown in Exp3 (refer to Supplementary Table 1 for genotype details). These lines are putatively heat-tolerant elite spring varieties (Sunstar, Sokoll, and Waagaan) parental lines and their pre-breeding materials (Cossani and Reynolds, 2015; Erena, 2018), as well as two lines from the United Kingdom included as controls.
Experimental Design and Heat Stress Treatment
Exp1
This experiment was conducted between June and September of 2020 and followed a split randomized complete block design (split RCBD) with four blocks/biological replicates. Fourteen genotypes (Supplementary Table 1) were randomly assigned to the whole plot within each block, and temperature treatments (control/CT and heat/HT) were assigned to the subplots within each whole plot. Two seeds were sown in separate pots filled with Rothamsted Standard compost (75% medium grade peat; 12% screened sterilized loam; 3% medium grade vermiculite; 10% 5 mm screened lime-free grit) and fertilized with Osmocote Exact for 3–4 months at the rate of 3.5 kg/m3. One week after sowing, seedlings were thinned down to one per pot and grown under natural light glasshouse conditions with a 16-h light period; they were supplemented with artificial light (230W LED; Kroptek Ltd., London, United Kingdom) if natural light intensity fell below 175 μmol/m2/s1. The temperature in the glasshouse was set at 21°C/15°C (day/night, actual value: 21.5 ± 0.4/16.3 ± 0.5°C) and the relative humidity (RH) was around 60/75% (day/night) (Supplementary Figure 1A). At the booting stage, plants with the primary tiller reached the targeted AIL (Supplementary Figure 2A) of 6 cm (actual value for each genotype in Supplementary Figure 2B) and were sequentially moved into Fitotron Modular Plant Growth Chambers (HGC1514; Weiss Technik UK Ltd., Loughborough, United Kingdom) for HT treatment (36.97 ± 0.03/26.95 ± 0.17°C, day/night). The light period was 16 h and the intensity was maintained around 600 μmol/m2/s1 at the plant level. The RH was maintained between 70 and 75% (Supplementary Figure 1D). Plants for CT treatment were kept in the glasshouse. After 5 days of HT treatment, heat-stressed plants were moved back to the glasshouse until final harvest.
Exp2
This experiment was conducted between August and December of 2020 and with the same set of 14 genotypes (Supplementary Table 1). The experiment design and plant cultivation conditions were similar as in Exp1. The temperature in the glasshouse was also set at 21°C/15°C (day/night, actual value: 20.6 ± 0.8/15.4 ± 0.7°C) and the RH was around 57/69% (day/night) (Supplementary Figure 1B). When the AIL of the primary tiller reached 2–3 cm (actual value for each genotype in Supplementary Figure 2C), plants for HT treatment were sequentially moved into the same growth chamber as Exp1 with the temperature of 37.02 ± 0.01/27.00 ± 0.01°C (day/night) (Supplementary Figure 1E). Plants for CT treatment were also moved into a similar growth chamber with the temperature of 21.01 ± 0.01/15.01 ± 0.01°C (Supplementary Figure 1E). The light period, intensity, and RH were similar between CT and HT treatments with settings of 16 h, 600 μmol/m2/s1, and 70–75% respectively. After 5 days in the growth chambers, plants were moved back to the glasshouse until final harvest.
Exp3
This experiment was conducted between January and May of 2021 and with 22 genotypes (Supplementary Table 1). The experiment design and plant cultivation conditions were the same as the previous two experiments. The temperature in the glasshouse was also set at 21°C/15°C (day/night, actual value: 21.7 ± 0.4/15.4 ± 0.3°C) and the RH was around 42/49% (day/night) (Supplementary Figure 1C). When the AIL of the primary tiller reached 2–3 cm (actual value in Supplementary Figure 2D), plants for HT treatment were sequentially moved into the growth chamber with a mild temperature of 34.00 ± 0.02/28.00 ± 0.01°C (day/night) (Supplementary Figure 1F). Plants for CT treatment were also moved into a similar growth chamber with the temperature of 21.01 ± 0.02/15.01 ± 0.01°C (Supplementary Figure 1F). The light period, intensity, and RH were similar between CT and HT treatments with settings of 16 h, 600 μmol/m/s, and 70–75% respectively. After 5 days in the growth chambers, plants were moved back to the glasshouse until final harvest.
Morphological, Phenological, and Physiological Measurements
On the day before (day 0) and after (day 6) HT treatment, the primary tiller of each plant was tagged and measured for AIL. Plant height (PH) was also recorded at these two time points in Exp2 and Exp3 (Supplementary Figure 1G). Chlorophyll SPAD and Fv/Fm (maximum potential quantum efficiency of Photosystem II) were measured at the same time as AIL and weekly thereafter (Supplementary Figure 1G). The SPAD measurement was performed with an MC-100 Chlorophyll Concentration Meter (Apogee Instruments, Inc., Logan, UT, United States). Fv/Fm was measured with a Pocket PEA (Hansatech Instruments Ltd., Norfolk, United Kingdom) after 15–20 min dark adaptation. For each plant, the mean SPAD value of measurements at the tip, middle, and bottom of flag leaf was obtained and one measurement of Fv/Fm was made in the middle of flag leaf. After the HT treatment, the date of each plant was recorded to calculate days to heading in Exp2 and Exp3. Physiological maturity of the spike on the tagged tiller was recorded as days to maturation. These measurements were conducted with four biological replicates of each genotype and treatment combination.
Measurement With Tagged Tillers/Spikes for Pollen Fertility and the Number of Grains Per Spike
During anthesis in Exp2 and Exp3, the fourth or fifth spikelet (counted from the bottom) was sampled from the tagged tiller. One anther from two florets at the bottom was photographed for a representative image and length of the anther was measured. In Exp3, the remaining five anthers from two florets at the bottom were pooled together for pollen viability analysis using staining with Lugol’s solution. Fully stained pollen was scored as viable, whereas partially stained or aberrant shaped pollen was scored as non-viable. At maturity, the number of filled grains of the tagged spike was counted and recorded as the number of grains per spike, and also, spike length (cm) and number of spikelets were measured. Four biological replicates were used for these measurements.
Measurement of Tillering Dynamics
In Exp1, development of extra young spikes after heat stress was observed. In Exp2 and Exp3, tiller number was therefore continuously counted for four biological replicates of CT and HT-treated plants of each genotype on the day (day 0) before and after (day 6) the 5-day HT treatment, and on weekly intervals thereafter until maximum tillering (Supplementary Figure 1G).
Yield-Related Measurements at Maturation
At maturity, the spikes per plant were distinguished into “old” spikes (labeled just before starting the HT treatment in Exp2 and Exp3) and “new” spikes, harvested separately, and then dried in oven at 40°C for 7 days prior to mechanical threshing and cleaning. The weight, number, length, and width were then determined for grain samples from old and new spikes separately with a scale and a MARViN digital seed analyzer (MARViTECH GmbH., Wittenburg, Germany). Grain yield per plant was calculated as the sum of grains from old and new spikes. The aboveground biomass for each plant was determined as the weight of all straw materials dried in an oven at 80°C for 48 h. The yield-related measurements were analyzed with four biological replicates.
Statistical Analysis
The data from the time course SPAD measurements were fitted using a generalized additive model (GAM) for each of the three experiments to estimate maximum SPAD (SPADmax), senescence onset (SenOnset), and senescence rate (SenRate) (Supplementary Figure 3). SPAD was predicted by a smooth function of time (days counted from stress initiation), with a separate smooth function fitted for each combination of genotype and treatment. The Exp1 model used eight basis functions, whereas Exp2 and Exp3 used seven basis functions. SPADmax was estimated from the fitted predicted model. SenOnset was calculated as the day that SPAD fell to 95% of the maximum SPAD. Senescence period was defined over 14 days from the onset or until the end of the measurement, whichever was shorter. SenRate was then calculated as the daily reduction of SPAD over the senescence period. GAMs were fitted in R package (version 3.6.1) using the “mgcv” package (version 1.8-35) (Wood, 2011).
All trait measurements and calculated parameters (Supplementary Table 2) were used for statistical analysis in R 4.0.3.1. First, descriptive statistics were summarized with the “describeBy” from the “psych” package (Revelle, 2020). The effects of genotype treatment and the interaction were obtained from ANOVA with the model fitted with “lmer” from the R package “lmerTest” (Kuznetsova et al., 2017); genotype, treatment, and their interaction were treated as fixed factors, while block and genotype nested in block were treated as random effects. Later, Tukey’s post hoc test was carried out for multiple test comparisons to identify genotypic variation. Estimated marginal means were calculated for each combination of genotype and treatment. Subsequently, for either CT or HT treatment, Pearson correlation coefficient table was calculated by using “tab_corr” from “sjPlot” package (Lüdecke, 2021) among measurements and pairwise-deletion method was used to account for missing data. For each experiment and temperature treatment, correlations among different traits were visualized as networks with the “qgraph” package (Epskamp et al., 2012).
Results
Heat-Impaired Pollen Fertility and Number of Grains Per Spike
To understand the effects of heat on pollen development and grain setting, the primary tiller of each plant was tagged and measured. In Exp1 and Exp2, the imposed severe heat treatment of 37°C/27°C caused nearly complete loss of grain setting for all genotypes, except for Paragon and Cadenza (Supplementary Figure 4). The anther morphology was also severely changed by the HT treatments indicating complete absence of viable pollen (Figures 1A,B). In Exp3, relatively mild heat stress (34°C/28°C) also significantly reduced anther length; however, this was less severe compared with Exp2 (Figures 1A,B) and pollen viability was therefore analyzed by staining with Lugol’s solution. The results showed considerable variation among genotypes, ranging from 0 to 60%. One line (SWBL1.1, a progeny between the cross of Sokoll and Weebill1) had the highest pollen viability (relative to control value), followed by SWES, SUN (Sunstar), and WBL1.2 (Weebill1) (Figure 1C). The number of grains of the tagged primary spike was also variable among the genotypes, with SUN showing the highest value relative to control value (Figure 1D). Further analysis found a positive correlation between pollen viability and the number of grains per spike under HT (Figure 1F), but not under CT treatment (Figure 1E).
Figure 1. The heat effects on wheat reproductive traits. Anther morphology after heat treatments in two representative genotypes is shown in (A, bar = 1 mm). Average anther length across all analyzed genotypes in Exp2 (37/27°C) and Exp3 (34/28°C) is shown in (B). Genotypic variation of pollen viability (C) and the number of grains per spike (D) under heat in Exp3. Representative pollen images for CT, HT-tolerant, and HT-sensitive were inserted (C, bar = 100 μm). The data were obtained from the primary tiller and calculated as relative to control values, which are shown on the top of each bar. Pearson correlation between pollen viability and number of grains per spike under control (E) and heat (F). Significance level: ***p < 0.001; **p < 0.01; *p < 0.05; *p < 0.05.
Heat-Stimulated Tillering/Spike Formation and Its Association With Yield
During the ripening stage of Exp1, the senescence status of tillers/spikes was clearly separated into two groups (Figure 2A) and tillers were therefore distinguished into old (pre-heat) and new (post-heat) spikes for each plant. About 1 week after heat treatment, new tiller outgrowth was noticed from the bottom of HT-stressed plants (Figures 2B,C). A final count of spikes found significantly more new spikes in the HT-treated plants compared with the CT plants in Exp1 (p < 0.001), Exp2 (p < 0.001), and Exp3 (p < 0.001) (Figures 2D–F), while the number of old spikes was similar between HT and CT conditions (Supplementary Figure 5). In addition, there was no significant interaction between treatment and genotype (Figures 2D–F), indicating that all genotypes responded similarly to the HT treatment in terms of new spike formation. The analysis of tillering dynamics in Exp2 and Exp3 showed that onset of new tiller development commenced at 2–3 weeks after the HT treatment, with a stronger effect observed in Exp2 (Figures 2G,H). Moreover, the more severe heat stress in Exp2 (37°C/27°C) also caused tiller retardation on day 6, 1 day after the end of the HT treatment (Figure 2G), but this was not observed under the milder heat stress condition in Exp3 (34°C/28°C) (Figure 2H).
Figure 2. The effects of heat treatments on tillering/spike formation. After heat (HT) treatment, new tillers/spikes (in yellow ellipse) were vigorously stimulated, while controls were not (CT, A). Tiller outgrowth from CT (B) and HT (C) around 1 week after stress treatment. A comparison of the number of new spikes between CT and HT in Exp1 (D), Exp2 (E), and Exp3 (F). The effects of genotype/G, treatment/T, interaction/G × T were indicated for each panel. Dynamic change of tillering before (day 0), after (day 6) the 5-day HT treatment (days 1–5), and weekly intervals in Exp2 (G) and Exp3 (H). Significance level: ***p < 0.001; **p < 0.01.
As new tillers developed after the HT treatment and extended the days to maturity of the plants, the aboveground biomass per plant (including both old and new tillers) was very similar between HT and CT treatments (Figures 3A–C). Nevertheless, the overall grain yield per plant was significantly reduced after the HT treatment in all three experiments (p < 0.001 for all) (Figures 3D–F). This was primarily due to heat-induced sterility in the old spikes (Figures 3G–I). However, heat-induced formation of new spikes gave rise to similar (Exp2, Figure 3K) or even significantly higher grain yield from new spikes in Exp1 and Exp3 (Figures 3J,L). The proportion of yield from new spikes after the HT treatment was therefore significantly higher than under CT conditions (Supplementary Figure 6). As sink size was reduced by limited seed setting of old spikes under HT condition, source supply became more than sufficient for the survived developing grains, and their width and length were significantly higher than grains of old spikes from control plants (Supplementary Table 2). By contrast, the grains from new spikes showed variable responses in terms of width and length (Supplementary Table 2).
Figure 3. Average values of wheat genotypes grown under high temperature (HT) or control (CT) across the three experiments are shown for biomass (A–C), grain yield per plant (D–F), grain yield from old spikes (G–I), and grain yield from new spikes (J–L). The effects of genotype (G), treatment (T), and interaction (G × T) are indicated for each panel.
Heat Effects on Plant Morphology, Phenology, and Chlorophyll Dynamics
When wheat plants were exposed to heat stress during the early booting stage, the increase in AIL (p < 0.001 for Exp1 and Exp2) and PH (p < 0.001 Exp2) during the 5-day treatments was significantly reduced by the HT of 37°C/27°C
compared with the CT of 21°C/15°C (Figures 4A,B,D). In contrast, the AIL (p = 0.065) and PH (p = 0.279) were marginally affected under the milder HT of 34°C/28°C in Exp3 (Figures 4C,E). HT treatments also changed plant phenology as indicated by the significantly reduced number of days to heading (DTH) (p < 0.001 for Exp2 and Exp3) (Figures 4F,G) and days to maturation (DTM) (p < 0.001 for Exp1, Exp2, and Exp3) (Figure 4H).
Figure 4. The heat effects on plant morphology and phenology. Comparison of the effect of heat (HT) and control (CT) treatments on the increase of auricle interval length (AIL, A–C), the increase in plant height (D,E), days to heading (F,G), and days to maturation (H).
To understand the physiological basis of changes in phenology, dynamic changes in SPAD and Fv/Fm were compared between CT and HT treatments. On day 6 (1 day after treatment), in comparison to the corresponding CT conditions, SPAD value was significantly reduced by the severe heat (37°C/27°C) in Exp1 (Figure 5A) and Exp2 (Figure 5B), but surprisingly increased slightly after the mild heat (34°C/28°C) in Exp3 and maintained a higher maximum SPAD value (Figures 5E,F). At later stages, however, an accelerated decrease in SPAD was observed under HT conditions in all three experiments, irrespective of heat stress intensity (Figures 5A,B,E). Based on the time course of SPAD measurements, GAMs were fitted to estimate maximum SPAD (SPADmax), senescence onset (SenOnset), and senescence rate (SenRate) for each combination of genotype and treatment. In Exp1 and Exp3, SenOnset from HT treatment was reproducibly and significantly advanced in comparison with CT conditions, whereas SenRate was similar between treatments (Figures 5C,F). By contrast, Exp2 showed an opposite response with similar SenOnset between treatments, but an increased SenRate under HT (Figure 5D). This variation between the three experiments may be due to variable intensities of natural sunlight. Even within the same experiment, some genotypes showed earlier SenOnset, while others showed faster SenRate under HT treatment (Supplementary Figures 7–9). In addition, the Fv/Fm value at day 6 was always significantly reduced by HT treatment in all three experiments indicating a negative effect of the HT on PSII (Supplementary Figure 10).
Figure 5. The effects of heat on SPAD (chlorophyll content index) and leaf senescence parameters. Comparison of the effect of heat (HT) and control (CT) treatments on the dynamic change of SPAD measured before (day 0) and 1 day after (day 6) a 5-day HT treatment, and in weekly intervals thereafter in Exp1 (A), Exp2 (B), and Exp3 (E). Senescence-related parameters, SPADmax (maximum SPAD value), SenOnset (Senescence onset day), and SenRate (Senescence rate), were compared between CT and HT for Exp 1 (C), Exp2 (D), and Exp3 (F). Significance level: ***p < 0.001; **p < 0.01; *p < 0.05.
Analysis of Trait Correlations From Different Experiments and Temperature Conditions
To understand the relationships among different traits across genotypes, correlations were calculated (Supplementary Table 3) and visualized as networks (Figure 6). Number of grain per spike (GpS) showed different correlations under control and HT conditions; In Exp1 and Exp2, there was no correlation between GpS and any other trait under HT, but under CT, it was strongly and positively correlated with the number of spikelets (SpikeletN) (r = 0.98*** for Exp1 and 0.81** for Exp2) and length (SpikeL) (r = 0.94*** for Exp1 and 0.79** for Exp2) of the tagged spike, as well as with biomass (r = 0.88***) and yield (r = 0.88***) in Exp1. In Exp3, GpS was also associated with different traits between CT and HT. The importance of induced new tillers and spikes after heat stress was corroborated by the reproducible strong positive correlations (r = 0.98***, 0.92***, 0.80*** for Exp1, Exp2, and Exp3, respectively) between grain yield of new tillers (GY.NT) and total grain yield (GY), observed in all three experiments (Figure 6 and Supplementary Table 3). This suggests a critical role of new spikes in mitigating heat-induced yield reduction. In addition, the morphological traits, increase in AIL and PH, were generally positively correlated with yield or biomass-related traits, regardless of temperature treatments. Ultimately, SPAD and Fv/Fm parameters were not consistently correlated with other traits from different experiments and treatments. In Exp1, SenOnset showed HT-specific weak positive correlations with a grain yield of old tillers (GY.OT) (r = 0.62*) and a spike number of new tillers (SpikeN.NT) (r = 0.69*); SenRate was closely related to yield traits in both Exp2 (r = 0.62* with GY.NT under CT; r = 0.79** with GY.NT and 0.69* with GY under HT) and Exp3 (r = 0.47* with GY and 0.53* with SpikeL under CT; r = 0.45* with GpS and 0.49* with SpikeL under HT), but not heat-specific; in Exp2, SPADmax was important, as it was strongly correlated with Spikelet (r = 0.66* under CT, 0.75** under HT) and SpikeL (r = 0.70* under CT, 0.74** under HT) of tagged spike (Figure 6 and Supplementary Table 3).
Figure 6. The effects of heat on trait relationships across different experiments. Correlation networks for Exp1 (A,B), Exp2 (C,D), and Exp3 (E,F). Only significant correlations are shown and the width of the edges indicated correlation r-value. Orange solid edges represent positive correlations, while blue dashed edges represent negative correlations. Trait abbreviations: SPADres: (SPAD at day 6—SPAD at day 0)/SPAD at day 0; Fv/Fmres: (Fv/Fm at day 6 -Fv/Fm at day 0)/(Fv/Fm at day 0); SPADmax, maximum value of SPAD; SenOnset, senescence onset time; SenRate, senescence rate; AILincrease, auricle interval length increase during the 5-day treatment; PHincrease, plant height increase during the 5-day treatment; PS, observed frequency of paired spikelet from all spikes; SR, observed frequency of sham ramification from all spikes; DTH, days to heading; DTM, days to maturation; SpikeL, spike length (tagged primary spike); SpikeletN, number of spikelets (tagged primary spike); GpS, number of grains per spike (tagged primary spike); SpikeN.OT, number of spikes from old tillers; SpikeN.NT, number of spikes from new tillers; GY.OT, grain yield of old tillers; GY.NT, grain yield of new tillers; GY, grain yield per plant (sum of old and new tillers); Biomass, the dry weight of all straw per plant; PV, pollen viability from the middle spikelet of tagged spike.
Discussion
Importance and Limitation of Pollen Viability as a Target Trait for Wheat Heat Research
In the present study, anther morphology was gradually affected under two levels of heat stress, 34°C/28°C (day/night) and 37°C/27°C, applied for 5 days during early booting stage coinciding with pollen development. The more severe heat stress condition in this study led to a complete loss of pollen viability, while results from a parallel study, in which the same 37°C/27°C heat treatment was applied that lasted for only 3 days (Erena, 2018), were similar to the 5-day, milder temperature (34°C/28°C) treatment, suggesting that both stress intensity and duration are critical to screening reproductive heat tolerance. Under the 34°C/28°C condition, pollen viability was considerably variable among genotypes. Two of the lines (SWB1.1 and SWES) with high pollen viability share one common parent, Sokoll, in their pedigree. Sokoll is an advanced wheat line derived from synthetic hexaploid wheat and has shown a yield advantage under terminal heat stress in other reports (Cossani and Reynolds, 2015; Thistlethwaite et al., 2020), although it did not show particularly high pollen viability after early booting-stage heat stress in this study. These results suggest stage-specific heat tolerance; therefore, it is necessary to pyramid tolerant traits across different developmental stages. Another parental line included in this study, Weebill1 (WBL1.1 and WBL1.2), has previously been reported to be tolerant to a wide range of variable environmental conditions (Singh et al., 2007). One of the most tolerant genotypes identified in this study was Sunstar, in agreement with data reported by Erena (2018) who also demonstrated the reproductive heat tolerance of Sunstar. These identified genotypes with heat tolerance during pollen development may be suitable donors for breeding and warrant further studies to understand the underlying genetic and molecular-physiological mechanisms. The importance of pollen viability is supported by its positive correlation with the number of grains per spike under heat stress. Interestingly, similar relationships have been reported in other crops (Xu et al., 2017; Shi et al., 2018) and abiotic stresses (Ji et al., 2010), indicating that pollen fertility is a general limiting factor for final grain number under suboptimal growth conditions. Therefore, it should be an important target trait for heat-related research and breeding. Nevertheless, the response of pollen viability to heat stress is highly dependent on the developmental stage when stress is applied (Saini and Aspinall, 1982; Prasad and Djanaguiraman, 2014) and it is thus important to consider genotypic differences and carefully target meiosis to microspore stage when applying heat stress to exclude confounding effects. Currently, the most widely used morphological marker for pollen developmental stage is AIL, which is also known as auricle distance (Ji et al., 2010; Erena, 2018; Bokshi et al., 2021). However, AIL corresponding to a specific pollen developmental stage varies among different genotypes (Erena, 2018) and must be determined for each genotype, which is laborious. Fortunately, progress has been made by non-destructive X-ray micro-computed tomography scanning (Fernández-Gómez et al., 2020), and integrating this with modeling could be a promising way to overcome difficulties with accurate identification of developmental stages of wheat pollen.
Utilizing Developmental Plasticity to Mitigate Heat Effects on Yield
The number of spikes per plant, interacting with spikelet number and floret fertility, determines grain number and thereby final yield. Our data show that a short episode of heat stress during early booting stage induced the development of new tillers and spikes, which is in agreement with other studies (Bányai et al., 2014; Chavan et al., 2019; Hütsch et al., 2019). Although tillering was initially inhibited under severe heat stress, new tillers started emerging at 2 weeks after recovery, corresponding to about 1 week after anthesis. This timing suggests that available photosynthates stored in vegetative tissue that cannot be translocated into grain due to spikelet sterility can be reallocated into the development of new tillers and spikes. Additional photo-assimilates for new tillers and spikes would be produced during recovery and this is reflected by its positive correlation with delayed onset of senescence (Figure 6C). The observed formation of new spikes after heat stress compensating for heat-induced biomass and yield losses under controlled environment conditions now needs to be corroborated under field conditions to ensure that it is a valid target trait for breeding. In addition, a higher frequency of paired spikelets (Boden et al., 2015) and sham ramification (Amagai et al., 2017) was observed in heat-treated plants and this may also be related to excessive source supply. Although these traits were not correlated with yield, they could contribute to understand mechanisms underlying such developmental abnormalities.
Accelerated Leaf Senescence After Brief Heat Stress During Early Booting Stage
Screening wheat for heat tolerance in the field is generally implemented by late-sowing to impose continuous terminal heat stress during grain filling, often resulting in accelerated leaf senescence (Bergkamp et al., 2018). In the present study, a similar stimulation of flag leaf senescence was observed after a brief episode of heat stress was applied during early booting. It is possible that plants are able to measure and memorize phenology or leaf age to program the senescence process (Woo et al., 2019). In our study, model fitting using SPAD time course data proved to be successful in identifying senescence parameters. Both earlier onset and faster senescence rate were identified and were closely related to accelerated leaf senescence, in agreement with similar results reported by Šebela et al. (2020). Heat-specific positive correlations between senescence onset (SenOnset), new spike formation (SpikeN.NT), and yield of old tillers (GY.OT) in Exp1 support the important role of late senescence. The observed positive associations between senescence rate (SenRate) and yield traits (grain yield/GY, number of grains per spike/GpS, spike length/SpikeL) in Exp2 and Exp3 suggest fast nutrient remobilization in high-yielding lines. Finally, both SPAD and Fv/Fm were reduced by heat immediately after the treatment (day6) in Exp1 and Exp2, but the mild temperature of 34°C/28°C only decreased Fv/Fm, not SPAD. These results indicate that Fv/Fm may be more sensitive and therefore a better parameter for heat tolerance evaluation (Cao et al., 2019). Therefore, these senescence-related parameters are useful for crop phenotyping, and integrating modeling with high-throughput imaging measurements will enable large-scale analysis.
Conclusion
In this study, a spring wheat panel, including heat-tolerant elite varieties and their pre-breeding lines, was dissected for reproductive, developmental, physiological, and yield responses along with their inter-relationships after a 5-day heat stress application during the early booting stage. In comparison with the control treatment, pollen viability from the tagged primary spike was significantly decreased by heat and subsequently reduced number of grains per spike. The heat stress, however, resulted in late tillering after the disruption of sink strength. Consequently, more new spikes were formed contributing to final yield and biomass, though an additional week was needed for the maturation of the late tillers. Flag leaf SPAD (Chlorophyll content index) and Fv/Fm (maximum potential quantum efficiency of Photosystem II) were reduced by heat stress. Model fitting with time course SPAD measurements showed accelerated leaf senescence by either earlier onset or faster senescence rate, and these parameters were associated with yield traits. Ongoing genomic and genetic studies will subsequently be used to dissect the mechanism of identified heat-tolerant genotypes (Sunstar, SWBL1.1). Taken together, these reproductive, developmental, and physiological traits could be further used as targets for understanding basic mechanisms and breeding heat-tolerant wheat.
Data Availability Statement
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author/s.
Author Contributions
SH conceived and supervised the project, together with EV-S and MP. JX designed and implemented the experiments and analyzed the data, with support from CL. MR and SD advised on the selection of genotypes included in this study and provided the seeds. SH, EV-S, SH-L, and MP held regular project planning discussions. JX wrote the manuscript, which was reviewed and edited by SH, CL, MR, SD, EV-S, SH-L, and MP. All authors contributed to the article and approved the submitted version.
Funding
This project was funded by the BBSRC UK-Mexico Newton Fund (BB/S012885/1 “Safeguarding Sonora’s Wheat from Climate Change”). SH and MP were supported by the Designing Future Wheat (DFW) Institute Strategic Programme (BB/P016855/1) and SH by NIAB.
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Acknowledgments
We would like to thank Tess Rose and Maria Oszvald for their support and suggestions on the manuscript, as well as Matthew Dale for advice on the experimental setup. We would also like to thank Fiona Gilzean, Jill Maple, and Jack Turner for taking excellent care of the plants and for their technical support, as well as Chris Hall for helping with the sample processing. Seeds were kindly provided by the CIMMYT and Nick Collins, University of Adelaide.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2022.886541/full#supplementary-material
Footnotes
References
Amagai, Y., Gowayed, S., Martinek, P., and Watanabe, N. (2017). The third glume phenotype is associated with rachilla branching in the spikes of tetraploid wheat (Triticum L.). Genet. Resour. Crop Evol. 64, 835–842. doi: 10.1007/s10722-017-0503-7
Azam, F. I, Chang, X., and Jing, R. (2015). Mapping QTL for chlorophyll fluorescence kinetics parameters at seedling stage as indicators of heat tolerance in wheat. Euphytica 202, 245–258. doi: 10.1007/s10681-014-1283-1
Bányai, J., Karsai, I., Balla, K., Kiss, T., Bedö, Z., and Láng, L. (2014). Heat stress response of wheat cultivars with different ecological adaptation. Cereal Res. Commun. 42, 413–425. doi: 10.1556/CRC.42.2014.3.5
Begcy, K., Weigert, A., Egesa, A. O., and Dresselhaus, T. (2018). Compared to australian cultivars, european summer wheat (Triticum aestivum) overreacts when moderate heat stress is applied at the pollen development stage. Agronomy 8:99. doi: 10.3390/agronomy8070099
Bergkamp, B., Impa, S. M., Asebedo, A. R., Fritz, A. K., and Jagadish, S. V. K. (2018). Prominent winter wheat varieties response to post-flowering heat stress under controlled chambers and field based heat tents. Field Crop Res. 222, 143–152. doi: 10.1016/j.fcr.2018.03.009
Bheemanahalli, R., Sunoj, V. S. J., Saripalli, G., Prasad, P. V. Vara, Balyan, H. S., Gupta, P. K., et al. (2019). Quantifying the impact of heat stress on pollen germination, seed set, and grain filling in spring wheat. Crop Sci. 59, 684–696. doi: 10.2135/CROPSCI2018.05.0292
Bhusal, N., Sharma, P., Sareen, S., and Sarial, A. K. (2018). Mapping QTLs for chlorophyll content and chlorophyll fluorescence in wheat under heat stress. Biol. Plant 62, 721–731. doi: 10.1007/s10535-018-0811-6
Boden, S. A., Cavanagh, C., Cullis, B. R., Ramm, K., Greenwood, J., Jean Finnegan, E., et al. (2015). Ppd-1 is a key regulator of inflorescence architecture and paired spikelet development in wheat. Nat. Plants 1, 1–6. doi: 10.1038/nplants.2014.16
Bokshi, A. I., Tan, D. K. Y., Thistlethwaite, R. J., Trethowan, R., and Kunz, K. (2021). Impact of elevated CO2and heat stress on wheat pollen viability and grain production. Funct. Plant Biol. 48, 503–514. doi: 10.1071/FP20187
Browne, R. G., Li, S. F., Iacuone, S., Dolferus, R., and Parish, R. W. (2021). Differential responses of anthers of stress tolerant and sensitive wheat cultivars to high temperature stress. Planta 254:4. doi: 10.1007/s00425-021-03656-7
Cai, C., Yin, X., He, S., Jiang, W., Si, C., Struik, P. C., et al. (2016). Responses of wheat and rice to factorial combinations of ambient and elevated CO2 and temperature in FACE experiments. Glob. Change Biol. 22, 856–874. doi: 10.1111/gcb.13065
Cao, Z., Yao, X., Liu, H., Liu, B., Cheng, T., Tian, Y., et al. (2019). Comparison of the abilities of vegetation indices and photosynthetic parameters to detect heat stress in wheat. Agric. For. Meteorol. 265, 121–136. doi: 10.1016/J.AGRFORMET.2018.11.009
Chavan, S. G., Duursma, R. A., Tausz, M., and Ghannoum, O. (2019). Elevated CO2 alleviates the negative impact of heat stress on wheat physiology but not on grain yield. J. Exp. Bot. 70, 6447–6459. doi: 10.1093/JXB/ERZ386
Cossani, C. M., and Reynolds, M. P. (2012). Physiological traits for improving heat tolerance in wheat. Plant Physiol. 160, 1710–1718. doi: 10.1104/pp.112.207753
Cossani, C. M., and Reynolds, M. P. (2015). Heat stress adaptation in elite lines derived from synthetic hexaploid wheat. Crop Sci. 55, 2719–2735. doi: 10.2135/cropsci2015.02.0092
Draeger, T., and Moore, G. (2017). Short periods of high temperature during meiosis prevent normal meiotic progression and reduce grain number in hexaploid wheat (Triticum aestivum L.). Theor. Appl. Genet. 130, 1785–1800. doi: 10.1007/s00122-017-2925-1
Dwivedi, S. K., Basu, S., Kumar, S., Kumari, S., Kumar, A., Jha, S., et al. (2017). Heat stress induced impairment of starch mobilisation regulates pollen viability and grain yield in wheat: study in Eastern Indo-Gangetic Plains. Field Crop Res. 206, 106–114. doi: 10.1016/J.FCR.2017.03.006
Epskamp, S., Cramer, A. O. J., Waldorp, L. J., Schmittmann, V. D., and Borboom, D. (2012). {qgraph}: network visualizations of relationships in psychometric data. J. Stat. Softw. 48, 1–18.
Erena, M. F. (2018). Genetic and Physiological Bases of Heat-Induced Floret Sterility in Wheat. Adelaide: University of Adelaide.
Erena, M. F., Lohraseb, I., Munoz-Santa, I., Taylor, J. D., Emebiri, L. C., and Collins, N. C. (2021). The WtmsDW locus on wheat chromosome 2B controls major natural variation for floret sterility responses to heat stress at booting stage. Front. Plant Sci. 12:376. doi: 10.3389/FPLS.2021.635397/BIBTEX
Fernández-Gómez, J., Talle, B., Tidy, A. C., and Wilson, Z. A. (2020). Accurate staging of reproduction development in Cadenza wheat by non-destructive spike analysis. J. Exp. Bot. 71, 3475–3484. doi: 10.1093/jxb/eraa156
Hassan, F. S. C., Solouki, M., Fakheri, B. A., Nezhad, N. M., and Masoudi, B. (2018). Mapping QTLs for physiological and biochemical traits related to grain yield under control and terminal heat stress conditions in bread wheat (Triticum aestivum L.). Physiol. Mol. Biol. Plants 24:1231. doi: 10.1007/S12298-018-0590-8
Hütsch, B. W., Jahn, D., and Schubert, S. (2019). Grain yield of wheat (Triticum aestivum L.) under long-term heat stress is sink-limited with stronger inhibition of kernel setting than grain filling. J. Agron. Crop Sci. 205, 22–32. doi: 10.1111/JAC.12298
Ji, X., Shiran, B., Wan, J., Lewis, D. C., Jenkins, C. L., Condon, A. G., et al. (2010). Importance of pre-anthesis anther sink strength for maintenance of grain number during reproductive stage water stress in wheat. Plant Cell Environ. 33, 926–942. doi: 10.1111/j.1365-3040.2010.02130.x
Kumar, P., Gupta, V., Singh, G., Dhakar, R., and Ramakrishnan, R. S. (2021). Assessment of terminal heat tolerance based on agro-morphological and stress selection indices in wheat. Cereal Res. Commun. 49, 217–226. doi: 10.1007/s42976-020-00112-2
Kuznetsova, A., Brockhoff, P. B., and Christensen, R. H. B. (2017). {lmerTest} package: tests in linear mixed effects models. J. Stat. Softw. 82, 1–26. doi: 10.18637/jss.v082.i13
Liu, B., Asseng, S., Müller, C., Ewert, F., Elliott, J., Lobell, D. B., et al. (2016). Similar estimates of temperature impacts on global wheat yield by three independent methods. Nat. Clim. Change 6, 1130–1136. doi: 10.1038/nclimate3115
Mirosavljević, M., Mikić, S., Župunski, V., Špika, A. K., Trkulja, D., Ottosen, C. O., et al. (2021). Effects of high temperature during anthesis and grain filling on physiological characteristics of winter wheat cultivars. J. Agron. Crop Sci. 207, 823–832. doi: 10.1111/jac.12546
Murchie, E. H., and Lawson, T. (2013). Chlorophyll fluorescence analysis: a guide to good practice and understanding some new applications. J. Exp. Bot. 64, 3983–3998. doi: 10.1093/jxb/ert208
Omidi, M., Siahpoosh, M. R., Mamghani, R., and Modarresi, M. (2014). The influence of terminal heat stress on meiosis abnormalities in pollen mother cells of wheat. Cytologia 79, 49–58. doi: 10.1508/cytologia.79.49
Pequeno, D. N. L., Hernández-Ochoa, I. M., Reynolds, M., Sonder, K., Molero-Milan, A., and Robertson, R. (2021). Climate impact and adaptation to heat and drought stress of regional and global wheat production. Environ. Res. Lett. 16:54070. doi: 10.1088/1748-9326/abd970
Pinto, R. S., Lopes, M. S., Collins, N. C., and Reynolds, M. P. (2016). Modelling and genetic dissection of staygreen under heat stress. Theor. Appl. Genet. 129, 2055–2074. doi: 10.1007/s00122-016-2757-4
Prasad, P. V. V., and Djanaguiraman, M. (2014). Response of floret fertility and individual grain weight of wheat to high temperature stress: sensitive stages and thresholds for temperature and duration. Funct. Plant Biol. 41, 1261–1269. doi: 10.1071/FP14061
Revelle, W. (2020). psych: Procedures for Psychological, Psychometric, and Personality Research Version 2.2.3.
Richardson, A. D., Duigan, S. P., and Berlyn, G. P. (2002). An evaluation of noninvasive methods to estimate foliar chlorophyll content. New Phytol. 153, 185–194. doi: 10.1046/j.0028-646X.2001.00289.x
Rieu, I., Twell, D., and Firon, N. (2017). Pollen development at high temperature: from acclimation to collapse. Plant Physiol. 173, 1967–1976. doi: 10.1104/pp.16.01644
Sade, N., Del Mar Rubio-Wilhelmi, M., Umnajkitikorn, K., and Blumwald, E. (2018). Stress-induced senescence and plant tolerance to abiotic stress. J. Exp. Bot. 69, 845–853. doi: 10.1093/jxb/erx235
Saini, H. S., and Aspinall, D. (1982). Abnormal sporogenesis in wheat (Triticum aestivum L.) induced by short periods of high temperature. Ann. Bot. 49, 835–846. doi: 10.1093/OXFORDJOURNALS.AOB.A086310
Šebela, D., Bergkamp, B., Somayanda, I. M., Fritz, A. K., and Krishna Jagadish, S. V. (2020). Impact of post-flowering heat stress in winter wheat tracked through optical signals. Agron. J. 112, 3993–4006. doi: 10.1002/AGJ2.20360
Sharma, D., Singh, R., Rane, J., Gupta, V. K., Mamrutha, H. M., and Tiwari, R. (2016). Mapping quantitative trait loci associated with grain filling duration and grain number under terminal heat stress in bread wheat (Triticum aestivum L.). Plant Breed. 135, 538–545. doi: 10.1111/pbr.12405
Shi, W., Li, X., Schmidt, R. C., Struik, P. C., Yin, X., and Jagadish, S. V. K. (2018). Pollen germination and in vivo fertilization in response to high-temperature during flowering in hybrid and inbred rice. Plant Cell Environ. 41, 1287–1297. doi: 10.1111/pce.13146
Shirdelmoghanloo, H., Cozzolino, D., Lohraseb, I., and Collins, N. C. (2016). Truncation of grain filling in wheat (Triticum aestivum) triggered by brief heat stress during early grain filling: association with senescence responses and reductions in stem reserves. Funct. Plant Biol. 43, 919–930. doi: 10.1071/FP15384
Singh, R. P., Huerta-Espino, J., Sharma, R., Joshi, A. K., and Trethowan, R. (2007). High yielding spring bread wheat germplasm for global irrigated and rainfed production systems. Euphytica 157, 351–363. doi: 10.1007/s10681-006-9346-6
Slafer, G. A., Savin, R., and Sadras, V. O. (2014). Coarse and fine regulation of wheat yield components in response to genotype and environment. Field Crop Res. 157, 71–83. doi: 10.1016/j.fcr.2013.12.004
Stone, P. J., and Nicolas, M. E. (1995). Effect of timing of heat stress during grain filling on two wheat varieties differing in heat tolerance. I. Grain growth. Aust. J. Plant Physiol. 22, 927–934. doi: 10.1071/PP9950927
Talukder, A. S. M. H. M., McDonald, G. K., and Gill, G. S. (2014). Effect of short-term heat stress prior to flowering and early grain set on the grain yield of wheat. Field Crop Res. 160, 54–63. doi: 10.1016/J.FCR.2014.01.013
Thistlethwaite, R. J., Tan, D. K. Y., Bokshi, A. I., Ullah, S., and Trethowan, R. M. (2020). A phenotyping strategy for evaluating the high-temperature tolerance of wheat. Field Crop Res 255:107905. doi: 10.1016/j.fcr.2020.107905
Touzy, G., Lafarge, S., Redondo, E., Lievin, V., Decoopman, X., Le Gouis, J., et al. (2022). Identification of QTLs affecting post-anthesis heat stress responses in European bread wheat. Theor. Appl. Genet. 135, 947–964. doi: 10.1007/S00122-021-04008-5
Trnka, M., Rötter, R. P., Ruiz-Ramos, M., Kersebaum, K. C., Olesen, J. E., Žalud, Z., et al. (2014). Adverse weather conditions for European wheat production will become more frequent with climate change. Nat. Clim. Change 4, 637–643. doi: 10.1038/nclimate2242
Woo, H. R., Kim, H. J., Lim, P. O., and Nam, H. G. (2019). Leaf senescence: systems and dynamics aspects. Annu. Rev. Plant Biol. 70, 347–376. doi: 10.1146/annurev-arplant-050718-095859
Xu, J., Wolters-Arts, M., Mariani, C., Huber, H., and Rieu, I. (2017). Heat stress affects vegetative and reproductive performance and trait correlations in tomato (Solanum lycopersicum). Euphytica 213:156. doi: 10.1007/s10681-017-1949-6
Yeh, C. H., Kaplinsky, N. J., Hu, C., and Charng, Y. Y. (2012). Some like it hot, some like it warm: phenotyping to explore thermotolerance diversity. Plant Sci. 195, 10–23. doi: 10.1016/j.plantsci.2012.06.004
Zampieri, M., Ceglar, A., Dentener, F., and Toreti, A. (2017). Wheat yield loss attributable to heat waves, drought and water excess at the global, national and subnational scales. Environ. Res. Lett. 12:064008. doi: 10.1088/1748-9326/aa723b
Keywords: heat stress, booting, pollen viability, tillering, SPAD and Fv/Fm, wheat
Citation: Xu J, Lowe C, Hernandez-Leon SG, Dreisigacker S, Reynolds MP, Valenzuela-Soto EM, Paul MJ and Heuer S (2022) The Effects of Brief Heat During Early Booting on Reproductive, Developmental, and Chlorophyll Physiological Performance in Common Wheat (Triticum aestivum L.). Front. Plant Sci. 13:886541. doi: 10.3389/fpls.2022.886541
Received: 28 February 2022; Accepted: 04 April 2022;
Published: 16 May 2022.
Edited by:
Pasala Ratnakumar, Indian Institute of Oilseeds Research (ICAR), IndiaReviewed by:
Renu Munjal, Chaudhary Charan Singh Haryana Agricultural University, IndiaAshok Singamsetti, Banaras Hindu University, India
Copyright © 2022 Xu, Lowe, Hernandez-Leon, Dreisigacker, Reynolds, Valenzuela-Soto, Paul and Heuer. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Sigrid Heuer, sigrid.heuer@niab.com