AtNF-YB1 induces early flowering in rice
The gene, AtNF-YB1, was amplified by PCR method (The primers were listed in Table S1) using cDNA from Arabiopspsis (ecotype Columbia) leaves and cloned into plant expression vector pCXUN (Chen et al. 2009). AtNF-YB1 was first transferred into rice Kasalath. The plants in T1 populations that contain the transgene detected in PCR showed early flowering; but the plants flowering lately as wild-type Kasalath were negative in the PCR examination. Further investigations were conducted using selected homozygous lines in the NLD condition (Chengdu, in summer) or NSD condition (Sanya, in winter). The transgenic Kasalath had fully flowered spikes with seed set; however, the wild-type was at its initial flowering stage (Fig. 1A). As a weak photoperiod-sensitive Indica landrace, Kasalath was flowering at about 105 days and 101 days in NLD and NSD, respectively. In contrast, the selected transgenic lines flowered early, at about 92 and 89 days in NLD and NSD, respectively.
To further explore the usage that AtNF-YB1 induces flowering early in breeding hybrid rice with shorter growth duration, the gene was transferred into rice Jinfeng B, a male sterile maintainer line of three-line hybrid rice. The gene was further introduced into Jinfeng A, the male sterile line, through backcrossing. As long as both the maintainer line and its sterile line containing AtNF-YB1as homozygote were ready, the two restorer lines, Chenghui 727 and Chuanhui 907, were respectively crossed with the sterile line to get two hybrids with the transgene, called transgenic Jinfeng X Chenghui 727 (T7) and transgenic Jinfeng X Chuanhui 907(T9).
The selected stable transgenic maintainer lines, Jinfeng B, and its sterile lines, Jinfeng A, were flowering early in the fields, about 12 days earlier under NLD conditions and four days under NSD. Both Jinfeng X Chenghui 727 (NT727) and Jinfeng X Chuanhui 907(NT907) were about 112 days to flowering under NLD. Under NSD, NT727 was 107 days and NT907 was 99 days to flowering. The transgenic hybrid T7 and T9 were both flowering at about 96 days after sowing under NLD in Chengdu. The transgenic plants were fully flowering (Figs. 1B and 1C), whereas the wild-type plants were at the flowering initial stages. Under NSD, the hybrids T7 and T9 were flowering at 95 and 88 days, respectively, in Sanya. All these indicate that the transgenic rice, either pure lines or hybrids, flower early under both NLD and NSD conditions.
AtNF-YB1 regulates shoot growth reducing plant height in rice
Although the height of the transgenic Kasalath plants seemed similar to that of the wild-type plants (Fig. 1A), height measurement at the flowering stage under NLD revealed a shorter phenotype of the transgenic plants (Fig. 1D). Despite few variations, the shorter statures of the transgenic lines were due to shorter length of each internode including the panicle (Fig. 1D). Shorter plant height among the transgenic lines was also observed during harvest under NLD (Table 1) and NSD (Table S2). The leaf numbers on the main culm were further investigated under NLD (Fig. 1E). Wild-type Kasalath had 12 or 13 leaves on the main culm with an average of 12.4, while the transgenic TK-01 had 11 or 12 leaves with an average of 11.3. These results indicate that AtNF-YB1 in Kasalath may result in reduced plant height with one less internode and shorter internode length of each internode, as well.
For hybrid rice, the transgenic plants of T7 and T9 were also shorter than its correspondent wild-type plants under NLD (Figs. 1F and 1G). However, the main panicle length was similar between transgenic and the correspondent wild-type plants for both hybrids. The first (I) and second (II) internode in the hybrids were even longer. Notably, the lengths of all the following internodes continually decreased. The final (VII) internode was almost invisible in both T7 and T9. Leaf numbers on the main culm in both T7 and T9 had about one less leaf than their correspondent wild-type hybrids (Fig. 1E). Under NSD, both hybrids with the transgene were shorter at harvest (Table S6). Therefore, expression of AtNF-YB1 in hybrid rice reduces plant height through regulation of the length of each internode; an internode is closer to the top, the length decrease is lesser. The top internode with panicle in the transgenic plants is even comparable to that in their wild-type plants.
AtNF-YB1 affects grain yield in rice
The transgenic Kasalath was cultivated in natural fields under NLD twice to investigate their agronomic traits (Table 1). Except for plant height, the main panicle length, the number of grains per panicle, and 1000-grain weight were reduced in the transgenic plants. However, the number of productive panicles per plant increased, decreased, or varied in different years or lines. The same phenomena were also observed under the NSD condition (Table S2). Thus, the transgenic Kasalath with AtNF-YB1 showed lower grain yields in both NLD and NSD conditions.
Several agronomic traits were investigated under NLD and NSD conditions for the two hybrids at harvest (Table 2). Plant heights were reduced; the main panicle lengths of transgenic and wild-type lines were similar in the mature plants under NLD conditions. The number of productive panicles was increased, but the 100-grain weight was reduced in the transgenic hybrids. The grain-setting rate between the transgenic and correspondent non-transgenic lines was similar, and the number of filled grains per panicle decreased in T7 but increased in T9. Thus, higher grain yields were recorded for the hybrids.
Under NSD conditions, the hybrids were tested twice in the winter season in Sanya. The reduced plant height with a shorter panicle was observed in both transgenic hybrids so that the filled grains on the panicle were reduced (Table S3). The 100-grain weight and the grain-setting rates decreased. However, the number of productive panicles per plant remained the same or increased. All of these resulted in a reduction of grain yields of both T7 and T9 in NSD.
AtNF-YB1 stimulates flowering pathways in rice
To explore the molecular basis of the early flowering, transcription levels of two florigens, Hd3a and RFT1, were measured first by a quantitative real-time PCR in Kasalath grown under NLD conditions. Both Hd3a (Fig. 2A) and RFT1 (Fig. 2B) were up-regulated in the transgenic line at the eighth week after sowing with peaks in the early morning or at midnight (Figs. 3A and 3B). The expressions of Hd1 (Figs. 2C and 3C) and Ehd1 (Figs. 2D and 3D), two major regulators of the florigen gene transcriptions (Sun et al. 2014), were induced in the sixth and seventh week before Hd3a and RFT1, respectively. Further analysis of the upstream genes in the flowering pathway revealed that Ghd7 (Figs. 2E and 3E) was clearly down-regulated before up-regulation of the florigen genes, but DTH8 (Figs. 2F and 3F) and PRR37 (Figs. 2G and 3G) had some up-or down-regulation after the flowering transition.
A diurnal expression experiment was then carried out under NSD conditions. Both Hd3a and RFT1 were up-regulated in the transgenic Kasalath (Figs. S1A and S1B). Hd1 and Ehd1 were also up-regulated (S1C and S1D). The expression of Ghd7 was decreased (Fig. S1E), but DTH8 and PRR37 were up- or down-regulated (Figs. S1F and S1G). These results are consistent with the observations under the NLD condition.
Hd1 in Kasalath genome is a non-functional allele (Yano et al. 2000). The observed up-regulation of Hd1transcriptions was thus not able to regulate the flowering induction. Ghd7, as only one among three upstream genes of Ghd7, PRR37, and DTH8 in the pathway of Ehd1-Hd3a/RFT1, was markedly down-regulated. Therefore, the suppression of the Ghd7 expression played a major role in promoting flowering in both NLD and NSD.
The genes in flowering-time pathways were also investigated in the two hybrids under NLD and NSD conditions. Under NLD, both T7 and T9 showed elevated expressions of both Hd3a (Figs. 2A and 3A) and RFT1 (Figs. 2B and 3B) at approximately the ninth week after sowing. The expressions of Hd1 (Fig.2C and 3C) and Ehd1 (Fig.2D and 3D) were also clearly up-regulated before the induction of the two florigen genes. Among the three upstream genes in the Ehd1- H3a/RFT1 pathway, Ghd7 (Figs. 2E and 3E) was the only dramatically down-regulated gene in the two transgenic hybrids before the up-regulation of the florigen genes. The transcript levels of DTH8 (Figs. 2F and 3F) and PRR37 (Figs. 2G and 3G) were only a little increased or decreased at the stage of flowering induction.
Under NSD conditions, expressions of Hd3a and RFT1 increased in both T7 and T9 (Fig. S1A, S1B). Ehd1 and Hd1 were also up-regulated (Fig. S1C, S1D). Ghd7 was dramatically down-regulated again (Fig. S1E). But, the transcript levels of DTH8 were unchanged (Fig. S1F), while PRR37 showed a little down-regulation (Fig. S1G). Thus, in the hybrids, suppressing the transcription of Ghd7 to activate the Ghd7-Ehd1-Hd3a/RFT1 pathway was also a major reason for the early flowering in both NLD and NSD conditions. However, Hd1 is a functional allele in Jinfeng B (GenBank, OP030521), so Hd1-Hd3a/RFT1 pathway could also impact the flowering induction in the hybrids.
AtNF-YB1 impacts multiple processes in rice
To further investigate the roles of AtNF-YB1 in rice, an RNA-seq experiment was carried out. Young spikes grow very fast after flowering induction in rice, playing a critical role in final grain yield. Therefore, spikes with a length of 0.5–1.0 cm were collected for the study. The experiment was designed to compare the transcriptions between the transgenic and corresponding wild-type plants, including Kasakath, the maintainer Jinfeng B, and two hybrids. The two hybrids were selected to be validated by quantitative real-time PCR with 16 genes (The primers were listed in Table S1). The PCR results were consistent with the results from RNA-Seq (Fig. S2).
Based on the criteria Fold Change (FC) ≥ 1.5 and P value < 0.01, 983 and 692 DEGs (Table S4) were identified in Kasalath (Table S5) and Jinfeng B (Table S6), respectively. In Kasalath, the up-regulated genes were in greater number than the down-regulated ones, whereas in Jinfeng B, the down-regulated genes were more than the up-regulated ones. KEGG (Kyoto Encyclopedia of Genes and Genomes) pathway enrichment analysis showed that the top five regulated pathways in Kasalath were protein processing in the endoplasmic reticulum, endocytosis, galactose metabolism, spliceosome, and plant hormone signal transduction (Fig. 4A). For Jinfeng B, the pathway of protein processing in endoplasmic reticulum was also the most significantly enriched pathway as in Kasalath. Interestingly, three pathways, photosynthesis-antenna proteins, photosynthesis, and carbon fixation in photosynthetic organisms, appeared among the top five pathways in Jinfeng B (Fig. 4B). Looking into the transcript levels of the significantly regulated genes in the pathway of photosynthesis - antenna proteins, six among the seven genes were also up-regulated in Kasalath (Table S9). In the photosynthesis pathway, seven out of the 12 genes were expressed higher in Kasalath (Table S10). In the pathway of carbon fixation in photosynthetic organisms, four out of nine genes were up-regulated in Kasalath (Table S11).
In the two hybrids, 471 and 308 DEGs (Table S5) were respectively identified in Jinfeng X Chenghui 727 (Table S7) and Jinfeng X Chuanhui 907 (Table S8). Further analysis of the data (Table S5) revealed that 101 DEGs were present in both hybrids. Up-regulated DEGs were up-regulated in both, and down-regulated DEGs were down-regulated in both with only one exception. The KEGG pathway enrichment analysis (Figs. 4C and 4D) revealed that the top five regulated pathways contained protein processing in endoplasmic reticulum, which was down-regulated in both Jinfeng X Chenghui 727 and Jinfeng X Chuanhui 907. Further, the pathways of photosynthesis-antenna proteins, photosynthesis, and carbon fixation in photosynthetic organisms were up-regulated in Jinfeng X Chuanhui 907 as observed in Jinfeng B. However, such pathways related to photosynthesis were not observed in Jinfeng X Chenghui 727. But, some of the genes, of course not all, showed also up-regulation (Table S12, S13, and S14) in the hybrid.
The data demonstrate that AtNF-YB1 in rice markedly alters the transcriptome. However, the actual actions may be different depending on the genetic background. The pathway of protein processing in endoplasmic reticulum might be one of the highly impacted pathways. Plant photosynthesis might be significantly affected, as well.
AtNF-YB1 impacts photosynthesis in rice
The discovery of alteration of the photosynthesis pathways further led to measuring photosynthetic traits. Both transgenic Jinfeng B and hybrid T9 showed markedly higher photosynthetic rates. However, transgenic Kasalath and T7 increased or decreased by 5%, respectively (Fig. 5A). These results are consistent with the transcript levels of the genes in the three pathways related to plant photosynthesis from the RNA-seq, indicating the influence of AtNF-YB1 on photosynthesis in rice may depend on the genetic background.
The chlorophyll contents in the leaves were then analyzed. Obtained data showed either chlorophyll a or b was a little increased in the transgenic Kasalath but decreased in the transgenic Jinfeng B (Fig. 5B). However, both hybrids showed a clear increase of both chlorophyll a and b in the transgenic plants (Fig. 5C). The data indicate that enhancement of the photosynthetic pigments in the transgenic hybrids, whereas, in pure line, the effects were quite small, and up-or down-regulation also depends on the genetic background.