
==== Front
eLife
Elife
eLife
eLife
2050-084X
eLife Sciences Publications, Ltd

96553
10.7554/eLife.96553
version of record
Research Article
Plant Biology
OsNF-YB7 inactivates OsGLK1 to inhibit chlorophyll biosynthesis in rice embryo
Yang Zongju https://orcid.org/0000-0003-4700-049X
12†
Bai Tianqi 1†
E Zhiguo 3†
Niu Baixiao 12
Chen Chen https://orcid.org/0000-0002-9748-3111
chenchen@yzu.edu.cn
12
1 https://ror.org/03tqb8s11 Jiangsu Key Laboratory of Crop Genomics and Molecular Breeding/ Zhongshan Biological Breeding Laboratory, Agricultural College of Yangzhou University Yangzhou China
2 https://ror.org/03tqb8s11 Jiangsu Co-Innovation Center for Modern Production Technology of Grain Crops/ Key Laboratory of Plant Functional Genomics of the Ministry of Education, Agricultural College of Yangzhou University Yangzhou China
3 https://ror.org/05szcn205 State Key Laboratory of Rice Biology and Breeding, China National Rice Research Institute Hangzhou China
Casati Paula Reviewing Editor https://ror.org/040me2j46 Center of Photosynthetic And Biochemical Studies (CEFOBI) Argentina

Kleine-Vehn Jürgen Senior Editor https://ror.org/0245cg223 University of Freiburg Germany

† These authors contributed equally to this work.

17 9 2024
2024
13 RP9655305 2 2024
This manuscript was published as a preprint.05 2 2024

This manuscript was published as a reviewed preprint.22 3 2024

The reviewed preprint was revised.29 8 2024

© 2024, Yang, Bai, E et al
2024
Yang, Bai, E et al
https://creativecommons.org/licenses/by/4.0/ This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.

As a master regulator of seed development, Leafy Cotyledon 1 (LEC1) promotes chlorophyll (Chl) biosynthesis in Arabidopsis, but the mechanism underlying this remains poorly understood. Here, we found that loss of function of OsNF-YB7, a LEC1 homolog of rice, leads to chlorophyllous embryo, indicating that OsNF-YB7 plays an opposite role in Chl biosynthesis in rice compared with that in Arabidopsis. OsNF-YB7 regulates the expression of a group of genes responsible for Chl biosynthesis and photosynthesis by directly binding to their promoters. In addition, OsNF-YB7 interacts with Golden 2-Like 1 (OsGLK1) to inhibit the transactivation activity of OsGLK1, a key regulator of Chl biosynthesis. Moreover, OsNF-YB7 can directly repress OsGLK1 expression by recognizing its promoter in vivo, indicating the involvement of OsNF-YB7 in multiple regulatory layers of Chl biosynthesis in rice embryo. We propose that OsNF-YB7 functions as a transcriptional repressor to regulate Chl biosynthesis in rice embryo.

rice
chlorophyll biosynthesis
embryo
LEC1
Research organism

Other
http://dx.doi.org/10.13039/501100001809 National Natural Science Foundation of China 32170344 Chen Chen http://dx.doi.org/10.13039/501100001809 National Natural Science Foundation of China 32300689 Yang Zongju Project of Zhongshan Biological Breeding Laboratory BM2022008-02 Chen Chen Natural Science Foundation of the Higher Education Institutions of Jiangsu Province 21KJB21003 Yang Zongju http://dx.doi.org/10.13039/501100002858 China Postdoctoral Science Foundation 2022M712702 Yang Zongju http://dx.doi.org/10.13039/501100002949 Jiangsu Province Government JBGS001 Chen Chen http://dx.doi.org/10.13039/501100012246 Priority Academic Program Development of Jiangsu Higher Education Institutions PAPD Chen Chen Independent Scientific Research Project Funds of the Jiangsu Key Laboratory of Crop Genomics and Molecular Breeding PLR202101 Chen Chen The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.Author impact statementOsNF-YB7, an ortholog of Arabidopsis LEC1, is a negative regulator of Chl biosynthesis in rice embryos by interacting with OsGLK1 to inhibit the transactivity of OsGLK1.
publishing-routeprc
==== Body
pmcIntroduction

Angiosperms can be divided into chloroembryophytes and leucoembryophytes, depending on the presence or absence of chlorophyll (Chl) in the embryo, respectively (Puthur et al., 2013; Smolikova and Medvedev, 2016). Some plant species such as Arabidopsis (A. thaliana) produce chloroembryos, which have photochemically active chloroplasts capable of producing assimilates that are further converted into reserve biopolymers (Simkin et al., 2020). Grass species, such as rice (Oryza sativa), are unable to produce Chl during embryo development. To the best of our knowledge, it remains completely unclear what determines the ability to biosynthesize Chl in plant embryo. However, studies have identified several genes that contribute to Chl degradation for chloroembryos, the mutation of which can lead to a stay-green phenotype in mature seeds such as Mendel’s green pea (Sato et al., 2007; Delmas et al., 2013; Wang et al., 2018; Li et al., 2017; Thomas and Ougham, 2014).

Golden 2-Like (GLK) genes encode GARP-type transcription factors (TFs), which are key components regulating chloroplast development and Chl biosynthesis in plants (Rossini et al., 2001; Fitter et al., 2002; Nakamura et al., 2009). GLK can recognize the CCAATC cis-element of the Chl biosynthesis- and photosynthesis-associated nuclear genes to trigger their expression (Waters et al., 2009). The G-box (CACGTG) is also enriched in the GLK-targeted genes in Arabidopsis, possibly due to GLK being able to interact with certain G-box binding factors (Tamai et al., 2002). A genetic study showed that GLK activates Chl biosynthesis in roots, in a manner dependent on the G-box binding TF Elongated Hypocotyl 5 (HY5) (Kobayashi et al., 2012).

There are two copies of GLK in the rice genome, designated as OsGLK1 and OsGLK2 (Fitter et al., 2002). They redundantly regulate a set of genes, such as rice Chlorophyllide A Oxygenase (OsCAO) and Protochlorophyllide Oxidoreductase A (OsPORA) responsible for Chl biosynthesis, and Light Harvesting Complex B1 (OsLHCB1) and OsLHCB4 responsible for photosynthesis (Nakamura et al., 2009; Wang et al., 2013; Sakuraba et al., 2017). OsGLK1 overexpression in rice leads to green calli and chloroplast development in the vascular bundles (Nakamura et al., 2009; Wang et al., 2017). Organelle development in rice vascular sheath cells is induced by ectopically expressed maize (Zea mays) GLK genes, mimicking a key step in the evolutionary transition from C3 to C4 plants (Wang et al., 2017). In accordance with this, the ectopic expression of maize GLKs in rice can boost biomass and grain yield by facilitating chloroplast development and photosynthesis (Li et al., 2020; Yeh et al., 2022). Moreover, the overexpression of maize GLKs in calli was shown to improve the ability of rice and maize to regenerate (Luo et al., 2023). A recent study found that OsGLK1 also participates in tapetum plastid development and programmed cell death, consequently affecting pollen fertility in rice (Zheng et al., 2022). The findings indicate that GLK plays multiple roles in relation of plant development.

Leafy Cotyledon 1 (LEC1), a member of the nuclear factor Y (NF-Y) TF family, is a central regulator controlling many aspects of seed development, such as Chl accumulation in the embryo (Meinke et al., 1994; Pelletier et al., 2017). LEC1 can also act as a pioneer TF to regulate flowering by reprogramming the embryonic chromatin state (Tao et al., 2017). Previous studies have reported that lec1 mutants have paler green embryos than wild-type (WT) at maturation in Arabidopsis (Meinke et al., 1994; West et al., 1994). LEC1 transcriptionally regulates the expression of genes that encode light-reaction components of photosystems I and II, as well as the expression of genes involved in chloroplast biogenesis in Arabidopsis and soybean embryos (Pelletier et al., 2017; Jo et al., 2020). These findings suggest that LEC1 is important for photosynthesis and chloroplast development during seed development. However, the underlying molecular mechanisms remain largely unclear.

It is still a mystery why plant species such as rice cannot synthesize chlorophyll in embryos, while species such as Arabidopsis can. There are two LEC1 homologs, OsNF-YB7 and OsNF-YB9, encoded by the rice genome (Zhiguo et al., 2018). OsNF-YB7 is restricted to the embryo, and defective OsNF-YB7 may result in seed lethality (Niu et al., 2021b; Zhang and Xue, 2013). Here, we found that OsNF-YB7 acts as a key inhibitor of Chl biosynthesis in rice embryo. It inactivates the transactivation activity of OsGLK1, at multiple regulatory layers, to inhibit Chl accumulation in the embryo of rice, explaining the achlorophyllous embryo produced in rice.

Results

Loss of function of OsNF-YB7 leads to chloroembryos

By observing seeds produced by the loss-of-function mutant of OsNF-YB7, we surprisingly found that the osnf-yb7 embryo was greenish (Figure 1A and B), suggesting that OsNF-YB7 plays a negative role in chloroplast biogenesis or Chl biosynthesis, or both, during embryogenesis. We therefore examined the phenotype of green embryos at various seed development stages. The results showed that WT had an achlorophyllous embryo throughout embryonic development, whereas osnf-yb7 embryo turned green at 5 days after fertilization (DAF) and the chloroembryo remained green until maturity (Figure 1—figure supplement 1A–C). In support of this, transmission electron microscopy (TEM) showed well-developed chloroplast in the scutellum tip of the mutant, but this was not seen in the WT sections (Figure 1C). By measuring Chl contents in the developing and mature embryos, we found that total Chl content in osnf-yb7 was consistently higher than that in the WT, both at 10 DAF and at maturation (30 DAF) (Figure 1D). As revealed by confocal laser scanning microscopy (CLSM), Chl autofluorescence was detectable as early as 5 DAF in osnf-yb7, indicating the initiation of Chl accumulation at this stage (Figure 1E). This is consistent with our previous finding that OsNF-YB7 was highly activated at 5 DAF in rice embryo (Niu et al., 2021b).

Figure 1. OsNF-YB7 negatively regulates Chlorophyll (Chl) biosynthesis in embryo.

(A, B) Morphologies of wild-type (WT) and osnf-yb7 detached embryos at 10 days after fertilization (DAF) (A) and longitudinally dissected embryos at maturation (B). Scale bars = 200 μm. (C) Transmission electron microscopy images of embryos from WT and osnf-yb7 at 10 DAF. Scale bars = 5 μm. (D) Chl levels in WT and osnf-yb7 embryos at 10 DAF and maturation. Data are means ± SD (n=3). **, p<0.01; Student’s t-test was used for statistical analysis. (E) Chl autofluorescence of WT and mutant embryos at 5 DAF. BF, bright field; Chl, Chl autofluorescence. Scale bars = 100 μm. (F) Photosynthesis-related pathways enriched among the differentially expressed genes (DEGs) identified from the 5- and 10-DAF-old osnf-yb7 embryos compared to WT, respectively. (G) Most of the photosynthesis-related genes were upregulated in the mutant embryos. (H) A heat map shows the expression of the Chl biosynthesis-related genes in the WT and osnf-yb7 embryos at 5 and 10 DAF. Reads per kilobase per million mapped reads (RPKM) is used to indicate the expression level.

Figure 1—figure supplement 1. Mutation of OsNF-YB7 leads to chloroembryo.

Figure 1—figure supplement 2. Light is required but not sufficient for Chl biosynthesis in the chloroembryo of osnf-yb7.

Figure 1—figure supplement 3. Mutation of OsNF-YB7 actives the expression of photosynthesis-related genes.

Figure 1—figure supplement 4. OsNF-YB7 negatively regulates the expression of light harvest and Chl biosynthesis-associated genes.

Light is a critical signal triggering Chl biosynthesis (Wietrzynski and Engel, 2021). To determine whether the Chl biosynthesis in osnf-yb7 is induced by light, we investigated the embryo phenotype of WT and osnf-yb7 in the dark, using aluminum foil to cover rice panicles prior to flowering. The result showed that, similar to WT, osnf-yb7 embryo was achlorophyllous in the dark, although the embryogenesis defects, such as degenerated epiblast and coleorhiza, and maldeveloped coleoptile, were still observable (Figure 1—figure supplement 2A–H). This indicated that Chl biosynthesis in the mutant is light dependent. Only a small amount of light can be perceived by rice embryos because the external hulls block light penetration (Simkin et al., 2020). We removed the hulls to directly expose the embryo to light; however, Chl accumulation still failed to occur in the WT embryo (Figure 1—figure supplement 2I–L), suggesting that there are internal signals that repress Chl biosynthesis in rice embryo. We suspected that the OsNF-YB7 mutation may attenuate the activity of such inhibitors; alternatively, OsNF-YB7 itself could be an inhibitor.

OsNF-YB7 negatively regulates photosynthesis- and Chl biosynthesis-related genes

