Research article | Open | Published:
Proteomic insight into fruit set of cucumber (Cucumis sativus L.) suggests the cues of hormone-independent parthenocarpy
BMC Genomicsvolume 18, Article number: 896 (2017)
Parthenocarpy is an excellent agronomic trait that enables crops to set fruit in the absence of pollination and fertilization, and therefore to produce seedless fruit. Although parthenocarpy is widely recognized as a hormone-dependent process, hormone-insensitive parthenocarpy can also be observed in cucumber; however, its mechanism is poorly understood. To improve the global understanding of parthenocarpy and address the hormone-insensitive parthenocarpy shown in cucumber, we conducted a physiological and proteomic analysis of differently developed fruits.
Physiological analysis indicated that the natural hormone-insensitive parthenocarpy of ‘EC1’ has broad hormone-inhibitor resistance, and the endogenous hormones in the natural parthenocarpy (NP) fruits were stable and relatively lower than those of the non-parthenocarpic cultivar ‘8419 s-1.’ Based on the iTRAQ technique, 683 fruit developmental proteins were identified from NP, cytokinin-induced parthenocarpic (CP), pollinated and unpollinated fruits. Gene Ontology (GO) analysis showed that proteins detected from both set and aborted fruits were involved in similar biological processes, such as cell growth, the cell cycle, cell death and communication. Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis revealed that ‘protein synthesis’ was the major biological process that differed between fruit set and fruit abortion. Clustering analysis revealed that different protein expression patterns were involved in CP and NP fruits. Forty-one parthenocarpy-specialized DEPs (differentially expressed proteins) were screened and divided into two distinctive groups: NP-specialized proteins and CP-specialized proteins. Furthermore, qRT-PCR and western blot analysis indicated that NP-specialized proteins showed hormone- or hormone-inhibitor insensitive expression patterns in both ovaries and seedlings.
In this study, the global molecular regulation of fruit development in cucumber was revealed at the protein level. Physiological and proteomic comparisons indicated the presence of hormone-independent parthenocarpy and suppression of fruit abortion in cucumber. The proteomic analysis suggested that hormone-independent parthenocarpy is regulated by hormone-insensitive proteins such as the NP-specialized proteins. Moreover, the regulation of fruit abortion suppression may be closely related to protein synthesis pathways.
Parthenocarpy is considered the most cost-effective solution for improving the fruit set rate when pollination or fertilization is suppressed by sub-optimum growth conditions, such as low temperature, weak light intensity or facility environments, and ensuring yields of vegetable and fruit crops that are self-sterile or gynoecious. Moreover, seedless fruits produced by parthenocarpy have a better texture, appearance and shelf life  and avoid yield loss caused by seed development [2,3,4].
Parthenocarpy is a widely recognized hormone-dependent biological process. Independent evidence indicates that auxins play a vital role in parthenocarpic fruit set. The genes linking the auxin signal transduction pathway to fruit set have been identified in recent decades. The involvement of IAA9, a member of the tomato Aux/IAA gene family of transcriptional regulators, was confirmed in tomato fruit set. Auxin dose-response assays showed that the down-regulation of IAA9 led to auxin hypersensitivity and resulted in parthenocarpy . Another auxin signaling component involved in fruit set is ARF8, which was identified as a candidate gene for two parthenocarpy QTLs in tomato . Based on the findings of Goetz et al. , a model was proposed for the mechanism of parthenocarpic induction. According to this model, ARF8 forms an inhibitory complex together with an AUX/IAA protein, possibly IAA9, to repress the transcription of the auxin response genes and consequently induce parthenocarpy. Wittwer et al.  showed that a second class of hormone, gibberellins (GAs), could also stimulate parthenocarpic fruit set. The only known gibberellin signaling component shown to be involved in fruit set is DELLA. The reduction in SlDELLA mRNA levels induces the formation of parthenocarpic tomato fruit . Null and loss-of-function recessive mutations in the DELLA genes of Arabidopsis provoke a constitutive GA-response phenotype, including parthenocarpy . Besides Auxin and GAs, Cytokinin (CK) is also involved in parthenocarpy, which accumulates to high levels in ovaries during fruit set [11,12,13,14,15]. Recent studies suggested that CKs may induce parthenocarpy partially through modulation of IAA and GA metabolisms [16,17,18]. Ethylene and abscisic acid (ABA) also play important roles in the regulation of fruit set and development. Ethylene is likely involved in the fruit set program by functioning coordinately with auxin [19,20,21]. ABA may acts as an antagonist of GA or auxin to induce and maintain the dormant state of ovaries, likely by repressing their transition to fruit . These studies demonstrate the complicated and confusing relationships among hormone responses during fruit set. However, the key integrating molecular players remain largely undiscovered, and a global understanding of the mechanisms underlying parthenocarpy has yet to be attained.
Genetic studies have suggested that a majority of natural parthenocarpic properties in crops are quantitative traits regulated by both genetic and environmental factors [6, 22,23,24,25,26]. Photoperiod, temperature, light intensity and nutritional conditions have considerable influences on parthenocarpy [27,28,29]. A hypothesis proposed in the 1930s suggested that plant developmental responses to environmental stimuli were due to the spatiotemporal variations in phytohormone synthesis and transport [30, 31]. Studies have suggested that short-daylight conditions could enhance parthenocarpy by increasing the activity of auxin, while high temperatures suppressed the parthenocarpy rate by inhibiting the synthesis of auxin and gibberellin in the ovary of cucumber [32, 33]. Kim et al.  also found that ovaries had twice the auxin content at 15 °C than that at 25 °C, resulting in a higher rate of parthenocarpy in cucumber. Anyhow, the agricultural application of parthenocarpy was limited by its environmental sensitivity. In practice, the excessive application of exogenous hormones was often used to overwhelm the environmental effects on hormone synthesis, thereby inducing environmentally stable parthenocarpy.
Cucumber (Cucumis sativus L.) is emerging as a model for plants in the Cucurbitaceae family because of its small and fully sequenced genome (2n = 2x = 14, 367 Mb genome) . The mechanism for sex determination, vascular system development, and typical pepo fruit are also well documented in cucumber. The rich parthenocarpic germplasm of cucumber offers an opportunity to investigate the coordination and communication of hormone signals and genes during parthenocarpic fruit set. Genetic studies of parthenocarpy in cucumber started in 1930s. Early studies suggested that parthenocarpy in cucumber is controlled by single genes [36,37,38,39,40]. While most recent studies confirmed that inheritance of parthenocarpy in cucumber is consistent with characteristics of quantitative traits [25, 26, 41,42,43]. Although cucumber is rich in parthenocarpic germplasm resources e.g. main branch specialized/lateral branch specialized parthenocarpy, temperature/photoperiods sensitive parthenocarpy and parthenocarpy with accelerated ovary expansion before anthesis, the agricultural application of parthenocarpic cucumber was limited by their environmental sensitivity . A serial of studies indicated that stable parthenocarpy of cucumber can be induced by auxin or auxin transport inhibitors [34, 44,45,46,47,48]. Meanwhile artificially increasing of endogenous auxin in the ovary by introducing the DefH9-iaaM auxin-synthesizing gene into cucumber might also stimulate parthenocarpy . Besides auxin, application of other hormones such as cytokinins, gibberenllins and brassinosteroids (BRs) could also promote parthenocarpy in cucumber [50, 51]. However, it reported that auxin and GAs had less potential to induce parthenocarpic fruit growth than CKs in cucumber [34, 51, 52]. Therefore, in practice application of exogenous cytokinin, particularly CPPU (N-(2-chloro-4-pyridyl)-N′-phenyl urea, a type of synthetic cytokinin), to induce parthenocarpy is widely used in cucumber production.
In previous studies, an excellent parthenocarpic cucumber cultivar, ‘EC1,’ was found, which showed environmentally stable parthenocarpy under different culture conditions [25, 43]. Previous transcriptome studies have demonstrated that the natural parthenocarpy (NP) of ‘EC1’ has many different aspects compared with cytokinin-induced parthenocarpy (CP) at the mRNA level . However, mRNA levels are not always in accordance with protein activity. To improve the global understanding of parthenocarpy and address the environmental stability of parthenocarpy in cucumber, we conducted a physiological analysis and an iTRAQ (isobaric tags for relative and absolute quantitation)-based proteomic analysis in the natural parthenocarpic fruits of ‘EC1’ and cytokinin-induced parthenocarpic fruits of ‘8419 s-1’ (a non-parthenocarpic variety).
Physiological comparison of the parthenocarpic and non-parthenocarpic cucumber cultivars
Experiments were conducted to investigate the physiological differences between parthenocarpic cultivar ‘EC1’ and non-parthenocarpic cultivar ‘8419 s-1’ during natural/cytokinin-induced parthenocarpy, pollinated fruit set and unpollinated fruit abortion. The detail information of the two cultivars was described in Material and Methods section. The longitudinal and radial growth of the natural parthenocarpic fruits of ‘EC1’ and CPPU -induced parthenocarpy of ‘8419 s-1’, pollinated and unpollinated fruits of ‘8419 s-1’ were measured (Fig. 1a). Our results showed that the length and diameter of the parthenocarpic and pollinated fruits linearly increased from 0 to 6 dpa (days post-anthesis), and the natural and CPPU-induced parthenocarpic fruits showed similar growth curves, wherein the fruit size was generally larger than the pollinated fruits. In contrast, the growth of the unpollinated fruits of ‘8419 s-1’ was blocked, and the length and diameter of the abortive fruits also decreased slightly.
Kim et al.  suggested that genetic factor for parthenocarpy in cucumber may be associated with high content of IAA in the ovaries at anthesis. In this study, endogenous auxins, cytokinins, and gibberellins were analyzed in the cucumber fruits noted above (Fig. 1b). The induction of both naturally occurring and hormone induced parthenocarpy is attributed to the presence of sufficient phytohormones in the ovaries [54,55,56,57,58]. However, parthenocarpic fruits of ‘EC1’ had relatively low and stable hormone levels compared with the fruits of ‘8419 s-1.’ Moreover, the auxin and gibberellin concentrations also decreased unexpectedly during the natural parthenocarpic fruit set of ‘EC1’ (Fig. 1b).
For a further comparison of the fruit developmental differences between ‘EC1’ and ‘8419 s-1,’ we conducted ovary treatment experiments. The ovaries of ‘EC1’ and ‘8419 s-1’ were treated with hormones, hormone inhibitors and pollen separately at anthesis. The weight, length and diameter measurement of the treated ovaries was conducted at 4dpa to reveal different phytohormone responses between ‘EC1’ and ‘8419 s-1’. In cucumber, etiolation of ovary tips is the principal identifying symbol to identify whether the fruits are set or aborted since 2 dpa. It was showed that the ovaries with etiolated tips did not grow or even wilt by comparing with the 0 dpa ovaries (Table 1). The non-etiolated tip phenotypes and growth of ovaries suggested that parthenocarpy of ‘8419 s-1’ could be induced by all of the exogenous hormones, including NAA, CPPU, GA3 and EBR; however, NAA, GA3 and EBR exerted weak effects on fruit growth, causing the treated ovaries to grow slightly in length and remain in a dormant-like state (Table 1, Additional file 1: Figure S1). Marcelis et al.  suggested that cell division in fruit of cucumber occurs about a week after anthesis whereas cell size increases markedly only after cell division begins to decline. Fruit size increases because of increase in both the number and size of cells. Therefore we speculated that NAA, GA3 and EBR might be primarily involved in cell division, rather than cell expansion during fruit development of cucumber, thus ovaries treated by these hormones are smaller than pollinated and CCPU-treated fruits. On the other hand, the weak effect of NAA, GA3 and EBR might be due to the low hormone concentration used in this study. Fu et al.  showed that 0.2 μM (20 mg/L) of EBR led to high efficiency of fruit growth in cucumber. But in this study we found 10 mg/L of EBR showed weak effect on fruit expansion. However, half concentration of CPPU (50 mg/L) referred to Fu’s method still lead to strong effect on fruit set and growth.
The strong effect of CPPU on fruit development of cucumber was observed and documented in many other reports [34, 51, 52]. Similar observation was also found in many other species such as watermelon, apple, kiwifruit and blueberry [60,61,62,63]. Early fruit development generally consists of three stages: fruit formation, cell division and cell expansion . The growth of cucumber fruit size is often mirrored by the increase in cell number and size . Previous cytological observation showed that the CPPU-treated cucumber fruits was initiated with an increase of cell numbers in the pericarp and placenta tissues, and the size of pericarp cells were bigger than natural parthenocarpic fruits, although the number of cell layers was similar . The findings suggested that CPPU might be involved in both processes of cell division and cell expansion. In addition, fruit growth is tightly related to the availability of carbohydrate, because fruit is a very strong metabolic sink. Many studies confirmed that CKs was demonstrated to regulate carbohydrate allocation in fruit [65, 66]. We thought that maybe another reason why CPPU induced parthenocarpic fruit was consistently bigger than the fruit induced by other PRGs.
Martínez et al.  have demonstrated that the inhibition of ethylene response (STS treatment) is sufficient to induce the set and early development of the fruit in absence of pollination in both the parthenocarpic and the non-parthenocarpic cultivar of zucchini squash. Coincidentally, it was showed that STS has stimulated parthenocarpy in the three non-parthenocarpic cucumber cultivars; however it had no effect on fruit development of parthenocarpic cultivar ‘EC1’ (Additional file 2: Table S6). Besides, diameter, length and weight of ethephon treated fruit of ‘EC1’ showed no significant differences to the natural parthenocarpic fruits (Additional file 2: Table S6). It indicated that neither ethephon nor STS (ethylene response inhibitor) could affect parthenocarpy of EC. Interestingly, the ethephon treated ovaries of non-parthenocarpic cucumbers displayed more severe atrophy of ovaries by comparing with their unpollinated ovaries (Additional file 2: Table S6), suggested that the fruit abortion of the non-parthenocarpic cultivars maybe accelerated by ethephon. Irradiated pollen treatment could promote stenospermocarpy in 8419 s-1, but the stenospermocarpic ovaries were much smaller than the active pollen-treated ovaries, implying that seed set may be essential for the fruit growth of non-parthenocarpic varieties. Interestingly, the seedless fruit of ‘EC1’ formed by parthenocarpy was much larger than its pollinated fruit. The pollination fruits of ‘8419 s-1’ were blocked by either a hormone inhibitor mixture or individual hormone inhibitors (Table 1). Etiolation was observed in the pollen and hormone inhibitor co-treated ovaries. In contrast, the natural parthenocarpic fruit set of ‘EC1’ could not be blocked by hormone inhibitors, but the growth of the ovaries was suppressed (Table 1).
iTRAQ-based proteomic study of differently developed cucumber fruits
Early fruit development generally consists of three stages: fruit formation, cell division and cell expansion . The earliest stage, in which the ovary is aborted or allowed to proceed with fruit development, is referred to as fruit set. In this study, we focused on the fruit set stage in cucumber, which is 0 to 2 days post-anthesis [53, 59, 64, 67, 68]. The natural parthenocarpic fruits of ‘EC1,’ cytokinin-induced parthenocarpic fruits of ‘8419 s-1’, the pollinated and unpollinated fruits of ‘8419 s-1’ were investigated (described in the M&M section). The proteomes of these differently developed fruits were analyzed by iTRAQ with two technical replicates per sample. The strategy for analysis is shown in Additional file 3: Figure S2.
Protein homolog identification was conducted by BLASTP against the cucumber Refseq database and the Arabidopsis thaliana Refseq database (E-value <1E-10). After redundancies were removed, 683 unique proteins were identified, including 359 fruit set-related proteins (natural/cytokinin-induced parthenocarpic fruit and pollinated fruit) and 377 fruit abortion-related proteins. Gene Ontology (GO) analysis showed that the proteins detected in both set and aborted fruits were involved in similar biological processes (Fig. 2a). Proteins related to cell growth (GO:0016049), cell cycle (GO:0007049), cell death (GO:0008219) and cell communication (GO:0007154) were detected that actively expressed during fruit development, that was consistent with our previous transcriptomic study , suggesting the cell growth, cell cycle, cell death and cell communication related genes were actively involved in fruit development at both mRNA and protein levels.
Identification and comparison of differentially expressed proteins involved in different fruit developmental processes
Large numbers of polymorphic SNPs (Single Nucleotide Polymorphisms), InDels (Insertion-Deletion Polymorphisms) and SVs (Structural Variations) were detected between ‘EC1’ and ‘8419 s-1’ by genome resequencing study . Besides, 84 differentially expressed proteins were also identified between 0dpa fruits of ‘EC1’ and ‘8419 s-1’ (Additional file 4: Figure S3). These findings demonstrated that the genetic background of the two varieties was significantly different with each other. In order to identifiy fruit developmental DEPs (Differentially Expressed Proteins), a proteome comparison strategy was employed (Additional file 3: Figure S2). Proteomic comparisons were not conducted between cultivars, therefore the differentially expression of identified DEPs were not caused by genetic variations but only correlated with development of fruit. The false discovery rate (FDR) method (FDR ≤ 0.05, |fold ≥ 1.5|, P-value < 0.05) was used to determine the significance of the differential protein expression. Consequently 138 DEPs were identified and four groups of DEPs were screened (Additional file 4: Figure S3). In pollination and CP fruits, most of the DEPs were up-regulated; in contrast, more down-regulated proteins were detected in the abortive fruit of ‘8419 s-1’ and the natural parthenocarpic (NP) fruit. Clustering analysis indicated that the DEPs in CP and pollinated fruit showed similar protein expression patterns, but the protein expression profiles of the CP and NP fruits were clustered into two groups, and the cluster distance between CP and NP fruits was greater than that between CP and the abortive fruit (Fig. 2b).
Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis of the DEPs was conducted using MapMan software , according to the biological pathway maps of Arabidopsis (http://mapman.gabipd.org/web/guest/home). Large numbers of DEPs were shown to be involved in the biological process of protein synthesis and were mainly up-regulated during fruit set and down-regulated during fruit abortion (Fig. 2c). Transcriptome analysis has shown that amino acid metabolism, glycolysis and TCA cycle-related genes were actively expressed during fruit set , corroborating the present proteomic analysis, which also showed that the proteins related to these biological processes were actively expressed (Fig. 2c).
The interactions between the DEPs were analyzed based on the reference proteome-wide binary protein-protein interaction (PPI) map of Arabidopsis . A core fruit developmental PPI network was revealed, which consisted of 30 DEPs and 19 ‘bridging’ interaction proteins (Fig. 3, Additional file 2: Table S1). The proteins within the PPI network were mainly involved in the processes of protein metabolism (GO:0019538), transport (GO:0006810) and signal transduction (GO:0007165). Interactions within the PPI network occurred more frequently in fruit abortion, while the smallest scale of interactions was involved in natural parthenocarpy (Additional file 5: Figure S4). Two interaction proteins exhibited specialized expression in natural parthenocarpy: CER9 (ECERIFERUM 9, Csa7M073540.1), which is involved in cuticle metabolism and the maintenance of plant water status , and PRL (PROLIFERA, Csa7M407650.1), which is specifically expressed in populations of dividing cells in the sporophytic tissues of the plant body . TOC159 (Csa1M229500.1, a chloroplast biogenesis-related protein), IMPA-6 (Csa1M597740.1, a nuclear import protein) and RS6 (Csa1M229500.1, a putative ovule development regulator) showed specialized up-regulation in CPPU-induced parthenocarpic fruit set (Additional file 2: Table S1).
The common and specific DEPs in the differently developed cucumber fruits are shown in the Venn diagram in Fig. 4. No DEPs were commonly expressed in NP, CP and pollination fruit. Twelve common proteins were identified in the CP and pollinated fruits and showed similar expression patterns (Additional file 2: Table S2). Most of these CP and pollination-specialized DEPs are closely related to pollen and seed development. Eleven DEPs were commonly expressed in NP and abortive fruit, most of which were protein metabolism-related proteins and mainly involved in the biological processes of pollen germination, gametophyte and endosperm development as well as root morphogenesis; however, these proteins showed opposite expression trends during parthenocarpy and fruit abortion (Additional file 2: Table S2).
Forty-one parthenocarpy-specialized DEPs were identified, but only 3 DEPs are commonly found in the NP and CP fruits (Additional file 2: Table S3). The remaining 39 parthenocarpy-specialized DEPs were divided into two groups, of which 19 were uniquely expressed in NP fruits, while the other 19 existed only in CP fruits (Table 2). These DEPs were mainly involved in the biological processes of pollen germination, seed and seedling development, cell proliferation and programmed cell death (Table 2). The NP-specialized DEPs Csa7M073540.1 and Csa7M450640.1, which were related to cell cycle and proliferation, and Csa1M025890.1 and Csa4M036590.1, the amino acid biosynthesis-related proteins, showed dramatically up-regulated expression during fruit set (expression fold >5) (Table 2). Moreover, Csa2M139820.1, which has the putative function of translational elongation, was the only DEP that was dramatically increased in the CP fruit (Table 2).
Expression analysis of NP- and CP-specialized proteins in response to phytohormones
The transcription profiles of the parthenocarpy-specialized DEPs in the hormone- and hormone inhibitor-treated fruits were compared (the treated fruits are described in the 1st part of the results section). Consistent with the result of iTRAQ, the specialized DEPs from both parthenocarpy groups showed inactive transcription during pollinated fruit set (Fig. 5, middle panel). The CP-specialized DEPs were actively expressed in the cytokinin-induced parthenocarpic fruits of ‘8419 s-1.’ However, in the parthenocarpic fruits of ‘EC1,’ including the hormone inhibitor-treated but not blocked parthenocarpic fruits, most proteins showed decreased expression (Fig. 5, left column of the panel). In contrast, the NP-specialized DEPs showed up-regulated expression in parthenocarpic fruits of ‘EC1’, whereas many of these DEPs were silenced in the cytokinin-induced parthenocarpic fruits (Fig. 5, right column of the panel). Although a few of the NP-specialized DEPs were differentially transcribed in cytokinin-induced parthenocarpic fruits, the expression levels of the genes in these fruits were the same as those in the parthenocarpic fruits of ‘EC1.’ Moreover, similar transcription patterns of NP-specialized DEPs between the natural parthenocarpic and unblocked natural parthenocarpic fruits were found, indicating that the transcription of these proteins could not be affected by hormone inhibitors (Fig. 5, top right of the panel). These findings, to some extent, indicated that the transcription of NP-specialized DEPs was not sensitive to hormones or hormone inhibitors.
iTRAQ showed that four NP-specialized proteins—Csa7M073540.1, Csa7M450640.1, Csa1M025890.1 and Csa4M036590.1—showed a high abundance in expression during parthenocarpy (expression >5-fold; Table 2, marked by stars). To further investigate the expression characteristics of the active parthenocarpy-specialized proteins, we conducted a western blot analysis. Consistent with the results of iTRAQ, the NP-specialized proteins were up-regulated during natural parthenocarpic fruit set but inactively expressed in the cytokinin-induced fruit, while the NP-specialized protein was not detected in the NP fruit but was actively expressed in the CP fruit (Fig. 6a). Auxins, cytokinins and gibberellins are essential fruit developmental phytohormones, which can induce parthenocarpy in cucumber. The expression patterns of the parthenocarpy-specialized proteins responding to these hormones were also investigated (Fig. 6b). Our results showed that the treatment with endogenous hormones NAA, CPPU and GA3 increased the expression of the CP-specialized protein Csa2M139820.1 but decreased the expression of Csa4M036590.1 and Csa7M450640.1 (NP-specialized proteins). Although the NP-specialized proteins Csa1M025890.1 and Csa7M073540.1 were up-regulated by GA3, they were not significantly changed by NAA and CPPU (Fig. 6b, Additional file 6: Figure S5A).
The protein expression pattern of Csa2M059750.1 was also analyzed, which was up-regulated during natural parthenocarpy and down-regulated in abortive fruit (Additional file 2: Table S3). Csa2M059750.1 was the unique DEP located in the chromosome region of the major parthenocarpic QTL Parth2.1 of ‘EC1’ . Western blot analysis showed that Csa2M059750.1 was degraded during fruit abortion but actively expressed during fruit set. However, in contrast to the natural parthenocarpic fruit set, the increasing expression of Csa2M059750.1 was delayed in the pollinated and CP fruits until 3 dpa (Fig. 7a). Western blot analysis also showed that the expression of Csa2M059750.1 was not significantly affected by hormones, indicating that Csa2M059750.1 may be a hormone-insensitive protein (Fig. 7b, Additional file 6: Figure S5B).
Although studies on parthenocarpy have been conducted for over 100 years, current understanding of parthenocarpy remains at a nascent stage. Cucumber is emerging as a model for the Cucurbitaceae family. The rich parthenocarpic germplasm of cucumber offers an opportunity to investigate the coordination and communication of environmental factors, hormone signals and genes during parthenocarpy. In this study, we investigated the proteomes of cucumber fruits to help improve the global understanding of parthenocarpy.
Post-translational regulation of fruit development in cucumber
Proteins are executors with a vast array of functions within organisms. The translational and post-translational regulations of proteins such as protein synthesis, proteolysis, glycosylation, phosphorylation and folding are essential for plant development. Studies on hormone-dependent biological processes have revealed that post-translational regulations frequently occur during fruit development. For instance, in the absence of auxins, the function of ARFs (auxin response factors) was inhibited via heterologous dimerization with Aux/IAA. However, in the presence of auxin, the ARFs dissociated from Aux/IAA proteins, whose targeting and degradation were mediated by the E3 ubiquitin ligase SCFTIR1/AFB, and consequently stimulated fruit set [5, 73,74,75,76,77]. The fruit developmental responses to ethylene were mediated by the SCFEBF1/EBF2-dependent proteolysis of EIN3 (Ethylene insensitive 3) . Furthermore, the function of EIN3 was regulated by the MAPK-dependent phosphorylation within the EPR1 domain of the transcriptional factor . In addition, our previous transcriptome study confirmed that glycosylation reactions were dramatically active throughout fruit development in cucumber . In the present proteomic study, the protein folding-related proteins were up-regulated in both pollinated and parthenocarpic fruits, including Csa1M255160.1, which was annotated as a TCP-1 (T-COMPLEX PROTEIN 1 ALPHA SUBUNIT), and was actively expressed in natural and cytokinin-induced parthenocarpic fruits (Additional file 2: Table S3). Moreover, Csa2M099450.1, also defined as a TCP-1-like protein, showed specialized expression in pollinated and cytokinin-induced parthenocarpic fruits (Additional file 2: Table S1).
Many lines of evidence indicate that protein synthesis/degradation may function during cell growth [80, 81]. The present proteomic study showed that over 30% of the differentially expressed proteins during the cucumber fruit development were related to the biological process of protein metabolism (Fig. 2c). Within the predicted IPP network, nearly half of the interaction proteins were protein metabolism-related proteins, such as the three TIF (TRANSLATION INITIATION FACTOR) proteins involved in the initiation phase of eukaryotic translation, four CSN (COP9 signalosome) multi-proteins that functioned in the ubiquitin–proteasome pathway, and three ribosomal proteins (Additional file 2: Table S1). Protein metabolism related proteins were commonly expressed in unpollinated and natural parthenocarpic fruit (Additional file 2: Table S2). The opposite expression patterns of the common proteins indicated the fate (set or abortion) of the mature ovary in cucumber, which was determined by protein metabolism pathways.
The cues of hormone-independent parthenocarpy in cucumber
Growth measurement showed that natural and cytokinin-induced parthenocarpic fruits presented similar growth curves (Fig. 1a). However, clustering analysis showed that the protein expression profiles of the CP and NP fruits were quite different from each other. Moreover, the cluster distance between the CP and NP fruits was greater than that between the CP fruit and the abortive fruit (Fig. 2b). The specialized proteins expressed in the parthenocarpic fruits were divided into two individual groups, which were separately involved in the NP and CP fruit set, as shown in the Venn diagram (Fig. 4; Table 2). These findings suggested that there may be individual parthenocarpic pathways in cucumber.
Gustafson [54, 55] proposed that plants produce parthenocarpic fruits because the ovary contains enough auxins to promote fruit initiation. Since then, many studies have confirmed that parthenocarpy is a phytohormone-dependent biological process. In cucumber, polar auxin transport-blocking experiments have shown that parthenocarpy could be triggered by the sufficient accumulation of auxin in the ovary . Moreover, the application of exogenous hormones such as auxins, cytokinins, gibberellins and brassinosteroids could induce parthenocarpy . In this study, hormone measurement showed that the endogenous hormone levels increased during fruit set but decreased during fruit abortion in ‘8419 s-1’ (Fig. 1b). However, the endogenous hormone levels were relatively low and remained stable during natural parthenocarpic fruit set in ‘EC1’ compared with ‘8419 s-1.’ Moreover, the NP fruits showed a broad resistance to hormone inhibitors (Fig. 1b; Table 1; Additional file 1: Figure S1), indicating the existence of a hormone-independent parthenocarpic mechanism in ‘EC1.’ This speculation was supported by expression analysis in the parthenocarpy-specialized proteins, whereby the NP-specialized proteins performed hormone-insensitive transcriptional and translational functions (Figs. 5, 6 and 7; Additional file 6: Figure S5).
Inhibiting the regulation of fruit abortion in cucumber
Dormant fruits, as a result of first-fruit inhibition or nutritional stress, can always be observed in the field [25, 26]. However, the dormant state of these fruits usually leads to fruit abortion in a short time (2 days at most). Although natural parthenocarpic fruits of ‘EC1’ could not be blocked by hormone inhibitors, these treated fruits stayed in a dormant state for a long time (more than 4 days) (Additional file 1: Figure S1). We speculated that inhibitory regulations of fruit abortion might exist in ‘EC1’, causing the fruits to maintain a dormant state. Coincidentally, the common proteins that detected in the NP fruits of ‘EC1’ and the abortive fruit of ‘8419 s-1,’ showed opposite expression trends. Most of these proteins were down-regulated during fruit abortion but up-regulated during NP fruit set (Additional file 2: Table S2). Besides, ethephon treating experiments showed that although fruit abortion of the non-parthenocarpic cultivars was accelerated by ethephon which had no effect on fruit development of ‘EC1’, further suggesting inhibitory regulations of fruit abortion might exist in ‘EC1’. Conversely, hormone inhibitor-induced dormant state in ‘EC1’ indicated that hormone stimuli might be required for fruit expansion in either parthenocarpic cultivars or non-parthenocarpic cultivars.
Based on the evidence provided in this study, a working hypothesis for the cucumber parthenocarpic fruit set was proposed (Fig. 8), whereby parthenocarpy in cucumber may be promoted by a ‘parallel switch,’ namely, hormone-dependent and hormone-independent pathways. During hormone-independent parthenocarpy, fruit set was promoted by hormone-insensitive regulatory proteins, such as the NP-specialized proteins in ‘EC1.’ In the presence of sufficient hormones, young fruits formed through both hormone-dependent and -independent pathways could continuously grow to maturity. In the absence of hormones, the development of hormone-sensitive fruits proceeds to fruit abortion, whereas the hormone-insensitive fruits remain in a dormant state because of the increasing expression of abortion-inhibiting proteins. However, the expansion of dormant fruits and their further promotion are unknown. Although the accurate regulation of parthenocarpy in cucumber remains unclear, our studies provide a theoretical framework for understanding the mechanism of parthenocarpy for its application in agricultural production.
Plant material and growth conditions
In this study, the cucumber cultivar ‘EC1’ was used as a parthenocarpic sample (Gynoecious inbred line, European Glasshouse type, parthenocarpic rate ≥ 95%) and ‘8419 s-1’ as a non-parthenocarpic sample (Monoecious inbred line, European Glasshouse type, the rare occurrence of parthenocarpy is occasionally observed in the senescence phase of the cultivar). Plants were grown in a greenhouse at Nanjing Agricultural University with a 14 h photoperiod, a mean daily air temperature of 28/20 °C (day/night).
Phytohormones were separately analyzed through ELISA using IAA, ZR and GA3 ELISA Kits (Sangon Biotech Company) based on Weiler’s method . The results are presented as the mean ± SE (n = 10) with three technical replicates.
The female flowers (at the 12-15th node of the main stem) of the above cucumber cultivars were previously trapped with bags in order to prevent pollen contamination on the day before anthesis. When anthesis, the trapped ovaries were treated separately: keeping trapping (unpollination), hand pollination  and CPPU treatment. CPPU (N-(2-chloro-4-pyridyl)-N′-phenyl urea) is a kind of synthetic cytokinin which could induce parthenocarpy in cucumber. For CPPU treatment, 20 μL CPPU solution (100 mg/L) was sprayed on the surface of the ovaries. All the treated ovaries were harvested at 0, 1, 2 and 3 dpa (days-post-anthesis). Thirty ovaries of each treatment were ground into powder with liquid nitrogen and mix as a sample pool for iTRAQ and western blot analysis.
The trapped ovaries of ‘EC1’ and ‘8419 s-1’ were also treated with phytohormones (NAA, 1-naphthaleneacetic acid, 50 mg/L; CPPU, 100 mg/L; GA3, Gibberellin A3, 50 mg/L; EBR, Epi-Brassinosteroids, 10 mg/L; Ethephon, 100 mg/L) and hormone inhibitors (TIBA, 2,3,5-triodobenzoic acid, inhibitor of auxin, 50 mg/L; Lovastatin, inhibitor of cytokinin, 50 mg/L; uniconazole, inhibitor of gibberellin 50 mg/L; Brz, Brassinazole, inhibitor of Brassinosteroids, 10 mg/L; STS, silver thiosulphate, inhibitor of ethylene, 0.25 mM). The exogenous phytohormones and hormone inhibitors were separately sprayed on the surface of the ovaries at 0dpa. Active pollens and irradiated pollens (γ-ray irradiation at a dose of 200Gy) also used to treat the 0dpa ovaries by hand-pollination. After spaying and hand-pollination, the ovaries were trapped again. The weight, length and diameter of the ovaries were measured at 4dpa. The experiments were repeated three times (n = 30). The treated ovaries were also harvested at 2dpa, of which the RNA was isolated for qRT-PCR analysis.
For western blotting, the seedlings of ‘8419 s-1’ (at the three true leaf stage) were also treated with exogenous phytohormones (50 μM, 10 μM, or 5 μM NAA; 10 μM CPPU and 10 μM GA3) by spraying the solutions on the surface of the true leaves. After growing in the growth chamber with a 14 h photoperiod and 25 °C for 24 h, in total 30 true leaves from five individual plants by same treatment were collected and mixed by grinding in liquid nitrogen, then stored at −80 °C before protein extraction.
Protein extraction and quantization
Approximately 1 g of powdered sample was mix with 3 mL extraction buffer [500 mM Tris-HCl (pH 7.5), 150 mM NaCl, 50 mM ethylene diaminetetraacetic acid (EDTA), 1% Triton-X-100, 2 mM dithiothreitol (DTT), 2 mM phenylmethanesulfonyl fluoride (PMSF)]. Protein extraction was performed using the methods described by Omar et al. . The protein precipitation was collected and washed with cold methanol containing 10 mM DTT three times, cold acetone containing 10 mM DTT twice and then dried by vacuum freeze. The extracted proteins were quantified by using the Bradford method .
Protein digestion and iTRAQ labeling
One hundred micrograms Proteins from each samples were precipitated with five volume of cold acetone at −20 °C for 1, centrifuged by 12,000 rpm for 15 min at 4 °C, and dried by vacuum freeze dryer (Thermo savant, USA). Pellets were dissolved in the dissolution buffer with reducing reagent described in iTRAQ Reagent 8-Plex kit (Applied Biosystems, USA), and alkylated by cysteine-blocking reagent according to the manufacturer’s instructions . After digestion with 50 μl of 50 ng/μl sequence grade modified trypsin (Promega, USA) solution overnight at 37 °C, the peptide samples were labeled. The samples were labeled with the iTRAQ tags as described in Additional file 1: Figure S1.
SCX chromatography and LC–MS/MS analysis
The vacuum dried iTRAQ labeled samples were re-suspended with 100 μl SCX (Strong cation exchange) buffer A (10 mM ammonium formate, 20% ACN (acetonitrile), pH 2.8) and fractionated using a Poly-SEA HPLC (High Performance Liquid Chromatography) column (2.0 × 150 mm, 5 μm particle size, 300 Ǻ pore size) using at a flow rate of 0.3 ml/min on the Agilent 1200 HPLC System (Agilent, USA). The 50 min HPLC gradient consisted of 100% buffer A (10 mM ammonium formate, 20%ACN, pH 2.8) for 5 min, 0–50% buffer B (500 mM ammonium formate, 20%ACN, pH 2.8) for 25 min, then 50–80% buffer B for 15 min, followed by 80–100% buffer B for 10 min, and lastly 100% buffer B for 15 min. Chromatograms were recorded at 215 and 280 nm. All the collected fractions were vacuum dried, and re-suspended with Nano- RPLC (Reversed Phase Liquid Chromatography) buffer A (0.1% FA, folic acid; 2%ACN). Samples were desalted with C18 nanoLC trap column (100 μm ID × 3 cm, 3 μm particle size, 150 Ǻ pore size) and Nano-RPLC buffer A (0.1%FA, 2%ACN) at 2 μl/min for 10 min for LC–MS/MS analysis.
The mass spectroscopy analysis was performed using a Triple TOF 5600 System (AB SCIEX, USA), coupled with the Eksigent nanoLC-Ultra™ 2D System (AB SCIEX, USA).The iTRAQ labeled peptides were separated using an analytical ChromXP C18 column (75 μm ID × 15 cm, 3 μm particle size, 120 Å pore size) (New Objectives, USA) with a nanospray emitter (2500 V, 30 PSI (pounds per square inch) curtain gas, 5 PSI nebulizer gas, 150 °C interface heater temperature) (New Objectives, USA), and analyzed by LC-MS/MS. A rolling collision energy setting was applied to all precursor ions for collision-induced dissociation (CID). For information dependent acquisition (IDA), survey scans were acquired in 250 ms and as many as 35 product ion scans were collected if they exceeded a threshold of 150 counts per second (counts/s) with a 2+ to 5+ charge-state. The total cycle time was fixed to 2.5 s. Dynamic exclusion was set for one-half of peak width (18 s), and then the precursor was refreshed off the exclusion list. The peak areas of the iTRAQ reporter ions reflect the abundance of the proteins in the samples.
Mass spectrometric data was processed with Protein Pilot Software v. 4.0 (AB SCIEX, USA) against Cucumber database using the Paragon algorithm, and further processed by a Pro Group algorithm where isoform-specific quantification was adopted to trace the differences between expressions of various isoforms which was applied to the peptide identification. Protein identification was performed with emphasis on biological modifications option. Database search parameters were the followings: instrument was TripleTOF 5600, iTRAQ 8-plex quantification, cysteine modified with iodoacetamide, biological modifications were selected as the ID focus, trypsin digestion. An automatic decoy database search strategy was employed to estimate the false discovery rate (FDR) using the Proteomics System Performance Evaluation Pipeline Software (PSPEP) t was integrated in the Protein Pilot Software. In this study, only protein quantification data with the value of global FDR ≤0.05 were chosen for further analysis, and proteins with a |fold change ≥1.5| were considered to be significantly differentially expressed.
Quantitative real time PCR
Proteins based on their differential expression patterns revealed by iTRAQ were selected for verification by Quantitative real-time PCR (qRT-PCR) with primers designed using Primer 5.0 software (Additional file 2: Table S4). Total RNA of the samples described above was extracted by Trizol (Invitrogen, USA). After extraction, total RNA was treated with DNase I (Fermentas, UK) according to the manufacturer’s protocol. First-strand cDNA synthesis was carried-out using the PrimeScript™ RT-PCR Kit (TaKaRa, Japan). The real-time qRT-PCR was accomplished in a thermal cycler and analyzed by an IQ5 multicolor Real-time PCR detection system (Bio-Rad, USA). To determine relative fold differences for each sample in each experiment, the CT values were normalized using Cs-actin as an internal control and calculated relative to a calibrator using the formula 2-△△Ct. The experiment was repeated three times.
Western blot analysis
Proteins extracted from ovaries and leaves as previously described were mixed with protein lysis buffer at 4 °C. Total protein lysis was boiled at 98 °C for 10 min and separated with 10% SDS-PAGE (Sodium dodecyl sulfate-polyacrylamide gel electrophoresis), then transferred to nitrocellulose (NC) membranes (GE Hybond, USA) with semi-dry approach. After 2 h of blocking with 5% milk in TBST (Tris Buffered Saline with Tween), membranes were incubated with polyclonal antibodies against the proteins of cucumber that were raised in rabbit by synthetic peptides (Lufei, P.R.China; Additional file 2: Table S5). Polyclonal antibody against cucumber beta-actin was used as internal control. The antibodies were used at 1:300 dilutions. Membranes were incubated with goat anti-rabbit IgG (Proteintech, USA) at 1:2000 dilutions in TBST for 2 h, after that membranes were washed with TBST for five times. Signals were detected by using enhanced electro-chemiluminescence (Beyotime, P.R. China).
Brassinazole, a kind of inhibitor of Brassinosteroids
Cytokinin induced parthenocarpy
N-(2-chloro-4-pyridyl)-N0-phenylurea, a diphenylurea-derived cytokinin
Differentially expressed protein
days post anthesis
False discovery rate
isobaric tags for relative and absolute quantitation
Quantitative trait locus
2,3,5-triodobenzoic acid, a kind of inhibitor of auxin
Pandolfini T. Seedless fruit production by hormonal regulation of fruit set. Nutrients. 2009;1(2):168–77.
Tiedjens VA. The relation of environment to shape of fruit in Cucumis Sativus L. and its bearing on the genetic potentialities of the plants. J Agr Res. 1928;36:794–809.
Denna DW. Effects of genetic parthenocarpy and gynecious flowering habit on fruit production and growth of cucumber Cucumis Sativus L. J Am Soc Hortic Sci. 1973;98(6):602–4.
Falavigna A, Soressi GP. Influence of the pat-sha gene on plant and fruit traits in tomato (L. esculentum mill.) In: Modern trends in tomato genetics and breeding; 1987. p. 128.
Wang H, Jones B, Li Z, et al. The tomato aux/IAA transcription factor IAA9 is involved in fruit development and leaf morphogenesis. Plant Cell. 2005;17(10):2676–92.
Gorguet B, Eggink PM, Ocana J, et al. Mapping and characterization of novel parthenocarpy QTLs in tomato. Theor Appl Genet. 2008;116(6):755–67.
Goetz M, Vivian-Smith A, Johnson SD, et al. AUXIN RESPONSE FACTOR8 is a negative regulator of fruit initiation in Arabidopsis. Plant Cell. 2006;18(8):1873–86.
Wittwer SH, Bukovac MJ, Sell HM, et al. Some effects of gibberellin on flowering and fruit setting. Plant Physiol. 1957;32(1):39.
Martí C, Orzáez D, Ellul P, et al. Silencing of DELLA induces facultative parthenocarpy in tomato fruits. Plant J. 2007;52(5):865–76.
Carrera E, Ruiz-Rivero O, Peres LEP, et al. Characterization of the procera tomato mutant shows novel functions of the SlDELLA protein in the control of flower morphology, cell division and expansion, and the auxin-signaling pathway during fruit-set and development. Plant Physiol. 2012;160(3):1581–96.
Mapelli S. Changes in cytokinin in the fruits of parthenocarpic and normal tomatoes. Plant Sci Lett. 1981;22(3):227–33.
Bohner J, Bangerth F. Effects of fruit set sequence and defoliation on cell number, cell size and hormone levels of tomato fruits (Lycopersicon esculentum mill.) within a truss. Plant Growth Regul. 1988;7(3):141–55.
Gillaspy G, Ben-David H, Gruissem W. Fruits: a developmental perspective. Plant Cell. 1993;5(10):1439.
Srivastava A, Handa AK. Hormonal regulation of tomato fruit development: a molecular perspective. J Plant Growth Regul. 2005;24(2):67–82.
Trueman SJ. Endogenous cytokinin levels during early fruit development of macadamia. Afr J Agric Res. 2010;5(24):3402–7.
Matsuo S, Kikuchi K, Fukuda M, et al. Roles and regulation of cytokinins in tomato fruit development. J Exp Bot. 2012;63(15):5569–79.
Ding J, Chen B, Xia X, et al. Cytokinin-induced parthenocarpic fruit development in tomato is partly dependent on enhanced gibberellin and auxin biosynthesis. PLoS One. 2013;8(7):e70080.
Li Y, Yu JQ. Photosynthesis and 14C-assimilate distribution as influenced by CPPU treatment on ovary. Acta Agric Nucleatae Sin. 2001;15:355–9. (in Chinese with English abstract)
Vriezen WH, Feron R, Maretto F, et al. Changes in tomato ovary transcriptome demonstrate complex hormonal regulation of fruit set. New Phytol. 2008;177(1):60–76.
Pascual L, Blanca JM, Cañizares J, et al. Transcriptomic analysis of tomato carpel development reveals alterations in ethylene and gibberellin synthesis during pat3/pat4 parthenocarpic fruit set. BMC Plant Biol. 2009;9(1):67.
Martínez C, Manzano S, Megías Z, et al. Involvement of ethylene biosynthesis and signalling in fruit set and early fruit development in zucchini squash (Cucurbita pepo L.). BMC Plant Biol. 2013;13(1):139.
Fos M, Nuez F, García-martínez JL. The gene pat-2, which induces natural parthenocarpy, alters the gibberellin content in unpollinated tomato ovaries. Plant Physiol. 2000;122(2):471–80.
Beraldi D, Picarella ME, Soressi GP, et al. Fine mapping of the parthenocarpic fruit (pat) mutation in tomato. Theor Appl Genet. 2004;108(2):209–16.
Miyatake K, Saito T, Negoro S, et al. Development of selective markers linked to a major QTL for parthenocarpy in eggplant (Solanum melongena L.). Theor Appl Genet. 2012;124(8):1403–13.
Wu Z, Zhang T, Li L, et al. Identification of a stable major-effect QTL (Parth 2.1) controlling parthenocarpy in cucumber and associated candidate gene analysis via whole genome re-sequencing. BMC Plant Biol. 2016;16(1):182.
Lietzow CD, Zhu H, Pandey S, et al. QTL mapping of parthenocarpic fruit set in north American processing cucumber. Theor Appl Genet. 2016;129(12):2387–401.
Smith O, Cochran HL. Effect of temperature on pollen germination and tube growth in the tomato. New York: The University; 1935.
Rylski I. Effects of season on parthenocarpic and fertilized summer squash (Cucumis pepo L.). Exp Agric. 1974;10(01):39–44.
Sun C, Li Y, Zhao W, et al. Integration of hormonal and nutritional cues orchestrates progressive corolla opening. Plant Physiol. 2016;171(2):1209–29.
Cholodny N. Beiträge zur hormonalen Theorie von Tropismen. Planta. 1928;6(1):118–34.
Went FW, Thimann KV. Phytohormones. NY: Macmillan; 1937.
Rudich J, Baker LR, Sell HM. Parthenocarpy in Cucumis sativus L. as affected by genetic parthenocarpy, thermo-photoperiod, and femaleness. J Am Soc Hort Sci 1977;102(2):225-8.
Matlob AN, Kelly WC. Growth regulator activity and parthenocarpic fruit production in snake melon and cucumber grown at high temperature. J Am Soc Hortic Sci. 1975;100:406–9.
Kim IS, Okubo H, Fujieda K. Endogenous levels of IAA in relation to parthenocarpy in cucumber (Cucumis sativus L.). Sci Hortic. 1992;52(1–2):1–8.
Huang S, Li R, Zhang Z, et al. The genome of the cucumber, Cucumis sativus L. Nat Genet. 2009;41(12):1275–81.
Hawthorn LR, Wellington R. Geneva, a greenhouse cucumber that develops fruit without pollination. N Y State Agric Exp Station. 1930;2:3–11.
Pike LM, Peterson CE. Inheritance of parthenocarpy in the cucumber (Cucumis sativus L.). Euphytica. 1969;18:101–5.
Kvasnikov BV, Rogova NT, Tarakanova SI, Ignatov SI. Methods of breeding vegetable crops under the covered ground. Trudy Prikl Bot Genet Selek. 1970;42:45–57.
Juldasheva L. Inheritance of the tendency towards parthenocarpy in cucumbers. Byull Vsesoyuznogo Ordena Lenina Inst Rastenievodstva Imeni NI Vavilova. 1973;32:58–9.
Meshcherov E, Juldasheva L. Parthenocarpy in cucumber. Trudy Prikl Bot Genet Selek. 1974;51:204–13.
Shawaf EI, Baker L. Inheritance of parthenocarpic yield in gynoecious pickling cucumber for once-over mechanical harvest by diallel analysis of six gynoecious lines. J Am Soc Hortic Sci. 1981;106:359–64.
Sun ZY, Lower RL, Staub JE. Analysis of generation means and components of variance for parthenocarpy in cucumber (Cucumis sativus L.). Plant Breed. 2006;125:277–80.
Yan LY, Lou LN, Li XL, et al. Evaluation of parthenocarpy in cucumber germplasm. Acta Hortic Sin. 2009;36:975–82. (in Chinese with English abstract)
Robinson RW, Cantliffe DJ, Shannon S. Morphactin-induced parthenocarpy in the cucumber. Science. 1971;171(3977):1251–2.
Cantliffe D. Parthenocarpy in the cucumber induced by some plant growth-regulating chemicals. Can J Plant Sci. 1972;52(5):781–5.
Quebedeaux B, Beyer E. Chemically-induced parthenocarpy cucumber by a new inhibitor of auxin transport. Hortscience. 1972;7:474–6.
Elassar GJ, Patevitch D, Kedar N. Induction of Parthenocarpic fruit development cucumber by growth regulators. Hortscience. 1974;9(3):238–9.
Takeno K, Ise H, Minowa H, et al. Fruit growth induced by benzyladenine in Cucumis Sativus L.: influence of benzyladenine on cell division, cell enlargement and indole-3-acetic acid content. J Jpn Soc Hortic Sci. 1992;60(4):915–20.
Yin Z, Malinowski R, Ziolkowska A, et al. The DefH9-iaaM-containing construct efficiently induces parthenocarpy in cucumber. Cell Mol Biol Lett. 2006;11(2):279–90.
Ogawa Y, Inoue N, Aoki S. Promotive effects of exogenous and endogenous gibberellins on the fruit development in Cucumis sativus L. J Jpn Soc Hortic Sci. 1989;58(2):327–31.
Fu FQ, Mao WH, Shi K, et al. A role of brassinosteroids in early fruit development in cucumber. J Exp Bot. 2008;59(9):2299–308.
Hikosaka S, Sugiyama N. Effects of exogenous plant growth regulators on yield, fruit growth, and concentration of endogenous hormones in Gynoecious Parthenocarpic cucumber (Cucumis sativus L.). Hortic J. 2015;84(4):342–9.
Li J, Wu Z, Cui L, et al. Transcriptome comparison of global distinctive features between pollination and parthenocarpic fruit set reveals transcriptional phytohormone cross-talk in cucumber (Cucumis sativus L.). Plant Cell Physiol. 2014;55(7):1325–42.
Gustafson FG. Further studies on artificial parthenocarpy. Am J Bot. 1938a;25:237–44.
Gustafson FG. Induced parthenocarpy. Bot Gaz. 1938b;99(4):840–4.
Schwabe WW. Hormones and parthenocarpic fruit set, a literature survey. Hortic Abstr. 1981;51:661–98.
Vivian-Smith A, Koltunow AM. Genetic analysis of growth-regulator-induced parthenocarpy in Arabidopsis. Plant Physiol. 1999;121(2):437–52.
Spena A, Rotino GL, Bhajwani SS, Soh WY. Parthenocarpy. State of the art, Current trends in the embryology of Angiosperms. Kluwer Academic Publishers; 2001. p. 435-450.
Marcelis LFM, Baan Hofman-Eijer LR. Cell division and expansion in the cucumber fruit. J Hortic Sci. 1993;68(5):665–71.
Huitrón MV, Diaz M, Diánez F, et al. Effect of 2, 4-D and CPPU on triploid watermelon production and quality. Hortscience. 2007;42(3):559–64.
Bangerth F, Schröder M. Strong synergistic effects of gibberellins with the synthetic cytokinin N-(2-chloro-4-pyridyl)-N-phenylurea on parthenocarpic fruit set and some other fruit characteristics of apple. Plant Growth Regul. 1994;15(3):293–302.
Lewis DH, Burge GK, Hopping ME, et al. Cytokinins and fruit development in the kiwifruit (Actinidia deliciosa). II. Effects of reduced pollination and CPPU application. Physiol Plant. 1996;98(1):187–95.
NeSmith DS. Response of rabbiteye blueberry (Vaccinium ashei Reade) to the growth regulators CPPU and gibberellic acid. Hortscience. 2002;37(4):666–8.
Boonkorkaew P, Hikosaka S, Sugiyama N. Effect of pollination on cell division, cell enlargement, and endogenous hormones in fruit development in a gynoecious cucumber. Sci Hortic. 2008;116(1):1–7.
Fang JB, Tian LL, Li SH, et al. Influence of CPPU on the sink and source of kiwifruit. Acta Horticulturae Sin. 2000;27(6):444–6. (in Chinese with English abstract)
Zeng H, Yang W, Lu C, et al. Effect of CPPU on carbohydrate and endogenous hormone levels in young macadamia fruit. PLoS One. 2016;11(7):e0158705.
Ando K, Grumet R. Transcriptional profiling of rapidly growing cucumber fruit by 454-pyrosequencing analysis. J Am Soc Hortic Sci. 2010;135(4):291–302.
Ando K, Carr KM, Grumet R. Transcriptome analyses of early cucumber fruit growth identifies distinct gene modules associated with phases of development. BMC Genomics. 2012;13(1):518.
Thimm O, Bläsing O, Gibon Y, et al. MAPMAN: a user-driven tool to display genomics data sets onto diagrams of metabolic pathways and other biological processes. Plant J. 2004;37(6):914–39.
Dreze M, Carvunis A-R, Charloteaux B, et al. Evidence for network evolution in an Arabidopsis interactome map. Science. 2011;333:601–7.
Lü S, Zhao H, Des Marais DL, et al. Arabidopsis ECERIFERUM9 involvement in cuticle formation and maintenance of plant water status. Plant Physiol. 2012;159(3):930–44.
Springer PS, Holding DR, Groover A, et al. The essential Mcm7 protein PROLIFERA is localized to the nucleus of dividing cells during the G (1) phase and is required maternally for early Arabidopsis development. Development. 2000;127(9):1815–22.
De Jong M, Wolters-Arts M, Feron R, et al. The Solanum Lycopersicum auxin response factor 7 (SlARF7) regulates auxin signaling during tomato fruit set and development. Plant J. 2009;57(1):160–70.
Goetz M, Hooper LC, Johnson SD, et al. Expression of aberrant forms of AUXIN RESPONSE FACTOR8 stimulates parthenocarpy in Arabidopsis and tomato. Plant Physiol. 2007;145(2):351–66.
Ren Z, Li Z, Miao Q, et al. The auxin receptor homologue in Solanum Lycopersicum stimulates tomato fruit set and leaf morphogenesis. J Exp Bot. 2011;62(8):2815–26.
Dharmasiri N, Dharmasiri S, Weijers D, et al. Plant development is regulated by a family of auxin receptor F box proteins. Dev Cell. 2005;9(1):109–19.
Kepinski S, Leyser O. The Arabidopsis F-box protein TIR1 is an auxin receptor. Nature. 2005;435(7041):446–51.
Guo H, Ecker JR. Plant responses to ethylene gas are mediated by SCF EBF1/EBF2-dependent proteolysis of EIN3 transcription factor. Cell. 2003;115(6):667–77.
Li J, Li Z, Tang L, et al. A conserved phosphorylation site regulates the transcriptional function of ETHYLENE-INSENSITIVE3-like1 in tomato. J Exp Bot. 2012;63(1):427–39.
Key JL. Ribonucleic acid and protein synthesis as essential processes for cell elongation. Plant Physiol. 1964;39(3):365.
Faurobert M, Mihr C, Bertin N, et al. Major proteome variations associated with cherry tomato pericarp development and ripening. Plant Physiol. 2007;143(3):1327–46.
Weiler EW, Jourdan PS, Conrad W. Levels of indole-3-acetic acid in intact and decapitated coleoptiles as determined by a specific and highly sensitive solid-phase enzyme immunoassay. Planta. 1981;153(6):561–71.
Omar AA, Song WY, Grosser JW. Introduction of Xa21, a Xanthomonas-resistance gene from rice, into ‘Hamlin’ sweet orange [Citrus sinensis (L.) Osbeck] using protoplast-GFP co-transformation or single plasmid transformation. J Hortic Sci Biotechnol. 2007;82(6):914–23.
Bradford MM. A rapid and sensitive method for the quantitation of microgram quantities of protein utilizing the principle of protein-dye binding. Anal Biochem. 1976;72(1–2):248–54.
Wisniewski JR, Zougman A, Nagaraj N, et al. Universal sample preparation method for proteome analysis. Nat Methods. 2009;6(5):359.
Wang Y, Zhang WZ, Song LF, et al. Transcriptome analyses show changes in gene expression to accompany pollen germination and tube growth in Arabidopsis. Plant Physiol. 2008;148(3):1201–11.
Rajjou L, Belghazi M, Huguet R, et al. Proteomic investigation of the effect of salicylic acid on Arabidopsis seed germination and establishment of early defense mechanisms. Plant Physiol. 2006;141(3):910–23.
Menges M, Hennig L, Gruissem W, et al. Cell cycle-regulated gene expression in Arabidopsis. J Biol Chem. 2002;277(44):41987–2002.
Hajduch M, Hearne LB, Miernyk JA, et al. Systems analysis of seed filling in Arabidopsis: using general linear modeling to assess concordance of transcript and protein expression. Plant Physiol. 2010;152(4):2078–87.
Hwang IS, Kim NH, Choi DS, et al. Overexpression of Xanthomonas campestris pv. Vesicatoria effector AvrBsT in Arabidopsis triggers plant cell death, disease and defense responses. Planta. 2012;236(4):1191–204.
Irshad M, Canut H, Borderies G, et al. A new picture of cell wall protein dynamics in elongating cells of Arabidopsis Thaliana: confirmed actors and newcomers. BMC Plant Biol. 2008;8(1):94.
We would like to thank Professor Robert N. Trigiano (University of Tennessee, USA) for critical reading of the manuscript.
This work was supported by the National Natural Science Foundation of China [31430075, 31672168]; Special Fund for Agro-Scientific Research in the Public Interest ; National Key Research and Development Program of China [2016YFD0101705-5]; Agricultural science and Technology Innovation Fund of Jiangsu Province [CX(15)1019]. 31430075, 31672168 were involved in iTRAQ analysis. 201403032, 2016YFD0101705-5 and CX(15)1019 were involved in sample collection. The funding bodies were not involved in the design of the study and collection, analysis, and interpretation of data and in writing the manuscript.
Availability of data and materials
The datasets generated during the current study are available from the corresponding author on reasonable request.
Ethics approval and consent to participate
The plant material used in this work has been generated for research. All the materials were grown in Jiangpu experimental station of Nanjing Agricultural University. The seeds are maintained in State Key Laboratory of Crop Genetics and Germplasm Enhancement (Nanjing Agricultural University).
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
The typical phenotypes of the treated ovaries of ‘EC1’ and ‘8419 s-1’ at 4 dpa. (DOCX 196 kb)
The annotation and expressional information of the proteins within the protein-protein interactions (PPIs) which presented in Fig. 5. Table S2. DEPs commonly expressed in pollinated and Cytokinin induced parthenocarpic fruit (I), and DEPs commonly existed in unpollinated and natural parthenocarpic fruit (II). Table S3. DEPs commonly expressed in natural and cytokinin induced parthenocarpic fruits. Table S4. The primers of natural and cytokinin induced parthenocarpy specialized protein encoding genes for qRT-PCR. Table S5. The antigenic peptides of parthenocarpy specialized proteins for producing rabbit polyclonal antibodies. Table S6. The measurement of ovary weight, length and diameter after treated by ethephon and ethylene inhibitors.The treatments were conducted at the anthesis day (0 dpa), and the measurements were conducted at 4 dpa. Means (±SE) of three independent experiments were calculated.Letters indicate differences between the treated ovaries with statistical significance at P ≤ 0.05 (t-test).Detail information of EC1 and 8419 s-1 was described in M&M section. Both cucumber cultivars CC3 and CCMC were non-parthenocarpic, monoecious inbred lines, Asia ecotype. (XLSX 35 kb)
The analytical strategy of the iTRAQ based proteome analysis of cucumber fruits. Total proteins of each sample were extracted and labeled separately with tags (113 to 118). Proteomic analysis was conducted by iTRAQ. The differentially expressed proteins (DEPs) were identified by comparing proteomes of 0 and 2 dpa ovaries of each treatment. CK1 and CK2: Control sample 1 and Control sample 2; NP: natural parthenocarpic fruits of EC1; Unp: Unpollination fruits of 8419 s-1 (fruit abortion); P: pollination fruits of 8419 s-1; CP: Cytokinin induced parthenocarpic fruits of 8419 s-1. Bar = 20 mm. (DOCX 357 kb)
Statistics of differentially expressed proteins during different fruit developmental processes of cucumber. NP: natural parthenocarpic fruits of EC1; Unp: Unpollination fruits of 8419 s-1 (fruit abortion); P: pollination fruits of 8419 s-1; CP: Cytokinin induced parthenocarpic fruits of 8419 s-1. (DOCX 247 kb)
Protein interactions occurred in different fruit developmental processes based on the predicted PPIs network. The interactions marked in red color means the interactions involved in Cytokinin induced parthenocarpy (A), natural parthenocarpy (B), pollination fruit set (C) and unpollination fruit abortion (D). (DOCX 349 kb)
Band intensity analysis of western blotting by using the software ImageJ. The relative expression fold of each parthenocarpy specialized protein was calculated by the formula: (band intensity after hormone treatments/band intensity without hormone treatment)/(band intensity of beta-actin in hormone treated sample/band intensity of beta-actin in untreated sample). Each value represents the mean ± SE of three Western blotting replicates. (DOCX 238 kb)