Cytokinins and gibberellins (GAs) play antagonistic roles in regulating reproductive meristem activity. Cytokinins have positive effects on meristem activity and maintenance. During inflorescence meristem development, cytokinin biosynthesis is activated via a KNOX-mediated pathway. Increased cytokinin activity leads to higher grain number, whereas GAs negatively affect meristem activity. The GA biosynthesis genes GA20oxs are negatively regulated by KNOX proteins. KNOX proteins function as modulators, balancing cytokinin and GA activity in the meristem. However, little is known about the crosstalk among cytokinin and GA regulators together with KNOX proteins and how KNOX-mediated dynamic balancing of hormonal activity functions. Through map-based cloning of QTLs, we cloned a GA biosynthesis gene, Grain Number per Panicle1 (GNP1), which encodes rice GA20ox1. The grain number and yield of NIL-GNP1TQ were significantly higher than those of isogenic control (Lemont). Sequence variations in its promoter region increased the levels of GNP1 transcripts, which were enriched in the apical regions of inflorescence meristems in NIL-GNP1TQ. We propose that cytokinin activity increased due to a KNOX-mediated transcriptional feedback loop resulting from the higher GNP1 transcript levels, in turn leading to increased expression of the GA catabolism genes GA2oxs and reduced GA1 and GA3 accumulation. This rebalancing process increased cytokinin activity, thereby increasing grain number and grain yield in rice. These findings uncover important, novel roles of GAs in rice florescence meristem development and provide new insights into the crosstalk between cytokinin and GA underlying development process.
Rice panicle architecture, a valuable composite agronomic trait that includes grain number per panicle (GNP), panicle length and so on, is strongly associated with rice grain yield. GNP is one of the most important agronomic characteristics of ideal plant architecture [1]. To improve rice grain yields to meet the needs of the rapidly growing population, numerous studies have focused on identifying and cloning genes/QTLs contributing to rice panicle architecture development. Many genes and pathways have recently been identified, including transcriptional and plant hormone regulators that contribute to the reproductive meristem activity maintenance processes. Cytokinins play a fundamental role in regulating reproductive meristem activity by promoting cell division [2]. Grain number 1a (Gn1a), a cytokinin metabolism-related gene, encodes a cytokinin oxidase/dehydrogenase (OsCKX2) that catalyzes the degradation of active cytokinins in reproductive meristems. Thus, a null allele of Gn1a leads to improved rice grain yield through increased active cytokinin levels and reproductive meristem activity [3]. Another gene, LONELY GUY (LOG) encodes a cytokinin nucleoside 5’-monophosphate phosphoribohydrolase. LOG transcripts are specifically enriched in the apical regions of vegetative and reproductive meristems. LOG functions in the activation of cytokinin, catalyzing the conversion of inactive cytokinins to biologically active forms. Reduced active cytokinin levels in the meristem due to malfunctioning of cytokinin activation is likely responsible for the defective meristem activity in the log mutant [4]. In additions, the zinc finger transcription factor DROUGHT AND SALT TOLERANCE (DST) directly induces the expression of OsCKX2 in the inflorescence meristems. The mutant allele DSTreg1 reduces OsCKX2 expression, thus increasing cytokinin levels in the inflorescence meristem, and therefore, the number of panicle branches and grains [5,6]. Gibberellins (GAs) are crucial for plant growth and developmental processes, such as seed germination [7], grain setting [8] and so on. However, unlike cytokinins, GAs are primarily associated with high yield rice breeding due to their roles in plant height promotion. Most mutants or RNAi transgenic lines of GA biosynthesis genes, including CPS, KS, KAO [9], KO [10], GA20oxs [11–13] and GA3oxs [14], show dwarfism phenotypes, which results in improved lodging resistance, a valuable trait for rice breeding under high inputs [15]. At the same time, transgenic-activated expression of GA catabolism genes, GA2oxs, also leads to dwarfism [16,17]. However, GA signals are also active in inflorescence meristems. OsGA20ox2, OsGA3ox2, Gα and SLR1 are highly expressed in inflorescence meristems and leaf primordia [18]. In maize, the expression domains of GA2ox1 and KN1 (a maize KNOX gene) overlap, mainly at the base of the shoot apical meristem. The KNOX gene KN1 directly induces GA2ox1 expression in reproductive meristems [19]. In tobacco and Arabidopsis, GA20ox expression could be directly excluded from the corpus of the shoot apical meristem [20,21]. These findings suggest that GAs are detrimental to meristem activity. Although the importance of GAs in meristem establishment and maintenance has been recognized, the GA biosynthesis and regulatory networks underlying this process are largely unknown, and it also remains to be determined whether certain GA biosynthesis and regulatory genes can be useful for increasing grain number and yield in rice. KNOX proteins are a class of homeodomain transcription factors that function in meristem establishment and maintenance. OSH1 (a rice KNOX gene) can directly activate the expression of other KNOX paralogs (OSH15, for example) and itself. The positive autoregulation of KNOX genes and activation by cytokinin are both essential for meristem maintenance [22]. In rice and Arabidopsis, KNOX proteins can activate cytokinin biosynthesis in the meristems through the induction of genes encoding adenosine phosphate isopentenyltransferase (IPT). IPTs are important enzymes that convert ATP, ADP and AMP to the iP riboside 5’-triphosphate (iPRTP), iP riboside 5’-diphosphate (iPRDP) and iP riboside 5’- moophosphate (iPRMP) forms [23,24]. As KNOX proteins reduce GA activity, they play an indispensable role in maintaining shoot apical meristem activity, probably by balancing cytokinin and GA activity in the meristems, increasing cytokinin levels and reducing GA levels [25,26]. Here, we report the identification and characterization of a QTL, Grain Number per Panicle1 (GNP1), which encodes rice GA biosynthetic protein OsGA20ox1. We propose that the upregulation of GNP1 in the inflorescence meristems may increase cytokinin activity via a KNOX-mediated feedback regulation loop and increase GA catabolism activity through inducing the expression of GA2oxs. This process would result in increased cytokinin activity, rebalancing cytokinin and GA activity and increasing grain number and grain yield. These results provide insights into the mechanism underlying KNOX-mediated cytokinin and GA crosstalk during rice inflorescence meristem development, and they suggest that GNP1 is a suitable target gene for high yield rice breeding.
To identify QTLs, we constructed two sets of reciprocal introgression lines (ILs) derived from a japonica rice variety Lemont (LT) and an indica variety Teqing (TQ), TQ-ILs and LT-ILs. In these two ILs, multiple QTLs for Grain Number per Panicle (GNP) were identified in Beijing and Sanya, respectively (S1 Table). Among these, QTLs affecting GNP in the RM227–RM85 region on chromosome 3 were detected in both TQ- and LT-ILs, suggesting that this QTL is stable for the grain number trait in rice. This QTL was designated Grain Number per Panicle1 (GNP1). From 201 LT-ILs, an IL named GG306 (BC3F4), containing chromosome segment RM227–RM85 from TQ and 92. 6% of the genetic background of LT, was selected (Fig 1A) and backcrossed twice to LT. Self-pollination of BC5F1 plants heterozygous for this fragment resulted in heterozygous near-isogenic lines (NILs) with almost all of the genetic background of LT except for the introgressed segment (Fig 1A). The BC5F2 was successively self-pollinated several times to obtain segregating NIL-F2 (BC5F3, BC5F4 and BC5F5) populations for fine mapping of GNP1 and construction of NILs, NIL-GNP1LT and NIL-GNP1TQ (Fig 1B). An analysis of a BC5F3 population of 163 individuals derived by self-pollination of the BC5F2 heterozygotes at the region RM227–RM85 showed that the trait segregated as a single locus with a Mendelian ratio, which was confirmed by data from BC5F4 families (S1 Fig and S2 Table). Through map-based cloning of GNP1, we narrowed the GNP1 locus down to a 33. 7 kb region between SL65 and SL54 (Fig 1C and S2 Fig). This region contains four predicted genes (LOC_Os03g63970, LOC_Os03g63980, LOC_Os03g63990 and LOC_Os03g63999, http: //rice. plantbiology. msu. edu/cgi-bin/gbrowse). To further investigate the effects of the GNP1 locus on grain number and other traits, we analyzed near-isogenic lines, NIL-GNP1LT and NIL-GNP1TQ, in the LT genetic background, which only differed in the ~66. 1 kb region containing GNP1 derived from LT and TQ (Fig 1A). We observed a significant increase in the total grain number per panicle (GNP; +56%), filled grain number per panicle (FGN; +28%) and secondary branch number (SBN) in NIL-GNP1TQ (Fig 1D, Fig 1E and S3F Fig, the same pattern in SBN between LT and TQ (S3G Fig) ), but only a small increase in plant height (+8%; Fig 1B and S3A Fig), a slight decrease in grain length (-4%; S3B Fig), grain width (-5%; S3C Fig) and 1,000-grain weight (-12%; S3D Fig) and no effect on panicle length (S3E Fig) and primary branch number (S3F Fig the same pattern in PBN between LT and TQ (S3G Fig) ) compared with the NIL-GNP1LT isogenic control in plants grown in Shanghai. These results indicate that the GNP1TQ locus in NIL-GNP1TQ has pleiotropic effects on rice development, primarily on inflorescence development, especially secondary branch number and grain number. To determine whether GNP1TQ affects grain yield, we evaluated the grain yields of NIL-GNP1TQ and the isogenic control (Lemont), together with other related traits. In different fields, the grain number was still substantially higher in NIL-GNP1TQ than in the control, leading to a significant increase in grain yield (5. 7–9. 6%) despite the slightly reduced grain weight (Table 1 and S3 Table). These results suggest that the GNP1TQ locus can potentially be used in high yield rice breeding. According to the mapping results, LOC_Os03g63980 and LOC_Os03g63990 are predicted to encode transposon and retrotransposon proteins, LOC_Os03g63999 encodes a small peptide with unknown function and LOC_Os03g63970 encodes GA 20-oxidase 1, which is thought to catalyze the conversion of GA12 to GA20 within a multi-step process. Therefore, LOC_Os03g63970 is the most likely candidate for the GNP1 locus. We sequenced the promoter (2 kb before ATG) and LOC_Os03g63970 in both TQ and LT. The two parents exhibited base differences at 21 positions in the promoter region, including 17 single-base substitutions, as well as two single-base and two multi-base insertions and deletions. The coding region contains two single-base substitutions, one of which leads to an amino acid substitution (S4 Fig). These results suggest that the sequence differences in the promoter and coding region of this gene might lead to changes in gene expression levels and protein function and may help increase grain number in NIL-GNP1TQ. To validate this hypothesis, we obtained the LOC_Os03g63970 T-DNA gain-of-function mutant gnp1-D from the Rice T-DNA Insertion Sequence Database. TAIR-PCR screening showed that the T-DNA was inserted at position -514 to -492 of the LOC_Os03g63970 promoter relative to the start codon ATG (Fig 2A), which constitutively induces the expression of LOC_Os03g63970 throughout the plant. We analyzed traits of the homozygous gnp1-D mutant and control via PCR with specific primers designed based on the insertion sequence (Fig 2A and Fig 2B), finding a significant increase in plant height (Fig 2C and Fig 2D) with increasing LOC_Os03g63970 expression in flag leaves (Fig 2E). Interestingly, a substantial increase in GNP (+51. 5%) and FGN (+71. 6%) were also observed (Fig 2F and Fig 2G). These results suggest that LOC_Os03g63970 is the gene for GNP1 and that the increased GNP1 expression in this mutant might influence GA biosynthesis during rice panicle meristem development. We then constructed a binary vector harboring the GNP1TQ coding sequence (CDS) driven by a CaMV 35S promoter, which we used to transform japonica rice (O. sativa L.) variety Zhonghua 11 (ZH11), whose GNP1 CDS matches that of LT. GNP1 was expressed at levels several hundred- to over a thousand-fold that of CK (transgenic negative control) in flag leaves (Fig 3A). Compared with CK, the GNP of line p35S: : GNP1TQ-3 increased by 36. 3%, accompanied with hugely increased height (S5A and S5B Fig) and greatly increased sterility, while lines p35S: : GNP1TQ-1 and p35S: : GNP1TQ-2 had significantly increased GNP (FGN) by 27. 8% (35. 5%) and 26. 5% (33. 4%) (Fig 3B and Fig 3C), and slightly increased height (S5A and S5B Fig). These results indicate that the expression disturbances associated with the promoter activity variations at the GNP1 locus are responsible for the phenotypic variation in GNP and plant height with a dose-dependent manner and a very high expression level of GNP1 may have a negative effect on seed setting rate. Then, in order to find out whether decreased expression of GNP1 could show some negative effect on grain number phenotype, we transformed ZH11 with the mimic artificial microRNA oligo sequence designed for GNP1 silencing driven by the CaMV 35S promoter. Interestingly, the grain number of six transgenic-positive independent lines increased (S6A Fig), which was negatively correlated with GNP1 expression (S6B Fig). These lines also had reduced plant height (S6C and S6D Fig). These results indicate that the reduced expression of GNP1 might contribute to attenuated GA biosynthesis activity, leading to reduced GA levels and partially reducing the negative effects of GAs on maintaining inflorescence meristem activity [26], which might be responsible for the higher grain number in these mimic artificial miRNA transgenic lines. To further confirm the function of GNP1LT CDS, we transformed NIL-GNP1LT with GNP1LT CDS driven by the GNP1 promoter from Lemont (pGNP1LT). Similar to gnp1-D gain-of-function mutant and GNP1TQ overexpression lines, as the expression level of GNP1 increased (up to nearly ten-fold compared to the control; Fig 3D), we observed an increase in GNP and FGN (Fig 3E), as well as plant height (S5C Fig). These results indicate that both GNP1LT and GNP1TQ could affect panicle development. These results indicate that the accumulation of GNP1LT or GNP1TQ transcripts (or both) in the plant has a positive effect on grain number and plant height. To determine whether the differences between the GNP1LT and GNP1TQ promoter regions (S4 Fig) influence GNP1 expression, and account for the differences in grain number, we analyzed the expression patterns of GNP1 between NIL-GNP1LT and NIL-GNP1TQ in different tissues during panicle initiation to the booting stage. GNP1 was mainly expressed in developing panicles and nodes (S7 Fig), which is consistent with effects of this gene on grain number and plant height. In addition, compared to NIL-GNP1LT, GNP1 transcripts were much more abundant in NIL-GNP1TQ tissues (S7 Fig). Meanwhile, GNP1 expression in seedling leaf sheaths was negatively correlated with the dose of GA3 used for treatment (Fig 4A and Fig 4C) and positively correlated with that of the GA biosynthesis inhibitor uniconazole-P (Fig 4B and Fig 4D), suggesting that GNP1 expression is controlled by biologically active GA levels. The GNP1LT allele was much more sensitive to uniconazole-P treatment and endogenous GA signal feedback regulation (Fig 4B and Fig 4D), probably due to the sequence variations among promoters. We also investigated GNP1 expression in the shoot apical meristems and inflorescence meristems. Similar to OSH1, a key factor in rice meristem maintenance and regulation, GNP1 was also expressed in the apical regions of these meristems (S8 Fig). OSH1 expression signal in NIL-GNP1TQ meristems is still strong and specific (S8 Fig), These results suggest that during NIL-GNP1TQ inflorescence meristem development, the sequence variations of the promoter might lead to a failure to maintain low GNP1 expression level, resulting in induced GNP1 expression in the panicle meristems of NIL-GNP1TQ. The above findings demonstrate that the variations in promoters leading to changes in GNP1 expression in the panicle meristems are the main contributor to the differences in grain number between NIL-GNP1TQ and NIL-GNP1LT. Moreover, the total GNP was positively correlated with the expression level of GNP1. In vitro, GNP1 (GA20ox1) directly catalyzes the biosynthesis of GA53, GA44, GA19 and GA20 in the early-13-hydroxylation pathway with various catalyzing efficiency for each steps [27]. GA20 is then used for GA1 and GA3 biosynthesis via catalyzing by GA3oxs (Fig 5A) [28]. We therefore measured the contents of five endogenous GA biosynthesis intermediates, finding that GA20 and GA12 accumulated preferentially in the panicle meristems of NIL-GNP1TQ, whereas GA44 levels were much lower and there were no changes in GA19 levels relative to NIL-GNP1LT (Fig 5B and Fig 5C), indicating that GA20 biosynthesis was accelerated. GNP1 mRNA levels were much higher in NIL-GNP1TQ, suggesting that the catalytic activity of GNP1 markedly increased as well, leading to higher accumulation of the GA biosynthesis intermediate GA20. The increased accumulation of GA12 suggests that GA biosynthesis activities including GA12 biosynthesis and previous steps might have been activated in this line. However, in the panicle meristems of NIL-GNP1TQ, bioactive GA1 and GA3 were not detected although they were detected in NIL-GNP1LT (Fig 5D), indicating that GA1 and GA3 levels in the NIL-GNP1TQ panicle meristems were too low to quantify. Consistent with this result, the GA signal transduction-related genes RGL3 and SLR1 were induced in this line (S9 Fig). RGL3 and SLR1 are DELLA proteins and negative regulators of GA signaling, whose degradation by GAs in collaboration with GID1 (gibberellin receptor) [29,30] and F-box protein is a key event in GA signaling activation [31–33]. Indeed, bioactive GA1 and GA3 levels were reduced in NIL-GNP1TQ panicle meristems. By contrast, most GA biosynthesis-related genes were upregulated, including OsKAO, OsKO, OsKS, OsCPS and OsGA3ox2 (S9 Fig), leading to increased GA12 levels (Fig 5B), likely due to feedback activation by reduced bioactive GA (GA1 and GA3) levels. At the same time, most bioactive GA catabolism genes, i. e. , OsGA2oxs (Fig 5E), were induced. As GA2oxs directly catalyze progressive catabolic processes that convert active GAs into inactive forms (Fig 5A), the increased catabolic activities in NIL-GNP1TQ panicle meristems regulate GA levels much more effectively, regardless of the activated GA biosynthesis process described above. Based on these findings, during NIL-GNP1TQ panicle meristem development, GA (GA1 and GA3) levels happened to be reduced, although the catabolic activities of GNP1 were enhanced. Cytokinins significantly affect reproductive meristem activity [2]. The abnormal GA metabolism in NIL-GNP1TQ observed in the current study might be caused by KNOX-mediated responses. To investigate this possibility, we analyzed the expression of five rice KNOX genes, including OSH1, OSH6, OSH15, OSH43 and OSH71. The expression of these genes significantly increased in the panicle meristems of NIL-GNP1TQ (Fig 6A). OsIPTs, which are directly regulated by KNOX proteins, were also upregulated in NIL-GNP1TQ, as was the cytokinin activating gene LOG (Fig 6B), perhaps leading to cytokinin accumulation. We also examined endogenous cytokinins levels in NIL-GNP1TQ, finding that the levels of several cytokinins and cytokinin biosynthesis intermediates increased in this line (Fig 6C to 6F), leading to increased expression of cytokinin signal response factors (Fig 6G). These results indicate that cytokinin activity was substantially enhanced in NIL-GNP1TQ panicle meristems, resulting in increased grain number compared to NIL-GNP1LT.
A previous in vitro study showed that recombinant OsGA20ox1 could catalyze the conversion of GA12 and GA53 to GA9 and GA20, but it acts more effectively on GA53 [27]. The present study shows that GNP1 encodes a rice OsGA20ox1 protein. OsGA20ox1 activity is induced via increased expression of GNP1, which increases GA20 levels in vivo. Moreover, GNP1 transcript levels in seedling leaf sheaths were positively correlated with the treatment dose of uniconazole-P and negatively correlated with that of GA3 (Fig 4C and 4D), suggesting that GNP1 expression is controlled by biologically active GA levels. Moreover, NIL-GNP1LT was much more susceptible to endogenous GA signal feedback regulation than NIL-GNP1TQ, likely due to the sequence variations among promoters leading to altered expression of GNP1. GNP1 transcripts were mainly detected in newly initiated panicles and in apical regions of meristems overlapping with OSH1 (a rice KNOX gene) expression (S8 Fig). This specific expression pattern implies that GNP1 also plays a fundamental role in regulating panicle meristem activity that is similar to that of cytokinin biosynthesis and signaling genes. The increased grain number of NIL-GNP1TQ due to enhanced expression of GNP1 supports this notion. Cytokinins positively regulate reproductive meristem activity [2], GAs are detrimental to meristem activity [20,21] and KNOX proteins play an irreplaceable role in balancing cytokinin and GA activity in the meristem [25]. We observed increased cytokinin activity in the panicle meristems of NIL-GNP1TQ, including KNOX-mediated induction of OsIPTs and increased levels of cytokinins and cytokinin biosynthesis intermediates, together with enhanced cytokinin responses. In additions, these plants failed to accumulate bioactive GA1 and GA3 and exhibited significantly increased KNOX transcript levels. Taken together, these results demonstrate that increased GNP1 activity positively induces the expression of KNOX genes via a feedback loop (Fig 7, red arrow). This promotion of KNOX gene expression leads to increased cytokinin activity through directly inducing OsIPT expression, as well as upregulation of GA2oxs, which negatively regulate GA biosynthesis, thereby reducing GA1 and GA3 levels. The activation of GA biosynthesis might be due to feedback regulation compensating for the defects in GA1 and GA3 accumulation, leading to increased accumulation of GA12. The tendency for activated GA biosynthesis may be much less effective than that for GA catabolism. This feedback mechanism rebalances cytokinin and GA activity, resulting in increased cytokinin levels and contributing to the higher GNP1 expression level of NIL-GNP1TQ. On the other hand, decreased expression of GNP1 could lead to lower GA1 and GA3 level in those positive GNP1 mimic artificial miRNA transgenic lines, which might eliminate the suppression effect of higher GA1 and GA3 level on meristem activities, and increase grain number in turn (Fig 7). We propose that during inflorescence meristem development and maintenance processes, increased expression of GNP1 in those NILs leads to promoted cytokinin activities and gives increased grain number and yield, while decreased expression of GNP1 in those mimic artificial miRNA transgenic lines most probably contributes to alleviation of the detrimental effect of gibberellins to meristem activity, according to those previous reports, which in turn also gives increased grain number. Numerous efforts aimed at increasing food production to sustain the growing population have focused on elucidating the mechanisms underlying the development of several important agronomic traits in rice, such as panicle architecture. In this study, we cloned a rice GA20ox1 gene, GNP1, whose expression strongly increases rice grain number. Increasing GNP1 expression may be useful for high yield rice breeding, as these GNP1 higher-expressed NILs exhibited increased grain number and grain yield, although they were also slightly taller than the controls. When we overexpressed GNP1 in ZH11, similar results were obtained, thus representing a new strategy for high yield rice breeding.
Two sets of reciprocal introgression lines (ILs) derived from a japonica rice (O. sativa L.) variety Lemont and an indica variety Teqing were used as materials for QTL mapping [34]. ZH11 and Lemont were used for the transgenic experiments. The gnp1-D T-DNA mutant line PFG_2D-41474. R was identified from the Rice Functional Genomic Express Database (RiceGE, http: //signal. salk. edu/cgi-bin/RiceGE) and obtained from the Rice T-DNA Insertion Sequence Database (RISD DB, http: //cbi. khu. ac. kr/RISD_DB. html) [35]. Oligo sequences used for genotyping the progeny of gnp1-D T-DNA insertional line are shown in S4 Table. For map-based cloning of GNP1, we performed genotyping of 5,500 BC5F3 individuals from five BC5F2 plants that were heterozygous only at the region RM227–RM85, harboring five markers. We identified 16 informative recombinants of four genotypes within this region. Using multiple comparisons of the homozygous recombinant BC5F4 lines for GNP with the non-recombinant controls, we localized GNP1 to a 309. 5kb region between SL13 and RM85. Further fine mapping using 9,500 BC5F4 plants with six new markers between SL13 and RM85 identified six informative recombinants and four genotypic classes in the target region. We localized GNP1 to a high-resolution linkage map by progeny testing of BC5F5 homozygous recombinant plants and narrowed the GNP1 locus down to a 33. 7 kb region between SL65 and SL54. Primers used for fine mapping are shown in S5 Table. GA3 and uniconazole-P treatment were carried out as previously described [36] with minor modifications. For GA3 treatment, manually dehulled seeds were sterilized with 75% ethanol for 1 min, washed three times with distilled water, sterilized with 2. 5% sodium hypochlorite for 35 min, washed five times with sterile distilled water and incubated on 1/2 MS medium at 4°C for 3 days in the dark. The germinated seeds were transferred to plastic containers containing 1% (w/v) agar with various concentrations of GA3 (63492-1G, Sigma-Aldrich). For uniconazole-P treatment, the seeds were incubated in distilled water with various concentrations of uniconazole-P (19701-25MG, Sigma-Aldrich) at 4°C for 24 h, followed by 26°C for an additional 24 h. The seeds were washed three times with distilled water and incubated for an additional 24 h in distilled water at 26°C. The germinated seeds were grown in 1% (w/v) agar in plastic containers. Seedlings were grown for 7 days under fluorescent light with a 12 h light/12 h dark photoperiod at 26°C. The second leaf sheath lengths of 48 seedlings per treatment were measured and analyzed. For qRT-PCR analysis, second leaf sheaths were also used, with six pooled replicates for each treatment. To produce the overexpression constructs, the full-length coding sequence of GNP1 was amplified from NIL-GNP1TQ and cloned into plant binary vector pCAMBIA1300 under the control of single CaMV 35S promoter. The artificial microRNA oligo sequences used for GNP1 silencing were designed as previously described [37] (http: //wmd3. weigelworld. org/cgi-bin/webapp. cgi? page=Home; project=stdwmd) and amplified using primer set G-11491 and G-11494. The oligo sequences were inserted into the XbaI and KpnI sites of pCAMBIA1300 containing one CaMV 35S promoter. Oligo sequences for three different target sites were independently used for construction and transformation. The overexpression and silencing plasmids were introduced into Agrobacterium tumefaciens strain EHA105 and transferred into the japonica variety ZH11. To produce the construct for the complementary test, 2. 2 kb promoter sequence with full-length coding sequences of GNP1 were amplified from NIL-GNP1LT. The sequences were then cloned into pCAMBIA1300, introduced into Agrobacterium tumefaciens strain EHA105 and used for transformation of NIL-GNP1LT. All constructs were confirmed by sequencing. The primer sets are shown in S6 Table, and plant transformation processes were carried out as previously described [38]. Total RNA was extracted from various plant tissues using TRIZOL Reagent (Invitrogen). Approximately 500 ng of total RNA was transcribed into first-strand cDNA using ReverTra Ace qPCR RT Master Mix with gDNA Remover (TOYOBO). Real-time PCR data were obtained using an ABI 7300 Real Time PCR System with Fast Start Universal SYBR Green Master Mix with ROX (Roche) and analyzed using the ΔΔCt method. The cycling parameters were 10 min at 95°C, followed by 40 cycles of amplification (95°C for 10 s and 60°C for 1 min). The ubiquitin and actin genes were used for normalization. The standard amplification slope for real-time PCR primer OsGA2ox1f/OsGA2ox1r was -3. 498971, which was used to calculate amplification efficiency. All analyses were repeated at least three times. Primer sets are shown in S7 Table. NIL-GNP1TQ and NIL-GNP1LT plants were grown in open fields for approximately 5 weeks. Freshly initiated panicles approximately 1 cm long were harvested, and ~1 g samples were used for measurements, with three independent biological repeats per sample. Quantification of endogenous GAs [39] and cytokinins [40] was performed as previously described. NIL-GNP1TQ plants were grown in open fields for approximately 3 weeks. Samples ~0. 5 cm in length including the meristem region were harvested and fixed in 4% (w/v) paraformaldehyde with 0. 1% Tween-20,0. 1% Triton-x-100 and 1% (v/v) 25% glutaraldehyde solution in 0. 1 M sodium phosphate buffer (pH 7. 4) overnight at 4°C. The samples were then dehydrated with a graded ethanol series followed by a dimethylbenzene series. The samples were then embedded in Paraplast Plus (Sigma, P3683), cut into 10 μm sections and mounted on pre-coated poly-prep slides (Sigma, P0425). Digoxigenin-labeled RNA probes were prepared following the instructions of the DIG RNA labeling kit (SP6/T7) (Roche, 11175025910). Hybridization and signal detection were performed as previously described [41]. The primer sets are shown in S8 Table. Yield and related traits for NIL-GNP1TQ and the isogenic control (Lemont) were evaluated at five locations: Beijing (40. 2°N, 116. 2°E); Nanning (22. 1°N, 107. 5°E), Guangxi province; Jingzhou (30. 3°N, 112. 2°E), Hubei province; Pingxiang (27. 6°N, 113. 9°E), Jianxi province and Sanya (18. 3°N, 109. 3°E), Hainan province, China. NIL-GNP1TQ and Lemont plants were grown in a randomized plot design with three replications per line. The area of each plot was 13. 2 m2, with a single plant transplanted per hill at 25 d after sowing and a spacing of 17 cm between hills and 25 cm between rows. As a basal dressing, 50 kg ha-1 each of N, P and K was applied the day before transplanting, and 30 kg ha-1 of N was applied twice as topdressing at 1 and 5 weeks after transplanting. At the heading stage, heading date (HD) and plant height (PH) were recorded when 30% of plants contained panicles in each line. At maturity, whole plots were harvested for yield measurements based on a 14% moisture content after air drying. Eight plants were sampled and dried in an oven at 70°C for 5 d for trait investigation, including panicle number per plant (PNP), panicle length (PL), filled grains per panicle (FGP), grain number per panicle (GNP), thousand grain weight (TGW), grain length (GL) and grain width (GW). QTLs affecting GNP were identified using IciMapping 3. 0 [42], combined with genotypic data for 157 SSRs and three morphological markers (Ph, gl-1 and C) for the ILs [34]. The permutation method was used to obtain empirical thresholds for claiming QTLs based on 1,000 runs in which the trait values were randomly shuffled [43].
Grain number per panicle, a valuable agronomic trait for rice yield improvement, is profoundly affected by reproductive meristem activity. This activity, in turn, is controlled by transcriptional and plant hormone regulators, especially KNOX proteins and cytokinins. However, little is known about the roles of GAs in these processes in rice and how the regulatory network functions due to the complexity of crosstalk between plant hormone regulators. In this study, we identify a novel GA biosynthesis gene in rice and demonstrate its role in improving grain number and grain yield. We also propose that the KNOX-mediated cytokinin-GA activity rebalancing mechanisms regulate inflorescence meristem development and maintenance processes, providing a possible tool for high-yield rice breeding.
plant anatomy
panicles
floral meristem
plant physiology
hormones
cereal crops
plant science
rice
model organisms
plant hormones
crops
inflorescences
plants
research and analysis methods
cytokinins
grasses
crop science
gibberellins
leaves
agriculture
biochemistry
plant biochemistry
plant and algal models
biology and life sciences
biosynthesis
meristems
organisms
The QTL GNP1 Encodes GA20ox1, Which Increases Grain Number and Yield by Increasing Cytokinin Activity in Rice Panicle Meristems