Using 5- and 10-DAF-old WT and osnf-yb7 embryos, we previously performed deep sequencing of the transcriptome (RNA-seq) to identify possible downstream genes of OsNF-YB7 (Niu et al., 2021a). As revealed by Mapman analysis, photosynthesis-related pathways, such as photosystem I, photosystem II, and the light reaction, were significantly enriched for the differentially expressed genes (DEGs) in the mutant embryos (Figure 1F, Figure 1—figure supplement 3A and B). Moreover, 96.9% (95/98) and 96.4% (106/110) of the photosynthesis-related DEGs that identified from 5- and 10–DAF-old embryos, respectively, were upregulated in the mutant (Figure 1G). To confirm this finding, we next examined the expression of OsLHCAs and OsLHCBs, which are primarily associated with photosystems I and II, respectively, via quantitative real-time PCR (RT-qPCR). It showed higher expression of all of the studied OsLHCAs and OsLHCBs in the osnf-yb7 embryo at 10 DAF (Figure 1—figure supplement 4A). Likewise, many of the genes participating in Chl biosynthesis, including rice Genomes Uncoupled 4 (OsGUN4), Mg-Chelatase H Subunit (OsCHLH), OsCHLI, OsCHLD, Copper Response Defect 1 (OsCRD1), OsPORA, OsPORB, and Divinyl Reductase (OsDVR), were significantly activated in the mutant (Figure 1H and Figure 1—figure supplement 4B). The findings indicated that OsNF-YB7 might act as a repressor of Chl biosynthesis and photosynthesis, which is opposite to the role of its homologue LEC1 in Arabidopsis.

Since OsNF-YB7 is a TF, we assumed that it may directly regulate the expression of genes related to Chl biosynthesis and photosynthesis, such as OsPORA, and OsLHCB4, which were significantly activated in the osnf-yb7 embryo (Figure 2A and B). We first generated transgenic lines that overexpressed OsNF-YB7, tagged with either green fluorescent protein (NF-YB7-GFP) or 3×Flag (NF-YB7-Flag) in the Zhonghua11 (ZH11, O. sativa ssp. geng/japonica) or Kitaake (O. sativa ssp. geng/japonica) background, respectively. As expected, OsPORA and OsLHCB4 were significantly downregulated in leaves of the OsNF-YB7 overexpressors (Figure 2C and D, and Figure 1—figure supplement 4C). Similar to previously reported (Zhang and Xue, 2013; Ito et al., 2011), the OsNF-YB7 overexpression lines displayed severe reproductive development defects, which prevented us from obtaining sufficient embryo tissues for subsequent experiments. Instead, using the NF-YB7-Flag seedling, we conducted a chromatin immunoprecipitation assay coupled with quantitative PCR (ChIP-qPCR). The results showed that OsNF-YB7 was highly enriched in the promoter regions of OsPORA and OsLHCB4 harboring the G-box motif (Figure 2E and F), a putative binding site of OsNF-YB7 (Guo et al., 2022). To confirm the ability of OsNF-YB7 to bind to the OsPORA and OsLHCB4 promoters, we next performed electrophoretic mobility-shift assays (EMSAs), using recombinant OsNF-YB7-His protein, and biotin-labeled subfragments of the OsPORA or OsLHCB4 promoters containing G-boxes. The results showed that OsNF-YB7-His was able to bind to the labeled probes, and the shifted band signals were substantially weakened upon application of the unlabeled cold probes or hot probes with a mutated G-box (Figure 2G and H).

Figure 2. OsNF-YB7 binds to the promoters of OsPORA and OsLHCB4 to regulate their expression.

(A, B) Quantitative real-time PCR (RT-qPCR) analysis of the transcription levels of OsPORA (A) and OsLHCB4 (B) in the embryos of WT and osnf-yb7 at 10 DAF. Data are means ± SD (n=3). **, p<0.01; Student’s t-test was used for statistical analysis. (C, D) Expression of OsPORA (C) and OsLHCB4 (D) in the leaves of WT and OsNF-YB7-overexpressing transgenic plants (NF-YB7-GFP). Data are means ± SD (n=3). **, p<0.01; Student’s t-test was used for statistical analysis. (E, F) Chromatin immunoprecipitation assay coupled with quantitative PCR (ChIP-qPCR) analyses showing the enrichment of OsNF-YB7 at the OsPORA (E) and OsLHCB4 (F) promoters in 14-day-old OsNF-YB7-Flag seedlings. Precipitated DNA was quantiﬁed by qPCR and DNA enrichment is displayed as a percentage of input DNA. Data are means ± SD (n=3). *, p<0.05; **, p<0.01; Student’s t-test was used for statistical analysis. ACTIN was used as a nonspeciﬁc target gene. Diagrams in the upper panel showing the promoter structures of OsPORA and OsLHCB4, and the PCR amplicons used for ChIP-qPCR. (G, H) Electrophoretic mobility-shift assays (EMSAs) showing that OsNF-YB7 directly binds to the promoters of OsPORA (G) and OsLHCB4 (H). Hot probes were biotin-labeled. The hot mProbes contain mutant nucleic acid from CACATG to AAAAAA. The arrow heads indicate the shift bands. (I) Schematic diagram displaying the constructs used in the dual luciferase reporter (DLR) assays of J. LUC, firefly luciferase; REN, Renilla luciferase. (J) DLR assays showing that OsNF-YB7 directly represses the promoter activities of OsPORA and OsLHCB4. Data are means ± SD (n=3). *, p<0.05; **, p<0.01; Student’s t-test was used for statistical analysis.

Figure 2—source data 1. Uncropped and labelled gels.

Figure 2—source data 2. Raw unedited gels.

To confirm that the binding of OsNF-YB7 represses target gene expression, dual luciferase (LUC) reporter (DLR) assays were performed. We first generated LUC reporters driven by the OsPORA or OsLHCB4 promoters (designated as proOsPORA:LUC and proOsLHCB4:LUC hereafter) (Figure 2I). When the reporters were coexpressed with OsNF-YB7, which acted as an effector, in rice protoplasts, we found that OsNF-YB7 significantly repressed the activity of both proOsPORA:LUC and proOsLHCB4:LUC (Figure 2J). Taken together, these results suggested that OsNF-YB7 directly binds to the promoters of photosynthesis- and Chl-biosynthesis-related genes and represses their transcription.

OsNF-YB7 represses OsGLK1 in the embryo

Several TFs that regulate Chl biosynthesis or chloroplast development have been identified in plants (Jarvis and López-Juez, 2013). Some of these, such as rice OsGLK1, HY5-like 1 (OsHY5L1), PIF-like 14 (OsPIL14), and GATA Nitrate-inducible Carbon-metabolism-involved (OsGNC), were found to be upregulated in the embryo of osnf-yb7 (Figure 3A).

Figure 3. OsNF-YB7 associates with the promoter of OsGLK1 and represses its expression.

(A) A heat map showing the expression of transcription factors associated with Chl biosynthesis and chloroplast development in the 5- and 10-DAF-old embryos of WT and osnf-yb7. The colored dots indicate log2(RPKM mean) of the genes in three biological replicates. (B) RT-qPCR analysis of OsGLK1, OsGNC, and OsPIL14 expression levels in the embryos from WT and osnf-yb7 at 5- and 10 DAF. Numbers represent fold changes of expression. Data are means ± SD (n=3). *, p<0.05; **, p<0.01; Student’s t-test was used for statistical analysis. (C) RT-qPCR analysis of OsGLK1 expression levels in leaves from WT and NF-YB7-GFP. Data are means ± SD (n=3). **, p<0.01; Student’s t-test was used for statistical analysis. (D) ChIP-qPCR analysis showing the enrichment of OsNF-YB7 at the OsGLK1 promoter in 14-day-old OsNF-YB7-Flag seedling. Chromatin of each sample was immunoprecipitated using anti-Flag or IgG antibodies. Precipitated DNA was quantiﬁed by qPCR and DNA enrichment is displayed as a percentage of input DNA. Data are means ± SD (n=3). *, p<0.05; **, p<0.01; Student’s t-test was used for statistical analysis. ACTIN was used as a nonspeciﬁc target gene. The experiment was performed three times with similar results. The diagram in the upper panel showing the promoter structure of OsGLK1 and PCR amplicons (P1, P2, P3, P4, and P5) used for ChIP-qPCR. (E) Schematic diagram displaying the constructs used in the DLR assays of F. LUC, firefly luciferase; REN, Renilla luciferase; UAS, upstream activating sequence. (F) DLR assays showing that OsNF-YB7 represses the promoter activity of OsGLK1. Data are means ± SD (n=3). *, p<0.05; Student’s t-test was used for statistical analysis.

Figure 3—figure supplement 1. OsNF-YB7 does not directly binds to the promoter of OsGLK1 in vitro.

Figure 3—figure supplement 1—source data 1. Uncropped and labelled gels.

Figure 3—figure supplement 1—source data 2. Raw unedited gels.

OsGLK1 was the most strikingly activated TF (Figure 3A and B), while it was significantly repressed in NF-YB7-GFP (Figure 3C). To test the idea that OsGLK1 is a direct downstream target of OsNF-YB7, a ChIP-qPCR assay was first carried out using NF-YB7-Flag transgenic seedlings. The results showed that the promoter segments P3 and P4 were significantly enriched in the immunoprecipitated chromatin (Figure 3D), suggesting that OsNF-YB7 was able to bind to the OsGLK1 promoter in vivo. However, we failed to validate the binding ability in vitro using the EMSA assay (Figure 3—figure supplement 1), suggesting that OsNF-YB7 requires other NF-Ys or TFs to form a TF complex to recognize the promoter, as many NF-Y members do for function (Laloum et al., 2013). A DLR assay was next performed to investigate the negative regulation of OsNF-YB7 on OsGLK1 transcription using rice protoplasts. The OsGLK1 promoter was inserted upstream of the 5×upstream activating sequence (UAS) as a reporter; OsNF-YB7 was fused C-terminally to the DNA binding domain (BD) of the yeast GAL4 and the herpes virus VP16 transactivation domain (VP16:OsNF-YB7), as an effector (Figure 3E). Coexpression of VP16:OsNF-YB7 with the reporter in rice protoplasts significantly decreased the transcriptional activity of VP16 (Figure 3F), indicating that OsNF-YB7 represses OsGLK1 promoter activity.

OsGLK1 is involved in OsNF-YB7-regulated Chl biosynthesis in embryo

To confirm the contribution of OsGLK1 for the production of chloroembryo in osnf-yb7, we first generated OsGLK1-overexpressing lines (OsGLK1-OX) driven by the rice ubiquitin promoter. In association with over-accumulated Chl in the glume and seed coat of the transformant, RT-qPCR and Western blot assays confirmed OsGLK1 activation in OsGLK1-OX (Figure 4—figure supplement 1A-C). As observed in osnf-yb7, green embryos were produced (Figure 4A). The Chl content in OsGLK-OX were higher than that in WT (Figure 4B and C), suggesting that OsGLK1 overexpression in rice embryo induces Chl biosynthesis.

Figure 4. Chl biosynthesis in osnf-yb7 embryo requires active OsGLKs.

(A) Embryo morphologies of WT and OsGLK1-OX detached embryos at 10 DAF. Scale bars = 2 mm. (B) Chl autofluorescence of the WT and GLK-OX embryos at 10 DAF. BF, bright field; Chl, Chl autofluorescence. Scale bars = 100 μm. (C) Chl levels in WT and OsGLK1-OX embryos at maturation. Data are means ± SD (n=3). **, p<0.01; Student’s t-test was used for statistical analysis. (D–G) Morphologies of the embryos produced by WT (D), osnf-yb7 (E), osnf-yb7;osglk2 double mutant (F) and osnf-yb7;osglk1;osglk2 triple mutant (G). Scale bars = 1 mm. (H) A Venn diagram showing overlaps of the DEGs identified from the embryos of osnf-yb7 and OsGLK1-OX at 10 DAF. (I) A pie chart showing similar transcriptional changes of the common DEGs identified from osnf-yb7 and OsGLK1-OX. (J) Gene Ontology (GO) analysis of the common DEGs identified from osnf-yb7 and OsGLK1-OX.

Figure 4—source data 1. Differentially expressed genes (DEGs) in the 10-DAF-old embryos of OsGLK1-OX.

Figure 4—source data 2. Common DEGs identified from the OsGLK1-OX and osnf-yb7 embryos at 10 DAF.

Figure 4—figure supplement 1. Overexpression of OsGLK1 induces chloroembryo in rice.

Figure 4—figure supplement 1—source data 1. Uncropped and labelled gels.

Figure 4—figure supplement 1—source data 2. Raw unedited gels.

Figure 4—figure supplement 2. Generation of osnf-yb7, osglk1, and osglk2 high-order mutants.

