Skip to main content

Transcriptomic study of pedicels from GA3-treated table grape genotypes with different susceptibility to berry drop reveals responses elicited in cell wall yield, primary growth and phenylpropanoids synthesis

Abstract

Background

Gibberellins (GA3) are the most sprayed growth regulator for table grape production worldwide, increasing berry size of seedless varieties through pericarp cell expansion. However, these treatments also exacerbate berry drop, which has a detrimental effect on the postharvest quality of commercialized clusters. Several studies have suggested that pedicel stiffening caused by GA3would have a role in this disorder. Nevertheless, transcriptional and phenotypic information regarding pedicel responses to GA3is minimal.

Results

Characterization of responses to GA3treatments using the lines L23 and Thompson Seedless showed that the former was up to six times more susceptible to berry drop than the latter. GA3also increased the diameter and dry matter percentage of the pedicel on both genotypes. Induction of lignin biosynthesis-related genes by GA3has been reported, so the quantity of this polymer was measured. The acetyl bromide method detected a decreased concentration of lignin 7 days after GA3treatment, due to a higher cell wall yield of the isolated fractions of GA3-treated pedicel samples which caused a dilution effect. Thus, an initial enrichment of primary cell wall components in response to GA3was suggested, particularly in the L23 background. A transcriptomic profiling was performed to identify which genes were associated with these phenotypic changes. This analysis identified 1281 and 1787 genes differentially upregulated by GA3in L23 and cv. Thompson Seedless, respectively. Concomitantly, 1202 and 1317 downregulated genes were detected in L23 and cv. Thompson Seedless (FDR < 0.05). Gene ontology analysis of upregulated genes showed enrichment in pathways including phenylpropanoids, cell wall metabolism, xylem development, photosynthesis and the cell cycle at 7 days post GA3application. Twelve genes were characterized by qPCR and striking differences were observed between genotypes, mainly in genes related to cell wall synthesis.

Conclusions

High levels of berry drop are related to an early strong response of primary cell wall synthesis in the pedicel promoted by GA3treatment. Genetic backgrounds can produce similar phenotypic responses to GA3, although there is considerable variation in the regulation of genes in terms of which are expressed, and the extent of transcript levels achieved within the same time frame.

Background

The exogenous application of gibberellic acid-3 (GA3) during key phases of fruit development is critical for seedless grape varieties, since it promotes berry enlargement through elongation of pericarp cells [1,2]. Although these treatments are beneficial in terms of meeting the market standards at expected fruit caliber, they also exacerbate berry drop (also known as shattering) during postharvest [3]. This disorder affects negatively the general quality of the fruit, generating relevant economic losses. Its extent can also increase over time depending on postharvest storage conditions [4].

Preventing decay through the application of systematic measures during the postharvest handling of table grape is critical to achieve the commercialization of healthy clusters [5]. Hence it is important that marketable clusters present a non-senescent aspect with green stems and pedicels attached to the berry. The pedicel has been proposed as a main factor involved in berry drop [6] and several varieties have been described to be more prone to this postharvest issue under GA3治疗的条件。Berry drop in cv. Kyoho is related to an increase in rachis hardness caused by expansion and lignification of this structure as a consequence of GA3treatments [7] . Other commercial cultivars have been reported as susceptible, such as Flame Seedless [5] and Ruby Seedless [3]. Thus, genetic background and the growth regulators applied to certain cultivars, particularly GA3, seem to play an important role in postharvest berry drop. However, there are also varieties such as cv: Thompson Seedless which can be considered as non-sensitive to berry drop as long as proper management conditions are achieved: with only three GA doses applied for berry size enlargement spanning at specific fruit developmental stages [4,6].

Reports have indicated several morphological changes that take place in the rachis and pedicel following GA3treatment [8,9]. It has been found that GA3induces changes in the pedicel by increasing the area of xylem and pith structures [8]. One study proposed that the loss of pedicel flexibility could be one of the main factors underlying berry drop in cv. Thompson Seedless [10]. Recently, analysis of transcripts related to lignin synthesis on genetic backgrounds with contrasting susceptibility to berry drop was reported, suggesting an enhancement of lignification in response to GA3treatments [11]. In other woody species, an increase in the concentration of bioactive gibberellins (GA) enhanced growth and biomass production [12,13]. In fact, transgenic lines have been engineered to increase bioactive GA and repress lignin biosynthesis in xylem tissue, resulting in enriched biomass destined to liquid fuel production [14,15]. Therefore, it is suggested that the expected effects of GA3treatments in vascular tissue could be on cell wall-related processes, which in turn impact on the flexibility of this structure.

Transcriptomic approaches have been very useful to characterize the underlying biology of berry development [16,17,18] and its responses to gibberellin treatments [19,20]. These studies have been widely conducted on samples obtained from fruit tissues. In this case, our focus is the study of the responses of pedicel tissue to GA3通过转录组platf治疗评估orm, and how this could contribute to understand a complex trait such as berry drop.

The availability of contrasting phenotypes for postharvest berry drop may serve as a useful reference to understand this phenomenon [11]. L23, a line obtained in the framework of INIA’s table grape breeding program, is a genotype with high susceptibility to berry drop compared to cv. Thompson Seedless, considered as a non-sensitive genotype in this study, under similar agroclimatic and management conditions. By capturing transcriptional variation in response to GA3treatments from two contrasting genotypes for berry drop, valuable information can be extracted to unveil the underlying differential response through the identification of differentially expressed genes. The objective of this study was to characterize the early phenotypic and transcriptional responses to GA3application of two genotypes that show contrasting susceptibility to berry drop and to identify the main biological processes involved in each condition.

Results

Genotype L23 is more susceptible to postharvest berry drop than cv. Thompson seedless under GA3treatment conditions

比较的微分响应L23和简历。Thompson Seedless in terms of postharvest berry drop, these genotypes were characterized across three seasons following GA3-treatment (Fig.1). The expected effect on berry enlargement after GA3application was confirmed (Fig.1a). Significant differences in berry drop among conditions by season followed the same trend in all cases; L23 showed higher levels of berry drop than cv. Thompson Seedless but only under GA3treatment conditions (Fig.1b). In the first season L23 was up to 5.97-fold more susceptible to berry drop. Values for the second and third seasons were 4.11 and 3.11-fold, respectively.

Fig. 1
figure1

GA3not only increases berry and pedicel dimensions but also exacerbates postharvest berry drop.aEffects of GA3application on bunches from cv. Thompson Seedless. The image illustrates the notorious increase on berry size produced by GA3compared to a non-treated cluster of the same vine.bL23 exhibits high susceptibility to postharvest berry drop under GA3治疗的条件。总重量的浆果环保总局rated from cap stems against total weight of bunches after shaking is reported. Mean and standard deviation are shown on each bar and error bar (replicates for Season 1, Season 2 and Season 3 were 15, 8, 8, respectively, on each point). Red abscissa line in the graph shows 3% berry drop, a recommended tolerance value for high quality table grape exportation. Data were analyzed with ANOVA and Tukey HSD tests to find significant differences between conditions (p < 0.05).cPedicel dimensions increase in GA3-treated samples. Scatter plots of variables show differences in dimension of berry and pedicel. Total soluble solids and firmness were used as indicators of phenological development. V + n: Number of weeks (n) passed since veraison (V). The numbers ‘0’, ‘15’, ‘30’ represent the number of days spent in cold storage after harvest time (n = 30)

GA3treatment for berry enlargement promotes pedicel growth and enhances primary cell wall yield

Since GA3treatment is associated to postharvest berry drop, phenotypic variables of treated plants were observed during several developmental stages. Results of the first season are shown in Fig.1c; two attributes considered as maturity indicators in the study were firmness and total soluble solids of berry samples. As can be seen, the expected effect on berry enlargement was observed under GA3treatment conditions for berry diameter and length. A similar response of the pedicel to GA3treatment was detected, particularly in the thickness of this tissue.

Changes in pedicel dimensions induced by GA3have already been reported [11], and were observed in this study (Additional file1: Figure S1A). In our study, enlargement of cortex cells and transcript levels of monolignol biosynthesis-related genes (4CL,CAD6andCCR1L) were upregulated by GA3treatment. Transcriptional induction of lignin biosynthesis genes along with increased cell size could be related to the greater stiffness of pedicels as proposed in [10]. To test if the amount of lignin in samples treated with GA3is actually higher in these groups, lignin concentration was determined by the acetyl bromide (AB) method [21]. Successful validation of this method adapted to pedicel tissue was performed using the same reference tissue as in the original research article (Additional file1: Figure S1B and C).

