Skip to main content
Download PDF

Gene expression profile indicates involvement of NO in Camellia sinensis pollen tube growth at low temperature

BMC Genomics volume 17, Article number: 809 (2016) Cite this article

Abstract

Background

Nitric oxide (NO) functions as a critical signaling molecule in the low-temperature stress responses in plants, including polarized pollen tube growth in Camellia sinensis. Despite this, the potential mechanisms underlying the participation of NO in pollen tube responses to low temperature remain unclear. Here, we investigate alterations to gene expression in C. sinensis pollen tubes exposed to low-temperature stress and NO using RNA-Seq technology, in order to find the potential candidate genes related to the regulation of pollen tube elongation by NO under low-temperature stress.

Results

Three libraries were generated from C. sinensis cv. ‘Longjingchangye’ pollen tubes cultured at 25 °C (CsPT-CK) and 4 °C (CsPT-LT) or with 25 μM DEA NONOate (CsPT-NO). The number of unigenes found for the three biological replications were 39,726, 40,440 and 41,626 for CsPT-CK; 36,993, 39,070 and 39,439 for CsPT-LT; and 39,514, 38,298 and 39,061 for CsPT-NO. A total of 36,097 unique assembled and annotated sequences from C. sinensis pollen tube reads were found in a BLAST search of the following databases: NCBI non-redundant nucleotide, Swiss-prot protein, Kyoto Encyclopedia of Genes and Genomes, Cluster of Orthologous Groups of proteins, and Gene Ontology. The absolute values of log2Ratio > 1 and probability > 0.7 were used as the thresholds for significantly differential gene expression, and 766, 497 and 929 differentially expressed genes (DEGs) were found from the comparison analyses of the CK-VS-LT, CK-VS-NO and LT-VS-NO libraries, respectively. Genes related to metabolism and signaling pathways of plant hormones, transcription factors (TFs), vesicle polarized trafficking, cell wall biosynthesis, the ubiquitination machinery of the ubiquitin system and species-specific secondary metabolite pathways were mainly observed in the CK-VS-LT and CK-VS-NO libraries.

Conclusion

Differentially expressed unigenes related to the inhibition of C. sinensis pollen tube growth under low temperature and NO are identified in this study. The transcriptomic gene expression profiles present a valuable genomic tool to improve studying the molecular mechanisms underlying low-temperature tolerance in pollen tube.

Background

Tea (Camellia sinensis (L.) O. Kuntze) is a popular evergreen beverage tree, grown worldwide in different agro-climatic between 45°N and 34°S latitude [1]. Low-temperature stress is one of the most critical abiotic stresses that negatively influence the growth, development and geographical distribution of tea plants [2]. Previous studies have noted that low temperature lead to a weakened photosynthetic capacity and strengthened respiration rates [3], and a total of 105 metabolic processes were changed in tea plants, such as carbohydrate metabolism and amino acid metabolism [4], which ultimately decreased the yield and quality of the tea plants. Thus, finding the mechanisms of improving the tolerance of C. sinensis to low temperature is of great importance. Recently, certain molecular mechanisms driving the responses of tea plants to low temperature, such as Ca2+ signaling, APETALA2/ethylene responsive factor (AP2/ERF), basic helix-loop-helix (bHLH), WRKY and MYB domain-containing protein (MYB) transcription factors (TFs) and zinc finger protein (ZFP), have been demonstrated, particularly in vegetative tissues [5]. These findings provided a solid foundation for further exploring the cold-related genes and improved the understanding of plant tolerance to low-temperature stress. However, previous reports revealed that gene expression in reproductive tissues was different from that in vegetative tissues during cold tolerance [6], implying that reproductive tissue responds to low-temperature stress by different mechanisms [7]. Furthermore, the pollen tube, as a key component of genetic breeding, plays a crucial role in the higher plants’ reproduction, and its growth is highly susceptible to low temperature. For example, certain reports indicate that low temperature clearly reduce pollen tube elongation and disrupt its morphology [8, 9]. Interestingly, Camellia pollen can germinate in vitro even at 3 °C, whereas lily pollen cannot germinate at low temperature [10], implying higher tolerance to low temperature in Camellia pollen. Our previous report also demonstrated that C. sinensis pollen could germinate and grow at low temperature, although the germination and elongation were both retarded by cold stress [11]. Thus, the basic molecular mechanisms regulating pollen germination and tube elongation in response to low-temperature stress in C. sinensis deserve further attention.

Recently, nitric oxide (NO) has gained recognition for its significance in various essential physiological progresses in plants [12], and much attention has been focused on its roles in the plant tolerance to diverse abiotic stresses, including high salt [13], drought [14], heat [15], heavy metal [16], ultraviolet radiation [17], and cold stress [18]. Evidence suggests that NO production regulates diverse responses during cold stress. For example, NO production was increased in Triticum aestivum roots under treatment with cold stress [19], whereas it was reduced in Capsicum annuum leaves under treatment with low temperature [20]. Furthermore, NO has been confirmed to participate in the anabolic pathway and signal transduction of plant hormones under low-temperature stress. Guo et al. [21] revealed that NO participates in mediating the cold-induced expression of an enzyme related to the ethylene precursor MfSAMS1 (S-adenosylmethionine synthetase) by interacting with abscisic acid (ABA) and H2O2 in Medicago sativa subsp. Falcate. Recent reports have shown that NADPH oxidase inhibitors and reactive oxygen species scavengers reduced 24-epibrassinolide (EBR)-induced NO production, whereas EBR-induced H2O2 synthesis was not sensitive to NO. Decreased NO content inhibited EBR-induced low-temperature tolerance and partly inhibited EBR-induced expression of genes coding for several antioxidant enzymes and their relative activities in cucumber, suggesting that NO plays key roles in the H2O2-dependent induction of plant tolerance to low-temperature stress by brassinosteroids (BRs) [22]. Additionally, accumulating evidence has suggested that TFs are modulated by NO during cold stress. For instance, the C repeat binding factor 1 (CBF1) and CBF3 TFs are accumulated because of NO during Arabidopsis responding to low temperature [23]. Similarly, application of the NO donor sodium nitroprusside (SNP) increased LeCBF1 expression, whereas nitro-arginine (a NO synthesis competitive inhibitor) reduced its expression in tomato fruits during cold tolerance [24]. A recent study has indicated that NO participates in ABA-dependent cold signaling by negatively regulating the DNA binding of MYB2 by S-nitrosylation [25]. However, NO’s participation in low-temperature response is controlled by sophisticated signaling crosstalk for which the potential molecular mechanisms are still undeciphered.

As shown in the previous reports, NO plays crucial roles in pollen tube tip growing, for example its participation in the pollen tube growth regulation and reorientation in lily and Arabidopsis [26, 27]. In addition, extracellular nucleotides suppress Arabidopsis pollen germination and tube growth, which is modulated partly through the NO signaling pathway [28]. Wang et al. [29] reported that NO affects the cell wall construction in pollen tubes through the alteration of extracellular Ca2+ influx and F-actin organization, which consequently affects pollen tube elongation in Pinus bungeana. More interestingly, our previous report revealed that NO production from a NO synthase (NOS)-like enzymatic reaction enhanced the inhibiting effect of low-temperature stress on pollen germination and tube elongation in C. sinensis, and decreased free proline accumulation, which was accomplished partly through the guanosine 3’, 5’-cyclic monophosphate (cGMP) signaling pathway [11]. However, the downstream signaling and function mechanisms behind the NO-induced reduction of pollen tube elongation under low-temperature stress in C. sinensis remain unclear.

High-throughput transcriptomic analysis techniques can provide an overview of expression changes in response to biotic and abiotic stresses. In particular, RNA-Seq has become a widely used method for the analysis of species whose genome sequence is not yet available [30]. To understand the potential mechanisms driving NO-modulated pollen tube growth during exposure to low-temperature stress, we utilized the Illumina HiSeq™2000 platform to sequence the transcriptomes of C. sinensis pollen tubes cultured at 25 °C and low temperature and exposed to exogenous NO treatment, and all of the experiments were replicated three times to guarantee the reliability of the results in this study. Furthermore, the differential gene expression from the transcriptomes of the three treatments was systematically examined, the differentially expressed genes (DEGs) related to plant hormone signaling pathways, TFs, vesicle polarized trafficking, cell wall biosynthesis, the ubiquitination machinery of the ubiquitin system and species-specific secondary metabolite pathways were identified and analyzed. The assembled, annotated transcriptome will present a useful genomic resource for facilitating the investigations on the molecular mechanisms of NO’s involvement in pollen tube tolerance to low-temperature stress.

Results

Illumina sequencing and de novo assembly

To study the transcriptome expression of pollen tubes, C. sinensis pollen tubes were cultured under the following conditions: 25 °C (CK), 4 °C (LT) and 25 μM DEA NONOate (NO). Three libraries (CsPT-CK, CsPT-LT, and CsPT-NO) were designed for RNA-Seq using an Illumina HiSeqTM2000 genome analyzer, and all of the experiments were replicated three times. After cleaning and performing quality checks, 55.5, 52.4, and 51.3 million clean reads were obtained for CsPT-CK1, CsPT-CK2 and CsPT-CK3 (Accession number SRR3270364, SRR3270376 and SRR3270829 for library CsPT-CK1, CsPT-CK2 and CsPT-CK3), respectively; 50.2, 51.8, and 51.4 million clean reads were obtained for CsPT-LT1, CsPT-LT2 and CsPT-LT3 (Accession number SRR3270928, SRR3270974 and SRR3270993 for library CsPT-LT1, CsPT-LT2 and CsPT-LT3), respectively; and 51.9, 54.5, and 53.7 million clean reads were obtained for CsPT-NO1, CsPT-NO2 and CsPT-NO3 (Accession number SRR3270997, SRR3271001 and SRR3271002 for library CsPT-NO1, CsPT-NO2 and CsPT-NO3), respectively. The Q20 percentage was > 97.77 % in each library as shown in Table 1. The clean reads were assembled into 60,611, 61,246 and 63,117 contigs with mean lengths of 333, 335 and 331 nt for CsPT-CK, respectively; 56,102, 58,622 and 60,547 contigs with mean lengths of 343, 340 and 333 nt for CsPT-LT, respectively; and 59,750, 58,884 and 59,550 contigs with mean lengths of 331, 335 and 333 nt for CsPT-NO, respectively. These contigs were further assembled into 39,726, 40,440 and 41,626 unigenes with mean lengths of 648, 650 and 660 nt for CsPT-CK, respectively; 36,993, 39,070 and 39,439 unigenes with mean lengths of 650, 664 and 657 nt for CsPT-LT, respectively; and 39,514, 38,298 and 39,061 unigenes with mean lengths of 647, 650 and 639 nt for CsPT-NO, respectively (Table 2). These unigenes size distribution was displayed in Additional file 1: Figure S1. As shown in Fig. 1, the Pearson product moment coefficient was > 99.65 % in each library. All above data indicated that the quality of throughput and sequencing was high enough for further analysis.

Table 1 Summary of the output in the transcriptomes of pollen tubes
Full size table
Table 2 Summary of the assembly quality of the transcriptomes of the pollen tubes
Full size table
Fig. 1
figure1

Pearson correlation analysis of the replicates from the control (CK), low-temperature and NO treatments on the Camellia sinensis pollen tubes. a, b and c are analyses of the CK treatment; d, e and f are analyses of the low-temperature (LT) treatment, and g, h and i are analyses of the NO donor DEA NONOate (NO) treatment.

Full size image

Annotation of predicted proteins

We annotated 36,097 unique sequences on the basis of a BLAST search of six public databases: NCBI non-redundant (NR) database, mysql-nt (NT) database, Swiss-Prot protein database, Kyoto Encyclopedia of Genes and Genomes (KEGG) database, Cluster of Orthologous Groups of proteins (COG) database and Gene Ontology (GO) database. As shown in Fig. 2, we annotated 33,716 unique sequences using the NR database as a reference. Based on these NR annotations, 48.4 % of the sequences had a very strong homology to available plant sequences (E-value < 10−60), 20.4 % showed a strong homology (10−60 < E-value < 10−30), and an additional 31.3 % showed homology (10−30 < E-value < 10−5) (Fig. 2a). Similarity distributions were comparable, with 31.5 % of the sequences presenting similarities higher than 80 % and 69.5 % presenting similarities of 17–80 % (Fig. 2b). Furthermore, we found 43.7 % of the sequences matching homologues from Vitis vinifera, with additional hits to Amygdalus persica (10.1 %), Lycopersicon esculentum (10.1 %), Ricinus communis (7.9 %), and Populus balsamiferas ubsp. trichocarpa (6.7 %) (Fig. 2c).

Fig. 2
figure2

Homology search of unigenes against the NR database. a and b show the E-value and similarity distribution of the top BLAST hits for each unique sequence, respectively; c shows species distribution of the top BLAST hits for every homologous sequence

Full size image

Functional classification by COG

We searched all unigenes against the COG database for functional predictions and classifications. Our query assigned 13,029 of the 33,716 sequences that presented NR hits to COG classifications (Fig. 3). The COG-annotated putative proteins were functionally classified into at least 25 molecular families, including cellular structure, biochemistry metabolism, molecular processing, and signal transduction. Among them, the cluster “General function prediction” represented the largest group (4100; 16.96 %), followed by “Transcription” (2073; 8.57 %), “Replication, recombination and repair” (2064; 8.54 %), “Posttranslational modification, protein turnover and chaperones” (1992; 8.24 %), “Signal transduction mechanisms” (1713; 7.09 %), “Translation, ribosomal structure and biogenesis” (1560; 6.45 %), and “Carbohydrate transport and metabolism” (1503; 6.22 %); Only several unigenes were assigned to “Extracellular structure”.

Fig. 3
figure3

COG functional classification of the transcriptomes. Among 33,716 hits in the NR database, 13,029 unigenes with significant homologies in the COG database (E-value < 0.00001) were classified into 24 COG categories

Full size image

DEGs and GO classification

We filtered the DEGs (Additional file 2: Table S1, Additional file 3: Table S2 and Additional file 4: Table S3) using the absolute values of log2Ratio > 1 and probability > 0.7 as the thresholds according to Zhao et al. [31]. Comparing the library CsPT-CK with CsPT-LT (CK-VS-LT) showed that there were 766 DEGs (243 up-regulated and 523 down-regulated genes, 243/523). Moreover, we found 497 (222/275) and 929 (545/384) DEGs in the comparison between the library CsPT-CK and CsPT-NO (CK-VS-NO) and between the library CsPT-LT and CsPT-NO (LT-VS-NO), respectively (Fig. 4a). Among these DEGs, 357, 247 and 496 genes were specifically expressed in the CK-VS-LT, CK-VS-NO and LT-VS-NO, respectively. In addition, we observed 101, 284 and 125 genes co-expressed in CK-VS-LT and CK-VS-NO, CK-VS-LT and LT-VS-NO and CK-VS-NO and LT-VS-NO, respectively, and 24 genes were co-expressed in the three comparisons (Fig. 4b).

Fig. 4
figure4

The number of DEGs identified in the comparisons of CK-VS-LT, CK-VS-NO and LT-VS-NO. CK: control; LT: low-temperature treatment; NO: NO donor DEA NONOate treatment. CK-VS-LT: comparison between CK and LT. CK-VS-NO: comparison between CK and NO. LT-VS-NO: comparison between LT and NO. The absolute values of log2Ratio > 1 and probability > 0.7 were used as threshold for assigning significance. a the number of up-regulated (blue) and down-regulated (red) DEGs in the comparisons of CK-VS-LT, CK-VS-NO and LT-VS-NO are shown, respectively. b the number of specifically and co-expressed DEGs from the comparisons of CK-VS-LT, CK-VS-NO and LT-VS-NO are shown, respectively

Full size image

GO classifies genes on the basis of three functional ontologies, that is, Molecular function, Cellular component and Biological process. We used Blast2GO [32] to obtain GO annotations and Web Gene Ontology Annotation Plot (WEGO) [33] to perform GO functional classifications. For the CK-VS-LT, CK-VS-NO and LT-VS-NO libraries, 262 out of 766 DEGs (262/766), 156 out of 497 DEGs (156/497) and 328 out of 929 DEGs (328/929) could be assigned a GO classification, respectively (Additional file 5: Table S4, Additional file 6: Table S5 and Additional file 7: Table S6). The 262, 156 and 328 DEGs were assigned to 1957, 1239 and 2474 GO classifications in the CK-VS-LT, CK-VS-NO and LT-VS-NO libraries, of which 952, 593 and 1219 as “Biological process”, 705, 459 and 920 as “Cellular component” and 300, 187 and335 as “Molecular function”, respectively (Fig. 5). As shown in Fig. 5, the majority of DEGs in the “Biological process” category were associated with cellular processes, metabolic processes and single-organism processes, most DEGs in the “Cellular component” category were associated with cells, cell parts and organelles, and most in the “Molecular function” category were associated with binding and catalytic activity. Furthermore, between the CK-VS-LT and CK-VS-NO libraries, certain CK-VS-LT DEGs annotated with the “Biological process” category were related to biological adhesion and locomotion, whereas several CK-VS-NO DEGs annotated with the “Molecular function” category were associated with enzyme regulator activity and protein binding transcription factor activity.

Fig. 5
figure5

GO functional classifications of DEGs found from the comparisons of CK-VS-LT, CK-VS-NO and LT-VS-NO are annotated in the categories of Biological process (blue), Cellular component (red) and Molecular function (green). a GO classifications of DEGs from comparison of CK-VS-LT; b: GO classifications of DEGs from comparison of CK-VS-NO; c: GO classifications of DEGs from comparison of LT-VS-NO

Full size image

Cluster analysis of DEGs co-expressed in CK-VS-LT and CK-VS-NO

Genes are usually functionally correlated to those with similar expression patterns. In order to understand the expression patterns of the 125 co-expressed DEGs (Additional file 8: Table S7) in the CK-VS-LT and CK-VS-NO libraries, a cluster analysis of the gene expression patterns between these two libraries were performed. We classified these genes into three groups (Fig. 6). The largest group (Group 1) contained 71 genes (56.8 %) (from Unigene 11102_All to Unigene 20419_All) (Additional file 8: Table S7) that were up-regulated in both CK-VS-LT and CK-VS-NO. This group mainly contained genes encoding proteins related to oxidation-reduction reactions and amino acid metabolism, such as cytochrome c oxidase, arginine/serine-rich protein and serine/threonine-protein kinase. The second largest group (Group 3) included 52 genes (41.6 %) (from Unigene 10170_All to Unigene 8447_All) that were down-regulated in both CK-VS-LT and CK-VS-NO. This group mainly contained genes encoding proteins related to plant hormone metabolism and signaling pathways and TFs, such as ethylene receptor, ABA 8’-hydroxylase (ABA8ox), MADS-box protein and CaM-binding transcription activator (CAMTA). Finally, two genes (Group 2) were down-regulated in CK-VS-LT but up-regulated in CK-VS-NO, and one of them (Unigene 11403_All) was annotated as a U-box domain-containing protein.

Fig. 6
figure6

Clustering analysis of the DEGs co-expressed in CK-VS-LT and CK-VS-NO. This analysis is based on the log2Ratio values of DEGs. Columns displays the two comparisons (CK-VS-LT and CK-VS-NO), and rows represent the 125 DEGs. Up- and down-regulated gene expressions are showed by red and green, respectively. The 125 genes exhibit three expression patterns

Full size image

DEGs identified from pollen tubes treated with low temperature and NO

In this study, numerous genes were observed that exhibited changed expression levels when pollen tubes were exposed to low temperature or NO. These genes included factors related to plant hormone metabolism and signaling pathways, such as ABA8ox and type 2C protein phosphatase (PP2C) which participate in the ABA metabolism and signaling network (Tables 3 and 4); TFs, such as MYB and ZFP transcription factors (Tables 3 and 4); vesicle polarized trafficking and cell wall biosynthesis (Additional file 9: Table S8 and Additional file 10: Table S9); the ubiquitination machinery of the ubiquitin system (Additional file 11: Table S10 and Additional file 12: Table S11); and species-specific secondary metabolite pathways (Additional file 13: Table S12 and Additional file 14: Table S13). Furthermore, many of these DEGs responded to both low temperature and NO.

Table 3 DEGs related to plant hormone metabolism and signaling pathways and transcription factors (TFs) between CK and LT (CK-VS-LT)
Full size table
Table 4 DEGs involved in the plant hormone metabolism and signaling pathways and transcription factors (TFs) between CK and NO (CK-VS-NO)
Full size table

Quantitative real-time PCR (qRT-PCR) validation of DEGs in CK-VS-LT and CK-VS-NO

To confirm the gene expression patterns shown in the sequencing data, we performed a qRT-PCR analysis on nine randomly selected genes in pollen tubes exposed to low-temperature stress, 25 μM DEA NONOate (NO donor) or to low temperature with 200 μM 2-(4-carboxyphenyl)-4,4,5,5-tetramethylimidazoline-1-oxyl-3-oxide (cPTIO) (NO scavenger). According to our preliminary analysis and previous reports, we speculated that Ca2+ signaling and TFs played an important role in this process [11, 25]. Therefore, we selected CAMTA (Unigenes 2215), MYB (Unigenes 17702), ZFP (Unigenes 21920), AP2 (CL2776.Contig1) and PP2C (CL580.Contig48) genes. Furthermore, Phospholipase D (PLD) accumulated highly when Chorispora bungeana was exposed to low temperature [34] and PLD activation in C. bungeana exposed to low-temperature stress was related to a drop in Ca2+ content in membrane fractions [35]. Moreover, TOPLESS (TPL) was involved in mutual adjustment between hormones via NO in Petunia hybrida [36]. Additionally, genes encoding a non-specific serine/threonine protein kinase were down-regulated expressed during response to low temperature in two grapevine varieties [37]. Recent research indicated that phytochrome B (PHYB) and PHYA antagonistically regulated tomato low-temperature tolerance, in a process involving far-red light-induced PHYA activation to induce ABA signaling and eventually lead to cold response [38]. Here, we also randomly selected PLD (Unigenes 22243), TPL (Unigenes 8699) and serine/threonine protein kinase (CL10.Contig2) in CK-VS-LT and CK-VS-NO, and the PHYB (CL2425.Contig1) in CK-VS-LT for our qRT-PCR assay. The expression of Unigene 8699, CL10.Contig2 and CL580.Contig48 was up-regulated, whereas Unigenes 2215, 17702, 21920 and 22243 and CL2425.Contig1 and CL2776.Contig1 were down-regulated expressed in CK-VS-LT and CK-VS-NO (Fig. 7). These analyses support our DEGs data, and the high confirmation rate further confirmed the reliability of the data. In addition, application of 200 μM cPTIO under low-temperature stress rescued the changes in the expression levels of the DEGs co-expressed in the three treatments, such as Unigenes 2215, 8699, 17702, 21920 and 22243; CL10.Contig2 and CL580.Contig48. These results indicated that the seven genes intervened in pollen tube response to low temperature, possibly via the regulation of NO signaling.

Fig. 7
figure7

QRT-PCR analysis confirms the expression patterns indicated by the sequencing. Nine genes were randomly selected for qRT-PCR analysis and their primers were presented in Additional file 16: Table S14. The C. sinensis 18S rRNA was employed as reference

Full size image

Discussion

According to our previous report, NO negatively regulated C. sinensis pollen tube growth under low-temperature stress partially through the cGMP signaling pathway [11]. However, the potential molecular mechanism during this process is still unknown. Here, we comprehensively analyzed the gene expression changes in C. sinensis pollen tubes treated with low temperature and exogenous NO by RNA-Seq technology to investigate the pollen tube’s defense response to low temperature. Genes related to metabolism and signaling pathways of plant hormones, TFs, vesicle polarized trafficking and cell wall biosynthesis, the ubiquitination machinery of the ubiquitin system and species-specific secondary metabolite pathways were the main focus of this study.

Comparison of unigenes between C. sinensis leaves and pollen tubes

According to previous reports, Illumina sequencing of leaves transcriptome from C. sinensis cv. ‘Longjing43’ revealed an average transcript size of 356 bp with a N50 of 529 bp [5]. In this study, we reported that the average transcript size was > 331 bp, and the N50 was > 559 bp in every transcriptome of pollen tubes from C. sinensis cv. ‘Longjingchangye’. In addition, the previous study showed the annotation of 53,201 unigenes based on BLASTx searches against NR (51,883 unigenes), UniRef90 (52,217 unigenes), the Arabidopsis Information Resource (TAIR, version 10) (42,969 unigenes), KEGG (15,086 unigenes) and COG (29,356 unigenes) databases in tea plant leaves [5]. Among 71,289 transcripts based on NR annotation, 30,115 transcripts matched 80,176 GO term annotations [5]. Most of the unigenes both from leaves and pollen tubes are categorized as the cell and cell part (cellular component), binding and catalytic activity (molecular function) and cellular process and metabolic process (biological process) (Additional file 15: Figure S2). Interestingly, envelope (cellular component), auxiliary transport protein, transcription regulator (molecular function) and anatomical structure formation, cell killing, death, pigmentation, viral reproduction (biological process) were only annotated in GO categories of tea leaves transcriptome (Additional file 15: Figure S2A). However, cell junction, extracellular matrix, extracellular matrix part, membrane, membrane part, nucleoid, symplast (cellular component), channel regulator, metallochaperone, nucleic acid binding transcription factor, nutrient reservoir, protein binding transcription factor, protein tag, receptor (molecular function) and regulation of biological (biological process) were only annotated in GO categories of pollen tube transcriptome (Additional file 15: Figure S2B). These differences between transcriptomes from leaves and pollen tubes might result from the diversity of reproductive and vegetative tissues of plant [6], further confirming that reproductive tissue responds to low temperature stress using distinct mechanisms from vegetative tissues [7].

Plant hormone metabolism and signaling pathways

Plant hormones play essential roles in response to abiotic stress through their involvement in complicated signaling networks and physiological processes. For example, ABA has been shown to accumulate as a protective factor against cold injury in maize [39], and previous investigations have demonstrated that NO enhances tolerance to drought stress via promoting ABA-induced stomata closure [40, 41]. In our study, we identified one ABA8ox gene (Unigene14044_All) and one PP2C gene (CL580.Contig48_All) in the CK-VS-LT (Table 3) and CK-VS-NO libraries (Table 4). As shown in previous reports, NO-induced ABA content reduction in Arabidopsis correlates with the modulation of ABA8ox protein expression [42]. Additionally, ABA signaling is negatively regulated by NO in germination and early seedling growth in Arabidopsis through the S-nitrosylation of sucrose non-fermenting 1 (SNF1)-related protein kinase 2.2 (SnRK2.2) and SnRK2.3 [43]. As a negative regulator, PP2C inactivates SnRK2, thereby regulating ABA signals [44]. These reports and our studies indicate that NO might participate in ABA catabolism and signaling to regulate pollen tube elongation under cold stress through the genes ABA8ox (Unigene14044_All) and PP2C (CL580.Contig48_All). In addition, the plant steroidal hormone BRs has been found to promote growth and cell elongation and protect plants against abiotic stresses. Reports have shown that BRs exert positive effects on the mitigation of oxidative damage by enhancing the antioxidant defense system in Chorispora bungeana under cold stress [45]. BRs are recognized through an active complex of a leucine-rich repeat (LRR) receptor-like kinase (RLK) brassinosteroid-insensitive 1 (BRI1) and BRI1-associated kinase 1 (BAK1). BRs bind to the extracellular LRR domain of BRI1 and induce a phosphorylation-mediated cascade to modulate gene expression [46]. Meanwhile, BAK1 is involved in the heterodimerization and endocytosis associating with BRI1 to promote BRs-induced signal transducing instead of involved in binding to BRs. [47]. We identified two BRI1 genes (Unigene9712_All and Unigene4038_All) and one BAK1 gene (CL4145.Contig1_All) in CK-VS-LT (Table 3) as well as one BRI1 gene (Unigene4038_All) in CK-VS-NO (Table 4). Current studies have revealed that NO concentration in the root cells is increased by BRs, which is required for BR-induced alterations of the root architecture in Arabidopsis [48]. This finding suggests that BRs are involved in pollen tube resistance to low temperature through Unigene9712_All, Unigene4038_All and CL4145.Contig1_All and indicate that NO might participate in this process through Unigene4038_All. Moreover, BRs could overcome the inhibition of seed germination induced by ABA and increase the tolerance of Arabidopsis seedlings to cold [49]. BR-induced NO production and NO-activated ABA biosynthesis are significant mechanisms for BR-increased tolerance to water-stress [50]. Combining these reports and our findings, BRs and ABA might interact in NO signal transduction during the pollen tube responding to low-temperature stress, although the underlying mechanisms require further research.

