Next Article in Journal
One, No One, and One Hundred Thousand: The Many Forms of Ribonucleotides in DNA
Next Article in Special Issue
Genetic and Global Epigenetic Modification, Which Determines the Phenotype of Transgenic Rice?
Previous Article in Journal
Systematic Review of Polygenic Risk Scores for Type 1 and Type 2 Diabetes
Previous Article in Special Issue
Rice OsAAA-ATPase1 is Induced during Blast Infection in a Salicylic Acid-Dependent Manner, and Promotes Blast Fungus Resistance
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJMS
Volume 21
Issue 5
10.3390/ijms21051704
ijms-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A RING-Type E3 Ubiquitin Ligase, OsGW2, Controls Chlorophyll Content and Dark-Induced Senescence in Rice

by
Kyu-Chan Shim
1,
Sun Ha Kim
1,
Yun-A Jeon
1,
Hyun-Sook Lee
1,
Cheryl Adeva
1,
Ju-Won Kang
2,
Hyun-Jung Kim
3,
Thomas H Tai
4,5 and
Sang-Nag Ahn
1,*
1
Department of Agronomy, Chungnam National University, Daejeon 34134, Korea
2
Department of Southern Area Crop Science, National Institute of Crop Science, RDA, Miryang 50424, Korea
3
LG Chemical, Ltd., Seoul 07796, Korea
4
USDA-ARS Crops Pathology and Genetics Research Unit, Davis, CA 95616, USA
5
Department of Plant Sciences, University of California, Davis, CA 95616, USA
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2020, 21(5), 1704; https://doi.org/10.3390/ijms21051704
Submission received: 22 January 2020 / Revised: 25 February 2020 / Accepted: 27 February 2020 / Published: 2 March 2020
(This article belongs to the Special Issue Molecular Research in Rice: Agronomically Important Traits)
Browse Figures
Versions Notes

Abstract

:
Leaf senescence is the final stage of plant development. Many internal and external factors affect the senescence process in rice (Oryza sativa L.). In this study, we identified qCC2, a major quantitative trait locus (QTL) for chlorophyll content using a population derived from an interspecific cross between O. sativa (cv. Hwaseong) and Oryza grandiglumis. The O. grandiglumis allele at qCC2 increased chlorophyll content and delayed senescence. GW2 encoding E3 ubiquitin ligase in the qCC2 region was selected as a candidate for qCC2. To determine if GW2 is allelic to qCC2, a gw2-knockout mutant (gw2-ko) was examined using a dark-induced senescence assay. gw2-ko showed delayed leaf senescence in the dark with down-regulated expression of senescence-associated genes (SAGs) and chlorophyll degradation genes (CDGs). The association of the GW2 genotype with the delayed senescence phenotype was confirmed in an F2 population. RNA-seq analysis was conducted to investigate 30-day-old leaf transcriptome dynamics in Hwaseong and a backcross inbred line—CR2002—under dark treatment. This resulted in the identification of genes involved in phytohormone signaling and associated with senescence. These results suggested that transcriptional regulation was associated with delayed senescence in CR2002, and RING-type E3 ubiquitin ligase GW2 was a positive regulator of leaf senescence in rice.

1. Introduction

Chlorophyll (Chl) is a photosynthetic pigment that is an essential component of the plant photosystem. It changes solar energy to chemical energy. Because of its photosynthetic ability, increasing Chl content in crops may be an effective way to increase grain yield and biomass production [1]. A number of studies have demonstrated that Chl content is controlled by quantitative trait locus (QTL) in various genetic backgrounds in rice [2,3,4,5,6].
Senescence or biological aging is the final stage of plant development [7]. During senescence, leaf color turns from green to yellow because of chlorophyll degradation [8]. Leaf yellowing is frequently used as a senescence indicator. Stay-green (non-yellowing) mutants maintain leaf greenness after the grain-ripening stage [9]. Stay-green mutants exhibit delayed leaf senescence and have been found in various plant species [8,10]. In rice, several stay-green genes have been cloned and characterized, including Stay-Green Rice (SGR), NON-YELLOW COLORING 1 (NYC1), NYC1-LIKE (NOL), Oryza sativa NAC-like, activated by Apetala3/Pistillata (OsNAP) [10,11,12,13,14]. SGR, a highly conserved senescence-associated gene, encodes a novel chloroplast protein, and the expression of SGR is up-regulated in both natural and dark-induced senescence. This gene interacts with the light-harvesting chlorophyll-binding protein (LHCP) [10,14]. NYC1 and NOL encode short-chain dehydrogenase/reductase (SDR) [11,12]. These two genes are co-localized in the thylakoid membrane and function as a chlorophyll b reductase. OsNAP contains a typical NAC structure at the N terminus [13]. This gene plays a role in regulating leaf senescence and acts as a key component, linking ABA signaling in rice. These genes have been used as ideal markers for the onset of the senescence process.
Ubiquitination is the post-translational modification of protein substrates [15]. The ubiquitin-proteasome pathway is known to play an important role in plant seed and organ size determination, such as DA1, DA2, EOD1/BB, SAMBA, and SOD/UBP15 genes in Arabidopsis [16,17,18,19,20]. A major function of E3 ubiquitin ligase is to regulate polyubiquitination and to degrade their target substrate proteins. E3 ubiquitin ligase BigBrother (BB) controls final organ size and seed size acting in parallel with the DA1 gene in Arabidopsis, and these two genes also positively control leaf senescence [17,19,20]. Delayed senescence has been confirmed by measuring the expression of AHK3, CRF6, and ARF2. The negative regulators of senescence—AHK3 and CRF6—are expressed higher in da1-1_bb/eod1-2 leaves, while the auxin repressor ARF2 is expressed at a lower level than wild type [19]. DA2, an ortholog of OsGW2, controls seed and organ size by interacting with DA1 [18].
OsGW2, a QTL on chromosome 2, controls grain width and weight in rice and encodes RING-type E3 ubiquitin ligase [21]. This gene might function in the degradation by the ubiquitin-proteasome pathway, and loss-of-function in OsGW2 increases cell numbers and grain size [21]. The OsGW2 mutant also shows increased transcript levels of OsPCR1 during the developing grain stage, leading to an increase in Zn concentration in the seed [22]. Yeast two-hybrid and in vitro pull-down assays have shown that OsGW2 directly interacts with expansin-like 1 (EXPLA1), chitinase 14 (CHT14), and phosphoglycerate kinase (PGK) [23,24]. Transcription activation activity has been found in the GW2-C terminus (205 to 260) [24]. Together, the findings that OsGW2 interacts with various proteins and that the mutation of OsGW2 has a pleiotropic effect in rice raises the possibility that GW2 is also involved in regulating leaf senescence.
The objective of this study was to identify QTL controlling chlorophyll content and leaf senescence and to characterize candidate genes associated with this trait. qCC2, a major QTL for chlorophyll content, was mapped using an introgression line, CR2002, developed from an interspecific cross between O. sativa cv. Hwaseong and the wild species O. grandiglumis. CR2002 harboring qCC2 from O. grandiglumis showed delayed leaf yellowing and a higher Fv/Fm value than Hwaseong. Endogenous expression levels of senescence-associated genes (SAGs) and chlorophyll degradation genes (CDGs) in CR2002 were lower than in Hwaseong. To determine if GW2 is allelic to qCC2 and thereby associated with delayed senescence, a gw2-knockout mutant (hereafter termed gw2-ko) was examined using dark-induced senescence (DIS) assay. gw2-ko showed delayed leaf senescence under dark conditions with down-regulated expression of SAGs and CDGs. Delayed senescence was confirmed by segregation analysis in the F2 generation. Comparative RNA-seq analysis was conducted to identify differentially-expressed transcripts, using 30-day-old leaves from Hwaseong and CR2002 that were subjected to dark treatment. This enabled the identification of which genotype-specific expressions were enriched in CR2002 and reduced in Hwaseong and vice versa, including genes involved in phytohormone biosynthesis and signaling, NAC transcription factors, and senescence-associated genes. Collectively, the results of this study suggested that OsGW2 controlled chlorophyll content and was a positive regulator of leaf senescence in rice.

2. Results

2.1. QTL Analysis for Chlorophyll Content

