Transcriptome Profiling of Petal Abscission Zone and Functional Analysis of an Aux/IAA Family Gene RhIAA16 Involved in Petal Shedding in Rose

Roses are one of the most important cut flowers among ornamental plants. Rose flower longevity is largely dependent on the timing of petal shedding occurrence. To understand the molecular mechanism underlying petal abscission in rose, we performed transcriptome profiling of the petal abscission zone during petal shedding using Illumina technology. We identified a total of 2592 differentially transcribed genes (DTGs) during rose petal shedding. Gene ontology term enrichment and pathway analysis revealed that major biochemical pathways the DTGs were involved in included ethylene biosynthesis, starch degradation, superpathway of cytosolic glycolysis, pyruvate dehydrogenase and TCA cycle, photorespiration and the lactose degradation III pathway. This suggests that alterations in carbon metabolism are an important part of rose petal abscission. Among these DTGs, approximately 150 genes putatively encoding transcription factors were identified in rose abscission zone. These included zinc finger, WRKY, ERF, and Aux/IAA gene families, suggesting that petal abscission involves complex transcriptional reprogramming. Approximately 108 DTGs were related to hormone pathways, of which auxin and ethylene related DTGs were the largest groups including 52 and 41 genes, respectively. These also included 12 DTGs related to gibberellin and 6 DTGs in jasmonic acid pathway. Surprisingly, no DTGs involved in the biosynthesis/signaling of abscisic acid, cytokinin, brassinosteroid, and salicylic acid pathways were detected. Moreover, among DTGs related to auxin, we identified an Aux/IAA gene RhIAA16 that was up-regulated in response to petal shedding. Down-regulation of RhIAA16 by virus-induced gene silencing in rose promoted petal abscission, suggesting that RhIAA16 plays an important role in rose petal abscission.


