Mutations on ent-kaurene oxidase 1 encoding gene attenuate its enzyme activity of catalyzing the reaction from ent-kaurene to ent-kaurenoic acid and lead to delayed germination in rice


Authors: Hui Zhang aff001;  Ming Li aff001;  Dongli He aff001;  Kun Wang aff003;  Pingfang Yang aff001
Authors place of work: State Key Laboratory of Biocatalysis and Enzyme Engineering, School of Life Sciences, Hubei University, Wuhan, China aff001;  National Key Laboratory of Crop Genetic Improvement and National Centre of Plant Gene Research (Wuhan), Huazhong Agricultural University, Wuhan, China aff002;  School of Life Sciences, Wuhan University, Wuhan, China aff003
Published in the journal: Mutations on ent-kaurene oxidase 1 encoding gene attenuate its enzyme activity of catalyzing the reaction from ent-kaurene to ent-kaurenoic acid and lead to delayed germination in rice. PLoS Genet 16(1): e1008562. doi:10.1371/journal.pgen.1008562
Category: Research Article
doi: 10.1371/journal.pgen.1008562

Summary

Rice seed germination is a critical step that determines its entire life circle, with seeds failing to germinate or pre-harvest sprouting both reduce grain yield. Nevertheless, the mechanisms underlying this complex biological event remain unclear. Previously, gibberellin has been shown to promote seed germination. In this study, a delayed seed germination rice mutant was obtained through screening of the EMS induced mutants. Besides of delayed germination, it also shows semi-dwarfism phenotype, which could be recovered by exogenous GA. Through re-sequencing on the mutant, wild-type and their F2 populations, we identified two continuous mutated sites on ent-kaurene oxidase 1 (OsKO1) gene, which result in the conversion from Thr to Met in the cytochrome P450 domain. Genetic complementary analysis and enzyme assay verified that the mutations in OsKO1 gene block the biosynthesis of GA and result in the defect phenotypes. Further analyses proved that OsKO1 could catalyze the reaction from ent-kaurene into ent-kaurenoic acid in GA biosynthesis mainly at seed germination and seedling stages, and the mutations decrease its activity to catalyze the step from ent-kaurenol to ent-kaurenoic acid in this reaction. Transcriptomic and proteomic data indicate that the defect on GA biosynthesis decreases its ability to mobilize starch and attenuate ABA signaling, therefore delay the germination process. The results provide some new insights into both GA biosynthesis and seed germination regulatory pathway in rice.

Keywords:

Gene expression – Rice – Phenotypes – Embryos – Biosynthesis – Seeds – Seed germination – Gibberellins

Highlights

  • Strong evidences at both genetic and biochemical levels proved the function of OsKO1 gene in GA biosynthesis;

  • OsKO1 gene might function mainly at seed germination and young seedling stages;

  • Different steps of the reaction from ent-kaurene to ent-kaurenol, and then to ent-kaurenoic acid in GA biosynthetic pathway might be catalyzed by different motifs of ent-kaurene oxidase in the intact enzyme context.

Introduction

Rice is one of the main cereal grains that feed the world population. Rice seed with its independent role as a new generation could provide itself with enough reserves until the establishment of new seedling [1]. Defect of seed germination at either rate or timing could impair grain yield in agriculture [2]. It is vital to disclose the molecular mechanisms that regulate rice seed germination, which might help to facilitate molecular design in rice breeding.

Seed germination, which starts with water absorption by the dry seed and ends at radicle protrusion, is controlled by the crosstalk among various phytohormones and environmental stimuli. The phytohormones include abscisic acid (ABA), gibberellic acid (GA), ethylene, brassinosteroids, salicylic acid, cytokinin, auxin, jasmonic acid and oxylipins, with the first two ABA and GA playing the major roles. ABA controls storage reserves accumulation, and induces the expression of late embryogenesis abundant proteins (LEAs) that accumulated during seed maturation to protect seeds against desiccation. ABA also inhibits the germination of seeds on mother plant [3]. Considerable evidences have proven that ABA functions in dormancy maintenance [46], in which pre-harvest sprouting mutants in rice (Oryza sativa L.) all show reduced ABA content [7]. Based on the studies in Arabidopsis, it is known that series transcription factors form a cascade in ABA signaling pathway to control seed germination in which ABI1 and ABI2 are two homologous phosphatases that negatively regulate ABI3. ABI3 is a transcription activator that will be degraded by the recognition of AIP2, an E3 ligase ABI3 interaction protein [810]. WRKY41 directly activates ABI3 together with ABA to maintain dormancy [11]. ABI5 is a bZIP transcription factor that functions at the downstream of ABI3 to arrest seed germination [12]. ABI4 contained an APETALA2 (AP2) DNA binding domain in a separate pathway to restrain germination [13].

In contrast to ABA, GA functions in promoting germination and radical protrusion. GA induces the expression of enzymes that mobilize storage reserves like starch, proteins, and lipids [14]. Rice delayed germination mutants show reduced expression of either GA synthetic genes or GA responsive genes or both [8,15]. In rice, GA functions through activating the expression of α-amylase in aleurone layer to degrade starch into soluble sugars in endosperm, and the soluble sugars in endosperm are transported into embryo [16]. Reduced germination rate and seedling growth are always accompanied by decreased expression of α-amylase [17,18]. Mutants with reduced germination rate and delayed germination phenotypes all show the repression of α-amylase genes [8,19]. Recently characterized rice germination delayed mutants include germination defective 1 (GD1), knock-down mutant of rice lectin receptor-like kinase (OslecRK), and Oryza sativa delayed seed germination 1 (OsDSG1). GD1 is B3 domain containing repressor, and OsDSG1 is a ring finger E3 ligase. Their delayed germination phenotype accompanied with reduced plant height, decreased endogenous GA levels, reduced expression of α-amylase genes, down-regulated GA biosynthetic genes and up-regulated GA catabolic genes [8,15,19]. Especially an endosperm-imposed dormancy prolonged mutant seed dormancy1-2 (qSD1-2) has mutation on OsGA20ox2 was observed with semi-dwarf plant, decreased GA levels, extended dehydration and seed maturity [20]. Taken together, the control of seed germination relies on the fine-tuning of GA and ABA levels.

Seven genes have been shown to be involved in the biosynthesis of active GAs such as GA1, GA3, and GA4. Four single genes include ent-copalyl diphosphate synthase (CPS), ent-kaurene synthase (KS), ent-kaurene oxidase (KO, CYP701A), ent-kaurenoic acid oxidase (KAO, CYP88A) catalyze early steps of GA biosynthesis from precursor geranylgeranyl diphosphate (GGDP) to GA12. Two terpene synthases CPS and KS catalyze the conversion of GGDP to ent-kaurene on plastid [21,22] followed by two cytochrome P450 monooxygenases KO and KAO. KO converts ent-kaurene to ent-kaurenoic acid, and KAO catalyzes the formation of GA12 from ent-kaurenoic acid [23,24]. GA12 is the common precursor of active GA1 and GA4, while GA 13-oxidases are needed in the converting of GA12 into GA53 from 13-hydroxylation pathway. These are followed by the conversion of GA53 into GA1 in 13-hydroxylation pathway and the conversion of GA12 to GA4 in the non-13-hydroxylation pathway catalyzed by GA 20-oxidase (GA20ox) and GA 3-oxidase (GA3ox) that with small gene families [25]. Various mutants of the seven genes that show dwarfism phenotype have been isolated. There are five OsKO genes which are contiguously arranged in a 120 kb highly linked region on the sixth chromosome (S1A Fig). Previous study has shown that OsKO2 is the right gene encoding the correct ent-kaurene oxidase in GA biosynthetic pathway in rice [26,27]. Their study showed that another OsKO (It is named as OsKO1 in this study, but updated as OsKO5 in the newly released rice genome annotation version; S1A Fig) gene might also involve in GA biosynthesis [26]. OsKO2 is a Cytochrome P450 monooxygenase that may locate on the outer membrane of plastid [22]. The kinesin-like protein GDD1 binds to the promoter of OsKO2 and activates its expression to maintain endogenous GA levels in rice [28]. Recently, a study showed that both OsKO1 and OsKO5 might also involve in the GA biosynthesis in cooperation with OsKO2 based on the phenotyping of the corresponding gene’s CRISPR/Cas9 edited rice lines, of which both OsKO1 and OsKO5 edited lines were dwarfism [29]. However, no further direct evidences verified their involvement in GA biosynthesis.