CR2002 exhibited a stay-green phenotype and a higher SPAD value (a parameter of leaf greenness) than Hwaseong (Figure 1b). CR2002, which had four O. grandiglumis introgressions, also showed differences in several agronomic traits, including grain weight (Table 1, Figure 1a, and Figure S1). To identify loci associated with chlorophyll content and the stay-green phenotype, QTL analysis was conducted using F3 and F4 populations derived from a cross between Hwaseong and CR2002. Chlorophyll contents were measured at the heading stage and one month after heading. A total of three significant QTLs located on chromosomes 1, 2, and 6 were detected in the F3 and F4 populations (Table 2 and Figure S2). A QTL for chlorophyll content (qCC2) was located on chromosome 2 between RM12813 and RM12983. This QTL was repeatedly detected in F3 and F4 generations at both stages, heading (Chlorophyll I) and one month after heading (Chlorophyll II). qCC2 explained 24.6% of the phenotypic variation, and the O. grandiglumis allele contributed to the increased chlorophyll content, indicating that qCC2 was a major QTL responsible for leaf greenness at heading and ripening stage.
The qCC2 region was confirmed and fine-mapped by substitution mapping (Figure S3). SPAD values were measured four times, from 95 to 125 days after transplanting (DAT), for four substitution lines and two parental lines. CR7036 and CR7039 showed significantly higher SPAD values at 125 DAT than CR7038 and CR7058, suggesting that qCC2 was located in the marker interval between RM3390 and RM7288 (about 1.4 Mbp).

2.2. Dark-Induced Senescence in Hwaseong and CR2002

To determine whether qCC2 is associated with senescence under dark condition, detached leaves of Hwaseong and CR2002 were incubated on 3 mM MES buffer at 27 °C under complete darkness (Figure 2). Hwaseong leaves turned yellow four days after incubation (DAI), while CR2002 remained green (Figure 2a). Although the chlorophyll content of CR2002 was higher than Hwaseong at 0 DAI, a larger difference was observed at 4 DAI between CR2002 and Hwaseong (Figure 2b). The efficiency of photosystem II (Fv/Fm ratio) of CR2002 was significantly higher than Hwaseong, indicating CR2002 showed delayed senescence in the DIS assay (Figure 2c).
To further examine the different senescence of Hwaseong and CR2002, transcript levels of SAGs (OsNAP, Osh36, and OsI57) and CDGs (SGR, RCCR1, PAO, NYC1, and NOL) were measured by qRT-PCR using detached leaf samples (Figure 3). During DIS, transcription levels of the SAGs and CDGs showed increases of 4.9-fold to 325-fold in Hwaseong, whereas increases of 4.7-fold to 200-fold were observed in CR2002. Transcript levels of SAGs and CDGs in CR2002 were significantly lower than Hwaseong at 4 DAI, with the exception of OsI57, indicating that the lower expression of SAGs and CDGs in CR2002 might result in the delayed senescence exhibited by CR2002 in the DIS assay.

2.3. Characterization of gw2-ko

The O. grandiglumis segment harboring qCC2 is approximately 1.4 Mb in size. In this region, GW2 was selected as a candidate gene based on its function as a RING-type E3 ubiquitin ligase. A 1-bp deletion at position 316 in the coding region of the O. grandiglumis GW2 allele, which causes a premature stop and leads to an increased grain width and thickness, was also observed in this study (Figure S4) [21]. In Arabidopsis, ubiquitin-activated peptidase DA1 regulates leaf senescence together with E3 ubiquitin ligase BB [18]. In addition, DA1 controls seed and final organ size with interacting partner DA2, which is orthologous to OsGW2. Due to the previously identified pleiotropic effect of GW2, the GW2 gene was selected as a candidate gene for the qCC2 QTL.
To determine if GW2 is associated with delayed senescence, a T-DNA insertional mutant (PFG_1B-10017.R) was selected from a T-DNA insertional mutant library [25]. The T-DNA insertion was confirmed by two sets of PCR primers, and homozygous T-DNA lines were selected by genomic PCR (Figure 4a,b). The OsGW2 knock-out mutant (gw2-ko) showed increased grain size, and the transcription of OsGW2 in wild type (WT Dongjin) and gw2-ko was evaluated using semi-quantitative RT-PCR (Figure 4c,d). The OsGW2 transcript level was examined during DIS in WT. The expression level was down-regulated at 2 DAI, whereas gene expression was up-regulated after 4 DAI, and a similar level of gene expression with 0 DAI was observed (Figure 4e). Consistent expression patterns of OsGW2 and OgGW2 (O. grandiglumis GW2) were found in Hwaseong and CR2002, but CR2002 showed lower transcript levels than Hwaseong at 0 and 4 DAI (Figure S5).
The gw2-ko mutant showed significantly higher SPAD values than Dongjin from 85 to 130 DAT under the field conditions (Figure 5a,b). To characterize the phenotypic difference between Dongjin and gw2-ko during DIS, fully expanded detached leaves were incubated in 3 mM MES buffer (pH 5.8) at 27 °C in the dark (Figure 5c). Leaves of Dongjin turned yellow at 6 DAI, while gw2-ko remained green. Total chlorophyll contents of Dongjin and gw2-ko at 6 DAI were also significantly different, despite the total chlorophyll contents of gw2-ko at 0 DAI being slightly higher than Dongjin (Figure 5d). The Fv/Fm ratio remained higher in gw2-ko than Dongjin at 5 DAI (Figure 5e). These results demonstrated that leaf senescence was delayed in gw2-ko.
An F2 population (n = 107) derived from a cross between Dongjin and gw2-ko was evaluated for DIS to determine if the T-DNA insertion in OsGW2 was responsible for delayed senescence. After five days of dark incubation, plants homozygous for gw2-ko genotype (KO/KO group) showed significantly higher total chlorophyll contents than those homozygous for Dongjin genotype (DJ/DJ group) (Figure 6a,b). This segregation analysis suggested that T-DNA insertion in OsGW2 caused the delayed senescence exhibited by gw2-ko.
Endogenous expression levels of SAGs and CDGs were examined (Figure 7). Among CDGs, SGR, NYC1, and NOL were significantly down-regulated in gw2-ko at 6 DAI, but PAO and RCCR1 showed no significant difference between Dongjin and gw2-ko. SAGs (OsNAP, Osh36, and OsI57) were also down-regulated in gw2-ko. The loss-of-function of GW2 led to the down-regulation of SAGs and CDGs during DIS.
The physiological phenotypes and endogenous transcript levels during DIS indicated that gw2-ko showed delayed leaf senescence compared to wild type and that GW2 positively regulated leaf senescence in rice.

2.4. Transcriptome Analysis using Hwaseong and CR2002

To gain a better understanding of the mechanism for delayed senescence under dark-treatment conditions, the transcriptome of leaves from Hwaseong and CR2002 plants before (control) and after dark-treatment for 3 days were investigated by RNA-seq. A total of 264 million reads were generated using the Illumina Hiseq 4000 platform, and adapter and low-quality read trimming were conducted (Table S2). A total of 262 million high-quality reads were mapped against the reference genome sequence, and an average mapping rate of 66% was observed (Table S2). A multi-dimensional scaling plot showed that three replicates were clustered, indicating samples were consistently generated (Figure S6).
Differentially expressed genes (DEGs) were identified in Hwaseong and CR2002 under dark treatment compared with control conditions. DEGs were determined based on the criteria of a greater than 2-fold change with significance at Benjamini-Hochberg false discovery rate adjusted p < 0.05. A total of 10,339 DEGs was detected in the comparison of 0 day and 3 days of dark treatment in Hwaseong, and 9885 DEGs were found in CR2002 genotype during the dark treatment (Figure 8a). A Venn diagram of the DEGs indicated that Hwaseong and CR2002 shared 7901 of them (Table S5). Among these 7901 genes, only 29 genes showed different regulation patterns between Hwaseong and CR2002. For example, Os09g0511600 was up-regulated in Hwaseong by 5.13 log2FC (fold-change), while down-regulated in CR2002 by -8.26 log2FC. In addition, a total of 2438 and 1984 genes were specifically identified as DEGs in Hwaseong and CR2002, respectively (Tables S3 and S4). For Hwaseong, 1240 and 1198 genes were specifically up- and down-regulated during the dark treatment, respectively (Figure 8a and Table S3). A total of 838 up-regulated genes and 1146 down-regulated genes were specifically identified during the dark treatment in CR2002 (Figure 8a and Table S4).
To know functional information of the DEGs specifically detected in Hwaseong and CR2002, gene ontology (GO) analysis was performed with agriGO v2.0. For the 1984 CR2002-specific DEGs, a total of 33 GO terms was significantly enriched (Figure 8b). Among 33 GO terms, biological process, cellular component, and molecular function ontologies included 8, 21, and 4 GO terms, respectively. For Hwaseong-specific DEGs, a total of 33 GO terms were also significantly overrepresented. Of these GO terms, 6, 23, and 4 were included in the biological process, cellular component, and molecular function ontologies, respectively (Figure S7). Hwaseong-specific DEGs included intracellular signaling cascade, organelle membrane, cytosol, cytoskeleton, Golgi membrane, endomembrane system, and Golgi apparatus GO terms, while iron ion transport, cytoplasmic part, envelope, cellular response to chemical stimulus, cell cycle, microtubule-associated complex, and organelle envelope were included only in CR2002-specific DEGs.
To understand the mechanism involved in delayed senescence in CR2002, the SAGs, CDGs, and genes involved in phytohormone biosynthesis and signaling were examined from the DEGs. A total of 109 genes were identified (Table S6). Among these genes, we focused on those specifically detected in Hwaseong or CR2002 (Table 3). During the dark treatment, OsNIT1, OsDWF1, and OsRR3 were up-regulated in CR2002, and OsERF5, OsABI1, OsJAZ1, and OsGLU were down-regulated. Phytohormone signaling and biosynthesis genes were mainly found in the CR2002 DEGs. When these seven DEGs were subjected to RT-qPCR using the same G0/G3 samples as for RNA-Seq, similar expression patterns were observed, suggesting the involvement of these DEGs in delayed senescence in CR2002 (Figure S8). For Hwaseong, 18 genes were identified from the DEGs, and 10 and 8 genes were down- and up-regulated, respectively. Hwaseong-specific DEGs included phytohormone biosynthesis and signaling genes, NAC transcription factors, and senescence-associated genes.

