Sequence Variation and Expression Analysis of Seed Dormancy- and Germination-Associated ABA- and GA-Related Genes in Rice Cultivars

Abscisic acid (ABA) and Gibberellic acid (GA) play key roles in regulating seed dormancy and germination. First, when examining germination of different rice cultivars, we found that their germination timing and dormancy status are rather distinct, coupled with different GA/ABA ratio. Second, we studied genomic sequences of ABA and GA dormancy- and germination-associated genes in rice and discovered single nucleotide polymorphisms and insertions/deletions (Indels) in both coding and regulatory sequences. We aligned all these variations to the genome assemblies of 9311 and PA64s and demonstrated their relevance to seed dormancy both quantitatively and qualitatively based on gene expression data. Third, we surveyed and compared differentially expressed genes in dry seeds between 9311 and PA64s to show that these differentially expressed genes may play roles in seed dormancy and germination.

to hormone content, the transition from the dormant to the nondormant state of many seeds is characterized by a decrease in ABA sensitivity and an increase in GA sensitivity (Chiwocha et al., 2005).
Recently, the major enzymes involved in ABA and GA metabolism pathways have been identified. ABA metabolic genes have been identified such as ABA biosynthetic gene NCED that encodes a 9-cis-epoxycarotenoid dioxygenase (NCED) and ABA catabolic gene CYP707A that encodes a ABA 8′-hydroxylase (CYP707A; Seo and Koshiba, 2002;Lefebvre et al., 2006;Okamoto et al., 2006;Yang and Guo, 2007). Genes associated with GA metabolism have also been identified such as GA biosynthetic gene GA20ox that encodes a GA 20-oxidase (GA20ox), GA3ox that encodes a GA 3-hydroxylase (GA3ox), and GA catabolic gene GA2ox that encodes a GA 2-oxidase (GA2ox; Yamauchi et al., 2004;Finch-Savage et al., 2007). However, no potential receptors that could perceive this increase in abscisic acid had been identified until recent reports of four abscisic acid binding proteins: the chloroplast protein magnesium protoporphyrin-IX chelatase H subunit (Shen et al., 2006), the membrane-associated protein G protein coupled receptor 2 (Liu et al., 2007;Verslues and Zhu, 2007), two GPCR-type G proteins (GTG1 and GTG2) (Ma et al., 2009;Pandey et al., 2009), and the PYR/PYL/RCAR family proteins (Ma et al., 2009;Park et al., 2009). Since there have been controversies over structural and functional classification of ABA receptors, we carried out analyses on the recently available recent receptor candidates in this study.

IntroductIon
Seed dormancy has been defined as a temporary failure of seed germination under conditions that favor germination (Bewley, 1997). The initiation of dormancy can either occur in the dry state (after-ripening) or be triggered by imbibition under defined conditions in mature seeds (Gubler et al., 2005). Pre-harvest sprouting of cereal grain leads to reduced grain yield and poor quality products, resulting in economic losses of cultivars. Strong dormancy also prevents grain such as barley from being re-planted immediately after harvest, causing delayed or non-uniform germination. Both situations would result in poor crop establishment and grain performance.
Seed dormancy and germination are regulated by nearly all plant hormones. Several studies have shown that ethylene, auxin, and brassinosteroids promote the germination of dormant seeds, but it is now generally accepted that ABA and GA are the leading regulators (Koornneef et al., 2002;Gubler et al., 2005). ABA deficiency during seed development is associated with the absence of primary dormancy in mature seeds, whereas over-expression of ABA biosynthesis genes can increase seed ABA content and enhance seed dormancy or delay germination. ABA plays important roles in many cellular processes including seed development, dormancy, germination, vegetative growth, and environmental stress responses. ABA level can increase significantly during seed maturation and in response to environmental stresses. In addition Sequence variation and expression analysis of seed dormancy-and germination-associated ABA-and GA-related genes in rice cultivars The identification of GA receptors has also improved our understanding of GA function in relation to germination (Nakajima et al., 2006).
Rice is one of the most important cereal crops, providing food for billions of people. In addition, it is also a model plant with various cultivars and mutants for the study of heterosis and domestic process. 9311 is a cultivar of Oryza sativa L. ssp. indica -the major rice subspecies grown in Asia-Pacific regions. Pei-Ai 64s (PA64s) has a major background of indica and a minor background of japonica and javanica, and the first two are two other commonly cultivated subspecies in China. 9311 and PA64s as parental lines show significant difference in their phenotype of seed dormancy. Their F1 offspring, an elite super hybrid rice LYP9, also shows significant heterosis including weak dormancy like its paternal line 9311. Compared to its wild ancestor O. rufipogon, cultivated rice typically exhibits reduced dormancy (Veasey et al., 2004;Sweeney and Mccouch, 2007). Therefore, to elucidate the molecular mechanisms regulating rice seed dormancy and germination is of importance for both plant biology and crop development. Toward this end, we identified genes encoding ABA and GA metabolic enzymes and found that there were many single nucleotide polymorphisms (SNPs) and Indels in both protein-coding and regulatory regions between different rice cultivars. We also found differentially expressed transcription units in the seed of our rice model 9311 and PA64s. We demonstrate that seed dormancy state in rice may be influenced by both gene sequence variation as well as expression patterns of ABA and GA metabolic genes.