Figure 4—figure supplement 3. Chl levels in WT, osnf-yb7 and osnf-yb7;osglk1;osglk2 embryos at 10 DAF.

Figure 4—figure supplement 4. Overexpression of OsGLK1 activates the expression of Chl biosynthesis- and photosynthesis-related genes.

Using a clustered regularly interspaced short palindromic repeats (CRISPR)/CRISPR-associated 9 (Cas9) cassette including three tandemly arrayed guide RNAs targeting OsNF-YB7, OsGLK1, and OsGLK2, respectively, we successfully obtained the osnf-yb7;osglk2 double mutant and the osnf-yb7;osglk1;osglk2 triple mutant (Figure 4—figure supplement 2A and B). By analyzing the embryos that the mutants produced, we found that significantly less Chl accumulated in osnf-yb7;osglk1;osglk2 than in osnf-yb7 and osnf-yb7;osglk2 (Figure 4D–G and Figure 4—figure supplement 3). In comparison to the achlorophyllous embryo of WT, the osnf-yb7;osglk1;osglk2 triple mutant still showed somewhat green coloration in the apical part of the embryos (Figure 4G). We hypothesized that this was at least partially due to that genes like OsLHCB4 and OsPORA can be induced by the mutation of OsNF-YB7, given the fact that OsNF-YB7 represses the genes’ expression independent of OsGLK1 (Figure 2J).

OsNF-YB7 and OsGLK1 regulate a common set of genes in the embryo

We next performed RNA-seq analysis to explore the transcriptomic changes in the chloroembryo of OsGLK1-OX at 10 DAF. More than 64.4% of the DEGs identified in OsGLK1-OX were overlapped with the ones identified in osnf-yb7, the vast majority of which were either activated (60.6%) or repressed (27.6%) in both OsGLK1-OX and osnf-yb7 (Figure 4H and I, and Figure 4—source data 1 and 2). As revealed by the Gene Ontology (GO) analysis, genes involved in Chl biosynthesis and photosynthesis, such as OsPORA and OsLHCB4, were markedly enriched among the common DEGs (Figure 4J and Figure 4—figure supplement 4), implying that OsNF-YB7 and OsGLK1 antagonistically regulate a common set of genes for Chl biosynthesis and photosynthesis.

The ChIP-qPCR results showed that OsGLK1 associated with the regions of the OsPORA and OsLHCB4 promoters to which OsNF-YB7 binds (Figure 5A and B). In agreement with this, the EMSA results suggested that OsGLK1 directly binds to the same DNA probes of OsPORA and OsLHCB4 in vitro (Figure 5C and D). In opposite to OsNF-YB7, when we coexpressed the OsGLK1 effector vector with the reporter vector proOsLHCB4:LUC or proOsPORA:LUC in rice protoplasts, OsGLK1 showed significant transactivation activity on OsPORA and OsLHCB4 (Figure 5—figure supplement 1A and B).

Figure 5. OsNF-YB7 and OsGLK1 regulate a common set of genes in the embryo.

(A, B) ChIP-qPCR analysis showing the enrichment of OsGLK1 on the OsPORA (A) and OsLHCB4 (B) promoters. PCR amplicons used for ChIP-qPCR are indicated in the schematic diagrams. OsGLK1-GFP was transiently expressed in protoplasts isolated from green tissues of the 14-day-old WT seedling. Chromatin of each sample was immunoprecipitated using anti-GFP or igG antibodies. Precipitated DNA was quantiﬁed by qPCR and DNA enrichment is displayed as a percentage of input DNA. Data are means ± SD (n=3). *, p<0.05; **, p<0.01; Student’s t-test was used for statistical analysis. ACTIN was used as a nonspeciﬁc target gene. The experiment was performed three times with similar results. (C, D) EMSA assays showing that OsGLK1 directly binds to the promoters of OsPORA (C) and OsLHCB4 (D). The arrow heads indicate the shift bands. (E) A Venn diagram showing overlaps of the target genes between OsGLK1 and OsNF-YB7. (F) Distribution of the distance between the OsNF-YB7 and OsGLK1 summits showing that OsNF-YB7 and OsGLK1 bind to proximal regions of their common targets.

Figure 5—source data 1. Common targets of OsGLK1 and OsNF-YB7.

Figure 5—source data 2. Uncropped and labelled gels.

Figure 5—source data 3. Raw unedited gels.

Figure 5—figure supplement 1. OsGLK1 significantly activates the promoter activities of OsPORA and OsLHCB4.

Figure 5—figure supplement 2. Many of the OsNF-YB7 and OsGLK1 common targets that activated in osnf-yb7 and OsGLK1-OX are involved in Chl biosynthesis and photosynthesis.

Figure 5—figure supplement 3. OsNF-YB7 and OsGLK1 bind to proximal regions of their common targets.

Recently, the putative binding sites of OsNF-YB7 and OsGLK1 were investigated at the whole genome scale (Guo et al., 2022; Tu et al., 2022). This allows us to test our hypothesis that OsNF-YB7 and OsGLK1 can target a common set of genes involved in Chl biosynthesis and photosynthesis. By reanalyzing the ChIP-seq data, we found that 91.4% (235/257) and 90.7% (167/184) of the OsGLK1-binding genes overlapped with the OsNF-YB7-binding genes in two replicates (Figure 5E and Figure 5—source data 1), although the peak number of OsNF-YB7 was much greater than that of OsGLK1 (Guo et al., 2022; Tu et al., 2022). A large number of the common targets were activated in the embryos of osnf-yb7 and OsGLK1-OX at 10 DAF (Figure 5—figure supplement 2A). GO analysis suggested that most of the common targets were genes involved in Chl biosynthesis or photosynthesis (Figure 5—figure supplement 2B). By retrieving the sequences of the OsNF-YB7 and OsGLK1 binding peaks in the common targets for MEME analysis, we found that the TFs probably recognize similar DNA motifs. A short sequence containing the G-box motif was the most significantly enriched (Figure 5—figure supplement 3A and B). There was also enrichment of another sequence containing a CCAAT motif recognized by the NF-Y TF complexes and a CCAATC motif recognized by GLKs, as previously reported (Waters et al., 2009; Pelletier et al., 2017; Figure 5—figure supplement 3A and B). The results suggested that OsNF-YB7 and OsGLK1 can bind to the same region of their common target. In agreement with this, approximately 75.0% and 76.5% of the peak summits of OsGLK1 were located proximally to the summit of OsNF-YB7 in two replicates, within an adjacent region no more than 100 bp away (Figure 5F and Figure 5—source data 1). Consistent with our biochemical evidence, the ChIP-seq results showed that OsNF-YB7 does bind to the same regions of the OsPORA and OsLHCB4 promoters as OsGLK1 binds to (Figure 5—figure supplement 3C).

OsNF-YB7 physically interacts with OsGLK1 to inhibit its transcriptional activity

To resolve how OsNF-YB7 and OsGLK1 bind to the same regions to regulate their common targets, we speculated that the TFs probably form a dimer in rice in order to exert their functions. We therefore carried out a yeast two-hybrid (Y2H) assay by transforming yeast cells with a bait construct expressing OsNF-YB7 fused with the GAL4 DNA-binding domain (BD), together with a prey construct expressing OsGLK1 fused with the yeast GAL4 activation domain (AD). The results showed that OsNF-YB7 interacts with OsGLK1 in yeast (Figure 6A and B). Furthermore, to determine the functional domains required for the interaction, we generated two truncated versions of OsNF-YB7, which contained the N- or C-terminal, and three truncated versions of OsGLK1, which contained the N-terminal, DNA-BD, or GCT box domain, as previously reported (Fitter et al., 2002; Zhang et al., 2021a; Figure 6A). We found that the C-terminus of OsNF-YB7 was sufficient for the interaction (Figure 6B). In addition, the GCT box domain of OsGLK1 strongly interacted with the full length or C-terminus of OsNF-YB7, while the DNA-BD of OsGLK1 showed a weak interaction (Figure 6B).

Figure 6. OsNF-YB7 interacts with OsGLK1 to regulate the expression of OsPORA and OsLHCB4.

(A) Schematic diagrams showing the protein structures of OsGLK1 and OsNF-YB7. (B) Yeast-two-hybrid (Y2H) assays showing the interaction between OsNF-YB7 and OsGLK1. AD and BD indicate the activation domain and binding domain of GAL4, respectively. The full length or truncated OsNF-YB7 and OsGLK1 were fused with BD and AD, respectively. The indicated combinations of constructs were cotransformed into yeast cells and grown on the nonselective medium SD/-L-W and selective medium SD/-L-W-H. (C) A split complementary luciferase (LUC) assay confirmed the interaction between OsNF-YB7 and OsGLK1. Coexpression of the fusion of OsGLK1 and the N-terminal half of LUC (nLUC-OsGLK1) and the fusion of the C-terminal half of LUC and OsNF-YB7 (cLUC-OsNF-YB7) in the epidermal cells of N. benthamiana leaves induced LUC activities, whereas the epidermal cells coexpressing nLUC-OsGLK1 and cLUC, nLUC and cLUC-OsNF-YB7, or nLUC and cLUC did not show LUC activities. (D) Bimolecular fluorescence complementation (BiFC) assays showed interactions between OsNF-YB7 and OsGLK1 in the nucleus. OsNF-YB7 was fused with the N-terminal of yellow fluorescent protein (nYFP-OsNF-YB7); OsGLK1 was fused with the C-terminal of yellow fluorescent protein (cYFP-OsGLK1). The indicated combinations of constructs were transiently coexpressed in the leaf epidermal cells of N. benthamiana. Scale bar = 10 μm. (E) Co-immunoprecipitation (Co-IP) assays showing that OsNF-YB7 interacts with OsGLK1 in vivo. 35 S::OsNF-YB7:GFP (OsNF-YB7-GFP) and 35 S::OsGLK1:3×Flag (OsGLK1-Flag) were coexpressed in rice protoplasts and were immunoprecipitated with an anti-GFP antibody, and the immunoblots were probed with anti-GFP and anti-Flag antibodies. 35 S::GFP (GFP) was used as a negative control. (F) DLR assays showing that OsNF-YB7 represses the transactivation activity of OsGLK1 on OsPORA and OsLHCB4. Protoplasts isolated from etiolated seedlings were used for the analyses. EV, empty vector. Data are means ± SD (n=3). *, p<0.05; **, p<0.01; Student’s t-test was used for statistical analysis. (G) EMSA assay indicated that OsNF-YB7 inhibits the DNA binding of OsGLK1 to the promoter of OsPORA. The black and white arrow heads indicate the OsGLK1- and OsNF-YB7-bound probes, respectively; “+” and “++” indicate that 2- and 4 µM recombinant proteins were used for the reactions. The GST-His was used as a negative control.

Figure 6—source data 1. Uncropped and labelled gels.

Figure 6—source data 2. Raw unedited gels.

Figure 6—figure supplement 1. OsNF-YB7 interacts with OsGLK2.

Figure 6—figure supplement 2. A proposed model of the OsNF-YB7-mediated suppression of Chl biosynthesis in rice embryo.

A split complementary LUC assay further confirmed the interaction between OsNF-YB7 and OsGLK1 in the epidermal cells of Nicotiana benthamiana (Figure 6C). As suggested by the bimolecular fluorescence complementation (BiFC) analysis, the interaction occurred exclusively in the nuclei of N. benthamiana epidermal cells (Figure 6D). Moreover, we transiently coexpressed OsNF-YB7 tagged with GFP (OsNF-YB7-GFP) and OsGLK1 tagged with 3×Flag (OsGLK1-Flag) in rice protoplasts, and co-immunoprecipitation (Co-IP) analysis showed that OsGLK1-Flag could be immunoprecipitated by the anti-GFP antibody (Figure 6E), indicating that the interactions do occur in vivo. In addition, both Y2H and split complementary LUC assays showed that OsGLK2 could interact with OsNF-YB7 (Figure 6—figure supplement 1A and B). These findings indicated that OsNF-YB7 interacts with OsGLKs, explaining why OsNF-YB7 and OsGLKs share a common set of targets in rice.

To explore the biological meaning of the interaction, we next performed DLR assays using rice protoplasts. Transient expression of OsGLK1 in rice protoplasts substantially activated the reporters, driven by either the OsPORA or the OsLHCB4 promoter; however, when we coexpressed OsNF-YB7 with OsGLK1, the transactivation ability of OsGLK1 was significantly repressed in rice protoplasts (Figure 6F). These findings suggested that OsNF-YB7/OsGLK1 dimerization reduces the transactivation ability of OsGLK1 for fine-tuning the Chl biosynthetic and photosynthetic genes, such as OsPORA and OsLHCB4. The EMSA assay showed that OsGLK1-MBP recombinant proteins could bind to the promoter of OsPORA; however, when we incubated the probe with OsGLK1-MBP and OsNF-YB7-His together, the binding ability of OsGLK1 substantially decreased (Figure 6G). As a control, when the probe was incubated with OsGLK1-MBP and GST-His, the binding ability of OsGLK1 remained unchanged (Figure 6G). These findings implied that the reduced transactivity of OsGLK1 is likely due to the formation of OsNF-YB7/OsGLK1 heterodimers inhibiting the binding of OsGLK1 to its downstream genes.