3. Discussion

3.1. qCC2 is Associated with Chlorophyll Content and Leaf Senescence

Chlorophyll content, one of the most important traits due to its role in photosynthesis, is regulated by polygenic loci. Several studies have identified a number of QTLs across various genetic resources and environments [3,4,5,6]. A number of studies reported QTLs on the short arm of chromosome 2. Four QTLs have been mapped for stay-green traits on chromosome 2 between RM145 and RM322 using 190 doubled haploid lines from a cross between ‘Zhenshan 97′ and ‘Wuyujing 2′ [2]. Jiang et al. (2012) reported QTL for chlorophyll content from heading to the maturity stage (qCHM2) located on the short arm of chromosome 2 between RM145 and RM29 [5]. In this study, three QTLs for chlorophyll content were detected using progeny from an interspecific cross between O. sativa and O. grandiglumis. qCC2 was located on chromosome 2 between RM3390 and RM7288, which were approximately 1.4 Mbp apart. This QTL was detected in F3 and F4 generations and at two different developmental stages, respectively, and confirmed by substitution mapping (Figure S3). qCC2 shared a similar location with QTLs reported in previous studies [2,5]. Although allelism tests are needed to clarify the relationships among these QTLs, these studies indicated the presence of genes associated with chlorophyll content and leaf senescence on the short arm region of chromosome 2.
When Hwaseong and CR2002 were tested for dark-induced senescence, CR2002 showed higher chlorophyll content with delayed leaf senescence than Hwaseong (Figure 2a,b). Fv/Fm values were also significantly different at 4 days after incubation (Figure 2c). Endogenous expression levels of SAGs and CDGs were severely down-regulated in CR2002 than Hwaseong, and these results indicated that qCC2 was associated with leaf senescence in CR2002.

3.2. GW2 Positively Regulates Leaf Senescence

GW2 encodes a RING-type E3 ubiquitin ligase and resides in the qCC2 region (1.4 Mbp). The O. grandiglumis GW2 allele in CR2002 has the 1-bp deletion in 4th exon, which causes a premature stop and truncated protein (Figure S4) [21]. Many studies have suggested ubiquitination as a candidate pathway for the regulation of senescence. SAUL1 encodes E3 ubiquitin ligase, which is required for suppression of premature senescence in Arabidopsis [26]. Also, spl11 encodes E3 ubiquitin ligase, which regulates cell death in rice [27]. E3 ubiquitin ligase Big Brother (BB) gene and ubiquitin receptor DA1 gene negatively regulate leaf size and promote senescence [19]. Mutations in DA1 and BB enhance leaf growth, an effect that is synergistically increased in the double mutant. A da1-1/eod1-2 double mutant especially exhibits a longer lifespan than wild type or the single mutants [19]. In addition, DA1 physically interacts with DA2, an ortholog of GW2, and DA2 acts synergistically with DA1 to regulate seed size in Arabidopsis [18]. Based on these previous results, GW2 was chosen as a candidate gene for the delayed senescence phenotype in CR2002. The transcript level of GW2 was down-regulated at 2 DAI and showed a similar expression level at 0 DAI and at 4 and 6 DAI in the dark (Figure 4e and Figure S5). This expression pattern was similar to the NYC4 gene, which functions in the degradation of chlorophyll and chlorophyll-protein complexes during DIS [28]. The effect of gw2-ko was examined in DIS (Figure 5c). gw2-ko showed delayed senescence and down-regulated expression of SAGs and CDGs. Compared to Hwaseong and Dongjin, most of the SAGs and CDGs genes were consistently down-regulated in CR2002 and gw2-ko during DIS. However, transcript levels of RCCR1 (red chlorophyll catabolite reductase 1), PAO (Pheophorbide a oxygenase), and OsI57 (putative 3-ketoacyl-CoA thiolase) did not show significant differences between Dongjin and gw2-ko or Hwaseong and CR2002 in DIS. These results suggested that RCCR1 and PAO might function down-stream in the chlorophyll degradation pathway [8]. Because the transcript levels of RCCR1 and PAO in gw2-ko were not significantly different from wild-type Dongjin, GW2 might affect genes up-stream in the chlorophyll degradation pathway.

3.3. Differentially Expressed Genes Associated with Leaf Senescence

Comparative RNA-seq analysis was conducted to identify differentially-expressed transcripts, using 30-day-old leaves from Hwaseong and CR2002 to gain a better understanding of the mechanism involved in delayed senescence. A total of 10,339 and 9885 DEGs were identified in Hwaseong and CR2002 under dark treatment compared with the control, respectively (Figure 8a). To examine the difference between Hwaseong and CR2002, Hwaseong- and CR2002-specific DEGs were extracted. In CR2002-specific DEGs, 33 GO terms were enriched, and these included iron ion transport, response to hormone stimulus, cell cycle, response to endogenous stimulus, and chloroplast. Transporter genes are closely associated with leaf senescence process. During senescence, remobilization of nutrients from senescing cells to developing tissues is mediated by transporters, and transporter genes show notable changes in expression [13,29]. Chloroplast GO term contains 15 CR2002-specific DEGs, including Os07g0462000 (OsSG1), Os07g0558500 (NYC4), Os02g0152400 (RbcS1) [28,30]. Os07g0462000 encodes glutamate-cysteine ligase and possibly controls chlorophyll content and stay-green phenotype [30]. Os07g0558500 (NYC4) is an ortholog of THF1 in Arabidopsis. The AtTHF1 is known as a multi-function protein involved in acclimation to high light, sugar sensing, and disease resistance, while nyc4-1 mutant shows the stay-green phenotype in rice [28]. Os02g0152400 encodes the rubisco small subunit 1, which reflects involvement in the regulation of rubisco catalytic activity.
Transcription factors (TFs) play an important role in regulating leaf senescence. NAC (NAM, ATAF1, and CUC2) transcription factor family is one of the large TF families, and many NAC TFs are up-regulated during leaf senescence [31]. In addition, several NAC TFs regulating leaf senescence have been characterized in various plant species [13,32,33,34,35,36]. Transcript levels of NAC transcription factors (OsNAP, ONAC120, and ONAC122) were specifically up-regulated in Hwaseong (Table 3). OsNAP is a positive regulator of leaf senescence, and its orthologs TtNAM-B1 and AtNAP have shown similar function in wheat and Arabidopsis [13,32,33]. OsNAC122 (or OsNAC10, Os11g0126900) over-expression and root-specific expression transgenic have shown significantly improved drought, high salinity, and low-temperature tolerance in rice [37].
Phytohormone signaling and biosynthesis-related genes were identified in Hwaseong- and CR2002-specific DEGs. During dark treatment, ethylene (OsERF5), ABA (OsABI1), and jasmonic acid (OsJAZ1) signaling genes were significantly down-regulated, while auxin synthesis (OsNIT2), brassinosteroid biosynthesis (OsDWF1), and cytokinin signaling (OsRR3) genes were up-regulated in CR2002 (Table 3). In contrast, ethylene biosynthesis and signaling (OsACS2 and OsERF2) and auxin signaling (OsARF2) genes were up-regulated in Hwaseong, whereas cytokinin signaling (OsRR6, OsRR9, and OsRR10), brassinosteroid signaling (OsBRI1), and gibberellin 20-oxidase (OsGA20ox1) were down-regulated (Table 3). These results indicated that delayed leaf senescence in CR2002 (gw2-ko) was associated with phytohormone signaling or biosynthesis pathways in rice. The rice stay-green mutant oself3.1 has shown down-regulated transcript levels of ABA-, ethylene-, JA-associated genes in microarray analysis [7]. The results from the present study were consistent with the report that OsELF3.1 promoted leaf senescence by modulating signaling pathways [7]. To further investigate this relationship with phytohormones, gene expression of GW2 in response to plant hormones was examined using the RiceXPro database. GW2 did not show large changes in response to plant hormones (data not shown) [38]. In durum wheat, a GW2 knock-down mutant has shown enhanced cytokinin dehydrogenase 1 (CKX1) transcript level and down-regulation of cytokinin dehydrogenase 2 (CKX2) and gibberellin oxidase 3 (GA3-ox) [39]. In addition, a wheat near-isogenic line (NIL) harboring the TaGW2-6A allelic variant has shown an increase in not only grain size but also endogenous cytokinin content (Z+ZR) [40]. Compared to the control (Chinese Spring), the expression level of cytokinin biosynthesis genes (TaIPT2, TaIPT3, TaIPT5, and TaIPT8) are found to be up-regulated, while cytokinin degradation genes (TaCKX1, TaCKX2, and TaCKX6) are found to be down-regulated in the NIL during the endosperm development [40]. These results supported the possibility that delayed senescence in CR2002 and gw2-ko was mediated by phytohormone-related pathways.