INTRODUCTION
Plant organ abscission is a crucial process that occurs throughout the life span of plants, and regulates the detachment of organs from main body (Roberts et al., 2002). This will benefit plants for recycling nutrients for continuous growth, propagating, facilitating reproduction, and preventing from disease infections (Addicott, 1982;. In particular, flower, fruit, and seed abscission is closely correlated with plant reproductive success, crop quality, and productivity (Addicott, 1982;Roberts et al., 2002;Estornell et al., 2013;. Initiation of flower organ abscission is triggered by both internal and external cues (Taylor and Whitelaw, 2001). As internal cues, interaction of auxin and ethylene plays a critical role in abscission initiation (Meir et al., 2010). Depletion of the polar flow of auxin passing through the abscission zone (AZ) makes the AZ sensitive to ethylene, which triggers the separation process (Abeles and Rubinstein, 1964;Addicott, 1982;Taylor and Whitelaw, 2001). Ethylene biosynthesis and signal transduction pathways are involved in the regulation of abscission. In Arabidopsis, ethylene-insensitive mutants etr1-1and ein2 inhibited flower organ shedding, suggesting the roles of ETR1 and EIN2 in abscission (Patterson and Bleecker, 2004). In tomato, repression of EIN3-like gene LeEIL expression retarded the flower pedicel abscission and fruit ripening processes (Tieman et al., 2001). Tomato never ripe (nr), sletr1-1, and sletr1-2 mutants affected ethylene receptor function and ethylene sensitivity, thereby delayed fruit ripening and organ abscission (Lanahan et al., 1994;Okabe et al., 2011). On the contrary, auxin as a negative regulator of abscission inhibits the cell separation process (Addicott, 1982;Estornell et al., 2013). The change in auxin flow results in the changes of transcript abundance of many genes involved in auxin biosynthesis, signal transduction, and transport. Functional studies of auxin response factors (ARFs) 1, 2, 7, and 19 demonstrated that these transcriptional regulators have functions in floral organ abscission (Ellis et al., 2005;Okushima et al., 2005). However, the roles of other gene families in the auxin pathway in the regulation of the petal abscission process are still largely unknown.
The perception and transduction of auxin signaling involve the cooperative action of several components. Among them, auxin/indole-3-acetic acid (Aux/IAA) proteins act as transcription repressors by dimerizing with auxin response factors (ARFs; Leyser, 2002;Woodward and Bartel, 2005). In presence of auxin, Aux/IAA proteins binding to the transport inhibitor response one/auxin signaling F-box (TIR1/AFB) cause degradation of Aux/IAA proteins, which then releases ARFs to trigger the expression of auxin responsive genes (Kepinski and Leyser, 2002;Woodward and Bartel, 2005). In Arabidopsis, analyses of Aux/IAA gain-of-function mutants revealed functional redundancy among Aux/IAA members (Overvoorde et al., 2005). However, distinctive expression patterns in different organs and tissues among Aux/IAA genes are displayed in several plant species such as rice (Jain et al., 2006), tomato (Audran-Delalande et al., 2012), Populus (Kalluri et al., 2007). In addition, functional analyses of Aux/IAA homologues in different plant species demonstrated the distinct and diverse roles of Aux/IAA proteins in plant development and growth (Audran-Delalande et al., 2012). In tomato, expression of an Aux/IAA family gene Sl-IAA3 is auxin-and ethylene-dependent. The phenotypes resulting from Sl-IAA3-silenced transgenic tomato suggested that Sl-IAA3 plays a role in interaction of auxin and ethylene in differential growth (Chaabouni et al., 2009a,b). However, the roles of Aux/IAAs in flower petal abscission are not well documented.
Roses are one of the most important cut flowers among ornamental plants. The opening and longevity of rose flower are major factors in determining the ornamental value of rose flower. Moreover, rose flower longevity is largely dependent on the timing of petal shedding occurrence. However, information on the molecular mechanism governing the rose petal abscission is scarce. To date, the regulatory genes in abscission signaling pathway, including IDA (Cho et al., 2008), NEVERSHED (Liljegren et al., 2009), EVERSHED (Leslie et al., 2010), were identified and characterized in model plants by genetic mapping of mutants. This forward genetic technique is difficult and time-consuming to identify and characterize the genes in nonmodel plant systems including rose. Next generation sequencing technology has become an effective method to investigate the regulatory network of abscission. Transcriptome studies of the flower abscission process were previously performed in several plant species including tomato (Meir et al., 2010;Wang et al., 2013;Ma et al., 2015a;Sundaresan et al., 2016), olive (Gil-Amado and Gomez-Jimenez, 2013), melon (Corbacho et al., 2013), apple (Botton et al., 2011;Zhu et al., 2011), and litchi . Most of these studies have focused on dissecting the regulatory mechanism of pedicel abscission that is the last phase of fruit development and ripening (Giovannoni, 2004). Combined with reverse genetic techniques, roles of several genes obtained from those transcriptome data related to abscission have been confirmed, including SlERF52 (Nakano et al., 2012(Nakano et al., , 2014 and KD1 (Ma et al., 2015a). In Arabidopsis, HAESA (HAE) and HAESA-LIKE2 (HSL2)-dependent pathways were revealed to be involved in petal AZ by the comparison of the transcriptomes of wild-type and hae hsl2 double mutant (Niederhuth et al., 2013). However, the mechanism regulating petal abscission in non-model plants is still not well understood.
Here we investigated the transcriptome dynamics of the petal AZ during petal shedding in rose by Illumina technology and dissected the transcriptional network governing the abscission process. Furthermore, we identified and characterized an Aux/IAA gene, RhIAA16, which we revealed to have an important role in controlling the timing of petal abscission in rose.

Plant Materials
Rose flowers (Rosa chinensis Jacq. cv. Gold Medal) were grown at a greenhouse on the campus of China Agricultural University, Beijing, China. Rose flower opening was divided into six stages: stage 1, partially opened bud; stage 2, completely opened bud; stages 3 and 4, partially opened flower; stage 5, fully opened flower with visible anthers; stage 6, fully opened flower with abscised petals (Figure 1). Rose flowers at different opening stages were harvested. The flower stems were placed immediately in water, and transported to the laboratory within 15 min. The flower stems were re-cut to 20 cm in length under water and placed in deionized water until further processing. The petal was shed at separation layer of abscission zone. Both distal and proximal sides of separation layer belong to abscission zone. Distal side attaches to petal organ, and proximal side attaches to receptacle (Figure 1). Therefore, we took sample of petal abscission zone by excising the base of petal (less than 1 mm in length) and the receptacle where petals locate (less than 1 mm in length). Since petals at stage 6 started to abscise, we only focused on the stages prior to petal shedding. Therefore, AZ samples at stages 1, 3, and 5 were collected and used for the transcriptome profiling. Three biological samples were collected for each stage.
Our preliminary tests demonstrated that R. hybrida cv. Samantha showed better responses to virus-induced gene silencing (VIGS) and much higher silencing efficiency than R. chinensis Jacq. cv. Gold Medal that was used for the transcriptome profiling (data not shown). In addition, the plantlets of R. hybrida cv. Samantha bloom as quickly as 40 days after rooting under our growth conditions. Therefore, rose plantlets of R. hybrida cv. Samantha were selected for VIGS. Rose plantlets were propagated by tissue culture. Rose shoots with at least 1 node and approximately 2 cm in length were used as explants and cultured on Murashige and Skoog (MS) medium supplemented with 1.0 mg/L 6-Benzyl aminopurine, 3 mg/L Gibberellic Acid, and 0.05 mg/L α-Naphthalene acetic acid for 30 days, then transferred to 1/2-strength MS medium supplemented with 0.1 mg/L NAA for 30 days for rooting.

Total RNA Extraction and RNA-Seq Library Preparation
Total RNA was extracted using the hot borate method according to previously described (Wan and Wilkins, 1994), and treated with RNase-free DNase I (Promega) to remove any contaminating genomic DNA. Three biological replicates were performed for each developmental stage (stages 1, 3, and 5). Strand-specific RNA libraries were constructed using the protocol described previously (Zhong et al., 2011), and sequenced on a HiSeq 2500 system according to the manufacturer's instructions. The raw reads were deposited into NCBI SRA database under accession no. PRJNA325324.

RNA-Seq Data Processing, De novo Assembly and Annotation
RNA-Seq reads were first processed to remove low quality and adaptor sequences using Trimmomatic (Bolger et al., 2014). Reads shorter than 40 bp were removed. The resulting RNA-Seq reads were then aligned to the ribosomal RNA (rRNA) database (Quast et al., 2013) using Bowtie with default parameters (Li and Durbin, 2010). Reads mapped to the rRNA database were discarded. The high-quality cleaned reads were assembled de novo into contigs using the Trinity program (Grabherr et al., 2011). To remove the redundancy of Trinity-assembled contigs, the contigs were further assembled de novo using iAssembler . The final assembled rose contigs were blasted against the UniProt (Swiss-Prot and TrEMBL; The UniProt Consortium, 2014) and Arabidopsis protein (version TAIR10) databases (Lamesch et al., 2012) with a cutoff e-value of 1e -5 . Based on FIGURE 1 | Flower opening stages in rose (Left), and petal abscission zone region used for RNA-Seq (Right). (Left) flower opening of rose was divided into six stages: stage 1, partially opened bud; stage 2, completely opened bud; stages 3 and 4, partially opened flower; stage 5, fully opened flower with visible anthers; stage 6, fully opened flower with abscised petals. (Right) abscission zone (AZ) sample used for RNA-Seq was the AZ at the base of petal (less than 1 mm in length) and AZ at the receptacle where petals locate (less than 1 mm in length).
the UniProt and Arabidopsis protein blast results, functional descriptions (human readable descriptions) were assigned to each unigene using AHRD 1 . Gene ontology (GO) terms were assigned to the rose assembled transcripts based on the GO terms annotated to their corresponding homologues in the UniProt database. Biochemical pathways were predicted from the rose transcripts using the Pathway Tools (Karp et al., 2002).
To identify differentially expressed genes, high-quality cleaned reads were aligned to rose contigs using Bowtie allowing up to two mismatches. Only the best alignments for each read were retained. Following alignments, raw read count for each rose contig in each sample was derived and normalized to reads per kilo base exon model per million mapped reads (RPKM). The significance of differential gene expression between different samples was determined using DESeq (Anders and Huber, 2010), and raw p-values of multiple tests were corrected using false discovery rate (Benjamini and Hochberg, 1995).

Quantitative RT-PCR
To analyze the transcript abundance of selected genes for confirmation of RNA-Seq data, quantitative RT-PCR (qRT-PCR) reactions were performed, as previously described (Ma et al., 2015b). Briefly, total RNAs were isolated from AZ samples with three biological repeats. 1 µl of the first strand cDNA was used as template with the Step One PlusTM real-time PCR system (Applied Biosystems) using KAPA TM SYBR R FAST quantitative PCR kits (Kapa Biosystems). RhActin5 was used as a reference gene (Pei et al., 2013). The primers used for determining transcript abundance were listed in Supplementary  Table S1.

Virus-Induced Gene Silencing
A 290 bp fragment in the 3 end of RhIAA16 gene was amplified by PCR from rose cDNAs to specifically silence RhIAA16. The fragment was then inserted into pTRV2 vector to generate the pTRV2-RhIAA16 construct. The primers used for amplifying RhIAA16 are listed in Supplementary Table S1.
Constructs of pTRV1, pTRV2, and pTRV2-RhIAA16 were transformed into Agrobacterium tumefaciens GV3101. A. tumefaciens were cultured in Luria-Bertani medium supplemented with 10 mM MES, 20 µM acetosyringone, 50 µg/ml gentamicin sulfate, and 50 µg/ml kanamycin. The cultured bacterium was collected by centrifuge at 4,000 rpm for 10 min, and re-suspended in infiltration buffer (10 mM MgCl 2 , 200 µM acetosyringone, 10 mM MES, pH 5.6) to OD 600 of ∼1.5 (Yin et al., 2015). Mixtures of cultures containing an equal ratio (v/v) of pTRV1 and pTRV2 or pTRV2-RhIAA16 were used for inoculation. VIGS of rose plantlets was performed as previously described (Tian et al., 2014). Rose plantlets, as shown in Supplementary Figure S1, were immersed in bacterial suspension solution and infiltrated under a vacuum at 0.7 MPa for 2 min. After infiltration, plantlets were washed in deionized water, and transplanted into pots 1 https://github.com/asishallab/AHRD-1 containing a mixture of 1:1 (v/v) of peat and vermiculite. The plantlets were immediately placed in dark at 8 • C for 3 days in a low temperature incubator (MIR-253, SANYO), and then grown in a culture room at 22 ± 1 • C, 40% relative humidity. Three independent experiments were performed. 30 plantlets were used for each experiment. Prior to the petal abscission test, we PCR-screened the plants and determined the transcript abundance of RhIAA16 in petals among the 30 plantlets. We found that the transcript levels of RhIAA16 in more than 30% plantlets were reduced, compared to uninoculated and empty vector controls. These plants with down-regulated RhIAA16 expression were used for petal abscission assay.

Sequencing and De novo Assembly of Petal Abscission Zone Transcriptome in Rose
To perform transcriptomic analysis, petal AZs of rose flowers (Figure 1) with three biological replicates at stage 1, 3, and 5 were collected to construct a total of nine RNA-Seq libraries. A total of 75,752,884 paired-end raw reads with length of 100 nucleotides (nt) were obtained. After further filtering and cleaning, a total of 57,312,389 clean read pairs were obtained. De novo assembly of these high-quality cleaned reads generated 80,226 unique transcripts with an average length of 743 bp ( Table 1). The size distribution indicated that the lengths of the 19,414 transcripts were more than 1000 bp (Supplementary Figure S2). Correlation coefficients of transcriptome profiles among the nine libraries and between the biological replicates were calculated (Supplementary Table S2). High correlation coefficients were obtained, suggesting the robustness of our RNA-Seq dataset.
To further validate the expression profiles of RNA-Seq data, four selected transcripts were analyzed by qRT-PCR. The results from qRT-PCR analysis were generally in agreement with the expression profiles obtained by RNA-Seq data (Figure 2).

Dynamic Transcriptome Profiles during Petal Abscission in Rose
Differentially transcribed genes (DTGs) were determined using a cutoff ratio of >2 or <0.5 when comparing their transcript  abundance in stage 3 to that in stage 1 (S3 vs. S1), and/or in stage 5 to that in stage 3 (S5 vs. S3). A total of 2592 DTGs were obtained (Supplementary Table S3). Based on the change in ratio of DTG transcript abundance, the number of DTGs at stage 3/1 and stage 5/3 was counted (Figure 3). Compared with stage 1, 782 DTGs were up-regulated and 300 DTGs were down-regulated in stage 3. Compared with stage 3, 1179 transcripts were increased and 1408 transcripts were decreased in stage 5 (Figure 3), suggesting that major transcriptional dynamic for petal abscission occurs just prior to petal shedding (stage 5).
To evaluate the potential functions of genes that showed transcriptional changes during petal abscission, we identified the GO terms of the DTGs in the biological process category (Figure 4). At stages 3 and 5, the biological processes were enriched in the metabolic process and defense responses including response to abiotic stimulus, external stimulus, organic substance (Figures 4A,B).
To identify the biochemical pathways involved in petal abscission, we analyzed DTGs using the Pathway Tools (Karp et al., 2002). The 108 DTGs at stage 3 and 261 DTGs at stage 5 were classified into 42 and 92 biochemical pathways, respectively (Supplementary Table S4). The major pathways both at stage 3 and 5 included ethylene biosynthesis, starch degradation, superpathway of cytosolic glycolysis, pyruvate dehydrogenase and TCA cycle, and photorespiration (Supplementary Table S4). In addition, one of the major pathways at stage 5 is related to the lactose degradation III pathway (Supplementary Table S4). These results suggested that alterations in carbon metabolism play an important part in rose petal abscission. FIGURE 4 | Gene ontology functional classification analysis of differentially transcribed genes (DTGs) during petal abscission. Histograms representing functional distributions (GO biological process class) of DTGs obtained from stage 3 compared to stage 1 (A), and from stage 5 compared to stage 3 (B). DTGs were determined using a cutoff ratio of >2 or <0.5 when comparing its expression in stage 3 to that in stage 1 (S3 vs. S1), and/or in stage 5 to that in stage 3 (S5 vs. S3).
Among DTGs related to auxin pathway, six auxin transporter genes were identified including five auxin efflux carrier genes ( Table 2). In addition, 14 Aux/IAA family members were obtained, of which six members were up-regulated in stage 5, and eight members were down-regulated in stage 5 compared to stage 3 ( Table 2). Among DTGs related to ethylene pathway, 15 DTGs encoded ethylene biosynthesis related 1-Aminocyclopropane-1-carboxylic acid oxidase (ACO). Transcript abundance of 10 ACO genes was accumulated in stage 5 compared to stage 3. Furthermore, 19 DTGs encoded ERFs were detected (Supplementary Table S6). The transcript abundance of 16 ERFs was increased in stage 5 compared to stage 3 ( Supplementary  Table S6). Overall, our results suggested that auxin and ethylene may play central roles in petal abscission of rose.

RhIAA16 Transcript Abundance Is Induced during Petal Abscission
To identify key regulators governing petal abscission, we initiated functional analysis of DTGs using VIGS. We primarily focused on the up-regulated DTGs, hypothesizing that VIGS-down regulation of these genes might lead to changes in the petal abscission processes. Given the potential important roles that auxin plays in the regulation of abscission, we first examined the functions of several Aux/IAA genes (RSA33069, RSA04500, RSA45030) that were up-regulated in the rose abscission zone using rose cut flowers. We found that VIGS-silencing of the contig of RSA33069 exhibited accelerated petal abscission phenotype. Analysis of the RSA33069 cDNA sequence revealed that it encoded a deduced protein of 253 amino acids with a 762 bp predicted open reading frame ( Figure 6A). The predicted amino acid sequence of RSA33069 showed that it belongs to the Aux/IAA family, and has the four canonical conserved domains known for this family (Figure 6A) (Overvoorde et al., 2005). In addition, phylogenetic tree analysis suggested that the protein has high degree of sequence homology to FvIAA16 from Fragaria vesca, therefore was designated as RhIAA16 ( Figure 6B). RT-PCR analysis demonstrated that the transcript abundance of RhIAA16 in petal AZ was significantly induced during flower development to peak at stage 5 ( Figure 7A). These results were consistent with the transcript abundance changes of RhIAA16 in the RNA-Seq data.

Reduction of RhIAA16 Expression Promotes Petal Abscission in Rose
To further confirm the potential role of RhIAA16 in petal abscission, we chose a fragment from RhIAA16-specific 3 untranslated region (UTR) to silence RhIAA16 in rose plantlets. qRT-PCR results showed that transcript abundance of RhIAA16 in RhIAA16-silenced (TRV2-RhIAA16) petal was significantly reduced compared to TRV2 control ( Figure 7B). Petal abscission was detected at 5 days after full opening in the RhIAA16-silenced plantlets in contrast to 9 days after full opening in TRV2 control plantlets (Figures 7C,D), suggesting that silencing of RhIAA16 accelerated the timing of initial petal abscission. At 9 days and 11 days after full opening, petals in ∼40.5 and 69.0% of flowers in RhIAA16-silenced plantlets had abscised whereas petals in only 17.8 and 37.9% of TRV2 control flowers were abscised, respectively ( Figure 7C).   The phenotypes of flowers were recorded and photographed every 2 days. The results were the means of three biological replicates with standard deviation. Letters indicated significant differences according to Duncan's multiple range test (P < 0.05), and asterisks indicated statistically significant differences (Student's t-test, P < 0.05).

DISCUSSION
In this study, high-throughput sequencing and de novo assembly strategies permitted us to dissect the transcriptome of rose petal AZ during petal shedding. Our results demonstrated that among DTGs related to hormones during petal abscission, most of them were associated with auxin and ethylene pathways, suggesting that auxin and ethylene play important roles in petal abscission in rose flowers ( Figure 5B). This conclusion is in good agreement with tomato abscission studies (Meir et al., 2010). Furthermore, functional characterization of RhIAA16 partially supports this notion (Figure 7). Auxin and ethylene as key hormones in the initiation of abscission have been demonstrated not only by physiological experiments, but also by transcriptome studies from different organ AZ. In tomato pedicel AZ, auxin depletion by auxin transport inhibitor, or flower removal stimulated pedicel abscission while ethylene action inhibitor treatment prevented the abscission induced by auxin depletion (Meir et al., 2010). In addition, the transcriptome of pedicel AZ demonstrated that acquisition of ethylene sensitivity in the AZ is associated with altered expression of auxin-regulated genes (Meir et al., 2010). Auxin homeostasis and signaling are usually modulated by the Aux/IAA genes (Song et al., 2009). Canonical Aux/IAA proteins function as transcriptional repressors of auxin-regulated genes (Tiwari et al., 2001(Tiwari et al., , 2004. In our study, the DTGs included 11 Aux/IAA family members ( Table 2). Among them, the transcript abundance of six members of Aux/IAA family was up-regulated during petal abscission ( Table 2). Intriguingly, VIGS-silencing RhIAA16, one of the up-regulated Aux/IAA genes, accelerated the petal abscission process (Figure 7), suggesting that RhIAA16 might be required for preventing premature abscission. In soybean, IAA16 has been reported as abscission-specific transcription factor by a transcriptome analysis of soybean leaf abscission, although the expression of IAA16 was down-regulated during soybean leaf abscission (Kim et al., 2016). In Arabidopsis, genetic evidences suggest that ARFs, which interact with Aux/IAA proteins, play regulatory roles in the petal abscission process. ARF1, ARF2, ARF7, and ARF19 were identified as regulators of abscission (Ellis et al., 2005). Given that ARF proteins interact with Aux/IAA in the auxin signal pathway (Leyser, 2002), further characterization of interactions between RhIAA16 and ARFs in rose may shed light on petal abscission activation.
It is worth pointing out that the transcript levels of RhIAA16 were not significantly changed in response to ethylene and ethylene action inhibitor 1-Methylcyclopropene (1-MCP) treatments (data not shown), suggesting that RhIAA16 is involved in either an ethylene independent pathway or up-stream of the ethylene pathway during abscission initiation.
In tomato, microarray assay showed that the expression of genes related to different steps of ethylene biosynthetic pathway, including S-adenosylmethionine (SAM) synthase, ACC synthase, and ACO genes, were altered during the pedicel abscission process (Meir et al., 2010). However, our transcriptome study revealed that multiple ACO genes in the final step of ethylene biosynthesis (Yang and Hoffman, 1984) were increased in AZ of stage 5 flower (Supplementary Table S6), indicating the critical role of ACO in controlling the rose petal abscission. Furthermore, the transcript abundance of many ERFs was induced in stage 5 compared with stage 3 (Supplementary Table S6). The importance of ERFs in flower abscission was recently demonstrated in tomato (Nakano et al., 2014). These researchers found that silencing of SlERF52, which is specifically expressed in the pedicel AZ, delayed tomato flower pedicel abscission.
In Arabidopsis, the transcriptomic analysis of petal abscission indicated that ethylene and ABA pathways were enriched during petal abscission (Niederhuth et al., 2013). In addition, JA signaling and biosynthesis genes were also involved in petal abscission (Niederhuth et al., 2013). Functional analysis showed that mutants of JA biosynthesis gene allene oxide synthase (AOS) retarded the petal abscission (Kim et al., 2013). This delayed abscission phenotype can be enhanced by ethylene insensitive mutant ein2 and ABA deficient mutant aba2, suggesting that ethylene, ABA and JA might synergistically regulate the petal abscission in Arabidopsis (Ogawa et al., 2009;Kim et al., 2013). In our study, the transcriptomic analysis suggested that among these three hormones, ethylene might play a major role, and JA might also be recruited in rose petal abscission, but not ABA (Figure 5). In GO analysis of our transcriptome showed that defense responses such as response to abiotic stimulus and external stimulus were enriched during petal abscission process (Figure 4), suggesting that the genes responsive to stress are involved in the activation of abscission. Indeed, abscission is considered to be a physiological process related to hormonemediated stress responses (Addicott, 1982;Estornell et al., 2013). In Arabidopsis, GO enrichment analysis of petal AZ transcriptome also demonstrated that biological processes that significantly enriched include defense response to abiotic and biotic stresses (Niederhuth et al., 2013). Similarly, in soybean, the biological processes of leaf AZ transcriptome significantly enriched include responses to endogenous and external stimuli (Kim et al., 2016).
Zinc finger transcription factors are a large and diverse family involved in many aspects of plant growth and development and play critical roles in cellular functions such as transcriptional regulation, RNA binding, and protein-protein interactions (Ciftci-Yilmaz and Mittler, 2008). Our results showed that zinc finger genes are the largest group among differentially transcribed transcription factors in rose petal AZ ( Figure 5A). In Arabidopsis, ZINC FINGER PROTEIN2 (ZFP2) encoding a ZFP, has been revealed in stamen AZ transcriptome profiling, which was increased during floral organ abscission process. Overexpression of ZFP2 exhibited delayed floral organ abscission phenotype (Cai and Lashbrook, 2008). Our data showed that 18 members of zinc finger gene were elevated in stage 5 (Supplementary Table S5), indicated the important roles of zinc finger gene in rose petal abscission. Functional analysis of these regulatory genes would be an important step toward elucidating their roles in petal abscission in the future.

CONCLUSION
Our RNA-Seq analysis has permitted us to dissect the rose AZ transcriptome during petal abscission, and reveal that auxin and ethylene are important hormones in the regulation of the abscission process. Furthermore, our data demonstrated that an Aux/IAA gene, RhIAA16, played an important role in rose petal abscission.

AUTHOR CONTRIBUTIONS
CM and C-ZJ conceived and supervised the study. YG, CL, XL, HX, and YL performed the experiments. YG and CL analyzed the data and drafted the manuscript. CM, NM, ZF, C-ZJ, and JG provided technical support, conceptual advice, analyzed the data, and participated in writing the manuscript. C-ZJ extensively revised the manuscript.

FUNDING
This work was supported by the National Science Foundation of China (Grants numbers 31520103913).

ACKNOWLEDGMENTS
We thank Dr. Bo Hong of China Agricultural University, China, for helpful comments.