As Fig.2shows, lignin concentration in treated groups was unexpectedly lower than in control samples; this difference was more evident at early stages of development, particularly following GA3application (2A). The difference between control and treated conditions in lignin was less accentuated in the véraison stage (2B), and at harvest the observed variation was too great to detect significant differences among conditions (2C). Measurement of samples obtained at a postharvest stage (15 days of cold storage), showed the same trend that at harvest stage (data not shown). Since the method quantifies the amount of soluble lignin obtained from a protein-free cell wall fraction, measurements are relative to the dry matter obtained in each sample.

Fig. 2
figure2

Lignin concentration in pedicels diminishes after GA3treatment. Quantification of soluble lignin by acetyl bromide method of pedicels sampled from three developmental stages. Mean and standard deviation are shown on each bar and error bar, respectively (n = 8). ANOVA followed by Tukey tests were performed to detect significant differences between conditions (p < 0.05). Stages of evaluation were:aSeven days after treatment,bveraison andcHarvest. Homoscedasticity of data was previously verified by the Levene test, withp-values of 0.6583 (a), 0.1799 (b) and 0.1252 (c)

Considering that pedicel dimensions differed between control and treated groups, changes in dry matter among conditions were recorded. Table1reports changes induced by GA3related to dry matter and cell wall accumulation. Dry matter increased in both genotypes; for cell wall yield, which is proportional to the protein-free cell wall fraction obtained from the AB protocol, GA3-treated samples showed a significant difference for the L23 genotype.

Table 1 Dry matter and cell wall deposition are increased in pedicel samples from GA3treated clusters

In summary, independent of the berry drop (shattering) amount associated to each genotype, the results suggest that size enlargement-related changes in the pedicel induced by GA3treatments could be related to promotion of primary cell wall yield and secondary cell wall modification-related processes.

Early transcriptional responses of pedicel to GA3are characterized by an enrichment of processes such as cell cycle, photosynthesis, cell wall modification, phenylpropanoid metabolism and xylem development

To characterize the early transcriptional responses of pedicel tissue to GA3treatment a transcriptomic profiling of both genotypes to this growth regulator was implemented. Information about treatments and timing of samples is summarized in Fig.3a (the pipeline followed to perform differential expression analysis is illustrated in Additional file2: Figure S2).

Fig. 3
figure3

Characterization of the early response to GA3-treatment in two backgrounds with different susceptibility to berry drop.aThe scheme illustrates the rationale behind the timing of pedicel samples since GA3-treatment is a major factor involved in berry drop incidence. Phenological-based sampling was critical since the L23 and cv. Thompson Seedless fruit development time frames differ slightly (2–3 weeks). The reference time frame was based on the study of [16].bUpregulated genes in response to GA3 treatments. The number of DE genes was 1787 and 1281 in cv. Thompson seedless and L23, respectively (FDR < 0.05).cRepressed genes in response to GA3 treatments. The number of DE genes was 1317 and 1202 in cv. Thompson seedless and L23 s, respectively (FDR < 0.05)

The number of differentially expressed genes was obtained as indicated in the methods section, after testing relative to a fold-change threshold [22] (Additional file3: Figure S3). Differential expression analysis revealed a total of 1281 and 1787 upregulated genes for L23 and cv. Thompson Seedless, respectively (Fig.3b). A total of 1202 and 1317 downregulated genes were detected for L23 and cv. Thompson Seedless (Fig.3b). This is the first report about differential gene expression in GA3treatments in genotypes with contrasting susceptibility to berry drop (See report for differential expression analysis in Additional file4: Data S1).

To test if the differentially expressed genes (DEGs) were significantly linked to specific biological processes, gene set enrichment analysis was conducted in gene ontology (GO) annotation terms. Significant GO terms revealed numerous functions associated with enhanced primary cell growth, cell wall modification, histone and chromatin modification, among many others. GO functions of repressed genes were enriched in stress and defense-related processes (See full report in Additional file5: Data S2).

To reduce redundancy and detect meaningful information among the GO terms obtained, the method described by Supek et al. [23] was followed. Enriched GO terms for upregulated genes in L23 is shown in Fig.4。The same analysis conducted on cv. Thompson Seedless can be seen in Additional file6: Figure S4 (a list with detailed terms and p-adjusted values is given in Additional file7: Table S1 for L23 and Additional file8: Table S2 for cv. Thompson Seedless). Several GO terms grouped to cell growth and DNA replication processes. Other relevant functions were related to photosynthesis, xylem development, phenylpropanoid and flavonoid metabolism and cell wall modification processes, among many others.

Fig. 4
figure4

Comparison between genotypes of GO terms enrichment analysis for the upregulated genes by GA3。A list of the 73 most prominent and least redundant GO accessions obtained from revigo analysis of L23 (Cutoff value: 0.7) is shown here; this set of GO accessions is compared to the results obtained in cv. Thompson seedless. The log10adjusted p-value is shown here; the size is relative to the number of query items matching the corresponding GO accession obtained from the single enrichment analysis. The fold-change of every DE gene matching each GO accession by genotype was averaged and is shown according to the scale color gradient depicted in the figure. Input of DE upregulated genes was 1281 and 1787 for L23 and cv. Thompson seedless, respectively

To corroborate data from transcriptome profiling and validate the results, transcript abundance detected by RNAseq analysis for 18 genes was compared to the measurements obtained by quantitative PCR (qPCR). Genes sampled from RNAseq analysis are given in Additional file9: Table S3 along with the primers designed for qPCR. Results for the correlation tests between these techniques are given in Additional file9: Figure S5; qPCR correlated significantly with RNAseq measurements (p = 1.368e-13,R2 = 0.80), validating the transcriptomic assay.

In summary, transcriptomic responses to GA3are mainly associated to the regulation of primary cell wall metabolism and cell cycle regulation; which could support primary growth and, by extension, be associated to cell wall fraction yield. Other features such as histone modification, positive regulation of photosynthesis-related genes, flavonoids and phenylpropanoids metabolism, along with xylem development were strongly regulated by GA3in pedicel tissue.

Genes related to phenylpropanoid, photosynthesis and cell wall pathways were identified in the highly susceptible L23 background as possible biomarkers of berry drop

To gain insight about the highly regulated genes in the susceptible L23 background, as an approach to associate berry drop and gene regulation, a list of 30 highly regulated genes was obtained from the transcriptomic data. For this, the DEGs grouping under functional categories were identified through the method proposed by [24]; Mapman results can be seen in Additional file10: Figure S6. Relevant categories were cell wall, along with lignin and some secondary metabolism-related genes. The main reason for the selection of these categories was their high representation as functional categories in terms of associated probability values; all significant categories identified by this method can be seen in Fig.5。Quantitative transcriptome values of gene expression were fitted to a linear model, considering lowest to highest levels of observed berry drop values (cv. Thompson Seedless untreated, L23 untreated, cv. Thompson Seedless treated, L23 treated) [25]. Regression analysis was performed, and 12 genes were selected (Fig.6).

Fig. 5
figure5

Enriched functional categories by Mapman analysis in L23 genes regulated by GA3。检测到的所有重要类别照办xon test with p-value adjusted by Benjamini-Yeuketeli are shown on the Cleveland dot plot. Red segmented vertical line depicts threshold significance value (p adjusted value< 0.05). Functional categories were consistent with the results given by the gene ontology analysis and were used as reference in the identification of possible biomarker genes according to the genetic background assessed for the berry drop trait

Fig. 6
figure6

Clustering and heatmap of transcriptomic data from potential candidate genes identified from RNAseq analysisPotential markers were extracted from DE genes. Data corresponds to fragments per kilobase of exon per million reads (FPKM). Below are depicted the libraries from RNAseq analysis (TC: Thompson Untreated, LC: L23 Untreated, TG: Thompson Treated, LG: L23 Treated. Numbers correspond to biological replicates)

Expression studies of L23 showed predominant regulation of flavonoids, xylem-related genes and secondary cell wall formation

The list of the 12 genes proposed to be studied at several developmental stages comparing the two genotypes is given in Table2, along with the primers used to measure the abundance of the transcripts by quantitative PCR. Additional file14shows that the expression of candidate genes in pedicel samples collected from several stages revealed genes related to secondary cell wall formation such as cellulose synthase subunit-7 (CESA7) and endoglucanase (GUN10)that are overexpressed under GA3treatment conditions, particularly in the L23 background. Phenylpropanoid-related genes such as phenylalanine ammonia lyase (PALY), caffeoyl CoA O-methyltransferase (CCoAOMT) and cinnamoyl CoA-reductase (CCR) were highly expressed in cv. Thompson Seedless and L23. These are also positively regulated by GA3(7DAT) as well as by cold treatment (30 days at 0 °C), most notoriously for thePALYgene.

Table 2 Primer list of candidate genes to characterize the underlying responses to GA3in pedicel samples from genotypes with contrasting susceptibility to berry drop