3.4. GW2 as an E3 Ubiquitin Ligase

The ubiquitination process is a post-translational modification. E3 ligases function in the regulation of polyubiquitination and are involved in target protein degradation. Several studies have reported interacting partners of GW2. Target proteins of GW2 have been revealed through yeast two-hybrid screening [24]. Yeast two-hybrid and pull-down assays have shown that expansin-like 1 (EXPLA1) directly interacts with GW2. In addition, in vitro assays have identified the ubiquitination of EXPLA1 by GW2 at lysine 279 (K279). Proteomic analysis has shown that chitinase 14 (CHT14) and phosphoglycerate kinase (PGK) directly interact with GW2 [23]. However, genes that are directly associated with leaf senescence have not been reported. Identification of GW2 substrate protein(s) responsible for leaf senescence would be key to understanding the post-translational molecular mechanism of leaf senescence in rice.

4. Materials and Methods

4.1. Plant Materials

CR2002 was derived from a cross between O. sativa Korean elite line ‘Hwaseong’ as a recurrent parent and wild species O. grandiglumis (2n = 48, CCDD, Acc. No. 101154) [41]. Four O. grandiglumis chromosome segments were found in CR2002 on chromosomes 1, 2, 6, and 10. CR2002 was backcrossed with Hwaseong to produce an F2 population. A total of 705 F2 plants were genotyped using SSR markers from the introgression regions. Fifty-eight F3 plants with different combinations of O. grandiglumis homozygous segments were selected and advanced to the F3. Among the 58 F3 lines, 17 that have different recombination break-points within the qCC2 region were selected and advanced to the F4 generation. QTL analysis was performed using two replications of the 58 F3 lines and 17 F4 lines, which were grown in the Chungnam National University paddy field during 2015 and 2016, respectively. For substitution mapping, four F4 recombinant lines in the qCC2 region were investigated in 2017. To determine whether OsGW2 is associated with leaf senescence, a gw2-knockout (gw2-ko) mutant (PFG_1B-10017.R) in the japonica rice ‘Dongjin’ background was used for evaluating leaf senescence [42,43]. An F2 population (n = 107) derived from a cross between Dongjin and gw2-ko was generated and used for DIS assay to confirm whether the T-DNA insertion in OsGW2 locus was responsible for delayed senescence.

4.2. Phenotypic Evaluation

A total of 10 traits, including chlorophyll content and three-grain morphology traits (Table 1), were evaluated in the parental lines (n= 10) using the methods, as described in Yoon et al. (2006) [41].
To evaluate chlorophyll content, three methods were used. First, fluorescence ratio F735/F700 was measured in F3 and F4 generation plants, and the relative chlorophyll content (mg/m2) was evaluated using CCM-300 (OPTI-SCIENCE, Hudson, NH, USA). This data was used for QTL analysis. Second, SPAD values were measured to compare parental lines (Hwaseong, CR2002, Dongjin, and gw2-ko) using a SPAD-502Plus (KONICA MINOLTA, Tokyo, Japan). Third, total chlorophyll was extracted from detached leaves and used for dark-induced senescence using pre-chilled 80% acetone [7]. The absorbance of supernatants was measured at 645 and 663 nm using a UV/VIS spectrophotometer (Hanson Tech., Seoul, Korea) [44]. For the dark-induced senescence (DIS) assay, fully expanded flag leaves were detached and placed in 3 mM MES buffer (pH 5.8) and incubated at 27 °C under complete darkness [45]. Fv/Fm ratio was measured using a Hansatech PEA MK2 (Hansatech, Norfolk, UK). The middle part of leaf samples from DIS was adapted in the dark for 30 min to complete the oxidation of QA, and then the Fv/Fm value was measured. A student’s t-test was performed for trait comparisons.

4.3. DNA Extraction and QTL Analysis

DNA was extracted from leaf tissues using the extraction method, as described in Causse et al. (1994) [46]. Genotyping was conducted in 58 F3 and 17 F4 lines using 11 polymorphic SSR markers between CR2002 and Hwaseong. PCR was conducted, as described in Shim et al. (2019) [47]. PCR products were separated on 3% metaphor agar stained with StaySafe Nucleic Acid Gel Stain (RBC, Taiwan) or 4% polyacrylamide denaturing gel stained with Silver Staining Kit (Bioneer, Korea), respectively.
QTLs were determined by single-marker analysis (SMA). SMA was conducted to establish the effect of each marker on each trait. In SMA, a QTL was declared when the phenotype was associated with SSR genotype at p < 0.05 by one-way analysis of variance (ANOVA). The observed phenotypic variation was estimated using the coefficient of determination (R2). For Tukey’s test, Minitab19 software was used. Student’s t-test was conducted using Microsoft Excel. Boxplots were drawn using R version 3.5.3.

4.4. RNA Isolation and Quantitative Real-Time PCR

Total RNA of detached leaves from DIS was isolated using NucleoSpin RNA (Macherey Nagel, Deuren, Germany), according to the manufacturer’s instructions. Following reverse-transcription into the first-strand cDNA with SmartGene Mixed cDNA synthesis kit (SJ Bioscience, Daejeon, Korea), real-time PCR was performed using a CFX Connect Real-Time System (Bio-Rad, Hercules, CA, USA). Each reaction contained 10 µL of 2× SYBR Green Master Mix reagent, 2 µL of diluted cDNA samples, 1 µL of 10 pmol gene-specific primers in a final volume of 20 µL. The real-time PCR protocol was: 15 min at 95 °C to denature and activate an enzyme, followed by 40 cycles at 95 °C for 20 s (denaturation), 60 °C for 40 s (annealing), 72 °C for 30 s (extension). qRT-PCR was conducted according to the 2-∆∆Ct method [48]. Rice UBIQUITIN5 (OsUBQ5) was used for normalization, and relative expression levels were compared by Student’s t-test. Primers used in this study are listed in Table S1.

4.5. RNA-Seq Analysis

