Skip to main content
  • Research article
  • Open access
  • Published:

MaMADS2 repression in banana fruits modifies hormone synthesis and signalling pathways prior to climacteric stage

Abstract

Background

While the role of ethylene in fruit ripening has been widely studied, the contributions of additional plant hormones are less clear. Here we examined the interactions between the transcription factor MaMADS2-box which plays a major role in banana fruit ripening and hormonal regulation. Specifically, we used MaMADS2 repressed lines in transcriptome and hormonal analyses throughout ripening and assessed hormone and gene expression perturbations as compared to wild-type (WT) control fruit.

Results

Our analyses revealed major differences in hormones levels and in expression of hormone synthesis and signaling genes mediated by MaMADS2 especially in preclimacteric pulp. Genes encoding ethylene biosynthesis enzymes had lower expression in the pulp of the repressed lines, consistent with reduced ethylene production. Generally, the expression of other hormone (auxin, gibberellins, abscisic acid, jasmonic acid and salicylic acid) response pathway genes were down regulated in the WT pulp prior to ripening, but remained high in MaMADS2 repressed lines. Hormone levels of abscisic acid were also higher, however, active gibberellin levels were lower and auxin levels were similar with MaMADS2 repression as compared to WT. Although abscisic level was higher in MaMADS2 repression, exogenous abscisic acid shortened the time to ethylene production and increased MaMADS2 mRNA accumulation in WT. Exogenous ethylene did not influence abscisic acid level. CRE - a cytokinin receptor, increased its expression during maturation in WT and was lower especially at prebreaker in the repressed line and zeatin level was lower at mature green of the repressed line in comparison to WT.

Conclusions

In addition to previously reported effects of MaMADS2 on ethylene, this transcription factor also influences other plant hormones, particularly at the pre-climacteric stage. The cytokinin pathway may play a previously unanticipated role via MaMADS2 in banana ripening. Finally, abscisic acid enhances MaMADS2 expression to promote ripening, but the transcription factor in turn auto inhibits ABA synthesis and signaling. Together, these results demonstrate a complex interaction of plant hormones and banana fruit ripening mediated by MaMADS2.

Background

Fruits of plants in the genus Musa including bananas and plantains represent staple foods for millions of people, especially in developing countries, and an important carbohydrate and nutrient source for billions more, the world over. Although multiple genetic components of ripening control have been discovered, especially in the model plant tomato [1,2,3], the mechanism in banana is less characterized, though critical transcription factors have recently been discovered [4, 5]. Banana, as typical of climacteric fruit, exhibit an increase in respiration and a surge in ethylene production concomitant with ripening [6, 7].

The involvement of ethylene in ripening has been established, but less so the role of other hormones: indole-3-acetic acid (IAA or auxin), cytokinin (CK), abscisic acid (ABA), gibberellins (GAs), jasmonic acid (JA) and salicylic acid (SA) [8,9,10]. Auxin is involved in early fruit development, but its role in ripening is less clear. Contradicting conclusions as a negative or positive regulator were reported based on its exogenous application [11], including in banana [12] and from its measurements or gene-related expression [13, 14]. By analyzing the function of ARF genes it was concluded that auxin may play a complex role in fruit ripening possibly affecting other hormones [15, 16]. Moreover, it has been shown in tomato that auxin may specifically intersect with the ethylene response pathway at the point of EIN2 [17]. So, while the role of auxin remains ambiguous and may vary based on development, tissue or species, available evidence, nevertheless, indicates involvement in fruit ripening.

Gibberellins (GAs) have been mostly associated with fruit growth [18], with little evidence regarding a specific role in fruit ripening. Interestingly, exogenous GA3 did not alter the respiration or the ethylene profile of banana, but delayed starch breakdown and sucrose accumulation, possibly by reducing levels of sucrose phosphate synthase [19].

ABA most likely play a positive role in ripening of climacteric and non-climacteric fruit [20]. Abscisic acid (ABA) was found to increase prior to the increase of ethylene in tomato [21]. ABA levels and the expression of 9-cis-epoxycarotenoid dioxygenase (NCED), encoding a critical enzyme in ABA biosynthesis, increased during fruit maturation, and declined during ripening as ethylene increased [22]. External application of ABA to tomato and banana promoted ripening [22, 23], and this effect is consistent with ethylene induction [8, 24]. Silencing tomato SlNCED decreased fruit softening and increased shelf life, but increased ethylene production [25], supporting a complex role in both ethylene synthesis and response.

Several studies suggest that jasmonates (JAs) might be positive regulators of fruit ripening through induced expression of ethylene synthesis pathway genes [18, 26, 27]. Salicylic acid (SA) levels have not been determined during ripening, however application of exogenous SA to different fruits including banana, reduced respiration and ethylene production and decreased cell wall deterioration [28, 29]. Finally, cytokinins (CKs) are usually associated with delayed senescence, cell death [30] and fruit ripening [31], and their levels have been reported to decrease with tomato development [18, 32]. However, tomato plants overexpressing the isopentenyl transferase (ipt) gene, encoding a key enzyme in the CKs biosynthesis, presented a sectored ripening phenotype [33] and in grapes, the levels of CK increased following veraison and remained high during ripening [34], Hence, the role of CKs in ripening is still obscure. While some evidence exists for multiple hormones having effects on ripening, the molecular mechanism(s) through which they coordinate and balance their effects and their interactions with the ripening-related transcription factors remain uncertain.

We have previously shown that transgenic MaMADS2 repressed fruits are delayed in ripening [4] and thus, may present a system for addressing hormone roles in the ripening process. Here we characterize in wild type (WT) and MaMADS2 repressed lines the transcriptomes in peel and pulp through fruit ripening and hormones levels in pulp. Our analysis revealed that hormone synthesis and signaling pathways, and hormone levels have been altered in response to MaMADS2 repression prior to normal ethylene induction, suggesting a role for MaMADS2 in coordinating hormone cross-talk at preclimacteric stage.

Methods

Plant material and sample collections

Banana (Musa acuminate, AAA Cavendish subgroup, cv. Grand Nain) WT and MaMADS2 repressed lines (antisense line #44 and RNAi line #32) were grown in the Jordan Valley, and fruits were collected for transcriptome, hormone analysis or ABA and ethylene treatment. Fruit were harvested at 75% of final cross-sectional filling and fruit maturity, determined by measuring the average of the locular angles of banana fruit and the ratio of peel to pulp in cross sections [35]. Hands of the first, second or third tiers, containing 10–30 fruit were separated from the bunch and stored at 20 °C and 95% RH. Samples of both control and transformed banana were collected every 2–3 days and color, ethylene (C2H4) production and carbon dioxide (CO2) emission were determent. Fruits of different ripening stages: mature green (MG; 1 day after harvest), prebreaker (PB; 5(WT), 9(#44) and 10(#32) days after harvest) just before the ethylene peak, breaker (B; 7(WT), 12(#44) and 17(#32) days after harvest) and ripe (R; 10(WT), 17(#44) and 20 (#32) days after harvest) with brown spots were analyzed in RNAseq. The stages were determined by the color of the peel using Minolta. The average of hue (o) from at least two replicates for WT, MaMADS2 antisense (#44) and MaMADS2 RNAi (#32), were: MG:118, 119.3, 119.3; PB:106.3, 108, 113.3; B:100.5, 96.1, 99.1; R: 86.2, 88.4, 89, respectively.

ABA and ethylene treatments

Slices of WT banana fingers taken from MG stage were treated with ABA. Slices were placed on a Whatman paper immersed with 10 μM ABA or 0.01% methanol in a petri dish. Ethylene production rate was followed periodically by transferring the slices to a 120 ml sealed jar, and mRNA was extracted from three days-treated and non-treated slices. For ethylene treatment, slices of WT banana fingers were taken from MG and PB and placed in a 120 ml sealed glass jars at 20 °C with 1 ppm ethylene for 24 h. ABA was extracted by acidic methanol extraction as described [36].

Ripening measurements

Ethylene and carbon dioxide production were determined by GC head-space gas analysis of individual banana fingers in 2-L sealed glass jars at 20 °C as described [5]. Peel color was determined from surface area of at least two individual banana fingers using a Minolta CR-300 colorimeter (Minolta Corporation, New Jersey, USA). Total soluble solids (TSS) of peel and pulp juice, resulting from freezing and thawing of the tissues, was determined using a handheld HSR-500 refractometer (Atago Co. Ltd., Japan).

Transcriptome analysis

The transcriptomes of WT, as well as, MaMADS2 repressed lines were determined by RNAseq analysis. Total RNA was extracted using the Spectrum™ Plant Total RNA Kit (Sigma Aldrich) and preparation of samples for RNAseq library construction is described along with procedures for sequencing and sequence data analysis with levels of expression determined as Reads per Kilo base of transcript per Million mapped reads (RPKM) [37]. As an independent secondary confirmation of RNAseq data, several genes were also analyzed by Quantitative reverse transcriptase PCR (qRT-PCR).

Comparative transcriptome analysis and differential gene expression

Gene expression was normalized and differential expression was calculated using the R Bioconductor [38] “DESeq” package [39]. The results are displayed in BaseMean (mean normalized counts, averaged over all replicates). The transcript reference of banana was downloaded from the banana genome hub website (http://banana-genome.cirad.fr/) [40] and were further annotated by blasting the transcripts against the NCBI non-redundant protein database [41] with the aid of blast2go software [42]. The translation of chromosome location of each of the genes to genes number is summarized in Additional file 1: Table S1. The correlation among all the conditions (peel/pulp, stages of ripening: MG/PB/B/R, genotypes: WT/#44/#32) is calculated based on the normalized count table with all the average expression of all the genes and not only the differentially expressed. The correlation takes into consideration the array of expression matrix of every gene in each condition and the distance matrix is calculated. The correlation between gene expressions in different developmental stages was calculated with “dist” function in R applying the “Canberra” method. The distance matrix was visualized using gplots R package. Functional analysis of the differentially modified genes (FDR < 0.1) was performed by Ontologizer (Parent-Child-Union, Benjamini-Hochberg) to detect enrichment of GO (Gene Ontology) http://ontologizer.de/ [43].

Hormones profiling

The hormonal profile of WT, as well as MaMADS2 repressed line (#32) fruit tissues, was determined by UHPLC-MS/MS as described [44]. Samples were collected at several developmental stages: MG (1 day after harvest), PB (5(WT) and 14(#32) days after harvest) and B (10(WT) and 20 (#32) days after harvest). The pulp of six banana fingers was analyzed for each stage. The stages were determined by the color of the peel using Minolta (see above). The average of hue (o) for WT and MaMADS2 repressed (#32), were: MG:120.3, 120.9; PB:108.9, 111; B:98.6, 99.2, respectively.

Hormonometer analysis

The Hormonometer was developed as a tool in Arabidopsis thaliana to identify the transcriptome footprint of several hormones and their temporal induction [45]. We have identified 10,606 genes with homology of more than 60% identity to Arabidopsis thaliana genes. This “Arabidopsis homologs gene pool” was used for the Hormonometer analysis. The fold change in expression between transgenic lines and WT, as well as, the p-value was loaded to Hormonometer tool (https://hormonometer.weizmann.ac.il/hormb/about) and the output gave the correlation values between the query gene pool and the transcriptome footprint for each hormone [45]. High positive values for a specific hormone indicates that there is high representation of the specific hormone-induced genes within the query gene pool.

Results

Peel and pulp transcriptome diverged during ripening and MaMADS2 repression delayed maturation

To better assess the role of MaMADS2 during ripening, we performed RNA-seq transcriptome analysis on WT and MaMADS2 repressed fruit at the same physiological stages based on peel color. This approach enabled assessment of whether the underlying transcriptomes were the same or different at each stage which would in turn, indicates whether most or a subset of ripening phenomena were influenced by MaMADS2. The transcriptomes of WT, and the MaMADS2 repressed lines; MaMADS2 RNAi (#32) and antisense (#44) lines were assessed throughout post-harvest ripening at the mature-green (MG; average H°: 118, 119.3, 119.3 respectively), pre-breaker (PB; average H°: 106.3, 113.3, 108 respectively), breaker (B; average H°: 100.5, 99.1, 96.1 respectively) and ripe (R; average H°: 86.2, 89, 88.4 respectively) stages.

Bananas of lines #32 and #44 (Fig. 1) were late in all ripening parameters examined: CO2, ethylene production, color and TSS. The ethylene peak in WT fruits occurred 7 days after harvest, and in lines #32 and #44 it was delayed by 10 and 5 days, respectively (Fig. 1A). CO2 levels increased 5 days after harvest in WT fruits, but it was delayed by 7 days in the transgenic fruits (Fig. 1B). As with ethylene levels, WT fruits became yellow 7 days after harvest (100.5 ± 3.5H°), while the yellow color was delayed by 10 and 5 days in lines #32 and #44, respectively (Fig. 1C).

Fig. 1
figure 1

Ethylene production, respiration rate, color and total soluble solids (TSS) in WT and MaMADS2 repressed banana fruits after harvest. Ethylene production rate (a), CO2 production rate (b), color expressed as hue angle (o) (c) and total soluble solids expressed as brix (%) (d) were measured in fruits of the second hand of WT, MaMADS2 RNAi line (#32) and MaMADS2 antisense line (#44) during ripening. For each time point at least three fruits were measured and average and standard error plotted. Arrow in panel A indicate the samples taken as PB in each of the lines

Fig. 2A presents a heat map and a hierarchical clustering of gene expression across tissues, developmental stages and genotypes. This analysis revealed major differences in PB between WT and the MaMADS2 repressed lines, but almost no differences at MG or breaker/ripe, supporting that the fruits from the repressed lines and WT were in a similar physiological stages at MG and B/R in both peel and pulp (Fig. 2B). Three different expression clusters are defined as I, II and III (Fig. 2A). Cluster I includes all the expression averages of samples taken immediately after harvest (MG) from WT and the two transgenic plants both peel and pulp. Cluster I also contains expression averages of samples taken at PB from the pulp and peel of lines #32 and #44, but not from PB of WT. Hence, this cluster represents the least physiologically mature tissues. Cluster II includes all the expression averages of samples taken from peel at B and at R and the expression averages of samples taken from WT peel at PB. Effectively, it contains all the ripening peel samples. Cluster III includes all the expression averages of samples taken from pulp at B and at R and the expression averages of samples taken from WT pulp at PB and hence it contains all the ripening pulp samples. Hence, the transcriptome profiles of WT during PB of peel and pulp are closer to the transcriptome profiles of these tissues at B, than to that of MG. In contrast, the transcriptome profiles during PB of line #32 and line #44 pulp and peel are closer to the transcriptome profiles of MG, than to that of B.

Fig. 2
figure 2

Transcriptome differences between WT and MaMADS2 repressed banana fruits during ripening. a. Heatmap of the normalized gene expression averages of all samples and hierarchical clustering of all samples. On right the different samples clustered into three groups (I, II and III). The analysis was performed on samples taken from peel (P) and pulp (U) of WT (WT) and repressed lines MaMADS2 antisense (#44) and MaMADS2 RNAi (#32) at harvest/mature green (H), prebreaker (PB), breaker (B) and ripe (R). the color scale shows the distances between different conditions; red show low distance while blue larger distance; b. Number of genes significantly changed (FDR < 0.05) between #32 or #44 lines and WT fruits in peel and pulp in the developmental stages of MG, PB, B and R

Transcriptome analysis revealed elevated expression of hormone-related genes at the pre-breaker stage in MaMADS2 repressed lines

Comparison of expression data of WT and MaMADS2 repressed lines #32 and #44 revealed a significant difference in gene expression levels for the PB stage in pulp and to a lesser degree in the peel, and no significant change between the WT and the transgenic plants in MG, B and R (Fig. 2B). This further supports the fact that all lines were collected at a comparable developmental stage and that transgenic fruits achieved full ripeness, although later.

In the peel of line 32 the expression of 162 genes was significantly changed at PB in comparison to WT (FDR < 0.05), and more than 90% of these genes showed the same direction of change in the peel of line #44 (increase or decrease), with larger FDR. This likely reflects the less severe phenotype of line #44 as compared to line #32. In the pulp of lines #32 and #44 the expression of 692 and 1374 genes, respectively, were significantly changed at PB (Additional file 1: Fig. S1A). Interestingly, although line #32 showed a more delayed phenotype than line #44, fewer genes were significantly changed in its pulp. Perhaps the parameters of fruit ripening reflecting mainly the events which occur in the peel (color change, ethylene production), which were the base for transcriptome analysis sampling, do not fully capture pulp ripening outcomes. Most of the genes that were changed in the pulp of #32 were also significantly changed in the pulp of line #44 and the rest exhibited the same directional change. Similarly, up to 99% of the genes that were altered in the pulp of line #44 showed the same directional change in #32 line. Based on this consistency, the restriction imposed on selection of the genes from line #32 was reduced (FDR < 0.1) and 1216 modified genes from this line were used in subsequent analysis (Additional file 1: Fig. S1B). Hence, the transcriptome analysis was performed on 1685 genes which constitute the union between the modified genes in #32 and #44 pulp at PB.

Out of 1685 modified genes, 545 and 1140 genes showed lower and higher expression, respectively, in both lines #44 and #32 compared to WT at PB. Gene Ontologizer (GO) was applied to identify functional groups (Additional file 1: Fig. S2A, B). Among the genes exhibiting lower expression in MaMADS2 repressed lines, there were groups related to protein assembly and to cytoskeleton and organelle movement, indicating that these processes are induced by MaMADS2 and likely associated with ripening (Additional file 1: Fig. S2A). Among the genes exhibiting higher expression, were groups related to RNA metabolism and gene expression, transport, oxidative stress, various biosynthetic pathways and hormonal signaling (Additional file 1: Fig. S2B), indicating that these processes are negatively influenced by MaMADS2. Further analysis of the hormonal signaling pathways was performed by KEGG on the genes included in the 1685 differentially expressed gene pool. This analysis identified putative proteins involved in the hormonal synthesis and response pathways of ethylene, ABA, IAA, CK, GA, SA and JA.

MaMADS2 influences ethylene synthesis and response pathways genes

Genes of the ethylene biosynthesis pathway have significantly lower expression in MaMADS2 repressed lines than in WT at PB in pulp (Fig. 3). One homolog of SlACS2 (GSMUA_Achr4T29150) involved in tomato ripening ethylene biosynthesis [46], is expressed in the pulp of WT at MG and its expression increased by three fold at PB (Fig. 3A), and a similar expression pattern was observed in the WT peel (Fig. 3A). In contrast, the level of GSMUA_Achr4T29150 in the pulp of MaMADS2 repressed lines, was significantly lower at MG and PB (less than 4% and about 30%, for lines #32 and #44, respectively), and also in the peel of line #32 at PB stage (Fig. 3A, B, Additional file 2: Table S2).

Fig. 3
figure 3

Ethylene biosynthesis and signaling pathway in WT and MaMADS2 repressed fruits. a. Heat map of the genes changed in PB between WT and the MaMADS2 repressed lines. Each gene is normalized to the expression level in WT at MG; b. The ethylene biosynthesis and signaling pathway scheme. The genes’ expression significantly higher (red) or lower (blue) in #32 or #44 fruits in comparison to WT at PB. Yellow color indicates no change. Detailed description of genes is provided at [47]

Fifteen ACO homologs were expressed in banana fruits with expression levels above 50 BaseMean. Two tomato SlACO1 homologs, which contribute to ethylene production during ripening [46], were identified in banana fruits and one, GSMUA_Achr1T18250, is also significantly lower at PB in MaMADS2 repressed lines in comparison to WT (Fig. 3A, B, Additional file 2: Table S2). Another SlACO1 homolog, GSMUA_Achr7T18420 (ACO1) is highly expressed (100 times more than the others ACO genes) and it is higher in WT than in the repressed lines only in the peel (Additional file 2: Table S2). Other ACO genes presented a diversity of expression patterns in both peel and pulp, however, most of the induced genes showed earlier induction in WT than in the repressed lines (Additional file 1: Fig. S3).

KEGG analysis also identified a banana homolog of the negative regulator of the tomato ethylene response pathway SlCTR1 [47]. CTR1 homologue GSMUA_Achr8T07600 is expressed in WT at MG in both pulp and peel but decreases sharply to 3% and 18% at PB, respectively (Fig. 3A, Additional file 2: Table S2). This gene is highly expressed in the pulp and the peel of the MaMADS2 repressed lines at MG, but also at PB and decrease at B (Additional file 2: Table S2). Downstream in the ethylene response pathway, a group of ERF (Ethylene Response Element) homologs were also highly expressed in the pulp of MaMADS2 repressed lines at PB. One of them, GSMUA_Achr8T06870 (ERF11), was previously demonstrate to represses the expression of ACO1 [48].

MaMADS2 influences ABA synthesis and response

ABA level is dependent on the rate of biosynthesis and catabolism. The 9-cis-epoxycarotenoid Dioxygenases (NCED) encodes an enzyme involved in ABA biosynthesis, while CYP707A, encodes a hydroxylase that converts ABA to 8′-OH-ABA and hence is involved in ABA degradation [49, 50]. In the transcriptome data there are two expressed NCED homologs and the expression of one, GSMUA_Achr4T22870, was low in pulp of WT and significantly higher in MaMADS2 repressed lines at PB (Fig. 4A, D, Additional file 2: Table S2). The expression of another NCED gene GSMUA_Achr5T02570 increased during ripening in WT, but oscillated in the MaMADS2 repressed lines with high levels at MG. In addition, one CYP707A was found to be highly expressed in the banana peel (GSMUA_Achr7T07200), but it had low expression in the repressed lines at PB, as well as, in the pulp of WT (Fig. 4A, Additional file 2: Table S2).

Fig. 4
figure 4

ABA biosynthesis and signaling pathway and ABA level in WT and MaMADS2 repressed fruits. a. Heat map of the genes changed in PB between WT and the MaMADS2 repressed lines. Each gene is normalized to the expression level in WT at MG; b. ABA level in WT and #32 pulps; c. The average of MaMADS2 expression levels in WT after ABA treatment (**: p < 0.01 standard error is presented); d. The ABA biosynthesis and signaling pathway scheme. The genes’ expression significantly higher (red) or lower (blue) in #32 or #44 fruits in comparison to WT at PB. Yellow color indicates no change. Detailed description of genes is provided by [50, 51]; e. Effect of ethylene treatment on ABA content in banana slices; f. A suggested model for the interactions of ABA with MaMADS2 and ethylene: rectangular: transcript, oval: protein, X: interaction that was rejected based on our results

In agreement with GSMUA_Achr4T22870 and GSMUA_Achr5T02570 expression, the ABA levels increase in pulp and peel (Fig. 4B and unpublished data) at the transition from MG (average H°: WT 120.3 and line #32 120.9) to PB (average H°: WT 108.9 and line ##32 111) in all lines prior to ethylene peak which occurred at B (Fig. 1) and were higher in the MaMADS2 repressed lines than in WT in pulp at all stages (Fig. 4B).

To elucidate the interactions between ethylene, ABA and MaMADS2-box, WT banana slices were treated with ABA and the effect on ethylene and MaMADS2 expression was examined. Ethylene production peak occurred on the fifth day following ABA treatment, and the level was 5.42 ± 0.75 μl/kg*h, while the peak in non-ABA treated slices occurred on the seventh day and the level was 3.75 ± 0.22 μl/kg*h. ABA also increased the expression of MaMADS2 (Fig. 4C). The reciprocal experiment showed, however, that ethylene treatment did not increase ABA levels (Fig. 4E).

Homologs of the ABA response pathway: ABA receptor, PYRABACTIN RESISTANCE1 (PYR/PYL) (GSMUA_Achr8T19910, GSMUA_AchrUn_randomT25050, and GSMUA_Achr3T03480); Type 2C protein phosphatases (PP2C) which activates the SnRK2 and is inactivated by binding to the receptor (GSMUA_Achr8T05700) and the sucrose non fermenting (SNF1)-related protein kinase (SnRK2) which phosphorylate proteins involved in ABA response (GSMUA_Achr6T25570 and GSMUA_Achr3T01710) [50, 51] were expressed in the pulp at WT MG but declined thereafter (Fig. 4A, Additional file 2: Table S2). A similar pattern was observed in the peel (Fig. 4A). The expression of these genes was significantly higher in both pulp and peel at PB in the repressed lines than in WT and was reduced later on (Fig. 4A).

In order to examine if there is an ABA-related transcriptome footprint in WT and the transgenic lines, the “Arabidopsis homologs gene pool” (see Methods) was analyzed by the Hormonometer [45]. The analysis showed high correlations between the query gene pool and ABA transcriptomic footprint of Arabidopsis thaliana (Fig. 9). This indicates that both transgenic lines have higher expression of ABA-induced genes at BP compared to WT.

MaMADS2 regulates Zeatin level and response

CK levels are determined by the biosynthesis and degradation enzymes Isopentenyl transferases and CK dehydrogenase (CKX), respectively [52]. Two homologs of CK dehydrogenase (CKX), GSMUA_Achr10T28080 and GSMUA_Achr9T17760, were significantly modified through ripening and in transgenic lines (Fig. 5A). GSMUA_Achr10T28080 was hardly expressed in WT in the pulp, but was expressed in the repressed lines. GSMUA_Achr9T17760, was expressed at MG in WT pulp and peel, and in pulp it decreased at PB, but in the repressed lines the decline was delayed (Fig. 5A, Additional file 2: Table S2).

Fig. 5
figure 5

CK signaling pathway and zeatin level in WT and MaMADS2 repressed fruits. a. Heat map of the genes changed in PB between WT and the MaMADS2 repressed lines. Each gene is normalized to the expression level in WT at MG; b. Zeatin level in WT and #32 pulps; c. The CK signaling pathway scheme. The genes’ expression significantly higher (red) or lower (blue) in #32 or #44 fruits in comparison to WT at PB. Yellow color indicates no change. Detailed description of genes is provided at [55]

Zeatin was the only detectable CK in our analysis (Fig. 5B). In WT there was no significant change during ripening, however, the levels of zeatin were lower in MaMADS2 suppressed line than in WT, and gradually increased during ripening in this line (Fig. 5B).

Plants respond to CK by using a phospho-relay mechanism which involves the membrane two component receptor, Cytokinin Response (CRE1) [53] and downstream transcription activators A/B-Arabidopsis Response Regulators (ARR) [54, 55]. Consistent with lower zeatin levels in MaMADS2 repressed lines than in WT, two CRE1 homologs (GSMUA_Achr8T22930 and GSMUA_AchrUn_randomT08400) clearly exhibited lower expression in the repressed lines than in WT pulp at PB (Fig. 5A, C and Additional file 2: Table S2). In addition, the expression levels of all other 7 CRE1 homologs expressed in ripening fruits was lower in the pulp of the repressed lines than in WT (more than twice the levels) at MG and/or PB (Additional file 1: Fig. S4).

On the other hand, the expression of A-ARR homologs (GSMUA_Achr1T17800 and GSMUA_Achr5T08230) was higher at PB in the pulp of the repressed lines in comparison to WT (Fig. 5A, C and Additional file 2: Table S2). It has been shown that the A-ARR are induced by CK and they repress their own expression [56], shutting down the CK response [57]. Moreover, the only ARR-B homologe, GSMUA_Achr8T30850, that was highly expressed in banana fruits was lower (less than 20%) in the pulp of repressed lines compared with WT at PB (Fig. 5A and Additional file 2: Table S2). Hormonometer analysis [45] showed high negative correlation between the query gene pool and a CK response footprint, indicating that the cytokinin-related transcriptomic footprint is lower in the transgenic lines compared to WT.

MaMADS2 has a dual effect on auxin response pathway(s)

There are no differences in IAA levels between WT and MaMADS2 repressed lines, and the levels of auxin was high at MG, declined at PB and remained low during maturation (Fig. 6B). Key elements of auxin transport and signaling pathways were modified in the repressed lines (Fig. 6A, C). Auxin transport is controlled by influx (AUX) and efflux carriers (PIN) and gene transcription is regulated through the TIR-AUX/IAA-ARF pathway (Fig. 6C) [58]. The transcript expression of auxin transport components and TIR-AUX/IAA-ARF signaling pathway GSMUA_Achr11T24200 (AUX1), GSMUA_Achr9T27160 (PIN1), GSMUA_Achr6T20580 (TIR1), GSMUA_Achr9T25550 (AUX/IAA) and GSMUA_Achr4T10590 and GSMUA_Achr6T36980 (ARF) genes were significantly higher in the pulp of repressed lines at PB in comparison to WT (Fig. 6A, C and Additional file 2: Table S2). This is primarily due to a delay in the reduction of these transcripts during ripening in the transgenic lines. All of these genes exhibited high expression levels at MG in WT pulp, and most of them also in peel, and their expression declined at PB to a very low levels, yet the decline was delayed to B in the repressed lines (Fig. 6A and Additional file 2: Table S2). In addition, all 10 homologs of TIR1 in banana follow this profile in both peel and pulp (Additional file 1: Fig. S5).

Fig. 6
figure 6

IAA signaling pathway and IAA level in WT and MaMADS2 repressed fruits. a. Heat map of the genes changed in PB between WT and the MaMADS2 repressed lines. Each gene is normalized to the expression level in WT at MG; b. IAA level in WT and #32 pulps; c. The IAA signaling pathway scheme. The genes’ expression significantly higher (red) or lower (blue) in #32 or #44 fruits in comparison to WT at PB. Yellow color indicates no change. Detailed description of genes is provided at [58, 59]

Recently, another poorly understood pathway involving protein translation has been suggested to contribute to auxin response [59] (Fig. 6C). This pathway is comprised of TARGET OF RAPAMYCIN (TOR) and Phospholipase D zeta (PLDz) which produce a lipid-related second messenger, phosphatidic acid (PA), that activates the PHOSPHOLIPID-DEPENDENT PROTEIN KINASE 1 (PDK1) involved in RIBOSOMAL S6 KINASE (S6K) regulation of protein translation [59]. Transcript homologs of several genes in this pathway, GSMUA_Achr9T16530 (PLDz), GSMUA_Achr10T10490 (PDK1) and GSMUA_Achr4T01370 (TOR), showed elevated expression during ripening and were much lower in the pulp of MaMADS2 repressed lines at PB (Fig. 6A, C and Additional file 2: Table S2). In the peel there were no differences in the expression of these genes between WT and the repressed lines.

MaMADS2 modifies GA synthesis and response

The compounds GA19, GA24, GA9, and GA20 are precursors of the active GA4/GA1 and GA3/GA7 [60]. Analysis of GA levels in the pulp of WT and transgenic plant revealed that, while GA9 level was similar in both WT and the repressed line in the pulp (Additional file 1: Fig. S6B), that of GA24 was higher and that of GA4 was lower at MG in the repressed line in comparison to WT (Fig. 7B). GA1, GA19, and GA20 could not be detected. These results coincide with the expression levels of genes encoding enzymes of biosynthesis and degradation of GA. GA 20-oxidase is responsible for the synthesis of the GA precursors [60] and the expression of its homolog GSMUA_Achr8T19120 was very low in WT pulp at all stages, and slightly higher in the peel reaching a peak at B stage (Fig. 7A and Additional file 2: Table S2). However, in the repressed line in both pulp and peel the expression was higher at MG and PB, possibly explaining the higher GA24 precursor in the repressed line. The expression of GSMUA_Achr10T21600, a homolog of GA 2-oxidase, which deactivates GA [60] decreased sharply after MG in the WT pulp. However, in both repressed lines the expression in pulp remained high through PB (Fig. 7A and Additional file 2: Table S2). This likely explains the low levels of GA4 observed in pulp of MaMADS2 suppress line.

Fig. 7
figure 7

GA biosynthesis and signaling pathway and GA levels in WT and MaMADS2 repressed fruits. a. Heat map of the genes changed in PB between WT and the MaMADS2 repressed lines. Each gene is normalized to the expression level in WT at MG; b. GA24 and GA4 levels in WT and #32 pulps; c. The GA biosynthesis and signaling pathway scheme. The genes’ expression significantly higher (red) or lower (blue) in #32 or 44 fruits in comparison to WT at PB. Yellow color indicates no change. Detailed description of genes is provided at [60, 61]

GA signaling includes GIBBERELLIN-INSENSITIVE DWARF1 (GID1/2) genes encoding gibberellins receptors which mediate the degradation of the negative regulator DELLA [61]. The expression of GSMUA_Achr3T24930 and GSMUA_Achr8T05910 (GID1), GSMUA_Achr9T22610 (GID2) and GSMUA_Achr8T18250 (DELLA) was higher in the MaMADS2 repressed lines pulp at MG and PB in comparison to WT (Fig. 7A, C and Additional file 2: Table S2). Moreover, out of six additional DELLA homologs three were expressed and their levels were higher in the pulp of repressed line at MG and PB than in WT (Additional file 1: Fig. S6A).

Combining high levels of DELLA homologs functioning as negative regulator(s) [61], with the low levels GA4 in the MaMADS2 suppressed lines, indicates that GA response is likely reduced in the transgenic lines. However, a Hormonometer analysis [45] of the query gene pool showed non-conclusive results; while there was negative correlation between banana query gene pool with 3 h GA3-induced response, there was a positive correlation with 6 and 9 h GA3-induced response (Fig. 9).

MaMADS2 influences JA and SA levels as well as their response pathways

The levels of JA and SA were significantly higher in MaMADS2 repressed lines than in WT in most ripening stages (Fig. 8B, E). The transcripts of the JA biosynthesis and JA and SA responses pathways, exhibited higher expression in the repressed lines in the pulp at PB (Fig. 8A, C, D, F and Additional file 2: Table S2). The expressions in WT pulp of GSMUA_Achr4T15630 (JAR1 homolog, involved in formation of an active JA conjugate) [62] and GSMUA_Achr9T14370 and GSMUA_Achr11T09290 (JAZ homologs, negative regulator of the pathway) [62] are mostly low in WT pulp throughout ripening, but elevated in the pulp of the repressed lines (Fig. 8A, C and Additional file 2: Table S2). The Hormonometer analysis showed a positive correlation between the query gene pool and the JA- transcriptome foot print, indicating a relatively high JA-related transcriptome foot print of the transgenic lines (Fig. 9).

Fig. 8
figure 8

JA and SA signaling pathway and JA and SA levels in WT and MaMADS2 repressed fruits. a. Heat map of the genes in JA signaling pathway changed in PB between WT and the MaMADS2 repressed lines. Each gene is normalized to the expression level in WT at MG; b. JA level in WT and #32 pulps; c. The JA signaling pathway scheme. The genes’ expression significantly higher (red) or lower (blue) in #32 or #44 fruits in comparison to WT at PB. Yellow color indicates no change; d. Heat map of the genes in SA signaling pathway changed in PB between WT and the MaMADS2 repressed lines. Each gene is normalized to the expression level in WT at MG; e. SA level in WT and #32 pulps; f. The SA signaling pathway scheme. The genes’ expression significantly higher (red) or lower (blue) in #32 or #44 fruits in comparison to WT at PB. Yellow color indicates no change. Detailed description of genes is provided at [62] for JA and at [63] for SA

Fig. 9
figure 9

Difference in hormonal footprint between MaMADS2 repressed lines and WT fruits. Hormonometer analysis of the “Arabidopsis homologs gene pool” (see methods) expressed in MaMADS2 repressed lines in comparison to WT. Dark brown indicates strong hormone footprint in MaMADS2 repressed line and blue indicates strong hormone footprint in WT.

The SA signaling pathway includes two major genes NON-EXPRESSOR OF PATHOGENESISRELATEDGENES1 (NPR1) and TGA that acts together to activate the SA-induced transcription [63]. Homologs of these genes GSMUA_Achr5T08160 (NPR1) and GSMUA_Achr2T06770 (TGA) had higher expression in PB in the pulp of the MaMADS2 repressed lines (Fig. 8D, F and Additional file 2: Table S2). The expression of these genes in WT pulp decreased immediately after MG and gradually during ripening in the peel. Nevertheless, the decline in expression occurs also in the pulp of the transgenic, but was delayed till B, while that in peel was similar to that in WT (Fig. 8D and Additional file 2: Table S2). The Hormonometer analysis showed a small positive correlation between the query gene pool to the transcriptome footprint of SA [45] (Fig. 9), suggesting a SA-transcriptome footprint in the transgenic lines.

Discussion

Transcriptome analysis of peel and pulp revealed similar pattern at MG which diverged during maturation to distinct peel and pulp transcriptome, indicating that different ripening programs exist in these parts. There were almost no differences between WT and the MaMADS2 repressed lines at MG or breaker/ripe supporting that fruits from the repressed lines and WT were in a similar physiological stages at MG and B/R which indicates that repressed lines reached full ripening (Fig. 2B). The majority of differences between the MaMADS2 repressed fruit and WT appeared in the pulp of PB. In addition, the time from PB to B was extended in the repressed lines, in comparison to WT which indicates that MaMADS2 has a role in enhancing the ripening progress. This further supports our previous conclusion that MaMADS2 determines the rate of ripening [4]. Since GO analysis identified functional groups related to hormonal regulation which were significantly changed in the pulp of MaMADS2 repressed lines in comparison to WT (Fig. 2D), it is suggested that MaMADS2 coordinate the hormone response of ethylene, ABA, auxin, cytokinin, gibberellins, jasmonic acid and salicylic acid prior to the breaker stage.

Ethylene biosynthesis was lower and delayed in the MaMADS2 repressed lines in comparison to WT (Fig. 1), and it was in agreement with the expression of a system II tomato SlACS2 homolog that was suppressed at PB in MaMADS2 repressed lines, but increased later on. It is also in agreement with lower expression in the pulp of the repressed lines of few of the ACO genes (Fig. 3). The expression of both ACS and ACO genes was lower also in tomato SEPALLATA SlMADS-RIN (rin) mutants [64, 65]. In tomato the SlMADS-RIN binds directly to promoters of ACS genes and not ACO genes, but it is not clear if MaMADS2 acts in a similar way. It has been demonstrated that the ERF11 in banana recruits histone deacetylase to the promoters of ACO1 gene and most likely deactivates its expression [48]. In agreement with lower ethylene production in MaMADS2 repressed lines, this gene is up regulated in these lines.

Another gene that was modified in tomato rin mutants is the negative regulator CTR1 [66]. Similar to rin mutant, in MaMADS2 repressed lines the expression of CTR1 homolog was higher in PB compared to WT. However, while the expression in banana is reduced towards breaker, in tomato it is not [67], and it is induced by ethylene [68].

ABA signaling response gene expression and hormone levels revealed a crosstalk between ABA, ethylene and MaMADS2. In the transgenic plants, low levels of ethylene and low MaMADS2 expression [4], coexist with high levels of ABA at all ripening stages (Figs. 1 and 4). To exclude the possibility that ethylene may affect ABA levels and the low levels of ethylene which occur in the transgenic plants might cause the ABA increase, ethylene application to MG banana fruits did not alter the levels of ABA (Fig. 4E). Similarly, treatment of tomato fruits with ACC, a precursor of ethylene, did not alter ABA levels [22], indicating there is no direct effect of ethylene on ABA synthesis (Fig. 4F). However, since application of ABA reduce the time to ethylene induction and increased ethylene production, and also increased the levels of MaMADS2 expression (Fig. 4C), we suggest that ABA acts via an increase in MaMADS2 expression to enhance ethylene biosynthesis and ripening. This mechanism may explain results in tomato and banana showing that ABA increased ethylene production [21, 22, 24].

Higher level of ABA was observed in MaMADS2 repressed fruit, suggests that ABA production in banana fruit is under autoinhibitory regulation via MaMADS2 (Fig. 4F). The expression of NECD encoding genes responsible for ABA synthesis may explain both the increase in ABA during ripening and the higher levels in repressed lines, since the expression of one gene was increased at PB in WT and another was higher in the repressed lines than in WT.

The genes within the ABA signaling pathway most likely are negatively regulated by MaMADS2-box gene, since the reduction in their expression is delayed in the repressed lines. Likewise, in strawberry, a non-climacteric fruit, the SnRK2 transcript level decreased during ripening [69]. Higher levels of ABA and higher expression of the genes within the signaling response pathway in the transgenic lines might explain the high ABA-transcriptomic footprint which exists in the MaMADS2 repressed lines (Fig. 9).

Most surprising are the results of cytokinin signal transduction and levels. Cytokinins action is usually associated with delayed senescence and fruit ripening [18], but our data showed the zeatin levels are lower in the repressed line than in WT at MG (Fig. 5B), suggesting that MaMADS2 enhances cytokinins levels prior to climacteric stage. Interestingly, CKs levels were highest in the red sectors (in comparison to green sectors) of transgenic plants overexpressing ipt gene involved in cytokinin biosynthesis [32]. In agreement with a positive role for cytokinin in ripening, The CRE genes, encoding cytokinin receptors, were higher in WT than in the repressed lines (Fig. 5A and Additional file 1: Fig. S4). In addition, the reduced expression in WT of A-ARR during maturation and the high levels in the transgenic plants further support a positive role of cytokinin in ripening. The A-ARR is induced by CK and overexpression of type A-ARRs represses their own expression [56], leading to shutting down of the CK response [57]. This mechanism might be augmented in the repressed lines, and as can be expected in this case, cytokinin-related transcriptomic footprint is low in those fruits (Fig. 9). It has been demonstrated in Arabidopsis that a components of the CK response pathway (ARR2) play a role in signaling of both CK and ethylene [57]. Hence, it is possible that higher CK action in WT leads to higher ethylene production by induction or shutdown of proteins leading to stabilization or degradation, respectively, of ACS homologous enzymes [70, 71]. Taken together, our results strongly suggest a positive role for CK in ripening.

Although accumulating evidence suggest a role for auxin in fruit ripening [16], the exact role is still not clear. Auxin level was high in MG pulp and declined at PB. This result is in agreement with reports in tomato and grapes where an auxin peak was observed prior to ripening and decreased during later maturation [18, 21, 72]. In contrast to banana MaMADS2 repressed lines, in the homologous tomato rin and apple MdMADS8/9 repressed lines, there were higher auxin levels in the than in WT [73, 74]. These differences might results from the sampling of more mature fruit in our study or due to species specific differences.

Repressing the expression of the MaMADS2 affected differently the two auxin signaling pathways; on one hand, it increased the expression of genes within the TIR-AUX/IAA-ARF pathway and on the other hand, it reduced the expression of genes in IAA-induced protein translation pathway. The lower expression of genes within the latest pathway in the repressed banana pulp might impair the translation of upstream open reading frames (uORF) containing mRNAs [75]. Interestingly, it has been suggested that the translation of few proteins in the ethylene response pathway are regulated by uORFs [76]. In summary, our results suggest that in the MaMADS2 repressed lines, auxin activity lingers on at PB, although levels decrease at this stage. This is consistent with the hypothesis developed based on data in grape and tomato that higher auxin levels precede ethylene induction and only when auxin levels and response is reduced, it enables ripening to proceed [8, 13, 14].

GA and SA were shown to delay ripening and JA to promote it [8, 18, 27, 33, 77]. Our data provide molecular support to a role for MaMADS2 in mediating the ripening-related alteration in the levels of these hormone, however, their role in banana fruit ripening remains unclear.

Conclusions

The majority of differences in WT versus MaMADS2 repressed fruit reflected slower acquisition of ripening phenotypes that together is manifested in delayed ripening. That is, transcriptome changes in the repressed lines generally follow WT but at additional days after harvest and to an overall lesser magnitude. The hormonal response-related transcriptome and hormonal profiles of banana fruits from WT and MaMADS2 repressed lines point to a major effect of MaMADS2 on multiple hormonal signaling pathways and suggest that MaMADS2 coordinate the response of the different hormones, besides ethylene, in the regulation of fruit ripening. Moreover, this study suggests an intricate crosstalk between ABA and MaMADS2 and a positive role for cytokinin in fruit ripening. Since the regulation of hormonal response during ripening was rarely investigated in fleshy fruit, this study open further investigations.

Abbreviations

ABA:

Abscisic acid

B:

Breaker

CK:

Cytokinin|

GA:

Gibberellin

IAA:

Indole acetic acid

JA:

Jasmonic acid

MG:

Mature Green

PB:

Prebreaker

SA:

Salicylic acid

References

  1. Giovannoni J, Nguyen C, Ampofo B, Zhong S, Fei Z. The epigenome and transcriptional dynamics of fruit ripening. Ann RevPlant Biol. 2017;68:61–84.

    Article  CAS  Google Scholar 

  2. Klee HJ, Giovannoni JJ. Genetics and control of tomato fruit ripening and quality attributes. Ann RevGen. 2011;45:41–59.

    Article  CAS  Google Scholar 

  3. Seymour GB, Østergaard L, Chapman NH, Knapp S, Martin C. Fruit Development and Ripening. Ann RevPlant Biol. 2013;64:219–41.

    Article  CAS  Google Scholar 

  4. Elitzur T, Yakir E, Quansah L, Zhangjun F, Vrebalov J, Khayat E, Giovannoni JJ, Friedman H. Banana MaMADS transcription factors are necessary for fruit ripening and molecular tools to promote shelf-life and food security. Plant Physiol. 2016;171:380–91.

    Article  CAS  Google Scholar 

  5. Elitzur T, Vrebalov J, Giovannoni JJ, Goldschmidt EE, Friedman H. The regulation of MADS-box genes expression during ripening of banana and their regulatory interaction with ethylene. J Exp Bot. 2010;61:1523–35.

    Article  CAS  Google Scholar 

  6. Seymour GB. Banana. In: Seymour GB, Taylor J, Tucker G, editors. Biochemistry of fruit ripening. London: Chapman and Hall; 1993. p. 83–106.

    Chapter  Google Scholar 

  7. Grierson D. Ethylene and the control of fruit ripening. In: Seymour GB, Poole M, Giovannoni JJ, Tucker GA, editors. The molecular biology and biochemistry of fruit ripening: Blackwell Publishing Ltd.; 2013. p. 43–73.

  8. Kumar R, Khurana A, Sharma AK. Role of plant hormones and their interplay in development and ripening of fleshy fruits. JExp Bot. 2014;65:4561–75.

    CAS  Google Scholar 

  9. Cherian S, Figueroa CR, Nair H. ‘Movers and shakers’ in the regulation of fruit ripening: a cross-dissection of climacteric versus non-climacteric fruit. J ExpBot. 2014;65:4705–22.

    CAS  Google Scholar 

  10. Osorio S, Scossa F, Fernie AR. Molecular regulation of fruit ripening. FrontPlant Scie. 2013;4:198.

    Google Scholar 

  11. Cohen JD. In vitro tomato fruit cultures demonstrate a role for indole-3-acetic acid in regulating fruit ripening. J Amer Soci Hort Scie. 1996;121:520–4.

    CAS  Google Scholar 

  12. Vendrell M. Reversion of senescence: effects of 2, 4-dichlorophenoxyacetic acid and indoleacetic acid on respiration, ethylene production, and ripening of banana fruit slices. AusJBiolScie. 1969;22:601–10.

    CAS  Google Scholar 

  13. Tatsuki M, Nakajima N, Fujii H, Shimada T, Nakano M, Hayashi K-i, Hayama H, Yoshioka H, Nakamura Y. Increased levels of IAA are required for system 2 ethylene synthesis causing fruit softening in peach (Prunus persica L. Batsch). J Exp Bot. 2013;64:1049–59.

    Article  CAS  Google Scholar 

  14. Trainotti L, Tadiello A, Casadoro G. The involvement of auxin in the ripening of climacteric fruits comes of age: the hormone plays a role of its own and has an intense interplay with ethylene in ripening peaches. J Exp Bot. 2007;58:3299–308.

    Article  CAS  Google Scholar 

  15. Breitel DA, Chappell-Maor L, Meir S, Panizel I, Puig CP, Hao Y, Yifhar T, Yasuor H, Zouine M, Bouzayen M. AUXIN RESPONSE FACTOR 2 intersects hormonal signals in the regulation of tomato fruit ripening. PLoS Gen. 2016;12:e1005903.

    Article  Google Scholar 

  16. Hao Y, Hu G, Breitel D, Liu M, Mila I, Frasse P, Fu Y, Aharoni A, Bouzayen M, Zouine M. Auxin response factor SlARF2 is an essential component of the regulatory mechanism controlling fruit ripening in tomato. PLoS Gen. 2015;11:e1005649.

    Article  Google Scholar 

  17. Bouzayen M, Latché A, Nath P, Pech JC. Mechanism of fruit ripening. In: Pua EC, Davey MR, editors. Plant developmental biology - biotechnological perspectives: Springer Berlin Heidelberg; 2010. p. 319–39.

    Google Scholar 

  18. Srivastava A, Handa A. Hormonal regulation of tomato fruit development: a molecular perspective. J Plant Gro Regul. 2005;24:67–82.

    Article  CAS  Google Scholar 

  19. Rossetto MRM, Purgatto E, JROd N, Lajolo FM, Cordenunsi BR. Effects of gibberellic acid on sucrose accumulation and sucrose biosynthesizing enzymes activity during banana ripening. Plant Gro Regul. 2003;41:207–14.

    Article  CAS  Google Scholar 

  20. Leng P, Yuan B, Guo Y. The role of abscisic acid in fruit ripening and responses to abiotic stress. J Exp Bot. 2014;65(16):4577–88.

    Article  CAS  Google Scholar 

  21. Buta J, Spaulding D. Changes in indole-3-acetic acid and abscisic acid levels during tomato (Lycopersicon esculentum mill.) fruit development and ripening. J Plant Gro Regul. 1994;13:163–6.

    Article  CAS  Google Scholar 

  22. Zhang M, Yuan B, Leng P. The role of ABA in triggering ethylene biosynthesis and ripening of tomato fruit. JExpBot. 2009;60:1579–88.

    CAS  Google Scholar 

  23. Mou W, Li D, Luo Z, Mao L, Ying T. Transcriptomic analysis reveals possible influences of ABA on secondary metabolism of pigments, flavonoids and antioxidants in tomato fruit during ripening. PLoS One. 2015;10:e0129598.

    Article  Google Scholar 

  24. Jiang Y, Joyce DC, Macnish AJ. Effect of abscisic acid on banana fruit ripening in relation to the role of ethylene. J Plant Gro Regul. 2000;19:106–11.

    Article  CAS  Google Scholar 

  25. Sun L, Sun Y, Zhang M, Wang L, Ren J, Cui M, Wang Y, Ji K, Li P, Li Q. Suppression of 9-cis-epoxycarotenoid dioxygenase, which encodes a key enzyme in abscisic acid biosynthesis, alters fruit texture in transgenic tomato. Plant Physiol. 2012;158:283–98.

    Article  Google Scholar 

  26. Fan X, Mattheis JP, Fellman JK. A role for jasmonates in climacteric fruit ripening. Planta. 1998;204:444–9.

    Article  CAS  Google Scholar 

  27. Peña-Cortés H, Barrios P, Dorta F, Polanco V, Sánchez C, Sánchez E, Ramírez I. Involvement of jasmonic acid and derivatives in plant response to pathogen and insects and in fruit ripening. J Plant Gro Regul. 2004;23:246–60.

    Google Scholar 

  28. Srivastava MK, Dwivedi UN. Delayed ripening of banana fruit by salicylic acid. Plant Sci. 2000;158:87–96.

    Article  CAS  Google Scholar 

  29. Asghari M, Aghdam MS. Impact of salicylic acid on post-harvest physiology of horticultural crops. TrenFood ScieTechnol. 2010;21:502–9.

    CAS  Google Scholar 

  30. Gan S, Amasino RM. Making sense of senescence. Molecular genetic regulation and manipulation of leaf senescence. Plant Physiol. 1997;113:313–9.

    Article  CAS  Google Scholar 

  31. Vendrell M, Palomer X. Hormonal control of fruit ripening in climacteric fruits. In: VIII international symposium on plant bioregulation in fruit production 463; 1997. p. 325–34.

    Google Scholar 

  32. Martineau B, Houck CM, Sheehy RE, Hiatt WR. Fruit-specific expression of the a. tumefaciens isopentenyl transferase gene in tomato: effects on fruit ripening and defense-related gene expression in leaves. Plant J. 1994;5(1):11–9.

    Article  CAS  Google Scholar 

  33. McAtee P, Karim S, Schaffer R, David K. A dynamic interplay between phytohormones is required for fruit development, maturation, and ripening. FrontPlant Scie. 2013;4:79.

    Google Scholar 

  34. Böttcher C, Burbidge CA, Boss PK, Davies C. Interactions between ethylene and auxin are crucial to the control of grape (Vitis vinifera L.) berry ripening. BMC Plant Biol. 2013;13:222.

    Article  Google Scholar 

  35. Dadzie BK, Orchard JE: Routine postharvest screening of banana/plantain hybrids: criteria and methods. In. Montpellier, France: INIBAP Technical Guidelines; 1997: pp75.

  36. Or E, Belausov E, Popilevsky I, Bental Y. Changes in endogenous ABA level in relation to the dormancy cycle in grapevines grown in a hot climate. J Hort Scie Biotechnol. 2000;75:190–4.

    Article  CAS  Google Scholar 

  37. Zhong S, Joung J-G, Zheng Y, Chen Y-r, Liu B, Shao Y, Xiang JZ, Fei Z, Giovannoni JJ. High-throughput Illumina Strand-specific RNA sequencing library preparation. Cold Spring Harb Protoc. 2011;2011(8) pdb.prot5652.

    Article  Google Scholar 

  38. Gentleman RC, Carey VJ, Bates DM, Bolstad B, Dettling M, Dudoit S, Ellis B, Gautier L, Ge Y, Gentry J. Bioconductor: open software development for computational biology and bioinformatics. Genome Biol. 2004;5:R80.

    Article  Google Scholar 

  39. Anders S. Analysing RNA-Seq data with the DESeq package. Mol Biol. 2010;43:1–17.

    Google Scholar 

  40. D’Hont A, Denoeud F, Aury J-M, Baurens F-C, Carreel F, Garsmeur O, Noel B, Bocs S, Droc G, Rouard M. The banana (Musa acuminata) genome and the evolution of monocotyledonous plants. Nature. 2012;488(7410):213–7.

    Article  Google Scholar 

  41. Pruitt KD, Tatusova T, Maglott DR. NCBI reference sequences (RefSeq): a curated non-redundant sequence database of genomes, transcripts and proteins. Nucleic Acids Res. 2007;35:D61–5.

    Article  CAS  Google Scholar 

  42. Conesa A, Götz S, García-Gómez JM, Terol J, Talón M, Robles M. Blast2GO: a universal tool for annotation, visualization and analysis in functional genomics research. Bioinformatics. 2005;21:3674–6.

    Article  CAS  Google Scholar 

  43. Bauer S, Grossmann S, Vingron M, Robinson PN. Ontologizer 2.0—a multifunctional tool for GO term enrichment analysis and data exploration. Bioinformatics. 2008;24:1650–1.

    Article  CAS  Google Scholar 

  44. Müller M, Munné-Bosch S. Rapid and sensitive hormonal profiling of complex plant samples by liquid chromatography coupled to electrospray ionization tandem mass spectrometry. Plant Methods. 2011;7:37.

    Article  Google Scholar 

  45. Volodarsky D, Leviatan N, Otcheretianski A, Fluhr R. HORMONOMETER: a tool for discerning transcript signatures of hormone action in the Arabidopsis transcriptome. Plant Physiol. 2009;150:1796–805.

    Article  CAS  Google Scholar 

  46. Barry CS, Llop Tous MI, Grierson D. The regulation of 1-aminocyclopropane-1-carboxylic acid synthase gene expression during the transition from system-1 to system-2 ethylene synthesis in tomato. Plant Physiol. 2000;123:979–86.

    Article  CAS  Google Scholar 

  47. Alexander L, Grierson D. Ethylene biosynthesis and action in tomato:a model for climacteric fruit ripening. J Exp Bot. 2002;53:2039–55.

    Article  CAS  Google Scholar 

  48. Han Y, Kuang J, Chen J, Liu X, Xiao Y, Fu C, Wang J, Wu K, Lu W. Banana transcription factor MaERF11 recruits histone deacetylase MaHDA1 and represses the expression of MaACO1 and Expansins during fruit ripening. Plant Physiol. 2016; http://www.plantphysiol.org/content/early/2016/04/05/pp.16.00301.

  49. Nambara E, Marion-Poll A. Abscisic acid biosynthesis and catabolism. Annual Review of Plant Biol. 2005;56:165–85.

    Article  CAS  Google Scholar 

  50. Finkelstein R. Abscisic acid synthesis and response. The Arabidopsis Book. 2013;11:e0166.

    Article  Google Scholar 

  51. Raghavendra AS, Gonugunta VK, Christmann A, Grill E. ABA perception and signalling. Tren Plant Scie. 2010;15:395–401.

    Article  CAS  Google Scholar 

  52. Sakakibara H. CYTOKININS: activity, biosynthesis, and translocation. Ann Rev Plant Biol. 2006;57:431–49.

    Article  CAS  Google Scholar 

  53. Inoue T, Higuchi M, Hashimoto Y, Seki M, Kobayashi M, Kato T, Tabata S, Shinozaki K, Kakimoto T. Identification of CRE1 as a cytokinin receptor from Arabidopsis. Nature. 2001;409(6823):1060–3.

    Article  CAS  Google Scholar 

  54. Perilli S, Moubayidin L, Sabatini S. The molecular basis of cytokinin function. Curr Opin Plant Biol. 2010;13:21–6.

    Article  CAS  Google Scholar 

  55. Werner T, Schmülling T. Cytokinin action in plant development. Curr Opin Plant Biol. 2009;12:527–38.

    Article  CAS  Google Scholar 

  56. Haberer G, Kieber JJ. Cytokinins. New insights into a classic phytohormone. Plant Physiol. 2002;128:354–62.

    Article  CAS  Google Scholar 

  57. Hass C, Lohrmann J, Albrecht V, Sweere U, Hummel F, Yoo SD, Hwang I, Zhu T, Schäfer E, Kudla J. The response regulator 2 mediates ethylene signalling and hormone signal integration in Arabidopsis. EMBO J. 2004;23:3290–302.

    Article  CAS  Google Scholar 

  58. Mockaitis K, Estelle M. Auxin receptors and plant development: a new signaling paradigm. Ann Rev Cell Develop Biol. 2008;24:55–80.

    Article  CAS  Google Scholar 

  59. Bögre L, Henriques R, Magyar Z. TOR tour to auxin. EMBO J. 2013;32:1069–71.

    Article  Google Scholar 

  60. Hedden P, Phillips AL, Rojas MC, Carrera E, Tudzynski B. Gibberellin biosynthesis in plants and fungi: a case of convergent evolution? J Plant Grow Regul. 2001;20:319–31.

    Article  CAS  Google Scholar 

  61. Lor VS, Olszewski NE. GA signalling and cross-talk with other signalling pathways. Essays In Biochem. 2015;58:49–60.

    Article  Google Scholar 

  62. Wasternack C, Hause B. Jasmonates: biosynthesis, perception, signal transduction and action in plant stress response, growth and development. An update to the 2007 review in annals of botany. Ann Bot. 2013;111:1021–58.

    Article  CAS  Google Scholar 

  63. Vlot AC, Dempsey DMA, Klessig DF. Salicylic acid, a multifaceted hormone to combat disease. Annu Rev Phytopathol. 2009;47:177–206.

    Article  CAS  Google Scholar 

  64. Kitagawa M, Ito H, Shiina T, Nakamura N, Inakuma T, Kasumi T, Ishiguro Y, Yabe K, Ito Y. Characterization of tomato fruit ripening and analysis of gene expression in F1 hybrids of the ripening inhibitor (rin) mutant. Physiol Plant. 2005;123:331–8.

    Article  CAS  Google Scholar 

  65. Martel C, Vrebalov J, Tafelmeyer P, Giovannoni JJ. The tomato MADS-box transcription factor RIPENING INHIBITOR interacts with promoters involved in numerous ripening processes in a COLORLESS NONRIPENING-dependent manner. Plant Physiol. 2011;157:1568–79.

    Article  CAS  Google Scholar 

  66. Lin Z, Alexander L, Hackett R, Grierson D. LeCTR2, a CTR1-like protein kinase from tomato, plays a role in ethylene signalling, development and defence. Plant J. 2008;54:1083–93.

    Article  CAS  Google Scholar 

  67. Adams-Phillips L, Barry C, Kannan P, Leclercq J, Bouzayen M, Giovannoni J. Evidence that CTR1-mediated ethylene signal transduction in tomato is encoded by a multigene family whose members display distinct regulatory features. Plant Mol Biol. 2004;54:387–404.

    Article  CAS  Google Scholar 

  68. Leclercq J, Adams-Phillips L, Zegzouti H, Jones B, Latche A. LECTR1, a tomato CTR1-like gene, demonstrates ethylene signaling ability in Arabidopsis and novel expression patterns in tomato. Plant Physiol. 2002;130:1132–42.

    Article  CAS  Google Scholar 

  69. Han Y, Dang R, Li J, Jiang J, Zhang N, Jia M, Wei L, Li Z, Li B, Jia W. SUCROSE NONFERMENTING1-RELATED PROTEIN KINASE2. 6, an ortholog of OPEN STOMATA1, is a negative regulator of strawberry fruit development and ripening. Plant Physiol. 2015;167(3):915–30.

    Article  CAS  Google Scholar 

  70. Wang KL-C, Yoshida H, Lurin C, Ecker JR. Regulation of ethylene gas biosynthesis by the Arabidopsis ETO1 protein. Nature. 2004;428(6986):945–50.

    Article  CAS  Google Scholar 

  71. Chae HS, Faure F, Kieber JJ. The eto1, eto2, and eto3 mutations and cytokinin treatment increase ethylene biosynthesis in Arabidopsis by increasing the stability of ACS protein. Plant Cell. 2003;15:545–59.

    Article  CAS  Google Scholar 

  72. Böttcher C, Keyzers RA, Boss PK, Davies C. Sequestration of auxin by the indole-3-acetic acid-amido synthetase GH3-1 in grape berry (Vitis vinifera L.) and the proposed role of auxin conjugation during ripening. J Exp Bot. 2010;61:3615–25.

    Article  Google Scholar 

  73. Hong S, Lee S. Changes in endogenous plant hormones during ripening of tomato fruits. Physiological Basis of Postharvest Technologies. 1992;343:220–4.

    Google Scholar 

  74. Schaffer RJ, Ireland HS, Ross JJ, Ling TJ, David KM. SEPALLATA1/2-suppressed mature apples have low ethylene, high auxin and reduced transcription of ripening-related genes. AoB Plants. 2013;5.

  75. Schepetilnikov M, Dimitrova M, Mancera-Martínez E, Geldreich A, Keller M, Ryabova LA. TOR and S6K1 promote translation reinitiation of uORF-containing mRNAs via phosphorylation of eIF3h. EMBO J. 2013;32:1087–102.

    Article  CAS  Google Scholar 

  76. Robles L, Stepanova A, Alonso J. Molecular mechanisms of ethylene–auxin interaction. Mol Plant. 2013;6:1734–7.

    Article  CAS  Google Scholar 

  77. Kondo S, Yamada H, Setha S. Effect of jasmonates differed at fruit ripening stages on 1-aminocyclopropane-1-carboxylate (ACC) synthase and ACC oxidase gene expression in pears. J Amer Socie Horti Scie. 2007;132:120–5.

    CAS  Google Scholar 

Download references

Acknowledgements

This is contribution no. 775/17 from Agriculture Research Center (ARO). We thank Ron Izchaki at Zemach Experiments for banana cultivation.

Funding

This work was supported by Israel-USA Bi-National Agricultural Research and Development (BARD fund No. IS-4371-10C and National Science Foundation grants IOS-1322714 (BSF NSF # 2012837) and IOS-1539831. The funding bodies were not involved in data interpretation nor were they involved in any way in writing the manuscript. Aspects of the work were presented as proposed activities in the proposals that were eventually funded.

Availability of data and materials

The datasets generated and/or analyzed during the current study are not publicly available due need for additional data mining, but are available from the corresponding author on reasonable request. Nevertheless, some of the data is presented on Additional file 2: Table S2.

Author information

Authors and Affiliations

Authors

Contributions

EY-experimental design, data mining, and preparation of the manuscript, FZ-RNA seq analysis, NS-bioinformatic analysis, VS-hormonometer analysis, AD and JRJ-RNA preparation, MM and SMB-hormonal analysis, YX and JV- RNA sample generation and RNA-seq data generation, JJG-discussion and critical analysis of results, HF-principal investigator. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Haya Friedman.

Ethics declarations

Ethics approval

Transgenic plants were cultivated according to Israel Plant Protection and Inspection Services (PPIS) regulations and instructions.

Consent for publication

All participants gave their consent for publication.

Competing interests

The authors declare that they have no competing interests.

Publisher’s Note

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

Additional files

Additional file 1:

Figure S1. Number of genes significantly change between #32 or #44 and WT pulp at PB. Fig. S2. Gene Ontology (GO) groups of significantly lower (A) or higher (B) at PB in MaMADS2 repressed fruits than in WT. Fig. S3. Expression of ACO homologs in banana fruits during ripening. Fig. S4. Expression of CRE1 homologs in banana fruits during ripening. Fig. S5. Expression of TIR homologs in banana fruits during ripening. Fig. S6. Expression of DELLA homologs (A) and levels of GA9 (B) in banana fruits during ripening. Table S1. Banana genes annotation and chromosome location. (PPTX 833 kb)

Additional file 2:

Table S2. BaseMean levels of the hormone-related genes in banana fruits during ripening. U-pulp, P-peel, H-harvest (MG), PB-prebreaker, B-breaker, R-ripe. (XLSX 21 kb)

Rights and permissions

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

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Yakir, E., Zhangjun, F., Sela, N. et al. MaMADS2 repression in banana fruits modifies hormone synthesis and signalling pathways prior to climacteric stage. BMC Plant Biol 18, 267 (2018). https://0-doi-org.brum.beds.ac.uk/10.1186/s12870-018-1480-5

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://0-doi-org.brum.beds.ac.uk/10.1186/s12870-018-1480-5

Keywords