Flavonoid 3′-hydroxylase (F3H), an enzyme belonging to the flavonoid biosynthesis pathway, is highly expressed in the L23 genotype and upregulated by GA3。The expression pattern of this gene in the susceptible genotype is strikingly different from cv. Thompson. Some laccases, also differentially expressed, such as Laccase-4 (LAC4) and Laccase-17 (LAC17) may exert an important role in lignin condensation, although their expression was relatively low. Other genes such as Nephentesin-1 (NEP1), Xylem serine proteinase (XSP1) and Chlorophyll a − b binding protein 40 chloroplastic (CB12) had higher basal levels in the L23 background than in cv. Thompson Seedless and were regulated positively by GA3。These sets of genes responsive to GA3showed distinctive expression patterns between genotypes and treatments across several stages of development.

Discussion

Gibberellins are the most used growth regulators in table grape production; particularly in seedless varieties, in which the endogenous generation of the hormone is reduced to basal levels and so the application of GA3is mandatory to obtain berries of marketable size [26]. However, there is also a detrimental effect of GA3, which increases postharvest berry drop, a problem well known for a long time for some cultivars such as Thompson Seedless [6]. Some authors have proposed that the rachis and pedicel could be one of the main structures influencing berry drop [9,10,11]. But scarce information exists about the transcriptional response of this tissue to GA3treatments.

This is the first study to characterize the transcriptomic response of table grape pedicel to GA3in two genetic backgrounds with contrasting susceptibility to berry drop. The main findings were: (i) GA3treatments induced major phenotypical changes in the pedicel, by increasing its size and modifying cell-wall properties; (ii) accumulation of lignin was negatively affected by an early response of cell wall augmentation, resulting in an initial dilution effect with lower concentration of the polymer. As a consequence, lignin concentration reached similar values for both treated and untreated groups at advanced stages such as harvest time; (iii) transcriptome profiling of samples after treatment provided cues on the main biological processes upregulated by GA3,可以概括为促进vegetative growth through positive regulation of cell growth, cell wall modification, xylem development, phenylpropanoid and flavonoid metabolism, along with downregulation of stress-related responses; and (iv) core responses are similar between the genotypes, one more tolerant and one more susceptible to postharvest berry drop, although there are a number of DEGs that are specific to each one. Genes which showed contrasting expression patterns were identified and studied at several stages; this set of genes, or some of them, could be used potentially as indicators or biomarkers of postharvest berry drop. However, to obtain robust indicators it is necessary to generate information in the context of segregating populations to understand the underlying genetics of the shattering trait and subtraits. With this information, comprehensive models can be inferred which could serve to predict the occurrence of this trait in wider genetic backgrounds.

Berry drop in seedless varieties has been linked to growth regulator treatments, most of them correlating with the application of the gibberellic acid-3 (GA3) isoform for berry size enlargement [3,5,6,27,28]. Researchers have reported the effects of this growth regulator on several tissues and conditions. In pedicel tissue, positive regulation by GA3on lignin monomers biosynthesis, cell expansion and xylem development has been detected [8,11,28]. The present study supports these findings in pedicel with the identification of upregulated genes with significant annotations for these processes. However, although gene expression data could suggest a possible upregulation of lignin deposition in this tissue, our results showed scarce alteration of lignin quantity after GA3treatment; as a matter of fact, its concentration was diluted in a cell wall matrix which was augmented in response to this plant growth regulator. Polymerization of lignin is a process which seems to be finely regulated by the cellular surroundings [29,30,31] and could be highly dependent of the analyzed tissue. For instance, a recent study showed that lignin deposition in berry skin is in fact diminished by GA3, along with the activities of enzymes which synthetize several lignin monomers [32]. This would suggest that pedicel stiffness is probably more related to a strong cell wall synthesis response elicited by GA3than to a differential response on lignin deposition, which was proposed elsewhere for the contrasting genotypes studied here [11]. However, some studies have shown that pedicel stiffness is not always the expected phenotypic response to GA3。For instance, in cv. Italia the berry detachment force (an indicator of pedicel strength) is unaffected by GA3[33]. In summary, berry drop, and the phenotypes associated to GA3, could have a wide arrange of responses depending on the genetic background assessed and the tissue analyzed, even under similar agroclimatic conditions. Besides pedicel stiffening, other factors would have taken into account to predict berry drop under postharvest conditions such as brush strength, varietal response to GA3and other relevant postharvest conditions.

Measurement of lignin was coherent with other techniques performed on similar tissues. For example, lignin content in grape stalks has been estimated to be 17.4% by the Klason method (corrected previously by tannin, polyphenols and other substances) [34]. In our study, the range for reported lignin content was 15.7–20.4%. We observed that this diminished content at early stages of GA3application could be explained by a strong response in primary cell wall synthesis that sustains the growth of cell components through creep processes [35]; however, to prove such statement other analyses must be conducted to measure the differences among the different components of the lignocellulosic matrix. This could be an interesting and challenging task to perform in this woody tissue. Nevertheless, recent methodology that facilitates the analysis of cell wall components in other woody species through non-destructive protocols could be critical to correlate cell wall changes with transcriptional data [36] and remains an interesting topic to be addressed in future studies.

A mild but significant difference in cell wall yield was observed between genotypes; however, as mentioned before, the magnitude alone of this parameter will not explain differences in berry drop between backgrounds since the phenomenon could be far more complex and other factors could be involved. The transcriptional profiles of a set of genes, which were mostly involved in primary growth, could correlate with the observed phenotypic changes elicited by GA3。The gene with the most different expression profile between genotypes was Cellulose synthase subunit-7 (CESA7). which has been described as a gene involved in secondary cell wall synthesis [37]. In this study, this gene was highly expressed in the pedicel of L23 background along with several genes encoding other CESA-related subunits (GSVIVG01023643001, GSVIVG01021248001, GSVIVG01028402001, GSVIVG01033278001; see Additional file3更详细的结果)。进一步analyses will be required to evaluate the specific effects of these genes on cell wall synthesis and its interactions with GA3。这也将是有趣的来确定re are sequence differences in these genes and their regulatory upstream regions between cv. Thompson Seedless and L23 backgrounds, to uncover possible regulation motifs or structural variations associated with the different transcription levels detected by transcriptomics and qPCR methods, along with the measurement of cell wall components to corroborate the effects suggested by transcriptional data .

Responses of seedless berry to GA3treatments are strongly oriented towards hormone crosstalk, as observed in cv. Centennial Seedless, which ultimately orchestrates a general response of cell wall relaxation that promotes berry size enlargement [20]. The pedicel, however, showed more enrichment in ‘effector-related’ responses than in hormone-related pathways in the same time frame of 7 days after GA3treatment (7DAT). Remnants of seed tissues in the berry may play a major role in the response to exogenous hormone applications, and since the pedicel is a vascular and supporting tissue of berry growth, the responses of the latter to GA3applications could be focused on downstream pathways, whereas the former is integrating the external stimuli (in this case, mostly GA3) to the signals given by the embryo remnants. Results of the gene set enrichment analysis of the pedicel to GA3(Additional file5: Data S2) showed positive responses on plant cell wall loosening (GO:0009828) and organization (GO:0071669), among other cell wall-related and cell growth processes, which can be attributed to tissue enlargement similar to what has been described in berry tissue. In contrast, hormone-related terms showed scarce enrichment; although, response to gibberellin (GO:0009739) was detected in positive DEGs. Other responses related to hormones in negative DEGs were related to jasmonic acid (GO:0071395), salicylic acid (GO:0009754, GO:0009863, GO:0071446), ABA (GO:009738) and ethylene (GO:0009723, GO:0009692, GO:0009693). It is worth mentioning that in the present study the response to GA3treatment was only evaluated at 7 days after application. As can be seen in [20], sequencing of samples from earlier stages of evaluation (one and 3 days after treatment) showed that hormone-related responses were richer. Thus, it is important to consider how the response to hormone signaling develops over time. Nevertheless, the aim of this study was to evaluate how different is the response to GA3between two genetic backgrounds with contrasting performance in postharvest berry drop.

Enrichment in processes such as photosynthesis, phenylpropanoids and flavonoids highlights interesting features on how pedicel tissue is regulated by GA3。Although this study did not include measurements of photosynthetic variables, transcriptomics showed that many photosynthesis-related genes were upregulated by GA3treatment. KEGG analysis based on KEGG Orthologs showed several genes involved in photosynthesis reactions (Additional file11: Figure S7A). Non-foliar photosynthesis has been discussed as an important alternative source of carbon acquisition [38]. The green shoots from annual woody species, known as stalks, diminish CO2efflux through photosynthesis in their chlorenchyma cells [39]. Therefore, it may be suggested that photosynthetic positive responses of pedicels to GA3could be favorable in the context of sustaining general growth, although the extent of CO2assimilation and its contribution to pedicel and berry growth remains to be further analyzed. The phenylpropanoids pathway showed positive regulation of genes related to the synthesis of these metabolites, except for those involved in the final steps of monolignol biosynthesis (Additional file11: Figure S7B). These results suggest that phenylpropanoids are upregulated towards the flavonoid pathway in response to GA3(Additional file11: Figure S7C). Polyphenolic compounds derived from plant secondary metabolism, such as phenylpropanoids and flavonoids, include a wide spectrum of structurally and functionally rich compounds, to which a high antioxidant capacity has been attributed [40]. The role of these compounds in the pedicel and its interaction with GA3似乎与一般植物应答se, and further metabolomics analysis could help to identify the variety of secondary metabolites present in this tissue. The abundance of these metabolites is suggested to be higher in treated samples, since polyphenolic compounds interfere with RNA extraction as was observed in this study prior to the optimization of extraction methods for pedicels (data not shown).

Open questions remain about the long-term response of GA3对待组织,如冷tr的影响eatment during storage conditions and their gene activity, which have not yet been addressed under postharvest conditions. Besides, even though a genotype-dependent factor has been suggested in berry drop, there are no genetic studies linking traits such as berry drop in segregating populations, due to the complexity of phenotyping this disorder. Several efforts have focused on increasing the precision of measurement of berry drop [41]. Hopefully, more studies will help to understand the genetic factors and conditions determining berry drop in table grapes.

Conclusions

This is the first study to investigate the effects of GA3treatment on the pedicel transcriptomic response of genotypes with different susceptibility to postharvest berry drop. This phenomenon in table grape is related to GA3由于治疗,增加茎维度to promotion of cell wall and phenylpropanoid metabolism. This is consistent with the differentially expressed genes and gene families found here by a transcriptomic approach. Molecular responses to GA3treatment varied according to the genetic background: genes involved in cell wall synthesis, phenylpropanoid, flavonoid and laccase-related processes were highly upregulated in the L23 susceptible background, which showed a higher incidence of berry drop compared to cv. Thompson Seedless. Positive regulation of these genes could partially explain the differences observed in cell wall yield-related results. Moreover, a strong response of photosynthesis-related genes was detected by gene ontology analysis, suggesting a possible role of the pedicel in non-foliar photosynthesis for both genotypes.

Considering that GA3is a widely used growth regulator for table grape production, the differentially expressed genes identified here could be valuable candidate genes for assisted selection in breeding programs. Future analyses will address the changes in cell wall-related components suggested by the transcriptional data gathered here and could be valuable to characterize the dynamics concerning to this plant growth regulator in the pedicel of table grape.

Methods

Plant material

This study was performed at La Platina Research Center that belongs to the Instituto de Investigaciones Agropecuarias (INIA), located in Santiago, Chile (33°34′20′′S; 70°37′32′′W; 630 m.a.s.l.), during the 2015–2016, 2016–2017 and 2017–2018 seasons. The plant material (plants of 5–10 years old grown on their own roots, conducted with the Spanish trellis system and managed under standard irrigation, fertilization and pest management programs) of two seedless table grape genotypes: cv. Thompson Seedless and L23(‘Line #23’:F1 from cv. Ruby Seedless x cv. Centennial Seedless) were used in this study. Six individual plants from each genotype were used for qPCR experiments; four biological replicates were used for RNAseq experiments. Each plant had treated (GA3) and untreated (Control) clusters to reduce any bias related to individual variation.

Timing of samples was based on the Einchorn-Lorenz phenology system adapted forVitis viniferaL. by Coombe, 1995 [42]. Stage(s) evaluated in each experiment are indicated in each figure.

Regarding the origin of the genotypes assessed in this study, cv. Thompson Seedless (also known as ‘Sultanina’) is a free, ancestral variety, native to west Asia [43]; it was introduced into Chile by the beginning of twentieth century from California, USA. On the other hand, L23 is a segregant belonging to the table grape breeding program of INIA, which was the result of the crossing of cv. Ruby Seedless x Centennial Seedless. Both parents of ‘L23’ are currently free resources originally bred and registered at California Agricultural Experiment Station (University of California) [44] . As with ‘Thompson seedless’, they were introduced into Chile by the mid-twentieth century and kept through vegetative propagation. Cv. Thompson Seedless and L23 are both available from and maintained at INIA grapevine collection.

Growth regulator treatment

Gibberellic acid (GA3, Pro-Gibb 40% Valent Biosciences Chile S.A., Santiago, Chile) treatments of the two genotypes were as follows: (1) for cv. Thompson Seedless, a total of five applications were carried out: 10 ppm GA3for bunch elongation at pre-bloom stage; 15 ppm for thinning at full bloom stage and 200 ppm for berry enlargement, distributed in three doses of 40, 100 and 60 ppm respectively, at EL27, EL29 and EL31 stages, according to Eichhorn-Lorenz phenology system adapted to grapevine by [42])。(2)对于L23,集群被喷three times, only for berry enlargement.:Applications of 40 ppm were at EL27, EL29, and EL31. Dosage for cv. Thompson Seedless and L23 was defined in previous assays, optimizing berry enlargement at the lowest possible value of berry drop, as has been described in [11].

Physiological parameters

The parameters measured in this study were: berry and pedicel diameter (mm), titratable acidity (g L− 1tartaric acid), soluble solids (% w/w g sucrose per 100 g solution), and firmness (g mm− 1) of 10 clusters randomly sampled from different vines. For each variable, 30 healthy (e.g. clusters non-affected by any disease such as powdery mildew) and homogenous berries attached with their cap stems (pedicels) were randomly sampled. The diameter of berries and pedicels was measured using a digital caliper. The soluble solid content of fruit was determined by a refractometer (ATC-1E, Atago, Tokyo, Japan), and was periodically measured to follow the fruit development. The data regarding firmness of berries (considering both skin and flesh) was obtained using a firmness tester (Firmtech II, BioWorks, Wamego, KS). Lastly, titratable acidity was determined by measuring the pooled juice of 10 berries per cluster by titration with 0.1 N NaOH (pH 8.2).

As a part of this evaluation, pedicels from clusters of each individual plant were randomly sampled and pooled (300 ~ 500 pedicels). Then, plant material was immediately frozen in liquid nitrogen and stored at − 80 °C for further molecular analyses (e.g. lignin measurement and transcriptional activity-related assays).

Berry drop

Clusters from each genotype were collected when the average soluble solids content reached 18° on Brix scale (average of ~ 30 berries randomly sampled). For postharvest storage, 10 clusters per treatment (control and GA3-treated) were stored immediately after harvest during 15, 30 and 45 days at 0 °C under standard postharvest management conditions (temperature control, gas control, and application of sulphur dioxide to prevent abiotic deterioration factors such as gray mold) [5),然后贝瑞测量。估计the percentage of berry drop (shattering), 10 clusters per condition were weighed before and after being shaken for 30 s; total weight of detached fruit was registered. Then, berry drop percentage was calculated as follows BD% = [(DB + SB)/TB] / 100. In this formula, DB was the total weight of dropped berries during storage time at 0 °C, SB was the total weight of dropped berries after shaking and TB represents the total weight of the cluster.

Determination of lignin

For lignin estimation, quantitative determination by the acetyl-bromide method [21] was followed with slight modifications. Pedicels were dried at 105 °C for 16 h in an oven and cooled in a vacuum desiccator until the next step. Dry samples (0.3 g) were homogenized in 50 mM potassium phosphate buffer (5 mL, pH 7.0) using an Ultra Turrax T25 homogenizer (IKA-Werke GmbH & Co., Staufen, Germany). The solution was centrifuged (1400×g, 5 min) and the pellet was washed by successive stirring and centrifugation as follows: twice with phosphate buffer (pH 7.0; 7 mL), three times with 1% (v/v) Triton X-100 in pH 50 mM potassium phosphate 7.0 buffer (7 mL), and six times with acetone (5 mL). All supernatants were monitored by spectrophotometric measurements at 280 nm to ensure no contamination with protein and UV-absorbing materials in downstream steps. The pellet was dried in an oven (60 °C, 24 h) and kept cool in a vacuum desiccator until further reaction with acetyl bromide. Samples and a blank solution were prepared as follows: dry matter obtained from the previous step (considered as a protein-free cell wall fraction) was placed in a screw cap 15 mL centrifuge tube (20 mg per sample) containing 0.5 mL 25% acetyl bromide (v/v in glacial acetic acid) and incubated at 70 °C for 30 min. After digestion, samples were cooled in an ice bath and mixed with 0.9 mL of 2 M NaOH and 0.1 mL of 5 M hydroxylamine-HCl. Samples and blank solution were solubilized in 8 mL of glacial acetic acid and absorbance of supernatant was measured at 280 nm immediately after centrifugation (1400×g, 5 mins). Dilution in glacial acetic acid was performed to avoid signal saturation. The standard curve was constructed using alkali lignin (Aldrich 37, 096–7) with a reported absorptivity value of 23.03 g− 1 L cm− 1(Additional file1: Figure S1B). Results were expressed as mg lignin g− 1cell wall as appears in the cited protocol [21].

Dry matter

For the determination of changes in pedicel dry matter accumulation elicited by GA3treatments, samples from different conditions at harvest were measured in an analytical balance. Samples taken from bunches when berries reached an average of 18 °Brix were pooled into a group of 30~60 pedicels and weighed for fresh weight determination (at least 0.3 g fresh weight to avoid technical error). Then replicates were dried at 105 °C in an oven for 16 h and dry weight was measured [21]. Percentage of dry matter and weight of pedicels was determined, considering the number of pedicels weighed in each batch.

RNA isolation

Total RNA was isolated from 1.0 g frozen tissue using a modified hot borate method [45], following all the indications listed in the cited protocol with the only exception that sample tissue was reduced to a third of the recommended amount and polyvinylpyrrolidone was doubled to avoid phenolic compound contamination before RNA isolation. The quantity and quality of RNA were assessed with Qubit® 2.0 fluorometer (Invitrogen™ by Life Technologies, Singapore). Spectrophotometric determination of A260/280and A260/230ratios and electrophoresis in 1.2% formaldehyde-agarose gels verified the quality and integrity of extracted RNA.

Library synthesis and sequencing

For the RNAseq experiment, pooled pedicel samples of individual plants were obtained 7 days after the last GA3application (corresponding to EL31 stage according to the Einhoch-Lorenz phenology system [42], plus 7 days—defined in this study as 7DAT). Two levels were used in the treatment factor, control and GA; and two levels were considered in the genotype factor, cv. Thompson Seedless and L23. Four biological replicates were used for each condition, summing to 16 libraries which were synthesized and sequenced.

Prior to library synthesis, total RNA was assessed by fragment analyzer PROSize® 2.0 version 1.3.1.1 (Advanced Analytical Technologies, Inc., Ames, IA, USA). Mean RNA quality in all samples was 7.35 (SD: 0.56), confirming the integrity of extractions. Then 2.5 μg aliquots were used to isolate poly(A) mRNA for preparation of libraries using TruSeq RNA Sample Prep Kit v2, following the manufacturer’s instructions described in the TruSeq RNA Sample Preparation v2 Guide, Part #15026495 Rev. F (Illumina, Inc.). Libraries were sequenced using the HiSeq 4000 platform (Illumina). Libraries were sequenced as paired-end data (2 × 100 bp).

Pre-processing of reads

Low-quality reads (Phred score Q < 25, nucleotides with undefined base assignmentN> 1 and read length < 50) were removed using the wrapper tool Trim Galore! [46]. Contamination with Illumina adapters was handled with the same tool, removing sequences matching the adapters.

Mapping and matrix count extraction

Alignment of clean reads to theV. viniferaL. PN40024 12X reference genome [47] was performed using STAR software, considering paired-end data default parameters [48]. Uniquely mapped reads were kept for generation of a count matrix using HTseq [49]. The gene model was extracted from the CRIBI database. Genes reported were based on Genoscope structural annotation. Summarized mapping results are given in Additional file12: Table S4.

Data from the count matrix of uniquely mapped reads were visualized to confirm major effects expected from the experimental design. Assessment using multidimensional scaling and hierarchical clustering based on counts per million reads showed no anomalies; differences were given mainly by genotype and treatment factors, with low variability of replicates within groups (for experimental details please see Additional file12Figure S8).

Differential expression analysis

Count matrices were analyzed with EdgeR package [50] under R 3.4.4 software. Correction for composition bias in each sample was handled with the CalcNormFactor option. Then limma-voom transformation of data with was performed with the limma package [25]. Since the number of DEGs was relatively high, an additional filter of log2-fold-change was established. To keep a low false discovery rate, the method described in [22] was followed, considering a prior log2-fold-change of 1. For differential expression analysis considered an adjustedp-value (FDR) < 0.05.

Gene set enrichment analysis and Revigo

Gene set enrichment analysis was conducted on the Agrigo platform v2.0 [51] using hypergeometric distribution of data;p-values were adjusted by the Benjamini-Yeukiteli method and significant terms (FDR < 0.05) were analyzed with the Revigo algorithm to reduce redundancy of terms based on semantic relatedness. The cutoff value used was C: 0.7 [23].

MapMan analysis and KEGG ortholog visualization

Mapping of differentially expressed genes as determined in 2.8 was performed using MapMan software (X4 version). Annotation ofV. viniferaL. was downloaded from BioMart of the Phytozome v12.1 resource page; amino acid sequences were fed into the Mercator 4 online tool, creating an up-to-date mapping file. The Wilcoxon test with adjusted p-value by the Benjamini-Yeukiteli method was used to detect significant differences among BIN categories provided by MapMan [24].

Visualization of KEGG orthologs was accomplished by the usage ofPathviewpackage under R enviroment [52], which is also available as an online utility [53].

cDNA synthesis

In the case of qPCR assays which required the synthesis of cDNA, the reverse transcription was accomplished in 2 μg of total RNA as template using a MMLV-RT reverse transcriptase with previously primed oligo dT (Promega, Madison, WI), according to the manufacturer instructions. The obtained cDNA was then stored at − 20 °C until further qPCR assays.

Real-time qPCR assays

Transcript abundance was assessed by qPCR with the LightCycler® 96 system from Roche (Roche Diagnostics, Mannheim, Germany), as described in [54]. Gene-specific primers (Table2使用在线软件Primer3p)设计lus [55] and were synthesized by IDT (Integrated DNA Technologies). The qPCR assays were performed on six biological samples in duplicate.

The reference gene was GSVIVG01011810001 – probable fructose-bisphosphate aldolase 3 chloroplastic (AFLC3). This gene was identified based on RNAseq results; all non-differentially expressed genes were filtered and ranked according to the lowest variance and coefficient of variation, considering the 16 libraries. Five genes were selected, and primers were synthetized (Additional file13: Table S5). AFLC3 gene values were used in downstream calculations related to relative expression.

Experimental design and statistical analyses

A fully randomized experimental design was used. The two main factors considered were Genotype factor (two levels: cv. Thompson Seedless and L23)and Treatment factor (two levels: Control and GA3).

Data were previously analyzed with Levene test to verify homogeneity of variance assumption (p > 0.05). Then, Analysis of Variance test (ANOVA) and Tukey’s Honest Significance Difference test (Tukey HSD) were conducted to find significant differences between conditions (p < 0.05). These analyses were applied to find significant differences on Fig.1b (analysis was performed on each season, separately. Since that ‘season’ effect could not be fitted to a linear model), Fig.2, Additional file14and Table1。The number of observations (n) can be found on the corresponding description of each figure.

To extract differential expressed genes from RNAseq experiments, data was contrasted against a negative bimomial data distribution through EdgeR package [50], obtainedp-values were applied with the Benjamini-Hochberg method to control false discovery ratio (FDR). As was explained before, since the number of differentially expressed genes (DEGs) was high. The method described by [22] on limma package [25] was applied to consider additional log2fold-change (FC) of the treated over control samples between each genotype with a prior log2FC = 1, p-values were corrected according to Benjamini-Hochberg method and a final FDR < 0.05 was considered to identify DEGs. Please refer to Additional file3to see results regarding this set of statistical analyses.

For gene set enrichment analysis on gene ontology (GO) terms, hypergeometric distribution of data was assumed, and calculated p-values were applied with Benjamini-Yekuiteli method to control false discovery at a rate of 0.05 (FDR < 0.05). From here, revigo analysis was conducted [23)(基于relativ去条款之间的距离e similiraty index). Cutoff value was considered 0.7 in cv. Thompson Seedless and 0.5 in L23.

For mapman [24] results, Wilcoxon test was conducted on assigned categories to detect significant differences p-values were corrected by Benjamini-Hochberg method (FDR < 0.05).

Hierarchal clustering method on heatmap figure (Fig.6)是通过计算欧几里得distance matrix among genes (rows) and among replicates (columns), the quantitative variable was the fragment per kilobase of exon per million reads (FPKM). This quantity was normalized to Z-score for the purpose of comparing the abundance on each cell generated by this method.

Data were mostly handled and analyzed under R software environment [56], unless indicated otherwise.

Availability of data and materials

The datasets used and/or analysed during the current study are available from the corresponding author on reasonable request.

Abbreviations

°C:

Degrees on Celsius scale

4CL:

4-coumarate-CoA ligase

7DAT:

Seven days after GA3-treatment

A260/230:

Ratio of the absorbance at 260 nm and 230 nm

A260/280:

Ratio of the absorbance at 260 nm and 280 nm

AB:

Acetyl bromide

AFLC3:

Probable fructose-bisphosphate aldolase 3 chloroplastic

ANOVA:

Analysis of variance

CAD6:

Cinnamyl-alcohol dehydrogenase 6

CB12:

Chlorophyll a-b binding protein 40 chloroplastic

CCoAOMT:

Caffeoyl-CoA O-methyltransferase

CCR:

Cinnamoyl-CoA reductase

CCR1L:

Cinnamoyl-CoA reductase 1-like

CESA7:

Cellulose synthase A catalytic subunit 7 [UDP-forming]

CHS:

Chalcone synthase

CO2:

Carbon dioxide

DEGs:

Differentially expressed genes

EL17:

Eichhorn-Lorenz phenological stage 17

F3H:

Flavonoid 3′-hydroxylase

FC:

Fold-change

FDR:

False discovery rate

g:

Gram

GA:

Gibberellic acid

GO:

Gene ontology

GUN10:

Endoglucanase 10

KEGG:

Kyoto Encyclopedia of Genes and Genomes

LAC17:

Laccase-17

LAC4:

Laccase-4

mg:

Milligram

mL:

Milliliter

mM:

Milimolar

mm:

Millimeter

MTDH:

Probable mannitol dehydrogenase

NEP1:

Aspartic proteinase nephentesin-1

nm:

Nanometer

PALY:

Phenylalanine ammonia-lyase

PMR5:

Powdery mildew resistance 5

ppm:

Parts per million

qPCR:

Quantitative Polymerase Chain Reaction

RNAseq:

Ribonucleic acids sequencing (whole transcriptome shotgun sequencing)

Ta:

Annealing temperature

TIP11:

Aquaporin tonoplast intrinsic protein 1–1

Tukey HSD test:

Tukey’s Honest Significant Difference test

v/v:

Volume of solute to volume of solution ratio

XSP1:

Xylem serine proteinase 1

References

  1. 1.

    Casanova L, Casanova R, Moret A, Agustí M. The application of gibberellic acid increases berry size of “Emperatriz” seedless grape. Spanish J Agric Res. 2009;7:919.https://doi.org/10.5424/sjar/2009074-1105

    Article谷歌学者

  2. 2.

    Rademacher W. Chemical regulators of gibberellin status and their application in plant production. In: Annual Plant Reviews, vol. 49. Chichester: Wiley; 2016. p. 359–404.https://doi.org/10.1002/9781119210436.ch12

    Chapter谷歌学者

  3. 3.

    Zoffoli JP, Latorre英航Naranjo p .收获前的达成ications of growth regulators and their effect on postharvest quality of table grapes during cold storage. Postharvest Biol Technol. 2009;51:183–92.https://doi.org/10.1016/j.postharvbio.2008.06.013

    CASArticle谷歌学者

  4. 4.

    Zoffoli JP, Latorre BA. Table grape (Vitis viniferaL.). In: Yahia EMBT-PB and T of T and SF, editor. Postharvest Biology and Technology of Tropical and Subtropical Fruits. Elsevier; 2011. 179–212. doi:https://doi.org/10.1533/9780857092885.179

    Chapter谷歌学者

  5. 5.

    Zoffoli JP. Postharvest handling of table grape. In: Acta Horticulturae; 2008. p. 415–9.

    谷歌学者

  6. 6.

    Ben-Tal Y. Effects of Gibberellin Treatments on Ripening and Berry Drop from Thompson Seedless Grapes. Am J Enol Vitic. 1990;41:142 LP–146http://www.ajevonline.org/content/41/2/142.abstract

    谷歌学者

  7. 7.

    Nakamura M, Hori Y. Postharvest berry drop of seedless berries produced by GA treatment in grape cultivar “Kyoho”. I. Relationship between postharvest berry drop and rachis hardness. Tohoku J Agric Res. 1981;32:1–13.

    CAS谷歌学者

  8. 8.

    Murcia G, Pontin M, Reinoso H, Baraldi R, Bertazza G, Gómez-Talquenca S, et al. ABA and GA 3 increase carbon allocation in different organs of grapevine plants by inducing accumulation of non-structural carbohydrates in leaves, enhancement of phloem area and expression of sugar transporters. Physiol Plant. 2016;156:323–37.https://doi.org/10.1111/ppl.12390

    CASArticlePubMed谷歌学者

  9. 9.

    Nakamura M, Hori Y. Postharvest berry drop of postharvest berry drop of seedless berries produced by GA treatment in grape cultivar “Kyoho” III. Effect of GA on xylem differentiation and lignification of rachis. Tohoku J Agric Res. 1984;34:55–64.

    CAS谷歌学者

  10. 10.

    Retamales J, Cooper T. Berry drop and fruit removal forces as related with GA3 applications in table grape. In: Seventh International Symposium on Plant Growth Regulators in Fruit Production: Jerusalem, Israel, June 14–19, 1992; 1992. p. 81–3.https://doi.org/10.13140/2.1.4192.7680

    Chapter谷歌学者

  11. 11.

    García-Rojas M, Meneses M, Oviedo K, Carrasco C, Defilippi B, González-Agüero M, et al. Exogenous gibberellic acid application induces the overexpression of key genes for pedicel lignification and an increase in berry drop in table grape. Plant Physiol Biochem. 2018;126:32–8.https://doi.org/10.1016/j.plaphy.2018.02.009

    CASArticlePubMed谷歌学者

  12. 12.

    Eriksson ME, Israelsson M, Olsson O, Moritz T. Increased gibberellin biosynthesis in transgenic trees promotes growth, biomass production and xylem fiber length. Nat Biotechnol. 2000;18:784–8.https://doi.org/10.1038/77355

    CASArticlePubMed谷歌学者

  13. 13.

    Biemelt S, Tschiersch H, Sonnewald U. Impact of altered gibberellin metabolism on biomass accumulation, lignin biosynthesis, and photosynthesis in transgenic tobacco plants. Plant Physiol. 2004;135:254–65.https://doi.org/10.1104/pp.103.036988

    CASArticlePubMedPubMed Central谷歌学者

  14. 14.

    Cho J, Jeon H, Kim M, Vo TK, Kim J, Park E, et al. Wood forming tissue-specific bicistronic expression of Pd GA 20ox1 and Ptr MYB 221 improves both the quality and quantity of woody biomass production in a hybrid poplar. Plant Biotechnol J. 2019;17:1048–57.https://doi.org/10.1111/pbi.13036

    CASArticlePubMed谷歌学者

  15. 15.

    Jeon HW, Cho JS, Park EJ, Han KH, Choi YI, Ko JH. Developing xylem-preferential expression of PdGA20ox1, a gibberellin 20-oxidase 1 from Pinus densiflora, improves woody biomass production in a hybrid poplar. Plant Biotechnol J. 2016;14:1161–70.

    CASArticle谷歌学者

  16. 16.

    Dal Santo S, Tornielli GB, Zenoni S, Fasoli M, Farina L, Anesi A, et al. The plasticity of the grapevine berry transcriptome. Genome Biol. 2013;14:r54.https://doi.org/10.1186/gb-2013-14-6-r54

    CASArticlePubMedPubMed Central谷歌学者

  17. 17.

    Fasoli M, Dal Santo S, Zenoni S, Tornielli GB, Farina L, Zamboni A, et al. The grapevine expression atlas reveals a deep Transcriptome shift driving the entire plant into a maturation program. Plant Cell. 2012;24:3489–505.https://doi.org/10.1105/tpc.112.100230

    CASArticlePubMedPubMed Central谷歌学者

  18. 18.

    Grimplet J, Tello J, Laguna N, Ibáñez J. Differences in Flower Transcriptome between Grapevine Clones Are Related to Their Cluster Compactness, Fruitfulness, and Berry Size. Front Plant Sci. 2017;8:1–17.https://doi.org/10.3389/fpls.2017.00632

    Article谷歌学者

  19. 19.

    Wang C, Xin M, Zhou X, Liu W, Liu D, Qin Z. Transcriptome profiling reveals candidate genes associated with sex differentiation induced by night temperature in cucumber. Sci Hortic (Amsterdam). 2018;232:162–9.https://doi.org/10.1016/j.scienta.2017.12.018

    CASArticle谷歌学者

  20. 20.

    Chai L, Li Y, Chen S, Perl A, Zhao F, Ma H. RNA sequencing reveals high resolution expression change of major plant hormone pathway genes after young seedless grape berries treated with gibberellin. Plant Sci. 2014;229:215–24.https://doi.org/10.1016/j.plantsci.2014.09.010

    CASArticlePubMed谷歌学者

  21. 21.

    Moreira-Vilar FC, Siqueira-Soares RDC, Finger-Teixeira A, de Oliveira DM, Ferro AP, da Rocha GJ, et al. The acetyl bromide method is faster, Simpler and Presents Best Recovery of Lignin in Different Herbaceous Tissues than Klason and Thioglycolic Acid Methods. PLoS One. 2014;9:e110000.https://doi.org/10.1371/journal.pone.0110000

    CASArticlePubMedPubMed Central谷歌学者

  22. 22.

    McCarthy DJ, Smyth GK. Testing significance relative to a fold-change threshold is a TREAT. Bioinformatics. 2009;25:765–71.https://doi.org/10.1093/bioinformatics/btp053

    CASArticlePubMedPubMed Central谷歌学者

  23. 23.

    Supek F, Bošnjak M, Škunca N, Šmuc T. REVIGO summarizes and visualizes long lists of gene ontology terms. PLoS One. 2011;6:e21800.https://doi.org/10.1371/journal.pone.0021800

    CASArticlePubMedPubMed Central谷歌学者

  24. 24.

    Usadel B, Poree F, Nagel A, Lohse M, Czedik-Eysenberg A, Stitt M. A guide to using MapMan to visualize and compare Omics data in plants: a case study in the crop species, maize. Plant Cell Environ. 2009;32:1211–29.https://doi.org/10.1111/j.1365-3040.2009.01978.x

    Article谷歌学者

  25. 25.

    Ritchie ME, Phipson B, Wu D, Hu Y, Law CW, Shi W, et al. Limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Res. 2015;43.https://doi.org/10.1093/nar/gkv007

    Article谷歌学者

  26. 26.

    Hedden P, Thomas SG. Annual plant reviews, volume 49. Chichester: Wiley; 2016.https://doi.org/10.1002/9781119210436

    Book谷歌学者

  27. 27.

    Zoffoli JP, Latorre BA. Table grape (Vitis viniferaL .) Cambridge: Woodhead Publishing Limited;

  28. 28.

    Nakamura M, Hori Y. Postharvest berry drop of seedless berries produced by GA treatment in grape cultivar “Kyoho” IV. Rachis hardness and phenol production of rachis in vitro following GA application as affected by the addition of β-Nitrostyrene derivatives. Tohoku J Agric Res. 1985;35:81–9.

    CAS谷歌学者

  29. 29.

    Barros J, Serk H, Granlund I, Pesquet E. The cell biology of lignification in higher plants. Ann Bot. 2015;115:1053–74.https://doi.org/10.1093/aob/mcv046

    CASArticlePubMedPubMed Central谷歌学者

  30. 30.

    Tobimatsu Y, Schuetz M. Lignin polymerization: how do plants manage the chemistry so well? Curr Opin Biotechnol. 2019;56:75–81.https://doi.org/10.1016/j.copbio.2018.10.001

    CASArticlePubMed谷歌学者

  31. 31.

    Donaldson LA. Lignification and lignin topochemistry — an ultrastructural view. Phytochemistry. 2001;57:859–73.https://doi.org/10.1016/S0031-9422(01)00049-8

    CASArticlePubMed谷歌学者

  32. 32.

    Hou X, Wei L, Xu Y, Khalil-Ur-rehman M, Feng J, Zeng J, et al. Study on russet-related enzymatic activity and gene expression in ‘shine Muscat’ grape treated with ga3 and cppu. J Plant Interact. 2018;13:195–202.https://doi.org/10.1080/17429145.2018.1455904

    CASArticle谷歌学者

  33. 33.

    Ferrara G, Mazzeo A, Netti G, Pacucci C, Matarrese AMS, Cafagna I, et al. Girdling, gibberellic acid, and forchlorfenuron: effects on yield, quality, and metabolic profile of table grape cv. Italia. Am J Enol Vitic. 2014;65:381–7.

    CASArticle谷歌学者

  34. 34.

    Prozil SO, Evtuguin DV, Silva AMS, Lopes LPC. Structural characterization of lignin from grape stalks ( Vitis vinifera L.). J Agric Food Chem. 2014;62:5420–8.https://doi.org/10.1021/jf502267s

    CASArticlePubMed谷歌学者

  35. 35.

    Cosgrove DJ. Growth of the plant cell wall. Nat Rev Mol Cell Biol. 2005;6:850–61.https://doi.org/10.1038/nrm1746

    CASArticle谷歌学者

  36. 36.

    Fernández JL, Sáez F, Castro E, Manzanares P, Ballesteros M, Negro MJ. Determination of the Lignocellulosic components of olive tree pruning biomass by near infrared spectroscopy. Energies. 2019;12:2497.https://doi.org/10.3390/en12132497

    CASArticle谷歌学者

  37. 37.

    Kim WC, Kim JY, Ko JH, Kang H, Han KH. Identification of direct targets of transcription factor MYB46 provides insights into the transcriptional regulation of secondary wall biosynthesis. Plant Mol Biol. 2014;85:589–99.

    CASArticle谷歌学者

  38. 38.

    Aschan G, Pfanz H. Non-foliar photosynthesis - a strategy of additional carbon acquisition. Flora. 2003;198:81–97.

    Article谷歌学者

  39. 39.

    Pfanz H, Aschan G, Langenfeld-Heyser R, Wittmann C, Loose M. Ecology and ecophysiology of tree stems: Corticular and wood photosynthesis. Naturwissenschaften. 2002;89:147–62.

    CASArticle谷歌学者

  40. 40.

    Quideau S, Deffieux D, Douat-Casassus C, Pouységu L. Plant polyphenols: chemical properties, biological activities, and synthesis. Angew Chemie Int Ed. 2011;50:586–621.https://doi.org/10.1002/anie.201000044

    CASArticle谷歌学者

  41. 41.

    Vinokur Y, Rodov V, Levi A, Kaplunov T, Zutahy Y, Lichter A. A method for evaluating fruit abscission potential of grapes and cherry tomato clusters. Postharvest Biol Technol. 2013;79:20–3.https://doi.org/10.1016/j.postharvbio.2012.12.012

    CASArticle谷歌学者

  42. 42.

    Coombe BG. Growth stages of the grapevine: adoption of a system for identifying grapevine growth stages. Aust J Grape Wine Res. 1995;1:104–10.https://doi.org/10.1111/j.1755-0238.1995.tb00086.x

    Article谷歌学者

  43. 43.

    Aradhya MK, Dangl GS, Prins BH, Boursiquot J-M, Walker MA, Meredith CP, et al. Genetic structure and differentiation in cultivated grape, Vitis vinifera L. Genet Res. 2003;81:179–92.

    CASArticle谷歌学者

  44. 44.

    Maul et al. Vitis International Variety Catalogue. 2019.www.vivc.de。Accessed 13 Dec 2019.

    谷歌学者

  45. 45.

    Gudenschwager O, González-Agüero M, Defilippi BG. A general method for high-quality RNA isolation from metabolite-rich fruits. South Afr J Bot. 2012;83:186–92.https://doi.org/10.1016/j.sajb.2012.08.004

    CASArticle谷歌学者

  46. 46.

    Martin M. Cutadapt removes adapter sequences from high-throughput sequencing reads. EMBnet J. 2011;17:10.https://doi.org/10.14806/ej.17.1.200

    Article谷歌学者

  47. 47.

    Jaillon O, Aury J-M, Noel B, Policriti A, Clepet C, Casagrande A, et al. The grapevine genome sequence suggests ancestral hexaploidization in major angiosperm phyla. Nature. 2007;449:463–7.https://doi.org/10.1038/nature06148

    CASArticlePubMed谷歌学者

  48. 48.

    Dobin A, Davis CA, Schlesinger F, Drenkow J, Zaleski C, Jha S, et al. STAR: ultrafast universal RNA-seq aligner. Bioinformatics. 2013;29:15–21.https://doi.org/10.1093/bioinformatics/bts635

    CASArticle谷歌学者

  49. 49.

    Anders S, Pyl PT, Huber W. HTSeq--a Python framework to work with high-throughput sequencing data. Bioinformatics. 2015;31:166–9.https://doi.org/10.1093/bioinformatics/btu638

    CASArticlePubMedPubMed Central谷歌学者

  50. 50.

    Robinson MD, McCarthy DJ, Smyth GK. edgeR: a bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics. 2010;26:139–40.https://doi.org/10.1093/bioinformatics/btp616

    CASArticlePubMedPubMed Central谷歌学者

  51. 51.

    Tian T, Liu Y, Yan H, You Q, Yi X, Du Z, et al. agriGO v2.0: a GO analysis toolkit for the agricultural community, 2017 update. Nucleic Acids Res. 2017;45:W122–9.https://doi.org/10.1093/nar/gkx382

    CASArticlePubMedPubMed Central谷歌学者

  52. 52.

    Luo W, Brouwer C. Pathview: an R/bioconductor package for pathway-based data integration and visualization. Bioinformatics. 2013;29:1830–1.https://doi.org/10.1093/bioinformatics/btt285

    CASArticlePubMedPubMed Central谷歌学者

  53. 53.

    Luo W, Pant G, Bhavnasi YK, Blanchard SG, Brouwer C. Pathview web: user friendly pathway visualization and data integration. Nucleic Acids Res. 2017;45:W501–8.https://doi.org/10.1093/nar/gkx372

    CASArticlePubMedPubMed Central谷歌学者

  54. 54.

    García-Rojas M, Gudenschwager O, Defilippi BG, González-Agüero M. Identification of genes possibly related to loss of quality in late-season ‘Hass’ avocados in Chile. Postharvest Biol Technol. 2012;73:1–7.https://doi.org/10.1016/j.postharvbio.2012.05.007

    CASArticle谷歌学者

  55. 55.

    Untergasser A, Nijveen H, Rao X, Bisseling T, Geurts R, Leunissen JAM. Primer3Plus, an enhanced web interface to Primer3. Nucleic Acids Res. 2007;35 Web Server:W71–4.https://doi.org/10.1093/nar/gkm306

    Article谷歌学者

  56. 56.

    R Core Team. R: A language and environment for statistical computing. 2019.https://www.r-project.org/

    谷歌学者