The 30-day-old Hwaseong and CR2002 plants were dark-treated to induce senescence for three days. From plants each before (day 0) and 3 days after dark-treatment (day 3), the 5th leaves were collected with three replicates, and total RNA was extracted, as described above. A total of 12 transcriptome libraries were constructed using TruSeq RNA Sample Prep Kit v2 (Illumina, San Diego, CA, USA), according to the manufacturer’s instructions. The library construction was conducted by GnC Bio Inc. (Daejeon, Korea). The constructed libraries were sequenced on a HiSeq 4000 instrument with 150 bp of paired-end reads. The sequence data generated from this study have been deposited at the National Center for Biotechnology Information Sequence Read Archive (SRA) under accession number PRJNA602373. Sequencing results were checked using FastQC (version 0.11.7). Raw sequences were processed to obtain clean reads by removing adapters and low-quality reads, using SCYTHE (https://github.com/vsbuffalo/scythe) and SICKLE (https://github.com/najoshi/sickle), respectively. High-quality reads were mapped on rice reference genome IRGSP-1.0, and transcript abundances were quantified using Salmon version 0.8.2 [49]. Differentially expressed genes (DEGs) were determined with more than 2-fold change with significance at Benjamini-Hochberg false discovery rate adjusted p < 0.05 using Limma-Voom package [50,51]. GO enrichment was carried out using agriGO v2.0 singular enrichment analysis [52].

5. Conclusions

A mapping population derived from an interspecific cross between O. sativa and O. grandiglumis facilitated the discovery of qCC2, a major QTL for chlorophyll content. The O. grandiglumis allele at qCC2 increased chlorophyll content and delayed senescence. Building on several previous studies, GW2, which encodes E3 ubiquitin ligase, was selected as the most probable candidate gene for this phenotype, and this was confirmed through genetic analysis of a knockout mutant. The gw2-ko mutant showed delayed leaf senescence under dark conditions with down-regulation of SAGs and CDGs, and the delayed senescence trait completely segregated with the gw2-ko gene. On the basis of this work, the RING-type ubiquitin ligase GW2 appeared to function as a positive regulator of leaf senescence in rice. This study also represented the first investigation of global gene expression in rice leaves subjected to complete darkness. The results of this whole-genome transcriptome analysis of Hwaseong and CR2002 highlighted the processes and genes that are likely to play a role in the mechanisms underlying leaf senescence in rice (e.g., carbohydrate metabolism, pyrophosphate-dependent energy conservation, and ethylene signaling). The QTL analysis and the genotype-specific regulation suggested that senescence was regulated at the transcriptional level, although the possibility of post-translational regulation by GW2 requires further investigation.

Supplementary Materials

Supplementary materials can be found at https://www.mdpi.com/1422-0067/21/5/1704/s1.

Author Contributions

K.-C.S. and S.H.K. performed experiments and analyzed data. K.-C.S. and S.-N.A. conceived the study, designed, and supervised the study. Y.-A.J., H.-S.L., J.-W.K., and C.A. investigated the traits. H.-J.K. and T.H.T. assisted in analyzing the RNA-seq data. K.-C.S., T.H.T., and S.-N.A. wrote and edited the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This work was carried out with the support of “Cooperative Research Program for Agriculture Science and Technology Development (Project No. PJ01321401)”, Rural Development Administration and the Basic Science Research Program through the National Research Foundation (NRF) of Korea funded by the Ministry of Education (NRF-2017R1A2B2007554).

Acknowledgments

The gw2-knockout (gw2-ko) mutant seeds were provided by G. An, Kyung Hee University. We thank Cheol Seong Jang and Ju Hee Kim, Kangwon National University, for their critical comments on this work and Young-Sook Kim for her technical assistance.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

Abbreviations

DISDark-induced senescence
SNPSingle nucleotide polymorphism
InDelInsertion and deletion
qRT-PCRQuantitative real-time PCR
SSRSimple sequence repeat
CDGChlorophyll degradation gens
QTLQuantitative trait loci
SAGSenescence-associated gene
WTWild type

References

  1. Wang, F.H.; Wang, G.X.; Li, X.Y.; Huang, J.L.; Zheng, J.K. Heredity, physiology and mapping of a chlorophyll content gene of rice (Oryza sativa L.). J. Plant Physiol. 2008, 165, 324–330. [Google Scholar] [CrossRef] [PubMed]
  2. Jiang, G.H.; He, Y.Q.; Xu, C.G.; Li, X.H.; Zhang, Q. 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. 2004, 108, 688–698. [Google Scholar] [CrossRef] [PubMed]
  3. Teng, S.; Qian, Q.; Zeng, D.L.; Kunihiro, Y.; Fujimoto, K.; Huang, D.N.; Zhu, L.H. QTL analysis of leaf photosynthetic rate and related physiological traits in rice (Oryza sativa L.). Euphytica 2004, 135, 1–7. [Google Scholar] [CrossRef]
  4. Takai, T.; Kondo, M.; Yano, M.; Yamamoto, T. A Quantitative Trait Locus for Chlorophyll Content and its Association with Leaf Photosynthesis in Rice. Rice 2010, 3, 172–180. [Google Scholar] [CrossRef] [Green Version]
  5. Jiang, S.K.; Zhang, X.J.; Zhang, F.M.; Xu, Z.J.; Chen, W.F.; Li, Y.H. Identification and Fine Mapping of qCTH4, a Quantitative Trait Loci Controlling the Chlorophyll Content from Tillering to Heading in Rice (Oryza sativa L.). J. Hered. 2012, 103, 720–726. [Google Scholar] [CrossRef] [Green Version]
  6. Abdelkhalik, A.F.; Shishido, R.; Nomura, K.; Ikehashi, H. QTL-based analysis of leaf senescence in an indica/japonica hybrid in rice (Oryza sativa L.). Theor. Appl. Genet. 2005, 110, 1226–1235. [Google Scholar] [CrossRef]
  7. Sakuraba, Y.; Han, S.H.; Yang, H.J.; Piao, W.L.; Paek, N.C. Mutation of Rice Early Flowering3.1 (OsELF3.1) delays leaf senescence in rice. Plant Mol. Biol. 2016, 92, 223–234. [Google Scholar] [CrossRef]
  8. Hortensteiner, S. Stay-green regulates chlorophyll and chlorophyll-binding protein degradation during senescence. Trends Plant Sci. 2009, 14, 155–162. [Google Scholar] [CrossRef]
  9. Thomas, H.; Smart, C.M. Crops That Stay Green. Ann. Appl. Biol. 1993, 123, 193–219. [Google Scholar] [CrossRef]
  10. Park, S.Y.; Yu, J.W.; Park, J.S.; Li, J.; Yoo, S.C.; Lee, N.Y.; Lee, S.K.; Jeong, S.W.; Seo, H.S.; Koh, H.J.; et al. The senescence-induced staygreen protein regulates chlorophyll degradation. Plant Cell 2007, 19, 1649–1664. [Google Scholar] [CrossRef] [Green Version]
  11. Kusaba, M.; Ito, H.; Morita, R.; Iida, S.; Sato, Y.; Fujimoto, M.; Kawasaki, S.; Tanaka, R.; Hirochika, H.; Nishimura, M.; et al. Rice NON-YELLOW COLORING1 is involved in light-harvesting complex II and grana degradation during leaf senescence. Plant Cell 2007, 19, 1362–1375. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Sato, Y.; Morita, R.; Katsuma, S.; Nishimura, M.; Tanaka, A.; Kusaba, M. Two short-chain dehydrogenase/reductases, NON-YELLOW COLORING 1 and NYC1-LIKE, are required for chlorophyll b and light-harvesting complex II degradation during senescence in rice. Plant J. 2009, 57, 120–131. [Google Scholar] [CrossRef] [PubMed]
  13. Liang, C.Z.; Wang, Y.Q.; Zhu, Y.N.; Tang, J.Y.; Hu, B.; Liu, L.C.; Ou, S.J.; Wu, H.K.; Sun, X.H.; Chu, J.F.; et al. OsNAP connects abscisic acid and leaf senescence by fine-tuning abscisic acid biosynthesis and directly targeting senescence-associated genes in rice. Proc. Natl. Acad. Sci. USA 2014, 111, 10013–10018. [Google Scholar] [CrossRef] [Green Version]
  14. Jiang, H.W.; Li, M.R.; Liang, N.B.; Yan, H.B.; Wei, Y.L.; Xu, X.; Liu, J.F.; Xu, Z.; Chen, F.; Wu, G.J. Molecular cloning and function analysis of the stay green gene in rice. Plant J. 2007, 52, 197–209. [Google Scholar] [CrossRef] [PubMed]
  15. Ardley, H.C.; Robinson, P.A. E3 ubiquitin ligases. Essays Biochem. 2005, 41, 15–30. [Google Scholar] [CrossRef]
  16. Li, N.; Li, Y.H. Signaling pathways of seed size control in plants. Curr. Opin. Plant Biol. 2016, 33, 23–32. [Google Scholar] [CrossRef]
  17. Li, Y.H.; Zheng, L.Y.; Corke, F.; Smith, C.; Bevan, M.W. Control of final seed and organ size by the DA1 gene family in Arabidopsis thaliana. Genes Dev. 2008, 22, 1331–1336. [Google Scholar] [CrossRef] [Green Version]
  18. Xia, T.; Li, N.; Dumenil, J.; Li, J.; Kamenski, A.; Bevan, M.W.; Gao, F.; Li, Y.H. The Ubiquitin Receptor DA1 Interacts with the E3 Ubiquitin Ligase DA2 to Regulate Seed and Organ Size in Arabidopsis. Plant Cell 2013, 25, 3347–3359. [Google Scholar] [CrossRef] [Green Version]
  19. Vanhaeren, H.; Nam, Y.J.; De Milde, L.; Chae, E.; Storme, V.; Weigel, D.; Gonzalez, N.; Inze, D. Forever Young: The Role of Ubiquitin Receptor DA1 and E3 Ligase BIG BROTHER in Controlling Leaf Growth and Development. Plant Physiol. 2017, 173, 1269–1282. [Google Scholar] [CrossRef] [Green Version]
  20. Disch, S.; Anastasiou, E.; Sharma, V.K.; Laux, T.; Fletcher, J.C.; Lenhard, M. The E3 ubiquitin ligase BIG BROTHER controls Arabidopsis organ size in a dosage-dependent manner. Curr. Biol. 2006, 16, 272–279. [Google Scholar] [CrossRef] [Green Version]
  21. Song, X.J.; Huang, W.; Shi, M.; Zhu, M.Z.; Lin, H.X. A QTL for rice grain width and weight encodes a previously unknown RING-type E3 ubiquitin ligase. Nat. Genet. 2007, 39, 623–630. [Google Scholar] [CrossRef]
  22. Song, W.Y.; Lee, H.S.; Jin, S.R.; Ko, D.; Martinoia, E.; Lee, Y.; An, G.; Ahn, S.N. Rice PCR1 influences grain weight and Zn accumulation in grains. Plant Cell Environ. 2015, 38, 2327–2339. [Google Scholar] [CrossRef]
  23. Lee, K.H.; Park, S.W.; Kim, Y.J.; Koo, Y.J.; Song, J.T.; Seo, H.S. Grain width 2 (GW2) and its interacting proteins regulate seed development in rice (Oryza sativa L.). Bot. Stud. 2018, 59, 1–7. [Google Scholar] [CrossRef] [Green Version]
  24. Choi, B.S.; Kim, Y.J.; Markkandan, K.; Koo, Y.J.; Song, J.T.; Seo, H.S. GW2 Functions as an E3 Ubiquitin Ligase for Rice Expansin-Like 1. Int. J. Mol. Sci. 2018, 19. [Google Scholar] [CrossRef] [Green Version]
  25. An, S.Y.; Park, S.; Jeong, D.H.; Lee, D.Y.; Kang, H.G.; Yu, J.H.; Hur, J.; Kim, S.R.; Kim, Y.H.; Lee, M.; et al. Generation and analysis of end sequence database for T-DNA tagging lines in rice. Plant Physiol. 2003, 133, 2040–2047. [Google Scholar] [CrossRef] [Green Version]
  26. Raab, S.; Drechsel, G.; Zarepour, M.; Hartung, W.; Koshiba, T.; Bittner, F.; Hoth, S. Identification of a novel E3 ubiquitin ligase that is required for suppression of premature senescence in Arabidopsis. Plant J. 2009, 59, 39–51. [Google Scholar] [CrossRef]
  27. Zeng, L.R.; Qu, S.H.; Bordeos, A.; Yang, C.W.; Baraoidan, M.; Yan, H.Y.; Xie, Q.; Nahm, B.H.; Leung, H.; Wang, G.L. Spotted leaf11, a negative regulator of plant cell death and defense, encodes a U-box/armadillo repeat protein endowed with E3 ubiquitin ligase activity. Plant Cell 2004, 16, 2795–2808. [Google Scholar] [CrossRef] [Green Version]
  28. Yamatani, H.; Sato, Y.; Masuda, Y.; Kato, Y.; Morita, R.; Fukunaga, K.; Nagamura, Y.; Nishimura, M.; Sakamoto, W.; Tanaka, A.; et al. NYC4, the rice ortholog of Arabidopsis THF1, is involved in the degradation of chlorophyll protein complexes during leaf senescence. Plant J. 2013, 74, 652–662. [Google Scholar] [CrossRef]
  29. Lee, S.; Jeong, H.; Lee, S.; Lee, J.; Kim, S.J.; Park, J.W.; Woo, H.R.; Lim, P.O.; An, G.; Nam, H.G.; et al. Molecular bases for differential aging programs between flag and second leaves during grain-filling in rice. Sci. Rep. 2017, 7, 1–16. [Google Scholar] [CrossRef] [Green Version]
  30. Zhao, Y.; Qiang, C.G.; Wang, X.Q.; Chen, Y.F.; Deng, J.Q.; Jiang, C.H.; Sun, X.M.; Chen, H.Y.; Li, J.; Piao, W.L.; et al. New alleles for chlorophyll content and stay-green traits revealed by a genome wide association study in rice (Oryza sativa). Sci. Rep. 2019, 9, 1–11. [Google Scholar] [CrossRef] [Green Version]
  31. Gregersen, P.L.; Holm, P.B. Transcriptome analysis of senescence in the flag leaf of wheat (Triticum aestivum L.). Plant Biotechnol. J. 2007, 5, 192–206. [Google Scholar] [CrossRef]
  32. Uauy, C.; Distelfeld, A.; Fahima, T.; Blechl, A.; Dubcovsky, J. A NAC gene regulating senescence improves grain protein, zinc, and iron content in wheat. Science 2006, 314, 1298–1301. [Google Scholar] [CrossRef] [Green Version]
  33. Guo, Y.F.; Gan, S.S. AtNAP, a NAC family transcription factor, has an important role in leaf senescence. Plant J. 2006, 46, 601–612. [Google Scholar] [CrossRef]
  34. Zhou, Y.; Huang, W.F.; Liu, L.; Chen, T.Y.; Zhou, F.; Lin, Y.J. Identification and functional characterization of a rice NAC gene involved in the regulation of leaf senescence. BMC Plant Biol. 2013, 13. [Google Scholar] [CrossRef] [Green Version]
  35. Mao, C.J.; Lu, S.C.; Lv, B.; Zhang, B.; Shen, J.B.; He, J.M.; Luo, L.Q.; Xi, D.D.; Chen, X.; Ming, F. A Rice NAC Transcription Factor Promotes Leaf Senescence via ABA Biosynthesis. Plant Physiol. 2017, 174, 1747–1763. [Google Scholar] [CrossRef] [Green Version]
  36. Sakuraba, Y.; Piao, W.; Lim, J.H.; Han, S.H.; Kim, Y.S.; An, G.; Paek, N.C. Rice ONAC106 Inhibits Leaf Senescence and Increases Salt Tolerance and Tiller Angle. Plant Cell Physiol. 2015, 56, 2325–2339. [Google Scholar] [CrossRef]
  37. Jeong, J.S.; Kim, Y.S.; Baek, K.H.; Jung, H.; Ha, S.H.; Do Choi, Y.; Kim, M.; Reuzeau, C.; Kim, J.K. Root-Specific Expression of OsNAC10 Improves Drought Tolerance and Grain Yield in Rice under Field Drought Conditions. Plant Physiol. 2010, 153, 185–197. [Google Scholar] [CrossRef] [Green Version]
  38. Sato, Y.; Takehisa, H.; Kamatsuki, K.; Minami, H.; Namiki, N.; Ikawa, H.; Ohyanagi, H.; Sugimoto, K.; Antonio, B.A.; Nagamura, Y. RiceXPro Version 3.0: Expanding the informatics resource for rice transcriptome. Nucleic Acids Res. 2013, 41, D1206–D1213. [Google Scholar] [CrossRef] [Green Version]
  39. Sestili, F.; Pagliarello, R.; Zega, A.; Saletti, R.; Pucci, A.; Botticella, E.; Masci, S.; Tundo, S.; Moscetti, I.; Foti, S.; et al. Enhancing grain size in durum wheat using RNAi to knockdown GW2 genes. Theor. Appl. Genet. 2019, 132, 419–429. [Google Scholar] [CrossRef]
  40. Geng, J.; Li, L.Q.; Lv, Q.; Zhao, Y.; Liu, Y.; Zhang, L.; Li, X.J. TaGW2-6A allelic variation contributes to grain size possibly by regulating the expression of cytokinins and starch-related genes in wheat. Planta 2017, 246, 1153–1163. [Google Scholar] [CrossRef]
  41. Yoon, D.B.; Kang, K.H.; Kim, H.J.; Ju, H.G.; Kwon, S.J.; Suh, J.P.; Jeong, O.Y.; Ahn, S.N. Mapping quantitative trait loci for yield components and morphological traits in an advanced backcross population between Oryza grandiglumis and the O. sativa japonica cultivar Hwaseongbyeo. Theor. Appl. Genet. 2006, 112, 1052–1062. [Google Scholar] [CrossRef]
  42. Jeon, J.S.; Lee, S.; Jung, K.H.; Jun, S.H.; Jeong, D.H.; Lee, J.; Kim, C.; Jang, S.; Lee, S.; Yang, K.; et al. T-DNA insertional mutagenesis for functional genomics in rice. Plant J. 2000, 22, 561–570. [Google Scholar] [CrossRef]
  43. Jeong, D.H.; An, S.; Park, S.; Kang, H.G.; Park, G.G.; Kim, S.R.; Sim, J.; Kim, Y.O.; Kim, M.K.; Kim, S.R.; et al. Generation of a flanking sequence-tag database for activation-tagging lines in japonica rice. Plant J. 2006, 45, 123–132. [Google Scholar] [CrossRef]
  44. Wu, G.; Lee, C.-H. Chlorophyll a/b Ratio as a Criterion for the Reliability of Absorbance Values Measured for the Determination of Chlorophyll Concentration. J. Life Sci. 2019, 29, 509–513. [Google Scholar]
  45. Oh, M.H.; Kim, Y.J.; Lee, C.W. Leaf senescence in a stay-green mutant of Arabidopsis thaliana: Disassembly process of photosystem I and II during dark-incubation. J. Biochem. Mol. Biol. 2000, 33, 256–262. [Google Scholar]
  46. Causse, M.A.; Fulton, T.M.; Cho, Y.G.; Ahn, S.N.; Chunwongse, J.; Wu, K.S.; Xiao, J.H.; Yu, Z.H.; Ronald, P.C.; Harrington, S.E.; et al. Saturated Molecular Map of the Rice Genome Based on an Interspecific Backcross Population. Genetics 1994, 138, 1251–1274. [Google Scholar]
  47. Shim, K.C.; Kim, S.; Le, A.Q.; Lee, H.S.; Adeva, C.; Jeon, Y.A.; Luong, N.H.; Kim, W.J.; Akhtamov, M.; Ahn, S.N. Fine Mapping of a Low-Temperature Germinability QTL qLTG1 Using Introgression Lines Derived from Oryza rufipogon. Plant Breed. Biotechnol. 2019, 7, 141–150. [Google Scholar] [CrossRef]
  48. Livak, K.J.; Schmittgen, T.D. Analysis of relative gene expression data using real-time quantitative PCR and the 2(T)(-Delta Delta C) method. Methods 2001, 25, 402–408. [Google Scholar] [CrossRef]
  49. Patro, R.; Duggal, G.; Love, M.I.; Irizarry, R.A.; Kingsford, C. Salmon provides fast and bias-aware quantification of transcript expression. Nat. Methods 2017, 14, 417. [Google Scholar] [CrossRef] [Green Version]
  50. Law, C.W.; Chen, Y.S.; Shi, W.; Smyth, G.K. voom: Precision weights unlock linear model analysis tools for RNA-seq read counts. Genome Biol. 2014, 15, R29. [Google Scholar] [CrossRef] [Green Version]
  51. Ritchie, M.E.; Phipson, B.; Wu, D.; Hu, Y.F.; Law, C.W.; Shi, W.; Smyth, G.K. limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Res. 2015, 43, e47. [Google Scholar] [CrossRef]
  52. Tian, T.; Liu, Y.; Yan, H.Y.; You, Q.; Yi, X.; Du, Z.; Xu, W.Y.; Su, Z. agriGO v2.0: A GO analysis toolkit for the agricultural community, 2017 update. Nucleic Acids Res. 2017, 45, W122–W129. [Google Scholar] [CrossRef]
Ijms 21 01704 g001 550
Figure 1. CR2002 showed a higher SPAD value and delayed senescence than Hwaseong in the field condition. Comparison of (a) plant morphology and (b) SPAD value (n = 30) of two parents. Error bars indicate the standard error. DAT: days after transplanting.
Figure 1. CR2002 showed a higher SPAD value and delayed senescence than Hwaseong in the field condition. Comparison of (a) plant morphology and (b) SPAD value (n = 30) of two parents. Error bars indicate the standard error. DAT: days after transplanting.
Ijms 21 01704 g001
Ijms 21 01704 g002 550
Figure 2. CR2002 showed delayed senescence under dark-induced senescence (DIS) conditions. Hwaseong and CR2002 were grown in a paddy field, and fully expanded flag leaves at the heading stage were used for DIS. (a) Detached leaves were incubated in 3 mM MES buffer (pH 5.8) at 27 °C under dark conditions. (b,c) Total chlorophyll contents (n= 6) and Fv/Fm ratio (n= 5) were compared between Hwaseong and CR2002. Error bars indicate a standard error, and more than three samples were used for each experiment. * and *** indicate significant difference at p < 0.05 and p < 0.001 based on Student’s t-test, respectively.
Figure 2. CR2002 showed delayed senescence under dark-induced senescence (DIS) conditions. Hwaseong and CR2002 were grown in a paddy field, and fully expanded flag leaves at the heading stage were used for DIS. (a) Detached leaves were incubated in 3 mM MES buffer (pH 5.8) at 27 °C under dark conditions. (b,c) Total chlorophyll contents (n= 6) and Fv/Fm ratio (n= 5) were compared between Hwaseong and CR2002. Error bars indicate a standard error, and more than three samples were used for each experiment. * and *** indicate significant difference at p < 0.05 and p < 0.001 based on Student’s t-test, respectively.
Ijms 21 01704 g002
Ijms 21 01704 g003 550
Figure 3. Expression of senescence-associated genes and chlorophyll degradation genes in Hwaseong and CR2002. qRT-PCR was conducted to determine the transcript level of genes. OsUBQ5 was used for normalization. Error bars indicate the standard deviation of three replications. *, **, and *** indicate significant difference at p < 0.05, p < 0.01, and p < 0.001 based on Student’s t-test, respectively. DAI: days after incubation.
Figure 3. Expression of senescence-associated genes and chlorophyll degradation genes in Hwaseong and CR2002. qRT-PCR was conducted to determine the transcript level of genes. OsUBQ5 was used for normalization. Error bars indicate the standard deviation of three replications. *, **, and *** indicate significant difference at p < 0.05, p < 0.01, and p < 0.001 based on Student’s t-test, respectively. DAI: days after incubation.
Ijms 21 01704 g003
Ijms 21 01704 g004 550
Figure 4. OsGW2 T-DNA insertion knock-out mutant (gw2-ko) in wild-type Dongjin genetic background. (a) Schematic diagram showing the location of T-DNA insertion of the PFG_1B-10017.R mutant. (b) T-DNA insertion was confirmed using two sets of primers indicated in (a), and lines homozygous for the T-DNA insertion were used for this experiment. (c) Comparison of grain shape between wild type and gw2-ko. (d) Confirmation of OsGW2 knock-out by semi-quantitative RT-PCR. OsUBQ5 was used as a loading control. (e) Gene expression pattern of OsGW2 during the dark-induced senescence. The transcript level of genes was determined by qRT-PCR. OsUBQ5 was used for normalization.
Figure 4. OsGW2 T-DNA insertion knock-out mutant (gw2-ko) in wild-type Dongjin genetic background. (a) Schematic diagram showing the location of T-DNA insertion of the PFG_1B-10017.R mutant. (b) T-DNA insertion was confirmed using two sets of primers indicated in (a), and lines homozygous for the T-DNA insertion were used for this experiment. (c) Comparison of grain shape between wild type and gw2-ko. (d) Confirmation of OsGW2 knock-out by semi-quantitative RT-PCR. OsUBQ5 was used as a loading control. (e) Gene expression pattern of OsGW2 during the dark-induced senescence. The transcript level of genes was determined by qRT-PCR. OsUBQ5 was used for normalization.
Ijms 21 01704 g004
Ijms 21 01704 g005 550
Figure 5. gw2-ko mutant showed delayed senescence. Comparison of (a) plant morphology (at 138 days after transplanting) and (b) SPAD value (n = 30) of wild-type Dongjin and gw2-ko. (c) Wild type and gw2-ko were grown in a paddy field, and fully expanded flag leaves at the heading stage were used for DIS. (d,e) Total chlorophyll contents (n = 6) and Fv/Fm ratio (n = 10) were compared between wild type and gw2-ko. Error bars indicate a standard error, and more than six samples were used for each experiment. * and *** indicate significant difference at p < 0.05 and p < 0.001 based on Student’s t-test, respectively.
Figure 5. gw2-ko mutant showed delayed senescence. Comparison of (a) plant morphology (at 138 days after transplanting) and (b) SPAD value (n = 30) of wild-type Dongjin and gw2-ko. (c) Wild type and gw2-ko were grown in a paddy field, and fully expanded flag leaves at the heading stage were used for DIS. (d,e) Total chlorophyll contents (n = 6) and Fv/Fm ratio (n = 10) were compared between wild type and gw2-ko. Error bars indicate a standard error, and more than six samples were used for each experiment. * and *** indicate significant difference at p < 0.05 and p < 0.001 based on Student’s t-test, respectively.
Ijms 21 01704 g005
Ijms 21 01704 g006 550
Figure 6. Boxplot of chlorophyll contents at (a) 0 day after incubation and (b) 5 days after incubation under the dark-induced senescence conditions in the F2 population derived from a cross between Dongjin and gw2-ko. Genotype was determined with two sets of primers, as indicated in Figure 4a—LP+RP and RB+RP. One-way ANOVA was conducted to determine the significant difference between genotypes. *DJ/DJ: homozygous for Dongjin, DJ/KO: heterozygous, KO/KO: homozygous for gw2-ko.
Figure 6. Boxplot of chlorophyll contents at (a) 0 day after incubation and (b) 5 days after incubation under the dark-induced senescence conditions in the F2 population derived from a cross between Dongjin and gw2-ko. Genotype was determined with two sets of primers, as indicated in Figure 4a—LP+RP and RB+RP. One-way ANOVA was conducted to determine the significant difference between genotypes. *DJ/DJ: homozygous for Dongjin, DJ/KO: heterozygous, KO/KO: homozygous for gw2-ko.
Ijms 21 01704 g006
Ijms 21 01704 g007 550
Figure 7. Expression of senescence-associated genes and chlorophyll degradation genes compared in Dongjin and gw2-ko. qRT-PCR was conducted to determine the transcript level of genes. OsUBQ5 was used for normalization. Error bars indicate the standard deviation of three replications. *, **, and *** indicate significant difference at p < 0.05, p < 0.01, and p < 0.001 based on Student’s t-test, respectively. DAI: days after incubation.
Figure 7. Expression of senescence-associated genes and chlorophyll degradation genes compared in Dongjin and gw2-ko. qRT-PCR was conducted to determine the transcript level of genes. OsUBQ5 was used for normalization. Error bars indicate the standard deviation of three replications. *, **, and *** indicate significant difference at p < 0.05, p < 0.01, and p < 0.001 based on Student’s t-test, respectively. DAI: days after incubation.
Ijms 21 01704 g007
Ijms 21 01704 g008 550
Figure 8. Summary of the numbers of differentially expressed genes (DEGs) upon incubation of leaves of two lines (Hwaseong and CR2002) in complete darkness. (a) A Venn diagram, showing the number of genes shared and distinct to each genotype. (b) Significantly enriched gene ontology (GO) terms for CR2002-specific DEGs. Blue and green bars indicate the input list of CR2002 and background/reference, respectively. H0, H3, G0, and G3 indicate Hwaseong day 0, Hwaseong day 3, CR2002 day 0, and CR2002 day 3 in dark condition, respectively.
Figure 8. Summary of the numbers of differentially expressed genes (DEGs) upon incubation of leaves of two lines (Hwaseong and CR2002) in complete darkness. (a) A Venn diagram, showing the number of genes shared and distinct to each genotype. (b) Significantly enriched gene ontology (GO) terms for CR2002-specific DEGs. Blue and green bars indicate the input list of CR2002 and background/reference, respectively. H0, H3, G0, and G3 indicate Hwaseong day 0, Hwaseong day 3, CR2002 day 0, and CR2002 day 3 in dark condition, respectively.
Ijms 21 01704 g008
Table
Table 1. Comparison of agronomic traits and chlorophyll content between Hwaseong and CR2002.
Table 1. Comparison of agronomic traits and chlorophyll content between Hwaseong and CR2002.
TraitHwaseongCR2002p-Value
Plant height (cm)97 ± 4.00 *101 ± 3.840.000
Stem diameter (mm)4.47 ± 0.245.26 ± 0.480.000
Panicle length (cm)21 ± 1.2819 ± 1.040.000
First internode diameter (mm)1.72 ± 0.121.99 ± 0.140.000
Grain length (mm)7.14 ± 0.457.46 ± 0.170.002
Grain width (mm)3.29 ± 0.133.83 ± 0.100.000
Grain thickness (mm)2.20 ± 0.062.56 ± 0.100.000
1000-grain weight (g)25.8 ± 0.6832.9 ± 1.040.000
Chlorophyll content I (mg/m2)440 ± 16.38457 ± 37.570.048
Chlorophyll content II (mg/m2)181 ± 50.57216 ± 50.860.038
* Data are presented as mean ± standard deviation. The student’s t-test was conducted for p-value. Chlorophyll content I and II were measured at the heading stage and 1 month after heading, respectively.
Table
Table 2. QTLs (quantitative trait loci) for chlorophyll content based on one-way ANOVA in the F3 and F4.
Table 2. QTLs (quantitative trait loci) for chlorophyll content based on one-way ANOVA in the F3 and F4.
Trait 1Gen.QTLChr.Markerp-ValueR2 (%)H/H 2G/G
Chlorophyll
content I		F3qCC11RM11302-RM113150.0023.48442452
F4qCC22RM7288-RM129830.02910.76546570
Chlorophyll
content II		F3qCC22RM128130.01810.68173189
F4qCC22RM12813-RM7288 0.00124.63497545
F3qCC66RM5840.01810.62454434
1 Chlorophyll content I and II were measured at heading stage and heading after 1 month, respectively. 2 H/H and G/G indicate the mean chlorophyll contents of homozygous for Hwaseong and O. grandiglumis genotypes, respectively.
Table
Table 3. List of differentially expressed genes (DEGs) associated with senescence, chlorophyll degradation, and phytohormone specifically detected in CR2002 and Hwaseong.
Table 3. List of differentially expressed genes (DEGs) associated with senescence, chlorophyll degradation, and phytohormone specifically detected in CR2002 and Hwaseong.
GeneCR2002SymbolDescription
log2FCAdj. p-Value
Os07g0410300−3.6760.037OsERF5Conserved hypothetical protein.
Os09g0532400−2.6530.031OsABI1Signal transduction response regulator, receiver region domain containing protein.
Os10g0392400−1.5770.044OsJAZ1Tify domain containing protein.
Os07g0658400−1.0130.012OsGLUFerredoxin-dependent glutamate synthase, Leaf senescence, and nitrogen remobilization.
Os02g06352001.0810.000OsNIT1Similar to Nitrilase 2.
Os10g03974001.3600.003OsDWF1Dim/dwf1 protein, Cell elongation protein DIMINUTO/Dwarf1, Brassinosteroid (BR) biosynthesis.
Os02g08302002.5920.003OsRR3A-type response regulator, Cytokinin signaling.
GeneHwaseongSymbolDescription
log2FCAdj. p-value
Os11g0143300−2.9660.011OsRR9A-type response regulator, Cytokinin signaling.
Os12g0139400−2.6710.006OsRR10A-type response regulator, Cytokinin signaling.
Os11g0143200−2.4230.047OsCPD1Similar to Cytochrome P450 90A1.
Os03g0856700−2.1260.027OsGA20ox1Gibberellin 20 oxidase 1.
Os04g0673300−1.9810.002OsRR6A-type response regulator, Cytokinin signaling.
Os02g0164900−1.5240.002OsARF6Similar to Auxin response factor 3.
Os01g0208600−1.4620.021OsSCAR1SCAR-like protein 2, Component of the suppressor of cAMP receptor/Wiskott-Aldrich syndrome protein family verprolin-homologous (SCAR/WAVE) complex, Actin organization, Panicle development, Regulation of water loss.
Os01g0718300−1.3280.005OsBRI1Brassinosteroid LRR receptor kinase, Similar to Brassinosteroid-insensitive 1.
Os01g0723100−1.3230.001 Senescence-associated family protein.
Os03g0265100−1.2290PLS2Glycosyl transferase, group 1 domain containing protein.
Os07g02090001.0960OsDGLDolichyl-diphosphooligosaccharide-protein glycosyltransferase 48 kDa subunit precursor, N-glycosylation.
Os02g03247001.1840.001 Similar to senescence-associated protein.
Os01g09276001.1920OsARF2Similar to Auxin response factor 2 (ARF1-binding protein).
Os01g07525001.3210.006OsERF2APETELA2/ethylene response factor (AP2/ERF) type transcription factor, Negative regulation of disease resistance, Negative regulation of salt tolerance.
Os03g03278001.8720.013OsNAPNAC Family transcriptional activator, Abiotic stress response, Positive regulator of leaf senescence.
Os10g04776001.9420.024ONAC120Similar to NAM / CUC2-like protein.
Os11g01269004.0260.048ONAC122NAC-domain protein, Drought tolerance.
Os04g05780005.5200.028OsACS2ACC synthase, Ethylene biosynthesis.

Share and Cite

MDPI and ACS Style

Shim, K.-C.; Kim, S.H.; Jeon, Y.-A.; Lee, H.-S.; Adeva, C.; Kang, J.-W.; Kim, H.-J.; Tai, T.H.; Ahn, S.-N. A RING-Type E3 Ubiquitin Ligase, OsGW2, Controls Chlorophyll Content and Dark-Induced Senescence in Rice. Int. J. Mol. Sci. 2020, 21, 1704. https://doi.org/10.3390/ijms21051704

AMA Style

Shim K-C, Kim SH, Jeon Y-A, Lee H-S, Adeva C, Kang J-W, Kim H-J, Tai TH, Ahn S-N. A RING-Type E3 Ubiquitin Ligase, OsGW2, Controls Chlorophyll Content and Dark-Induced Senescence in Rice. International Journal of Molecular Sciences. 2020; 21(5):1704. https://doi.org/10.3390/ijms21051704

Chicago/Turabian Style

Shim, Kyu-Chan, Sun Ha Kim, Yun-A Jeon, Hyun-Sook Lee, Cheryl Adeva, Ju-Won Kang, Hyun-Jung Kim, Thomas H Tai, and Sang-Nag Ahn. 2020. "A RING-Type E3 Ubiquitin Ligase, OsGW2, Controls Chlorophyll Content and Dark-Induced Senescence in Rice" International Journal of Molecular Sciences 21, no. 5: 1704. https://doi.org/10.3390/ijms21051704

APA Style

Shim, K. -C., Kim, S. H., Jeon, Y. -A., Lee, H. -S., Adeva, C., Kang, J. -W., Kim, H. -J., Tai, T. H., & Ahn, S. -N. (2020). A RING-Type E3 Ubiquitin Ligase, OsGW2, Controls Chlorophyll Content and Dark-Induced Senescence in Rice. International Journal of Molecular Sciences, 21(5), 1704. https://doi.org/10.3390/ijms21051704

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop