# Generation and Analysis of Expressed Sequence Tags fromChimonanthus praecox (Wintersweet) Flowers for Discovering Stress-Responsive and Floral Development-Related Genes

Authors: Shunzhao Sui; Jianghui Luo; Jing Ma; Qinlong Zhu; Xinghua Lei; Mingyang Li
Journal: Comparative and Functional Genomics
 (2012)

---

## Abstract

A complementary DNA library was constructed from the flowers ofChimonanthus praecox, an ornamental perennial shrub blossoming in winter in China. Eight hundred sixty-seven high-quality expressed sequence tag sequences with an average read length of 673.8 bp were acquired. A nonredundant set of 479 unigenes, including 94 contigs and 385 singletons, was identified after the expressed sequence tags were clustered and assembled. BLAST analysis against the nonredundant protein database and nonredundant nucleotide database revealed that 405 unigenes shared significant homology with known genes. The homologous unigenes were categorized according to Gene Ontology hierarchies (biological, cellular, and molecular). By BLAST analysis and Gene Ontology annotation, 95 unigenes involved in stress and defense and 19 unigenes related to floral development were identified based on existing knowledge. Twelve genes, of which 9 were annotated as “cold response,” were examined by real-time RT-PCR to understand the changes in expression patterns under cold stress and to validate the findings. Fourteen genes, including 11 genes related to floral development, were also detected by real-time RT-PCR to validate the expression patterns in the blooming process and in different tissues. This study provides a useful basis for the genomic analysis of C. praecox.


---

## Body

## 1. Introduction

Chimonanthus praecox(L.) Link, wintersweet, belongs to the Calycanthaceae family. It is a perennial deciduous shrub and blossoms in winter, from late November to March. Its unique flowering time and long blooming period make it one of most popular ornamental plants in China [1]. C. praecox is mainly a garden plant that also provides cut flowers. The flower is strongly fragrant and may be used as a source of essential oil, which has received much attention in New Zealand [2].C. praecox thrives in cold environments and blooms in low-temperature seasons with little rainfall. The plant is assumed to be rich in genes related to floral development and adversities, especially those responding to environmental stress factors. However, the molecular mechanism that regulates floral development and copes with stresses in C. praecox flowers remains unclear.Expressed sequence tags (ESTs) have been proven to be an efficient and rapid means to identify novel genes (and proteins) induced by environmental changes or stresses [3–7]. Genes related to flower form, longevity, and scent from roses, Phalaenopsis equestris, and Pandanus fascicularis were identified by ESTs [8–10]. The present study used transcriptomic analysis of C. praecox flowers to identify novel genes induced by environmental changes or related to floral development and ultimately to better understand the physiological and genetic basis of cold acclimation in flowers of woody plants.


## 2. Materials and Methods

### 2.1. Complementary DNA (cDNA) Library Construction and Sequencing

The process ofC. praecox blossoming includes the following: Stage 1, sprout period; Stage 2, flower-bud period; Stage 3, display-petal period; Stage 4, initiating bloom period; Stage 5, bloom period; Stage 6, wither period [11]. C. praecox flower buds or flowers at the six stages of development (Figure 1) were collected from the nursery at Southwest University, Chongqing, China, for cDNA library construction. The samples were immediately frozen with liquid nitrogen and then refrigerated at −80°C until RNA isolation. Total RNA samples were extracted from these flower buds or flowers using RNA isolation kits (W6771; Watson Biotechnologies Inc.), and RNA quality was detected by BioSpec-mini (Shimadzu). The final RNA sample for the cDNA library construction was bulked by pooling equal amounts of total RNA from each stage.Figure 1
Stages ofC. praecox blooming.  Stage 1, sprout period: bud scales loosen; Stage 2, flower-bud period: flower buds turn green; Stage 3, display-petal period: flower buds enlarge and turn yellow; Stage 4, initiating bloom period: fragrance emerges; Stage 5, bloom period: flowers fully open with strong fragrance; Stage 6, wither period: petals begin withering.cDNA synthesis from the mixed total RNA and library construction through directional cloning of cDNAs into theλTriplEx2 vector was performed using a SMART cDNA Library Construction Kit (K1051-1) according to the manufacturer’s instructions. The size of the insert fragment and the recombinant rate were measured by PCR, as described by Gao and Hu [12], using a random selection of 50 clones. All clones were sequenced from the 5′ end using ABI 3700 (Applied Biosystems).


### 2.2. Sequence Clustering, Annotation, and Functional Categorization