Plant MaterIals
Mature seeds from fresh harvest of 9311, PA64s, IR24, SHOEMED, AZUCENA, DOURADO AGULHA, MAKALIOKA34, and IR36 were harvested at the maturing stage then dried and stored at room temperature. The seeds of 9311 and PA64s were provided by the Chinese National R&D Center on Hybrid Rice, Changsha, China. The seeds of IR24, SHOEMED, AZUCENA, DOURADO AGULHA, MAKALIOKA34, and IR36 were provided by the International Rice Research Institute (IRRI). Six months after harvesting the seeds were used for the experiments. Mature rice embryos were separated manually from their seeds using scalpels, and then stored in −80°C immediately before RNA extraction.

GerMInatIon exPerIMents
For the germination experiment, we started with 50 seeds in a Petri dish with multiple layers of wet filter papers (sterile double-distilled water) at 37°C. The germination rate was calculated daily based on radical emergence. Experiments were performed in triplicate for each cultivar examined.

deterMInatIon of aBa and Ga levels
Mature rice embryos were homogenized in 80% (v/v) acetone. After adding internal standards 13 C-abscisic acid (ABA) and 3 H-labled GA 3 , the homogenate was shaken for 10 h on ice in darkness and then centrifuged at 2,000g for 30 min. The precipitate was then reextracted, and the combined supernatant was evaporated to remove residual acetone. After a series of organic extractions, the extracts were purified through C18 column. ABA was methylated with diazomethane, whereas GA 3 were trimethylsilylated with BSTFA at 100°C for 60 min. Gas chromatography-electron impact ionization mass spectrometry was carried out to determine ABA and GA concentrations. The following mass-to-charge ratio peaks were used for quantification: for ABA,192 (labeled) and 190 (endogenous); and for GA 3 , 506 (labeled) and 504 (endogenous). ABA and GA concentrations of 9311 and PA64s from each individual assay of the triplicate are provided in Data Sheet S1 in Supplementary Material.

teMPlate PreParatIon and real-tIMe Pcr valIdatIon
Total RNA was purified from each sample using Trizol (Invitrogen) according to the manufacturer's instructions. Eight microgram total RNA of each sample was used for first-strand cDNA synthesis in 25 μl reaction containing 5 μl 5 × RT buffer, 2.5 μl 10 mM dNTP, 50 ng random primer, 50 ng oligoDT(15), 2.5 μl RNase inhibitor on the t-test in GraphPad Prism 5. Wilcoxon rank sum test was performed to compare SNPs and Indels between coding regions and regulatory regions, also between paralogs and orthologs.

GerMInatIon test and aBa/Ga concentratIon MeasureMent of rIce seeds froM varIous cultIvars
We initiated our study by carrying out a germination experiment using seeds from various cultivars (9311, PA64s, IR24, SHOEMED, AZUCENA, DOURADO AGULHA, MAKALIOKA34, and IR36) with diverse genetic backgrounds (Figure 1). Although all seeds started to germinate after 1 day rehydration, we observed significant differences among the cultivars. For instance, IR24 showed the highest germination rate (82%) and about half of the seeds from 9311, SHOEMED, IR36, and AZUCENA germinated within 36-h imbibition. In contrast, only 29% of MAKALIOKA34, 14% of DOURADO, and less than 10% of PA64s seeds germinated at this time period.
Since rice germination rate is highly associated with the ABA/GA ratio (Nonogaki et al., 2010), we quantified ABA and GA in 9311 and PA64s dry seeds (see Materials and Methods). The endogenous ABA level of 9311 was obviously lower than those of PA64s and the GA/ABA ratio was much higher in 9311 than in PA64s, which is consistent with our experimental results that PA64s exhibited stronger dormancy than 9311 (Figure 2).

IdentIfIcatIon of rIce aBa and Ga MetaBolIc Genes
Because ABA and GA metabolic genes and their receptors are not wellannotated in the rice genome assemblies, we extracted well-annotated Arabidopsis genes from the public databases and searched for orthologs and paralogs from the rice databases 1 . To confirm the predicted protein (20 U/μl), 4 μl reverse transcriptase (50 U/μl; Invitrogen), 2.5 μl DTT. Reverse transcription was performed at 42°C for 60 min with a final denaturation at 70°C for 15 min.
Specific primers for Real-time PCR in the experiments are listed as follows: OsNCED1 Real-time PCR was performed with Quant SYBR Green PCR kit (BIO-RAD). The data were quantitated by referencing to a PCR product of house-keeping gene of GAPDH under the same reaction condition (94°C for 4 min, then 40 cycles of 94°C for 15 s, 60°C for 20 s, and 68°C for 30 s). Melting curves for each PCR reaction were carefully monitored to avoid non-specific amplifications. Gene expression quantification was performed with the comparative Ct formula normalized with the expression of GADPH.

statIstIc analysIs
We used the unpaired Student's t-test to analyze expression levels of ABA and GA, transcripts of ABA-and GA-related genes in dry seeds when comparing the two groups of 9311 and PA64s, and Ka/ Ks ratios. Statistical significance (p < 0.01) was determined based each other (80% identity) than any other reported NCED proteins from Arabidopsis, were also closely related to AtNCED2, AtNCED3, AtNCED5, AtNCED6, and AtNCED9 (59-63% identities).
Since Arabidopsis CYP707A has four paralogs (AtCYP707A1, AtCYP707A2, AtCYP707A3, and AtCYP707A4), we tried really hard to find each of their orthologs in rice. However, we only identified three rice CYP707A genes; one of them is identical to what was reported previously (Saika et al., 2007). OsCYP707A5, OsCYP707A6, and OsCYP707A7 reside on different chromosomes: 2, 8, and 9, respectively. We observed independent gene duplication events in both rice and Arabidopsis ( Figure 3A).
As to GA3ox in rice, we identified only two previously reported genes , matching to four GA3ox orthologs in Arabidopsis (Hedden et al., 2001). OsGA3ox1 and OsGA3ox2 are located on chromosomes 5 and 1, respectively. The amino acid sequence of OsGA3ox1 and OsGA3ox2 displayed about 37% similarities to Arabidopsis GA3ox proteins. Phylogenetic analysis ( Figure 3B) revealed that GA3ox genes duplicated separately after the monocot and dicot divide.
We identified four OsGA20ox genes; two were previously reported in rice (Sasaki et al., 2002) and five were reported in Arabidopsis (Hedden et al., 2001). In rice, OsGA20ox1, OsGA20ox2, OsGA20ox3, and OsGA20ox4 are located on chromosomes 3, 1, 12, and 5, respectively. The amino acid sequence of OsGA20ox1 has the highest homology with Arabidopsis GA20ox proteins (52-57% identities). In comparison, OsGA20ox2 and OsGA20ox4 have 46-51% identities, respectively. Phylogenetic analysis ( Figure 3B) revealed that GA20ox proteins were subdivided into two groups. OsGA20ox2 and OsGA20ox4 were more closely related (64% identity) than the other OsGA20ox proteins, and they formed one separate branch, but then OsGA20ox1 and OsGA20ox3 were grouped with AtGA20ox1, AtGA20ox2, and AtGA20ox3 to form another.
We identified five GA2ox genes in rice as compared to seven previously reported Arabidopsis counterparts (Schomburg et al., 2003); four (OsGA2ox1, OsGA2ox2, OsGA2ox3, and OsGA2ox4) of them were reported previously (Sakamoto et al., 2001;Sakai et al., 2003). OsGA2ox1, OsGA2ox4, and OsGA2ox5 are located on chromosome 5, and OsGA2ox2 and OsGA2ox3 are on chromosome 1. Although some OsGA2ox genes are localized on the same chromosomes, they are well separated with large phylogenetic distance that suggests early duplication events. OsGA2ox3 showed highest similarity to Arabidopsis GA2ox (49-51% identities), whereas the other OsGA2ox showed identities between 35 and 48% with the Arabidopsis counterparts. Phylogenetic analysis ( Figure 3B) revealed that OsGA2ox1 has lower identities with those of Arabidopsis counterparts, constituting a separate clade.

coMParatIve analysIs of aBa and Ga MetaBolIc Genes aMonG rIce cultIvars
Using a threshold of nucleotide identity over 95% and coverage over 80%, we defined orthologous gene pairs in 9311, PA64s, and Nipponbare (also including a 3-kb upstream sequence containing promoters and UTRs) for comparative analysis ( Table 2). Although 9311, PA64s, and Nipponbare share a common ancestor of O. sativa, there are many sequence variations between orthologs that may be a result of artificial selection. In general, there are significantly more SNPs and Indels in coding regions between paralogs than between orthologs (Wilcoxon rank sum test: p < 0.01). Compared to coding regions, regulatory regions have significantly higher numbers of sequences, we used PFAM database 2 to identify conserved domains and found that rice has five NCED, three CYP707A, four GA20ox, two GA3ox, and five GA2ox gene family members. They are distributed on chromosome 1 (OsGA20ox2, OsGA3ox2, OsGA2ox2, and OsGA2ox3), chromosome 2 (OsNCED1 and OsCYP707A5), chromosome 3 (OsNCED3 and OsGA20ox1), chromosome 5 (OsGA20ox4, OsGA3ox1, OsGA2ox1, OsGA2ox4, and OsGA2ox5), chromosome 7 (OsNCED4), chromosome 8 (OsCYP707A6), chromosome 9 (OsCYP707A7), and chromosome 12 (OsNCED2, OsNCED5, and OsGA20ox3). Our annotations were further confirmed in GBROWSE 3 ( Table 1).