Download references

Acknowledgements

Authors sincerely thanks Macarena Bastías for her assistance on library synthesis. MM thanks all the INIA – La Platina Department of Postharvest and Biotechnology for its critical role on planning and discussing experimental design and techniques. Authors also thank CM research group for their expertise on the analysis of transcriptomics data. MM thanks to Anibal Riveros for his advice on transcriptomics-related figures. MM also thanks Pedro Martínez Gómez for his guidance on the thesis related to this research paper.

Funding

Funding to support this research was received from Genoma-Chile, FONDEF program, grant G131–0003. MM acknowledges Programa de Formación de Capital Humano Avanzado CONICYT-PCHA/Doctorado Nacional/2016–21160793 fellowship for the financial support on the PhD program. Funding bodies did not take participation in the design of the study, collection, analysis and interpretation of data nor in the writing of this manuscript.

Author information

Affiliations

Authors

Contributions

MM, MGR, PH and BD planned and designed the research. MM and MGR performed the experiments and analyzed the data concerning phenotyping, lignin measurement and qPCR. MM and TCV analyzed the transcriptomic data. MGA assisted on qPCR-related experiments. MM and PH wrote the manuscript. CM, CME, RI, PH contributed to the critical discussion, interpretation of data and edited the manuscript. All authors have read and approved the final manuscript.

Corresponding author

Correspondence toPatricio Hinrichsen

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary information

Additional file 1: Figure S1.

Validation of AB method [21] for analysis of lignin content in grapevine pedicel.APedicels increases its fresh and dry matter in response to GA3-treatment. Representative images from pedicels sampled of L23 genotype at harvest time (17.4° Brix). Pedicels were cut and imaged before and after being dried at 105 °C by 16 h on oven to illustrate differences between treated and non-treated groups regarding fresh and dried condition.BCalibration curve of standard Lignin Alkali. Catalog No. Aldrich 37, 096–7.CPreliminary report of soybean seed coat performed on this research compared to the value given in literature (n = 16).

Additional file 2: Figure S2.

Pipeline for RNAseq analysisScheme followed for extraction of differentially expressed genes and gene ontology analysis from paired-end sequencing Illumina reads. Green hexagons show programs and some relevant features and/or filters considered on each step.

Additional file 3: Figure S3.

Differentially expressed genes before and after testing relative to a threshold method [22]. Mean-Difference plots are shown on each figureAResponse to GA3in L23 genotype (FDR < 0.05)BResponse to GA3in L23 genotype, considering an additional parameter of prior log-fold-change equal to 1 (FDR < 0.05)CResponse to GA3in cv. Thompson Seedless (FDR < 0.05)DResponse to GA3in cv. Thompson Seedless, considering an additional parameter of prior log-fold-change equal to 1 (FDR < 0.05). The number of DE genes is specified on the top right section of each graph.

Additional file 4: Data S1.

Results for differential expression analysis conducted on L23 and cv. Thompson Seedless.

Additional file 5: Data S2.

Results for single enrichment analysis on GO categories for L23 and cv. Thompson Seedless.

Additional file 6: Figure S4.

Reduction and visualization of data from gene set enrichment analysis on overexpressed genes from Thompson Seedless. List of significantly enriched GO terms (FDR < 0.05) was used as input for revigo analysis [23]. Distances between terms was based on simrelindex. Cutoff value: 0.5.

Additional file 7: Table S1.

GO terms obtained from L23 filtered according to revigo method proposed by [23].

Additional file 8: Table S2.

GO terms obtained from cv. Thompson Seedless filtered according to revigo method proposed by [23].

Additional file 9: Table S3.

Primer list specifying the genes used to corroborate results given by the transcriptomic approachEighteen primers were designed to obtain log2-fold-change values. The list details sequences for each gene, gene ID and the corresponding description based on BLAST results of the annotated sequence.Figure S5.Correlation between RNAseq and qPCR methods on the measurements of transcript abundance. A significant correlation based on Pearson correlation test (p = 1.368e-13) was detected, which was assayed on pedicel samples obtained from 7DAT stage. Eighteen genes were analyzed for both genotypes (cv. Thompson Seedless and L23), dots show the directional response considering fold-changes of averaged values of treated samples over control samples.

Additional file 10: Figure S6.

Mapman figure of DE genes from L23. Mapping of differentially expressed genes from L23 genotype under GA treatment condition is significantly related to promotion of cell wall-related processes. The DE genes fromL23were visualized through MapMan software. Blue dots show overexpressed genes while red dots account repressed genes by GA treatment. Grey dots correspond to filtered genes with non-differential expression. The templates for displaying genes were taken from Mapman [24].

Additional file 11: Figure S7.

Prominent pathways observed to be elicited by GA3on L23 genotype susceptible to berry drop. Three KEGG maps based on KEGG orthologs forVitis viniferaL. are presented in this figure: PhotosynthesisA, phenylpropanoidsB, flavonoidsC。Only transcriptional data is shown and is based on the observed log2fold-change ratios of GA3treated over control samples of differentially expressed genes. The templates for displaying pathways were taken from Pathview [52,53].

Additional file 12: Table S4.

Summarized results obtained from mapping to reference genome PN40024.Figure S8.Effects on transcriptomic profiles were given mainly by genotype and treatment factors. Two unsupervised methods were performed to identify grouping patterns of samples based on the counts per million reads data obtained after trimming step, before the differential expression analysis.AMultidimensional scaling plot from the Euclidean distance matrix calculated from count per million reads by library is shown. Each dot is a library (n = 16 libraries).BHierarchical clustering of libraries with heatmap of the 500 most variable genes across samples. Genes were ranked according to its variance from log-counts data, the 500 most variable were extracted and hierarchical clustering method was performed according to Euclidean distance matrix calculated from such set. Rows represents genes and columns represents libraries. Orange corresponds to cv. Thompson Seedless and purple corresponds to L23 genotype. On topleft section, color key and a histogram of gene profile expression is shown. Regarding bottom labels for each column: TC (cv. Thompson Seedless-Control), TG (cv. Thompson Seedless-GA), LC (L23-Control), LG (L23-GA), the numbers is the identifier assigned to each biological replicate.

Additional file 13: Table S5.

Primer list of reference genes assayed on this study to measure transcriptional expression in pedicel.

Additional file 14:

Gene expression changes in genes proposed as candidate biomarkers reveal strong regulation by GA3 treatment. Mean and standard deviation are shown on each bar and error bar, respectively (n= 6). The reference gene used in this study was AFLC3 (GSVIVG01011810001 – possible fructose-bisphosphate aldolase 3 chloroplastic), identified from RNAseq data. Common letters on top of each bar indicate no significant differences between conditions and were identified based on Tukey’s post-hoc test on log2-transformed expression data (p< 0.05).

Rights and permissions

Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/),允许无限制的使用、分配nd reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Meneses, M., García-Rojas, M., Muñoz-Espinoza, C.et al.Transcriptomic study of pedicels from GA3-treated table grape genotypes with different susceptibility to berry drop reveals responses elicited in cell wall yield, primary growth and phenylpropanoids synthesis.BMC Plant Biol20,66 (2020). https://doi.org/10.1186/s12870-020-2260-6

Download citation

Keywords

  • Gibberellic acid
  • Berry drop
  • RNAseq
  • Vitis viniferaL.
  • Shattering
  • Postharvest, pedicel