Poor-quality sequences, sequences with less than 100 bases, and vector sequences were trimmed from the raw sequences using SeqMan II (DNASTAR, Inc., Madison, WI) and manually. The trimmed cDNA sequences were assembled into clusters using the assembly program within SeqMan II set to default parameters. The assembly parameters were set to require a minimum match of 80% over 12 bp to initiate the assembly process. A consensus sequence for unigenes was exported from SeqMan II.Unigenes (contigs and singletons) were annotated using BLASTX against the NCBI nonredundant protein database with a cut-offE value of the best hit of ≤10−5 [13]. Sequences without a reliable match (>10−5) were subsequently compared with the NCBI nonredundant nucleotide database by performing BLASTN (score >100) for complementary annotation [14]. All well-annotated unigenes were then further classified and mapped to the three Gene Ontology (GO) categories (biological, cellular, and molecular) via AmiGO (http://amigo.geneontology.org/) [15].


### 2.3. Expression Analysis of Selected Genes Using Quantitative Real-Time RT-PCR

For real-time expression studies,C. praecoxseeds were kept under a 16 h light/8 h dark cycle in a growth chamber at 25°C. Small seedlings were subsequently transferred to plastic pots containing a mixture of soilrite and vermiculite (1 : 1). The plantlets were raised to the sixth leaf stage and then subjected to cold treatment (4°C) for 15 min, 1 h, and 6 h. Control plantlets were maintained at 25°C. The tissues were harvested and snap-frozen in liquid nitrogen and kept at −80°C until further use.The total RNA of flowers at the different developmental stages as well as floral organs dissected from Stage 4 flowers, roots, stems, and leaves were extracted with the above-mentioned method. Equal amounts of DNA-free RNA (5μg) from different tissues were reverse-transcribed using a PrimeScript RT Reagent Kit with gDNA Eraser (TaKaRa). The primers used for real-time PCR (Table 1) were designed by Primer Premier 5.0 (PREMIER Biosoft, USA). Real-time PCRs were performed in triplicate on 10 μL reactions using 5 μL of SsoFast EvaGreen Supermix (Bio-Rad), 0.5 μL of each primer (final concentration of 500 nM), 3.5 μL of water, and 0.5 μL of cDNA template. Table 1 shows the primer pair-specific temperatures. PCR was carried out using a Bio-Rad CFX96 real-time system machine based on the manufacturer’s instructions under the following conditions: denaturation at 95°C for 30 s; 40 cycles of 95°C for 5 s; 59°C for 5 s. Dissociation kinetics was performed using the real-time PCR system at the end of the experiment (60 to 95°C; continuous fluorescence measurement) to check the annealing specificity of the oligonucleotides. A comparative Ct (threshold cycles) method of relative quantification was used to analyze the real-time quantitative RT-PCR data using Bio-Rad CFX Manager Software Version 1.6. Actin and tubulin were used as housekeeping genes for the calculation of relative transcript abundance. The sizes of the amplified products were confirmed via gel electrophoresis. Negative controls with no templates were carried out concurrently.Table 1
Forward (F) and reverse (R) primers used in real-time PCR.
Gene codeAnnotation (sequence homology)Primer sequence (5′→3′)Tm (°C) aCyclesProduct lengthCp88Membrane channel proteinF: ATGTGGACTTTTTGCTGGAGTTTTT594090R: GCAATACAAGCATTTACCCAGTCCTCp173Low temperature and salt responsive proteinF: CTGACGCTCTTTGGATGGCTACC5940172R: ACGAATCAACGGTCACAAACACGCp215Putative abscisic acid-induced proteinF: TGCTCTTTACTTGCTGGCTCGTCT5940193R: TCTCGCCCATTCTCTTCTATGTATTTCCp274Glutathione S-transferaseF: ATGTGCTAAATCTCATTCCAGTCTTCC594085R: GCCGTGATATCCTGCGAAACTAGCp3591,4-Alpha-glucan-maltohydrolaseF: TATTATCTGGCAAGGTTCCACTCGTT594086R: TGTGTTGTAGAACCCAGCAGTCAGCCp364Cold acclimation Protein WCOR413-like protein beta formF: TCCAATAAGCCGAATCAAAATACAGT594084R: GTCTGTCTTCATCGCCAAATACTCCCp3753-Oxoacyl-[acyl-carrier-protein] synthaseF: GTAACTGCGCCGAAACGAGAAAAAG594078R: CCCAAAGACAGAAACCAGCCCCCp440Putative multiple stress-responsive zinc-finger proteinF: GTTTCCCGATTTTGTGTGGCG5940183R: AGTTGTTGATGCAGAGAATTGGTCCTCp465CatalaseF: CATCGTTCGCTTCTCTACCGTTATT5940118R: TGTTTCCGACCAGATCAAAGTTACCCp197AGL9.2F: TGAGAATAAGATAAATCGGCAGGTGAC5940128R: TTTTCCTCTGCTGGAGAAGACAACGCp203LLP-B3 proteinF: GCGATGTGAAACTTGTCAGTAGCCC5940171R: AGTTTGGCACAATCAGGACGAGGTTCp268Caffeoyl-CoAO-methyl-transferase-like proteinF: AACACAAATGGAGAAGAGAAAACCC5940193R: CATCGGCAGAGGTCGTCATAAGACp297SRG1-like proteinF: TAAGAACACCGTCCCATCTCGCTAC5940143R: GAGTGCAGTCTCTCCAATTCATCAGCp328STYLOSA proteinF: GGGGAACCAACAAGAACAACAGAT5940161R: GCCCAGTGGCACCATTTAGAAGATCp330RAC-like G-protein Rac1F: ACCGCTTCTACCGCCTTTTCTCC5940166R: AGGTCTTTCCAACAGCACCATCTCCCp360Peroxisomal fatty acid betao-xidation multifunctional protein AIM1F: TATGCTGGTGTTTTGAAAGAGGCG5940193R: CTATGCCAGACCCCATCAAACCTCCp383Secondary cell wall-related glycosyltransferase family 47F: CAGGCGACACCCCCTCCTCTAAC5940150R: TCAAGGCATCAGAGTTACGCACAAAGCp423Putative polyketide synthaseF: GAGTTCAGACACTGTATGCCCTTGG5940165R: TTGCTTGTAAAATTCGTCCTTCGTTCp436Allergen-like protein BRSn20F: CTCTTATCGCTCTCTGCTTCCTTCC5940176R: GTTGCTGTTAACGTCTCTGCATTCCCp458MADS-box protein 9F: TGGTTGAGTATGATGTAGAGAAGCGAC594098R: ACCCATCATCTGAAGGTGTGCTATTCp82Dormancy related proteinF: TGGAATAGAAAGAGCAAATACAGCG5940172R: TACAAAAGGCTCAATGGCGTCCCp24No hitsF: GCAGTTTACATACTACGGGAGAGGCT5940106R: CGGCTTACGGAATCGTCATCACCp64Putative proteinF: CCAATCACTCTCCCTGAGGATGTAT5940159R: TTCACCGACTCCTTGTTCTTTTAGCInternal controlActinF: GTTATGGTTGGGATGGGACAGAAAG5940199R: GGGCTTCAGTAAGGAAACAGGAInternal controlTublinF: TAGTGACAAGACAGTAGGTGGAGGT5940139R: GTAGGTTCCAGTCCTCACTTCATCaAnnealing temperature.


## 2.1. Complementary DNA (cDNA) Library Construction and Sequencing

The process ofC. praecox blossoming includes the following: Stage 1, sprout period; Stage 2, flower-bud period; Stage 3, display-petal period; Stage 4, initiating bloom period; Stage 5, bloom period; Stage 6, wither period [11]. C. praecox flower buds or flowers at the six stages of development (Figure 1) were collected from the nursery at Southwest University, Chongqing, China, for cDNA library construction. The samples were immediately frozen with liquid nitrogen and then refrigerated at −80°C until RNA isolation. Total RNA samples were extracted from these flower buds or flowers using RNA isolation kits (W6771; Watson Biotechnologies Inc.), and RNA quality was detected by BioSpec-mini (Shimadzu). The final RNA sample for the cDNA library construction was bulked by pooling equal amounts of total RNA from each stage.Figure 1
Stages ofC. praecox blooming.  Stage 1, sprout period: bud scales loosen; Stage 2, flower-bud period: flower buds turn green; Stage 3, display-petal period: flower buds enlarge and turn yellow; Stage 4, initiating bloom period: fragrance emerges; Stage 5, bloom period: flowers fully open with strong fragrance; Stage 6, wither period: petals begin withering.cDNA synthesis from the mixed total RNA and library construction through directional cloning of cDNAs into theλTriplEx2 vector was performed using a SMART cDNA Library Construction Kit (K1051-1) according to the manufacturer’s instructions. The size of the insert fragment and the recombinant rate were measured by PCR, as described by Gao and Hu [12], using a random selection of 50 clones. All clones were sequenced from the 5′ end using ABI 3700 (Applied Biosystems).


## 2.2. Sequence Clustering, Annotation, and Functional Categorization

Poor-quality sequences, sequences with less than 100 bases, and vector sequences were trimmed from the raw sequences using SeqMan II (DNASTAR, Inc., Madison, WI) and manually. The trimmed cDNA sequences were assembled into clusters using the assembly program within SeqMan II set to default parameters. The assembly parameters were set to require a minimum match of 80% over 12 bp to initiate the assembly process. A consensus sequence for unigenes was exported from SeqMan II.Unigenes (contigs and singletons) were annotated using BLASTX against the NCBI nonredundant protein database with a cut-offE value of the best hit of ≤10−5 [13]. Sequences without a reliable match (>10−5) were subsequently compared with the NCBI nonredundant nucleotide database by performing BLASTN (score >100) for complementary annotation [14]. All well-annotated unigenes were then further classified and mapped to the three Gene Ontology (GO) categories (biological, cellular, and molecular) via AmiGO (http://amigo.geneontology.org/) [15].


## 2.3. Expression Analysis of Selected Genes Using Quantitative Real-Time RT-PCR

For real-time expression studies,C. praecoxseeds were kept under a 16 h light/8 h dark cycle in a growth chamber at 25°C. Small seedlings were subsequently transferred to plastic pots containing a mixture of soilrite and vermiculite (1 : 1). The plantlets were raised to the sixth leaf stage and then subjected to cold treatment (4°C) for 15 min, 1 h, and 6 h. Control plantlets were maintained at 25°C. The tissues were harvested and snap-frozen in liquid nitrogen and kept at −80°C until further use.The total RNA of flowers at the different developmental stages as well as floral organs dissected from Stage 4 flowers, roots, stems, and leaves were extracted with the above-mentioned method. Equal amounts of DNA-free RNA (5μg) from different tissues were reverse-transcribed using a PrimeScript RT Reagent Kit with gDNA Eraser (TaKaRa). The primers used for real-time PCR (Table 1) were designed by Primer Premier 5.0 (PREMIER Biosoft, USA). Real-time PCRs were performed in triplicate on 10 μL reactions using 5 μL of SsoFast EvaGreen Supermix (Bio-Rad), 0.5 μL of each primer (final concentration of 500 nM), 3.5 μL of water, and 0.5 μL of cDNA template. Table 1 shows the primer pair-specific temperatures. PCR was carried out using a Bio-Rad CFX96 real-time system machine based on the manufacturer’s instructions under the following conditions: denaturation at 95°C for 30 s; 40 cycles of 95°C for 5 s; 59°C for 5 s. Dissociation kinetics was performed using the real-time PCR system at the end of the experiment (60 to 95°C; continuous fluorescence measurement) to check the annealing specificity of the oligonucleotides. A comparative Ct (threshold cycles) method of relative quantification was used to analyze the real-time quantitative RT-PCR data using Bio-Rad CFX Manager Software Version 1.6. Actin and tubulin were used as housekeeping genes for the calculation of relative transcript abundance. The sizes of the amplified products were confirmed via gel electrophoresis. Negative controls with no templates were carried out concurrently.Table 1
Forward (F) and reverse (R) primers used in real-time PCR.
Gene codeAnnotation (sequence homology)Primer sequence (5′→3′)Tm (°C) aCyclesProduct lengthCp88Membrane channel proteinF: ATGTGGACTTTTTGCTGGAGTTTTT594090R: GCAATACAAGCATTTACCCAGTCCTCp173Low temperature and salt responsive proteinF: CTGACGCTCTTTGGATGGCTACC5940172R: ACGAATCAACGGTCACAAACACGCp215Putative abscisic acid-induced proteinF: TGCTCTTTACTTGCTGGCTCGTCT5940193R: TCTCGCCCATTCTCTTCTATGTATTTCCp274Glutathione S-transferaseF: ATGTGCTAAATCTCATTCCAGTCTTCC594085R: GCCGTGATATCCTGCGAAACTAGCp3591,4-Alpha-glucan-maltohydrolaseF: TATTATCTGGCAAGGTTCCACTCGTT594086R: TGTGTTGTAGAACCCAGCAGTCAGCCp364Cold acclimation Protein WCOR413-like protein beta formF: TCCAATAAGCCGAATCAAAATACAGT594084R: GTCTGTCTTCATCGCCAAATACTCCCp3753-Oxoacyl-[acyl-carrier-protein] synthaseF: GTAACTGCGCCGAAACGAGAAAAAG594078R: CCCAAAGACAGAAACCAGCCCCCp440Putative multiple stress-responsive zinc-finger proteinF: GTTTCCCGATTTTGTGTGGCG5940183R: AGTTGTTGATGCAGAGAATTGGTCCTCp465CatalaseF: CATCGTTCGCTTCTCTACCGTTATT5940118R: TGTTTCCGACCAGATCAAAGTTACCCp197AGL9.2F: TGAGAATAAGATAAATCGGCAGGTGAC5940128R: TTTTCCTCTGCTGGAGAAGACAACGCp203LLP-B3 proteinF: GCGATGTGAAACTTGTCAGTAGCCC5940171R: AGTTTGGCACAATCAGGACGAGGTTCp268Caffeoyl-CoAO-methyl-transferase-like proteinF: AACACAAATGGAGAAGAGAAAACCC5940193R: CATCGGCAGAGGTCGTCATAAGACp297SRG1-like proteinF: TAAGAACACCGTCCCATCTCGCTAC5940143R: GAGTGCAGTCTCTCCAATTCATCAGCp328STYLOSA proteinF: GGGGAACCAACAAGAACAACAGAT5940161R: GCCCAGTGGCACCATTTAGAAGATCp330RAC-like G-protein Rac1F: ACCGCTTCTACCGCCTTTTCTCC5940166R: AGGTCTTTCCAACAGCACCATCTCCCp360Peroxisomal fatty acid betao-xidation multifunctional protein AIM1F: TATGCTGGTGTTTTGAAAGAGGCG5940193R: CTATGCCAGACCCCATCAAACCTCCp383Secondary cell wall-related glycosyltransferase family 47F: CAGGCGACACCCCCTCCTCTAAC5940150R: TCAAGGCATCAGAGTTACGCACAAAGCp423Putative polyketide synthaseF: GAGTTCAGACACTGTATGCCCTTGG5940165R: TTGCTTGTAAAATTCGTCCTTCGTTCp436Allergen-like protein BRSn20F: CTCTTATCGCTCTCTGCTTCCTTCC5940176R: GTTGCTGTTAACGTCTCTGCATTCCCp458MADS-box protein 9F: TGGTTGAGTATGATGTAGAGAAGCGAC594098R: ACCCATCATCTGAAGGTGTGCTATTCp82Dormancy related proteinF: TGGAATAGAAAGAGCAAATACAGCG5940172R: TACAAAAGGCTCAATGGCGTCCCp24No hitsF: GCAGTTTACATACTACGGGAGAGGCT5940106R: CGGCTTACGGAATCGTCATCACCp64Putative proteinF: CCAATCACTCTCCCTGAGGATGTAT5940159R: TTCACCGACTCCTTGTTCTTTTAGCInternal controlActinF: GTTATGGTTGGGATGGGACAGAAAG5940199R: GGGCTTCAGTAAGGAAACAGGAInternal controlTublinF: TAGTGACAAGACAGTAGGTGGAGGT5940139R: GTAGGTTCCAGTCCTCACTTCATCaAnnealing temperature.


## 3. Results

### 3.1. General Characteristics of the cDNA Library and the Ests

A cDNA library constructed fromC. praecox flowers at different stages of development was used as a source of ESTs. The titers of the primary library and amplified library were 1.4×106 and 1.0×1010 pfu/mL, respectively, with a recombinant rate of 96% for the original library. The sizes of the inserts ranged from 0.5 to 2.5 kb, and the average insert size was estimated to be 1.1 kb by PCR amplification of inserts from 50 randomly selected clones. These results indicate that our cDNA library was qualified.In total, 896 random cDNA clones were successfully sequenced to generate ESTs. Trimming of the short sequences (<100 bp), vector sequences, and poor-quality sequences resulted in 867 high-quality ESTs, constituting a total of 584,201 bases in theC. praecox sequence. The average read length of these ESTs was 673.8 bp. All 867 sequences were deposited in GenBank under accession numbers DW222667 to DW223533. The clustering of ESTs generated 94 contigs (containing 2 or more ESTs) and 385 singletons (containing only 1 EST), yielding 479 unigenes. The redundancy of the library was calculated as 44.8% [(1 − Number of Unigenes/Number of ESTs) × 100%]. Figure 2 shows the distribution of ESTs in unigenes after clustering. Forty-two contigs had more than 2 ESTs, with the largest one containing 77 ESTs.Figure 2
Distribution ofC. praecox ESTs among unigenes.


### 3.2. Functional Annotation and Classification ofC. praecox Unigenes

The 479 unigenes were compared with the nonredundant protein and nucleotide sequences database in NCBI using BLAST. Four hundred five unique sequences, corresponding to 84.6% of all the unigenes, shared significant homology with sequences in the public databases. Of these, 266 were similar to genes of known functions, whereas 139 were similar to putative uncharacterized proteins (Table2). The remaining 74 unigenes (15.4%) had only very weak or no matches and were considered as novel genes in C. praecox flowers.Table 2
Overview ofC. praecox flowers ESTs.
ItemsNumberTotal sequenced cDNA896Total number of ESTs analyzed867Total reading valid length (bp)584,201Average EST length (bp)673.8Unigenes479Contigs94Singletons385Redundancy44.8Table3 summarizes the highly expressed genes that contained more than 5 ESTs in one contig. The first and third most abundant ESTs were homologous to lipid transfer protein (LTP); the second most abundant ESTs encoded protein related to the adenine nucleotide translocator and then seed-specific protein, mannose-specific lectin, and LEA III protein (Table 3). These ESTs possibly corresponded to the most abundantly expressed genes in C. praecox flowers. Two hundred eighty-one ESTs were found to have more than five copies (Table 3), 385 had single copies of ESTs, and the remaining ESTs contained two to four copies in C. praecox flowers.Table 3
The 15 most abundant ESTs in theC. praecoxflowers cDNA library.
Gene codeLength(bp)Putative functionE-valueNumber of ESTsCp11319lipid transfer protein1.00E-3477Cp22124adenine nucleotide translocator9.00E-8047Cp41024lipid transfer protein2.00E-3537Cp61144hypothetical protein1.00E-3419Cp7710seed specific protein Bn15D18B4.00E-1918Cp20896mannose specific lectin5.00E-2817Cp16593putative LEA III protein isoform 26.00E-1010Cp10900glutathione S-transferase GST 224.00E-808Cp11928palmitoyl-acyl carrier protein thioesterase7.00E-228Cp141531hypothetical protein2.00E-798Cp17717stearoyl acyl carrier protein desaturase2.00E-317Cp261074aquaporin8.00E-977Cp23696S28 ribosomal protein4.00E-206Cp24574No hitsNULL6Cp27652proline-rich protein1.00E-306Total1488215281The database-matched 405 unigenes were found to have BLAST hits with 99 organisms, among which the highest number was fromArabidopsis thaliana (32.6%; 132 unigenes), followed by Oryza sativa (24.0%; 97 unigenes) and Nicotiana tabacum (4.2%; 17 unigenes) (Table 4). The remaining 67 unigenes (16.5%) had BLAST hits with a single organism only. The extensive distribution of the matched organisms may be attributed to the fact that the genome of C. praecox significantly differed from the genomes of model plants, as well as the fact that the relative plants’ genomes have not yet been widely studied.Table 4
Statistics of organism origin of matched-function homologs.
Organism nameNo. of unigenesPercentage (%)aArabidopsis thaliana13232.6%Oryza sativa (japonica cultivar-group)9724.0%Nicotiana tabacum174.2%Glycine max82.0%Lycopersicon esculentum71.7%Gossypium hirsutum61.5%Cucumis sativus51.2%Capsicum annuum51.2%Pisum sativum41.0%Solanum tuberosum41.0%Petunia x hybrida41.0%Nicotiana attenuata30.7%Helianthus annuus30.7%Beta vulgaris30.7%Calycanthus floridus var. glaucus30.7%Ipomoea batatas30.7%Hyacinthus orientalis30.7%Persea americana30.7%Ricinus communis30.7%Lily mottle virus30.7%Medicago sativa20.5%Zea mays20.5%Cicer arietinum20.5%Triticum aestivum20.5%Hevea brasiliensis20.5%Asparagus officinalis20.5%Panax ginseng20.5%Malus x domestica20.5%Elaeis oleifera20.5%Pinus taeda20.5%Spinacia oleracea20.5%The other 67 organisms116.5%aProportion from unigenes with BLAST hits.The initial annotations were further simplified into plant-specific annotations (plant GO slim;http://amigo.geneontology.org/) to obtain additional insights into the putative functions of unigenes. Of the 479 C. praecox unigenes, 364 were assigned GO terms in any category (biological, cellular, and molecular). Figures 3, 4, and 5 classify the unigenes according to terms in the biological process ontology, molecular function ontology, and cellular component ontology, respectively.Figure 3
GO classification of the ESTs based on their biological functions in theC. praecox flower cDNA library.Figure 4
GO classification of the ESTs based on their molecular functions in theC. praecox flower cDNA library.Figure 5
GO classification of the ESTs based on their cellular components in theC. praecox flower cDNA library.


### 3.3. Sequences Related to Stress and Defense

Among the ESTs that matched genes with known or putative functions, approximately 95 unigenes (291 ESTs) involved in stress and defense accounted for 19.8% (95/479) of all unigenes and 33.6% (291/867) of all ESTs. Table5 shows the nonredundant ESTs that share similarities with genes related to defense and stress response according to GO classifications and previously published data. As expected, the ESTs involved in stress and defense were highly abundant in the library because the C. praecox thrives and blossoms in winter, thereby confirming previous reports of related transcripts with higher levels of defense in developing flowers [16–20].Table 5
Sequences related to stress and defense.
Gene codePutative functionE-valueMatches no.Number of ESTsCp465catalase1.00E-148BAC79443.11Cp346transporter1.00E-138AAP13421.12Cp382peroxidase1.00E-125AAK52084.11Cp374GMPase1.00E-122AAT58365.11Cp333putative aldolase1.00E-122AAM64281.11Cp397chitinase1.00E-105AAF04454.11Cp452cysteine synthase3.00E-99BAA05965.11Cp314tryptophan synthase beta chain2.00E-94AAN15574.11Cp357Stromal cell-derived factor 2-like protein2.00E-93XP_481233.11Cp191glutathione S-transferase GST 235.00E-92AAG34813.12Cp25dehydroascorbate reductase2.00E-89AAL71857.13Cp438ras-related protein RAB8-39.00E-89BAB84324.11Cp459GTP-binding protein1.00E-86AAN31076.11Cp340putative 3-Hydroxyisobutyryl-coenzyme A hydrolase9.00E-84AAN41356.11Cp379cell-autonomous heat shock cognate protein 701.00E-82gAAN86276.11Cp39auxin-binding protein3.00E-80AAB51240.14Cp10glutathione S-transferase GST 224.00E-80AAG34812.18Cp304putative 2-Nitropropane dioxygenase2.00E-79AAL34288.13Cp21glutathione S-transferase GST 222.00E-78AAG34812.13Cp311GTP-binding protein6.00E-73AAM12880.11Cp353patatin-like protein1.00E-73CAB16788.12Cp274glutathione S-transferase1.00E-72AAF61392.11Cp3591,4-Alpha-glucan-maltohydrolase1.00E-71CAH60892.11Cp59BTF3b-like transcription factor2.00E-68AAT67244.13Cp76NDKB_FLABI Nucleoside diphosphate kinase B2.00E-63P479201Cp166putative tropinone reductase6.00E-59AAM62552.11Cp313putative glutathione S-transferase T35.00E-56AAG16758.11Cp156protolysis and peptidolysis6.00E-56CAA50022.11Cp120putative pyruvate dehydrogenase E1 alpha subunit8.00E-56XP_467697.11Cp80sadtomato protein7.00E-55AAS77347.11Cp440putative multiple stress-responsive zinc-finger protein6.00E-54BAD35553.11Cp175epoxide hydrolase6.00E-50AAB02006.11Cp364cold acclimation protein WCOR413-like protein beta form4.00E-49AAG13394.11Cp471pathogenesis-related protein 4B1.00E-48eCAA41438.11Cp145Thioredoxin4.00E-44CAA77847.11Cp406auxin-responsive family protein1.00E-43NP_565113.11Cp38basic PR-1 protein precursor3.00E-43AAU20808.14Cp215putative abscisic acid-induced protein4.00E-41BAD37454.11Cp66GLRX_RICCO Glutaredoxin5.00E-40sp|P55143|2Cp235histone H2A5.00E-38CAA07234.11Cp111topoisomerase 6 subunit B5.00E-37CAC24690.11Cp232cysteine proteinase inhibitor-like protein2.00E-36BAB03156.12Cp220glyoxalase I family protein2.00E-35NP_973579.12Cp214acyl-CoA-binding protein3.00E-35CAA70200.11Cp27proline-rich protein1.00E-30AAF78903.16Cp20lectin5.00E-28AAD45250.117Cp134metallothionein-like protein6.00E-27CAB52585.11Cp127Gip1-like protein2.00E-24CAD10106.11Cp296bZIP transcription factor2.00E-24AAN61914.12Cp3753-Oxoacyl-[acyl-carrier-protein] synthase2.00E-23AAC78479.11Cp104metallothionein-like protein type 24.00E-18CAB77242.12Cp372allyl alcohol dehydrogenase5.00E-17BAA89423.11Cp36metallothionein-like protein type 21.00E-13AAV97748.15Cp173low temperature and salt responsive protein1.00E-12BAC23051.11Cp472cystatin2.00E-07AAO18638.11Cp88membrane channel protein3.00E-06CAB83138.12Also related to developmentCp391actin1.00E-166AAN40685.11Cp402xyloglucan endotransglycosylase precursor1.00E-110AAC09388.11Cp26aquaporin8.00E-97CAE53881.17Cp3814-Coumarate-CoA ligase-like protein2.00E-89AAP03021.11Cp461glyceraldehyde-3-phosphate dehydrogenase1.00E-86CAC88118.11Cp2913-Ketoacyl-CoA thiolase6.00E-84CAA47926.11Cp378probable ubiquitin-conjugating enzyme E22.00E-83AAC32141.11Cp400auxin response factor 45.00E-81BAD19064.11Cp288expansin4.00E-80AAO15998.12Cp226calmodulin3.00E-79AAB68399.14Cp323putative calmodulin3.00E-73CAC84563.12Cp298putative adrenal gland protein AD-0042.00E-65XP_477670.11Cp326UDP-glucose:salicylic acid glucosyltransferase2.00E-65AAF61647.11Cp58actin-depolymerizing factor 11.00E-63AAK72617.12Cp167ABC transporter4.00E-60AAP80385.11Cp286putative CCR4-associated factor4.00E-60AAN13153.11Cp368MYB8 protein1.00E-60CAD87008.11Cp324actin depolymerizing factor2.00E-59AAD23407.11Cp450TAF91.00E-58AAR28026.11Cp405myb family transcription factor-like2.00E-52BAD29385.11Cp77putative actin-depolymerizing factor 14.00E-48XP_475079.11Cp121TAF103.00E-41AAR28030.11Cp347putative transcription factor2.00E-38XP_479103.112Cp4lipid-transfer protein2.00E-35AAS13435.137Cp1lipid transfer protein1.00E-34CAA63340.177Cp275lipid transfer protein precursor4.00E-30AAL27855.12Cp9lipid transfer protein 1 precursor3.00E-25AAT45202.11Cp216lipid-transfer protein3.00E-24NP_915960.14Cp90zinc finger homeobox family protein3.00E-22NP_565088.11Cp41lipid transfer protein3.00E-18BAC77694.11Cp13lipid transfer protein isoform 41.00E-16AAO33394.11Cp338polyprotein2.00E-16CAD92110.11Cp82dormancy related protein, putative6.00E-15AAG51119.11Cp8lipid-transfer protein7.00E-14AAS13435.11Cp141lipid-transfer protein4.00E-13AAS13435.11Cp211neutral/alkaline invertase 12.00E-12AAV28809.11Cp451putative MYB family transcription factor8.00E-11AAD23043.11Cp34Lea5 protein3.00E-10CAA86851.14Cp16putative LEA III protein isoform 26.00E-10CAC39110.110Twelve unigenes (22 ESTs) related to cold stress tolerance were found in the library; these were classified as “response to cold” according to GO terms [e.g.,β-amylase, acyl-CoA-binding protein, 3-hydroxyisobutyryl-coenzyme A hydrolase, catalase, low-temperature and salt-responsive protein, glutathione S-transferase (GST), membrane channel protein, and abscisic-acid-induced protein] and “cold acclimation” (e.g., fatty acid biosynthesis 1 and WCOR413-like protein). In this class, the most abundant sequences encode GST. Five GSTs encoded by 16 ESTs were identified in the library. GST genes exhibited a diverse range of responses to jasmonates, salicylic acid, ethylene, as well as oxidative stress in Arabidopsis [21], and were induced by heavy metals and hypoxic stress in rice roots [22]. In the current study, however, not enough cold-related unigenes were obtained, as expected. Some unknown functional genes related to cold stress tolerance likely exist in this library.Another class of genes involved in “response to absence of light,” “low light intensity,” and “response to red or blue light” according to GO terms was also represented in the library. Eight unigenes (10 ESTs) were identified, including acyl-CoA-binding protein, sadtomato protein, and catalase, among others.Of the stress- and defense-related unigenes, 39 were possibly related to development. Nine types of LTPs or LTP precursors were encoded by the highest number of ESTs in this study (Table5). Research has suggested different functions for LTPs in the physiology of plants, such as cutin synthesis, β-oxidation, somatic embryogenesis, pollen development, allergenics, plant signaling, and plant defense [23–25], but the true physiological role of LTPs in C. praecox flowers has yet to be determined. Fourteen ESTs encoding two kinds of LEA proteins, groups III and V, were identified. The presence of LEA proteins correlates well with freezing, water deficit, and salt stress [26–28], probably through the prevention of enzyme aggregation [29], and likely plays a similar role in C. praecox. Some other development-related unigenes were also found, such as MYB, calmodulin, actin, CCR4-associated factor, ubiquitin-conjugating enzyme E2, and ABC transporter, which are involved in transcription factor activity, signal transduction, cell structure, nucleotide metabolism, protein metabolic process, and transporter activity, respectively.


### 3.4. Sequences Related to Floral Development

Table6 shows the 19 unigenes related to floral development. Five of these were homologous to MADS box transcription factor genes. Plant life critically depends on the function of MADS box genes encoding MADS domain transcription factors, which are present to a limited extent in nearly all major eukaryotic groups but constitute a large gene family in land plants [30].Table 6
Sequences related to floral development.
Gene codePutative functionE-valueMatches no.Number of ESTsCp360Peroxisomal fatty acid beta-oxidation multifunctional protein AIM11.00E-108XP_464920.11Cp268caffeoyl-CoAO-methyl-transferase-like protein1.00E-107CAB80122.12Cp237MADS box protein1.00E-106AAQ83835.12Cp116MADS box transcription factor AP3-12.00E-95AAF73928.11Cp330RAC-like G-protein Rac11.00E-93AAD47828.21Cp217beta-D-glucosidase2.00E-92AAQ17461.11Cp413leucine-rich repeat transmembrane protein kinase6.00E-86NP_192248.21Cp408isoflavone reductase related protein3.00E-85AAC24001.11Cp423putative polyketide synthase8.00E-76NP_919053.11Cp335AGL 63.00E-73AAY25580.12Cp395isopentenyl pyrophosphate:dimethyllallyl pyrophosphate isomerase2.00E-70BAB09611.11Cp297SRG1-like protein1.00E-61CAB81342.13Cp383secondary cell wall-related glycosyltransferase family 473.00E-49AAX33321.11Cp197AGL9.24.00E-42AAX15924.11Cp436allergen-like protein BRSn201.00E-37AAM62935.11Cp407caffeateO-methyltransferase8.00E-23AAV36331.11Cp203LLP-B3 protein1.00E-22AAN76546.11Cp328STYLOSA protein1.00E-15CAF18245.11Cp458MADS-box protein 99.00E-15AAQ72497.11MADS box genes control diverse developmental processes in flowering plants—ranging from root to flower and fruit development—and they especially control the processes of the transition from vegetative to reproductive development and establishment of floral organ identity [31]. The present study has also identified six unigenes related to secondary metabolism that are probably involved in floral color and fragrance.The role of chilling temperature in dormancy in vegetative buds and induction of flowering has been investigated in many temperate-region species, particularly inA. thaliana and Populus spp. [32, 33]. The physiological processes of dormancy release and induction of flowering competence rely on longer-term chilling temperature and a period of vernalization, respectively [34, 35]. The current study has identified dormancy and vernalization-related genes; however, only one unigene (Cp82; E-value = 6.00E-15) was annotated as a dormancy-related protein in the library The possible reason was only a small-scale sequencing in our study or the processes of dormancy release and induction of flowering competence in C. praecox only last a short period.


### 3.5. Expression Analysis of Cold-Responsive and Floral Development-Related Genes

Real-time RT-PCR was performed for 12 unigenes, including 9 selected from the GO Slim annotation belonging to “response to cold” (Cp88, Cp173, Cp215, Cp274, Cp359, Cp364, Cp375, Cp440, and Cp465), 1 annotated as a dormancy-related protein (Cp82), and 2 without functional annotation (Cp24 and Cp64), to analyze the changes in their expression due to cold stress (Figure6). The data revealed that Cp82 responded to cold stress immediately after treatment and reached peak expression levels as early as 15 min. Cp24 (no annotation), Cp88 (membrane channel protein), Cp274 (GST), Cp364 (cold acclimation protein WCOR413-like protein), and Cp440 (multiple stress-responsive zinc-finger protein) were upregulated after 15 min of treatment and reached their expression peaks at 1 h. Cp173 (low-temperature and salt-responsive protein), Cp215 (abscisic acid-induced protein), Cp359 (1,4-alpha-glucan-maltohydrolase), Cp375 (3-oxoacyl-[acyl-carrier-protein] synthase), and Cp465 (catalase) displayed later responses and reached their peak expression levels at 6 h. Cp64 (no annotation) exhibited minor changes in its transcript level.Figure 6
Expression analysis of 12 genes under cold stress (4°C).Real-time RT-PCR was also applied to validate the expression patterns of 11 genes related to floral development (Cp197, Cp203, Cp268, Cp297, Cp328, Cp330, Cp360, Cp383, Cp423, Cp436, and Cp458), Cp24, Cp64, and Cp82 (Figure7). The expression patterns of these 14 genes were analyzed in roots, stems, leaves, outer tepals, middle tepals, inner tepals, stamina, and pistils. The results showed clear differences in their expression. The genes related to floral development, except for Cp203 (LLP-B3 protein) and Cp297 (SRG1-like protein), increased by more than twofold in middle tepals. Cp24 and Cp203 presented an active expression in stamina but showed a very slight accumulation in other tissues. Cp268 (caffeoyl-CoA O-methyl-transferase-like protein) and Cp458 (MADS box protein 9) presented a peak in middle tepals but were only slightly or not expressed at all in roots, stems, leaves, and pistils. Cp197 (AGL9.2) was higher in all reproductive organs and presented its peak in middle tepals but was not detected in roots, stems, and leaves. The expression of Cp328 (STYLOSA protein) and that of Cp383 (secondary cell wall-related glycosyltransferase family 47) were not detected in roots; moreover, Cp328 increased by more than fourfold in middle tepals, whereas Cp383 did in pistils and stems. The accumulation of Cp82 (dormancy-related protein) transcripts was higher in stems, leaves, stamina, and pistils. The other genes expressed in all the detected tissues, but Cp64 and Cp297 (SRG1-like protein) presented lower levels of expression. Cp360 (peroxisomal fatty acid β-oxidation multifunctional protein AIM1) increased by more than twofold in roots, middle tepals, and pistils.Figure 7
Expression analysis of 14 genes in different tissues ofC. praecox.Quantitative real-time PCR methods were used to validate the transcript levels of the 14 genes further during the blooming process inC. praecox (Figure 8). The results showed that all these genes were not detected or very slightly expressed in Stage 1 and that Cp64, Cp82, Cp268, Cp297, Cp330, Cp360, and Cp383 had almost no accumulation in Stage 2. The transcript accumulations of Cp203 and Cp458 were sharply elevated to the highest level in Stage 2 but dramatically decreased at the subsequent stages of floral development. The expressions of Cp24, Cp82, and Cp436 presented a peak in Stage 4. The transcript accumulations of the other 9 genes increased during the six developmental stages and reached their peak in Stage 6. The expression of Cp297(SRG1-like protein) was significantly high in Stage 6 but very low in the other stages, and it was associated with flower senescence. The SRG1 gene is reportedly expressed in senescing organs of Arabidopsis plants [36].Figure 8
Expression analysis of 14 genes duringC. praecox flowering.The current study found few references about bud dormancy inC. praecox. The expression pattern of Cp82 was very attractive. However, the results indicate that Cp82 is not certainly related to flower-bud dormancy. The dormancy-related protein CAA93825.1, which matched Cp82 (E-value = 1e−15), has been reported to play a role in dormancy breaking and in the germination of Trollius ledebourii seeds [37]. These data warrant further research into the relativity of Cp82 to dormancy breaking and germination in C. praecox seeds.


## 3.1. General Characteristics of the cDNA Library and the Ests

A cDNA library constructed fromC. praecox flowers at different stages of development was used as a source of ESTs. The titers of the primary library and amplified library were 1.4×106 and 1.0×1010 pfu/mL, respectively, with a recombinant rate of 96% for the original library. The sizes of the inserts ranged from 0.5 to 2.5 kb, and the average insert size was estimated to be 1.1 kb by PCR amplification of inserts from 50 randomly selected clones. These results indicate that our cDNA library was qualified.In total, 896 random cDNA clones were successfully sequenced to generate ESTs. Trimming of the short sequences (<100 bp), vector sequences, and poor-quality sequences resulted in 867 high-quality ESTs, constituting a total of 584,201 bases in theC. praecox sequence. The average read length of these ESTs was 673.8 bp. All 867 sequences were deposited in GenBank under accession numbers DW222667 to DW223533. The clustering of ESTs generated 94 contigs (containing 2 or more ESTs) and 385 singletons (containing only 1 EST), yielding 479 unigenes. The redundancy of the library was calculated as 44.8% [(1 − Number of Unigenes/Number of ESTs) × 100%]. Figure 2 shows the distribution of ESTs in unigenes after clustering. Forty-two contigs had more than 2 ESTs, with the largest one containing 77 ESTs.Figure 2
Distribution ofC. praecox ESTs among unigenes.


## 3.2. Functional Annotation and Classification ofC. praecox Unigenes

The 479 unigenes were compared with the nonredundant protein and nucleotide sequences database in NCBI using BLAST. Four hundred five unique sequences, corresponding to 84.6% of all the unigenes, shared significant homology with sequences in the public databases. Of these, 266 were similar to genes of known functions, whereas 139 were similar to putative uncharacterized proteins (Table2). The remaining 74 unigenes (15.4%) had only very weak or no matches and were considered as novel genes in C. praecox flowers.Table 2
Overview ofC. praecox flowers ESTs.
ItemsNumberTotal sequenced cDNA896Total number of ESTs analyzed867Total reading valid length (bp)584,201Average EST length (bp)673.8Unigenes479Contigs94Singletons385Redundancy44.8Table3 summarizes the highly expressed genes that contained more than 5 ESTs in one contig. The first and third most abundant ESTs were homologous to lipid transfer protein (LTP); the second most abundant ESTs encoded protein related to the adenine nucleotide translocator and then seed-specific protein, mannose-specific lectin, and LEA III protein (Table 3). These ESTs possibly corresponded to the most abundantly expressed genes in C. praecox flowers. Two hundred eighty-one ESTs were found to have more than five copies (Table 3), 385 had single copies of ESTs, and the remaining ESTs contained two to four copies in C. praecox flowers.Table 3
The 15 most abundant ESTs in theC. praecoxflowers cDNA library.
Gene codeLength(bp)Putative functionE-valueNumber of ESTsCp11319lipid transfer protein1.00E-3477Cp22124adenine nucleotide translocator9.00E-8047Cp41024lipid transfer protein2.00E-3537Cp61144hypothetical protein1.00E-3419Cp7710seed specific protein Bn15D18B4.00E-1918Cp20896mannose specific lectin5.00E-2817Cp16593putative LEA III protein isoform 26.00E-1010Cp10900glutathione S-transferase GST 224.00E-808Cp11928palmitoyl-acyl carrier protein thioesterase7.00E-228Cp141531hypothetical protein2.00E-798Cp17717stearoyl acyl carrier protein desaturase2.00E-317Cp261074aquaporin8.00E-977Cp23696S28 ribosomal protein4.00E-206Cp24574No hitsNULL6Cp27652proline-rich protein1.00E-306Total1488215281The database-matched 405 unigenes were found to have BLAST hits with 99 organisms, among which the highest number was fromArabidopsis thaliana (32.6%; 132 unigenes), followed by Oryza sativa (24.0%; 97 unigenes) and Nicotiana tabacum (4.2%; 17 unigenes) (Table 4). The remaining 67 unigenes (16.5%) had BLAST hits with a single organism only. The extensive distribution of the matched organisms may be attributed to the fact that the genome of C. praecox significantly differed from the genomes of model plants, as well as the fact that the relative plants’ genomes have not yet been widely studied.Table 4
Statistics of organism origin of matched-function homologs.
Organism nameNo. of unigenesPercentage (%)aArabidopsis thaliana13232.6%Oryza sativa (japonica cultivar-group)9724.0%Nicotiana tabacum174.2%Glycine max82.0%Lycopersicon esculentum71.7%Gossypium hirsutum61.5%Cucumis sativus51.2%Capsicum annuum51.2%Pisum sativum41.0%Solanum tuberosum41.0%Petunia x hybrida41.0%Nicotiana attenuata30.7%Helianthus annuus30.7%Beta vulgaris30.7%Calycanthus floridus var. glaucus30.7%Ipomoea batatas30.7%Hyacinthus orientalis30.7%Persea americana30.7%Ricinus communis30.7%Lily mottle virus30.7%Medicago sativa20.5%Zea mays20.5%Cicer arietinum20.5%Triticum aestivum20.5%Hevea brasiliensis20.5%Asparagus officinalis20.5%Panax ginseng20.5%Malus x domestica20.5%Elaeis oleifera20.5%Pinus taeda20.5%Spinacia oleracea20.5%The other 67 organisms116.5%aProportion from unigenes with BLAST hits.The initial annotations were further simplified into plant-specific annotations (plant GO slim;http://amigo.geneontology.org/) to obtain additional insights into the putative functions of unigenes. Of the 479 C. praecox unigenes, 364 were assigned GO terms in any category (biological, cellular, and molecular). Figures 3, 4, and 5 classify the unigenes according to terms in the biological process ontology, molecular function ontology, and cellular component ontology, respectively.Figure 3
GO classification of the ESTs based on their biological functions in theC. praecox flower cDNA library.Figure 4
GO classification of the ESTs based on their molecular functions in theC. praecox flower cDNA library.Figure 5
GO classification of the ESTs based on their cellular components in theC. praecox flower cDNA library.


## 3.3. Sequences Related to Stress and Defense

Among the ESTs that matched genes with known or putative functions, approximately 95 unigenes (291 ESTs) involved in stress and defense accounted for 19.8% (95/479) of all unigenes and 33.6% (291/867) of all ESTs. Table5 shows the nonredundant ESTs that share similarities with genes related to defense and stress response according to GO classifications and previously published data. As expected, the ESTs involved in stress and defense were highly abundant in the library because the C. praecox thrives and blossoms in winter, thereby confirming previous reports of related transcripts with higher levels of defense in developing flowers [16–20].Table 5
Sequences related to stress and defense.
Gene codePutative functionE-valueMatches no.Number of ESTsCp465catalase1.00E-148BAC79443.11Cp346transporter1.00E-138AAP13421.12Cp382peroxidase1.00E-125AAK52084.11Cp374GMPase1.00E-122AAT58365.11Cp333putative aldolase1.00E-122AAM64281.11Cp397chitinase1.00E-105AAF04454.11Cp452cysteine synthase3.00E-99BAA05965.11Cp314tryptophan synthase beta chain2.00E-94AAN15574.11Cp357Stromal cell-derived factor 2-like protein2.00E-93XP_481233.11Cp191glutathione S-transferase GST 235.00E-92AAG34813.12Cp25dehydroascorbate reductase2.00E-89AAL71857.13Cp438ras-related protein RAB8-39.00E-89BAB84324.11Cp459GTP-binding protein1.00E-86AAN31076.11Cp340putative 3-Hydroxyisobutyryl-coenzyme A hydrolase9.00E-84AAN41356.11Cp379cell-autonomous heat shock cognate protein 701.00E-82gAAN86276.11Cp39auxin-binding protein3.00E-80AAB51240.14Cp10glutathione S-transferase GST 224.00E-80AAG34812.18Cp304putative 2-Nitropropane dioxygenase2.00E-79AAL34288.13Cp21glutathione S-transferase GST 222.00E-78AAG34812.13Cp311GTP-binding protein6.00E-73AAM12880.11Cp353patatin-like protein1.00E-73CAB16788.12Cp274glutathione S-transferase1.00E-72AAF61392.11Cp3591,4-Alpha-glucan-maltohydrolase1.00E-71CAH60892.11Cp59BTF3b-like transcription factor2.00E-68AAT67244.13Cp76NDKB_FLABI Nucleoside diphosphate kinase B2.00E-63P479201Cp166putative tropinone reductase6.00E-59AAM62552.11Cp313putative glutathione S-transferase T35.00E-56AAG16758.11Cp156protolysis and peptidolysis6.00E-56CAA50022.11Cp120putative pyruvate dehydrogenase E1 alpha subunit8.00E-56XP_467697.11Cp80sadtomato protein7.00E-55AAS77347.11Cp440putative multiple stress-responsive zinc-finger protein6.00E-54BAD35553.11Cp175epoxide hydrolase6.00E-50AAB02006.11Cp364cold acclimation protein WCOR413-like protein beta form4.00E-49AAG13394.11Cp471pathogenesis-related protein 4B1.00E-48eCAA41438.11Cp145Thioredoxin4.00E-44CAA77847.11Cp406auxin-responsive family protein1.00E-43NP_565113.11Cp38basic PR-1 protein precursor3.00E-43AAU20808.14Cp215putative abscisic acid-induced protein4.00E-41BAD37454.11Cp66GLRX_RICCO Glutaredoxin5.00E-40sp|P55143|2Cp235histone H2A5.00E-38CAA07234.11Cp111topoisomerase 6 subunit B5.00E-37CAC24690.11Cp232cysteine proteinase inhibitor-like protein2.00E-36BAB03156.12Cp220glyoxalase I family protein2.00E-35NP_973579.12Cp214acyl-CoA-binding protein3.00E-35CAA70200.11Cp27proline-rich protein1.00E-30AAF78903.16Cp20lectin5.00E-28AAD45250.117Cp134metallothionein-like protein6.00E-27CAB52585.11Cp127Gip1-like protein2.00E-24CAD10106.11Cp296bZIP transcription factor2.00E-24AAN61914.12Cp3753-Oxoacyl-[acyl-carrier-protein] synthase2.00E-23AAC78479.11Cp104metallothionein-like protein type 24.00E-18CAB77242.12Cp372allyl alcohol dehydrogenase5.00E-17BAA89423.11Cp36metallothionein-like protein type 21.00E-13AAV97748.15Cp173low temperature and salt responsive protein1.00E-12BAC23051.11Cp472cystatin2.00E-07AAO18638.11Cp88membrane channel protein3.00E-06CAB83138.12Also related to developmentCp391actin1.00E-166AAN40685.11Cp402xyloglucan endotransglycosylase precursor1.00E-110AAC09388.11Cp26aquaporin8.00E-97CAE53881.17Cp3814-Coumarate-CoA ligase-like protein2.00E-89AAP03021.11Cp461glyceraldehyde-3-phosphate dehydrogenase1.00E-86CAC88118.11Cp2913-Ketoacyl-CoA thiolase6.00E-84CAA47926.11Cp378probable ubiquitin-conjugating enzyme E22.00E-83AAC32141.11Cp400auxin response factor 45.00E-81BAD19064.11Cp288expansin4.00E-80AAO15998.12Cp226calmodulin3.00E-79AAB68399.14Cp323putative calmodulin3.00E-73CAC84563.12Cp298putative adrenal gland protein AD-0042.00E-65XP_477670.11Cp326UDP-glucose:salicylic acid glucosyltransferase2.00E-65AAF61647.11Cp58actin-depolymerizing factor 11.00E-63AAK72617.12Cp167ABC transporter4.00E-60AAP80385.11Cp286putative CCR4-associated factor4.00E-60AAN13153.11Cp368MYB8 protein1.00E-60CAD87008.11Cp324actin depolymerizing factor2.00E-59AAD23407.11Cp450TAF91.00E-58AAR28026.11Cp405myb family transcription factor-like2.00E-52BAD29385.11Cp77putative actin-depolymerizing factor 14.00E-48XP_475079.11Cp121TAF103.00E-41AAR28030.11Cp347putative transcription factor2.00E-38XP_479103.112Cp4lipid-transfer protein2.00E-35AAS13435.137Cp1lipid transfer protein1.00E-34CAA63340.177Cp275lipid transfer protein precursor4.00E-30AAL27855.12Cp9lipid transfer protein 1 precursor3.00E-25AAT45202.11Cp216lipid-transfer protein3.00E-24NP_915960.14Cp90zinc finger homeobox family protein3.00E-22NP_565088.11Cp41lipid transfer protein3.00E-18BAC77694.11Cp13lipid transfer protein isoform 41.00E-16AAO33394.11Cp338polyprotein2.00E-16CAD92110.11Cp82dormancy related protein, putative6.00E-15AAG51119.11Cp8lipid-transfer protein7.00E-14AAS13435.11Cp141lipid-transfer protein4.00E-13AAS13435.11Cp211neutral/alkaline invertase 12.00E-12AAV28809.11Cp451putative MYB family transcription factor8.00E-11AAD23043.11Cp34Lea5 protein3.00E-10CAA86851.14Cp16putative LEA III protein isoform 26.00E-10CAC39110.110Twelve unigenes (22 ESTs) related to cold stress tolerance were found in the library; these were classified as “response to cold” according to GO terms [e.g.,β-amylase, acyl-CoA-binding protein, 3-hydroxyisobutyryl-coenzyme A hydrolase, catalase, low-temperature and salt-responsive protein, glutathione S-transferase (GST), membrane channel protein, and abscisic-acid-induced protein] and “cold acclimation” (e.g., fatty acid biosynthesis 1 and WCOR413-like protein). In this class, the most abundant sequences encode GST. Five GSTs encoded by 16 ESTs were identified in the library. GST genes exhibited a diverse range of responses to jasmonates, salicylic acid, ethylene, as well as oxidative stress in Arabidopsis [21], and were induced by heavy metals and hypoxic stress in rice roots [22]. In the current study, however, not enough cold-related unigenes were obtained, as expected. Some unknown functional genes related to cold stress tolerance likely exist in this library.Another class of genes involved in “response to absence of light,” “low light intensity,” and “response to red or blue light” according to GO terms was also represented in the library. Eight unigenes (10 ESTs) were identified, including acyl-CoA-binding protein, sadtomato protein, and catalase, among others.Of the stress- and defense-related unigenes, 39 were possibly related to development. Nine types of LTPs or LTP precursors were encoded by the highest number of ESTs in this study (Table5). Research has suggested different functions for LTPs in the physiology of plants, such as cutin synthesis, β-oxidation, somatic embryogenesis, pollen development, allergenics, plant signaling, and plant defense [23–25], but the true physiological role of LTPs in C. praecox flowers has yet to be determined. Fourteen ESTs encoding two kinds of LEA proteins, groups III and V, were identified. The presence of LEA proteins correlates well with freezing, water deficit, and salt stress [26–28], probably through the prevention of enzyme aggregation [29], and likely plays a similar role in C. praecox. Some other development-related unigenes were also found, such as MYB, calmodulin, actin, CCR4-associated factor, ubiquitin-conjugating enzyme E2, and ABC transporter, which are involved in transcription factor activity, signal transduction, cell structure, nucleotide metabolism, protein metabolic process, and transporter activity, respectively.


## 3.4. Sequences Related to Floral Development

Table6 shows the 19 unigenes related to floral development. Five of these were homologous to MADS box transcription factor genes. Plant life critically depends on the function of MADS box genes encoding MADS domain transcription factors, which are present to a limited extent in nearly all major eukaryotic groups but constitute a large gene family in land plants [30].Table 6
Sequences related to floral development.
Gene codePutative functionE-valueMatches no.Number of ESTsCp360Peroxisomal fatty acid beta-oxidation multifunctional protein AIM11.00E-108XP_464920.11Cp268caffeoyl-CoAO-methyl-transferase-like protein1.00E-107CAB80122.12Cp237MADS box protein1.00E-106AAQ83835.12Cp116MADS box transcription factor AP3-12.00E-95AAF73928.11Cp330RAC-like G-protein Rac11.00E-93AAD47828.21Cp217beta-D-glucosidase2.00E-92AAQ17461.11Cp413leucine-rich repeat transmembrane protein kinase6.00E-86NP_192248.21Cp408isoflavone reductase related protein3.00E-85AAC24001.11Cp423putative polyketide synthase8.00E-76NP_919053.11Cp335AGL 63.00E-73AAY25580.12Cp395isopentenyl pyrophosphate:dimethyllallyl pyrophosphate isomerase2.00E-70BAB09611.11Cp297SRG1-like protein1.00E-61CAB81342.13Cp383secondary cell wall-related glycosyltransferase family 473.00E-49AAX33321.11Cp197AGL9.24.00E-42AAX15924.11Cp436allergen-like protein BRSn201.00E-37AAM62935.11Cp407caffeateO-methyltransferase8.00E-23AAV36331.11Cp203LLP-B3 protein1.00E-22AAN76546.11Cp328STYLOSA protein1.00E-15CAF18245.11Cp458MADS-box protein 99.00E-15AAQ72497.11MADS box genes control diverse developmental processes in flowering plants—ranging from root to flower and fruit development—and they especially control the processes of the transition from vegetative to reproductive development and establishment of floral organ identity [31]. The present study has also identified six unigenes related to secondary metabolism that are probably involved in floral color and fragrance.The role of chilling temperature in dormancy in vegetative buds and induction of flowering has been investigated in many temperate-region species, particularly inA. thaliana and Populus spp. [32, 33]. The physiological processes of dormancy release and induction of flowering competence rely on longer-term chilling temperature and a period of vernalization, respectively [34, 35]. The current study has identified dormancy and vernalization-related genes; however, only one unigene (Cp82; E-value = 6.00E-15) was annotated as a dormancy-related protein in the library The possible reason was only a small-scale sequencing in our study or the processes of dormancy release and induction of flowering competence in C. praecox only last a short period.


## 3.5. Expression Analysis of Cold-Responsive and Floral Development-Related Genes

Real-time RT-PCR was performed for 12 unigenes, including 9 selected from the GO Slim annotation belonging to “response to cold” (Cp88, Cp173, Cp215, Cp274, Cp359, Cp364, Cp375, Cp440, and Cp465), 1 annotated as a dormancy-related protein (Cp82), and 2 without functional annotation (Cp24 and Cp64), to analyze the changes in their expression due to cold stress (Figure6). The data revealed that Cp82 responded to cold stress immediately after treatment and reached peak expression levels as early as 15 min. Cp24 (no annotation), Cp88 (membrane channel protein), Cp274 (GST), Cp364 (cold acclimation protein WCOR413-like protein), and Cp440 (multiple stress-responsive zinc-finger protein) were upregulated after 15 min of treatment and reached their expression peaks at 1 h. Cp173 (low-temperature and salt-responsive protein), Cp215 (abscisic acid-induced protein), Cp359 (1,4-alpha-glucan-maltohydrolase), Cp375 (3-oxoacyl-[acyl-carrier-protein] synthase), and Cp465 (catalase) displayed later responses and reached their peak expression levels at 6 h. Cp64 (no annotation) exhibited minor changes in its transcript level.Figure 6
Expression analysis of 12 genes under cold stress (4°C).Real-time RT-PCR was also applied to validate the expression patterns of 11 genes related to floral development (Cp197, Cp203, Cp268, Cp297, Cp328, Cp330, Cp360, Cp383, Cp423, Cp436, and Cp458), Cp24, Cp64, and Cp82 (Figure7). The expression patterns of these 14 genes were analyzed in roots, stems, leaves, outer tepals, middle tepals, inner tepals, stamina, and pistils. The results showed clear differences in their expression. The genes related to floral development, except for Cp203 (LLP-B3 protein) and Cp297 (SRG1-like protein), increased by more than twofold in middle tepals. Cp24 and Cp203 presented an active expression in stamina but showed a very slight accumulation in other tissues. Cp268 (caffeoyl-CoA O-methyl-transferase-like protein) and Cp458 (MADS box protein 9) presented a peak in middle tepals but were only slightly or not expressed at all in roots, stems, leaves, and pistils. Cp197 (AGL9.2) was higher in all reproductive organs and presented its peak in middle tepals but was not detected in roots, stems, and leaves. The expression of Cp328 (STYLOSA protein) and that of Cp383 (secondary cell wall-related glycosyltransferase family 47) were not detected in roots; moreover, Cp328 increased by more than fourfold in middle tepals, whereas Cp383 did in pistils and stems. The accumulation of Cp82 (dormancy-related protein) transcripts was higher in stems, leaves, stamina, and pistils. The other genes expressed in all the detected tissues, but Cp64 and Cp297 (SRG1-like protein) presented lower levels of expression. Cp360 (peroxisomal fatty acid β-oxidation multifunctional protein AIM1) increased by more than twofold in roots, middle tepals, and pistils.Figure 7
Expression analysis of 14 genes in different tissues ofC. praecox.Quantitative real-time PCR methods were used to validate the transcript levels of the 14 genes further during the blooming process inC. praecox (Figure 8). The results showed that all these genes were not detected or very slightly expressed in Stage 1 and that Cp64, Cp82, Cp268, Cp297, Cp330, Cp360, and Cp383 had almost no accumulation in Stage 2. The transcript accumulations of Cp203 and Cp458 were sharply elevated to the highest level in Stage 2 but dramatically decreased at the subsequent stages of floral development. The expressions of Cp24, Cp82, and Cp436 presented a peak in Stage 4. The transcript accumulations of the other 9 genes increased during the six developmental stages and reached their peak in Stage 6. The expression of Cp297(SRG1-like protein) was significantly high in Stage 6 but very low in the other stages, and it was associated with flower senescence. The SRG1 gene is reportedly expressed in senescing organs of Arabidopsis plants [36].Figure 8
Expression analysis of 14 genes duringC. praecox flowering.The current study found few references about bud dormancy inC. praecox. The expression pattern of Cp82 was very attractive. However, the results indicate that Cp82 is not certainly related to flower-bud dormancy. The dormancy-related protein CAA93825.1, which matched Cp82 (E-value = 1e−15), has been reported to play a role in dormancy breaking and in the germination of Trollius ledebourii seeds [37]. These data warrant further research into the relativity of Cp82 to dormancy breaking and germination in C. praecox seeds.


## 4. Conclusions

A cDNA library was constructed to generate an EST collection fromC. praecox, thereby providing a preliminary view into the genomic properties of this species. This collection of high-quality ESTs represents the first EST data set for C. praecox. Eight hundred sixty-seven valid EST sequences were generated, and 479 unigenes were assembled, among which 266 unigenes (55.53%) were identified according to their significant similarities with proteins of known functions. The EST sequences have been deposited in GenBank under accession numbers DW222667 to DW223533. Stress response genes and floral development-related genes were also identified. This study evaluated the expression patterns of 23 genes, including 2 novel ones, using real-time RT-PCR. Further investigations in this direction would help in the discovery of promising candidates with a key role in the development of stress tolerance for woody plants to bloom.



---

*Source: 134596-2012-03-21.xml*