PhyloGenIc analysIs of aBa and Ga MetaBolIc Genes
We identified five NCED genes in rice and nine in Arabidopsis (Schwartz et al., 2003;Tan et al., 2003;Kushiro et al., 2004). OsNCED2 and OsNCED5 are both located on chromosome 12 and OsNCED1, OsNCED3, and OsNCED4 are located on chromosomes 2, 3, and 7, respectively. The amino acid sequence of Arabidopsis AtNCED4 showed higher similarities with OsNCED1 and OsNCED2 (51-57% identities) but relatively lower identities with OsNCED3 to OsNCED5 (ranging between 33 and 36%). Phylogenetic analysis (Figure 3A) revealed that OsNCED3 and OsNCED5, which were more similar to We essentially identified two classes of differences among the three cultivars: SNPs and Indels ( Table 2). Some genes are obviously more conserved (by and large due to selection) than others, such as OsGA3ox2 and OsGA2ox3, who have no Indel or SNP discovered in either 9311 or PA64s (Table 2). One possible reason is that these genes play major functional roles and therefore are strongly selected. Another explanation is that these genes were recently duplicated and therefore less mutation was accumulated. Other genes in 9311, PA64s, and Nipponbare appear high variable, such as OsGA20ox4 (41 SNPs and 2 Indels with a total length of 540 bp) and OsGA2ox5 (65 SNPs and 9 Indels with a total length of 154 bp) ( Table 2). One possible reason for the variability is that they are relative older or less conserved due to weaker selection.
We have further scrutinized some SNPs and Indels of protein domains of the related genes to look for potential changes in protein structure and function. Ninety-seven candidate SNPs and 618 candidate Indels involving amino acid changes were found to lie in domains of RPE membrane protein, cytochrome P450, 2OG-Fe(II) oxygenase superfamily, lanthionine synthetase C-like protein, and alpha/beta hydrolase fold, as well as other residues associated with specific functions predicted by pfam (Table 3). We found that 2OG-Fe(II) oxygenase of OsGA3ox family is highly conserved but alpha/beta hydrolase fold of OsGID1 has 24 SNPs and 172 amino acid Indels (Figure 5). Compared to the sequences of PA64s and substitutions and Indels (Wilcoxon rank sum test: p < 0.01), which may also collectively contribute to the differential expression of these genes.
Based on the degree of sequence variations, we can readily identify the gene family members and their relationship with regard to obvious inheritances from their ancestor subspecies, indica or japonica. For instance, in NCED family, three genes (OsNCED1, OsNCED3, and OsNCED5) of PA64s and 9311 were inherited from the same indica ancestor as shown in Figure 4, which is consistent with lower variation between 9311 and PA64s ( Table 2). In another gene family, OsCYP707A6 and OsCYP707A7 of PA64s were from indica ( Figure 4B). In GA20ox family, PA64s inherited OsGA20ox2 from indica, whereas other genes were not obviously separated into indica or japonica ( Figure 4D). In GA3ox family, both OsGA3ox1 and OsGA3ox2 were from indica background, and in GA2ox family, PA64s inherited OsGA2ox4 from its indica background (Figures 4E  and 4F). ABA receptor (OsGPCR) of PA64s was most likely from japonica or javanica rather than from indica background. Both ABA and GA receptors were found with large amount of variations between 9311 and PA64s due to frequent hybridization with complex genetic background during rice domestication and breeding. We concluded that the dormancy trait of PA64s was mainly inherited from indica, though it also maintained a partial background of japonica and javanica. The mixed background may contribute to its significant difference in germination when compared with 9311. Arabidopsis and rice. We calculated non-synonymous and synonymous substitution rate (Ka and Ks) for each pair (Figure 6). As indicated by Ka/Ks < 1, all of pairwise comparisons appear to undergo purifying selection. Interestingly, the Ka/Ks ratio of OsGA20ox2 is significantly different from other members of OsGA20ox family (t-test: p < 0.01), indicating that it may evolve faster than others. This difference may be either a result of positive selection or relaxed purifying selection. Our results suggest that genes with multiple copies tend to evolve in different patterns, consistent with the view that one of the duplicates may undergo purifying selection after gene duplication while the other may enjoy more relaxed selective pressure. Furthermore, we also found that some of these genes (such as OsNCED1 and OsGA3ox2) had different Ka/Ks values between different rice cultivars. In imbibing seeds, ABA and GA contents are regulated by ABA biosynthetic enzyme NCED, ABA catabolic enzyme CYP707A, GA biosynthetic enzymes GA20ox and GA3ox, and GA catabolic enzyme GA2ox. We detected expression levels of these genes which are important for the regulation of ABA and GA contents Welsch et al., 2008). We also performed real-time PCR to confirm the expression levels of 19 genes in 9311 and PA64s, including 6 genes (OsNCED3, OsCYP707A5, OsGA20ox2, OsGA3ox2, OsGA2ox1, and OsGA2ox5) high-abundance and 13 (OsNCED1,OsNCED2,OsNCED4,OsNCED5,OsCYP707A6,OsCYP707A7,OsGA20ox1,OsGA20ox4,OsGA20ox3,OsGA3ox1,OsGA2ox2,OsGA2ox3,OsGA2ox4) low-abundance transcripts. Our results showed that the expression of ABA and GA metabolic genes was relative low and that there were significant differences between 9311 and PA64s (Figure 7). For instance, the expression of ABA synthetic gene OsNCED3 was 3.9-folds higher (t-test: p < 0.01) and catabolic gene OsCYP707A5 was 3.5-folds lower (t-test: p < 0.01) in PA64s than in 9311; this result implied a lower ABA level in 9311 (Barrero et al., 2006). The GA synthetic gene OsGA3ox2 was expressed at a much higher level (2.8-folds) in 9311 than in PA64s. In contrast, both of GA catabolic genes OsGA2ox1 and OsGA2ox5 were expressed at a much higher level in PA64s than in 9311; this result suggested higher GA level or less dormant seeds in 9311. Indeed, a correlation between the expression level of GA metabolism genes and the content of active GA has been reported (Busov et al., 2003;Oh et al., 2006). The differential expression among gene family members also confirmed the complexity of the regulatory network, also footnoted by higher variations in regulatory regions.