To further explore the regulatory mechanisms underlying rice seed germination, we screened a rice mutant showing both delayed germination and semi-dwarfism phenotype. MutMap sequencing of two F2 progeny bulks determined two SNP changes on the eighth exon of OsKO1 gene resulting in a Thr to Met transition. Further biochemical and molecular genetic analyses proved that OsKO1 could function as well as OsKO2 in the catalyzing of ent-kaurene into ent-kaurenoic acid in GA biosynthesis during germination and early seedling growth in rice. The defect on this gene leads to the decreasing of GA content, and then disorder the balance between GA and ABA signaling, which delays seed germination in rice. Our results may help to obtain further understanding about both GA biosynthesis and its regulation on seed germination in rice.

Results

Phenotypic characterization of the germination delayed dwarfism mutant

Through screening the EMS mutagenesis mutants with Nipponbare background, a rice seed germination delayed mutant was obtained. To further characterize the phenotype, we reproduced it for four years, and conducted germination assay along with the Nipponbare wild-type to see its stable inheritance. The mutant seeds germinated much slower than wild-type with a 45% decrease in germination rate (Fig 1A). Obvious difference on the germination appeared at 36 h after imbibition (Fig 1A). Besides of the delayed germination phenotype, mutant also displayed smaller seeds, dwarfism, and dark green leaves (Fig 1A and 1C). Further studies through microscopy analysis showed that the dwarfism of mutant is resulted from its reduced cell size in internode (n ≥ 20, p = 0.002) and leaf sheath (n ≥ 20, p = 3.7 × 10−5) (Fig 1C). All these are typical GA related phenotype, which indicates that it may be a GA deficient or insensitive mutant based on previous studies [3033].

Phenotypic characterization of the mutant.
Fig. 1. Phenotypic characterization of the mutant.
(A) Comparison of seed germination among wild-type, mutant, and mutant treated with 10−5 M GA3. (B) Elongation of the second leaf sheath of mutant and wild-type rice in responding to gradient concentration of GAs. (C) Morphological phenotypes of mutant in comparison with wild type. Left and middle panels show the seed size and seedlings, respectively; and the right panel shows the comparison of different internodes. Data are means from n≥25 plants. Scale bars = 1 cm. (D) Microscope images (left panels) and quantitative analysis (right panels) of the second internode (upside panels) and leaf sheath (bottom panels) cell length. “m” stands for the mutant, and “wt” stands for the wild type. Data are means±SE from at least 4 plant sections with at least 6 cells in each section. Significant difference was determined by Student’s t-test (*p < 0.05, **p < 0.01). Scale bars = 10 μm.

Therefore, we treated the mutant with exogenous GA3 to detect it response. Similar with previous report [28], the length of second leaf sheath of mutant was comparable with that of wild-type only with up to 10−6–10−5 M GA3 treatment (Fig 1B) [28]. Then, the mutant seeds were germinated with the presence of 10−5 M GA3, which showed almost the same germination speed and rate with the wild-type (Fig 1A). These results showed that the mutant phenotype could be rescued by exogenous GA3, and it is GA responsive (Fig 1A and 1B). Altogether, it may be a GA deficiency mutant.

Cloning and analyses of the mutated genes through re-sequencing strategy

To verify our hypothesis and further understand the mechanism underlying the delayed germination phenotype of the mutant, we applied whole-genome resequencing and mutmap method to characterize the mutations. The mutant was crossed with Nipponbare wild-type. All individuals of F1 population showed the same phenotype as the wild-type, which indicates that this is a homozygous recessive mutant. In F2 population, we obtained 5342 normal and 798 dwarf individuals showing a segregation ratio of 6.7:1. Considering the fact that there is 45% decrease of germination rate in the mutant, the ratio could be normalized to 3.69:1, indicating the mutant might be caused by a single gene or locus.

As described in M&M, two bulks of DNA from 60 individuals with mutant phenotype and 60 individuals with wild-type phenotype of F2 progeny were sequenced separately. Delta index analysis of the sequencing data showed a peak in a 3.05 Mb region on chromosome six (Fig 2), in which there are 20 genes including three OsKO genes, OsKO1, OsKO2 and OsKO5 (S1 Table). Since we have known that it might be a GA deficient mutant, we then sequenced these three OsKO genes including 2 kb sequence in their promoter regions again. The sequencing results were compared with those of the wild-type from RAP-DB. A short sequence replacement of ‘GACT’ to ‘AATG’ was found in the eighth exon of OsKO1 (starting from +1362) (Fig 2C) without any other detected mutations in the other two OsKO genes. The G to A replacement results in a silent mutation, while the CT to TG replacement leads to a Thr to Met missense mutation (Fig 2C). To confirm the mutation, Taqman probe with 5’ labeled with FAM fluorescent group and 3’ end labeled with BHQ-1 quenching group that completely matched with 18 bp region on mutant CDS containing three mutations (AATG genotype) was used to check the linkage of the mutations to mutant phenotype. The selected 625 mutant individuals in F2 population all matched with probe (AATG genotype), while wild-type individuals (GACT genotype) of the parental lines all did not match with probe, which verified the tight linkage between the mutation and phenotype (Table 1).

Characterization of the mutations through whole-genome resequencing and mutmap methods.
Fig. 2. Characterization of the mutations through whole-genome resequencing and mutmap methods.
(A) flowchart of the sequencing-based method. (B) Whole-genome resequencing of mutant and Nipponbare wild-type segregant F2 population to get casual mutations. Scrolling window plots show SNP index of sequencing results from wild-type pool, mutant pool and SNP index differences with their p values by chi-squared test and fisher test on 12 chromosomes. (C) Structure and mutation sites of OsKO1 gene and its corresponding protein. Black boxes indicate exons, white boxes are UTRs, red fonts are mutated nucleotids and corresponding amino acids.
Tab. 1. Taqman probe detection of linkage between SNP mutations and the mutant phenotype.
Taqman probe detection of linkage between SNP mutations and the mutant phenotype.

To confirm the mutant phenotype is the resultant of the mutations on OsKO1 gene, we conducted genetic complementation analyses in the mutant background. For complementary experiment, two types of transgenic plants were generated. The first one was generated with the vector containing a maize ubiquitin promoter and full length of OsKO1 cDNA from wild-type. The second one was conducted through the vector containing the native promoter (~2 kb) and full ORFs sequence of the gene. More than ten positive transgenic plants were obtained for each type of transformation, of which all rescued the mutant phenotype (Fig 3A–3D). The expressions of OsKO1 in all the transgenic lines were enhanced (Fig 3E). Interestingly, the expression of OsKO1 in the mutant is also higher in mutant than that in wild type (Fig 3E). Altogether, mutations on OsKO1 gene resulted in the mutant phenotype.

Genetic complementation analyses of the mutant.
Fig. 3. Genetic complementation analyses of the mutant.
(A) Germination assay of wild type, mutant, and different transgenic lines. Upside and downside panels are the data for the overexpression (left) and native promoter lines (right) of OsKO1 and OsKO2, respectively. (B) Images showing the germinating seeds from 0 to 48 h of different rice lines. (C) Rescue of the mutant phenotype in OsKO1overexpression lines. (D) Rescue of the mutant phenotype in OsKO1 transgenic lines driven by its native promoter. (E) Expression of OsKO1 gene in different rice lines. (F) Rescue of the mutant phenotype in OsKO2 overexpression lines. (G) Rescue of the mutant phenotype in OsKO2 transgenic lines driven by OsKO1 native promoter. (H) Expression of OsKO2 gene in different rice lines. “wt”, wild-type; m, mutant; 1OE-1,2,3, OsKO1 gene overexpression transgenic lines; 1CP-1,6,7, OsKO1 gene native promoter transgenic lines; 2OE-2,3,4, OsKO2 gene overexpression transgenic lines; 2CP-1,2,3, OsKO2 gene native promoter transgenic lines. Rice OsActin1 was used as internal reference for the qRT-PCR analysis. Data are means±SE (n≥3). Significant difference was determined by Student’s t-test (*p < 0.05, **p < 0.01).

