Abstract
Breeders agree that leaf senescence is a favorable process for wheat seed yield improvement due to the remobilization of leaf nutrients. However, several studies have suggested that staying green may be an important strategy for further increasing wheat yields. In this study, we performed a comparative transcriptome analysis between wheat cultivars CN17 and CN19 after heading and also measured photosynthetic parameters, chlorophyll (Chl) contents, and antioxidant enzyme activities at various time points after heading. The physiological and biochemical indexes revealed that CN17 exhibited a functionally stay-green phenotype while CN19 did not. We identified a total of 24,585 and 34,410 differential expression genes between genotypes at two time-points and between time-points in two genotypes, respectively, and we also found that 3 (37.5%) genes for leaf senescence, 46 (100%) for photosynthesis – antenna protein, 33 (70.21%) for Chl metabolism and 34 (68%) for antioxidative enzyme activity were upregulated in CN17 compared with CN19 during leaf senescence, which could be regulated by the differential expression of SAG39 (senescence-associated gene 39), while 22 (100%) genes for photosynthesis – antenna proteins, 6 (46.15%) for Chl metabolism and 12 (80%) for antioxidative enzyme activity were upregulated in CN17 compared with CN19 before the onset of leaf senescence. Here, we further clarified the expression profiles of genes associated with a functional stay-green phenotype. This information provides new insight into the mechanism underlying delayed leaf senescence and a new strategy for breeders to improve wheat yields.
Introduction
From 1910 to 2005, the global net primary production (HANPP) doubled, but the global population increased fourfold (). Thus, crop production must be further improved to meet the rising demand of the rapidly growing population (; ). As the most widely grown crop, wheat provides a great deal of calories and proteins for humans (). In the past, considerable improvements in wheat yield have been achieved through resistance improvement (Petersen et al., 2015; Ma et al., 2019), by the utilization of dwarfing genes (Thomas, 2017), and via morphological or type breeding (Pecetti and Annicchiarico, 1998). However, substantially increasing crop yields further via traditional breeding strategies has become increasingly difficult.
Leaf senescence, the final phase of leaf development, also has a large influence on crop yield (). Hence, promoting the favorable onset and speed of leaf senescence could constitute a potential strategy to further improve crop yields. Some mutants referred to as “stay-green mutants,” which exhibit a normal green phenotype even after the onset of leaf senescence, have been identified in several species, such as rice (; ), wheat () and maize (; ). Stay-green genotypes have been classified into two categories, functional stay-green genotypes and cosmetic stay-green genotypes, the difference of which is based on whether they have effective photosynthesis during the process of leaf senescence (). Many breeders usually treat the stay-green phenotype as an obviously negative indictor of yield traits (; ), which could partly result from the stay-green phenotype of many mutants being only cosmetic rather than functional, and thus, stay-green mutants have been largely ignored.
However, functional stay-green mutants might be helpful to further increase crop yields because they can maintain photosynthesis during the later stages of growth, thus leading to the accumulation of increased amounts of photosynthetic products (Luo et al., 2006; Wang et al., 2021). Various studies have suggested that the cosmetic stay-green phenotype usually results from mutations in genes involved in leaf senescence, photosynthesis, chlorophyll (Chl) metabolism and antioxidative enzyme activity (; Zhang et al., 2016). Unlike for cosmetic stay-green genotypes, accurate identification of the phenotypes of functional stay-green genotypes is difficult mainly because there is continuous variation in the characteristics caused by the complex genetic makeup (Yoo et al., 2007) and in part because of the influences of other factors, especially those under field conditions, such as available water (), temperature (), drought () and diseases (; ). Notably, because of these difficulties, developing functional stay-green commercial cultivars is the objective of only a few select breeders.
Fortunately, three functional commercial stay-green wheat cultivars, CN12, CN17 and CN18, with high yield and desirable agronomic traits were first developed via chromosome engineering, and they have been widely planted in southwestern China for approximately 20 years (Ren et al., 2003, 2004b). Previous studies have demonstrated that these cultivars are functional stay-green genotypes (Luo et al., 2006, 2009; ; ), and their stay-green-related functions could be explained by the chloroplast ultrastructure regeneration hypothesis (Luo et al., 2013), of which the core idea was that the photosynthetic apparatus could be regenerated in the senescing leaf when the stacked grana thylakoid of PSII was retained in the early stage of leaf senescence. However, there is still no strong evidence at the gene expression level supporting this hypothesis, although some differentially expressed sequence tag (EST) sequences have been obtained by suppression subtractive hybridization (SSH; Luo et al., 2013). In addition, the non-stay-green control wheat cultivar MY11 employed in a previous study was also not a contemporaneous counterpart, and thus, the evidence from both gene expression and key pathways supporting the views of the chloroplast ultrastructure regeneration hypothesis are still unclear.
To identify the key genes and core metabolic pathways underlying delayed leaf senescence and provide more solid evidence for the stay-green function, the stay-green cultivar CN17 and its counterpart control, the non-stay-green cultivar CN19, were used for a comparative transcriptome analysis both before and after the onset of leaf senescence. In this study, we were interested in two classes of genes. The first genes (type I) were differentially expressed between CN19 and CN17 before leaf senescence (constitutive differential expression), with no difference in expression or opposite differential expression after leaf senescence; this differential expression accompanied a different tendency of gene expression changes from before leaf senescence to after leaf senescence in CN19 and CN17, respectively. The other class of genes (type II) were not differentially expressed between CN19 and CN17 before leaf senescence (non-constitutive differential expression), and the differential expression of these genes after leaf senescence occurred alongside a different tendency of expression changes from before leaf senescence to after leaf senescence in CN19 and CN17, respectively. In addition, photosynthetic parameters, Chl contents and antioxidant enzyme activities were also measured throughout the whole process.
Materials and methods
Plant materials and growing environment
The wheat cultivar CN17 (stay green) and the control cultivar CN19 (non-stay green), which were developed and released in 2003 (Ren et al., 2003, 2004a) and are high yielding, were employed in the research. CN17 was selected from the cross (91S-23/A302), while CN19 was selected from the cross (Q1104A/R935). CN17 and CN19 have similar leaf area, plant height, tiller number and developmental progress including both heading and flowering time because they are contemporaneous cultivars. These plant materials were sown by hand in clay soil on October 31st, 2016 at the agricultural research station of Sichuan Agricultural University (SAU) in the Wenjiang district (latitude 30°43′N, longitude 103°52′E), Chengdu city, Sichuan Province, China. That area has a temperate climate, and during the experimental period average air temperatures were approximately 12°C and the experiment received approximately 390 mm of rain. However, during the period from anthesis (March15th) to maturity, the average air temperatures were 20.5°C and 347 mm rainfall was received. During the growing season, 6 g and 14 g N ammonium nitrate were applied at the one node and heading stage, respectively. Additionally, diseases and pests were effectively controlled with fungicides and pesticides as required. The experiment was conducted in a randomized complete block design with three replications. Each block between two 0.5 m width alleys consisted of twelve 3 m-long rows and rows were spaced at 0.33 m apart, in which the interval between two adjacent plants was 0.15 m. Thirty-one plants of each cultivar with similar development progress were randomly selected from the inside and marked in field plots for subsequent experiments, of which 10 plants of each genotype were chosen for measuring the photosynthetic indexes at 0, 10, 20, 30, and 40 days after heading (DAH). Similarly, 15 plants of each genotype were used to determine the biochemical indexes at the same five time points, and six plants of each genotype (three biological replicates for each) at 20 DAH (before the onset of leaf senescence, the key time point for determining the stay-green phenotype or not) and 40 DAH (during leaf senescence, the key time point for determining the final difference) were selected for RNA sequencing (RNA-seq).
Measurement of photosynthetic indexes
Photosynthetic indexes, including the net photosynthetic rate (Pn), intercellular carbon dioxide concentration (Ci), stomatal conductance (Gs), the transpiration rate (Tr), and the ratio of internal to atmospheric CO2 concentration (Ci/Ca), were measured with a Li-6400-02B photosynthetic apparatus (LI-COR, Lincoln, NE, United States) as previously described (). Measurements of the middle of flag leaves were made between 9.00 AM and 11.00 AM under the following conditions of 0.55–0.65 kPa vapor pressure deficit, 1000 μmol/m2/s actinic light intensity and 21–24°C air temperature. A measurement was logged 1 min after readings stabilized. The mean of three independent readings of each plant was taken to represent the value of the given plant, while the mean of ten plants was used as the parameter for the genotype.
Determination of Chl contents
Flag leaves of each genotype (three replications) at each time point were used to determine the Chl contents according to the methods of .
Measurements of antioxidant enzyme activity
Flag leaves of CN17 and CN19 (three replications) collected every ten DAH were used for measuring the activity of antioxidant enzymes, including superoxide dismutase (SOD), catalase (CAT) and peroxidase (POD), according to the methods described by ; Pigeolet et al. (1990), and , respectively.
Data statistics and analysis
The data were analyzed via SPSS Statistics software version 23 (IBM, Armonk, NY, United States). The means were calculated, and the differences between two cultivars and at each of the five time points were compared by analysis of variance (ANOVA) and independent-sample t-tests.
Ribonucleic acid extraction and RNA sequencing
In studied region, leaf senescence of wheat plants grown under natural field conditions usually occurs at 26–28 DAH, and therefore, three flag leaves of each genotype at 20 DAH were selected as representative samples before the onset of leaf senescence and at 40 DAH as representative samples after the onset of leaf senescence for RNA extraction. TRIzol reagent (Invitrogen, Carlsbad, CA, United States) was used to extract RNA in accordance with the manufacturer’s instructions. The RNA of the three flag leaves of each cultivar at each time point (constituting three biological replicates) was subjected to RNA-seq at Biomarker Technologies Corporation based on the Illumina HiSeq high-throughput sequencing platform. The raw data have been deposited in the BIG Data Center under accession number CRA007229.
Alignment of RNA sequencing data, differential gene expression analysis and functional annotations
The raw sequence data were screened by removing adapter and primer sequences, reads containing more than 10% N bases and reads of low quality. Then, the cleaned were aligned to the wheat reference genome IWGSC RefSeq v1.0 by using TopHat2 software to obtain the location information on the reference genome, as well as information concerning the specific characteristics of the sample sequences (). The expression level of all the genes was calculated using StringTie (Pertea et al., 2015) software. Then, DESeq was used to analyze the genes that were differentially expressed between samples (Wang et al., 2009). A fold change ≥ 2 and a false discovery rate (FDR) < 0.01 were used as screening criteria for the detection of differentially expressed genes (DEGs). The sequences of the obtained DEGs were aligned to the sequence information within the Kyoto Encyclopedia of Genes and Genomes (KEGG) database (), Pfam database (), SwissProt database (), and non-redundant (NR) database () to annotate gene functions via the BLAST method. KEGG pathway enrichment analysis of the DEGs was performed using the KEGG Orthology Based Annotation System (KOBAS; Mao et al., 2005). When we identified the constitutive DEGs truly responsible for leaf senescence prior to leaf senescence occurring, the genes with the same expression tendency in both CN19 and CN17 from before leaf senescence to after leaf senescence and the similar differential expression between these cultivars after leaf senescence were treated only as controls and were not included as true constitutive DEGs. Likewise, when we identified the induced DEGs truly responsible for leaf senescence, the genes with the same expression tendency in both CN19 and CN17 from before leaf senescence to after leaf senescence and with similar differential expression between them before leaf senescence were also treated only as controls and were not included in the genes with truly differential expression.
Single-nucleotide polymorphisms of differentially expressed genes detected via comparison of expressed sequences
Genome Analysis Toolkit (GATK; McKenna et al., 2010) software was used to detect SNPs, and SnpEff () software was used to annotate the variants and predict the effects of mutations. Reliable SNPs were screened from among those with no more than 3 consecutive single-base mismatches in 35 bp and for whose the quality after sequential depth standardization was greater than 2.0. To identify the key genes responsible for the stay-green phenotype, the expressed sequences of the DEGs were compared between samples; the same SNPs in the reference genome and CN19 were used as controls. In addition, three biological replicates of each genotype at each time point for the same SNP were selected to avoid false-positive SNPs.
Results
Phenotypic observations
With respect to CN17, a stay-green wheat cultivar, a large portion of leaves remained green past 40 DAH, whereas the flag leaves of the control cultivar (CN19) began to turn yellow at 30–32 DAH, which indicated that the onset of leaf senescence in CN19 usually occurred before 30 DAH, which was a normal senescence timeframe. However, CN17 exhibited a delayed onset of senescence; the visible change from green to yellow occurred at 44–45 DAH, and all of its leaves were yellow at 52–54 DAH.
Changes in photosynthetic parameters
The Pn of CN17 was significantly higher than that of CN19 at both 0 and 40 DAH at the P = 0.01 level, and the change in the Pn exhibited the opposite tendency between CN17 and CN19 from 0 to 10 DAH, during which there was a decrease in the Pn in CN17 and an increase in CN19. In addition, the rate of decrease in the Pn was lower for CN17 than for CN19 from 30 to 40 DAH, and therefore, the Pn at 40 DAH was maintained at more than 60% and less than 10% of the Pn at 0 DAH in CN17 and CN19, respectively (Figure 1A).
FIGURE 1
The other photosynthetic parameters, including the Ci, Gs, Tr and Ci/Ca, of both CN17 and CN19 tended to change similarly during the whole process (Figures 1B–E), but there was usually a lower degree of change in these parameters (excluding Gs) from 10 to 30 DAH in CN17 compared with CN19, especially from 30 to 40 DAH. Thus, compared with CN19, CN17 presented consequently higher Gs and Tr values as well as Pn values (Figures 1A,C,D) and lower Ci and Ci/Ca values (Figures 1B,E) at 40 DAH.
Chl components and contents
Compared with CN19, CN17 had significantly higher total Chl, chlorophyll a (Chl a) and chlorophyll b (Chl b) contents at all time points except for 20 DAH (Figures 2A–C). The contents of Chl a and Chl b in CN17 at 40 DAH were 68 and 84% of those at 0 DAH, respectively, while the contents of Chl a and Chl b in CN19 at 40 DAH were only 4.8 and 8.5% of those at 0 DAH (Figures 2A,B). Following a faster decrease, the Chl a/b ratio in CN19 significantly increased from 10 to 20 DAH, while a similar increase occurred during the later stage from 30 to 40 DAH in CN17 (Figure 2D).
FIGURE 2
Antioxidant enzyme activity
From 0 to 10 DAH, SOD activity significantly increased (P ≤ 0.05) in CN17 but significantly decreased (P ≤ 0.01) in CN19, and both CN19 and CN17 exhibited the same SOD activity change trend from 10 to 40 DAH, with a more rapid rate of change in CN17 than in CN19 from 20 to 40 DAH (Figure 3A). There was a similar change tendency of CAT activity in both CN19 and CN17 from 0 to 20 DAH (Figure 3B). After 20 DAH, the CAT activity increased from 20 to 30 DAH and slightly decreased from 30 to 40 DAH in CN17, while in CN19, it rapidly continuously decreased from 20 to 40 DAH (Figure 3B); thus, the CAT activity in CN17 was significantly higher than that in CN19 at 40 DAH (Figure 3B). The opposite change trend for POD activity occurred between CN17 and CN19 from 10 to 20 DAH and from 30 to 40 DAH, respectively (Figure 3C).
FIGURE 3
Analysis of differentially expressed genes
On average, 9.58 Gb of data for each RNA sample were obtained from the leaves of both CN17 and CN19 at 20 and 40 DAH (Supplementary Table 1). Clean reads of each sample were mapped to the wheat reference genome, with an average alignment rate of approximately 80.44% (Supplementary Table 1). At 20 DAH, the number of genes differentially expressed between CN19 and CN17 was 8,160, with 4,057 upregulated and 4,103 downregulated in CN17 compared with CN19 (Table 1). At 40 DAH, the number of DEGs between them reached 16,425, of which 7,792 were upregulated and 8,633 were downregulated in CN17 (Table 1). Further analysis revealed that there were only 9,395 genes differentially expressed in CN17 between 20 and 40 DAH, while there were up to 25,015 in CN19 (Table 1).
TABLE 1
| DEG set | Number of DEGs | Upregulated | Downregulated |
| CN17–20 DAH vs. CN17–40 DAH | 9,395 | 4,278 | 5,117 |
| CN19–20 DAH vs. CN17–20 DAH | 8,160 | 4,057 | 4,103 |
| CN19–20 DAH vs. CN19–40 DAH | 25,015 | 11,771 | 13,244 |
| CN19–40 DAH vs. CN17–40 DAH | 16,425 | 7,792 | 8,633 |
Comparison of DEGs in CN17 and CN19 at 20 DAH and 40 DAH.
DEG, differentially expressed gene; DAH, Days after heading. N = 3 three biological replications.
Classification of differentially expressed genes involved in the pathways possibly related to the stay-green phenotype
A graph of the statistical analysis results of KEGG pathway enrichment between CN19 and CN17 at 40 DAH showed that photosynthesis – antenna proteins were ranked first (Figure 4), and further analysis revealed that there were 22 (green ellipse) and 46 (red ellipse) DEGs involved in this pathway between CN19 and CN17 at 20 and 40 DAH, respectively, of which only 5 DEGs were detected both before and after leaf senescence. Moreover, there were 55 (yellow ellipse) and 69 (blue ellipse) DEGs in CN17 and CN19, respectively, between 20 DAH and at 40 DAH, of which 52 were common to both cultivars (Figure 5A). These results suggested that the expression of genes involved in the photosynthesis pathway exhibited obvious differences between the stay-green and non-stay green plants, and therefore, the pathways closely related to photosynthesis, including leaf senescence, Chl metabolism and reactive oxygen species (ROS) metabolism, were also analyzed with similar methods.
FIGURE 4
FIGURE 5
A total of 9 DEGs were involved in leaf senescence, which is usually considered antagonistic to staying green, and there were only 2 (green ellipse) and 8 (red ellipse) genes differentially expressed between CN19 and CN17 at 20 and 40 DAH, respectively, of which the two genes differentially expressed between CN19 and CN17 at 20 DAH were also detected between them at 40 DAH (Figure 5B). In addition, there were 5 common DEGs following the leaf senescence process in CN17 and CN19 (Figure 5B).
Of the 95 DEGs involved in the Chl metabolism pathway, 13 (green ellipse) and 47 (red ellipse) DEGs were identified between CN19 and CN17 at 20 DAH and 40 DAH, respectively, and the number of DEGs between 20 and 40 DAH in CN17 (yellow ellipse) was lower than that in CN19 (blue ellipse; Figure 5C). Further analysis found that the genes differentially expressed between CN19 and CN17 were mainly involved in Chl biosynthesis and that there were more DEGs between these cultivars at 40 DAH than at 20 DAH (Supplementary Table 2). Similarly, among the 68 DEGs involved in antioxidative enzyme pathways, the number of DEGs between genotypes at 40 DAH (red ellipse) was greater than that at 20 DAH (green ellipse), whereas the number of genes that were involved in this pathway as well as Chl metabolism and that were differentially expressed before and after leaf senescence in CN17 (yellow ellipse) was also lower than that in CN19 (blue ellipse; Figure 5D).
Identification of genes involved in the functional stay-green phenotype
Based on the foundational knowledge about wheat leaf senescence, two types of DEGs were interesting. Type I genes were differentially expressed between CN19 and CN17 before leaf senescence, and their differences in expression between genotypes after leaf senescence were not significant or opposite, which was possibly due to the different tendencies of expression changes during leaf development. Likewise, type II DEGs were identified between CN19 and CN17 after leaf senescence; these genes were not differentially expressed between genotypes before leaf senescence, which also resulted from the different responses to leaf senescence. In type I, 2, 6 and 3 genes were involved in photosynthesis – antenna proteins, Chl metabolism and antioxidative enzymes, respectively, and no DEGs involved in leaf senescence were classified as type I (Table 2). In contrast, 16, 17, 27 and 1 DEG involved in photosynthesis – antenna protein, Chl metabolism, antioxidative enzyme activity and leaf senescence, respectively, were classified as type II (Table 2). Interestingly, all 16 DEGs involved in photosynthesis – antenna proteins were found to encode Chl a-b binding protein.
TABLE 2
| Expression pattern | Classes | Types of DEGs | CN19–20 DAH vs. CN17–20 DAH | CN19–40 DAH vs. CN17–40 DAH | CN17–20 DAH vs. CN17–40 DAH | CN19–20 DAH vs. CN19–40 DAH | Porphyrin and chlorophyll metabolism | Photosynthesis –antenna proteins | Antioxidative enzymes | Senescence-associated genes |
| 95 | 73 | 68 | 9 | |||||||
| Constitutive expression | A1 | DEGs with similar expression trends from before leaf senescence to after leaf senescence in both genotypes exhibited similar expression between genotypes both before and after leaf senescence | Down | Down | / | / | 1 | 2 | ||
| Up | Up | Down | Down | 1 | 3 | 1 | ||||
| Up | Up | / | / | 1 | 4 | |||||
| Down | Down | Up | Up | 2 | ||||||
| A2 | DEGs with similar expression trends from before leaf senescence to after leaf senescence in both genotypes with no differential expression between genotypes after leaf senescence | Down | / | Down | Down | 1 | ||||
| Up | / | Down | Down | 1 | 15 | |||||
| Up | / | / | / | 2 | ||||||
| A3 | DEGs with different expression trends from before leaf senescence to after leaf senescence in both genotypes, with similar expression between genotypes both before and after leaf senescence | Down | Down | / | Down | 1 | ||||
| Up | Up | / | Down | 1 | 1 | 2 | ||||
| Up | Up | Up | Down | 1 | ||||||
| Up | Up | Up | / | 1 | ||||||
| A4 (Type I) | DEGs with different gene expression change trends from before leaf senescence to after leaf senescence in CN19 and CN17, with no differential expression or opposite expression between them after leaf senescence compared with before leaf senescence | Down | / | / | Down | 4 | 1 | |||
| Up | / | Down | / | 2 | 2 | 2 | ||||
| Non-constitutive expression | B1 | DEGs with similar expression trends from before leaf senescence to after leaf senescence in both genotypes, with similar expression between genotypes both before and after leaf senescence | / | / | Down | Down | 13 | 9 | 2 | |
| / | / | Up | Up | 4 | 1 | |||||
| B2 | DEGs with similar expression tendency from before leaf senescence to after leaf senescence in both genotypes exhibited differential expression between genotypes after leaf senescence, respectively | / | Up | Down | Down | 14 | 25 | 6 | ||
| / | Down | Up | Up | 9 | 4 | 3 | ||||
| / | Up | / | / | 1 | 3 | 2 | ||||
| / | Down | / | / | 1 | ||||||
| B3 | DEGs with different expression trends from before leaf senescence to after leaf senescence in both genotypes, with no differential expression between genotypes after leaf senescence | / | / | / | Down | 19 | 7 | 1 | ||
| / | / | / | Up | 2 | 3 | |||||
| / | / | Down | / | 2 | 1 | |||||
| B4 (Type II) | DEGs with different gene expression changetrends from before leaf senescence to after leaf senescence in CN19 and CN17, with differential expression between them after leaf senescence | / | Up | / | Down | 15 | 16 | 15 | 1 | |
| / | Down | / | Up | 2 | 9 | |||||
| / | Down | Down | / | 1 | ||||||
| / | Up | Up | / | 2 |
Numbers of DEGs of different classes in each pathway.
Bold text indicates gene of interest. “Up” and “Down” represent upregulated and downregulated DEGs in the latter compared with the former, respectively.
Based on the sequences of the DEGs, we found 5 SNPs between CN19 and CN17. Only one synonymous SNP belonged to type I and was in TraesCS1B01G059100 on chromosome 1B, encoding glutathione S-transferase. The other four non-synonymous SNPs (nsSNPs) were identified in DEGs of type II, and two nsSNPs were discovered in TraesCS2B01G447300, which encodes oxygen-dependent coproporphyrinogen-III oxidase, which participates in Chl biosynthesis, located at loci 640104224 and 640111777 of chromosome 2B. Similarly, a nsSNP was discovered in TraesCS7B01G192500 on chromosome 7B, which is associated with photosynthesis – antenna proteins, and a SNP downstream from TraesCS7A01G090400, which encodes SOD, on chromosome 7A was also detected.
Discussion
CN17 exhibits a true functional stay-green phenotype and delayed leaf senescence
The advantages of delayed leaf senescence in terms of crop breeding include sustained photosynthesis during the later grain filling stage (). Previous studies have suggested that wheat leaf senescence usually occurs at approximately 21 days post-anthesis and at approximately 25 DAH; in CN17, a widely planted wheat cultivar with a 1BL/1RS translocated chromosome, leaf senescence is delayed by approximately 14 days and the senescence speed is slower compared with that needed for normal leaf senescence (Luo et al., 2006, 2009, 2013). In the present study, a decrease in Pn occurred in both CN17 and CN19 from 20 DAH, while the rate of decrease was lower for CN17 than for CN19, especially during the period from 30 to 40 DAH (Figure 1A). The parallel change between the Pn and Gs in CN17 from 10 to 40 DAH (Figures 1A,C) indicated that the photosynthetic ability of CN17 during senescence was mainly regulated by stomatal factors rather than by the activity of the photosynthetic apparatus. In addition, the lower Tr accompanying the equivalent or larger Gs of CN17 before 40 DAH (Figures 1C,D) suggested that, the water-use efficiency of CN17 is also high, the phenomenon of which is usually associated with drought resistance and the stay-green trait (). Further comprehensive analysis of the changes in the Pn, Gs, and Ci suggested that the persistent increase in Ci in CN19 after 20 DAH possibly resulted mainly from the loss of photosynthetic ability rather than the decrease in Gs (Figures 1A–C). Consequently, the Pn of CN17 was significantly higher than that of CN19 at 40 DAH.
Previous studies have suggested that Chl b may be important for maintaining the stability of the light-harvesting Chl a/b-protein (LHCP) complex and the photosynthetic apparatus because it is primarily localized within the LHCP, while Chl a is a common component of all Chl-containing proteins (; Luo et al., 2013). In this study, the total Chl, Chl a and Chl b contents were significantly higher in CN17 than in CN19 after leaf senescence (Figures 2A–C). This could provide some explanations in terms of apparatus for increased photosynthetic potential in CN17.
The hallmark of the functional stay-green phenotype is the prolonged high level of photosynthesis to provide increased amounts of carbohydrates, ultimately leading to the production of larger grains (). Previous studies have shown that CN17 produced larger grains than did the control MY11 (Luo et al., 2006), which was partially contributed to the smaller grains resulted from the variety degeneration of MY11, which was released in 1983. Some studies have shown that CN17 and CN19 have similar disease resistance and high yield (Luo et al., 2005, 2008, 2009). In this study, the rapid changes in the amounts of both Chl a and Chl b in CN17 (Figures 2A,B) ruled out the possibility that the stay-green phenotype was caused by the blocking of Chl degradation. Therefore, the delayed leaf senescence of CN17 was a true functional stay-green trait and would be helpful for further yield improvement; this phenotype is governed by a high photosynthetic ability that is possibly associated with both Chl content and antioxidant enzyme activity especially SOD and CAT after 20 DAH, which was also indirectly supported by the previous studies suggesting that both enhanced oxidative stress and reduced antioxidant activity may be a crucial factor to accelerate the process of leaf senescence (Prochazkova et al., 2001; Wang et al., 2016).
The functional stay-green phenotype is associated with almost all pathways related to leaf senescence
The most noticeable difference between the functional stay-green phenotype and the cosmetic stay-green phenotype is whether these stay-green plants can maintain a high photosynthetic ability during leaf senescence (); the cosmetic stay-green phenotypes are related to only or mainly Chl degradation pathways (Pružinská et al., 2005; ).
In the present study, KEGG pathway enrichment analysis revealed that 46 genes differentially expressed between CN17 and CN19 after leaf senescence were mostly involved in photosynthesis – antenna proteins (Figures 4, 5A), and for 22 genes that were differentially expressed between the cultivars before leaf senescence, this pathway also ranked second (Figure 5A and Supplementary Figure 1). These findings indicated that the stay-green phenotype of CN17 could result from photosynthesis-related pathways such as leaf senescence, Chl biosynthesis and antioxidant enzyme activity because photosynthesis – antenna proteins play a key role in photosynthesis. We further identified 8, 47, and 50 genes that were differentially expressed between CN17 and CN19 at 40 DAH and involved in leaf senescence, porphyrin and Chl metabolism, and antioxidant enzymes, respectively (Figures 5B–D). In addition, the different change trends and values of the various photosynthetic parameters (Figure 1), Chl components and contents (Figure 2), and antioxidant enzyme activity (Figure 3) also indirectly suggested that these genes were differentially expressed.
There were also some genes that were involved in the abovementioned metabolic pathways and that were differentially expressed between the two genotypes at 20 DAH, but the number of DEGs and the change trend of some were obviously different. It is reasonable that the functional stay-green phenotype is associated with almost all the pathways related to leaf senescence, although we did not completely exclude the possibility of DEGs due to the differences between genotypes.
Identification of key genes responsible for the functional stay-green phenotype
Previous studies have demonstrated that unlike the cosmetic stay-green phenotype, the functional stay-green phenotype usually involves multiple metabolic pathways and is controlled by many loci (; Ren et al., 2022). The difference in leaf senescence between CN17 and CN19 is very conspicuous, especially in the later grain filling stage; therefore, despite their genetic background being different, as demonstrated by a pedigree analysis, these cultivars were used to detect DEGs ().
There are usually more DEGs detected when two genotypes have differences in their genetic background, and therefore, we need to pay more attention to screening and identifying the DEGs related to target traits. In this study, to identify the key genes responsible for the functional stay-green phenotype, we used the genes that were differentially expressed before leaf senescence as controls; these DEGs could have resulted from the differences in genetic background. At the same time, we further eliminated the common DEGs exhibiting similar expression change trends from 20 to 40 DAH between CN17 and CN19, which possibly resulted from the similar developmental processes of other traits, such as grain filling. Based on this idea, only type I and type II DEGs were interesting to us. There were 11 type I and 61 type II DEGs identified in the study (Table 2), which indicated that non-constitutive differential expression was the predominant mechanisms responsible for the stay-green phenotype, although we did not completely ignore the contribution of the constitutive DEGs involved in the stay-green phenotype.
Further analysis revealed that only one senescence-associated gene (SAG39) out of 72 interesting DEGs belonged to type II genes that we concerned (Table 2 and Supplementary Table 3); interestingly, a previous study confirmed that SAG could negatively regulate leaf senescence (). In this study, SAG39 was upregulated in CN17 compared with CN19 at 40 DAH, and there was no difference in expression between cultivars at 20 DAH. In addition, SAG39 was downregulated at 40 DAH compared with 20 DAH in CN19, while the expression change from before leaf senescence to after leaf senescence was not different in CN17 (Table 2 and Supplementary Table 3). These findings suggested that SAG39 negatively regulates leaf senescence and delays leaf senescence mainly by regulating the stability of Chl-protein complexes and that the upregulated expression of SAG39 in CN17 at 40 DAH possibly inhibited leaf senescence. Both 2 type I and all 16 type II DEGs involved in photosynthesis – antenna proteins encoded Chl a-b binding proteins (Table 2 and Supplementary Table 4); these proteins retained Chl in Chl-protein complexes and represent elegant control points of Chl degradation () and regulation of ROS balance (Xu et al., 2012).
Interestingly, 17 type II DEGs were involved in Chl metabolism, of which 15 DEGs involved in the Chl biosynthesis were upregulated in CN17 compared with CN19 at 40 DAH; however, following the leaf senescence process, they were downregulated in CN19 (Table 2 and Supplementary Table 2). In contrast, the other two DEGs involved in Chl degradation were downregulated in CN17 compared with CN19 at 40 DAH, while following the leaf senescence process, they were upregulated in CN19 (Table 2 and Supplementary Table 2). One of these DEGs encodes the Chl b reductase NYC1, which directly regulates Chl degradation (Teng et al., 2021), and the other encodes heme oxygenase 1, which putatively and negatively regulates Chl biosynthesis by enhancing the other branch pathways of porphyrin and inhibiting the number of porphyrins available for Chl biosynthesis (Mahawar and Shekhawat, 2018). In addition, 6 type I DEGs involved in Chl metabolism were detected, of which two (Table 2 and Supplementary Table 2) DEGs encoding protochlorophyllide reductase A, were upregulated in CN17 compared with CN19 at 20 DAH, while following leaf senescence, they were downregulated in CN17, and the expression change showed no difference in CN19. The measured Chl component and content parameters were well agreed with the differential expression of the genes involved in Chl metabolism (Figure 2).
The changes in antioxidant enzyme activity after heading (Figure 3) indicated that the genes encoding these enzymes had very versatile role in leaf senescence (Table 2 and Supplementary Table 5), and the 3 type I and 27 type II DEGs exhibiting many patterns of expression changes between CN17 and CN19 from before leaf senescence to after leaf senescence (Table 2 and Supplementary Table 5) demonstrated the multiple roles of the DEGs in regulating leaf senescence. However, the four genes encoding CAT or CAT isozymes showed consistent expression changes; these genes were upregulated in CN17 compared with CN19 at 40 DAH, while following leaf senescence, CN19 was downregulated and CN17 did not differ in terms of expression (Table 2 and Supplementary Table 5). Therefore, the predominant role of genes encoding CAT or CAT isozymes positively regulated the stay-green phenotype, especially after leaf senescence, possibly by inhibiting the degradation of Chl-containing proteins (Zavaleta-Mancera et al., 2007) and by defending against cell oxidative damage from ROS ().
Single nucleotide polymorphisms in the key genes possibly enhance functionally stay-green phenotype
Usually, nsSNPs are important for producing new genes with functions that differ from those of old genes. In this study, we identified four nsSNPs from three genes. TraesCS2B01G447300, which encodes oxygen-dependent coproporphyrinogen-III oxidase, on 2B, TraesCS7B01G192500, which encodes Chl a-b binding protein 1B-21, on 7B, and a downstream SNP from TraesCS7A01G090400, which encodes SOD, on 7A were also detected.
TraesCS2B01G447300 has 8 exons, within which two nsSNPs were located in exons 1 and 5, inducing a mutation from aspartic acid to asparagine and from alanine to aspartic acid, respectively, in the encoded protein. Coproporphyrinogen III oxidase (CPOX, EC 1.3.3.3) is a crucial enzyme in Chl biosynthesis, catalyzing coproporphyrinogen III to form protoporphyrinogen IX. Previous studies have indicated that mutations in the porphyrinogen III oxidase gene promoted early senescence in Arabidopsis thaliana ().
There were 4 exons in TraesCS7B01G192500, in which a nsSNP was present in the first exon, causing a mutation from threonine to alanine in the encoded protein. This gene encodes a light-harvesting Chl protein complex (LHC) protein, commonly known as an antenna, which absorbs and transports light energy to the reaction center of the photosystem (). In addition, LHCs also play crucial roles in photoprotection, dissipating excess excitation energy as heat and reducing the formation of ROS, thus reducing light-induced oxidative damage (; ). The remaining SNP was found in TraesCS7A01G090400, which encodes SOD and is involved in delaying senescence.
None of these three genes were differentially expressed before leaf senescence, but all three were substantially highly differentially expressed after leaf senescence between CN17 and CN19. In addition, different change tendencies during leaf senescence were detected between the genotypes. These findings indicated that these genes may indeed be differentially expressed as a result of the different responses to leaf senescence. These SNPs, especially the nsSNPs within the key genes, may contribute to the functional stay-green phenotype.
Conclusion
Altogether, CN17 was confirmed to exhibit a true functional stay-green phenotype with high photosynthetic potential. This feature could mainly be attributed to the following processes or gene expression differences, as shown in Figure 6. Firstly, the upregulation of SAG39 in CN17 could maintain the stability of chlorophyll-protein complexes by negatively regulating the leaf senescence process. Secondly, the upregulation of Chl synthesis-related genes and downregulation of Chl degradation-related genes may play a key role in maintaining high contents of chlorophyll and thus keeping the integrity of the chlorophyll-proteins. Thirdly, high expression of the genes encoding catalase or catalase isozyme could be helpful to maintain the perfect structure of chlorophyll-proteins.
FIGURE 6
Statements
Data availability statement
The original contributions presented in this study are included in the article/Supplementary material. The sequencing data are publicly available at the BIG Data Center (https://ngdc.cncb.ac.cn) under accession number: CRA007229. Further inquiries can be directed to the corresponding author.
Author contributions
PL, FT, and CC conceived the research. PL, TR, ZL, QL, and HY designed the experiments. QL, HY, JG, QH, and SZ conducted the experiments. QL, HY, JG, and PL analyzed the experimental data. QL, HY, and PL wrote the manuscript. All authors contributed to the article and approved the submitted version.
Funding
This work was financially supported by the Provincial Science and Technology Foundation of Sichuan, China (2022ZDZX0016), and the Science and Technology Research Program of Chongqing Municipal Education Commission of China (KJQN201903602).
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.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2022.971927/full#supplementary-material
References
1
AntoniettaM.FanelloD. D.AcciaresiH. A.GuiametJ. J. (2014). Senescence and yield responses to plant density in stay green and earlier-senescing maize hybrids from Argentina.Field Crops Res.155111–119. 10.1016/j.fcr.2013.09.016
2
ArchanaR.VinodK. K.KrishnanS. G.VadhanaD. C.BhowmickP. K.SinghV. J.et al (2021). Stay green behaviour of a novel mutant, PSG16 shows complex inheritance and functional relations with grain yield in rice.Indian J. Genet. Plant Breed.81495–504. 10.31742/IJGPB.81.4.1
3
BairochA.ApweilerR. (2000). The SWISS-PROT protein sequence database and its supplement TrEMBL in 2000.Nucleic Acids Res.2845–48. 10.1093/nar/28.1.45
4
BarryC. S. (2009). The stay-green revolution: Recent progress in deciphering the mechanisms of chlorophyll degradation in higher plants.Plant Sci.176325–333. 10.1016/j.plantsci.2008.12.013
5
BeauchampC.FridovichI. (1971). Superoxide dismutase: Improved assays and an assay applicable to acrylamide gels.Anal. Biochem.44276–287. 10.1016/0003-2697(71)90370-8
6
BellemareG.BartlettS. G.ChuaN. H. (1982). Biosynthesis of chlorophyll a/b-binding polypeptides in wild type and the chlorina f2 mutant of barley.J. Biol. Chem.2577762–7767. 10.1016/s0021-9258(18)34446-6
7
BestwickC. S.BrownI. R.MansfieldJ. W. (1998). Localized changes in peroxidase activity accompany hydrogen peroxide generation during the development of a nonhost hypersensitive reaction in Lettuce.Plant Physiol.1181067–1078. 10.1104/pp.118.3.1067
8
BlandinoM.PilatiA.ReyneriA. (2009). Effect of foliar treatments to durum wheat on flag leaf senescence, grain yield, quality and deoxynivalenol contamination in North Italy.Field Crops Res.114214–222. 10.1016/j.fcr.2009.08.008
9
BorrellA. K.HammerG. L.DouglasA. C. L. (2000b). Does maintaining green leaf area in Sorghum improve yield under drought? I. Leaf growth and senescence.Crop Sci.401026–1037. 10.2135/cropsci2000.4041026x
10
BorrellA. K.HammerG. L.HenzellR. G. (2000a). Does maintaining green leaf area in Sorghum improve yield under drought? II. Dry matter production and yield.Crop Sci.401037–1048. 10.2135/cropsci2000.4041037x
11
CalderiniD. F.CastilloF. M.Arenas-MA.MoleroG.ReynoldsM. P.CrazeM.et al (2021). Overcoming the trade-off between grain weight and number in wheat by the ectopic expression of expansin in developing seeds leads to increased yield potential.New Phytol.230629–640. 10.1111/nph.17048
12
ChenJ.LiangY.HuX.WangX.TanF.ZhangH.et al (2010). Physiological characterization of ‘stay green’ wheat cultivars during the grain filling stage under field growing conditions.Acta Physiol. Plant.32875–882. 10.1007/s11738-010-0475-0
13
ChibaneN.CaicedoM.MartinezS.MarcetP.RevillaP.OrdásB. (2021). Relationship between delayed leaf senescence (Stay-green) and agronomic and physiological characters in maize (Zea mays L.).Agronomy11:276. 10.3390/agronomy11020276
14
ChristopherJ. T.ManschadiA. M.HammerG. L.BorrellA. K. (2008). Developmental and physiological traits associated with high yield and stay-green phenotype in wheat.Aust. J. Agric. Res.59354–364. 10.1071/ar07193
15
CingolaniP.PlattsA.WangL. L.CoonM.NguyenT.WangL.et al (2012). A program for annotating and predicting the effects of single nucleotide polymorphisms.SnpEff. Fly (Austin)680–92. 10.4161/fly.19695
16
Clifton-BrownJ. C.LewandowskiI.BangerthF.JonesM. B. (2002). Comparative responses to water stress in stay-green, rapid- and slow senescing genotypes of the biomass crop, Miscanthus.New Phytol.154335–345. 10.1046/j.1469-8137.2002.00381.x
17
CookJ. P.AcharyaR. K.MartinJ. M.BlakeN. K.KhanI. J.HeoH. Y.et al (2021). Genetic analysis of stay-green, yield, and agronomic traits in spring wheat.Crop Sci.61383–395. 10.1002/csc2.20302
18
Crafts-BrandnerS. J.BelowF. E.WittenbachV. A.HarperJ. E.HagemanR. H. (1984). Differential senescence of maize hybrids following ear removal.Plant Physiol.74368–373. 10.1104/pp.74.2.368
19
CrepinA.CaffarriS. (2015). The specific localizations of phosphorylated Lhcb1 and Lhcb2 isoforms reveal the role of Lhcb2 in the formation of the PSI-LHCII supercomplex in Arabidopsis during state transitions.Biochim. Biophys. Acta (BBA) Bioenerg.18471539–1548. 10.1016/j.bbabio.2015.09.005
20
DengY. Y.LiJ. Q.WuS. F.ZhuY.ChenY. W.HeF. C. (2006). Integrated nr database in protein annotation system and its localization.Comput. Eng.3271–74.
21
FinnR. D.BatemanA.ClementsJ.CoggillP.EberhardtR. Y.EddyS. R.et al (2013). Pfam: The protein families database.Nucleic Acids Res.42D222–D230. 10.1093/nar/gkt1223
22
FrónaD.SzenderákJ.Harangi-RákosM. (2019). The challenge of feeding the world.Sustainability11:5816. 10.3390/su11205816
23
FuJ. D.YanY. F.KimM. Y.LeeS. H.LeeB. W. (2011). Population-specific quantitative trait loci mapping for functional stay-green trait in rice (Oryza sativa L.).Genome54235–243. 10.1139/g10-113
24
GuoC. Y.WuG. H.XingJ.LiW. Q.TangD. Z.CuiB. M. (2013). A mutation in a coproporphyrinogen III oxidase gene confers growth inhibition, enhanced powdery mildew resistance and powdery mildew-induced cell death in Arabidopsis.Plant Cell Rep.32687–702. 10.1007/s00299-013-1403-8
25
GuoY.RenG.ZhangK.LiZ.MiaoY.GuoH. (2021). Leaf senescence: Progression, regulation, and application.Mol. Hortic.1:5. 10.1186/s43897-021-00006-9
26
HörtensteinerS. (2009). Stay-green regulates chlorophyll and chlorophyll-binding protein degradation during senescence.Trends Plant Sci.14155–162. 10.1016/j.tplants.2009.01.002
27
JanssonS. (1999). A guide to the Lhc genes and their relatives in Arabidopsis.Trends Plant Sci.4236–240. 10.1016/s1360-1385(99)01419-3
28
JiangG. H.HeY. Q.XuC. G.LiX. H.ZhangQ. (2004). The genetic basis of stay-green in rice analyzed in a population of doubled haploid lines derived from an indica by Japonica cross.Theor. Appl. Genet.108688–698. 10.1007/s00122-003-1465-z
29
JoshiA. K.KumariM.SinghV. P.ReddyC. M.KumarS.RaneJ.et al (2007). Stay green trait: Variation, inheritance and its association with spot blotch resistance in spring wheat (Triticum aestivum L.).Euphytica15359–71. 10.1007/s10681-006-9235-z
30
KanehisaM.GotoS. (2000). KEGG: Kyoto encyclopedia of genes and genomes.Nucleic Acids Res.2827–30. 10.1093/nar/28.1.27
31
KimD.PerteaG.TrapnellC.PimentelH.KelleyR.SalzbergS. L. (2013). TopHat2: Accurate alignment of transcriptomes in the presence of insertions, deletions and gene fusions.Genome Biol.14:R36. 10.1186/gb-2013-14-4-r36
32
KrausmannF.ErbK. H.GingrichS.HaberlH.BondeauA.GaubeV.et al (2013). Global human appropriation of net primary production doubled in the 20th century.Proc. Natl. Acad. Sci. U.S.A.11010324–10329. 10.1073/pnas.1211349110
33
KumariM.SinghV. P.TripathiR.JoshiA. K. (2007). “Variation for staygreen trait and its association with canopy temperature depression and yield traits under terminal heat stress in wheat,” in Wheat production in stressed environments, edsBuckH. T.NisiJ. E.SalomónN. (Dordrecht: Springer Netherlands), 357–363.
34
Kura-HottaM.SatohK.KatohS. (1987). Relationship between photosynthesis and chlorophyll content during leaf senescence of rice seedlings.Plant Cell Physiol.281321–1329. 10.1093/oxfordjournals.pcp.a077421
35
LacanD.BaccouJ. C. (1998). High levels of antioxidant enzymes correlate with delayed senescence in nonnetted muskmelon fruits.Planta204377–382. 10.1007/s004250050269
36
LiQ.ZhongS.SunS.FatimaS. A.ZhangM.ChenW.et al (2017). Differential effect of whole-ear shading after heading on the physiology, biochemistry and yield index of stay-green and non-stay-green wheat genotypes.PLoS One12:e0171589. 10.1371/journal.pone.0171589
37
LiX.LiuT.ChenW.ZhongS.ZhangH.TangZ.et al (2015). Wheat WCBP1 encodes a putative copper-binding protein involved in stripe rust resistance and inhibition of leaf senescence.BMC Plant Biol.15:239. 10.1186/s12870-015-0612-4
38
LiZ.PengJ.WenX.GuoH. (2012). Gene network analysis and functional studies of senescence-associated genes reveal novel regulators of Arabidopsis leaf senescence.J. Integr. Plant Biol.54526–539. 10.1111/j.1744-7909.2012.01136.x
39
LichtenthalerH. (1987). Chlorophylls and carotenoids: Pigments of photosynthetic biomembranes.Methods Enzymol.148350–382. 10.1016/0076-6879(87)48036-1
40
LiguoriN.PerioleX.MarrinkS. J.CroceR. (2015). From light-harvesting to photoprotection: Structural basis of the dynamic switch of the major antenna complex of plants (LHCII).Sci. Rep.5:15661. 10.1038/srep15661
41
LuoP. G.DengK. J.HuX. Y.LiL. Q.LiX.ChenJ. B.et al (2013). Chloroplast ultrastructure regeneration with protection of photosystem II is responsible for the functional ‘stay-green’ trait in wheat.Plant Cell Environ.36683–696. 10.1111/pce.12006
42
LuoP. G.RenZ. L.ZhangH. Q.ZhangH. Y. (2005). Identification, chromosome location, and diagnostic markers for a new gene (YrCN19) for resistance to wheat stripe rust.Phytopathology951266–1270. 10.1094/phyto-95-1266
43
LuoP. G.ZhangH. Y.ShuK.WuX. H.ZhangH. Q.RenZ. L. (2009). The physiological genetic effects of 1BL/1RS translocated chromosome in “stay green” wheat cultivar CN17.Can. J. Plant Sci.891–10. 10.4141/cjps07072
44
LuoP. G.ZhangH. Y.ShuK.ZhangH. Q.LuoH. Y.RenZ. L. (2008). Stripe rust (Puccinia striiformis f. sp. tritici) resistance in wheat with the wheat-rye 1BL/1RS chromosomal translocation.Can. J. Plant Pathol.30254–259. 10.1080/07060661.2008.10540540
45
LuoP.RenZ.WuX.ZhangH.ZhangH.FengJ. (2006). Structural and biochemical mechanism responsible for the stay-green phenotype in common wheat.Chin. Sci. Bull.512595–2603. 10.1007/s11434-006-2175-0
46
MaH.ZhangX.YaoJ.ChengS. (2019). Breeding for the resistance to Fusarium head blight of wheat in China.Front. Agric. Sci. Eng.6:251–264. 10.15302/j-fase-2019262
47
MahawarL.ShekhawatG. S. (2018). Haem oxygenase: A functionally diverse enzyme of photosynthetic organisms and its role in phytochrome chromophore biosynthesis, cellular signalling and defence mechanisms.Plant Cell Environ.41483–500. 10.1111/pce.13116
48
MaoX.CaiT.OlyarchukJ. G.WeiL. (2005). Automated genome annotation and pathway identification using the KEGG Orthology (KO) as a controlled vocabulary.Bioinformatics213787–3793. 10.1093/bioinformatics/bti430
49
McKennaA.HannaM.BanksE.SivachenkoA.CibulskisK.KernytskyA.et al (2010). The genome analysis toolkit: A mapreduce framework for analyzing next-generation DNA sequencing data.Genome Res.201297–1303. 10.1101/gr.107524.110
50
PecettiL.AnnicchiaricoP. (1998). Agronomic value and plant type of Italian durum wheat cultivars from different eras of breeding.Euphytica999–15. 10.1023/a:1018346901579
51
PerteaM.PerteaG. M.AntonescuC. M.ChangT. C.MendellJ. T.SalzbergS. L. (2015). StringTie enables improved reconstruction of a transcriptome from RNA-seq reads.Nat. Biotechnol.33290–295. 10.1038/nbt.3122
52
PetersenS.LyerlyJ. H.WorthingtonM. L.ParksW. R.CowgerC.MarshallD. S.et al (2015). Mapping of powdery mildew resistance gene Pm53 introgressed from Aegilops speltoides into soft red winter wheat.Theor. Appl. Genet.128303–312. 10.1007/s00122-014-2430-8
53
PigeoletE.CorbisierP.HoubionA.LambertD.MichielsC.RaesM.et al (1990). Glutathione peroxidase, superoxide dismutase, and catalase inactivation by peroxides and oxygen derived free radicals.Mech. Ageing Dev.51283–297. 10.1016/0047-6374(90)90078-t
54
ProchazkovaD.SairamR. K.SrivastavaG. C.SinghD. V. (2001). Oxidative stress and antioxidant activity as the basis of senescence in maize leaves.Plant Sci.161765–771. 10.1016/S0168-9452(01)00462-9
55
PružinskáA.TannerG.AubryS.AndersI.MoserS.MüllerT.et al (2005). Chlorophyll breakdown in senescent Arabidopsis leaves. characterization of chlorophyll catabolites and of chlorophyll catabolic enzymes involved in the degreening reaction.Plant Physiol.13952–63. 10.1104/pp.105.065870
56
RenT.FanT.ChenS.ChenY.OuX.JiangQ.et al (2022). Identification and validation of quantitative trait loci for the functional stay green trait in common wheat (Triticum aestivum L.) via high-density SNP-based genotyping.Theor. Appl. Genet.1351429–1441. 10.1007/s00122-022-04044-9
57
RenZ.JiangH.ZhangH.TanF.YanB.ZhangH.et al (2004b). Wheat new cultivar chuannong 19 with high yields, high quality and high resistance to disease.J. Triticeae Crops24:159.
58
RenZ.ZhangH.FuT.TanF. (2003). Wheat new cultivars chuannong 12 and chuannong 17 with high yields and high resistance to disease and senescence.J. Triticeae Crops23:97.
59
RenZ.ZhangH.TanF.JiangH.YanB.ZhangH.et al (2004a). Wheat new cultivar chuannong 18 with high yields, high quality, resistance to disease and senescence.J. Triticeae Crops24:158.
60
TengK.TanP.GuanJ.DongD.LiuL.GuoY.et al (2021). Functional characterization of the chlorophyll b reductase gene NYC1 associated with chlorophyll degradation and photosynthesis in Zoysia japonica.Environ. Exp. Bot.191:104607. 10.1016/j.envexpbot.2021.104607
61
ThomasS. G. (2017). Novel Rht-1 dwarfing genes: Tools for wheat breeding and dissecting the function of DELLA proteins.J. Exp. Bot.68354–358. 10.1093/jxb/erw509
62
WangF.LiuJ.ZhouL.PanG.LiZ.ChengF. (2016). Senescence-specific change in ROS scavenging enzyme activities and regulation of various SOD isozymes to ROS levels in psf mutant rice leaves.Plant Physiol. Biochem.109248–261. 10.1016/j.plaphy.2016.10.005
63
WangL.FengZ.WangX.WangX.ZhangX. (2009). DEGseq: An R package for identifying differentially expressed genes from RNA-seq data.Bioinformatics26136–138. 10.1093/bioinformatics/btp612
64
WangP.HouS.WenH.WangQ.LiG. (2021). Chlorophyll retention caused by STAY-GREEN (SGR) gene mutation enhances photosynthetic efficiency and yield in soybean hybrid Z1.Photosynthetica5937–48. 10.32615/ps.2020.076
65
XuY. H.LiuR.YanL.LiuZ. Q.JiangS. C.ShenY. Y.et al (2012). Light-harvesting chlorophyll a/b-binding proteins are required for stomatal response to abscisic acid in Arabidopsis.J. Exp. Bot.631095–1106. 10.1093/jxb/err315
66
YooS. C.ChoS. H.ZhangH.PaikH. C.LeeC. H.LiJ.et al (2007). Quantitative trait loci associated with functional stay-green SNU-SG1 in rice.Mol. Cells2483–94.
67
Zavaleta-ManceraH. A.López-DelgadoH.Loza-TaveraH.Mora-HerreraM.Trevilla-GarcíaC.Vargas-SuárezM.et al (2007). Cytokinin promotes catalase and ascorbate peroxidase activities and preserves the chloroplast integrity during dark-senescence.J. Plant Physiol.1641572–1582. 10.1016/j.jplph.2007.02.003
68
ZhangJ.LiH.XuB.LiJ.HuangB. (2016). Exogenous melatonin suppresses dark-induced leaf senescence by activating the superoxide dismutase-catalase antioxidant pathway and down-regulating chlorophyll degradation in excised leaves of perennial ryegrass (Lolium perenne L.).Front. Plant Sci.7:1500. 10.3389/fpls.2016.01500
Summary
Keywords
wheat (Triticum aestivum L.), leaf senescence, photosynthesis, chlorophyll, antioxidant enzyme, transcriptome
Citation
Li Q, Yang H, Guo J, Huang Q, Zhong S, Tan F, Ren T, Li Z, Chen C and Luo P (2022) Comparative transcriptome analysis revealed differential gene expression involved in wheat leaf senescence between stay-green and non-stay-green cultivars. Front. Plant Sci. 13:971927. doi: 10.3389/fpls.2022.971927
Received
17 June 2022
Accepted
04 August 2022
Published
26 August 2022
Volume
13 - 2022
Edited by
Tesfamichael Kebrom, Prairie View A&M University, United States
Reviewed by
Ziming Wu, Jiangxi Agricultural University, China; Barbara George-Jaeggli, The University of Queensland, Australia; Thomas Roach, University of Innsbruck, Austria
Updates
Copyright
© 2022 Li, Yang, Guo, Huang, Zhong, Tan, Ren, Li, Chen and Luo.
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: Peigao Luo, lpglab@sicau.edu.cn
†These authors have contributed equally to this work
This article was submitted to Plant Physiology, a section of the journal Frontiers in Plant Science
Disclaimer
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.