dIscussIon
Studying the plasticity of dormancy and germination is of necessity in plant biology and crop breeding (Roberts, 1961). We engaged our study on seed germination and discovered that dormancy states are different among rice cultivars, consistent with previous publications (Seshu and Dadlani, 1991;Gu et al., 2004Gu et al., , 2006Veasey et al., 2004). Aligning sequences in coding and regulatory regions, we observed large amount of Indels and SNPs in genes that are related to regulate the ABA and GA metabolic pathways. For instance, two Indels of 306 bp in the coding sequence and two Indels of 260 bp in the upstream regulatory region of OsCYP707A5 between 9311 Nipponbare, 9311 has a large insertion in the alpha/beta hydrolase fold (Ollis et al., 1992), and this insertion may lead to enhanced sensitivity to GA. Sequence polymorphisms in the functional domains that differ between the three subspecies may be one of the direct causes for functional variation, although further functional tests are required to confirm their precise roles.

dIfferent functIonal constraInts on aBa and Ga MetaBolIc Genes Between 9311, PA64s, and NiPPoNbAre
Following a gene duplication event, the two duplicates may be subjected to different selective constraints and even new functions. To investigate whether or not different rice cultivars have undergone different types of artificial selection, we identified reciprocal best matches -orthologs -among ABA and GA metabolic genes from  the result suggests that they are most likely undergoing purifying selection despite the fact that most of them belong to multi-gene families. Compared to other OsNCED family members, OsNCED3, due to their high expression levels, may play a major role in controlling dormancy and germination (Lefebvre et al., 2006;Hwang et al., 2010), which is consistent with that the ectopic expression of OsNCED3 in Arabidopsis leads to a delay in seed germination (Hwang et al., 2010). Expression analysis of CYP707A genes suggests that OsCYP707A5 may play a key role in seed dormancy and germination (Kushiro et al., 2004;Millar et al., 2006;Okamoto et al., 2006). One SNP and several Indels with a total length of 59 amino acids were identified in functional regions of OsCYP707A5, indicating its potential involvement in regulating ABA concentration.
and PA64s may affect the function of OsCYP707A5 and the binding of its regulatory transcriptional factors. In addition, we compared genes differentially expressed in dry seed between 9311 and PA64s, and hypothesized that different cellular concentrations of ABA and GA may reflect the genetic variations among these genes, especially those related to the ABA and GA pathways Welsch et al., 2008). As rice cultivars have been under strong artificial selection during domestication and the loss of seed dormancy is one of the most important traits in the domestication syndrome, dormancy-and germination-associated genes are absolutely the targets of selection. Thus, we further calculated Ka/Ks values for the ABA and GA metabolism-related genes and found that all were lower than 0.5;  rice cultivars exhibit different dormancy states due to differential expression of ABA and GA metabolic genes. Environmental conditions, such as high temperature and light (Toh et al., 2008), can enhance expression of specific ABA biosynthesis genes (OsNCED family) and GA catabolic genes (OsGA2ox family), resulting in seed dormancy through induction of ABA dominance. Low temperature (Finch-Savage et al., 2007) can also enhance expression of specific ABA catabolic genes (OsCYP707A family) and GA biosynthesis genes (OsGA20ox and OsGA3ox families). In order to further illustrate molecular mechanisms of embryonic dormancy in rice, additional studies should be carried out, focusing on comparative sequence analysis and expression profiling of dormancyassociated genes among large samples of rice cultivars to construct gene regulation networks covering all the stages of seed dormancy and germination.
In OsGA20ox and OsGA3ox families, OsGA20ox2 and OsGA3ox2 may play major roles in the regulation of GA synthesis (Calvo et al., 2004); both are highly conserved as no SNP/Indel was found in their functional regions. Therefore, their differential expression may control the cellular concentration of GA (Mitchum et al., 2006;Rieu et al., 2008). By the same token, OsGA2ox1 and OsGA2ox5 of OsGA2ox family may play major roles in regulating GA catabolism (Mitchum et al., 2006) so do OsGPCR and OsGID1 (Aleman et al., 2008;Hirano et al., 2008;Murase et al., 2008). A balanced ABA/GA ratio governed by dynamics of hormone synthesis and catabolism regulates dormancy and germination through changing seed sensitivity to the external environment (Cadman et al., 2006;Seo et al., 2006). We propose a model extending the ABA/GA balance theory to address changes in gene function and expression (Figure 8). According to this model, different  FiGure 7 | real-time PCr based analyses of OsNCED3, OsCYP707A5, OsGA20ox2, OsGA3ox2, OsGA2ox1, and OsGA2ox5 discovered in 9311 and PA64s seeds. Transcripts were normalized with GAPDH as an internal control. Measurements were averaged over three replicates with error bars representing SE.
FiGure 8 | A schematic model for the regulation of dormancy and germination by ABA and GA in response to the environment in rice seeds. Under the same environmental conditions, different rice cultivars exhibit different dormancy states due to their differential expression of ABA and GA metabolic genes. Environmental conditions such as higher temperature enhance expression of specific ABA biosynthesis genes (OsNCED family) and GA catabolic genes (OsGA2ox family), which result in dormancy by inducing dominance of ABA. Moreover, whether a seed remains dormancy or germinates may depend on a balanced ABA and GA metabolism influenced by genetic and epigenetic mechanisms, a collective outcome determined by sequence and expression variations of ABA and GA metabolic genes and the function of ABA/GA receptors.

conclusIon
In this study, we identified rice ABA and GA metabolism-related genes based on comparative and phylogenetic analyses among rice cultivars. We found that there are many SNPs and Indels in the coding and regulatory sequences of these genes due to their redundancy and diverse genetic background. We further traced the variations on the genome sequences of 9311 and PA64s and showed that sequence variations may lead to functional variations and variable seed dormancy states regulated by ABA and GA. In addition, we also surveyed and compared differentially expressed genes between 9311 and PA64s and demonstrated that differential expression of related genes may also play roles in the variable dormancy and germination states in rice. Although precise correlation between sequence variations and dormancy states need further experimental verification, sequence and expression analyses of ABA and GA metabolism-related genes among rice cultivars pave a way to study molecular mechanism of seed dormancy and germination.

acknowledGMents
The work was supported by grants from the National Natural Science Foundation of China (30221004)