Since OsKOs are P450 family enzymes, the C-terminal region of this family is the heme-binding domain, which is quite conserved [34]. To see the conservative of the mutated site, we than conducted alignment analysis with the sequences from the five OsKO proteins and KO proteins involving in GA biosynthesis from Arabidopsis, wheat and Nelumbo nucifera. The result showed that sequences of these eight KO proteins are highly conserved with more than 70% similarity (S2 Fig). Previous study has shown that OsKO1, OsKO2 and OsKO5 were associated with the dwarfism phenotype, and have closer relationship to each other than with OsKO3 and OsKO4 [29]. Although the mutated site is not highly conserved, it prefers to be an–OH containing amino acid residue (S2 Fig). The replacement of T by M might lead to minor structural changes in the center of enzyme (S1B Fig). Based on these, full-length cDNA sequence of OsKO2 was also applied to conduct the complementation analysis under the same systems as OsKO1. It could also rescue the mutant phenotype. Interestingly, although its rescuing on the germination is weaker than OsKO1 (Fig 3A and 3B), it has a stronger rescuing on the semi-dwarfism phenotype, especially under the native promoter driving (Fig 3F and 3G). The expression of OsKO2 is also higher in the mutant, but with irregularly changes in different transgenic lines (Fig 3H).

OsKO1 contains ent-Kaurene oxidase enzyme activity and involved in GA biosynthesis in rice

Since the exogenous GA could rescue the mutant phenotype, we analyzed the contents of GA1, GA3 and their precursor GA53 in the wild-type and mutant. All of them were almost undetectable in the mutant, much less than that in the wild-type (Fig 4A). Consistent with the genetic complementation analyses, the GA contents could also be recovered in the transgenic mutant lines of both OsKO1 and OsKO2 genes (Fig 4A). Previous studies proved the ent-kaurene oxidase activity of OsKO2, which catalyzes the reaction from ent-kaurene to ent-kaurenol then to ent-kaurenoic acid in GA biosynthetic pathway [27]. Combined all these data, we hypothesize that OsKO1 might also catalyze the same reaction, and the mutation might attenuate its enzyme activity. To verify this hypothesis, transgenic yeast strains expressing wild-type and mutated OsKO1 on pESC-His vector were generated and used to conduct the enzyme assay as described in M&M. In the wild-type system, both ent-kaurenol and ent-kaurenoic acid could be produced, with more ent-kaurenoic acid detected (Fig 4B, S3 Fig). However, much less ent-kaurenoic acid was produced in the mutant system (Fig 4B and 4C). We then analyzed the endogenous content of ent-kaurenoic acid in rice tissues. Although there were undetectable ent-kaurenoic acid in the germinating seeds, it is detectable in the one-month-old seedling. Surprisingly, the content of ent-kaurenoic acid in the mutant is a little bit higher than that in the wild type (Fig 4D). Together, these results indicated that OsKO1 contains the enzyme activity to catalyze the producing of ent-kaurenol and ent-kaurenoic acid in GA biosynthetic pathway. The mutation from Thr to Met did not affect its activity to catalyze the reaction from ent-kaurene to ent-kaurenol, but decreased its activity to catalyze the ent-kaurenol to ent-kaurenoic acid.

GA contents and OsKO1 enzyme activity analyses.
Fig. 4. GA contents and OsKO1 enzyme activity analyses.
(A) Endogenous GA1, GA3 and GA53 contents in the embryos of different rice lines. (B) Yeast expression and GC-MS analysis of enzyme activity of normal OsKO1 (left panel) and mutated OsKO1 (middle panel), with empty vector pESC-His transgenic yeast as negative control (right panel). GC-MS chromatograms at single ion 257 showed products ent-kaurenol (peak 1) and ent-kaurenoic acid (peak 2) catalyzed by OsKO1 from ent-kaurene. The right panel. (C) Quantitative analysis on the production of ent-kaurenoic acid between normal OsKO1 and mutated OsKO1 in the transgenic yeast system. Data are means±SE (n = 3). (D) Endogenous content of ent-kaurenoic acid in one-month-old seedlings of different rice lines. Wild-type (wt), mutant (m) and transgenic lines in mutant background. 1OE, OsKO1 gene overexpression lines; 2OE, OsKO2 overexpression lines; 1CP, OsKO1 gene native promoter transgenic lines; 2CP, OsKO2 gene native promoter transgenic lines. Significant difference was determined by Student’s t-test (*p < 0.05, **p < 0.01).

Subcellular localization and expressional pattern of OsKO1

Previous study in Arabidopsis has shown that the AtKO1 protein localizes on the out envelope membrane of plastids [22]. To determine the subcellular localization of OsKO1, we applied transient transfection assay in Arabidopsis protoplast in parallel with OsKO2. When the full length CDS of OsKO1 and OsKO2 genes were fused with GFP respectively, no GFP signal could be detected in both constructs, although there was strong signal in the cell for the empty pGFP2 vector control (S4 Fig). Previous study has shown that transit peptides at the N-termini for OsCPS1 and OsCPS2 are enough for their transport into plastid [35], we then used the N-terminal 150 amino acids of OsKO1 and N-terminal 140 of OsKO2 proteins containing transit peptides to fuse with GFP, respectively. The signal of OsKO1 could be detected at both plastid and other membrane system, whereas, OsKO2 signal could be only detected at plastid membrane (Fig 5A). This is generally consistent with the result in Arabidopsis [22]. The expressional patterns of OsKO1 and OsKO2 were also analyzed through real-time PCR. Although there are tiny difference, both of them are ubiquitous expressed among different tissues detected in this study with the highest level in leaf (Fig 5B). The expression of OsKO1 was higher than that of OsKO2 in the germinating embryo and stem from young seedlings (Fig 5B). Since OsKO1, OsKO2 and OsKO5 were all shown to be involved in GA biosynthesis, we checked the expression of these three genes in the germinating seeds of both wild-type and mutant at 24 h after imbibition. It is shown that all of them expressed higher in mutant than wild-type, with OsKO1 having the highest expressional level (Fig 5C).

Subcellular localization of OsKO1 and OsKO2, and expression analyses on their encoding genes.
Fig. 5. Subcellular localization of OsKO1 and OsKO2, and expression analyses on their encoding genes.
(A) Subcellular localization of OsKO1 and OsKO2. OsKO1-N150 and OsKO2-N140 indicates the N-terminal 150 and 140 amino acids of OsKO1 and OsKO2, respectively. The N-terminal regions containing transit peptides of OsKO1 and OsKO2 were fused to GFP, and transiently expressed in Arabidopsis mesophyll protoplasts, respectively. Scale bars = 2.5 μm. (B) Expressional patterns of OsKO1 and OsKO2. Rice OsActin1 was used as internal reference for the qRT-PCR analysis. Data are means±SE (n = 3). The 0, 12, 24, 36 and 48 h stand for the embryos from the rice seeds at different time after imbibition, respectively. (C) Comparison of the expressional levels of OsKO1, OsKO2 and OsKO5 in the rice seed at 24 after imbibition between mutant and wild-type.

Pathways affected by the mutation that lead to the delayed germination

To further explore the pathways through which OsKO1 regulate seed germination in rice, RNA-seq was used to profile the differentially expressed genes (DEGs) in seeds at 36 h after imbibition among wild-type, mutant, and mutant treated with 10−5 M GA (S2 Table, S3 Table). For selecting DEGs, the threshold was set at p < 0.05. Totally, 281 genes were differentially expressed between wild-type and mutant, including 174 up-regulated and 107 down-regulated genes in wild-type compared to mutant (S4A Fig); 497 differentially expressed genes were detected between mutant treated with exogenous GA and mutant, including 259 up-regulated and 238 down-regulated genes in mutant treated with exogenous GA compared to mutant (S5A Fig). The reliability of eight freely selected DGEs were further confirmed by qRT-PCR (S5B Fig). Unsurprisingly, the ABA-responsive genes, including those encoding seed storage proteins, seed maturation proteins, and LEAs proteins, had higher expression in mutant compared to both wild-type and mutant treated with GA. To the contrary, genes encoding MYB type transcription factors, sugar and lipid transporters, hydrolase and cytochrome P450 family proteins were down-regulated in mutant (S3 Table). As we have reported, regulation at translational level is more important than that at transcriptional level for rice seed germination [36]. We then compared the proteome profiles of the germinating embryos between wild-type and mutant through a gel-based proteomic strategy. Among the 33 significantly changed proteins, 8 up-regulated proteins in mutant were seed storage proteins (U50, T44, T46) and ABA responsive late embryogenesis abundant proteins (U2, U32, U47, U49, U51) (S5C Fig, S3 Table). We then compared the DEGs list and the changed protein list. Although the number of overlapped genes was very limited, the ABA responsive LEAs and seed maturation proteins existed in both lists, and showed higher expression at both mRNA and proteins levels in mutant than wild-type (S3 Table, S4 Table).