Discussion

In crops such as soybean and canola, the presence of green embryos is considered as a valuable trait due to its association with increased photosynthetic capacity, which consequently promotes fatty acid biosynthesis (Ruuska et al., 2004). This highlights the potential value of the green embryo. Nevertheless, the degradation of chlorophyll in mature seeds is necessary to prevent adverse effects on seed viability and meal quality (Chung et al., 2006).

By surveying 1094 species from 666 genera and 182 families, Yakovlev and Zhukova, 1980 found that 428 angiosperms produce embryos with the presence of chlorophyll. Seeds with chlorophyllous embryo are scattered throughout the angiosperms and have evolved in many unrelated groups, but they are more prevalent in non-endospermous seeds (Dahlgren, 1980). Currently, the evolutionary force driving the divergence of chloroembryophytes and leucoembryophytes remains largely unknown. Embryo photosynthesis contributes a large amount of oxygen to fuel energy-generating pathways in seed (Simkin et al., 2020). The presence of chlorophyll in the embryo facilitates photosynthesis at early developmental stages, potentially leading to improved seedling growth and vigor (Smolikova and Medvedev, 2016). In many chloroembryophytes, such as Arabidopsis, the embryo occupies a large proportion of the seed. From an evolutionary perspective, the presence of chlorophyll in the embryo may promote adaptation in such chloroembryophytes, as more reserves can be accumulated in the seed through active photosynthesis, better supporting the embryo development and subsequent seedling growth (Sela et al., 2020). On the other hand, some leucoembryophytes, such as rice, have a persistent endosperm rich in storage reserves to nourish embryo development (Liu et al., 2022). The acquisition of the ability to accumulate chlorophyll in the embryo is not necessary in such species. Although chlorophyllous embryos are rare among the ‘primitive’ angiosperm superorders, they are observed in some (but not all) Nymphaeales species that are surrounded by some amount of endosperm (Yakovlev and Zhukova, 1980). As there is a paucity of knowledge regarding the evolution of chlorophyllous embryos, further comprehensive studies are necessary.

Some algae and gymnosperm species have evolved an ability to synthesize Chl in the dark (Myers, 1940; Ranade et al., 2019; Bogorad, 1950). However, it has remained largely unclear whether access to light accounts for the induction of chloroembryos (Dahlgren, 1980; Periasamy and Vivekanandan, 1981; Liu et al., 2017). Here, we found that removing the tissues covering an embryo failed to produce chlorophyllous embryos in the WT; however, light avoidance inhibited Chl accumulation in osnf-yb7 embryos (Figure 1—figure supplement 2A–H). These results suggest that light is necessary but insufficient to trigger chloroplast biogenesis in rice. Because light itself does not induce Chl accumulation in rice embryo (Figure 1—figure supplement 2I–L), we inferred that there are intrinsic cues to repress chloroplast development in rice embryo. However, to the best of our knowledge, the underlying mechanism that determines an embryo’s ability to synthesize Chl is completely unknown. Here, we showed that OsNF-YB7, a LEC1 homolog of rice, acts as an inhibitor to repress Chl accumulation in the embryo. In line with this, a recent independent study also showed that the OsNF-YB7 null mutant accumulates Chl in the embryo, although the cause of this remains unresolved (Guo et al., 2022).

Several TFs involved in chloroplast development and photomorphogenesis in rice have been identified (Nakamura et al., 2009; Li et al., 2019; Bai et al., 2019; Hudson et al., 2013). Interestingly, we found that many of them were upregulated in the osnf-yb7 embryos. For example, OsGLK1 was shown to be the most activated (Figure 3A and B). Our genetic and biochemical evidence suggests that OsGLK1 is involved in the OsNF-YB7-mediated repression of Chl biosynthesis in rice embryo. OsGLK1 overexpression mimicked the chloroembryo phenotype in WT, while knockout of OsGLK1 and OsGLK2, simultaneously, suppressed Chl accumulation in osnf-yb7 embryo (Figure 4A–G). The biochemical experiments suggested that OsNF-YB7 associates with the promoter of OsGLK1, in vivo, to transcriptionally inactivate OsGLK1 (Figure 3D–F), indicating that OsGLK1 is a downstream target of OsNF-YB7. However, OsNF-YB7 alone failed to bind to the promoter of OsGLK1 in vitro (Figure 3—figure supplement 1), presumably due to some as-yet-unidentified TFs being recruited by OsNF-YB7, assisting in recognizing OsGLK1 for transcriptional regulation. An NF-Y TF complex usually consists of three subunits (NF-YA/B/C) in order to exert its function (Laloum et al., 2013). NF-YA is responsible for DNA binding, while NF-YB and NF-YC are primarily responsible for transactivation (Chaves-Sanjuan et al., 2021; Gnesutta et al., 2017). OsNF-YB7 possibly forms a heterotrimer with NF-YA and NF-YC in the embryo to recognize OsGLK1. In addition, previous studies showed that Arabidopsis LEC1 can interact with different TFs for developmental regulation (Huang et al., 2015b; Yamamoto et al., 2009; Fatihi et al., 2016; Boulard et al., 2018; Huang et al., 2015a). Identifying such TFs may reveal how OsNF-YB7 recognizes OsGLK1 in order to exert its function.

OsNF-YB7 could recognize a common set of genes involved in Chl biosynthesis and photosynthesis that are recognized by OsGLK1 (Figure 2E–H and Figure 5—source data 1). The biochemical results suggested that OsNF-YB7 can directly repress OsPORA and OsLHCB4, which are activated by OsGLK1 (Figure 2J, and Figure 5—figure supplement 1). OsNF-YB7 probably hinders OsGLK1 to access the target genes by forming an OsGLK1/OsNF-YB7 heterodimer (Figure 6A–E), or by occupying the motif OsGLK1 recognized in the promoter, given the binding sites of OsGLK1 and OsNF-YB7 are likely overlapped (Figure 5F, Figure 5—figure supplement 3 and Figure 5—source data 1). The first of these hypotheses is more plausible because, upon co-incubation of OsGLK1 and OsNF-YB7 with the OsPORA promoter in vitro, the shifted band signal of OsGLK1 was substantially decreased (Figure 6G), indicating that the protein-protein interactions overwhelm the protein-DNA interactions. Thus, when OsNF-YB7 is expressed, it interacts with OsGLK1, then OsGLK1 is less available to activate downstream targets. The findings suggested that OsNF-YB7 plays a dual role in regulating Chl biosynthesis in rice embryo: first, it represses the downstream genes, achieved via its function as a transcriptional inactivator; second, OsNF-YB7 can interact with OsGLKs to disturb their abilities to transactivate genes related to Chl biosynthesis and photosynthesis (Figure 6—figure supplement 2). In addition to OsNF-YB7, a recent study showed that Deep Green Panicle 1, a plant-specific protein with a conserved TIGR01589 domain, can interact with OsGLKs to suppress OsGLK-mediated transcription (Zhang et al., 2021b). The findings suggested that, as the central regulator responsible for chloroplast development, GLK is tightly regulated at the post-translational level to fine-tune Chl biosynthesis in plants.

Previous studies suggested that Arabidopsis LEC1 is a positive regulator of Chl biosynthesis (Pelletier et al., 2017), given that the mature embryos of lec1 were paler than the WT (Meinke, 1992; West et al., 1994). However, the null mutation of OsNF-YB7 activated Chl biosynthesis, implying that the LEC1-type gene acts as a negative regulator in rice. By surveying the literature, we noted that Meinke reported that the cotyledons of the lec1 mutant remained green unusually late in development (Meinke, 1992). Although there was no significant difference in Chl content upon using whole seeds for quantification, Parcy et al. did observe that the tip of lec1 cotyledons accumulated more Chl (Parcy et al., 1997). Moreover, the lec1;abi3 double mutant embryos produced much more Chl than the abi3 single mutant (Parcy et al., 1997). Pelletier et al. recently reported that a cluster of LEC1 targets is enriched in genes related to photosynthesis and chloroplast development, and many of which are down-regulated in the embryo of lec1 (Pelletier et al., 2017); however, we noticed that in their dataset, the photosynthesis and Chl biosynthesis related genes were more strikingly enriched in the up-regulated genes of lec1, at either the mature green or postmature green stage. These observations challenge the concept that LEC1 positively regulates Chl biosynthesis and photosynthesis in Arabidopsis. Most studies in Arabidopsis have emphasized the importance of LEC1 in embryo development at the maturation stage. However, LEC1 is activated within 24 hr after fertilization (Lotan et al., 1998), but its role in the early embryo developmental stages for Chl biosynthesis is still unknown. Previously, we have expressed OsNF-YB7 in the lec1-1 background, driven by the native promoter of Arabidopsis LEC1 (Niu et al., 2021b). Given that OsNF-YB7 could rescue embryo morphogenesis defects in Arabidopsis (Niu et al., 2021b), we postulated that OsNF-YB7 in rice plays a similar role to that of LEC1 in Arabidopsis. In order to ascertain whether LEC1 can fully restore osnf-yb7, it is necessary to ectopically express LEC1 driven by the native OsNF-YB7 promoter in the osnf-yb7 background in the future, due to the possibility of functional divergence between the genes with regard to the regulation of chlorophyll biosynthesis in the embryo. It would also be worthwhile to carry out similar tests in other grass organisms to better understand the regulatory mechanism in rice.

In addition to producing chlorophyllous embryo, the osnf-yb7 mutants display an array of developmental defects, including abnormal embryogenesis, reduced dormancy, and desiccation intolerance, similar to those found in Arabidopsis lec1 mutants (Niu et al., 2021a). Mutation of osglk1 and osglk2 in osnf-yb7 could recover the chloroembyo phenotype, but did not alleviate the other embryo defects (Figure 4D–G). The findings suggest that OsGLKs specifically function in Chl biosynthesis, but that OsNF-YB7 is responsible for many aspects of embryo development. In agreement with this, the number of OsNF-YB7-targeted genes is far greater than that of OsGLK-targeted genes, and the DEGs in osnf-yb7 embryo is far greater than that in OsGLK1-OX embryo (Figure 4H and Figure 5E). The mechanisms underlying the OsNF-YB7-regulated multiple embryo developmental processes require further investigation.

Materials and methods

Plant materials and growth conditions

The osnf-yb7 mutant lines used in this study were previously generated in our laboratory (Niu et al., 2021b). Rice cultivars Zhonghua11 (ZH11) and Kitaake (Kit) were used for gene transformation. Rice plants were grown in a paddy field in Yangzhou, Jiangsu Province, China, or in a climate-controlled room under long-day conditions with a photocycle of 14 hr of light (32 °C) and 10 hr of darkness (28 °C), at 50% humidity. The N. benthamiana plants were grown in a growth chamber at 22 °C with long-day conditions (16 hr light/8 hr dark). To determine whether the chloroembryo of osnf-yb7 is light-dependent, the emerging panicles of the mutants were covered with aluminum foil until seed maturation. The lemmas of the WT were carefully removed with forceps at 1–2 DAF to expose the WT embryo to light in a climate-controlled room.

Vector construction and plant transformation

To generate the overexpression plants, full-length coding sequences (CDSs) of OsNF-YB7 or OsGLK1 (the stop codon removed) were cloned into pCAMBIA1300-Flag or pUN1301-GFP under the control of ubiquitin promoter, using the ClonExpress II One Step Cloning Kit (Vazyme). The constructs were transformed into rice calli through an Agrobacterium‐mediated strategy, as described previously (Chen et al., 2016). The higher order mutants of OsNF-YB7, OsGLK1 and OsGLK2 mutants, were generated using a previously described method for multiple gene editing in ZH11 (Cheng et al., 2021). The primers used for vector construction are listed in Supplementary file 1.

RNA extraction and RT-qPCR

Embryos or seedlings were collected by ﬂash freezing in liquid nitrogen and stored at −80 °C until processing. Samples were finely ground using a mortar and pestle with liquid nitrogen. Total RNA was isolated using the RNA-easy Isolation Reagent (Vazyme, R701-01). The experiments were performed with at least three biological replicates. The relative expression levels of the tested genes were normalized using the rice Actin gene and calculated by the 2ΔCt method. The primers used for the RT-qPCR are listed in Supplementary file 1.

RNA-sequencing and data analysis

Ten-DAF-old embryos of the Wild-type (WT) and OsGLK1-OX were used for RNA-sequencing. Two biological replicates were set. RNA extraction, library preparation, and high-throughput sequencing of the collected samples were outsourced to BGI Genomics Co., Ltd., Shenzhen, China. CLC Genomics Workbench 12.0 software was used for RNA-seq data analysis, as previously reported (Xu et al., 2021). The thresholds of fold change >2 and Bonferroni-corrected FDR <0.05 were used for defining a DEG. The previously generated ChIP-seq data of OsNF-YB7 and OsGLK1 (Guo et al., 2022; Tu et al., 2022) were reanalyzed using CLC Genomics Workbench 12.0 software for peak calling. Enriched motifs were identified by the online tool MEME-ChIP (https://meme-suite.org/meme/tools/meme-chip) with default parameter set. The software Mapman was used for pathway analysis (Usadel et al., 2009). The online tool AgriGO 2.0 was used for GO analysis (Du et al., 2010). The Venn diagrams were drawn using an online tool (https://bioinfogp.cnb.csic.es/tools/venny/index.html). TBtools was used for heat map generation (Chen et al., 2020).

Dual luciferase reporter (DLR) assays

The OsGLK1, OsLHCB4 (2 kb upstream of translation start site), and OsPORA (1.5 kb upstream of translation start site) promoter sequences were amplified from ZH11 genomic DNA and cloned into the vector pGreenII 0800-LUC (Hellens et al., 2005), as reporters; the 35 S::OsGLK1 and 35 S::OsNF-YB7 constructs were used as effectors. The reporters and effectors were transfected into rice etiolated protoplasts in different combinations and incubated overnight. Fireﬂy LUC and Renilla luciferase (REN) activities were measured using the DLR Assay Kit (Vazyme), following the manufacturer’s instructions, and the LUC:REN ratios were calculated for analysis. The primers used for generating these constructs are listed in Supplementary file 1.

Chromatin immunoprecipitation (ChIP) assays

The CDS (stop codons removed) of OsGLK1 was cloned into the pJIT163-GFP vector driven by a 35 S promoter to generate the OsGLK1-GFP construct. Protoplasts transformed with OsGLK1-GFP and OsNF-YB7-Flag transgenic lines were used for the ChIP assays, in accordance with a previously described method (Zhao et al., 2020). Briefly, protoplasts or 0.2 g of 14-day-old seedlings were harvested and crosslinked with 1% formaldehyde for 15 min, followed by neutralization using 0.125 M glycine for an additional 5 min. The seedlings were then ground into powder in liquid nitrogen. The nuclei were isolated and lysed using Buffer S (50 mM HEPES-KOH [pH 7.5], 150 mM NaCl, 1 mM ethylenediaminetetraacetic acid [EDTA], 1% Triton X-100, 0.1% sodium deoxycholate, 1% SDS) and Buffer F (50 mM HEPES-KOH [pH 7.5], 150 mM NaCl, 1 mM EDTA, 1% Triton X-100, 0.1% sodium deoxycholate). The chromatin was then sonicated with the segment size ranging from 200 to 600 bp. The lysates were then immunoprecipitated by anti-GFP (abcam no. ab290) and anti-Flag (Sigma no. F3165) antibodies, respectively. Immunocomplexes were washed with low-salt ChIP buffer (50 mM HEPES-KOH, 150 mM NaCl, 1 mM EDTA, 1% Triton X-100, 0.1% sodium deoxycholate, 0.1% SDS), high-salt ChIP buffer (low-salt ChIP buffer but replacing 150 mM NaCl with 350 mM NaCl), ChIP wash buffer (10 mM Tris-HCl pH 8.0, 250 mM LiCl, 0.5% NP-40, 1 mM EDTA, 0.1% sodium deoxycholate), and TE buffer (10 mM Tris-HCl, pH 8.0, and 1 mM EDTA). The protein-DNA complexes were eluted from beads using ChIP Elution buffer (50 mM Tris-HCl pH 7.5, 10 mM EDTA, 1% SDS) for 15 min at 65 °C and the crosslinking was reversed by incubation overnight with proteinase K. The fragment DNA was extracted with phenol:chloroform:isoamyl alcohol (25:24:1), precipitated with ethanol, and resuspended in TE buffer. The immunoprecipitated DNA was used as a template for qPCR. The primers used here are listed in Supplementary file 1.

Yeast two-hybrid assays

The CDSs of OsGLK1/2 and OsNF-YB7 were cloned into pGADT7 and pGBKT7, respectively. The constructs were cotransformed into yeast strain AH109 using Frozen-EZ Yeast Transformation II kit (Zymo), in accordance with the manufacturer’s instructions. The empty pGADT7 and pGBKT7 vectors were cotransformed in parallel as negative controls. The transformants were first selected on synthetic dropout medium (SD/-Leu-Trp) plates. We tested protein-protein interactions using selective SD/-Leu-Trp-His dropout medium. Interactions were observed after 3 days of incubation at 28 °C. The primers used for generating these constructs are listed in Supplementary file 1.

Split complementary LUC assays

Split complementary LUC assays were performed as previously described (Niu et al., 2020). The CDSs of OsGLK1/2 and OsNF-YB7 were cloned into JW771 and JW772 vectors to generate nLUC-OsGLK1/2 and cLUC-OsNF-YB7, respectively. The constructs were introduced into Agrobacterium tumefaciens strain GV3101 and then co-infiltrated into N. benthamiana leaves, and the LUC activities were analyzed after 48 hr of infiltration using Tanon Imaging System (5200 Multi; Tanon). The primers used for vector construction are shown in Supplementary file 1.

Bimolecular fluorescence complementation assays

The CDSs of OsNF-YB7 and OsGLK1 were cloned into the vector pSPYNE (nYFP) and pSPYCE (cYFP), respectively. The prepared plasmids were transformed into Agrobacterium strain GV3101, and the indicated transformant pairs were infiltrated into N. benthamiana leaves. Forty-eight hours after infiltration, the fluorescence signal of yellow fluorescent protein (YFP) was observed with a confocal laser scanning microscope (CLSM) (Carl Zeiss, LSM 710). Images were captured at 514 nm laser excitation and 519–620 nm emission for YFP. The primers used for vector construction are shown in Supplementary file 1.

Co-immunoprecipitation (Co-IP) assays

Co-IP assays were performed using rice protoplast as described previously (Zhang et al., 2011). The CDSs (stop codons removed) of OsGLK1 and OsNF-YB7 were cloned into the vectors pUC19-35S-FLAG-RBS and pJIT163-GFP driven by a 35 S promoter, respectively. Ten micrograms of plasmid DNA (OsGLK1-GFP, GFP, and OsNF-YB7-Flag) was transformed or cotransformed into 200 μl of protoplasts and incubated in WI buffer (0.5 M mannitol, 20 mM KCl, and 4 mM MES at pH 5.7) for 12 hr. The protoplasts were collected and lysed in 500 μl of lysis buffer (50 mM Tris-HCl, 150 mM NaCl, 5 mM EDTA [pH 8.0], 1% NP-40, 0.1 mM PMSF). The extracts were incubated with GFP-Trap agarose beads at 4 °C for 3 hr and washed three times with washing buffer. Samples were boiled in SDS protein loading buffer. Immunoblots were detected using corresponding primary antibodies (anti-GFP, ABclonal no. AE012; anti-Flag, Sigma no. F3165). The primers used for vector construction are shown in Supplementary file 1.

Chl measurement and confocal imaging

Approximately one hundred micrograms of embryos of the indicated genotypes were extracted in 1 ml of 100% dimethyl sulfoxide (DMSO) with incubation at 65 °C for 1 hr. Then, the absorbance values at wavelengths of 648.2 and 664.9 nm were measured by spectrophotometry and total Chl content was calculated as reported previously (Barnes et al., 1992). Chl autoﬂuorescence signal was detected by CLSM (Carl Zeiss, LSM 710), with excitation and emission wavelengths of 633 and 625–730 nm, respectively.

Electrophoretic mobility-shift assays (EMSAs)

OsNF-YB7 CDS was amplified by PCR and cloned into pET-28a vector to generate the OsNF-YB7-His construct. The full-length CDS of OsGLK1 was cloned into pMAL-c5X vector to generate OsGLK1-MBP construct. All constructs were transformed into E. coli strain BL21 to produce recombinant proteins. The promoter subfragments of OsPORA (42 bp, from –278––237) and OsLHCB4 (38 bp, from –259––222) were PCR amplified and labeled with biotin at the 3’ hydroxyl end of the double strands using EMSA Probe Biotin Labeling Kit (Beyotime, GS008). EMSA was performed using EMSA/Gel-Shift kit (Beyotime, GS009), in accordance with the manufacturer’s instructions. The labeled probes were detected in accordance with the instructions provided with the EMSA/Gel-Shift kit. All oligonucleotides used to generate the biotin-labeled probes are listed in Supplementary file 1.

Transmission electron microscopy (TEM)

TEM analysis was performed as described previously (Cheng et al., 2021). Brieﬂy, embryos of WT and osnf-yb7 were fixed overnight at 4 °C in 2.5% glutaraldehyde and 0.1 M PBS. The samples were subsequently washed three times with 0.1 M PBS and then fixed with 1% osmic acid for 4 hr. The samples were dehydrated in a series of ethanol and embedded in acrylic resin at 37 °C for 12 hr and at 60 °C for 48 hr. The samples were sectioned at 100 nm and observed by TEM (TECNAI 12).

Funding Information

This paper was supported by the following grants:

http://dx.doi.org/10.13039/501100001809 National Natural Science Foundation of China 32170344 to Chen Chen.

http://dx.doi.org/10.13039/501100001809 National Natural Science Foundation of China 32300689 to Zongju Yang.

Project of Zhongshan Biological Breeding Laboratory BM2022008-02 to Chen Chen.

Natural Science Foundation of the Higher Education Institutions of Jiangsu Province 21KJB21003 to Zongju Yang.

http://dx.doi.org/10.13039/501100002858 China Postdoctoral Science Foundation 2022M712702 to Zongju Yang.

http://dx.doi.org/10.13039/501100002949 Jiangsu Province Government JBGS001 to Chen Chen.

http://dx.doi.org/10.13039/501100012246 Priority Academic Program Development of Jiangsu Higher Education Institutions PAPD to Chen Chen.

Independent Scientific Research Project Funds of the Jiangsu Key Laboratory of Crop Genomics and Molecular Breeding PLR202101 to Chen Chen.

Acknowledgements

We thank Prof. Hengxiu Yu and Dr. Chao Zhang of Yangzhou University for kindly providing the positive and negative control vectors for dual luciferase reporter assays. This research was supported by grants from the National Natural Science Foundation of China (32170344, 32300689), the Project of Zhongshan Biological Breeding Laboratory (BM2022008-02), the Jiangsu Province Government (JBGS001), the Natural Science Foundation of the Higher Education Institutions of Jiangsu Province (21KJB21003), the China Postdoctoral Science Foundation (2022M712702), the Independent Scientific Research Project Funds of the Jiangsu Key Laboratory of Crop Genomics and Molecular Breeding (PLR202101), the Priority Academic Program Development of Jiangsu Higher Education Institutions (PAPD).

Additional information

Competing interests

Author contributions

Additional files

Supplementary file 1. Primers used in this study.

MDAR checklist

Data availability

Sequencing data have been deposited in NCBI BioProject under the accession code PRJNA998591.

The following dataset was generated:

Yang Z Niu B Chen C 2023 OsNF-YB7 inactivates OsGLK1 to prevent chlorophyll biosynthesis in rice embryo NCBI BioProject PRJNA998591

The following previously published datasets were used:

Zhong S Tu X 2022 GLK ChIP-seq and ATAC-seq data in Arabidopsis, Tobacco, Tomato, Rice and Maize NCBI Gene Expression Omnibus GSE220115

Guo F Bian H Zhang P Wu Y 2022 ChIP-seq of OsLEC1 NCBI Gene Expression Omnibus GSE179596

10.7554/eLife.96553.3.sa0
eLife assessment
Casati Paula Reviewing Editor Center of Photosynthetic And Biochemical Studies (CEFOBI) Argentina

Convincing
Important
This paper provides important insights into the role of rice OsNF-YB7, an ortholog of Arabidopsis LEC1, in chlorophyll biosynthesis, uncovering the genetic and molecular basis for negative regulation of chlorophyll production in the rice embryo. Mutational analysis, gene expression profiles and protein interaction combine for convincing evidence that OsNF-YB7 represses chlorophyll biosynthesis.

10.7554/eLife.96553.3.sa1
Reviewer #1 (Public review):
Reviewer
Summary:

This manuscript investigates the regulation of chlorophyll biosynthesis in rice embryos, focusing on the role of OsNF-YB7. The rigorous experimental approach, combining genetic, biochemical, and molecular analyses, provides a robust foundation for these findings. The research achieves its objectives, offering new insights into chlorophyll biosynthesis regulation, with the results convincingly supporting the authors' conclusions.

Strengths:

The major strengths include the detailed experimental design and the findings regarding OsNF-YB7's inhibitory role.

Weaknesses:

However, the manuscript's discussion on the practical implications for agriculture and the evolutionary analysis of regulatory mechanisms could be expanded.

10.7554/eLife.96553.3.sa2
Reviewer #2 (Public review):
Reviewer
Summary:

The authors set out to establish the role of the rice LEC1 homolog OsNF-YB7 in embryo development, especially as it pertains to the development of photosynthetic capacity, with chlorophyll production as a primary focus.

Strengths:

The results are well-supported and each approach used complements each other. There are no major questions left unanswered and the central hypothesis is addressed in every figure.

Weaknesses:

There are a handful of sections which could use clarifying for readers, but overall this is a solidly composed manuscript.

The authors clearly achieved their aims; the results compellingly establish a disparity between how this system operates in rice and Arabidopsis. Conclusions are thoroughly supported by the provided data and interpretations. This work will force a reconsideration of the value of Arabidopsis as a model organism for embryo chlorophyll biosynthesis and possibly photosynthesis during embryo maturation more broadly, as rice is a major crop organism and it very clearly does not follow the Arabidopsis model. It will thus be useful to carry out similar tests in other organisms rather than relying on Arabidopsis and attempt to more fully establish the regulatory mechanism in rice.

10.7554/eLife.96553.3.sa3
Reviewer #3 (Public review):
Reviewer
Summary:

In this study, the authors set out to understand the mechanisms behind chlorophyll biosynthesis in rice, focusing in particular on the role of OsNF-YB7, an ortholog of Arabidopsis LEC1, which is a positive regulator of chlorophyll (Chl) biosynthesis in Arabidopsis. They showed that OsNF-YB7 loss-of-function mutants in rice have chlorophyll-rich embryos, in contrast to Arabidopsis LEC1 loss-of-function mutants. This contrasting phenotype led the authors to carry out extensive molecular studies on OsNF-YB7, including in vitro and in vivo protein interaction studies, gene expression profiling and protein-DNA interaction assays. The evidence provided well supported the core arguments of the authors, emphasising that OsNF-YB7 is a negative regulator of Chl biosynthesis in rice embryos by mediating the expression of OsGLK1, a transcription factor that regulates downstream Chl biosynthesis genes. In addition, they showed that OsNF-YB7 interacts with OsGLK1 to negatively regulate the expression of OsGLK1, demonstrating the broad involvement of OsNF-YB7 in rice Chl biosynthetic pathways.

Strengths:

This study clearly demonstrated how OsNF-YB7 regulates its downstream pathways using several in vitro and in vivo approaches. For example, gene expression analysis of OsNF-YB7 loss-of-function and gain-of-function mutants revealed the expression of selected downstream chl biosynthetic genes. This was further validated by EMSA on the gel. The authors also confirmed this using luciferase assays in rice protoplasts. These approaches were used again to show how the interaction of OsNF-YB7 and OsGLK1 regulates downstream genes. The main idea of this study is very well supported by the results and data.

Weaknesses:

It would be interesting to see how two similar genes have come to play opposite roles in Arabidopsis and rice. Interspecies complementation might help to understand this point.

10.7554/eLife.96553.3.sa4
Author response
Yang Zongju Author Yangzhou University Yangzhou China

Bai Tianqi Author Yangzhou University Yangzhou China

E Zhiguo Author China National Rice Research Institute Hangzhou China

Niu Baixiao Author Yangzhou University Yangzhou China

Chen Chen Author Yangzhou University Yangzhou China

The following is the authors’ response to the original reviews.

eLife assessment

This is an important study on the regulation of chlorophyll biosynthesis in rice embryos. It provides insights into the genetic and molecular interactions that underlie chlorophyll accumulation, highlighting the inhibition of OsGLK1 by OsNF-YB7 and the broader implications for understanding chloroplast development and seed maturation in angiosperms. The results presented, including mutation analysis, gene expression profiles, and protein interaction studies, provide convincing evidence for the function of OsNF-YB7 as a repressor in the chlorophyll biosynthesis pathway.

Thank you very much for your positive assessment of our manuscript. We have carefully revised the manuscript according to the reviewers’ valuable suggestions and comments. For more details, please see the point-to-point response to the reviewers below.

Public Reviews:

Reviewer #1 (Public Review):

Summary:

This manuscript investigates the regulation of chlorophyll biosynthesis in rice embryos, focusing on the role of OsNF-YB7. The rigorous experimental approach, combining genetic, biochemical, and molecular analyses, provides a robust foundation for these findings. The research achieves its objectives, offering new insights into chlorophyll biosynthesis regulation, with the results convincingly supporting the authors' conclusions.

Strengths:

The major strengths include the detailed experimental design and the findings regarding OsNF-YB7's inhibitory role.

Weaknesses:

However, the manuscript's discussion on the practical implications for agriculture and the evolutionary analysis of regulatory mechanisms could be expanded.

Thank you for your insightful comments and suggestions. In the revised manuscript, we discussed the potential application of the chlorophyllous embryo (please see line 270-274). The presence of chlorophyll in the embryo facilitates photosynthesis at early developmental stages, potentially leading to improved seedling growth and vigor (Smolikova and Medvedev, 2016). In crops such as soybean and canola, green embryo is considered as a valuable trait due to its association with enhanced photosynthetic capacity, which consequently promotes fatty acid biosynthesis (Ruuska et al., 2004). However, chlorophyll degradation must be carefully managed during seed maturation to avoid negative effects on seed viability and meal quality (Chung et al., 2006). Interestingly, the green embryo of lotus (Nelumbo nucifera) is widely used as a food ingredient in Asian, Australia, and North America. It is employed in herbal medicine to treat nervous disorders, insomnia, and other conditions (Zhu et al., 2017; Ha et al., 2022), highlighting the significant potential value of the green embryo.

In many chloroembryophytes, such as Arabidopsis, the embryo occupies a large proportion of the seed. From an evolutionary perspective, the presence of chlorophyll in the embryo may promote adaptation in such chloroembryophytes because more reserves can be accumulated in the seed through active photosynthesis, better supporting the embryo development and subsequent seedling growth (Sela et al., 2020). On the other hand, some leucoembryophytes, such as rice, have persistent endosperm rich in storage reserves to nourish embryo development (Liu et al., 2022). Gaining the ability to accumulate chlorophyll in the embryo is unnecessary for such species. In agreement with this hypothesis, cholorophyllous embryos are more prevalent in non-endospermous seeds (Dahlgren, 1980). However, we would like to emphasize that the evolutionary force driving the divergence of chloroembryophytes and leucoembryophytes is currently almost completely unknown and deserves in-depth investigation in the future. We discussed the possible evolution of the ability to accumulate chlorophyll in the embryo, please find the details in Line 276-295.

Reviewer #2 (Public Review):

Summary:

The authors set out to establish the role of the rice LEC1 homolog OsNF-YB7 in embryo development, especially as it pertains to the development of photosynthetic capacity, with chlorophyll production as a primary focus.

Strengths:

The results are well-supported and each approach used complements each other. There are no major questions left unanswered and the central hypothesis is addressed in every figure.

Weaknesses:

There are a handful of sections that could use clarifying for readers, but overall this is a solidly composed manuscript.

The authors clearly achieved their aims; the results compellingly establish a disparity between how this system operates in rice and Arabidopsis. Conclusions are thoroughly supported by the provided data and interpretations. This work will force a reconsideration of the value of Arabidopsis as a model organism for embryo chlorophyll biosynthesis and possibly photosynthesis during embryo maturation more broadly, as rice is a major crop organism and it very clearly does not follow the Arabidopsis model. It will thus be useful to carry out similar tests in other organisms rather than relying on Arabidopsis and attempting to more fully establish the regulatory mechanism in rice.

Thank you very much for your positive comments. We have carefully revised the manuscript according to your and the other reviewers’ comments and suggestions. Particularly, we emphasized the necessary to carry out similar tests in other organisms rather than relying on Arabidopsis to better understand the regulatory mechanism in rice.

Reviewer #3 (Public Review):

Summary:

In this study, the authors set out to understand the mechanisms behind chlorophyll biosynthesis in rice, focusing in particular on the role of OsNF-YB7, an ortholog of Arabidopsis LEC1, which is a positive regulator of chlorophyll (Chl) biosynthesis in Arabidopsis. They showed that OsNF-YB7 loss-of-function mutants in rice have chlorophyll-rich embryos, in contrast to Arabidopsis LEC1 loss-of-function mutants. This contrasting phenotype led the authors to carry out extensive molecular studies on OsNF-YB7, including in vitro and in vivo protein interaction studies, gene expression profiling, and protein-DNA interaction assays. The evidence provided well supported the core arguments of the authors, emphasising that OsNF-YB7 is a negative regulator of Chl biosynthesis in rice embryos by mediating the expression of OsGLK1, a transcription factor that regulates downstream Chl biosynthesis genes. In addition, they showed that OsNF-YB7 interacts with OsGLK1 to negatively regulate the expression of OsGLK1, demonstrating the broad involvement of OsNF-YB7 in rice Chl biosynthetic pathways.

Strengths:

This study clearly demonstrated how OsNF-YB7 regulates its downstream pathways using several in vitro and in vivo approaches. For example, gene expression analysis of OsNF-YB7 loss-of-function and gain-of-function mutants revealed the expression of selected downstream chl biosynthetic genes. This was further validated by EMSA on the gel. The authors also confirmed this using luciferase assays in rice protoplasts. These approaches were used again to show how the interaction of OsNF-YB7 and OsGLK1 regulates downstream genes. The main idea of this study is very well supported by the results and data.

Weaknesses:

From an evolutionary perspective, it is interesting to see how two similar genes have come to play opposite roles in Arabidopsis and rice. It would have been more interesting if the authors had carried out a cross-species analysis of AtLEC1 and OsNF-YB7. For example, overexpressing AtLEC1 in an osnf-yb7 mutant to see if the phenotype is restored or enhanced. Such an approach would help us understand how two similar proteins can play opposite roles in the same mechanism within their respective plant species.

We appreciate your insightful comments and suggestions. It is a very interesting question whether AtLEC1 can fully restore osnf-yb7, given the possible functional divergence between the genes in terms of regulation of chlorophyll biosynthesis in the embryo. We have previously expressed OsNF-YB7 in the lec1-1 background in Arabidopsis, driven by the native promoter of LEC1 (Niu et al., 2021). We found that OsNF-YB7 could almost completely rescue the embryo defects in Arabidopsis, indicating that OsNF-YB7 plays a resemble role in rice as the LEC1 does in Arabidopsis (Niu et al., 2021). We sought to determine whether AtLEC1 can complement the chlorophyll defect in osnf-yb7. However, given the fact that osnf-yb7 shows severe callus induction defect, which is not surprising, because many studies have shown that LEC1 is indispensable for somatic embryo development in various plant species, we are struggling to obtain the genetic materials for analysis. We have to transform OsNF-YB7pro::AtLEC1 into the WT background first, and then cross the transformant with the osnf-yb7 mutant. This is a time-consuming process in rice, but hopefully we will able to isolate a line expressing OsNF-YB7pro::AtLEC1 in the osnf-yb7 background from the resulting segregating population.

Recommendations for the authors:

Reviewer #1 (Recommendations For The Authors):

A minor comment regarding the chlorophyll contents quantification in the study. Line 87: "The results showed that WT had an achlorophyllous embryo throughout embryonic development,...." In the TEM result, chloroplast was not observed in the WT embryo sections, indicating a lack of chlorophyll-containing structures, contrary to what was found in the osnf-yb7 embryos where chloroplasts were observed.

The authors stated that the embryo morphologies and Chl autofluorescence data showed that WT had an achlorophyllous embryo throughout embryonic development. However, the quantification of Chl levels in Figure 1D and Figure 4C showed that WT does produce some chlorophylls, albeit at lower levels than osnf-yb7 or OSGLK-OX embryos (WT values in the two figures are slightly different). This discrepancy warrants clarification to ensure consistency and accuracy in the manuscript's findings.

We re-evaluated the Chl content in the embryos of WT and OsGLK1-OX mature seeds. The result confirmed our previous finding that WT embryos produce a small amount of chlorophyll (please see the updated Fig. 4C). Notably, we observed that the dark-grown etiolated plants still have measurable chlorophyll content as reported in many studies (for example, Wang et al., 2017; Yoo et al., 2019), suggesting that there is potential bias in measuring chlorophyll content using an absorbance-based approach. We assume this possibly explains the concern you have raised.

Reviewer #2 (Recommendations For The Authors):

Mild editing for grammar is needed throughout, e.g. line 73, "It is still a mysterious why plant species".

We have carefully edited the grammar.

As a minor point, the placement of figure panels, such as in Figure 1, is not always intuitive.

Thank you for your suggestion. This figure has been revised as suggested. Please see the updated Fig. 1.

What is the significance of the two GFP mutants in Figures 2C and 2D? Is one of those the mislabeled Flag mutant?

The lines showed in Fig. 2C and D were not mislabeled. They were two independent transgenic events, both of which showed that OsNF-YB7 inhibited the expression of OsPORA and OsLHCB4 in rice. The transgenic lines overexpressing OsNF-YB7 tagging with the 3× Flag (NF-YB7-Flag) were also used for this experiment. In agreement, OsPORA and OsLHCB4 were significantly downregulated in the three independent NF-YB7-Flag lines (Fig. S4C), confirming the results showed in Fig. 2C and D.

In Figures 2G and 2H, what is that enormous band at the bottom of the gel?

The bands at the bottom of the gel were free probes. We indicated this in the revised figure.

Not until the Materials and Methods section did I realize that any of this study was being done in tobacco; the Introduction implies it's rice vs. Arabidopsis and it might be a good idea to mention the organism of study somewhere before Figure 6.

We apologize for any confusion caused by our previous writing. While the majority of this study was performed with rice plants or protoplasts, the split complementary LUC assays and BiFC assays were performed with tobacco. We have specified these in the revised manuscript as suggested.

Reviewer #3 (Recommendations For The Authors):

It would be nice if the author could show what the phenotype is in AtLEC1 OX in osnf-yb7 and also OsNF-YB7 OX in atlec1 mutants.

Thank you for your suggestion. We have previously expressed OsNF-YB7 in the lec1-1 background of Arabidopsis, driven by the native promoter of Arabidopsis LEC1 (Niu et al., 2021). Since OsNF-YB7 could rescue the embryo morphogenesis defects in Arabidopsis (Niu et al., 2021), we assumed that OsNF-YB7 plays a similar role in rice as the LEC1 does in Arabidopsis. However, it remains unknown whether expression of LEC1 in osnf-yb7 may restore the chlorophyllous embryo phenotype in rice. As the generation of genetic material is time-consuming, and especially given the fact that osnf-yb7 has a severe callus induction defect, we are struggling to obtain the complementary line for analysis. We have to transform OsNF-YB7pro::AtLEC1 in a WT background first, and then cross the transformant with the osnf-yb7 mutant. Hopefully, we will be able to isolate a line expressing OsNF-YB7pro::AtLEC1 in osnf-yb7 background, from the derived segregating population. We discussed the reviewer’s concern in the revised manuscript, please see Line 369-376.

Line 46, I think it is vague to mention that 'Like most plant species'. Some species might have different copy numbers, for example, a single GLK in liverwort M. polymorpha.

The statement has been revised. Please see Line 46.

Figures 2F and 5B, why was only one promoter region used for OsLHCB4? It would be better to have more regions like OsPORA.

Thank you for your comments. Here, we have examined more promoter regions (P1, P2 and P3) in the revised manuscript as suggested, among which, the previously selected promoter region (P3) contains both the G-box and CCAATC motifs that can be potentially recognized by GLK1. Consistent to our previous report, the results showed that OsNF-YB7 (left) and OsGLK1 (right) were associated with the P3 region, but showed no significant differences in the other probes. Please see the results in Fig. 2F and Fig. 5B of the revised manuscript.

Legend of Figures 2G, H, OsPORA (I), and OsLHCB (J) should be (G) and (H) respectively.

Corrected.

References

Chung, D.W., Pruzinska, A., Hortensteiner, S., and Ort, D.R. (2006). The role of pheophorbide a oxygenase expression and activity in the canola green seed problem. Plant Physiol 142, 88-97.

Ha, T., Kim, M.S., Kang, B., Kim, K., Hong, S.S., Kang, T., Woo, J., Han, K., Oh, U., Choi, C.W., and Hong, G.S. (2022). Lotus Seed Green Embryo Extract and a Purified Glycosyloxyflavone Constituent, Narcissoside, Activate TRPV1 Channels in Dorsal Root Ganglion Sensory Neurons. J Agric Food Chem 70, 3969-3978.

Liu, J., Wu, M.W., and Liu, C.M. (2022). Cereal Endosperms: Development and Storage Product Accumulation. Annu Rev Plant Biol 73, 255-291.

Niu, B., Zhang, Z., Zhang, J., Zhou, Y., and Chen, C. (2021). The rice LEC1-like transcription factor OsNF-YB9 interacts with SPK, an endosperm-specific sucrose synthase protein kinase, and functions in seed development. Plant J 106, 1233-1246.

Ruuska, S.A., Schwender, J., and Ohlrogge, J.B. (2004). The capacity of green oilseeds to utilize photosynthesis to drive biosynthetic processes. Plant Physiol 136, 2700-2709.

Sela, A., Piskurewicz, U., Megies, C., Mene-Saffrane, L., Finazzi, G., and Lopez-Molina, L. (2020). Embryonic Photosynthesis Affects Post-Germination Plant Growth. Plant Physiol 182, 2166-2181.

Smolikova, G.N., and Medvedev, S.S. (2016). Photosynthesis in the seeds of chloroembryophytes. Russ J Plant Physl+ 63, 1-12.

Wang, Z., Hong, X., Hu, K., Wang, Y., Wang, X., Du, S., Li, Y., Hu, D., Cheng, K., An, B., and Li, Y. (2017). Impaired Magnesium Protoporphyrin IX Methyltransferase (ChlM) Impedes Chlorophyll Synthesis and Plant Growth in Rice. Front Plant Sci 8, 1694.

Yoo, C.Y., Pasoreck, E.K., Wang, H., Cao, J., Blaha, G.M., Weigel, D., and Chen, M. (2019). Phytochrome activates the plastid-encoded RNA polymerase for chloroplast biogenesis via nucleus-to-plastid signaling. Nat Commun 10, 2629.

Zhu, M., Liu, T., Zhang, C., and Guo, M. (2017). Flavonoids of Lotus (Nelumbo nucifera) Seed Embryos and Their Antioxidant Potential. J Food Sci 82, 1834-1841.

No competing interests declared.

Conceptualization, Data curation, Funding acquisition, Validation, Investigation, Visualization, Methodology, Writing - original draft.

Validation, Investigation, Visualization, Methodology.

Software, Investigation.

Investigation.

Conceptualization, Supervision, Funding acquisition, Visualization, Writing - original draft, Project administration, Writing - review and editing.
==== Refs
References

Bai B Lu N Li Y Guo S Yin H He Y Sun W Li W Xie X 2019 OsBBX14 promotes photomorphogenesis in rice by activating OsHY5L1 expression under blue light conditions Plant Science 284 192 202 10.1016/j.plantsci.2019.04.017 31084872
Barnes JD Balaguer L Manrique E Elvira S Davison AW 1992 A reappraisal of the use of DMSO for the extraction and determination of chlorophylls A and b in lichens and higher plants Environmental and Experimental Botany 32 85 100 10.1016/0098-8472(92)90034-Y
Bogorad L 1950 Factors associated with the synthesis of chlorophyll in the dark in seedlings of pinus jeffreyi Botanical Gazette 111 221 241 10.1086/335591
Boulard C Thévenin J Tranquet O Laporte V Lepiniec L Dubreucq B 2018 LEC1 (NF-YB9) directly interacts with LEC2 to control gene expression in seed Biochimica et Biophysica Acta. Gene Regulatory Mechanisms 1861 443 450 10.1016/j.bbagrm.2018.03.005 29580949
Chaves-Sanjuan A Gnesutta N Gobbini A Martignago D Bernardini A Fornara F Mantovani R Nardini M 2021 Structural determinants for NF-Y subunit organization and NF-Y/DNA association in plants The Plant Journal 105 49 61 10.1111/tpj.15038 33098724
Chen C Begcy K Liu K Folsom JJ Wang Z Zhang C Walia H 2016 Heat stress yields a unique MADS box transcription factor in determining seed size and thermal sensitivity Plant Physiology 171 606 622 10.1104/pp.15.01992 26936896
Chen C Chen H Zhang Y Thomas HR Frank MH He Y Xia R 2020 TBtools: an integrative toolkit developed for interactive analyses of big biological data Molecular Plant 13 1194 1202 10.1016/j.molp.2020.06.009 32585190
Cheng X Pan M Zhiguo E Zhou Y Niu B Chen C 2021 The maternally expressed polycomb group gene OsEMF2a is essential for endosperm cellularization and imprinting in rice Plant Communications 2 100092 10.1016/j.xplc.2020.100092 33511344
Chung DW Pruzinská A Hörtensteiner S Ort DR 2006 The role of pheophorbide a oxygenase expression and activity in the canola green seed problem Plant Physiology 142 88 97 10.1104/pp.106.084483 16844830
Dahlgren R 1980 The taxonomic significance of chlorophyllous embryos in angiosperm seeds Botaniska Notiser 133 337 341
Delmas F Sankaranarayanan S Deb S Widdup E Bournonville C Bollier N Northey JGB McCourt P Samuel MA 2013 ABI3 controls embryo degreening through Mendel’s I locus PNAS 110 E3888 E3894 10.1073/pnas.1308114110 24043799
Du Z Zhou X Ling Y Zhang Z Su Z 2010 agriGO: a GO analysis toolkit for the agricultural community Nucleic Acids Research 38 W64 W70 10.1093/nar/gkq310 20435677
Fatihi A Boulard C Bouyer D Baud S Dubreucq B Lepiniec L 2016 Deciphering and modifying LAFL transcriptional regulatory network in seed for improving yield and quality of storage compounds Plant Science 250 198 204 10.1016/j.plantsci.2016.06.013 27457996
Fitter DW Martin DJ Copley MJ Scotland RW Langdale JA 2002 GLK gene pairs regulate chloroplast development in diverse plant species The Plant Journal 31 713 727 10.1046/j.1365-313x.2002.01390.x 12220263
Gnesutta N Saad D Chaves-Sanjuan A Mantovani R Nardini M 2017 Crystal structure of the Arabidopsis thaliana L1L/NF-YC3 histone-fold dimer reveals specificities of the LEC1 Family of NF-Y subunits in plants Molecular Plant 10 645 648 10.1016/j.molp.2016.11.006 27871811
Guo F Zhang P Wu Y Lian G Yang Z Liu W Buerte B Zhou C Zhang W Li D Han N Tong Z Zhu M Xu L Chen M Bian H 2022 Rice leafy cotyledon1 hinders embryo greening during the seed development Frontiers in Plant Science 13 887980 10.3389/fpls.2022.887980 35620685
Hellens RP Allan AC Friel EN Bolitho K Grafton K Templeton MD Karunairetnam S Gleave AP Laing WA 2005 Transient expression vectors for functional genomics, quantification of promoter activity and RNA silencing in plants Plant Methods 1 13 10.1186/1746-4811-1-13 16359558
Huang M Hu Y Liu X Li Y Hou X 2015a Arabidopsis leafy cotyledon1 controls cell fate determination during post-embryonic development Frontiers in Plant Science 6 955 10.3389/fpls.2015.00955 26579186
Huang M Hu Y Liu X Li Y Hou X 2015b Arabidopsis leafy cotyledon1 mediates postembryonic development via interacting with phytochrome-interacting factor4 The Plant Cell 27 3099 3111 10.1105/tpc.15.00750 26566918
Hudson D Guevara DR Hand AJ Xu Z Hao L Chen X Zhu T Bi YM Rothstein SJ 2013 Rice cytokinin GATA transcription Factor1 regulates chloroplast development and plant architecture Plant Physiology 162 132 144 10.1104/pp.113.217265 23548780
Ito Y Thirumurugan T Serizawa A Hiratsu K Ohme-Takagi M Kurata N 2011 Aberrant vegetative and reproductive development by overexpression and lethality by silencing of OsHAP3E in rice Plant Science 181 105 110 10.1016/j.plantsci.2011.04.009 21683874
Jarvis P López-Juez E 2013 Biogenesis and homeostasis of chloroplasts and other plastids Nature Reviews. Molecular Cell Biology 14 787 802 10.1038/nrm3702 24263360
Jo L Pelletier JM Hsu S-W Baden R Goldberg RB Harada JJ 2020 Combinatorial interactions of the LEC1 transcription factor specify diverse developmental programs during soybean seed development PNAS 117 1223 1232 10.1073/pnas.1918441117 31892538
Kobayashi K Baba S Obayashi T Sato M Toyooka K Keränen M Aro E-M Fukaki H Ohta H Sugimoto K Masuda T 2012 Regulation of root greening by light and auxin/cytokinin signaling in Arabidopsis The Plant Cell 24 1081 1095 10.1105/tpc.111.092254 22415275
Laloum T De Mita S Gamas P Baudin M Niebel A 2013 CCAAT-box binding transcription factors in plants: Y so many? Trends in Plant Science 18 157 166 10.1016/j.tplants.2012.07.004 22939172
Li Z Wu S Chen J Wang X Gao J Ren G Kuai B 2017 NYEs/SGRs‐mediated chlorophyll degradation is critical for detoxification during seed maturation in Arabidopsis The Plant Journal 92 650 661 10.1111/tpj.13710 28873256
Li Z Mo W Jia L Xu YC Tang W Yang W Guo YL Lin R 2019 Rice fluorescent1 is involved in the regulation of chlorophyll Plant and Cell Physiology 60 2307 2318 10.1093/pcp/pcz129 31290959
Li X Wang P Li J Wei S Yan Y Yang J Zhao M Langdale JA Zhou W 2020 Maize GOLDEN2-LIKE genes enhance biomass and grain yields in rice by improving photosynthesis and reducing photoinhibition Communications Biology 3 151 10.1038/s42003-020-0887-3 32238902
Liu H Wang X Ren K Li K Wei M Wang W Sheng X 2017 Light deprivation-induced inhibition of chloroplast biogenesis does not arrest embryo morphogenesis but strongly reduces the accumulation of storage reserves during embryo maturation in arabidopsis Frontiers in Plant Science 8 1 14 10.3389/fpls.2017.01287 28220127
Liu J Wu M-W Liu C-M 2022 Cereal endosperms: development and storage product accumulation Annual Review of Plant Biology 73 255 291 10.1146/annurev-arplant-070221-024405 35226815
Lotan T Ohto M Yee KM West MAL Lo R Kwong RW Yamagishi K Fischer RL Goldberg RB Harada JJ 1998 Arabidopsis leafy cotyledon1 is sufficient to induce embryo development in vegetative cells Cell 93 1195 1205 10.1016/s0092-8674(00)81463-4 9657152
Luo W Tan J Li T Feng Z Ding Z Xie X Chen Y Chen L Liu YG Zhu Q Guo J 2023 Overexpression of maize GOLDEN2 in rice and maize calli improves regeneration by activating chloroplast development Science China. Life Sciences 66 340 349 10.1007/s11427-022-2149-2 35982378
Meinke DW 1992 A homoeotic mutant of Arabidopsis thaliana with leafy cotyledons Science 258 1647 1650 10.1126/science.258.5088.1647 17742538
Meinke DW Franzmann LH Nickle TC Yeung EC 1994 Leafy cotyledon mutants of arabidopsis The Plant Cell 6 1049 1064 10.1105/tpc.6.8.1049 12244265
Myers J 1940 A study of the pigments produced in darkness by certain green algae Plant Physiology 15 575 588 10.1104/pp.15.4.575 16653658
Nakamura H Muramatsu M Hakata M Ueno O Nagamura Y Hirochika H Takano M Ichikawa H 2009 Ectopic overexpression of the transcription factor OsGLK1 induces chloroplast development in non-green rice cells Plant & Cell Physiology 50 1933 1949 10.1093/pcp/pcp138 19808806
Niu B Deng H Li T Sharma S Yun Q Li Q Zhiguo E Chen C 2020 OsbZIP76 interacts with OsNF-YBs and regulates endosperm cellularization in rice (Oryza sativa) Journal of Integrative Plant Biology 62 1983 1996 10.1111/jipb.12989 32621654
Niu B Yang Z Bai T Zhiguo E Xun X Yun Q Zhang J Zhang Z Lu X Li Q Liu QQ Chen C 2021a Rice Leafy Cotyledon1 Regulates Embryonic Envelope Development and Chlorophyll Biogenesis in Embryo Research Square 10.21203/rs.3.rs-888827/v1
Niu B Zhang Z Zhang J Zhou Y Chen C 2021b The rice LEC1-like transcription factor OsNF-YB9 interacts with SPK, an endosperm-specific sucrose synthase protein kinase, and functions in seed development The Plant Journal 106 1233 1246 10.1111/tpj.15230 33721364
Parcy F Valon C Kohara A Miséra S Giraudat J 1997 The abscisic acid-insensitive3, fusca3, and leafy cotyledon1 loci act in concert to control multiple aspects of arabidopsis seed development The Plant Cell 9 1265 1277 10.1105/tpc.9.8.1265 9286105
Pelletier JM Kwong RW Park S Le BH Baden R Cagliari A Hashimoto M Munoz MD Fischer RL Goldberg RB Harada JJ 2017 LEC1 sequentially regulates the transcription of genes involved in diverse developmental processes during seed development PNAS 114 E6710 E6719 10.1073/pnas.1707957114 28739919
Periasamy K Vivekanandan M 1981 Photosynthesis in the chloroembryo of cyamopsis tetragonaloba Taub Annals of Botany 47 793 797 10.1093/oxfordjournals.aob.a086078
Puthur JT Shackira AM Saradhi PP Bartels D 2013 Chloroembryos: A unique photosynthesis system Journal of Plant Physiology 170 1131 1138 10.1016/j.jplph.2013.04.011 23706538
Ranade SS Delhomme N García-Gil MR 2019 Global gene expression analysis in etiolated and de-etiolated seedlings in conifers PLOS ONE 14 e0219272 10.1371/journal.pone.0219272 31276530
Rossini L Cribb L Martin DJ Langdale JA 2001 The maize golden2 gene defines a novel class of transcriptional regulators in plants The Plant Cell 13 1231 1244 10.1105/tpc.13.5.1231 11340194
Ruuska SA Schwender J Ohlrogge JB 2004 The capacity of green oilseeds to utilize photosynthesis to drive biosynthetic processes Plant Physiology 136 2700 2709 10.1104/pp.104.047977 15347783
Sakuraba Y Kim EY Han SH Piao W An G Todaka D Yamaguchi-Shinozaki K Paek NC 2017 Rice Phytochrome-Interacting Factor-Like1 (OsPIL1) is involved in the promotion of chlorophyll biosynthesis through feed-forward regulatory loops Journal of Experimental Botany 68 4103 4114 10.1093/jxb/erx231 28922754
Sato Y Morita R Nishimura M Yamaguchi H Kusaba M 2007 Mendel’s green cotyledon gene encodes a positive regulator of the chlorophyll-degrading pathway PNAS 104 14169 14174 10.1073/pnas.0705521104 17709752
Sela A Piskurewicz U Megies C Mène-Saffrané L Finazzi G Lopez-Molina L 2020 Embryonic photosynthesis affects post-germination plant growth Plant Physiology 182 2166 2181 10.1104/pp.20.00043 32060052
Simkin AJ Faralli M Ramamoorthy S Lawson T 2020 Photosynthesis in non-foliar tissues: implications for yield The Plant Journal 101 1001 1015 10.1111/tpj.14633 31802560
Smolikova GN Medvedev SS 2016 Photosynthesis in the seeds of chloroembryophytes Russian Journal of Plant Physiology 63 1 12 10.1134/S1021443715060163
Tamai H Iwabuchi M Meshi T 2002 Arabidopsis GARP transcriptional activators interact with the Pro-rich activation domain shared by G-box-binding bZIP factors Plant & Cell Physiology 43 99 107 10.1093/pcp/pcf011 11828027
Tao Z Shen L Gu X Wang Y Yu H He Y 2017 Embryonic epigenetic reprogramming by a pioneer transcription factor in plants Nature 551 124 128 10.1038/nature24300 29072296
Thomas H Ougham H 2014 The stay-green trait Journal of Experimental Botany 65 3889 3900 10.1093/jxb/eru037 24600017
Tu X Ren S Shen W Li J Li Y Li C Li Y Zong Z Xie W Grierson D Fei Z Giovannoni J Li P Zhong S 2022 Limited conservation in cross-species comparison of GLK transcription factor binding suggested wide-spread cistrome divergence Nature Communications 13 7632 10.1038/s41467-022-35438-4 36494366
Usadel B Poree F Nagel A Lohse M Czedik-Eysenberg A Stitt M 2009 A guide to using MapMan to visualize and compare Omics data in plants: A case study in the crop species, Maize Plant, Cell & Environment 32 1211 1229 10.1111/j.1365-3040.2009.01978.x 19389052
Wang P Fouracre J Kelly S Karki S Gowik U Aubry S Shaw MK Westhoff P Slamet-Loedin IH Quick WP Hibberd JM Langdale JA 2013 Evolution of GOLDEN2-LIKE gene function in C(3) and C (4) plants Planta 237 481 495 10.1007/s00425-012-1754-3 22968911
Wang P Khoshravesh R Karki S Tapia R Balahadia CP Bandyopadhyay A Quick WP Furbank R Sage TL Langdale JA 2017 Re-creation of a key step in the evolutionary switch from C3 to C4 Leaf Anatomy Current Biology 27 3278 3287 10.1016/j.cub.2017.09.040 29056456
Wang M Li W Fang C Xu F Liu Y Wang Z Yang R Zhang M Liu S Lu S Lin T Tang J Wang Y Wang H Lin H Zhu B Chen M Kong F Liu B Zeng D Jackson SA Chu C Tian Z 2018 Parallel selection on a dormancy gene during domestication of crops from multiple families Nature Genetics 50 1435 1441 10.1038/s41588-018-0229-2 30250128
Waters MT Wang P Korkaric M Capper RG Saunders NJ Langdale JA 2009 GLK transcription factors coordinate expression of the photosynthetic apparatus in Arabidopsis The Plant Cell 21 1109 1128 10.1105/tpc.108.065250 19376934
West MAL Yee KM Danao J Zimmerman JL Fischer RL Goldberg RB Harada JJ 1994 Leafy cotyledon1 is an essential regulator of late embryogenesis and cotyledon identity in arabidopsis The Plant Cell 6 1731 1745 10.1105/tpc.6.12.1731 12244233
Wietrzynski W Engel BD 2021 Chlorophyll biogenesis sees the light Nature Plants 7 380 381 10.1038/s41477-021-00900-6 33875836
Xu X Zhiguo E Zhang D Yun Q Zhou Y Niu B Chen C 2021 OsYUC11 -mediated auxin biosynthesis is essential for endosperm development of rice Plant Physiology 185 934 950 10.1093/plphys/kiaa057 33793908
Yakovlev MS Zhukova GY 1980 Chhlorophyll in embryos of angiosperm seeds, a review Botaniska Notiser 133 323 336
Yamamoto A Kagaya Y Toyoshima R Kagaya M Takeda S Hattori T 2009 Arabidopsis NF-YB subunits LEC1 and LEC1-LIKE activate transcription by interacting with seed-specific ABRE-binding factors The Plant Journal 58 843 856 10.1111/j.1365-313X.2009.03817.x 19207209
Yeh S-Y Lin H-H Chang Y-M Chang Y-L Chang C-K Huang Y-C Ho Y-W Lin C-Y Zheng J-Z Jane W-N Ng C-Y Lu M-Y Lai I-L To K-Y Li W-H Ku MSB 2022 Maize Golden2-like transcription factors boost rice chloroplast development, photosynthesis, and grain yield Plant Physiology 188 442 459 10.1093/plphys/kiab511 34747472
Zhang Y Su J Duan S Ao Y Dai J Liu J Wang P Li Y Liu B Feng D Wang J Wang H 2011 A highly efficient rice green tissue protoplast system for transient gene expression and studying light/chloroplast-related processes Plant Methods 7 30 10.1186/1746-4811-7-30 21961694
Zhang JJ Xue HW 2013 OsLEC1/OsHAP3E participates in the determination of meristem identity in both vegetative and reproductive developments of rice Journal of Integrative Plant Biology 55 232 249 10.1111/jipb.12025 23230849
Zhang D Tan W Yang F Han Q Deng X Guo H Liu B Yin Y Lin H 2021a A BIN2-GLK1 signaling module integrates brassinosteroid and light signaling to repress chloroplast development in the dark Developmental Cell 56 310 324 10.1016/j.devcel.2020.12.001 33357403
Zhang C Zhang J Tang Y Liu K Liu Y Tang J Zhang T Yu H 2021b Deep green panicle1 suppresses GOLDEN2-LIKE activity to reduce chlorophyll synthesis in rice glumes Plant Physiology 185 469 477 10.1093/plphys/kiaa038 33721900
Zhao L Xie L Zhang Q Ouyang W Deng L Guan P Ma M Li Y Zhang Y Xiao Q Zhang J Li H Wang S Man J Cao Z Zhang Q Zhang Q Li G Li X 2020 Integrative analysis of reference epigenomes in 20 rice varieties Nature Communications 11 2658 10.1038/s41467-020-16457-5 32461553
Zheng S Dong J Lu J Li J Jiang D Yu H Ye S Bu W Liu Z Zhou H Ding Y Zhuang C 2022 A cytosolic pentatricopeptide repeat protein is essential for tapetal plastid development by regulating OsGLK1 transcript levels in rice The New Phytologist 234 1678 1695 10.1111/nph.18105 35306663
Zhiguo E Li T Zhang H Liu Z Deng H Sharma S Wei X Wang L Niu B Chen C 2018 A group of nuclear factor Y transcription factors are sub-functionalized during endosperm development in monocots Journal of Experimental Botany 69 2495 2510 10.1093/jxb/ery087 29514259