Ethylene (ET) is a gaseous plant hormone, and its signaling pathways participate in plant responses to low temperature [51, 52]. Previous reports indicated that ET signaling negatively regulates Arabidopsis tolerance to cold stress by inhibiting the expression of CBF and type-A ARR genes [53]. Ethylene response sensor 1 (ERS1) has been identified as one of the membrane-located receptor proteins in the ET signaling pathway [54]. The level of ET is usually low in normal conditions, and its response is suppressed by ERS1 via activating the downstream negative regulator CTR1 (Constitutive Triple Response 1). Furthermore, EIN2 (ethylene insensitive 2) protein acting as a critical positive regulator of ET signaling has been found to be inhibited by CTR1 in Arabidopsis [55]. As downstream of EIN2, the plant-specific transcription factors EIN3 and EIN3-like1 (EIL1) are sufficient and necessary to stimulate ethylene-response genes expression and to modulate the ethylene-related responses of plants [56]. EIN3-Binding F-box-1 (EBF1) and EBF2 tightly regulate the stability of EIN3 at the protein level [57], and EBF2 transcription is in turn activated by EIN3, thereby forming a negative feedback loop [58]. Interestingly, recent work has shown that ERS1 and EBF1/2 are differentially expressed in Papaya fruit ripening disorders caused by cold injury [59]. In our study, two ERS1 genes (Unigene4486_All and Unigene22534_All) and two EBF1/2 genes (CL2610.Contig1_All and Unigene10222_All) were down-regulated between CK and LT (CK-VS-LT) (Table 3). These findings imply that ERS1 and EBF1/2 are involved in pollen tube resistance to low temperature. Although few reports on the interaction between the two gaseous hormones NO and ET during low-temperature stress are available, studies have suggested that NO and ET have a complicated relationship in plants during water logging stress [60], salt stress [61], UV-B stress [62] and hypersensitive responses [63]. One ERS1 gene (Unigene22534_All) was also identified in the CK-VS-NO library (Table 4). It is tempting to speculate that ERS1 (Unigene22534_All) responds to low temperature via the NO signaling pathway and the other ERS1 genes in C. sinensis pollen tubes (Unigene4486_All and EBF1/2 (CL2610.Contig1_All and Unigene10222_All) respond to low temperature in a NO-independent manner.

Auxin is an essential morphogenetic signal participating in the regulation of cell identity during all the developmental processes of plants, and auxin signaling pathways constitute a critical component of mechanisms for plant tolerance to abiotic stresses. Substantial evidence has demonstrated a tight link between auxin responses and cold stress [64, 65]. For example, cold stress affects auxin polar transport via the selective inhibition of the intracellular trafficking of numerous proteins in Arabidopsis, including the auxin efflux carriers [66]. Two auxin-repressed protein genes (CL2941.Contig1_All and CL2941.Contig2_All), one auxin response factor gene (Unigene7793_All), one indole-3-acetic acid (IAA)-amido synthetase gene (Unigene13496_All) and one IAA-induced protein gene (Unigene21651_All) were identified in the CK-VS-LT library in this study (Table 3). Interestingly, as shown in a recent report, auxin participates in regulating pollen tube tip growth [67]. Therefore, it is reasonable to speculate that auxin intervenes in pollen tube responding to low-temperature stress in C. sinensis. Furthermore, Markus et al. [68] reported that the release of IAA and NO from peroxisomes is a spatially and temporally coordinated process in root formation in maize and Arabidopsis, thereby indicating that auxin and NO interact closely in plant developmental processes. We identified one IAA-amido synthetase gene (Unigene17743_All) and one IAA-induced protein gene (Unigene21651_All) in CK-VS-NO (Table 4). Unigene21651_All was simultaneously detected in both CK-VS-LT and CK-VS-NO libraries, indicating that IAA-induced protein (Unigene21651_All) might be related to the NO-mediated pathway in C. sinensis pollen tube growth in responding to low-temperature stress.

Additionally, gibberellins (GAs) are another class of plant hormones that play critical roles in modulating plant growth and development in response to environmental cues. Exposure of Arabidopsis seedlings to cold stress triggers a reduction in bioactive GA, which restrains root growth [69]. Previous studies have shown that the nuclear-localized growth-repressing DELLA proteins are central components of the GA-signaling pathway. These proteins accumulate through a posttranslational mechanism mediated by the decrease in bioactive GA via up-regulating the transcript level of GA 2-oxidase (GA2ox) gene, and this in turn increases the transcript levels of GA20ox and GA3ox gene by a feedback mechanism [70, 71]. The binding of bioactive GA to GA insensitive dwarf 1 (GID1), a pocket of the GA receptor, promotes the interaction between GID1 and DELLAs and forms the GA–GID1–DELLA complex [72]. This model has been employed to reveal the molecular mechanisms for GA perception. In this study, two DELLA genes (Unigene15959_All and Unigene10398_All), one GA2ox gene (Unigene11055_All), one GA20ox gene (Unigene1092_All) and one GID1 gene (Unigene9129_All) showed a significantly different expression in CK-VS-LT (Table 3). These findings suggest that GA participates in pollen tube response to low-temperature stress, whereas none of the GA signals or synthetic metabolism genes was detected in the DEGs of CK-VS-NO. In addition, María et al. [73] speculated that high levels of NO content could influence the increase of DELLA activity through the disappearance of PIN-FORMED 1 (PIN1) and subsequently inhibit cell elongation in the elongation-differentiation zone in Arabidopsis. However, the PIN1 gene is not detected in CK-VS-NO. Therefore, we hypothesize that NO might regulate C. sinensis pollen tube growth during responses to low temperature stress through a GA-independent signaling network.

Transcription factors (TFs)

Plants have evolved complex and highly efficient regulatory networks to respond and adapt to cold stress, and TFs play crucial roles in these regulatory processes. For example, the rice R2R3-type MYB gene, OsMYB2, plays a positive regulatory role in low-temperature stress tolerance [74]. Similarly, the soybean typical Cys2/His2-type (C2H2-type) zinc finger GmZF1 enhances cold stress tolerance by regulating gene expression in transgenic Arabidopsis [75]. In our study, two MYB genes and three ZFP genes were identified in CK-VS-LT as shown in Table 3. The results reveal that MYB and ZFP are involved in pollen tube elongation in response to low-temperature stress. Besides, NO action inhibits the DNA-binding of AtMYB2 through the S-nitrosylation of Cys53 [76]. Furthermore, the thiols of zinc-sulfur clusters can be S-nitrosated by NO, which reversibly disrupts zinc finger structures and presents a molecular mechanism for regulating ZFP TFs [77]. One MYB gene (Unigene17702_All) and one ZFP gene (Unigene21920_All) were detected among the comparisons of the CK-VS-LT and CK-VS-NO libraries (Tables 3 and 4), implying that MYB (Unigene17702_All) and ZFP (Unigene21920_All) are related to the NO signaling pathway during low-temperature stress responses in C. sinensis pollen tubes. Moreover, NO negatively regulates the DNA binding of MYB2 by S-nitrosylation in the ABA-dependent cold signaling pathway [25]. Combining this information with the results for ABA and NO, we speculate that MYB (Unigene17702_All) might participate in C. sinensis pollen tube growth during low-temperature stress via an ABA-dependent pathway that acts as a NO signaling pathway downstream regulator.

Pollen tube elongation is a tip growing process that can be affected by TFs. For example, the pollen-specific MIKC* class of MADS-box TFs is necessary for pollen maturation and tube elongation and controls a transcriptional switch that directs pollen maturation in Arabidopsis [78]. Furthermore, a large number of reports reveal that calmodulin (CaM) is a key element in pollen tube elongation and orientation and dependent on a complicated crosstalk process between multiple pathways [79–81]. The CaM-interacting protein group CAMTA has been identified in the plant defense signaling cascade [82]. Moreover, MADS-box [83] and CAMTA [84] TFs take part in the regulation of plant responses to low temperature. We detected five MADS-box genes and one CAMTA gene in CK-VS-LT (Table 3). These studies and our results indicate that CAMTA and MADS-box TFs participate in pollen tube elongation under cold stress. In CK-VS-NO, one MADS gene and two CAMTA genes were also detected, and in CK-VS-LT and CK-VS-NO, one MADS gene (CL1553.Contig1_All) and one CAMTA gene (Unigene2215_All) were differentially expressed (Tables 3 and 4). Previous experiments have demonstrated that NO modulates Ca2+ signaling in lily and Arabidopsis pollen tubes through directly imaging the concentration of cytosolic free Ca2+ during NO-induced tube re-orientation [26, 27]. We speculate that the CAMTA gene (Unigene2215_All) might participate in the regulation of Ca2+ by NO during pollen tube elongation under cold stress. Additionally, although directly relevant reports are not available on the crosstalk between MADS-box and NO during pollen tube growth, arginine is the substrate of NO synthases [85], and the metabolism of arginine is regulated by MADS-box [86]. Thus, it is tempting to suggest that the MADS gene (CL1553.Contig1_All) might influence C. sinensis pollen tube elongation under low temperature via the NO synthesis pathway.

In addition to the above TFs, we also identified two ERF genes, one WRKY gene and two bHLH genes in CK-VS-LT (Table 3) as well as one ERF gene and two bHLH genes in CK-VS-NO (Table 4). Using forward and reverse genetics, the ERF [51], WRKY [87] and bHLH [88] TFs were isolated and identified as the molecular switches of cold-responsive signaling networks in plants. Although these three TF family genes were found in CK-VS-LT and CK-VS-NO, none were simultaneously found in both CK-VS-LT and CK-VS-NO libraries. Thus, the complex networks and mechanisms of NO, ERF, WRKY, bHLH and low-temperature stress require further research.

Vesicle polarized trafficking and cell wall biosynthesis related genes

Pollen tube growth is a highly polarized process which is dependent on the polarized trafficking of vesicles containing cell wall materials to the tip region to establish a restricted growth zone [89]. Previous reports have shown that secretory vesicle formation in root hairs is impaired in the double phosphatidylinositol-4-kinase (PI4K) Arabidopsis mutant pi4kβ1/β2 because of the disruptions to phosphatidylinositol-4-phosphate (PI4P) gradient and polarized vesicle trafficking [90]. Furthermore, PI 4-phosphate 5-kinase (PIP5K) isoforms can catalyze the PI 4, 5-bisphosphate (PIP2) synthesis, which eventually regulates vesicle trafficking [91]. Additionally, Feng et al. [92] reported that vesicle fusion and recycling in the tip region of pollen tube apex might be mediated by vesicle-associated membrane protein 726 (VAMP726) in Petunia inflata. Reports have indicated that the vesicle trafficking process is accelerated in response to low temperature stress in Arabidopsis [93]. Here, we identified four PI4K genes, one PIP5K gene and two VAMP genes in CK-VS-LT (Additional file 9: Table S8), implying that these genes modulate pollen tube elongation during the response to cold stress. We also found that one of the PI4K genes (Unigene13861_All) is also differentially expressed in CK-VS-NO (Additional file 10: Table S9). Interestingly, NO involvement in pollen tube growth regulation has been reported to occur via the modulation of vesicle polarized trafficking [29], and NO is also required for vesicle formation and trafficking in root hair tip growth in Arabidopsis [94]. These results indicate that NO may modulate C. sinensis pollen tube elongation under low-temperature stress by mediating vesicle formation and trafficking via the PI4K gene (Unigene13861_All). Additionally, cell wall construction is among the most important mechanisms in pollen tube tip growth [95] and could be regulated by NO in Pinus bungeana [29]. COBRA-like genes are important factors for secondary cell wall development in the root growth zone of Arabidopsis [96]. We also identified one COBRA-like gene (Unigene10170_All) in the DEGs between CK-VS-LT and CK-VS-NO (Additional file 9: Table S8 and Additional file 10: Table S9). This COBRA-like gene (Unigene10170_All) might participate in cell wall construction in C. sinensis pollen tubes during low temperature stress through the NO signaling pathway.

Ubiquitination machinery of the ubiquitin system

The protein ubiquitination pathway is one of the most effective processes for performing versatile post-translational modifications that mediate the growth and development of all eukaryotic species. Ubiquitylated proteins are targeted for degradation by the 26S proteasome which is a proteolytic machine that degrades the target and allows ubiquitin moieties to be reused. Ubiquitin–26S proteasome system (UPS) widely intervenes in plant growth and development during tolerance to stresses by modulating proteins activity or affecting their localization or stability [97]. The UPS depends on the activity of ubiquitin activating enzymes (E1s), ubiquitin conjugating enzymes (E2s) and ubiquitin ligases (E3s). E3s are recognized as the critical factors that define substrate specificity [98] and are divided into four major types: Really Interesting New Gene (RING)-type, cullin–RING ligases (CRLs), Homology to E6-Associated Carboxyl-Terminus (HECT)-type and U-box-type. Furthermore, U-box-type E3s transfer ubiquitin tags directly from E2s-Ubs to target proteins [99]. Additionally, CRLs can consist of a substrate-recruiting protein, such as F-box [100] and broad complex, tramtrack, bric-a-brac/pox virus and zinc finger (BTB/POZ) [101]. In our study, we identified one E2 gene, four E3 genes, one 26S proteasome gene, four U-box genes, nine F-box genes and three BTB/POZ genes in the CK-VS-LT library (Additional file 11: Table S10). UPS modulates the levels of stress hormone biosynthesis and secondary messengers as well as the abundance of regulatory proteins that may be accumulated resulting from exposure to abiotic stress [102], which implies that the UPS plays a fundamental role in C. sinensis pollen tube elongation under low temperature stress. Moreover, most E3s have been linked to abiotic stress responses through their modulation of various stress signaling processes [98]. Previous reports have shown that multiple E3s regulate ABA-dependent stress signaling [103]. Combined with our findings on ABA and cold stress, this information indicates that the UPS is involved in an ABA-dependent signaling pathway and might influence other hormones and TFs through ABA-dependent signaling. Additionally, certain genes encoding E3 proteins are involved in NO-dependent modulation in response to cryptogein in Arabidopsis [104], indicating that NO and UPS interact closely in plants. In CK-VS-NO, four E3 genes, one 26S proteasome gene, two U-box genes, one F-box gene and one BTB/POZ protein gene were identified (Additional file 12: Table S11). Additionally, we also identified one E3s (Unigene6536_All), one 26S proteasome (Unigene21945_All), one U-box (Unigene11403_All) and one BTB/POZ (Unigene17396_All) gene co-expressed in the DEGs of CK-VS-LT and CK-VS-NO. These results indicate that the E3s (Unigene6536_All), 26S proteasome (Unigene21945_All), U-box (Unigene11403_All) and BTB/POZ (Unigene17396_All) might participate in an NO-dependent signaling pathway in C. sinensis pollen tubes under low temperature stress through their influence on the modulation of the UPS, which alters the regulatory proteins levels.

Species-specific secondary metabolites pathways

Major secondary metabolites are fundamental components that define the flavor, taste and health benefits of tea, such as flavonoids, caffeine and theanine. Flavonoids contain flavones, flavonols, isoflavones, flavanones, flavanols, and anthocyanidins, which have various functions in plant tolerance to stresses, including defense against phytopathogens, protection from ultraviolet light and antioxidant activity [105]. Furthermore, the accumulation of catechins in C. sinensis is very susceptible to low temperature [106]. In addition, caffeine is a purine alkaloid widely used as stimulant and pharmaceutical component, and mainly synthesized in young leaves in C. sinensis via a typical biosynthetic pathway that includes steps of purine synthesis and modification [107]. Theanine is emerged as a unique free amino acid that confers the unique “umami” taste of tea, and it occupies about half of the total free amino acids in tea [108]. It has been revealed that theanine biosynthesis starts from glutamine and pyruvate and includes further synthetic steps in buds, leaves and roots from C. sinensis [109, 110]. In the present study, some genes related to flavonoids, caffeine and theanine biosynthesis pathways were down-regulated in CK-VS-LT (Additional file 13: Table S12). In addition, we also identified some genes in CK-VS-NO (Additional file 14: Table S13). And one caffeine synthase 1 (Unigene1362_All) and glutamate receptor gene (CL4694.Contig1_All) were co-expressed in the comparisons of CK-VS-LT and CK-VS-NO. These results reveal that caffeine and theanine mediate cold stress tolerance via NO signal pathway by caffeine synthase 1 (Unigene1362_All) and glutamate receptor (CL4694.Contig1_All).

Conclusions

We produced a dataset containing 45,432 unigenes from the C. sinensis pollen tube transcriptomes using paired-end and de novo sequencing with the Illumina HiSeq™2000 platform. Among these unigenes, 36,097 were annotated with descriptions from the NR, NT, Swiss-Prot, KEGG, COG and GO databases. Our dataset includes comprehensive analysis in sequence and DEGs profiling data that provide a dynamic perspective on transcriptomic variations caused by NO and low-temperature in C. sinensis pollen tube elongation. Many DEGs were simultaneously found in both the CK-VS-LT and CK-VS-NO libraries in different functional pathways, such as plant hormone signaling (ABA, BRs, ET and auxin), TFs (MYB, ZFP, CAMTA and MADS-box), vesicle polarized trafficking and cell wall biosynthesis (PI4K and COBRA-like), the ubiquitination machinery of the ubiquitin system (E3s, 26S proteasome, U-box and BTB/POZ) and caffeine and theanine biosynthesis (caffeine synthase and glutamate receptor). Combining our findings and previous reports, we propose a crosstalk network which may indicate the probable mechanisms underlying NO mediating C. sinensis pollen tube growth during the response to low-temperature stress (Fig. 8) [102, 111]. Plant hormone signaling pathways (ABA, BRs, ET and auxin) respond to pollen germination and tube elongation at low temperature, and then TFs (MYB, ZFP, CAMTA and MADS-box) modulated relevant proteins expression. Moreover, the UPS (E3s, 26S proteasome, U-box and BTB/POZ) probably involves in affecting pollen germination and tube growth in response to low-temperature stress by modulating the activity, localization or stability of different proteins. Additionally, Ca2+ gradient, vesicle polarized trafficking as well as cell wall biosynthesis might participate in pollen tube elongation under low-temperature stress through the NO signaling pathway by CAMTA, PI4K and COBRA-like. Taken together, these findings offer a more profound understanding of the candidate genes employed by C. sinensis pollen tubes to integrate low temperature signals and finely tune their tolerance responses, present a fundamental transcriptomic resource to further explore the molecular mechanisms of reproductive tissue tolerance to low-temperature.

Fig. 8
figure8

Model of the potential signaling pathway of NO in C. sinensis pollen tube response to low temperature. This summarized model was based on the results reported by Wang et al. [111] and Lyzenga ea al. [102]. Plant hormone signaling pathways (ABA, BRs, ET and auxin) respond to pollen germination and tube elongation at low temperature, and then TFs (MYB, ZFP, CAMTA and MADS-box) modulated relevant proteins expression. Moreover, the UPS (E3s, 26S proteasome, U-box and BTB/POZ) probably involves in affecting pollen germination and tube growth in response to low-temperature stress by modulating the activity, localization or stability of different proteins. Additionally, Ca2+ gradient, vesicle polarized trafficking as well as cell wall biosynthesis might participate in pollen tube elongation under low-temperature stress through the NO signaling pathway by CAMTA, PI4K and COBRA-like

Full size image

Methods

Plant materials

Mature pollen was collected from C. sinensis (L.) O. Kuntze cv. ‘Longjingchangye’ growing in Dr. Sun Yat-sen’s Mausoleum Tea Factory, Nanjing, China (http://www.yhc1958.com/index.html). Then the pollen was pre-cultured in a control liquid medium (5 % sucrose, 0.05 % Ca(NO3)2  · 4H2O, 0.01 % H3BO3, 5 % PEG 4000 and 30 mM MES, pH 6.0) at 25 °C in darkness for 30 min [111]. The pre-cultured samples were incubated at 25 °C for 1 h as control (CK treatment). The pre-cultured pollen was incubated at 4 °C for 1 h as low-temperature treatment (LT treatment). In addition, adding 25 μM NO donor DEA NONOate to the pre-cultured pollen and incubated it at 25 °C for 1 h as NO treatment (NO treatment). For NO scavenger treatment, the pre-cultured pollen was incubated in the medium with 200 μM cPTIO at 4 °C for 1 h (LT with cPTIO treatment). All experiments were replicated three times. Pollen tubes snap frozen in liquid nitrogen after collection and then stored at −80 °C before use.

RNA extraction, library preparation and RNA-Seq

Nine pollen tube samples from three independent experiments were collected after incubated under control, low temperature and NO treatment. For Illumina sequencing, total RNA was extracted and treated with RNAiso Plus (TaKaRa, Japan) and RNase-free DNase I (Takara Biotechnology, China), respectively, according to the manufacturer’s instructions. We confirmed the integrity, quality and quantity of total RNA using Agilent 2100 Bioanalyzer (Agilent, Santa Clara, CA, USA) and a spectrophotometer (NanoDrop, Wilmington, DE), and used TruSeq RNA Sample Prep Kit v2 (Illumina) to purify 200 ng total RNA for each sample by oligo-dT. Then the total RNA was fragmented using Elute, Prime and fragmentation buffer. We produced the first-strand cDNA using First Strand Master Mix and Super Script II (Invitrogen) and created the second-strand cDNA using Second Strand Master Mix. The double-stranded DNA was then purified by QiaQuik PCR purification kit (Qiagen, Beijing, China) before end repaired by incubation with End Repair Mix at the 30 °C for 30 min. Subsequently, the end-repaired cDNA was added single A base followed by adapter ligation. We performed PCR with PCR Primer Cocktail and PCR Master Mix to amplify the cDNA. The final library was quantified determining the average molecule length measured by the Agilent 2100 bioanalyzer (Agilent DNA 1000 Reagents) and via QPCR (TaqMan Probe), the cDNA library was eventually sequenced on the HiSeq 2000 System (TruSeq SBS KIT-HS V3,Illumina) in the Beijing Genomics Institute (Shenzhen, China; http://www.genomics.cn/index). Above all, the RNA extraction, library preparation and RNA-Seq were carried out according to descriptions showed by Wang et al. [5] and Wang et al. [111].

De novo assembly and data filtering

Raw reads produced directly from the sequencing platform include dirty reads containing unknown or low-quality bases and adapters. After removing the reads with adaptors, the reads with more than 5 % unknown nucleotides and the low-quality reads (>20 % read rate with a quality value ≤ 10), the remaining clean reads were used for the subsequent analysis [112]. Trinity [112] software was then used to assemble short strips. The sequences resulting from Trinity are called unigenes. Using TGICL [113] to further account for splicing and redundancy, these sequences were assembled into homologous transcript clusters. The same process was applied to all three samples to obtain the longest possible non-redundant unigenes. Based on the clustering of the homologous transcripts, the unigenes were divided into clusters composed of several high-similarity (>70 %) unigenes (starting with CL) and singletons (starting with unigenes). Finally, the unigenes were aligned by BLASTX (E-value < 0.00001) to protein databases (in order of priority) NR, Swiss-Prot, KEGG and COG. The proteins showed the highest ranks in the results of BLAST search were used to determine the coding region sequences of the unigenes. We used ESTS software [114] to determine the sequence direction, if a unigene could not align with any of the above databases. Thereby the nucleotide sequence in the direction from 5’ to 3’ and the amino sequence of the predicted coding region are produced. We have deposited all sequencing data at the sequence read archive of NCBI (Accession number SRR3270364, SRR3270376, SRR3270829, SRR3270928, SRR3270974, SRR3270993, SRR3270997, SRR3271001 and SRR3271002).

DEGs selection

We used Tags per Million (TPM) to normalize clean tags of each library to obtain the normalized gene expression levels. DEGs were determined using the NOISeq-bio method as previously reported by Tarazona et al. [115]. We estimated the fold changes (log2Ratio) on the basis of the normalized gene expression levels in each sample. We used the absolute values of log2Ratio > 1 and probability > 0.7 as the thresholds criteria for significant differences in gene expression [31]. In addition, DEGs related to plant hormone metabolism and signaling, TFs and the vesicle polarized trafficking, cell wall biosynthesis, the ubiquitination machinery of the ubiquitin system and species-specific secondary metabolites pathways were mainly analyzed.

Unigene annotation and classification

Bioinformatics procedures were employed to classify the annotation of unigenes. BLAST was used against the NT, NR, KEGG, Swiss-Prot and COG (E-value ≤ 1.0E-5) databases for annotation, and subsequently the number of unigenes annotated within each database was counted. Based on the NR annotation, GO functional annotations can be obtained. First, the Blast2GO [32] to get the GO annotation of the unigenes, while the WEGO [33] was used to analyze GO functional classifications for all unigenes and determine the distribution of species gene functions at the macro level. The KEGG annotation was employed to get pathway annotations for our unigenes.

Clustering of gene expression profiles

We performed a hierarchical cluster analysis using the log2Ratio of the 125 co-expressed DEGs in the CK-VS-LT and CK-VS-NO expression libraries by cluster [116] and Java Treeview [117] software.

qRT-PCR analysis

The qRT-PCR analysis was carried out according to the description shown by Wang et al. [111]. The primers used in this assay were shown in Additional file 16: Table S14, and the C. sinensis 18S rRNA was employed as the reference.

Statistical analyses

Statistical analyses were performed as previous reported by Wang et al. [111]. Briefly, all experiments were replicated at least three times, and all data were represented as the means ± standard deviations (SD). Group differences were tested using one-way ANOVA and Duncan’s test, and significant differences are represented by different letters (P < 0.05). Data analysis was performed using SPSS 20 software.

Abbreviations

ABA:

Abscisic acid

ABA8ox:

ABA 8’-hydroxylase

AP2/ERF:

APETALA2/ethylene responsive factor

BAK1:

BRI1-associated kinase 1

bHLH:

Basic helix-loop-helix

BRI1:

Brassinosteroid-insensitive 1

BRs:

Brassinosteroids

BTB/POZ:

Broad complex tramtrack, bric-a-brac/pox virus and zinc finger

CaM:

Calmodulin

CAMTA:

CaM-binding transcription activator

CBF:

C repeat binding factor

cGMP:

Guanosine 3’ 5’-cyclic monophosphate

COG:

Cluster of Orthologous Groups of proteins database

cPTIO:

2-(4-carboxyphenyl)-4,4,5,5-tetramethylimidazoline-1-oxyl-3-oxide

CRLs:

Cullin–RING ligases

CTR1:

Constitutive Triple Response 1

DEGs:

Differentially expressed genes

E1s:

Ubiquitin activating enzymes

E2s:

Ubiquitin conjugating enzymes

E3s:

Ubiquitin ligases

EBF:

EIN3-Binding F-box

EBR:

Epibrassinolide

EIL1:

EIN3-like1

EIN:

Ethylene insensitive

ERS1:

Ethylene response sensor 1

ET:

Ethylene

GA20x:

GA 20oxidase

GA2ox:

GA 2-oxidase

GA3x:

GA 3oxidase

GAs:

Gibberellins

GID1:

GA insensitive dwarf 1

GO:

Gene Ontology database

HECT:

Homology to E6-Associated Carboxyl-Terminus

IAA:

Indole-3-acetic acid

KEGG:

Kyoto Encyclopedia of Genes and Genomesdatabase

LRR:

Leucine-rich repeat

MYB:

MYB domain-containing protein

NO:

Nitric oxide

NOS:

NO synthase

NR:

NCBI non-redundant database

NT:

Mysql-nt database

PHYA:

Phytochrome A

PHYB:

Phytochrome B

PI4K:

Phosphatidylinositol-4-kinase

PI4P:

Phosphatidylinositol-4-phosphate

PIN1:

PIN-FORMED 1

PIP2:

PI 4,5-bisphosphate

PIP5K:

PI 4-phosphate 5-kinase

PLD:

Phospholipase D

PP2C:

Type 2C protein phosphatase

qRT-PCR:

Real-time quantitative polymerase chain reaction

RING:

Really Interesting New Gene

RLK:

LRR receptor-like kinase

SAMS:

S-adenosylmethionine synthetase

SNF1:

S-nitrosylation of sucrose non-fermenting 1

SNP:

Sodium nitroprusside

SnRK:

SNF1-related protein kinase

TAIR:

The Arabidopsis Information Resource

TFs:

Transcription factors

TPL:

TOPLESS

TPM:

Tags per Million

UPS:

Ubiquitin–26S proteasome system

VAMP726:

Vesicle-associated membrane protein 726

WEGO:

Web Gene Ontology Annotation Plot

ZFP:

Zinc finger protein

References

  1. 1.

    Jain NK. Global advances in tea science. New Delhi: Aravali Books International Pvt Ltd.; 1999.

  2. 2.

    Vyas D, Kumar S. Tea (Camellia sinensis (L.) O. Kuntze) clone with lower period of winter dormancy exhibits lesser cellular damage in response to low temperature. Plant Physiol Bioch. 2005;43(4):383–8.

    CAS  Article  Google Scholar 

  3. 3.

    Aoki S. Interaction of light and low-temperature in depression of photosynthesis in tea leaves. Jpn J Crop Sci. 1986;55(4):496–503.

    Article  Google Scholar 

  4. 4.

    Shen J, Wang Y, Chen C, Ding Z, Hu J, Zheng C, Li Y. Metabolite profiling of tea (Camellia sinensis L.) leaves in winter. Sci Hortic-Amsterdam. 2015;192:1–9.

    CAS  Article  Google Scholar 

  5. 5.

    Wang X, Zhao Q, Ma C, Zhang Z, Cao H, Kong Y, Yue C, Hao X, Chen L, Ma J, et al. Global transcriptome profiles of Camellia sinensis during cold acclimation. BMC Genomics. 2013;14:415.

  6. 6.

    Lee JY, Lee DH. Use of serial analysis of gene expression technology to reveal changes in gene expression in Arabidopsis pollen undergoing cold stress. Plant Physiol. 2003;132(2):517–29.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  7. 7.

    Chinnusamy V, Zhu J, Zhu JK. Gene regulation during cold acclimation in plants. Physiol Plantarum. 2006;126(1):52–61.

    CAS  Article  Google Scholar 

  8. 8.

    Srinivasan A, Saxena NP, Johansen C. Cold tolerance during early reproductive growth of chickpea (Cicer arietinum L.): genetic variation in gamete development and function. Field Crop Res. 1999;60(3):209–22.

    Article  Google Scholar 

  9. 9.

    Gao YB, Wang CL, Wu JY, Zhou HS, Jiang XT, Wu J, Zhang SL. Low temperature inhibits pollen tube growth by disruption of both tip-localized reactive oxygen species and endocytosis in Pyrus bretschneideri Rehd. Plant Physiol Bioch. 2014;74:255–62.

    CAS  Article  Google Scholar 

  10. 10.

    Sekiya J, Yamashita K, Nakagawa S, Shibata Y, Hatanaka A. Phospholipids of Tea Pollen. Agr Biol Chem Tokyo. 1988;52(1):243–7.

    CAS  Google Scholar 

  11. 11.

    Wang YH, Li XC, Zhu-Ge Q, Jiang X, Wang WD, Fang WP, Chen X, Li XH. Nitric oxide participates in cold-inhibited Camellia sinensis pollen germination and tube growth partly via cGMP in vitro. PLoS One. 2012;7(12):e52436.

  12. 12.

    Arasimowicz M, Floryszak-Wieczorek J. Nitric oxide as a bioactive signalling molecule in plant stress responses. Plant Sci. 2007;172(5):876–87.

    CAS  Article  Google Scholar 

  13. 13.

    Bai XG, Yang LM, Yang YQ, Ahmad P, Yang YP, Hu XY. Deciphering the protective role of nitric oxide against salt stress at the physiological and proteomic levels in Maize. J Proteome Res. 2011;10(10):4349–64.

    CAS  PubMed  Article  Google Scholar 

  14. 14.

    Arasimowicz-Jelonek M, Floryszak-Wieczorek J, Kubis J. Interaction between polyamine and nitric oxide signaling in adaptive responses to drought in cucumber. J Plant Growth Regul. 2009;28(2):177–86.

    CAS  Article  Google Scholar 

  15. 15.

    Bouchard JN, Yamasaki H. Heat stress stimulates nitric oxide production in Symbiodinium microadriaticum: a possible linkage between nitric oxide and the coral bleaching phenomenon. Plant Cell Physiol. 2008;49(4):641–52.

    CAS  PubMed  Article  Google Scholar 

  16. 16.

    Cui XM, Zhang YK, Chen XL, Jin H, Wu XB. Effects of exogenous nitric oxide protects tomato plants under copper stress, 2009 3rd International Conference on Bioinformatics and Biomedical Engineering, Vols 1–11. 2009. p. 4417–23.

    Google Scholar 

  17. 17.

    Wang YB, Feng HY, Qu Y, Cheng JQ, Zhao ZG, Zhang MX, Wang XL, An LZ. The relationship between reactive oxygen species and nitric oxide in ultraviolet-B-induced ethylene production in leaves of maize seedlings. Environ Exp Bot. 2006;57(1–2):51–61.

    CAS  Article  Google Scholar 

  18. 18.

    Zhao MG, Chen L, Zhang LL, Zhang WH. Nitric reductase-dependent nitric oxide production is involved in cold acclimation and freezing tolerance in Arabidopsis. Plant Physiol. 2009;151(2):755–67.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  19. 19.

    Majlath I, Szalai G, Soos V, Sebestyen E, Balazs E, Vankova R, Dobrev PI, Tari I, Tandori J, Janda T. Effect of light on the gene expression and hormonal status of winter and spring wheat plants during cold hardening. Physiol Plantarum. 2012;145(2):296–314.

    CAS  Article  Google Scholar 

  20. 20.

    Airaki M, Leterrier M, Mateos RM, Valderrama R, Chaki M, Barroso JB, Del Rio LA, Palma JM, Corpas FJ. Metabolism of reactive oxygen species and reactive nitrogen species in pepper (Capsicum annuum L.) plants under low temperature stress. Plant Cell Environ. 2012;35(2SI):281–95.

    CAS  PubMed  Article  Google Scholar 

  21. 21.

    Guo Z, Tan J, Zhuo C, Wang C, Xiang B, Wang Z. Abscisic acid, H2O2 and nitric oxide interactions mediated cold-induced S-adenosylmethionine synthetase in Medicago sativa subsp falcata that confers cold tolerance through up-regulating polyamine oxidation. Plant Biotechnol J. 2014;12(5):601–12.

    CAS  PubMed  Article  Google Scholar 

  22. 22.

    Cui JX, Zhou YH, Ding JG, Xia XJ, Shi K, Chen SC, Asami T, Chen ZX, Yu JQ. Role of nitric oxide in hydrogen peroxide-dependent induction of abiotic stress tolerance by brassinosteroids in cucumber. Plant Cell Environ. 2011;34(2):347–58.

    CAS  PubMed  Article  Google Scholar 

  23. 23.

    Cantrel C, Vazquez T, Puyaubert J, Reze N, Lesch M, Kaiser WM, Dutilleul C, Guillas I, Zachowski A, Baudouin E. Nitric oxide participates in cold-responsive phosphosphingolipid formation and gene expression in Arabidopsis thaliana. New Phytol. 2011;189(2):415–27.

    CAS  PubMed  Article  Google Scholar 

  24. 24.

    Zhao RR, Sheng JP, Lv SN, Zheng Y, Zhang J, Yu MM, Shen L. Nitric oxide participates in the regulation of LeCBF1 gene expression and improves cold tolerance in harvested tomato fruit. Postharvest Biol Tec. 2011;62(2):121–6.

    CAS  Article  Google Scholar 

  25. 25.

    Sehrawat A, Gupta R, Deswal R. Nitric oxide-cold stress signalling cross-talk, evolution of a novel regulatory mechanism. Proteomics. 2013;13(12–13):1816–35.

    CAS  PubMed  Article  Google Scholar 

  26. 26.

    Prado AM, Colaco R, Moreno N, Silva AC, Feijo JA. Targeting of pollen tubes to ovules is dependent on nitric oxide (NO) signaling. Mol Plant. 2008;1(4):703–14.

    CAS  PubMed  Article  Google Scholar 

  27. 27.

    Prado AM, Porterfield DM, Feijo JA. Nitric oxide is involved in growth regulation and re-orientation of pollen tubes. Development. 2004;131(11):2707–14.

    CAS  PubMed  Article  Google Scholar 

  28. 28.

    Reichler SA, Torres J, Rivera AL, Cintolesi VA, Clark G, Roux SJ. Intersection of two signalling pathways: extracellular nucleotides regulate pollen germination and pollen tube growth via nitric oxide. J Exp Bot. 2009;60(7):2129–38.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  29. 29.

    Wang YH, Chen T, Zhang CY, Hao HQ, Liu P, Zheng MZ, Baluska F, Samaj J, Lin JX. Nitric oxide modulates the influx of extracellular Ca2+ and actin filament organization during cell wall construction in Pinus bungeana pollen tubes. New Phytol. 2009;182(4):851–62.

    CAS  PubMed  Article  Google Scholar 

  30. 30.

    Wang Z, Gerstein M, Snyder M. RNA-Seq: a revolutionary tool for transcriptomics. Nat Rev Genet. 2009;10(1):57–63.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  31. 31.

    Zhao YQ, Li J, Liu HJ, Xi Y, Xue M, Liu WH, Zhuang ZH, Lei MG. Dynamic transcriptome profiles of skeletal muscle tissue across 11 developmental stages for both Tongcheng and Yorkshire pigs. BMC Genomics. 2015;16.

  32. 32.

    Conesa A, Gotz S, Garcia-Gomez JM, Terol J, Talon M, Robles M. Blast2GO: a universal tool for annotation, visualization and analysis in functional genomics research. Bioinformatics. 2005;21(18):3674–6.

    CAS  PubMed  Article  Google Scholar 

  33. 33.

    Ye J, Fang L, Zheng HK, Zhang Y, Chen J, Zhang ZJ, Wang J, Li ST, Li RQ, Bolund L, et al. WEGO: a web tool for plotting GO annotations. Nucleic Acids Res. 2006;34(SI):W293–7.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  34. 34.

    Yang N, Yue XL, Chen XL, Wu GF, Zhang TG, An LZ. Molecular cloning and partial characterization of a novel phospholipase D gene from Chorispora bungeana. Plant Cell Tiss Org. 2012;108(2):201–12.

    CAS  Article  Google Scholar 

  35. 35.

    Yang N, Yue XL, Zhang H, Wu GF, Ding FX, Zhang TG, An LZ. Characterization of phospholipase D from Chorispora bungeana callus in response to freezing stress. Biologia Plantarum. 2013;57(1):113–20.

    CAS  Article  Google Scholar 

  36. 36.

    Druege U, Franken P, Hajirezaei MR. Plant hormone homeostasis, signaling, and function during adventitious root formation in cuttings. Front Plant Sci. 2016;7:381.

  37. 37.

    Kim SA, Ahn SY, Yun HK. Transcriptome analysis of grapevine shoots exposed to chilling temperature for four weeks. Hortic Environ Biote. 2016;57(2):161–72.

    Article  Google Scholar 

  38. 38.

    Wang F, Guo ZX, Li HZ, Wang MM, Onac E, Zhou J, Xia XJ, Shi K, Yu JQ, Zhou YH. Phytochrome A and B function antagonistically to regulate cold tolerance via abscisic acid-dependent Jasmonate signaling. Plant Physiol. 2016;170(1):459–71.

    CAS  PubMed  Article  Google Scholar 

  39. 39.

    Janowiak F, Maas B, Dorffling K. Importance of abscisic acid for chilling tolerance of maize seedlings. J Plant Physiol. 2002;159(6):635–43.

    CAS  Article  Google Scholar 

  40. 40.

    Mata CG, Lamattina L. Nitric oxide induces stomatal closure and enhances the adaptive plant responses against drought stress. Plant Physiol. 2001;126(3):1196–204.

    CAS  PubMed Central  Article  Google Scholar 

  41. 41.

    Neill SJ, Desikan R, Clarke A, Hancock JT. Nitric oxide is a novel component of abscisic acid signaling in stomatal guard cells. Plant Physiol. 2002;128(1):13–6.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  42. 42.

    Liu Y, Zhang J. Rapid accumulation of NO regulates ABA catabolism and seed dormancy during imbibition in Arabidopsis. Plant Signal Behav. 2009;4(9):905–7.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  43. 43.

    Wang P, Zhu JK, Lang Z. Nitric oxide suppresses the inhibitory effect of abscisic acid on seed germination by S-nitrosylation of SnRK2 proteins. Plant Signal Behav. 2015;10(6), e1031939.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  44. 44.

    Hauser F, Waadt R, Schroeder JI. Evolution of abscisic acid synthesis and signaling mechanisms. Curr Biol. 2011;21(9):R346–55.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  45. 45.

    Liu YJ, Zhao ZG, Si J, Di CX, Han J, An LZ. Brassinosteroids alleviate chilling-induced oxidative damage by enhancing antioxidant defense system in suspension cultured cells of Chorispora bungeana. Plant Growth Regul. 2009;59(3):207–14.

    CAS  Article  Google Scholar 

  46. 46.

    She J, Han ZF, Kim TW, Wang JJ, Cheng W, Chang JB, Shi SA, Wang JW, Yang MJ, Wang ZY, et al. Structural insight into brassinosteroid perception by BRI1. Nature. 2011;474(7352):472–U496.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  47. 47.

    Russinova E, Borst JW, Kwaaitaal M, Cano-Delgado A, Yin Y, Chory J, de Vries SC. Heterodimerization and endocytosis of Arabidopsis brassinosteroid receptors BRI1 and AtSERK3 (BAK1). Plant Cell. 2004;16(12):3216–29.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  48. 48.

    Tossi V, Lamattina L, Cassia R. Pharmacological and genetical evidence supporting nitric oxide requirement for 2,4-epibrassinolide regulation of root architecture in Arabidopsis thaliana. Plant Signal Behav. 2013;8(7), e24712.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  49. 49.

    Divi UK, Krishna P. Overexpression of the Brassinosteroid biosynthetic gene AtDWF4 in Arabidopsis seeds overcomes abscisic acid-induced inhibition of germination and increases cold tolerance in transgenic seedlings. J Plant Growth Regul. 2010;29(4):385–93.

    CAS  Article  Google Scholar 

  50. 50.

    Zhang AY, Zhang J, Zhang JH, Ye NH, Zhang H, Tan MP, Jiang MY. Nitric oxide mediates brassinosteroid-induced ABA biosynthesis involved in oxidative stress tolerance in maize leaves. Plant Cell Physiol. 2011;52(1):181–92.

    CAS  PubMed  Article  Google Scholar 

  51. 51.

    Ciardi JA, Deikman J, Orzolek MD. Increased ethylene synthesis enhances chilling tolerance in tomato. Physiol Plantarum. 1997;101(2):333–40.

    CAS  Article  Google Scholar 

  52. 52.

    Zhang Z, Huang R. Enhanced tolerance to freezing in tobacco and tomato overexpressing transcription factor TERF2/LeERF2 is modulated by ethylene biosynthesis. Plant Mol Biol. 2010;73(3):241–9.

    CAS  PubMed  Article  Google Scholar 

  53. 53.

    Shi YT, Tian SW, Hou LY, Huang XZ, Zhang XY, Guo HW, Yang SH. Ethylene signaling negatively regulates freezing tolerance by repressing expression of CBF and type-A ARR genes in Arabidopsis. Plant Cell. 2012;24(6):2578–95.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  54. 54.

    Hua J, Chang C, Sun Q, Meyerowitz EM. Ethylene insensitivity conferred by Arabidopsis ERS gene. Scince. 1995;269(5231):1712–4.

    CAS  Article  Google Scholar 

  55. 55.

    Roman G, Lubarsky B, Kieber JJ, Rothenberg M, Ecker JR. Genetic-analysis of ethylene signal-transduction in Arabidopsis-Thaliana - 5 novel mutant loci integrated into a stress-response pathway. Genetics. 1995;139(3):1393–409.

    CAS  PubMed  PubMed Central  Google Scholar 

  56. 56.

    Chao QM, Rothenberg M, Solano R, Roman G, Terzaghi W, Ecker JR. Activation of the ethylene gas response pathway in Arabidopsis by the nuclear protein ETHYLENE-INSENSITIVE3 and related proteins. Cell. 1997;89(7):1133–44.

    CAS  PubMed  Article  Google Scholar 

  57. 57.

    Potuschak T, Lechner E, Parmentier Y, Yanagisawa S, Grava S, Koncz C, Genschik P. EIN3-dependent regulation of plant ethylene hormone signaling by two Arabidopsis F box proteins: EBF1 and EBF2. Cell. 2003;115(6):679–89.

    CAS  PubMed  Article  Google Scholar 

  58. 58.

    Konishi M, Yanagisawa S. Ethylene signaling in Arabidopsis involves feedback regulation via the elaborate control of EBF2 expression by EIN3. Plant J. 2008;55(5):821–31.

    CAS  PubMed  Article  Google Scholar 

  59. 59.

    Zou Y, Zhang L, Rao S, Zhu XY, Ye LL, Chen WX, Li XP. The relationship between the expression of ethylene-related genes and papaya fruit ripening disorder caused by chilling injury. Plos One. 2014;9(12), e116002.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  60. 60.

    Sairam RK, Dharmar K, Chinnusamy V, Lekshmy S, Joshi R, Bhattacharya P. The role of non-symbiotic haemoglobin and nitric oxide homeostasis in waterlogging tolerance in Vigna species. Biologia Plantarum. 2012;56(3):528–36.

    CAS  Article  Google Scholar 

  61. 61.

    Li XL, Pan YJ, Chang BW, Wang YC, Tang ZH. NO promotes seed germination and seedling growth under high salt may depend on EIN3 protein in Arabidopsis. Front Plant Sci. 2016;6:1203.

  62. 62.

    He JM, Zhang Z, Wang RB, Chen YP. UV-B-induced stomatal closure occurs via ethylene-dependent NO generation in Vicia faba. Funct Plant Biol. 2011;38(4):293–302.

    CAS  Article  Google Scholar 

  63. 63.

    Mur LAJ, Laarhoven LJJ, Harren FJM, Hall MA, Smith AR. Nitric oxide interacts with salicylate to regulate biphasic ethylene production during the hypersensitive response. Plant Physiol. 2008;148(3):1537–46.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  64. 64.

    Jain M, Khurana JP. Transcript profiling reveals diverse roles of auxin-responsive genes during reproductive development and abiotic stress in rice. Febs J. 2009;276(11):3148–62.

    CAS  PubMed  Article  Google Scholar 

  65. 65.

    Rahman A. Auxin: a regulator of cold stress response. Physiol Plantarum. 2013;147(1):28–35.

    CAS  Article  Google Scholar 

  66. 66.

    Shibasaki K, Uemura M, Tsurumi S, Rahman A. Auxin response in Arabidopsis under cold stress: underlying molecular mechanisms. Plant Cell. 2009;21(12):3823–38.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  67. 67.

    Kovaleva LV, Voronkov AS, Zakharova EV. Role of auxin and cytokinin in the regulation of the actin cytoskeleton in the in vitro germinating male gametophyte of petunia. Russ J Plant Physl+. 2015;62(2):179–86.

    CAS  Article  Google Scholar 

  68. 68.

    Schlicht M, Ludwig-Muller J, Burbach C, Volkmann D, Baluska F. Indole-3-butyric acid induces lateral root formation via peroxisome-derived indole-3-acetic acid and nitric oxide. New Phytol. 2013;200(2):473–82.

    CAS  PubMed  Article  Google Scholar 

  69. 69.

    Achard P, Gong F, Cheminant S, Alioua M, Hedden P, Genschik P. The cold-inducible CBF1 factor-dependent signaling pathway modulates the accumulation of the growth-repressing DELLA proteins via its effect on gibberellin metabolism. Plant Cell. 2008;20(8):2117–29.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  70. 70.

    Colebrook EH, Thomas SG, Phillips AL, Hedden P. The role of gibberellin signalling in plant responses to abiotic stress. J Exp Biol. 2014;217(Pt 1):67–75.

    CAS  PubMed  Article  Google Scholar 

  71. 71.

    Thomas SG, Phillips AL, Hedden P. Molecular cloning and functional expression of gibberellin 2-oxidases, multifunctional enzymes involved in gibberellin deactivation. Proc Natl Acad Sci U S A. 1999;96(8):4698–703.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  72. 72.

    Wang YJ, Deng DX. Molecular basis and evolutionary pattern of GA-GID1-DELLA regulatory module. Mol Genet Genomics. 2014;289(1):1–9.

    CAS  PubMed  Article  Google Scholar 

  73. 73.

    Fernandez-Marcos M, Sanz L, Lorenzo O. Nitric oxide: an emerging regulator of cell elongation during primary root growth. Plant Signal Behav. 2012;7(2):196–200.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  74. 74.

    Yang A, Dai X, Zhang WH. A R2R3-type MYB gene, OsMYB2, is involved in salt, cold, and dehydration tolerance in rice. J Exp Bot. 2012;63(7):2541–56.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  75. 75.

    Yu GH, Jiang LL, Ma XF, Xu ZS, Liu MM, Shan SG, Cheng XG. A soybean C2H2-type zinc finger gene GmZF1 enhanced cold tolerance in transgenic Arabidopsis. Plos One. 2014;9(10), e109399.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  76. 76.

    Serpa V, Vernal J, Lamattina L, Grotewold E, Cassia R, Terenzi H. Inhibition of AtMYB2 DNA-binding by nitric oxide involves cysteine S-nitrosylation. Biochem Bioph Res Co. 2007;361(4):1048–53.

    CAS  Article  Google Scholar 

  77. 77.

    Kroncke KD. Zinc finger proteins as molecular targets for nitric oxide-mediated gene regulation. Antioxid Redox Signal. 2001;3(4):565–75.

    CAS  PubMed  Article  Google Scholar 

  78. 78.

    Verelst W, Twell D, de Folter S, Immink R, Saedler H, Munster T. MADS-complexes regulate transcriptome dynamics during pollen maturation. Genome Biol. 2007;8(11):R249.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  79. 79.

    Jiang XT, Gao YB, Zhou HS, Chen JQ, Wu JY, Zhang SL. Apoplastic calmodulin promotes self-incompatibility pollen tube growth by enhancing calcium influx and reactive oxygen species concentration in Pyrus pyrifolia. Plant Cell Rep. 2014;33(2):255–63.

    CAS  PubMed  Article  Google Scholar 

  80. 80.

    Chen T, Wu X, Chen Y, Li X, Huang M, Zheng M, Baluska F, Samaj J, Lin J. Combined proteomic and cytological analysis of Ca2+-calmodulin regulation in Picea meyeri pollen tube growth. Plant Physiol. 2009;149(2):1111–26.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  81. 81.

    Rato C, Monteiro D, Hepler PK, Malho R. Calmodulin activity and cAMP signalling modulate growth and apical secretion in pollen tubes. Plant J. 2004;38(6):887–97.

    CAS  PubMed  Article  Google Scholar 

  82. 82.

    Yang TB, Peng H, Whitaker BD, Conway WS. Characterization of a calcium/calmodulin-regulated SR/CAMTA gene family during tomato fruit development and ripening. BMC Plant Biol. 2012;12:19.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  83. 83.

    Saha G, Park JI, Jung HJ, Ahmed NU, Kayum MA, Chung MY, Hur Y, Cho YG, Watanabe M, Nou IS. Genome-wide identification and characterization of MADS-box family genes related to organ development and stress resistance in Brassica rapa. BMC Genomics. 2015;16:178.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  84. 84.

    Doherty CJ, Van Buskirk HA, Myers SJ, Thomashow MF. Roles for Arabidopsis CAMTA transcription factors in cold-regulated gene expression and freezing tolerance. Plant Cell. 2009;21(3):972–84.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  85. 85.

    Crawford NM. Mechanisms for nitric oxide synthesis in plants. J Exp Bot. 2006;57(3):471–8.

    CAS  PubMed  Article  Google Scholar 

  86. 86.

    Jamai A, Dubois E, Vershon AK, Messenguy F. Swapping functional specificity of a MADS box protein: residues required for Arg80 regulation of arginine metabolism. Mol Cell Biol. 2002;22(16):5741–52.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  87. 87.

    Zou C, Jiang W, Yu D. Male gametophyte-specific WRKY34 transcription factor mediates cold sensitivity of mature pollen in Arabidopsis. J Exp Bot. 2010;61(14):3901–14.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  88. 88.

    Lee CM, Thomashow MF. Photoperiodic regulation of the C-repeat binding factor (CBF) cold acclimation pathway and freezing tolerance in Arabidopsis thaliana. Proc Natl Acad Sci U S A. 2012;109(37):15054–9.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  89. 89.

    Moscatelli A, Idilli AI. Pollen tube growth: a delicate equilibrium between secretory and endocytic pathways. J Integr Plant Biol. 2009;51(8):727–39.

    CAS  PubMed  Article  Google Scholar 

  90. 90.

    Kang B, Nielsen E, Preuss ML, Mastronarde D, Staehelin LA. Electron tomography of RabA4b-and PI-4 K beta 1-labeled trans Golgi network compartments in Arabidopsis. Traffic. 2011;12(3):313–29.

    CAS  PubMed  Article  Google Scholar 

  91. 91.

    Weernink PAO, Schmidt M, Jakobs KH. Regulation and cellular roles of phosphoinositide 5-kinases. Eur J Pharmacol. 2004;500(1–3):87–99.

    Article  CAS  Google Scholar 

  92. 92.

    Guo F, McCubbin AG. The pollen-specific R-SNARE/longin PiVAMP726 mediates fusion of endo- and exocytic compartments in pollen tube tip growth. J Exp Bot. 2012;63(8):3083–95.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  93. 93.

    Deng SR, Sun J, Zhao R, Ding MQ, Zhang YN, Sun YL, Wang W, Tan YQ, Liu DD, Ma XJ, et al. Populus euphratica APYRASE2 enhances cold tolerance by modulating vesicular trafficking and extracellular ATP in Arabidopsis Plants. Plant Physiol. 2015;169(1):530. +.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  94. 94.

    Lombardo MC, Lamattina L. Nitric oxide is essential for vesicle formation and trafficking in Arabidopsis root hair growth. J Exp Bot. 2012;63(13):4875–85.

    CAS  PubMed  Article  Google Scholar 

  95. 95.

    Chen T, Teng N, Wu X, Wang Y, Tang W, Samaj J, Baluska F, Lin J. Disruption of actin filaments by latrunculin B affects cell wall construction in Picea meyeri pollen tube by disturbing vesicle trafficking. Plant Cell Physiol. 2007;48(1):19–30.

    CAS  PubMed  Article  Google Scholar 

  96. 96.

    Schindelman G, Morikami A, Jung J, Baskin TI, Carpita NC, Derbyshire P, McCann MC, Benfey PN. COBRA encodes a putative GPI-anchored protein, which is polarly localized and necessary for oriented cell expansion in Arabidopsis. Genes Dev. 2001;15(9):1115–27.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  97. 97.

    Vierstra RD. The ubiquitin-26S proteasome system at the nexus of plant biology. Nat Rev Mol Cell Biol. 2009;10(6):385–97.

    CAS  PubMed  Article  Google Scholar 

  98. 98.

    Lee JH, Kim WT. Regulation of abiotic stress signal transduction by E3 ubiquitin ligases in Arabidopsis. Mol Cells. 2011;31(3):201–8.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  99. 99.

    Yee D, Goring DR. The diversity of plant U-box E3 ubiquitin ligases: from upstream activators to downstream target substrates. J Exp Bot. 2009;60(4):1109–21.

    CAS  PubMed  Article  Google Scholar 

  100. 100.

    Lechner E, Achard P, Vansiri A, Potuschak T, Genschik P. F-box proteins everywhere. Curr Opin Plant Biol. 2006;9(6):631–8.

    CAS  PubMed  Article  Google Scholar 

  101. 101.

    Gingerich DJ, Gagne JM, Salter DW, Hellmann H, Estelle M, Ma L, Vierstra RD. Cullins 3a and 3b assemble with members of the broad complex/tramtrack/bric-a-brac (BTB) protein family to form essential ubiquitin-protein ligases (E3s) in Arabidopsis. J Biol Chem. 2005;280(19):18810–21.

    CAS  PubMed  Article  Google Scholar 

  102. 102.

    Lyzenga WJ, Stone SL. Abiotic stress tolerance mediated by protein ubiquitination. J Exp Bot. 2012;63(2):599–616.

    CAS  PubMed  Article  Google Scholar 

  103. 103.

    Yu F, Wu Y, Xie Q. Ubiquitin-proteasome system in ABA signaling: from perception to action. Mol Plant. 2016;9(1):21–33.

    CAS  PubMed  Article  Google Scholar 

  104. 104.

    Kulik A, Noirot E, Grandperret V, Bourque S, Fromentin J, Salloignon P, Truntzer C, Dobrowolska G, Simon-Plas F, Wendehenne D. Interplays between nitric oxide and reactive oxygen species in cryptogein signalling. Plant Cell Environ. 2015;38(2):331–48.

    CAS  PubMed  Article  Google Scholar 

  105. 105.

    Dixon RA, Pasinetti GM. Flavonoids and isoflavonoids: from plant biology to agriculture and neuroscience. Plant Physiol. 2010;154(2):453–7.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  106. 106.

    Lillo C, Lea US, Ruoff P. Nutrient depletion as a key factor for manipulating gene expression and product formation in different branches of the flavonoid pathway. Plant Cell Environ. 2008;31(5):587–601.

    CAS  PubMed  Article  Google Scholar 

  107. 107.

    Ashihara H, Kubota H. Patterns of adenine metabolism and caffeine biosynthesis in different parts of tea seedlings. Physiol Plantarum. 1986;68(2):275–81.

    CAS  Article  Google Scholar 

  108. 108.

    Yamaguchi S, Ninomiya K. Umami and food palatability. J Nutr. 2000;130(4):921s–6.

    CAS  PubMed  Google Scholar 

  109. 109.

    Shi CY, Yang H, Wei CL, Yu O, Zhang ZZ, Jiang CJ, Sun J, Li YY, Chen Q, Xia T, et al. Deep sequencing of the Camellia sinensis transcriptome revealed candidate genes for major metabolic pathways of tea-specific compounds. BMC Genomics. 2011;12.

  110. 110.

    Li CF, Zhu Y, Yu Y, Zhao QY, Wang SJ, Wang XC, Yao MZ, Luo D, Li X, Chen L, et al. Global transcriptome and gene regulation network for secondary metabolite biosynthesis of tea plant (Camellia sinensis). BMC Genomics. 2015;16:560.

  111. 111.

    Wang W, Sheng X, Shu Z, Li D, Pan J, Ye X, Chang P, Li X, Wang Y. Combined cytological and transcriptomic analysis reveals a nitric oxide signaling pathway involved in cold-inhibited Camellia sinensis pollen tube growth. Front Plant Sci. 2016;7:456.

    PubMed  PubMed Central  Google Scholar 

  112. 112.

    Grabherr MG, Haas BJ, Yassour M, Levin JZ, Thompson DA, Amit I, Adiconis X, Fan L, Raychowdhury R, Zeng QD, et al. Full-length transcriptome assembly from RNA-Seq data without a reference genome. Nat Biotechnol. 2011;29(7):130–644.

    Article  CAS  Google Scholar 

  113. 113.

    Pertea G, Huang XQ, Liang F, Antonescu V, Sultana R, Karamycheva S, Lee Y, White J, Cheung F, Parvizi B, et al. TIGR gene indices clustering tools (TGICL): a software system for fast clustering of large EST datasets. Bioinformatics. 2003;19(5):651–2.

    CAS  PubMed  Article  Google Scholar 

  114. 114.

    Iseli C, Jongeneel CV, Bucher P. ESTScan: a program for detecting, evaluating, and reconstructing potential coding regions in EST sequences. In: Lengauer T, Schneider R, Bork P, Brutlag D, Glasgow J, Mewes HW, Zimmer R, editors. Proceedings of the seventh international conference on intelligent systems for molecular biology. Menlo Park: AAAI Press; 1999. p.138–148.

  115. 115.

    Tarazona S, Garcia-Alcalde F, Dopazo J, Ferrer A, Conesa A. Differential expression in RNA-seq: A matter of depth. Genome Res. 2011;21(12):2213–23.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  116. 116.

    de Hoon MJ, Imoto S, Nolan J, Miyano S. Open source clustering software. Bioinformatics. 2004;20(9):1453–4.

    PubMed  Article  CAS  Google Scholar 

  117. 117.

    Saldanha AJ. Java treeview--extensible visualization of microarray data. Bioinformatics. 2004;20(17):3246–8.

    CAS  PubMed  Article  Google Scholar 

Download references

Acknowledgements

Not applicable.

Funding

This work was supported by the National Natural Science Foundation of China (No. 31370014), the Opening Project of State Key Laboratory of Tea Plant Biology and Utilization (No. SKLTOF20150113) and the Fundamental Research Funds For the Central Universities (No. KYZ201509).

Availability of data and materials

RNA-Seq data used in this study have been deposited into the NCBI Sequence Read Archive (SRA, http://www.ncbi.nlm.nih.gov/sra/) under accession number SRR3270364, SRR3270376, SRR3270829, SRR3270928, SRR3270974, SRR3270993, SRR3270997, SRR3271001 and SRR3271002.

Authors’ contributions

JP, YW and WW conceived the study. JP, WW carried out bioinformatics data analysis. DL, ZS and XY performed qPCR analysis and participated in the sequence alignment. JP and YW wrote the manuscript. YW, WW and PC critically reviewed the manuscript. All authors read and approved the final manuscript.

Competing interest

The authors declare that they have no competing interests.

Consent for publication

Not applicable.

Ethics approval and consent to participate

Not applicable.

Author information

Affiliations

  1. College of Horticulture, Nanjing Agricultural University, Nanjing, 210095, China

    Junting Pan, Weidong Wang, Dongqin Li, Zaifa Shu, Xiaoli Ye, Pinpin Chang & Yuhua Wang

Authors
  1. Junting Pan
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Weidong Wang
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Dongqin Li
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Zaifa Shu
    View author publications

    You can also search for this author in PubMed Google Scholar

  5. Xiaoli Ye
    View author publications

    You can also search for this author in PubMed Google Scholar

  6. Pinpin Chang
    View author publications

    You can also search for this author in PubMed Google Scholar

  7. Yuhua Wang
    View author publications

    You can also search for this author in PubMed Google Scholar

Corresponding author

Correspondence to Yuhua Wang.

Additional files

Additional file 1: Figure S1.

Assessment of assembly quality. Distribution of mapped reads within the assembled unigenes determined unigene assembly quality. (TIF 2384 kb)

Additional file 2: Table S1.

DEGs identified from the comparison between control (CsPT-CK) and 4 °C-treated (CsPT-LT) pollen tbues. All of the samples were replicated three times. CK and LT FPKM: fragments per kb per million reads for each unigene in the CK and LT libraries, respectively. The log2Ratio (LT/CK): ratio between the FPKM of LT and CK. The absolute values of log2Ratio > 1 and probability > 0.7 were used as threshold for assigning significance. Annotation of DEGs against NR, NT, Swiss-Prot protein, KEGG, COG and GO were all reported in the tables. “-”: no hit. (XLS 381 kb)

Additional file 3: Table S2.

DEGs identified from the comparison between control (CsPT-CK) and NO-treated (CsPT-NO) pollen tubes. All of the samples were replicated three times. CK and NO FPKM: fragments per kb per million reads for each unigene in the CK and NO libraries, respectively. The log2Ratio (NO/CK): ratio between the FPKM of NO and CK. The absolute values of log2Ratio > 1 and probability > 0.7 were used as threshold for assigning significance. Annotation of DEGs against NR, NT, Swiss-Prot protein, KEGG, COG and GO were all reported in the tables. “-”: no hit. (XLS 250 kb)

Additional file 4: Table S3.

DEGs identified from the comparison between 4 °C- (CsPT-LT) and NO-treated (CsPT-NO) pollen tubes. All of the samples were replicated three times. LT and NO FPKM: fragments per kb per million reads for each unigene in the LT and NO libraries, respectively. The log2Ratio (LT/NO): ratio between the FPKM of LT and CK. The absolute values of log2Ratio > 1 and probability > 0.7 were used as threshold for assigning significance. Annotation of DEGs against NR, NT, Swiss-Prot protein, KEGG, COG and GO were all reported in the tables. “-”: no hit. (XLS 451 kb)

Additional file 5: Table S4.

GO classification of DEGs identified from the comparison between control (CK) and 4 °C-treated (LT) pollen tubes. (XLS 52 kb)

Additional file 6: Table S5.

GO classification of DEGs identified from the comparison between control (CK) and NO-treated (NO) pollen tubes. (XLS 40 kb)

Additional file 7: Table S6.

GO classification of DEGs identified from the comparison between pollen tubes exposed to 4 °C (LT) and pollen tubes treated with NO donor (NO). (XLS 59 kb)

Additional file 8: Table S7.

The 125 DEGs identified from the comparison between CK-VS-LT and CK-VS-NO. (XLS 89 kb)

Additional file 9: Table S8.

DEGs involved in vesicle polarized trafficking and cell wall biosynthesis between CK and LT (CK-VS-LT). The absolute values of log2Ratio (LT/CK) > 1 and probability > 0.7 were used as threshold for assigning significance. CK: control; LT: 4 °C treatment. (DOC 37 kb)

Additional file 10: Table S9.

DEGs involved in vesicle polarized trafficking and cell wall biosynthesis between CK and NO (CK-VS-NO). The absolute values of log2Ratio (NO/CK) > 1 and probability > 0.7 were used as threshold for assigning significance. CK: control; NO: NO treatment. (DOC 31 kb)

Additional file 11: Table S10.

DEGs involved in the ubiquitination machinery of the ubiquitin system between CK and LT (CK-VS-LT). The absolute values of log2Ratio (LT/CK) > 1 and probability > 0.7 were used as threshold for assigning significance. CK: control; LT: 4 °C treatment. (DOC 55 kb)

Additional file 12: Table S11.

DEGs involved in the ubiquitination machinery of the ubiquitin system between CK and NO (CK-VS-NO). The absolute values of log2Ratio (NO/CK) > 1 and probability > 0.7 were used as threshold for assigning significance. CK: control; NO: NO treatment. (DOC 37 kb)

Additional file 13: Table S12.

DEGs involved in flavonoids, caffeine and theanine biosynthesis pathways between CK and LT (CK-VS-LT). The absolute values of log2Ratio (LT/CK) > 1 and probability > 0.7 were used as threshold for assigning significance. CK: control; LT: 4 °C treatment. (DOC 36 kb)

Additional file 14: Table S13.

DEGs involved in flavonoids, caffeine and theanine biosynthesis pathways between CK and NO (CK-VS-NO). The absolute values of log2Ratio (NO/CK) > 1 and probability > 0.7 were used as threshold for assigning significance. CK: control; NO: NO treatment. (DOC 34 kb)

Additional file 15: Figure S2.

GO functional classifications in transcriptome data of C. sinensis leaves and pollen tubes under low temperature stress. A: GO functional classifications in transcriptome data of C. sinensis leaves under cold stress by Wang et al. [5]. B: GO functional classifications in transcriptome data of C. sinensis pollen tubes under cold stress. (TIF 3893 kb)

Additional file 16: Table S14.

Primers used in the quantitative real-time PCR (qRT-PCR) analysis of DEGs. (DOC 33 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Pan, J., Wang, W., Li, D. et al. Gene expression profile indicates involvement of NO in Camellia sinensis pollen tube growth at low temperature. BMC Genomics 17, 809 (2016). https://doi.org/10.1186/s12864-016-3158-4

Download citation

Keywords

  • Pollen tube
  • Low temperature
  • Nitric oxide
  • Camellia sinensis
  • RNA-Seq
  • Transcriptome

BMC Genomics

ISSN: 1471-2164

Contact us