To further confirm the involvement of GA and ABA signaling, we then compared the expression of GA- and ABA-responsive and metabolic genes in the germinating embryos between wild-type and mutant through RT-PCR. Although most of the GA biosynthetic genes were up-regulated, the GA-responsive genes, such as α-amylase and GAMYB, were down-regulated in mutant related to wild-type (Fig 6A), which was similar with that in gd1 and dsg1 mutants [8,15]. GID1, a positive regulator of GA signaling, was also up-regulated in mutant (Fig 6A). Interestingly, ABA biosynthetic genes, OsNCED2 and OsNCED3, were down-regulated, and ABA catabolic gene (OsABA8ox2), playing key roles in ABA catabolism [3740], was up-regulated, although ABA-responsive gene ARAG1 was up-regulated in the mutant (Fig 6A). Considering that the GA induced expression of α-amylase is very important for the degradation of starch in cereal seed, we then detected starch granule in the germinating rice seed through microscopy. Obviously, there were more starch granules in the mutant than wild-type and mutant treated with exogenous GA (Fig 6A). Altogether, it seems the deficiency of GA biosynthesis lead to the delay of ABA signaling attenuation and starch degradation, which might result in the delayed radicle protrusion phenotype in the mutant.

Analyses on the expression of GA and ABA metabolism and responsive genes and the degradation of starch.
Fig. 6. Analyses on the expression of GA and ABA metabolism and responsive genes and the degradation of starch.
(A) Comparison on the expression of GA and ABA metabolism and responsive genes between the embryos from wild-type and mutant through qRT-PCR. Data are means±SE (n = 3). Significant difference was determined by Student’s t-test (*p < 0.05, **p < 0.01). (B) Microscope observation of the starch granules in the paraffin sections of rice embryos from wt, m and m treated with exogenous GA3 (m+GA). The down-side panels are the enlaged rectangle areas of their corresponding up-side panels. Bars = 15μm. All the embryos were from the seeds at 36 h after imbibition.

Discussion

Analysis of the germination delayed dwarfism rice mutant

Seed germination is an important agronomic trait for rice production. Seeds fail to germinate or germinate on mature rice plant (pre-harvest sprouting) will both decrease its yield. As has been reported in Arabidopsis [4143], seed germination is mainly controlled by the homostasis of plant hormones GA and ABA, which has not been fully verified in rice. Recently rice seed germination mutants were characterized based on either T-DNA insertion mutant cloning or map-based cloning methods [8,15,20]. In this study, we identified an mutant showing delayed seed germination and semi-dwarfism phenotype through screening the EMS induced rice mutants. Using whole-genome re-sequencing on F2 progeny bulks, two continuous mutations on OsKO1 gene were detected to be responsible for the delayed germination and semi-dwarfism of mutant. The osko2 mutant was previously reported as d35 mutant [26], with single mutation on the fifth exon leading to semi-dwarfism and an insertion in the 4th exon leading to serious dwarfism or even death [26]. However, involvement of OsKO2 in seed germination was from an indirect study of OsLOL1 [44]. In our study, the two continuous mutations on the eighth exons of OsKO1 lead to a transition from Thr to Met, which results in the delayed seed germination and semi-dwarfism phenotype. The complementary experiments showed that either OsKO1 or OsKO2 gene could rescue the mutant phenotype (Fig 3). Complementary expression of each of these two genes could recover the endogenous GA content as well (Fig 4). In addition, exogenous GA could also rescue the mutant phenotype. Together, it could be concluded that OsKO1 gene could affect seed germination through their function in GA biosynthesis.

Function verification of OsKO1 gene

There are five similar KO genes in rice genome. Phylogenic analysis sorted them into two clades, with OsKO1, OsKO2 and OsKO5 genes belonging to clade I, and OsKO3 and OsKO4 belonging to clade II [26,29]. Among them, OsKO3 and OsKO4 might involve in phytoalexin biosynthesis, and hence contribute to disease resistance [26,29,45]. OsKO2 has been proven to be involved in GA biosynthesis [26]. Studies showed that OsKO1 and OsKO5 might be also involved in GA biosynthesis [26,29,46]. Chen et al. [29] proposed that OsKO1, OsKO2 and OsKO5 might form into heterotrimer to catalyze the reaction in GA biosynthesis. However, previous report [27] and the yeast system enzyme assay in this study showed that both OsKO1 and OsKO2 could independently catalyze the reaction. Furthermore, there is only one KO gene in Arabidopsis genome, it is impossible to function in heterotrimer. Expression analyses of OsKO1 and OsKO2 from different tissues, including root, stem, leaf, flower, embryo at 0, 12, 24, 36, and 48 h after imbibition, showed that OsKO1 is highly expressed at the early stage of seed germination and in the stem of young seedling, which is different with that of OsKO2 (Fig 5). Except for the delayed germination, phenotyping also showed that the semi-dwarfism is mainly ascribed to the shortening of the first and second internodes (Fig 1B). Combining with the fact with OsKO2 rescuing weakly on the germination, it seems that OsKO1 is the major one that functions in GA biosynthesis during seed germination and young seedling stage in rice.

In GA biosynthetic pathway, KO catalyzes the reaction from ent-Kaurene to ent-Kaurenoic acid via ent-Kaurenol and ent-Kaurenal [23,24]. The enzyme activity of OsKO2 has been confirmed in previous study [27]. Consistent with the fact that OsKO1 could rescue the mutant phenotype (Fig 3), it could also catalyze the reaction from ent-Kaurene to ent-Kaurenoic acid in-vitro (Fig 4B and 4C). Since GA content in the OsKO1 complementary transgenic rice was recovered as well (Fig 4A), it could be concluded that OsKO1 also contains the same enzyme activity as OsKO2 and involves in GA biosynthesis. However, the endogenous ent-Kaurenoic acid in mutant is higher than that in the wild type in the one-month-old seedling. It seems that there is compensation in the mutant by OsKO2 after the early seedling growth stage. The mutation of Thr to Met in this mutant resulted in the decrease of its catalyzing efficiency to produce ent-Kaurenoic acid, with no effect on the producing of ent-Kaurenol (Fig 4B). Sequence alignment showed that the amino acids around the mutated site were highly conserved with the mutated site preferring to be Thr (S2 Fig). It seems that each of the reactions from ent-Kaurene to ent-Kaurenol then to ent-Kaurenal and finally to ent-Kaurenoic acid is catalyzed by different motif of KO. Based on the current study, the last conserved motif might not involve in catalyzing the producing of ent-Kaurenol, but in the producing of either ent-Kaurenal or ent-Kaurenoic acid. Further studies are still needed to determine the exact functional motif for each reaction. Taken together, OsKO1 is another functional OsKO gene on early GA biosynthesis pathway especially in ent-kaurenol into ent-kaurenoic acid in rice.

GA and ABA on the control of rice seed germination

GA and ABA are two antagonistic phytohormones, the defect in GA biosynthesis will enhance ABA signaling [47]. It is believed that the balance between GA and ABA rather than their absolute contents controls the seed germination process [4]. LEAs are typical ABA-responsive proteins and involved in many abiotic stress response [48,49]. In rice seed, LEAs along with seed maturation and seed storage proteins are accumulated in responding to the enhancing ABA signal during seed maturation, which will help the seed to survive from desiccation, and maintain its longevity [50,51]. At the early stage of seed germination, these proteins will be degraded to provide precursors and energy for seedling establishment [52]. RNA-seq data showed that the transcripts of these genes were much higher in mutant than wild-type, suggesting the enhanced ABA signaling, which is also supported by the higher level of ARAG1 transcript (Fig 6A) in the mutant. In addition, the degradation of LEAs, seed maturation and seed storage proteins were much slower than that in wild-type (S5C Fig, S4 Table), indicating the weakening of GA signaling, since the mobilization of reserves is promoted by GA. In cereal seed germination, starch degradation catalyzing by α-amylase is a typical activity [16,53], which is obviously weakened in the mutant (Fig 6B). This weakening could be ascribed to the down-regulation on the expression of α-amylase at both mRNA and protein levels (S3 Table, S4 Table), which might also be the result of up-regulation of ARAG1 [54]. Meanwhile, GAMYB is also down-regulated in mutant relative to wild type (Fig 6A, S3 Table). Together, it could be concluded that mutations on OsKO1 and OsKO2 lead to defect in GA biosynthesis, which results in the weakening of GA signaling and enhancing of ABA signaling, and hence delay the germination process. To the contrast of GA and ABA signaling results, most of the GA biosynthesis genes, ABA catabolic gene and GID2, a positive GA regulator [5557] were up-regulated, whereas, ABA biosynthesis genes were down-regulated in the mutant (Fig 6A). It seems there is feedback regulation to re-balance the GA and ABA in the mutant during germination.

Materials & methods

Plant materials, germination condition and GA treatment

The rice germination delayed and dwarfism mutant (hereafter named as mutant) was derived from EMS mutagenesis of Oryza sativa L. cv. Nipponbare (hereafter named as wild-type, and abbreviated as wt) and planted in paddy field in Wuhan, China. For measuring germination rates under GA treatment, 50 seeds of wt and mutant harvested at late September from three biological replicates were germinated on wet filter paper treated with sterilized water or 10−5 M GA (SIGMA, USA) at 25 oC constant temperature incubator with 12 h light (3 Lux)/ 12 h dark cycle, 70% relative humidity. To test whether the mutation is on GA biosynthesis pathway and GA responsive, 20 seeds of wt and mutant were surface sterilized with 3% NaClO for 30 min, then washed with sterilized water for three times. After sterilization, wt and mutant seeds were separately placed on 1% agar treated with 10−12–10−4 M GA. Length of the second leaf sheath was measured after one week as previous described [28].

Whole-genome resequencing of F2 progeny bulks DNA

F2 progenies derived from crossing then self-crossing of mutant with its parental line Nipponbare were bulk sequenced. In details, 60 individuals with mutant phenotype and 60 individuals with wild-type character from the F2 population were sequenced as previously described [58]. A total of 23.13 Gb parental line sequence reads and 24.28 Gb reads of mutant from Illumine Hiseq 2000 PE100 sequencer were controlled by fastqc software to get high quality reads. Clean reads were aligned to reference genome Os-Nipponbare-Reference-IRGSP-1.0 from the Rice Annotation Project by bwa 0.5.9-r16 [59,60]. Then the bam file after alignment was sorted and removed duplications by software samtools 0. 1. 18 [61]. After removing duplications, the file was realigned to reference genome to reduce false positive through GATK Indel Realigner. GATK Unified Genotyper was used to find SNPs and indels, then the SNPs and InDels were filtered through GATK VariantFiltration software to get 200508 filtered SNPs between mutant and Nipponbare [62]. SNP index was calculated according to previous study [63]. For checking linkage of mutations with mutant phenotype, taqman probe 5’-FAM-TGCAGGCAATGCACATCG-BHQ-1-3’ and primer pair 5’-ATGAACAGGAAGGAGTGGGAGT-3’, 5’-TGGAGCTTGTAGGCGGTGAG-3’ (S5 Table) were used. Functional concentration of probe is 1 pmol and was detected by CFX96 Real Time System (BIO-RAD, USA) at 95°C 10 min, 92°C 15 sec, 60°C 1 min, 40 cycles according to previous methods with modifications [64,65].

Sequence analysis

