Of UniGenes doi:10.1371/journal.pone.0057715.t001 72,688,546 65,561,528 91,193 47,Total Nucleotides (nt) ?5,900,537,520 ??Mean BI-78D3 biological activity Length (nt) ??414N50 ??868Transcriptome Analysis of Gerbera hybridaFigure 2. Figures of Nr classification. (A) E-value distribution. (B) MedChemExpress 1948-33-0 Similarity distribution. (C) Species distribution. doi:10.1371/journal.pone.0057715.gsteps of GA metabolism more precisely regulate concentrations of bioactive GA. GA20ox and GA3ox are the key enzymes in GA biosynthesis. Both GA20ox and GA3ox were identified in five transcripts of G.hybrida ray florets transcriptome (Table 3, Table S2). In Arabidopsis, GA3ox1 and GA3ox3 in stamen filaments and flower receptacles play major roles in anther development and petal development [31,32]. Emasculation of stamens in petunia (PetuniaTable 2. Summary of the annotations of Gerber hybrida ray floret UniGenes.Number of blasted UniGenes All UniGenes Unigenes of exogenous contaminated species All cleaned UniGenes UniGenes blasted against plant Nr UniGenes blasted against plant Nt UniGenes blasted against Swiss-Prot UniGenes blasted against KEGG UniGenes blasted against GO UniGenes blasted against COG All annotated UniGenes doi:10.1371/journal.pone.0057715.t002 47,104 223 46,881 36,693 28,245 23,040 20,375 15,721 13,239 37,Ratio ??100.00 78.27 60.25 49.15 43.46 33.53 28.23 79.75Transcriptome Analysis of Gerbera hybridaFigure 3. GO categories of the UniGenes. The UniGenes were annotated in three categories: biological processes, cellular components and molecular functions. doi:10.1371/journal.pone.0057715.gFigure 4. COG function classification of UniGenes. doi:10.1371/journal.pone.0057715.gTranscriptome Analysis of Gerbera hybridaTable 3. Statistics of GA metabolism related genes 
in G. hybrida ray florets.SymbolNumber of ECCount of transcriptsDistribution of corresponding hits by local BLASTNG. hybrida `Terra Regina’ A. annyaCPS KS KO KAO GA20ox GA3ox GA2ox 5.5.1.13 4.2.3.19 1.14.13.78 1.14.13.79 1.14.11.12 1.14.11.15 1.14.11.13 1 4 3 3 5 5 6 ?????1 ?5 5 1 ?1 6C. tinctorius??44 ?5 18H. annuus??14 5 8 12 ?doi:10.1371/journal.pone.0057715.thybrida) arrests corolla growth, which can be rescued by exogenous GAs [33]. The expression of GA20ox displayed slight up-regulation at stage 3 and stage 4, which was possibly caused by the petal elongation and expansion (Figure 6). Because of the aborted stamens in ray florets, we speculated that the bioactive GAs are transported from the tiny receptacle under every sole floret or from other places in G. hybrida to the petal to promote its development. However, why the hermaphrodite disc florets located at the same capitulum have extremely short instead of long petals remains unclear. An overdose of GA results in a damaged flower opening and fruit ripening [34]. GA2ox as the major deactivation enzyme are essential for precisely sustaining the optimal bioactive GA 
concentration (Figure 5). Six transcripts of GA2ox were identified in our experiment. The expression of GA2ox displayed tiny upregulation at stage 1 and stage 2. We selected the paralogs from four EST databases, G. hybrida `Terra Regina’, A. annya, C. tinctorius and H. annuus using local BLASTN. The results demonstrated that only one hit of GA3ox was found in G. hybrida `Terra Regina’ and a different number of hits were filtered from the other three species (Table 3). Thus, further study is required on the genes associated with GA biosynthesis.Candidate Genes Related to GA Signal TransductionTh.Of UniGenes doi:10.1371/journal.pone.0057715.t001 72,688,546 65,561,528 91,193 47,Total Nucleotides (nt) ?5,900,537,520 ??Mean Length (nt) ??414N50 ??868Transcriptome Analysis of Gerbera hybridaFigure 2. Figures of Nr classification. (A) E-value distribution. (B) Similarity distribution. (C) Species distribution. doi:10.1371/journal.pone.0057715.gsteps of GA metabolism more precisely regulate concentrations of bioactive GA. GA20ox and GA3ox are the key enzymes in GA biosynthesis. Both GA20ox and GA3ox were identified in five transcripts of G.hybrida ray florets transcriptome (Table 3, Table S2). In Arabidopsis, GA3ox1 and GA3ox3 in stamen filaments and flower receptacles play major roles in anther development and petal development [31,32]. Emasculation of stamens in petunia (PetuniaTable 2. Summary of the annotations of Gerber hybrida ray floret UniGenes.Number of blasted UniGenes All UniGenes Unigenes of exogenous contaminated species All cleaned UniGenes UniGenes blasted against plant Nr UniGenes blasted against plant Nt UniGenes blasted against Swiss-Prot UniGenes blasted against KEGG UniGenes blasted against GO UniGenes blasted against COG All annotated UniGenes doi:10.1371/journal.pone.0057715.t002 47,104 223 46,881 36,693 28,245 23,040 20,375 15,721 13,239 37,Ratio ??100.00 78.27 60.25 49.15 43.46 33.53 28.23 79.75Transcriptome Analysis of Gerbera hybridaFigure 3. GO categories of the UniGenes. The UniGenes were annotated in three categories: biological processes, cellular components and molecular functions. doi:10.1371/journal.pone.0057715.gFigure 4. COG function classification of UniGenes. doi:10.1371/journal.pone.0057715.gTranscriptome Analysis of Gerbera hybridaTable 3. Statistics of GA metabolism related genes in G. hybrida ray florets.SymbolNumber of ECCount of transcriptsDistribution of corresponding hits by local BLASTNG. hybrida `Terra Regina’ A. annyaCPS KS KO KAO GA20ox GA3ox GA2ox 5.5.1.13 4.2.3.19 1.14.13.78 1.14.13.79 1.14.11.12 1.14.11.15 1.14.11.13 1 4 3 3 5 5 6 ?????1 ?5 5 1 ?1 6C. tinctorius??44 ?5 18H. annuus??14 5 8 12 ?doi:10.1371/journal.pone.0057715.thybrida) arrests corolla growth, which can be rescued by exogenous GAs [33]. The expression of GA20ox displayed slight up-regulation at stage 3 and stage 4, which was possibly caused by the petal elongation and expansion (Figure 6). Because of the aborted stamens in ray florets, we speculated that the bioactive GAs are transported from the tiny receptacle under every sole floret or from other places in G. hybrida to the petal to promote its development. However, why the hermaphrodite disc florets located at the same capitulum have extremely short instead of long petals remains unclear. An overdose of GA results in a damaged flower opening and fruit ripening [34]. GA2ox as the major deactivation enzyme are essential for precisely sustaining the optimal bioactive GA concentration (Figure 5). Six transcripts of GA2ox were identified in our experiment. The expression of GA2ox displayed tiny upregulation at stage 1 and stage 2. We selected the paralogs from four EST databases, G. hybrida `Terra Regina’, A. annya, C. tinctorius and H. annuus using local BLASTN. The results demonstrated that only one hit of GA3ox was found in G. hybrida `Terra Regina’ and a different number of hits were filtered from the other three species (Table 3). Thus, further study is required on the genes associated with GA biosynthesis.Candidate Genes Related to GA Signal TransductionTh.