CDS sequence of mutated gene was cloned and aligned to sequence from RAP-DB (http://rapdb.dna.affrc.go.jp). Conserved domain annotation of cytochrome P450 genes were referenced from EMBL-EBI (http://www.ebi.ac.uk/interpro/search/sequence/).

Function identification of OsKO1 through yeast expression system and GC-MS analysis

The ORF of OsKO1 gene from wild-type rice and mutant rice were separately cloned into pESC-His vector under the restriction sites EcoR I and Spe I. The two plasmids OsKO1-wt-pESC-His, OsKO1-m-pESC-His were transformed into WAT11 yeast strain that express the Arabidopsis thaliana reductase NADPH-CYP under the control of a Gal1 promoter using previously reported lithium acetate method. After incubating the single colony of transgenic strain in synthetic dropout medium that lack His amino acid contained 2% glucose to OD 0.8, the expression of OsKO genes were induced by diluted the transgenic strain cultures in the same SD medium contained 2% galactose. The substrate ent-kaurene was added to a final concentration of 100 μM and co-incubated with OsKO or pESC-His empty vector transgenic strains for 72 hours. After extracted with ethyl acetate, the product was dried with N2 gas and then trimethylsilylated with N-methyl-N-trimethylsilyltrifluoroacetamide (MSTFA) and analyzed by GC-MS. Three independent biological replicates with different colonies were conducted for the whole procedure. For quantification of the first product ent-kaurenol, 0.1 μg, 0.5 μg, 1 μg, 5 μg, 10 μg ent-kaurenol standards and their relative GC-MS area were used to depict standard curve, then to calculate product content. To quantify another product ent-kaurenoic acid, 0.1 μg, 0.5 μg, 1 μg, 5 μg, 10 μg ent-kaurenoic acid standards were used and standard curve were depict based on the same method with ent-kaurenol.

GC-MS analysis was carried out on 7890 A gas chromatograph with triple-axis detector coupled with 5975 C mass spectrometer (Agilent, USA). Parameters of mass spectrometer and gas chromatograph were set as: ionization voltage, 70 eV; ion source temperature, 80°C; solvent delay, 4 min; mass to charge ratio range, 40 to 400 amu; flow rate of helium carrier gas, 1 mL/min. Samples were injected into a HP-5 MS column (30 m × 0.25 mm × 0.25 μm film thickness, Agilent) in the splitless mode with outlet temperature 250°C. The column temperature was programmed as: hold at 80°C for 1 min, and a linear gradient of 30°C min-1 to 200°C, from 200°C to 280°C at 5°C min-1 hold at 280°C for 2 min.

Construction of complementation, overexpression vectors for rice transformation

The full length of OsKO1 and OsKO2 CDS under the control of 2075 bp native promoter with 5’ UTR was inserted into binary vector pCAMBIA 1301 by HindIII, KpnI digestion and then transformed into mutant according to previous method with modification [26] to conduct the complementary experiment. For overexpression of OsKO1 and OsKO2 in mutant, cDNAs were cloned into modified pCAMBIA 1301 under the control of maize ubiquitin promoter. Plant transformation was with Agrobacterium strain EHA105 and according to previous study [66]. After using MS medium contained 30 μg/mL hygromycin to select homozygous transgenic plants from T1 generation, phenotypes on the transgenic plants were observed.

Analysis of gibberellins and ent-kaurenoic acid contents

Nipponbare wild-type and mutant seeds freshly harvested at the same time were surface sterilized with 3% NaClO, then germinated on MS medium at 25°C, 70% humidity with 12 h light (3 Lux)/ 12 h dark cycle constant incubator for 36 h. A total of 1 g embryos of wild-type and mutant seeds were excised and snap frozen in liquid nitrogen for GA measurement. Endogenous GA contents of 3 g leaves of OsKO1 overexpression and complementary transgenic lines grow for 64 d selected by 30 μg/mL hygromycin were compared to Nipponbare and mutant. Endogenous GA quantification method was according to previous study [67].

For ent-kaurenoic acid measurement, 100 mg tissues were ground into fine powder in liquid nitrogen, and extracted with methanol containing 1% FA through vortexing 2 min and 12 h at 4°C in darkness. After centrifugation at 12000 g for 10 min, the supernatant was harvested, and subjected to C18 matrix, and then dried. Resuspended it with methanol containing 1% FA, and centrifuged again. The supernatant was harvested and subjected to UHPLC (Thermo Scientific Ultimate 3000) and MS (Thermo Scientific TSQ Quantiva) analysis.

RNA isolation, quantitative real time PCR and digital gene expression profile

Total RNAs of root, stem, leaf, flower, and embryos at 0, 12, 24, 36 and 48 h after germination were extracted with TRIzol Reagent (Invitrogen, Carlsbad, USA) to check the expression profiles of OsKO1 and OsKO2 genes. For qRT-PCR and differentially expressed genes (DEGs) analysis with digital gene expression (DGE) profile, embryos of wild-type, mutant and mutant with 10−5 M GA germinated for 36 hours were used [68,69]. Three independent biological replicates were performed for qRT-PCR and DGE assay, respectively. For qRT-PCR, cDNAs were obtained through Reverse Tra Ace-α-First Strand cDNA synthesis Kit (TOYOBO, Osaka, JAPAN) with 2 μg of total RNA after DNase I digestion. qRT-PCR was performed on CFX96 Real Time System (BIO-RAD, USA) with SYBR green fluorescence according to manufacturer’s instructions. OsACTIN1 was used to normalize the expression ratio of each gene. For DGE sequencing (Illumina Hiseq platform), a total amount of 3 μg RNA was used for each library construction. Reads were mapped to the reference genome with software TopHat v2.0.12. For quantification of gene expression level and differentially expressed genes analysis HTSeq v0.6.1 and DESeq 1.10.1 were used. Quantification of genes used expected fragments per kilobase of transcript per million fragments sequenced (FPKM) based on gene length and reads count mapped to the gene [70]. Genes with p value < 0.05 found in DESeq analysis were assigned as differentially expressed.

Paraffin section and starch granule staining

Leaf sheath and internode of wt and mutant were fixed in FAA solution contained 5% glacial acetic acid, 5% formaldehyde, 70% ethanol until use. The samples were then undergone dehydration, infiltration and embedded in paraffin. After cutting into sections, they were preserved by neutral balsam. For starch granule staining, embryos of wt, mutant and mutant with GA treatment seeds germinated for 36 h were fixed and then embedded in paraffin. Section cutting method, Periodic Acid Schiff staining and specimen observation was according to previous study [36,71]. At least three replicates were fixed and observed for each sample.

Two-dimensional electrophoresis (2-DE)

Four independent biological replicates of mutant and wild-type rice seeds at 36 h after imbibition were sliced to get embryos. After extracting by acetone precipitation method, 1 mg proteins were separated by 17 cm, PH 4–7 IEF gels (BIO-RAD, USA) with isoelectric focusing and SDS-PAGE. Four independent replicate gels were run and analyzed in parallel for each sample. PDQuest 2-DE Analysis Software (BIO-RAD) was used to analyze and select proteins with significant changes (folds change>2, p<0.05). Changed protein spots were digested and desalted, and then peptide solution was identified by a MALDI-TOF mass spectrometer voyager STR (AB, Milwaukee, WI, USA) and searched against NCBI rice database (contains 63553 sequences and 22360987 residues) with Mascot software (http://www.matrixscience.com). Search parameters were set as: type of search, peptide mass fingerprint; enzyme, trypsin; fixed modifications, carbamidomethyl (C); protein mass, unrestricted; mass values, monoisotopic; peptide mass tolerance, ±200 ppm; peptide charge state, 1+; max missed cleavages, 1. Protein score is -10*Log(P), protein scores greater than 61 are significant (p < 0.05). Protein with confidence level > 95% were credible. Protein extraction, 2-DE, in gel digestion and desalting methods were all according to previous study [72].

Subcellular localization assay

ORFs encoding full length of OsKO1 and OsKO2 with its stop codon changed into Arg were fused to GFP under the control of CaMV 35S promoter into pGFP2 vector through XhoI, KpnI sites for OsKO1. The cDNA fragments encoding N-terminal 150 and 140 amino acids of OsKO1 and OsKO2 was cloned into pGFP2 vector through the same restriction endonuclease recognition sites with vectors for ORFs, respectively. Each constructed plasmid was introduced into four-week old mesophyll protoplasts of Arabidopsis wild-type (Col-0) through PEG-calcium mediated transformation to transiently express OsKO-GFP fusion proteins according to reliable guidelines [73]. GFP signal was detected by Leica TCS SP8 laser scanning confocal microscope (Leica, Germany) with same exciting light wavelength 488 nm, detection light wavelength 500 to 530 nm for GFP and 600 to 630 nm for chloroplast.

Supporting information

S1 Table [pdf]
Genes that exist in the region associated with mutant phenotype.

S2 Table [pdf]
Summary of RNA-seq data generated by Illumina Hiseq platform.

S3 Table [pdf]
Significantly differentially expressed genes.

S4 Table [xlsx]
MALDI TOF/TOF MS identified significantly changed proteins.

S5 Table [pdf]
Primers used in this study.

S1 Fig [a]
Distribution and sequence alignment analyses of genes.

S2 Fig [pdf]
Sequence alignment analysis of KO proteins from rice, Arabidopsis, wheat and Nelumbo nucifera.

S3 Fig [a]

S4 Fig [pdf]
Subcellular localization analysis of OsKO1 and OsKO2 with their full length CDS fused with GFP protein.

S5 Fig [a]
RNA-seq and 2-DE analyses on the mutant and wild type rice.


Zdroje

1. Bewley J. Seed germination and dormancy. Plant Cell. 1997; 9: 1055–1066. doi: 10.1105/tpc.9.7.1055 12237375

2. Wu K, Wang J, Kong Z, Ma ZQ. Characterization of a single recessive yield trait mutant with elevated endogenous ABA concentration and deformed grains, spikelets and leaves. Plant Sci. 2011; 180: 306–312. doi: 10.1016/j.plantsci.2010.10.001 21421375

3. Rajjou L, Duval M, Gallardo K, Catusse J, Bally J, Job C, Job D. Seed germination and vigor. Annu Rev Plant Biol. 2012; 63: 507–533. doi: 10.1146/annurev-arplant-042811-105550 22136565

4. Finkelstein R, Reeves W, Ariizumi T, Steber C. Molecular aspects of seed dormancy. Annu Rev Plant Biol. 2008; 59: 387–415. doi: 10.1146/annurev.arplant.59.032607.092740 18257711

5. Gubler F, Millar AA, Jacobsen JV. Dormancy release, ABA and pre-harvest sprouting. Curr Opin Plant Biol. 2005; 8: 183–187. doi: 10.1016/j.pbi.2005.01.011 15752999

6. Nambara E, Marion-Poll A. ABA action and interactions in seeds. Trends Plant Sci. 2003; 8: 213–217. doi: 10.1016/S1360-1385(03)00060-8 12758038

7. Fang J, Chai C, Qian Q, Li C, Tang J, Sun L, et al. Mutations of genes in synthesis of the carotenoid precursors of ABA lead to pre-harvest sprouting and photo-oxidation in rice. Plant J. 2008; 54: 177–189. doi: 10.1111/j.1365-313X.2008.03411.x 18208525

8. Park GG, Park JJ, Yoon J, Yu SN, An G. A RING finger E3 ligase gene, Oryza sativa Delayed Seed Germination 1 (OsDSG1), controls seed germination and stress responses in rice. Plant Mol Biol. 2010; 74: 467–478. doi: 10.1007/s11103-010-9687-3 20878348

9. Rodriguez P, Benning G, Grill E. ABI2, a second protein phosphatase 2C involved in abscisic acid signal transduction in Arabidopsis. FEBS Lett. 1998; 421: 185–190. doi: 10.1016/s0014-5793(97)01558-5 9468303

10. Zhang X, Garreton V, Chua N-h. The AIP2 E3 ligase acts as a novel negative regulator of ABA signaling by promoting ABI3 degradation. Gene Dev. 2005; 19: 1532–1543. doi: 10.1101/gad.1318705 15998807

11. Ding ZJ, Yan JY, Li GX, Wu ZC, Zhang SQ, Zheng SJ. WRKY41 controls Arabidopsis seed dormancy via direct regulation of ABI3 transcript levels not downstream of ABA. Plant J. 2014; 79: 810–823. doi: 10.1111/tpj.12597 24946881

12. Lopez-Molina L, Mongrand S, Mclachlin DT, Chait BT, Chua N-h. ABI5 acts downstream of ABI3 to execute an ABA-dependent growth arrest during germination. Plant J. 2002; 32: 317–328. doi: 10.1046/j.1365-313x.2002.01430.x 12410810

13. Finkelstein R, Wang M, Lynch T, Rao S, Goodman H. The Arabidopsis abscisic acid response locus ABI4 encodes an APETALA2 domain protein. Plant Cell. 1998; 10: 1043–1054. doi: 10.1105/tpc.10.6.1043 9634591

14. Peng J, Harberd NP. The role of GA-mediated signaling in the control of seed germination. Curr Opin Plant Biol. 2002; 5: 376–381. doi: 10.1016/s1369-5266(02)00279-0 12183174

15. Guo X, Hou X, Fang J, Wei P, Xu B, Chen M, et al. The rice GERMINATION DEFECTIVE 1, encoding a B3 domain transcriptional repressor, regulates seed germination and seedling development by integrating GA and carbohydrate metabolism. Plant J. 2013; 75: 403–416. doi: 10.1111/tpj.12209 23581288

16. Kaneko M, Itoh H, Ueguchi-Tanaka M, Ashikari M, Matsuoka M. The alpha-amylase induction in endosperm during rice seed germination is caused by gibberellin synthesized in epithelium. Plant Physiol. 2002; 128: 1264–1270. doi: 10.1104/pp.010785 11950975

17. Kashem MA, Itoh K, Iwabuchi S, Hori H, Mitsui T. Possible involvement of phosphoinositide-Ca2+ signaling in the regulation of alpha-amylase expression and germination of rice seed (Oryza sativa L.). Plant Cell Physiol. 2000; 41: 399–407. doi: 10.1093/pcp/41.4.399 10845452

18. Tamiru M, Undan JR, Takagi H, Abe A, Yoshida K, Undan JQ, et al. A cytochrome P450, OsDSS1, is involved in growth and drought stress responses in rice (Oryza sativa L.). Plant Mol Biol. 2015; 88: 85–99. doi: 10.1007/s11103-015-0310-5 25800365

19. Chen X, Wu Y, Guo J, Du B, Chen R, Zhu L, He G. A rice lectin receptor-like kinase that is involved in innate immune responses also contributes to seed germination. Plant J. 2013; 76: 687–698. doi: 10.1111/tpj.12328 24033867

20. Ye H, Feng J, Zhang L, Zhang J, Mispan MS, Cao Z, et al. Map-based cloning of seed dormancy1-2 identified a gibberellin synthesis gene regulating the development of endosperm-imposed dormancy in rice. Plant Physiol. 2015; 169: 2152–2165. doi: 10.1104/pp.15.01202 26373662

21. Aach H, Bode H, Robinson DG, Graebe JE. ent-Kaurene synthase is located in proplastids of meristematic shoot tissues. Planta. 1996; 202: 211–219.

22. Helliwell CA, Sullivan JA, Mould RM, Gray JC, Peacock WJ, Dennis ES. A plastid envelope location of Arabidopsis ent-Kaurene oxidase links the plastid and endoplasmic reticulum steps of the gibberellin biosynthesis. Plant J. 2001; 28: 201–208. doi: 10.1046/j.1365-313x.2001.01150.x 11722763

23. Sakamoto T, Miura K, Itoh H, Tatsumi T, Ueguchi-Tanaka M, Ishiyama K, et al. An overview of gibberellin metabolism enzyme genes and their related mutants in rice. Plant Physiol. 2004; 134: 1642–1653. doi: 10.1104/pp.103.033696 15075394

24. Yamaguchi S. Gibberellin metabolism and its regulation. Annu Rev Plant Biol. 2008; 59: 225–251. doi: 10.1146/annurev.arplant.59.032607.092804 18173378

25. Magome H, Nomura T, Hanada A, Takeda-Kamiya N, Ohnishi T, Shinma Y, et al. CYP714B1 and CYP714B2 encode gibberellin 13-oxidases that reduce gibberellin activity in rice. Proc Natl Acad Sci USA. 2013; 110: 1947–1952. doi: 10.1073/pnas.1215788110 23319637

26. Itoh H, Tatsumi T, Sakamoto T, Otomo K, Toyomasu T, Kitano H, et al. A rice semi-dwarf gene, Tan-Ginbozu (D35), encodes the gibberellin biosynthesis enzyme, ent-kaurene oxidase. Plant Mol Biol. 2004; 54: 533–547. doi: 10.1023/B:PLAN.0000038261.21060.47 15316288

27. Ko K-W, Lin F, Katsumata T, Sugai Y, Miyazaki S, Kawaide H, et al. Functional identification of a rice ent-kaurene oxidase, OsKO2, using the Pichia pastoris expression system. Biosci Biotech Bioch. 2014; 72: 3285–3288.

28. Li J, Jiang J, Qian Q, Xu Y, Zhang C, Xiao J, et al. Mutation of Rice BC12/GDD1, Which Encodes a Kinesin-Like Protein That Binds to a GA Biosynthesis Gene Promoter, Leads to Dwarfism with Impaired Cell Elongation. Plant Cell. 2011; 23: 628–640. doi: 10.1105/tpc.110.081901 21325138

29. Chen X, Tian X, Xue L, Zhang X, Yang S, Traw MB, Huang J. CRISPR-Based Assessment of Gene Specialization in the Gibberellin Metabolic Pathway in Rice. Plant Physiol. 2019; 180: 2091–2105. doi: 10.1104/pp.19.00328 31160507

30. Helliwell CA, Chandler PM, Poole A, Dennis ES, Peacock WJ. The CYP88A cytochrome P450, ent-kaurenoic acid oxidase, catalyzes three steps of the gibberellin biosynthesis pathway. Proc Natl Acad Sci USA. 2001; 98: 2065–2070. doi: 10.1073/pnas.041588998 11172076

31. Magome H, Yamaguchi S, Hanada A, Kamiya Y, Oda K. Dwarf and delayed-flowering 1, a novel Arabidopsis mutant deficient in gibberellin biosynthesis because of overexpression of a putative AP2 transcription factor. Plant J. 2004; 37: 720–729. doi: 10.1111/j.1365-313x.2003.01998.x 14871311

32. Margis-Pinheiro M, Zhou XR, Zhu QH, Dennis ES, Upadhyaya NM. Isolation and characterization of a Ds-tagged rice (Oryza sativa L.) GA-responsive dwarf mutant defective in an early step of the gibberellin biosynthesis pathway. Plant Cell Rep. 2005; 23: 819–833. doi: 10.1007/s00299-004-0896-6 15668792

33. Sponsel VM, Schmidt FW, Porter SG, Nakayama M, Kohlstruk S, Estelle M. Characterization of new gibberellin-Responsive semidwarf mutants of Arabidopsis. Plant Physiol. 1997; 115: 1009–1020. doi: 10.1104/pp.115.3.1009 9390435

34. Denisov IG, Makris TM, Sligar SG, Schlichting I. Structure and chemistry of cytochrome P450. Chem Rev. 2005; 105: 2253–2278. doi: 10.1021/cr0307143 15941214

35. Toyomasu T, Usui M, Sugawara C, Kanno Y, Sakai A, Takahashi H, et al. Transcripts of two ent-copalyl diphosphate synthase genes differentially localize in rice plants according to their distinct biological roles. J Exp Bot. 2014; 66: 369–376. doi: 10.1093/jxb/eru424 25336684

36. He D, Han C, Yao J, Shen S, Yang P. Constructing the metabolic and regulatory pathways in germinating rice seeds through proteomic approach. Proteomics. 2011; 11: 2693–2713. doi: 10.1002/pmic.201000598 21630451

37. Lefebvre V, North H, Frey A, Sotta B, Seo M, Okamoto M, et al. Functional analysis of Arabidopsis NCED6 and NCED9 genes indicates that ABA synthesized in the endosperm is involved in the induction of seed dormancy. Plant J. 2006; 45: 309–319. doi: 10.1111/j.1365-313X.2005.02622.x 16412079

38. Nambara E, Marion-Poll A. Abscisic acid biosynthesis and catabolism. Annu Rev Plant Biol. 2005; 56: 165–185. doi: 10.1146/annurev.arplant.56.032604.144046 15862093

39. Saika H, Okamoto M, Miyoshi K, Kushiro T, Shinoda S, Jikumaru Y, et al. Ethylene promotes submergence-induced expression of OsABA8ox1, a gene that encodes ABA 8' -hydroxylase in rice. Plant Cell Physiol. 2007; 48: 287–298. doi: 10.1093/pcp/pcm003 17205969

40. Song SY, Dai XY, Zhang W-H. A rice F-box gene, OsFbx352, is involved in glucose-delayed seed germination in rice. J Exp Bot. 2012; 63: 5559–5568. doi: 10.1093/jxb/ers206 22859682

41. Debeaujon I, Koornneef M. Gibberellin requirement for Arabidopsis seed germination is determined both by testa characteristics and embryonic abscisic acid. Plant Physiol. 2000; 122: 415–424. doi: 10.1104/pp.122.2.415 10677434

42. Ogawa M, Hanada A, Yamauchi Y, Kuwahara A, Kamiya Y, Yamaguchi S. Gibberellin biosynthesis and response during Arabidopsis seed germination. Plant Cell. 2003; 15: 1591–1604. doi: 10.1105/tpc.011650 12837949

43. Piskurewicz U, Jikumaru Y, Kinoshita N, Nambara E, Kamiya Y, Lopez-Molina L. The gibberellic acid signaling repressor RGL2 inhibits Arabidopsis seed germination by stimulating abscisic acid synthesis and ABI5 activity. Plant Cell. 2008; 20: 2729–2745. doi: 10.1105/tpc.108.061515 18941053

44. Wu J, Zhu C, Pang J, Zhang X, Yang C, Xia G, et al. OsLOL1, a C2C2-type zinc finger protein, interacts with OsbZIP58 to promote seed germination through the modulation of gibberellin biosynthesis in Oryza sativa. Plant J. 2014; 80: 1118–1130. doi: 10.1111/tpj.12714 25353370

45. Wang Q, Hillwig ML, Wu Y, Peters RJ. CYP701A8: a rice ent-kaurene oxidase paralog diverted to more specialized diterpenoid metabolism. Plant Physiol. 2012; 158: 1418–1425. doi: 10.1104/pp.111.187518 22247270

46. Zhu S, Gao F, Cao X, Chen M, Ye G, Wei C, Li Y. The rice dwarf virus P2 protein interacts with ent-kaurene oxidases in vivo, leading to reduced biosynthesis of gibberellins and rice dwarf symptoms. Plant Physiol. 2005; 139: 1935–1945. doi: 10.1104/pp.105.072306 16299167

47. White CN, Rivin CJ. Gibberellins and Seed Development in Maize. II. Gibberellin synthesis inhibition enhances Abscisic Acid signaling in cultured embryos. Plant Physiol. 2000; 122: 1089–1097. doi: 10.1104/pp.122.4.1089 10759504

48. Khurana P, Vishnudasan D, Chhibbar AK. Genetic approaches towards overcoming water deficit in plants-special emphasis on LEAs. Physiol Mol Biol Plants. 2008; 14: 277–297. doi: 10.1007/s12298-008-0026-y 23572894

49. Lu G, Gao C, Zheng X, Han B. Identification of OsbZIP72 as a positive regulator of ABA response and drought tolerance in rice. Planta. 2009; 229: 605–615. doi: 10.1007/s00425-008-0857-3 19048288

50. Nguyen TP, Cueff G, Hegedus DD, Rajjou L, Bentsink L. A role for seed storage proteins in Arabidopsis seed longevity. J Exp Bot. 2015; 66: 6399–6413. doi: 10.1093/jxb/erv348 26184996

51. Yazaki J, Kikuchi S. The Genomic View of Genes Responsive to the Antagonistic Phytohormones, Abscisic Acid, and Gibberellin. Vitam Horm. 2005; 72: 1–30. doi: 10.1016/S0083-6729(05)72001-X 16492467

52. Yang P, Li X, Wang X, Chen H, Chen F, Shen S. Proteomic analysis of rice (Oryza sativa) seeds during germination. Proteomics. 2007; 7: 3358–3368. doi: 10.1002/pmic.200700207 17849412

53. Fincher GB. Molecular and cellular biology associated with endosperm mobilization in germinating cereal grains. Annu Rev Plant Physiol Plant Mol Biol. 1989; 40: 305–346.

54. Zhao L, Hu Y, Chong K, Wang T. ARAG1, an ABA-responsive DREB gene, plays a role in seed germination and drought tolerance. Ann Bot-London. 2010; 105: 401–409.

55. Dai C, Xue HW. Rice early flowering1, a CKI, phosphorylates DELLA protein SLR1 to negatively regulate gibberellin signalling. EMBO J. 2010; 29: 1916–1927. doi: 10.1038/emboj.2010.75 20400938

56. Hirano K, Kouketu E, Katoh H, Aya K, Ueguchi-Tanaka M, Matsuoka M. The suppressive function of the rice DELLA protein SLR1 is dependent on its transcriptional activation activity. Plant J. 2012; 71: 443–453. doi: 10.1111/j.1365-313X.2012.05000.x 22429711

57. Sasaki A, Itoh H, Gomi K, Ueguchi-Tanaka M, Ishiyama K, Kobayashi M, et al. Accumulation of phosphorylated repressor for gibberellin signaling in an F-box mutant. Science. 2003; 299: 1896–1898. doi: 10.1126/science.1081077 12649483

58. Takagi H, Uemura A, Yaegashi H, Tamiru M, Abe A, Mitsuoka C, et al. MutMap-Gap: whole-genome resequencing of mutant F2 progeny bulk combined with de novo assembly of gap regions identifies the rice blast resistance gene Pii. New Phytol. 2013; 200: 276–283. doi: 10.1111/nph.12369 23790109

59. Li H, Durbin R. Fast and accurate short read alignment with Burrows-Wheeler Transform. Bioinformatics. 2009; 25: 1754–1760. doi: 10.1093/bioinformatics/btp324 19451168

60. Kawahara Y, Bastide M, Hamilton J, Kanamori H, Mccombie W, Ouyang S, et al. Improvement of the Oryza sativa Nipponbare reference genome using next generation sequence and optical map data. Rice. 2013; 6: 1–10. doi: 10.1186/1939-8433-6-1

61. Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N, et al. The sequence alignment/map format and samtools. Bioinformatics. 2009; 25: 2078–2079. doi: 10.1093/bioinformatics/btp352 19505943

62. McKenna A, Hanna M, Banks E, Sivachenko A, Cibulskis K, Kernytsky A, et al. The genome analysis toolkit: a mapreduce framework for analyzing next-generation DNA sequencing data. Genome Res. 2010; 20: 1297–1303. doi: 10.1101/gr.107524.110 20644199

63. Abe A, Kosugi S, Yoshida K, Natsume S, Takagi H, Kanzaki H, et al. Genome sequencing reveals agronomically important loci in rice using MutMap. Nat Biotechnol. 2012; 30: 174–178. doi: 10.1038/nbt.2095 22267009

64. Billard A, Laval V, Fillinger S, Leroux P, Lachaise H, Beffa R, Debieu D. The allele-specific probe and primer amplification assay, a new real-time PCR method for fine quantification of single-nucleotide polymorphisms in pooled DNA. Appl Environ Microb. 2012; 78: 1063–1068.

65. Hui L, DelMonte T, Ranade K. Genotyping Using the Taqman Assay. Curr Protoc Hum Genet. 2008; 56: 2.10.1–2.10.8.

66. Hiei Y, Ohta S, Komari T, Kumashiro T. Efficient transformation of rice (Oryza sativa L.) mediated by Agrobasterium and sequence analysis of boundaries of the T-DNA. Plant J. 1994; 6: 271–282. doi: 10.1046/j.1365-313x.1994.6020271.x 7920717

67. Chen ML, Fu XM, Liu JQ, Ye TT, Hou SY, Huang YQ, et al. Highly sensitive and quantitative profiling of acidic phytohormones using derivatization approach coupled with nano-LC-ESI-Q-TOF-MS analysis. J Chromatogr B. 2012; 905: 67–74.

68. Liu Y, Liu M, Li X, Cao B, Ma X. Identification of differentially expressed genes in leaf of Reaumuria soongorica under PEG-induced drought stress by digital gene expression profiling. PloS one. 2014; 9: e94277. doi: 10.1371/journal.pone.0094277 24736242

69. Yang C, Jiang M, Wen H, Tian J, Liu W, Wu F, Gou G. Analysis of differential gene expression under low-temperature stress in Nile tilapia (Oreochromis niloticus) using digital gene expression. Gene. 2015; 564: 134–140. doi: 10.1016/j.gene.2015.01.038 25617524

70. Trapnell C, Williams BA, Pertea G, Mortazavi A, Kwan G, van Baren MJ, et al. Transcript assembly and quantification by RNA-Seq reveals unannotated transcripts and isoform switching during cell differentiation. Nat Biotechnol. 2010; 28: 511–515. doi: 10.1038/nbt.1621 20436464

71. Li M, Wang K, Wang X, Yang P. Morphological and proteomic analysis reveal the role of pistil under pollination in Liriodendron chinense (Hemsl.) Sarg. PloS one 2014; 9: e99970. doi: 10.1371/journal.pone.0099970 24924488

72. Li M, Sha A, Zhou X, Yang P. Comparative proteomic analyses reveal the changes of metabolic features in soybean (Glycine max) pistils upon pollination. Sex Plant Reprod. 2012; 25: 281–291. doi: 10.1007/s00497-012-0197-0 22968406

73. Yoo SD, Cho YH, Sheen J. Arabidopsis mesophyll protoplasts: a versatile cell system for transient gene expression analysis. Nat Protoc. 2007; 2: 1565–1572. doi: 10.1038/nprot.2007.199 17585298

Štítky
Genetika Reprodukčná medicína

Článok vyšiel v časopise

PLOS Genetics


2020 Číslo 1

Najčítanejšie v tomto čísle

Tejto téme sa ďalej venujú…


Prihlásenie
Zabudnuté heslo

Nemáte účet?  Registrujte sa

Zabudnuté heslo

Zadajte e-mailovú adresu, s ktorou ste vytvárali účet. Budú Vám na ňu zasielané informácie k nastaveniu nového hesla.

Prihlásenie

Nemáte účet?  Registrujte sa