Next Article in Journal
Metformin Reduces Potassium Currents and Prolongs Repolarization in Non-Diabetic Heart
Next Article in Special Issue
Unique and Shared Proteome Responses of Rice Plants (Oryza sativa) to Individual Abiotic Stresses
Previous Article in Journal
Podocyte-Related Mechanisms Underlying Survival Benefit of Long-Term Angiotensin Receptor Blocker
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Metabolic Pathways for Observed Impacts of Crop Load on Floral Induction in Apple

1
Agriculture Victoria, AgriBio, Centre for AgriBioscience, Bundoora, VIC 3083, Australia
2
Manjimup Horticulture Research Centre, Department of Primary Industries and Regional Development, Manjimup, WA 6258, Australia
3
Tatura SmartFarm, Agriculture Victoria, Tatura, VIC 3616, Australia
4
Centre for Agricultural Innovation, University of Melbourne, Parkville, VIC 3010, Australia
5
School of Applied Systems Biology, La Trobe University, Bundoora, VIC 3083, Australia
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2022, 23(11), 6019; https://doi.org/10.3390/ijms23116019
Submission received: 7 April 2022 / Revised: 20 May 2022 / Accepted: 25 May 2022 / Published: 27 May 2022
(This article belongs to the Special Issue State-of-the-Art Molecular Plant Sciences in Australia)

Abstract

:
The triggers of biennial bearing are thought to coincide with embryonic development in apple and occurs within the first 70 days after full bloom (DAFB). Strong evidence suggests hormonal signals are perceived by vegetative apple spur buds to induce flowering. The hormonal response is typically referred to as the floral induction (FI) phase in bud meristem development. To determine the metabolic pathways activated in FI, young trees of the biennial bearing cultivar ‘Nicoter’ and the less susceptible cultivar ‘Rosy Glow’ were forced into an alternate cropping cycle over five years and an inverse relationship of crop load and return bloom was established. Buds were collected over a four-week duration within 70 DAFB from trees that had maintained a four-year biennial bearing cycle. Metabolomics profiling was undertaken to determine the differentially expressed pathways and key signalling molecules associated with biennial bearing. Marked metabolic differences were observed in trees with high and low crop load treatments. Significant effects were detected in members of the phenylpropanoid pathway comprising hydroxycinnamates, salicylates, salicylic acid biosynthetic pathway intermediates and flavanols. This study identifies plant hormones associated with FI in apples using functional metabolomics analysis.

Graphical Abstract

1. Introduction

Apple is globally one of the most valued fruit crops and is of high economic importance to the horticultural industry [1]. Apple production, however, is significantly impacted by biennial bearing, a phenomenon that has been studied and reported since the early twentieth century [2]. The cause of alternate bearing is still largely unknown [3] and typically occurs in a tree or branch that does not yield a consistent crop load every year but alternates between a heavy and extremely light one. This common irregularity occurs in both deciduous and evergreen trees, and is reported in nuts, temperate fruits, tropical fruits and forest trees [2]. Fruit growers require a consistent crop load with good quality fruit that attracts a high market value. Several different horticultural practices are implemented to achieve a viable industry. These include crop management strategies such as application of nutrients, pruning and seasonal chemical or hand thinning. Yet, even with the best practices, trees can enter a biennial bearing cycle. A high crop load tree (“ON”) can result in poor quality fruit. Fruit size, color appearance, firmness, sugars and acidity can correlate negatively with high crop load levels [4]. These trees would be required to be adequately thinned to improve fruit quality for the current season but also yield a consistent crop load for the following season. Light crop load trees (“OFF”) are associated with bigger fruit sizes that can have storage disorders [5]. Both phases of the cycle cause a disruption in cropping levels resulting in serious economic losses for the apple industry and can cause significant financial strain on fruit growers.
Histological analysis has shown that the initial development of a flower bud coincides with a developing embryo or fruitlet [6]. The developing seed in young fruitlets is believed to repress flower induction in vegetative buds via phytohormones and removal of fruitlets at the 3–18 mm fruit size, through thinning practices have been shown to lead to a more consistent flower formation (return bloom) for the following season [3,7,8]. There have been investigations of other factors that influence flower formation in apple such as stress associated with temperature, photoperiod [9,10], water deficit as well as internal factors associated with carbon-nitrogen ratio, hormones and interaction with other organs (leaves, terminal shoot growth, and fruit) [3]. A recent investigation showed that return bloom can also be influenced by cultivar and rootstock, that differentially expressed sugars and hormones, suggesting that some cultivars and rootstocks are more susceptible to biennial bearing [8].
Flowering in apple involves four stages: floral induction (FI), flower bud initiation, differentiation, and floral bloom [3,6,11]. FI occurs when a chemical and/or environmental signal causes a vegetative bud to become floral. FI corresponds to a specific time when seeds from developing fruit are believed to emit signals to adjacent bud meristem; the time point at which this occurs can differ in apple cultivars [12]. Floral initiation is the period of histological changes that occur in the flower bud and the subsequent flower differentiation stage describes the visible morphological changes occurring in buds [6,11]. In trees with high crop load, flower bud development is hypothesised to be inhibited by chemical signals, leading to a biennial bearing cycle that is difficult to alter once the cycle is established [13,14].
Regulation of flowering was initially thought to be driven by nutritional competition between flower bud initiation and concurrent fruit formation [11]. An in vitro investigation showed that flower bud development in Plumbago indica was stimulated by disaccharides (sucrose, cellobiose or maltose) and certain plant hormone mixtures (cytokinins, adenine and low levels of auxin), and inhibited by amino acids (glutamine, asparagine), riboflavin and gibberellins [9,15]. More recently, genes [3], transcripts [6,16] and proteins [6] associated with biennial bearing were identified and studies indicated that a combination of increased carbohydrates and activation of FI genes, likely mediated by phytohormones, either stimulated or inhibited FI [6]. Hormone-related genes were likely candidates involved in biennial bearing compared to flowering genes. Moreover, transcriptomics studies showed that hormone responses were also differentially expressed. Plant hormone response to stress is linked to redox hub activity, such that ascorbic acid is either required for biosynthesis of plant hormones or low-levels promote accumulation [16]. A multi-omics investigation also indicated that thiamine, chlorogenic acid and an adenine derivative are involved in the metabolic pathway promoting early flower bud development in apple [6]. Increased levels of flavonoids such as kaempferol derivatives were also identified in low crop load trees.
To determine the metabolic pathways involved in FI, an untargeted metabolomics analysis of buds was utilised for this study. Buds were collected from young apple trees of the biennial bearing cultivar ‘Nicoter’ (marketed as Kanzi®) and the less susceptible to biennial bearing cultivar ‘Rosy Glow’ (marketed as Pink Lady®). The collection of buds was undertaken over eight weeks in the key time frame associated with floral induction as identified by Milyaev et al. [6], who found floral initiation to occur 75 and 97 days after flower bloom (DAFB) for ‘Fuji’ and ‘Gala’ with FI hypothesised to occur at least two weeks earlier. Polar extracts of the buds were analysed using high-resolution mass spectrometry and metabolite expression was compared across bud samples with various crop load treatments using multivariate and univariate analysis to identify key chemical determinants and associated disrupted metabolic pathways involved in inhibition or promotion of FI.
The main objective of the present study was to contribute towards the largely unknown physiological mechanisms of biennial bearing in apple, thereby to better understand the underlying pathways and triggers of FI that might facilitate intervention opportunities for controlling apple crop load and thus ensuring stable apple production.

2. Results

Strong inverse relationships between crop load in the 2018/2019 season and the number of flower clusters in the 2019/2020 season were found for variable crop load trees of both ‘Nicoter’ and ‘Rosy Glow’, as shown in Figure 1. For the constant crop load treatments, ‘Rosy Glow’ showed a similarly strong correlation to the variable treatment trees whereas ‘Nicoter’ showed a poor correlation. This was thought to be due to poor yield despite high crop load treatments (150% and 200%) and replacement trees due to disease at the beginning of the 2017/2018 and 2018/2019 seasons. For this reason, this study only focused on the variable crop load trees.
To investigate the physiological mechanisms of biennial bearing in apples, variable crop load treatments were applied and apple buds collected once a week over 8 weeks post-treatment were analysed using ultra-high performance liquid chromatography–high resolution mass spectrometry (UHPLC-HRMS).
Evaluation of the metabolomic data of ‘Nicoter’ and ‘Rosy Glow’ revealed a total of 908 compounds in positive ionisation mode and 668 compounds in negative ionization mode. Putative (Level 3) identification of 436 metabolites in the positive mode and 265 in the negative mode, provided valuable information on the complex composition of the apple bud, including plant hormones, lipids, amino acids, vitamins and phenols.
Prior to statistical analysis of UHPLC-HRMS results, principal component analysis (PCA) plots were generated to assess the reproducibility of the pooled biological quality control (PBQC) samples for intensity drifts along the batch and no corrections were required. PCA plots of ‘Rosy Glow’ and ‘Nicoter’ (Figure 2) apple bud extracts revealed separation of the two cultivars in the positive (Figure 2A) and negative (Figure 2B) ionisation modes, indicating that the metabolomes of the individual cultivars are distinct. Although, crop load regulation is thought to be similar in cultivar [6], the PCA model indicates that the metabolome between the two cultivars may have differences and were thus explored individually. Furthermore, there may be differences in stress-induced pathways associated with apple trees treated with a variable crop load.
Initially, RV and NV datasets revealed no clear separation between treatment groups for RV (RVHIGH, RVMID, RVLOW) and NV (NVHIGH, NVMID, NVLOW) (Figure S1); however, some separation was observed between high and low treatment groups for each cultivar, particularly in the strong biennial cultivar ‘Nicoter’.
Thus, an orthogonal projection to latent structure discriminant analysis (OPLS-DA) model was applied to ‘Nicoter’ and ‘Rosy Glow’ high and low treatment groups for positive and negative whole datasets. The OPLS-DA score plot revealed discrimination between treatment groups in the positive and negative mode datasets.
The positive mode evaluation of model performance was well described with good predictive performance for NVLOW and NVHIGH (Q2 = 0.601, R2Y = 0.724), as shown in Figure 3A. The model was significant, as indicated by 100 different model permutations (Q2 = 0.807, p < 0.01 and R2Y = 0.983). The negative mode evaluation for NVLOW and NVHIGH, shown in Figure 3B, revealed good predictive performance (Q2 = 0.568, R2Y = 0.709). The model was significant indicated by 100 different model permutations (p value < 0.01, Q2 = 0.816 and R2Y = 0.934). Compounds significantly responsible for the separation in OPLS-DA models were ascertained using the variable importance in projection (VIP) scores. A VIP > 1.5 was used as a cut-off for variable selection, resulting in a total of 72 significant metabolites in the positive mode and 69 metabolites in the negative mode. To further confirm the discrimination observed between the treatment groups, a linear model (y (metabolite response) ~ return bloom) was applied to the ‘Nicoter’ variable (NVHIGH, NVMID, NVLOW) for both positive and negative datasets. The return bloom data were treated as a covariate for each cultivar. The Benjamini–Hochberg (BH) correction was used to adjust the significance (p value) of each of the variables for the ‘Nicoter’ dataset, and the subsequent adjusted p value was referred to as a Q-value. The majority of significant metabolites in the apple spur buds showed higher levels in NVHIGH treatments than NVLOW, as indicated in Table 1 and Table 2.
The OPLS-DA score plot of ‘Rosy Glow’ Variable UHPLC-HRMS ESI+ data with an associated 95% confidence ellipses demonstrates distinction between RVHIGH and RVLOW with lower predictive performance (Q2 = 0.231, R2Y = 0.647), as shown in Figure 4A. The model was significant indicated by 100 different model permutations (p value < 0.01, Q2 = 0.757, and R2Y = 0.985). The OPLS-DA score plot of ESI- with an associated 95% confidence ellipses demonstrating distinction between RVHIGH and RVLOW, as shown in Figure 4B, revealed poor predictive performance (Q2 = 0.071, R2Y = 0.593). The model was significant indicated by 100 different model permutations (p value < 0.01, Q2 = 0.785 and R2Y = 0.983). A VIP > 1.5 cut-off was also applied to the ‘Rosy Glow’ variable and resulted in a total of 109 significant metabolites in the positive mode and 54 metabolites in the negative mode. A linear model (y (metabolite response) ~ return bloom) was applied to ‘Rosy Glow’ variable (RVHIGH, RVMID, RVLOW) to confirm the significance of the compounds. The return bloom data were treated as a continuous variable (covariate) for each cultivar. p values were only reported for the ‘Rosy Glow’ datasets as no significant Q values were obtained. The majority of significant metabolites in the apple spur buds showed higher levels in RVHIGH treatments than RVLOW, as indicated in Table 3 and Table 4.
MS or MSn fragmentation of parent ion confirmed the identification of differentially expressed metabolites for ‘Nicoter’ (Table 1 and Table 2) and ‘Rosy Glow’ (Table 3 and Table 4) variable treatments. Level 3 identification or above is required for putative identification of compounds in accordance with the Metabolomics Standards Initiative and Schrimpe-Rutledge et al. [17,18]. Most metabolites in Table 1, Table 2, Table 3 and Table 4 had level 2 identification—i.e., compounds that have matching fragmentation pattern with metabolite MS/MS libraries. Level 4 identification requires a unique molecular formula and level 3 requires the precursor m/z to match with a metabolite database. Thus, chlorogenic acid derivatives were assigned level 4 identification as some of the MS/MS ions matched the chlorogenic acid fragmentation data, although the precursor m/z did not match a database.
The majority of the differentially expressed metabolites were elevated in buds that were harvested from trees with high return bloom compared to low return bloom in both ‘Nicoter’ and ‘Rosy Glow’. In ‘Nicoter’, many of the phenylpropanoid pathway intermediates increased, comprising chlorogenic acid and its derivatives (I-VI) and precursor molecule quinic acid, coumarins (7-hydroxycoumarin, coumaranone), salicylates (3,4 dihydroxybenzaldehyde, salicylaldehyde) and a glycosidic derivative (4-acetyl-3-hydroxy-5-methylphenyl β-D-glucopyranoside), flavanols (afzelechin-7-apioside, naringenin, kaempferol) and flavanol glycosides (7-hydroxy-2-(4-hydroxyphenyl)-4-oxo-3,4-dihydro-2H-chromen-5-yl β-D-glucopyranoside; kaempferol-3-O-α-L-rhamnopyranoside). Few compounds including the flavanol (R)-shinanolone, glycosides of hydroxycinnamic acid derivatives 1-O-feruloylglucose and 3-O-acetyl-2-O-coumaroyl-hexopyranose and a methyl salicylate glycoside derivative showed significant decreases in ‘Nicoter’. Similarly, ‘Rosy Glow’ showed increased levels of phenylpropanoid pathway intermediates chlorogenic acid and its derivatives (VIII, IX and 3-caffeoyl-1,5-quinolactone) and precursor molecule quinic acid; coumarins and related derivatives (hydroxycoumarin; coumaric acid; dihydro-trans-o-coumaric acid 2-glucoside, coumaranone), glycosidic derivative of salicylate (2,5-dihydroxybenzoic acid 2-O-β-D-glucoside) and related phenolic glycoside (phlorisobutanophenone glycoside) and flavanol (kaempferol). Sugar alcohol D-(-)-mannitol also increased. The amino acids L-histidine and L-aspartate decreased in ‘Rosy Glow.’
Pathway enrichment analysis showed that compounds benzoate, 1-O-feruloyl-β-D-glucose, kaempferol and trans-5-O-caffeoyl-D-quinate (chlorogenic acid) belong to the phenylpropanoid derivative biosynthesis and this pathway is significantly disrupted (p = 0.005) (Figure 5). Increased levels of chlorogenic acid indicate significant disruption in the chlorogenic acid biosynthesis II pathway (p = 0.04) and the identification of unknown chlorogenic acid derivatives further increases its relevance.
Plant hormones including cytokinins, gibberellins and auxins are believed to participate in floral development in many plant species including apple [19,20]. To investigate the effect of crop load on selected phytohormones and their structural derivatives, relative abundances were measured for the cytokinin precursors, adenine and adenosine, auxins (2-oxindole-3-acetic acid (OxIAA), indole-3-acetic acid (IAA) indole-3-acetonitrile, methyl-indole-3-acetic acid (MeIAA) and precursors tryptophan (Trp) and tryptamine, gibberellic acid (GA), salicylic acid (SA), methyl jasmonate and abscisic acid (ABA). A line graph showing the effect of individual treatments on the compounds is shown for ‘Nicoter’ (Figure S2) and ‘Rosy Glow’ (Figure S3).
A linear model (y (metabolite response) ~ return bloom) was applied to ‘Rosy Glow’ (RVHIGH, RVMID, RVLOW) and ‘Nicoter’ (NVHIGH, NVMID, NVLOW) together with t-tests and fold change in high and low treatments. However, statistics tests revealed there was no significant differences for the selected compounds in ‘Nicoter’ (Table S1) and ‘Rosy Glow’ (Table S2).
Together these data show that increased biosynthesis of phenylpropanoid pathway intermediates, including hydroxycinnamates (chlorogenic acid, coumarates, ferulates), coumarins, salicylates and flavanols, increased in response to low crop loads.

3. Discussion

This study investigated the metabolic pathways of FI in apple using ESI LCMS metabolic profiling on apple buds collected within the critical 70 DAFB time period. In this investigation, levels of hydroxycinnamic acid derivatives, flavonoids and salicylates were significantly increased in buds collected from trees with low crop load (“OFF”) that expressed high return bloom in the following season. Although previous multi-omics investigations have identified differentially expressed genes associated with plant hormone signal transduction and significant metabolites such as flavanols and chlorogenic acid in “OFF” trees, there are no direct reports on the type of plant hormones regulating early flower bud development in apple [3,6,16].
To elucidate the metabolic pathways and key repressors and promotors of FI, young apple trees were forced into biennial bearing. Return bloom response for the 2019/2020 season was coupled to metabolomics profiling of apple spur buds collected in the previous 2018/2019 season using positive and negative ESI LCMS analysis to determine the physiological pathways disrupted during bud meristem development. Our results show that trees exhibiting an “OFF” season triggered metabolome changes that not only corroborated with previous studies but also identified novel candidates for FI belonging to the salicylate group of phytohormones and metabolites within the SA biosynthetic pathway. Biosynthetic production of SA occurs via isochorismate or phenylalanine in the model plant Arabidopsis thaliana [21]. Several metabolites were significantly disrupted by biennial bearing treatments imposed in this study, particularly those represented in Figure 5, providing evidence that the phenylalanine ammonia-lyase (PAL) pathway is likely triggered in the biosynthesis of SA and modification can render SA inactivate (e.g. glycosylation) or confer complementary properties such as activation of stress response via hydroxylation [21].
Guitton et al. described that the differentially expressed transcripts of apple trees in “ON” and “OFF” years showed that the redox and hormonal statuses are likely to contribute to FI in apple trees [16]. Disruption in the redox status of the plant is largely caused by reactive oxygen species (ROS) which are known to be toxic by-products of metabolic processes in normal conditions [22,23]. Under abiotic and biotic stress, there is a boost in production of ROS species which is believed to activate defense genes and the biosynthesis of SA for the activation of a defense response [22,24,25]. Redox homeostasis is maintained by enzymes and metabolites, and ascorbic acid is in the main line of defense for the detoxification of ROS species and leads to the biosynthesis of plant hormones, including SA [24]. Moreover, flavonoids are also considered as radical scavengers due to their antioxidative properties. In this study, increased levels of flavonoids kaempferol, naringenin and afzelechin were observed in “OFF” trees of the biennial bearing cultivar ’Nicoter’. Increased levels of kaempferol were also observed in ’Rosy Glow’. The sugar alcohol D-Mannitol can also function as a radical scavenger and was augmented exclusively in “OFF” trees of ’Rosy Glow‘. Sugar alcohols are produced in plants in response to abiotic and biotic stress and can provide plants with salinity tolerance, efficient growth and pathogen resistance [26].
The regulation of redox potential via antioxidants is involved in the signalling networks in both spatial and temporal dimensions of plant growth and development. ROS are strongly linked to the response to environmental factors, particularly stress and crosstalk with plant hormonal signalling pathways such as salicylates which are known to regulate floral transition [27,28]. Induction of early flowering ensures survival under stress [29].
The high-levels of hydroxycinnamic acid derivatives represented as chlorogenic acid, coumarate and ferulate derivatives in “OFF” trees indicated disruption in the phenylpropanoid biosynthesis pathway which was consistent with a previous study conducted by Milyaev et al. [6]. A gene homologue to caffeic acid metabolism, a likely metabolic precursor to the chlorogenic acid pathway, were also down-regulated in “ON” olive trees in a transcriptomic study [30]. Levels of salicylates were also observed to increase accordingly in the biennial cultivar ‘Nicoter’, resulting from decarboxylation of cinnamates. The elevated levels of chlorogenic acid and its derivatives suggest it is likely a key intermediate in FI, although it may not directly induce flowering [31]. Increased levels of salicylates are known to impact the flowering process in model plants such as duckweeds [27,28,31,32] and Arabidopsis thaliana [33,34], suggesting that a stress activation pathway is likely exerted in the bud, although benzoic acid was shown to elicit a more effective flower induction response than SA in duckweed [27]. Although SA (2-hydroxybenzoic acid) was not differentially expressed in this study, the closely related derivatives 3,4 dihydroxybenzaldehyde and salicylaldhehyde showed significantly increased levels in the biennial bearing cultivar ‘Nicoter’ and 2,5-dihydroxybenzoic acid 2-O-β-D-glucoside in ‘Rosy Glow’.
The metabolite 4-hydroxycoumarin is also known to promote flowering in duckweed, and both cultivars ‘Nicoter’ and ‘Rosy Glow’ showed a significant increase in a hydroxycoumarin (unknown position of hydroxy group) [35]. Although there are no reports of flower-inducing activity of the closely related compound coumaranone, it may play an important role in FI, due to its significance and effect size reported in both cultivars.
‘Rosy Glow’—the less susceptible cultivar to biennial bearing—showed low levels of amino acids histidine and aspartate in “OFF” trees, which are likely consumed to produce the dipeptide aspartyl-histidine which subsequently increased. In Arabidopsis thaliana, Histidine-to-Aspartate phosphorelays are involved in signal transduction induced by cytokinin and other environmental stimuli [36]. Cytokinins are a group of phytohormones that influence growth and the stimulation of cell division and are abundant in the root tip, shoot apex and all plant tissue. Although, cytokinins were not differentially expressed in this study, the significantly increased levels of aspartyl-histidine suggest that the His-Asp phosphorelays are likely a downstream effect of the perceived cytokinin signal [36]. Cytokinins are also responsible for activating the biosynthetic pathway for SA in plants [37] and when applied to apple buds post bloom at the time of fruit bud initiation, an increase in the amount of return bloom occurs [20,38]. The participation of the cytokinin class of phytohormones also corroborates with gene expression studies conducted by Milyaev et al. [6]. It is worth noting that, although no significant effects were observed in the auxin indole acetic acid (IAA) and gibberellic acid (GA), it is possible that phytohormone signalling cross-talk is occurring through mobile signals from the developing fruitlet seed or stem. Studies have indicated that fruit presence stimulates IAA transport in citrus and olive stem [39] and possibly induces GA biosynthesis. In apple, fruitlet seeds are rich in GA’s and application of GA is known to inhibit flowering [12,40].
Using ESI LCMS profiling, and MSMS fragmentation techniques, only 14% of the differentially expressed metabolites in ‘Nicoter’ and ‘Rosy Glow’ cultivars were identified. Albeit these were the most significant, further annotation of the “unknown” metabolites would be required to characterise the full extent of the metabolic pathways driving biennial bearing.
In conclusion, to our knowledge, this is the first functional metabolomics analysis to identify plant hormones associated with FI in apple. This study showed that crop load treatments exerted significant effects on members of the phenylpropanoid pathway comprising hydroxycinnamates, salicylates, salicylic acid biosynthetic pathway intermediates and flavanols. Our findings provide evidence that the PAL salicylic acid biosynthetic pathway was activated in response to an “OFF” year during a biennial bearing cycle. Although no significant changes were exhibited in salicylic acid, its biosynthetic derivatives exhibited distinct increases in “OFF” trees of ‘Nicoter’ and ‘Rosy Glow’ apples. An additional mechanism of cytokinin involvement via histidine-aspartate phosphorelays’, which is also known to activate a defense response in trees, was observed in ‘Rosy Glow’. This study for the first time identifies the participation of the salicylate group of plant hormones in FI in apple.

4. Materials and Methods

The experiment was conducted in a commercial farm at Three Bridges, Yarra Valley, Victoria, Australia. Three-year-old trees of the cultivars ‘Nicoter’ (marketed as Kanzi®) and ‘Rosy Glow’ (marketed as Pink Lady®), trained on Open Tatura trellis, were used. Trees were managed according to the standard local practice and commercial operations. Five crop load treatments were first applied during the 2015–2016 growing season with six replicates, for a total of 30 trees per cultivar. Crop load treatments consisted of 1%, 50%, 100%, 150% and 200% of normal grower practice (6 fruit/cm2 for ‘Nicoter’, 8 fruit/cm2 for ‘Rosy Glow’), based on the tree’s trunk cross-sectional area (TCSA), which was measured 25 cm above the grafting union at the beginning of each growing season. In 2016–2017 and subsequent seasons, three replicates of each crop load treatment alternated between corresponding low and high treatments (e.g., 1% became 200%, 150% became 50%), thus forcing a biennial-type cropping behaviour on those trees (NV = ‘Nicoter’ variable crop load; RV = ‘Rosy Glow’ variable crop load).
During each season at full bloom (i.e., 80% open flowers), the number of flower clusters on each tree were counted manually to determine return bloom (RB). Flower clusters were then thinned by hand in an attempt to obtain the required number of fruit per tree. Only one flower per cluster was retained, except where insufficient clusters meant higher fruit numbers per cluster were needed to achieve the desired number of fruit per tree. Thinning was completed within 4 weeks of full bloom to minimize the chance of excess fruit affecting the following year’s return bloom.
The continuous return bloom variable was converted to a categorical dataset for ‘Nicoter’ and ‘Rosy Glow’ treatments denoted as: NVHIGH (1–2 fruit/cm2 TCSA; 300–400 RB; n = 14) and RVHIGH (1–4 fruit/cm2 TCSA; 190–240 RB; n = 13) = low crop load treatments, eliciting a high RB referred to as “OFF” trees. NVLOW (6–8 fruit/cm2 TCSA; 20–41 RB; n = 14) and RVLOW (11–15 fruit/cm2 TCSA; 32–70 RB; n = 14) = high crop load trees, eliciting a low return bloom, referred to as “ON” trees. Moderate treatments include NVMID (4.43 fruit/cm2 TCSA; 139.7 RB; n = 14) and RVMID (8.77 fruit/cm2 TCSA; 138 RB; n = 14).

4.1. Bud Selection and Preparation

In this study, buds were collected from the apple trees after thinning, in late spring and early summer of the 2018/2019 growing season. Weekly collection of one bud per tree began 4 weeks after full bloom and continued for 8 weeks. Buds were selected on spurs growing on at least two-year-old wood on lateral branches of the trees. Buds were prepared for molecular analysis by removing the scales from the bud using a scalpel then excising the growing tip containing the active meristem and immediately freezing it in dry ice in situ. Buds from each treatment were pooled into one 2 mL polypropylene mini centrifuge tube to ensure sufficient material for analysis. On return to the laboratory, these samples were stored at −80 °C until analysis.

4.2. Extraction of Metabolites from Buds

Samples collected in the first week were used for optimising the grinding and extraction methodology. Apple buds were lyophilised and subsequently one large (3.5–4.1 mm) and two small (2.8–3.2 mm) YTZP (yttria zirconia) beads were added to each tube. Samples were placed in 24-well cryo-blocks on a Geno/Grinder 2010 (SPEX Sample Prep, Metuchen, NJ, USA) and buds were ground at 1200 rpm for 1 min. The samples were extracted with 80% methanol/water (v/v), with extraction volumes adjusted proportionally to the weight of the lyophilised bud. Samples were centrifuged at 13,000 rpm for 2 min and 200 μL of the supernatant was transferred into a HPLC tube and stored at −20 °C until ready for LCMS analysis.

4.3. LCMS Methods for Untargeted and Targeted Analysis

For untargeted metabolite profiling, a Vanquish ultra-high performance liquid chromatography (UHPLC) system (Thermo Fisher Scientific, Bremen, Germany) with a binary pump, autosampler and temperature-controlled column compartment, coupled with a QExactive (QE) Plus mass spectrometer (Thermo, Bremen, Germany) with electrospray (ESI) probe operating in both positive and negative modes, was used. For MS data acquisition, positive and negative ion data were captured over a mass range of 80–1200 m/z, with a mass resolution set at 35,000. Samples were randomised, and blanks (80% methanol) injected every five samples. A pooled biological quality control (PBQC) was run every 10 samples. For MS2, data were acquired in full-scan MS/data-dependent MS2 (ddMS2) mode on positive and negative ionisation modes on selected samples. MS cycles were composed of 1 Full MS and up to 10 ddMS2. Ions within the inclusion list detected in the full MS survey scan (intensity threshold 1.6 × 105) triggered a MS2 event at the peak apex with an isolation window of 0.4 m/z. A 5.0 s delay was required for the same ion to trigger a new MS2 event (dynamic exclusion). Full MS scans were acquired from m/z 100 to 1500 for the positive ionisation mode and 80 to 1200 for the negative ionisation mode with a resolution of 35,000 (full width at half maximum, FWHM, at m/z 200); automatic gain control (AGC) target was 3 × 106; maximum injection time (IT) 200 ms. Scans (ddMS2) were acquired at a resolution of 17,500, the AGC target was 1 × 105 and the maximum IT was 50 ms. Ions were fragmented with stepped collision energy (20, 40 and 60%).
Prior to data acquisition, the system was calibrated with Pierce LTQ Velos ESI Positive and Negative Ion Calibration Solution (Thermo Fisher Scientific). Mass spectrometry data were acquired using Thermo Xcalibur V. 2.1 (Thermo Fisher Scientific Inc., Waltham, MA, USA). Nitrogen was used as the sheath, auxiliary and sweep gases at flow rates of 28, 15 and 4 L/min, respectively. Spray voltage was set at 4000 V (positive and negative).
A Thermo Fisher Scientific Hypersil Gold 1.9 μm, 100 mm × 2.1 mm column with a gradient mobile phase consisting of 0.1% formic acid in H2O (A) and 0.1% formic acid in acetonitrile (B), at a flow rate of 0.3 mL/min was used. The gradient began at 2% B, increasing to 100% B over 11 min; followed by 4 min at 100% B before a 5 min equilibration with 2% B.

4.4. Data Processing and Statistical Analyses

The data files obtained following LCMS analyses were processed in the Refiner MS module of Genedata Expressionist® 12.0 with the following parameters: (1) chromatogram chemical noise subtraction with removal of peaks with less than 4 scans, chromatogram smoothing using moving average estimator over 5 scans and 30% quantile over 151 scans for noise subtraction, (2) intensity thresholding using a clipping method and a threshold of 100,000, (3) selection of positive mode data only, (4) chromatogram RT alignment using a pairwise alignment-based tree and a maximum RT shift of 1 min, (5) chromatogram peak detection using a 5 scan summation window, a minimum peak size of 0.1 min, a maximum merge distance of 0.05 Da, a boundary merge strategy and a maximum gap/peak ratio of 70% with moving average smoothing over 10 scans for peak RT splitting, (6) a chromatogram isotope clustering using RT and m/z tolerance of 0.05 min and 0.05 Dalton, respectively, with a maximum charge of 2 and finally (7) an adduct detection using mainly M + H and allowable adducts (M + 2H, M + K, M + Na, M − H2O + H).
Statistical analyses were performed using the Analyst module of Genedata Expressionist® 12.0. Principal component analyses (PCA) were performed to identify tissue and treatment differences. Overlay of the PBQC and samples allowed for the validation of the high-quality dataset by ensuring RT variation, mass error and sensitivity changes throughout. Identification of metabolites was performed by searching experimental MS1 data through the following databases: Plant Metabolic Network (PMN) https://plantcyc.org (accessed on 23 January 2022); Human Metabolome DataBase (HMDB) (http://hmdb.ca) (accessed on 10 June 2021); ChemSpider (http://chemspider.com) (accessed on 13 June 2021); and Lipid Maps® (http://www.lipidmaps.org) (accessed on 20 June 2021). MS2 data were searched on MzCloud (https://www.mzcloud.org) (accessed on 30 June 2021). Identified significant compounds were inputted into a SmartTable in Plant Metabolic Network (PMN) https://plantcyc.org (accessed on 23 January 2022) for pathway enrichment analysis specific to Malus domestica.
A linear model (y (metabolite response) ~ return bloom) applied to the variable crop load treatments of the individual ‘Nicoter’ and ‘Rosy Glow’ cultivars, revealed significant metabolites. The Benjamini–Hochberg (BH) correction criteria were used to adjust the significance (p value) of each of the variables, and the subsequent adjusted p value is referred to as a Q-value.
Prior to multivariate analysis, the missing value imputation was applied and features with >10% missing values were removed and remaining missing values imputed by k-nearest neighbour (KNN) for each feature. Subsequently, a cube root transformation and autoscaling was applied to the data to achieve normality and homoscedasticity. An OPLSDA model was applied to each dataset using MetaboAnalyst 3.053 41 and Q2 value ≥ 0.4 indicates a model with good predictability.
For targeted analysis of plant hormones and structurally related compounds as described in Table 5, an extracted ion chromatogram (EIC) with a 5 ppm tolerance of [M + H] + was utilised to obtain the relative abundances from the MS spectra, in LCQUAN™ Quantitative Software (Thermo Fisher Scientific). All compounds were purchased from Sigma-Aldrich (St. Louis, MI, USA). Standards were prepared in 80% water/methanol.
Targeted analytes which resulted in p < 0.05 from a two-sided Student t-test and exhibited a log10 fold change of >0.5 in high and low treatments of ‘Nicoter’ and ‘Rosy Glow’ were regarded as significant.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/ijms23116019/s1.

Author Contributions

P.R., T.P., D.S., I.G. and S.R. conceptualised, supervised and resourced the research; P.R., T.P., V.E. and D.S. carried out the main experiments; P.R., T.P., V.E., A.S., I.G. and S.R. carried out the writing, editing and review of the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This study was a component of the apple and pear industry’s PIPS3 (Productivity, Irrigation, Pests and Soils) program of research and development funded by Hort Innovation (AP19003), using the Hort Innovation Apple and Pear research and development levy, contributions from the Australian Government and co-investment from Agriculture Victoria (RCM 6472). Hort Innovation is the grower-owned, not-for-profit research and development corporation for Australian horticulture.

Institutional Review Board Statement

Not Applicable.

Data Availability Statement

The datasets generated during the current study are available from the corresponding author upon reasonable request.

Acknowledgments

The authors would like to thank Henryk Flachowsky and the late Jens Wunsche for allowing us to participate in the project “Physiological, Metabolic and Molecular Basis of Biennial Bearing in Apple” from which the buds samples were sourced, and the Sanders brothers for allowing the experiment in their orchard. We would also like to thank the Molecular Phenomics group at Agriculture Victoria for their assistance in sample preparation and analysis.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. O’Rourke, D. Economic importance of the world apple industry. In The Apple Genome; Korban, S.S., Ed.; Springer International Publishing: Cham, Switzerland, 2021; pp. 1–18. [Google Scholar] [CrossRef]
  2. Monselise, S.P.; Goldschmidt, E.E. Alternate bearing in fruit trees. Hortic. Rev. 1982, 4, 128–173. [Google Scholar] [CrossRef]
  3. Guitton, B.; Kelner, J.J.; Velasco, R.; Gardiner, S.E.; Chagné, D.; Costes, E. Genetic control of biennial bearing in apple. J. Exp. Bot. 2012, 63, 131–149. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Embree, C.G.; Myra, M.T.D.; Nichols, D.S.; Wright, A.H. Effect of blossom density and crop load on growth, fruit quality, and return bloom in ‘Honeycrisp’ apple. HortScience 2007, 42, 1622–1625. [Google Scholar] [CrossRef]
  5. Atay, E.; Crété, X.; Loubet, D.; Lauri, P.-E. Effects of different crop loads on physiological, yield and fruit quality of ‘JoyaTM’ apple trees: High crop load decreases maximum daily trunk diameter and does not affect stem water potential. Int. J. Fruit Sci. 2021, 21, 955–969. [Google Scholar] [CrossRef]
  6. Milyaev, A.; Kofler, J.; Klaiber, I.; Czemmel, S.; Pfannstiel, J.; Flachowsky, H.; Stefanelli, D.; Hanke, M.-V.; Wünsche, J.-N. Toward systematic understanding of flower bud induction in apple: A multi-omics approach. Front. Plant Sci. 2021, 12, 286. [Google Scholar] [CrossRef]
  7. Arseneault, M.H. Crop Load Management of Apple by Chemical Fruitlet Thinning, Preharvest Drop Management, and Investigation of Abscission Physiology. Ph.D. Thesis, University of Guelph, Guelph, ON, Canada, 2016. [Google Scholar]
  8. Kviklys, D.; Samuolienė, G. Relationships among the rootstock, crop load, and sugar hormone signaling of apple tree, and their effects on biennial bearing. Front. Plant Sci. 2020, 11, 1213. [Google Scholar] [CrossRef]
  9. Nitsch, C.; Nitsch, J.P. The induction of flowering in vitro in stem segments of Plumbago indica L.—I. The production of vegetative buds. Planta 1967, 72, 355–370. [Google Scholar] [CrossRef]
  10. Nitsch, C.; Nitsch, J.P. The induction of flowering in vitro in stem segments of Plumbago indica L.—II. The production of reproductive buds. Planta 1967, 72, 371–384. [Google Scholar] [CrossRef]
  11. Hanke, M.-V.; Henryk, b.; Flachowsky, H.; Peil, A.; Hättasch, C. No flower no fruit—Genetic potentials to trigger flowering in fruit trees. Genes Genomes Genom. 2007, 1, 1–20. [Google Scholar]
  12. Dennis, F.G.; Neilsen, J.C. Physiological factors affecting biennial bearing in tree fruit: The role of seeds in apple. HortTechnology 1999, 9, 317–322. [Google Scholar] [CrossRef] [Green Version]
  13. Harley, C.P.; Magness, J.; Masure, M.; Fletcher, L.; Degman, E. Investigations on the Cause and Control of Biennial Bearing of Apple Trees; United States Department of Agriculture: Washinton, DC, USA, 1942.
  14. Ramírez, H.; Benavides, A.; Robledo, V.; Alonso, R.; Gómez, J. Gibberellins and cytokinins related to fruit bud initiation in apple. Acta Hortic. 2004, 636, 409–413. [Google Scholar] [CrossRef]
  15. Hoad, G.V. The role of seed derived hormones in the control of flowering in apple. Acta Hortic. 1978, 80, 93–104. [Google Scholar] [CrossRef]
  16. Guitton, B.; Kelner, J.-J.; Celton, J.-M.; Sabau, X.; Renou, J.-P.; Chagné, D.; Costes, E. Analysis of transcripts differentially expressed between fruited and deflowered ‘Gala’adult trees: A contribution to biennial bearing understanding in apple. BMC Plant Biol. 2016, 16, 55. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  17. Schrimpe-Rutledge, A.C.; Codreanu, S.G.; Sherrod, S.D.; McLean, J.A. Untargeted metabolomics strategies—Challenges and emerging directions. J. Am. Soc. Mass Spectrom. 2016, 27, 1897–1905. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  18. Sumner, L.W.; Amberg, A.; Barrett, D.; Beale, M.H.; Beger, R.; Daykin, C.A.; Fan, T.W.M.; Fiehn, O.; Goodacre, R.; Griffin, J.L.; et al. Proposed minimum reporting standards for chemical analysis: Chemical Analysis Working Group (CAWG) Metabolomics Standards Initiative (MSI). Metabolomics 2007, 3, 211–221. [Google Scholar] [CrossRef] [Green Version]
  19. Samuolienė, G.; Čeidaitė, A.; Sirtautas, R.; Duchovskis, P.; Kviklys, D. Effect of crop load on phytohormones, sugars, and biennial bearing in apple trees. Biol. Plant. 2016, 60, 394–400. [Google Scholar] [CrossRef]
  20. Ramírez, H. Endogenous Cytokinins and Fruit Bud Formation in Apple. pp. 245–248. Available online: https://www.actahort.org/books/514/514_28.htm (accessed on 24 May 2021).
  21. Dempsey, D.M.A.; Vlot, A.C.; Wildermuth, M.C.; Klessig, D.F. Salicylic acid biosynthesis and metabolism. Arab. Book 2011, 9, e0156. [Google Scholar] [CrossRef] [Green Version]
  22. Das, K.; Roychoudhury, A. Reactive oxygen species (ROS) and response of antioxidants as ROS-scavengers during environmental stress in plants. Front. Environ. Sci. 2014, 2, 53. [Google Scholar] [CrossRef] [Green Version]
  23. Asada, K. Production and scavenging of reactive oxygen species in chloroplasts and their functions. Plant Physiol. 2006, 141, 391–396. [Google Scholar] [CrossRef] [Green Version]
  24. Barth, C.; Moeder, W.; Klessig, D.F.; Conklin, P.L. The timing of senescence and response to pathogens is altered in the ascorbate-deficient Arabidopsis mutant vitamin c-1. Plant Physiol. 2004, 134, 1784–1792. [Google Scholar] [CrossRef] [Green Version]
  25. Levine, A.; Tenhaken, R.; Dixon, R.; Lamb, C. H2O2 from the oxidative burst orchestrates the plant hypersensitive disease resistance response. Cell 1994, 79, 583–593. [Google Scholar] [CrossRef]
  26. Stoop, J.M.H.; Williamson, J.D.; Pharr, D.M. Mannitol metabolism in plants: A method for coping with stress. Trends Plant Sci. 1996, 1, 139–144. [Google Scholar] [CrossRef]
  27. Khurana, J.P.; Cleland, C.F. Role of salicylic acid and benzoic acid in flowering of a photoperiod-insensitive strain, lemna paucicostata LP6. Plant Physiol. 1992, 100, 1541–1546. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  28. Fujioka, S.; Yamaguchi, I.; Murofushi, N.; Takahashi, N.; Kaihara, S.; Takimoto, A.; Cleland, C.F. The role of benzoic acid and plant hormones in flowering of Lemna gibba G3. Plant Cell Physiol. 1985, 26, 655–659. [Google Scholar] [CrossRef]
  29. Kocsy, G.; Tari, I.; Vanková, R.; Zechmann, B.; Gulyás, Z.; Poór, P.; Galiba, G. Redox control of plant growth and development. Plant Sci. 2013, 211, 77–91. [Google Scholar] [CrossRef] [PubMed]
  30. Dastkar, E.; Soleimani, A.; Jafary, H.; de Dios Alche, J.; Bahari, A.; Zeinalabedini, M.; Salami, S.A. Differential expression of genes in olive leaves and buds of ON- versus OFF-crop trees. Sci. Rep. 2020, 10, 15762. [Google Scholar] [CrossRef]
  31. Hatayama, T.; Takeno, K. The metabolic pathway of salicylic acid rather than of chlorogenic acid is involved in the stress-induced flowering of Pharbitis nil. J. Plant Physiol. 2003, 160, 461–467. [Google Scholar] [CrossRef]
  32. Wada, K.C.; Takeno, K. Stress-induced flowering. Plant Signal. Behav. 2010, 5, 944–947. [Google Scholar] [CrossRef]
  33. Martínez, C.; Pons, E.; Prats, G.; León, J. Salicylic acid regulates flowering time and links defence responses and reproductive development. Plant J. 2004, 37, 209–217. [Google Scholar] [CrossRef]
  34. Jin, J.B.; Jin, Y.H.; Lee, J.; Miura, K.; Yoo, C.Y.; Kim, W.Y.; Van Oosten, M.; Hyun, Y.; Somers, D.E.; Lee, I.; et al. The SUMO E3 ligase, AtSIZ1, regulates flowering by controlling a salicylic acid-mediated floral promotion pathway and through affects on FLC chromatin structure. Plant J. 2008, 53, 530–540. [Google Scholar] [CrossRef] [Green Version]
  35. Kaihara, S.; Takimoto, A. Induction of flowering in lemna paucicostata by dicoumarol and 4-hydroxycoumarin. Plant Cell Physiol. 1985, 26, 1465–1471. [Google Scholar]
  36. Oka, A.; Sakai, H.; Iwakoshi, S. His-Asp phosphorelay signal transduction in higher plants: Receptors and response regulators for cytokinin signaling in Arabidopsis thaliana. Genes Genet. Syst. 2002, 77, 383–391. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. Hare, P.D.; Cress, W.A.; Van Staden, J. The involvement of cytokinins in plant responses to environmental stress. Plant Growth Regul. 1997, 23, 79–103. [Google Scholar] [CrossRef]
  38. Nichols, D.; Embree, C.; Cline, J.; Ju, H.Y. Blossom and fruitlet thinners affect crop load, fruit weight, seed number, and return bloom of ‘Northern Spy’ apple. HortScience 2004, 39, 1309–1312. [Google Scholar] [CrossRef] [Green Version]
  39. Haim, D.; Shalom, L.; Simhon, Y.; Shlizerman, L.; Kamara, I.; Morozov, M.; Albacete, A.; Rivero, R.M.; Sadka, A. Alternate bearing in fruit trees: Fruit presence induces polar auxin transport in citrus and olive stem and represses IAA release from the bud. J. Exp. Bot. 2021, 72, 2450–2462. [Google Scholar] [CrossRef] [PubMed]
  40. Tu, Y.; Hirst, P.; Coolbaugh, R.; Pharis, R. 532 Endogenous Gibberellins in Developing Apple Seeds in Relation to Biennial Bearing. HortScience 2000, 35, 487B-487. [Google Scholar] [CrossRef] [Green Version]
  41. Xia, J.; Sinelnikov, I.V.; Han, B.; Wishart, D.S. MetaboAnalyst 3.0-making metabolomics more meaningful. Nucleic Acids Res. 2015, 43, W251–W257. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Flower clusters counted on trees in the 2019/2020 season against the crop loads (fruit/cm2 tree’s trunk cross sectional area (TCSA)) of those trees in the 2018/2019 season for: (A) ‘Nicoter’; (B) ‘Rosy Glow’. Each point is the average of three trees with standard error bars.
Figure 1. Flower clusters counted on trees in the 2019/2020 season against the crop loads (fruit/cm2 tree’s trunk cross sectional area (TCSA)) of those trees in the 2018/2019 season for: (A) ‘Nicoter’; (B) ‘Rosy Glow’. Each point is the average of three trees with standard error bars.
Ijms 23 06019 g001
Figure 2. Principal component analysis (PCA) score plot showing ‘Nicoter’ (n = 35) and ‘Rosy Glow’ (n = 33) apple spur bud extracts evaluated on (A) ESI+ UHPLC-HRMS and (B) ESI- UHPLC-HRMS with a pooled biological quality control (PBQC) and categorised by (A1,B1) variety and (A2,B2) collection date.
Figure 2. Principal component analysis (PCA) score plot showing ‘Nicoter’ (n = 35) and ‘Rosy Glow’ (n = 33) apple spur bud extracts evaluated on (A) ESI+ UHPLC-HRMS and (B) ESI- UHPLC-HRMS with a pooled biological quality control (PBQC) and categorised by (A1,B1) variety and (A2,B2) collection date.
Ijms 23 06019 g002
Figure 3. Orthogonal partial least squares discriminant analysis (OPLS-DA) of NVLOW and NVHIGH extracts acquired in UHPLC-HRMS. NVLOW (+) and NV HIGH (Δ) OPLS-DA score plot for (A) ESI+ mode with an associated 95% confidence ellipses (Q2 = 0.601, R2Y = 0.724) and (B) ESI- mode with an associated 95% confidence ellipses (Q2 = 0.568, R2Y = 0.709). Both models were significant indicated by 100 different model permutations for (A) (p < 0.01, Q2 = 0.807 and R2Y = 0.983) and (B) (p < 0.01, Q2 = 0.816 and R2Y = 0.934).
Figure 3. Orthogonal partial least squares discriminant analysis (OPLS-DA) of NVLOW and NVHIGH extracts acquired in UHPLC-HRMS. NVLOW (+) and NV HIGH (Δ) OPLS-DA score plot for (A) ESI+ mode with an associated 95% confidence ellipses (Q2 = 0.601, R2Y = 0.724) and (B) ESI- mode with an associated 95% confidence ellipses (Q2 = 0.568, R2Y = 0.709). Both models were significant indicated by 100 different model permutations for (A) (p < 0.01, Q2 = 0.807 and R2Y = 0.983) and (B) (p < 0.01, Q2 = 0.816 and R2Y = 0.934).
Ijms 23 06019 g003
Figure 4. Orthogonal partial least squares discriminant analysis (OPLS-DA) of RVLOW and RVHIGH extracts acquired in UHPLC-HRMS. RVLOW (+) and RVHIGH (Δ) OPLS-DA score plot of (A) ESI+ mode with an associated 95% confidence ellipses (Q2 = 0.231, R2Y = 0.647) and (B) ESI- mode with an associated 95% confidence ellipses (Q2 = 0.071, R2Y = 0.593). Both models were significant indicated by 100 different model permutations for (A) (p < 0.01, Q2 = 0.757 and R2Y = 0.985) and (B) (p < 0.01, Q2 = 0.785 and R2Y = 0.983).
Figure 4. Orthogonal partial least squares discriminant analysis (OPLS-DA) of RVLOW and RVHIGH extracts acquired in UHPLC-HRMS. RVLOW (+) and RVHIGH (Δ) OPLS-DA score plot of (A) ESI+ mode with an associated 95% confidence ellipses (Q2 = 0.231, R2Y = 0.647) and (B) ESI- mode with an associated 95% confidence ellipses (Q2 = 0.071, R2Y = 0.593). Both models were significant indicated by 100 different model permutations for (A) (p < 0.01, Q2 = 0.757 and R2Y = 0.985) and (B) (p < 0.01, Q2 = 0.785 and R2Y = 0.983).
Ijms 23 06019 g004
Figure 5. Compounds associated with the PAL (Phenylalanine Ammonia-Lyase) pathway to salicylic acid, including (A) hydroxycinnamates, (B) salicylates and (C) flavanols, show significant effects in ‘Nicoter’ and ‘Rosy Glow’ variable treatments associated with meristem bud formation.
Figure 5. Compounds associated with the PAL (Phenylalanine Ammonia-Lyase) pathway to salicylic acid, including (A) hydroxycinnamates, (B) salicylates and (C) flavanols, show significant effects in ‘Nicoter’ and ‘Rosy Glow’ variable treatments associated with meristem bud formation.
Ijms 23 06019 g005
Table 1. Metabolites identified in UHPLC-HRMS ESI+ data that were significant (VIP > 1.5) in the OPLSDA model between NVHIGH compared to NVLOW apple spur buds with associated effect size. Benjamini–Hochberg adjusted p values (Q-values) indicate those metabolites that are significant in the linear model with the (y ~ return bloom) for NVLOW, NVMID and NVHIGH.
Table 1. Metabolites identified in UHPLC-HRMS ESI+ data that were significant (VIP > 1.5) in the OPLSDA model between NVHIGH compared to NVLOW apple spur buds with associated effect size. Benjamini–Hochberg adjusted p values (Q-values) indicate those metabolites that are significant in the linear model with the (y ~ return bloom) for NVLOW, NVMID and NVHIGH.
IdentityRetention Time (min)Mass (m/z) [M + H] +Molecular FormulaMass
Error (ppm)
VIP ScoreEffect Size*Q-Value (BH
Adjusted p Value)
MS2 IonsMetabolite Level
chlorogenic acid3.68355.1019C16H18O9−1.412.52.11.4 × 10−6163.0390, 145.0287, 135.0443, 117.03392
hydroxycoumarin3.68163.0389C9H6O3−0.432.52.01.4 × 10−6135.0442, 107.0495, 95.0498, 79.03912
coumaranone3.68135.0440C8H6O2−0.412.52.01.5 × 10−6117.0339, 107.0491, 89.03862
chlorogenic acid derivative I3.67728.1708--2.32.08.9 × 10−5374.0763, 551.1243, 747.1460, 163.0393,4
chlorogenic acid derivative II3.67559.1093--2.31.92.5 × 10−5188.0710, 163.0394, 145.0287, 135.04444
chlorogenic acid derivative III3.67382.0613--2.31.72.0 × 10−5163.0393, 135.0446, 89.03924
chlorogenic acid derivative IV3.71645.1810--2.32.52.5 × 10−4291.0858, 163.0392, 139.03934
hydroxibenzoisochromanquinone derivative4.15230.0572C13H9O4−0.702.21.49.3 × 10−5147.0442, 119.0495, 91.05453
chlorogenic acid derivative V3.68551.1230--2.21.71.3 × 10−4163.0394, 145.0289, 135.0445, 117.0339, 89.03914
afzelechin 7-apioside5.59407.1336C20H22O9−0.142.11.34.1 × 10−3257.0814, 205.0495, 181.0498, 107.04972
chlorogenic acid derivative VI3.68217.5414- 2.11.8 1.6 × 10−4163.0392, 135.0442, 117.0337, 89.03884
3-O-acetyl-2-O-coumaroyl-hexopyranose4.66369.1180C17H20O9−0.022.11.69.0 × 10−3119.0494, 91.0547, 85.02892
(R)-shinanolone4.66193.0860C11H12O30.402.01.56.3 × 10−3193.0861, 91.0546, 57.03402
7-hydroxy-2-(4-hydroxyphenyl)-4-oxo-3,4-dihydro-2H-chromen-5-yl β-D-glucopyranoside5.43435.1284C21H22O10−0.391.81.31.7 × 10−2271.0615, 151.0028, 119.0493, 107.01282
1-O-feruloylglucose4.19357.1176C16H20O9−1.141.71.63.5 × 10−2177.0546, 147.0443, 137.0598, 119.04952
3,4 dihydroxybenzaldehyde4.04139.0389C7H6O3−0.511.71.23.4 × 10−2111.0444, 93.0339, 83.04962
naringenin5.43273.0757C15H12O5−0.181.61.33.6 × 10−2231.0662, 153.0184, 147.0443, 119.04952
kaempferol-3-O-α-L-rhamnopyranoside5.30433.1125C21H20O10−0.981.61.31.1 × 10−2287.0556, 271.0603, 85.0289, 71.04972
kaempferol5.24287.0547C15H10O6−1.101.61.32.9 × 10−2213.0549, 153.0182, 121.02842
salicylaldehyde4.04123.0440C7H6O2−0.461.61.24.8 × 10−295.0496, 77.0391, 67.0547, 53.03812
* ↑ = up-regulated in NVHIGH, ↓ = down-regulated in NVHIGH.
Table 2. Metabolites identified in UHPLC-HRMS ESI- data that were significant (VIP > 1.5) in the OPLSDA model between NVHIGH compared to NVLOW apple spur buds with associated effect size. Benjamini–Hochberg adjusted p values (Q-values) indicate those metabolites that are significant in the linear model with the (y ~ return bloom) for NVLOW, NVMID and NVHIGH.
Table 2. Metabolites identified in UHPLC-HRMS ESI- data that were significant (VIP > 1.5) in the OPLSDA model between NVHIGH compared to NVLOW apple spur buds with associated effect size. Benjamini–Hochberg adjusted p values (Q-values) indicate those metabolites that are significant in the linear model with the (y ~ return bloom) for NVLOW, NVMID and NVHIGH.
IdentityRetention Time (min)Mass (m/z) [M − H] −Molecular FormulaMass
Error (ppm)
VIP ScoreEffect Size*Q-Value (BH
Adjusted p Value)
MS2 IonsMetabolite Level
quinic acid3.68191.0556C7H12O6−2.682.42.41.4 × 10−8191.0556, 173.0449, 127.0391, 82.02842
chlorogenic acid3.68353.0882C16H18O91.122.32.44.8 × 10−8191.0556, 353.0886, 173.04502
4-acetyl-3-hydroxy-5-methylphenyl β-D-glucopyranoside3.52327.1092C15H20O82.022.21.88.3 × 10−7165.0548, 147.0449, 163.0392, 121.0646, 119.04932
chlorogenic acid derivative VII3.71643.1680--2.22.87.7 × 10−6191.0567, 353.08784
quinic acid isomer4.15191.0556C7H12O6−2.682.22.16.8 × 10−6191.0556, 173.0449, 127.0391, 82.02842
Methyl salicylate glycoside derivative4.33461.1670C20H30O121.192.11.86.1 × 10−4191.0555, 251.0773, 149.0455, 131.0340, 415.16112
* ↑ = up-regulated in NVHIGH, ↓ = down-regulated in NVHIGH.
Table 3. Metabolites identified in UHPLC-HRMS ESI+ data that were significant (VIP > 1.5) in the OPLSDA model between RVHIGH compared to RVLOW apple spur buds with associated effect size. p values indicate those metabolites that are significant in the linear model with the (y ~ returnbloom) for RVLOW, RVMID and RVHIGH.
Table 3. Metabolites identified in UHPLC-HRMS ESI+ data that were significant (VIP > 1.5) in the OPLSDA model between RVHIGH compared to RVLOW apple spur buds with associated effect size. p values indicate those metabolites that are significant in the linear model with the (y ~ returnbloom) for RVLOW, RVMID and RVHIGH.
IdentityRetention Time (min)Mass (m/z) [M + H] +Molecular FormulaMass Error (ppm)VIP ScoreEffect Size*p Value
Linear Model
MS2 IonsMetabolite Level
D-(-)-mannitol1.23183.0862C6H14O6−0.652.31.35.2 × 10−3181.0712, 101.0239, 89.02338, 71.01282
chlorogenic acid derivative VIII3.68775.1544- 2.21.24.5 × 10−1709.0716, 532.0876, 421.0622, 163.03934
chlorogenic derivative IX3.68551.1230- 1.91.65.9 × 10−2585.1142, 374.0770, 255.0183, 163.03944
kaempferol5.24287.0547C15H10O6−1.091.91.32.9 × 10−2213.0549, 153.0182, 121.02842
3-caffeoyl-1,5-quinolactone4.17337.0915C16H16O8−0.871.81.53.2 × 10−3163.0411, 145.0304, 135.0460, 117.03512
chlorogenic derivative X3.67728.1708- 1.81.52.8 × 10−1551.1243, 374.0763, 163.03934
chlorogenic acid3.68355.1019C16H18O9−1.291.61.41.3 × 10−1163.0390, 145.0287, 135.0443, 117.03392
hydroxycoumarin3.68163.0389C9H6O3−0.431.61.31.1 × 10−1135.0442, 107.0495, 95.0498, 79.03912
coumaranone3.68135.0440C8H6O2−0.411.51.31.4 × 10−1117.0339, 107.0491, 89.03862
* ↑ = up-regulated in RVHIGH ↓ = down-regulated in RVHIGH.
Table 4. Metabolites identified in UHPLC-HRMS ESI- data that were significant (VIP > 1.5) in the OPLSDA model between RVHIGH compared to RVLOW apple spur buds with associated effect size. p values indicate those metabolites that are significant in the linear model with the (y ~ returnbloom) for RVLOW, RVMID and RVHIGH.
Table 4. Metabolites identified in UHPLC-HRMS ESI- data that were significant (VIP > 1.5) in the OPLSDA model between RVHIGH compared to RVLOW apple spur buds with associated effect size. p values indicate those metabolites that are significant in the linear model with the (y ~ returnbloom) for RVLOW, RVMID and RVHIGH.
IdentityRetention Time (min)Mass (m/z) [M − H] −Molecular FormulaMass Error (ppm)VIP ScoreEffect Size*p ValueMS2 IonsMetabolite Level
L-histidine1.14154.0613C6H9N3O25.843.12.81.7 × 10−3154.0614, 137.0348, 93.0448, 80.06912
L-aspartate1.20132.0292C4H7NO47.812.01.35.5 × 10−2-3
aspartyl-histidine1.21269.0881C10H14N4O53.882.71.32.1 × 10−3-3
2,5-dihydroxybenzoic acid 2-O-β-D-glucoside2.96315.0730C13H16O92.701.81.43.2 × 10−2-3
dihydro-trans-o-coumaric acid 2-glucoside3.52327.1092C15H20O82.022.21.22.6 × 10−2181.0714, 165.0550, 145.0287, 119.04942
quinic acid3.68191.0556C7H12O62.681.51.41.1 × 10−1191.0556, 173.0449, 127.0391, 82.02842
chlorogenic acid3.68353.0882C16H18O91.121.81.47.9 × 10−2191.0556, 353.0886, 173.04502
phlorisobutanophenone glycoside3.70357.1193C16H22O90.542.01.24.5 × 10−2311.0544, 289.0723, 195.06592
1-O-feruloylglucose3.97355.1041C16H20O91.821.61.49.5 × 10−2177.0546, 147.0443, 137.0598, 119.04952
* ↑ = up-regulated in RVHIGH ↓ = down-regulated in RVHIGH.
Table 5. The molecular formula, accurate mass, [M + H] + and retention time of individual plant hormones and structural derivatives.
Table 5. The molecular formula, accurate mass, [M + H] + and retention time of individual plant hormones and structural derivatives.
NameMolecular FormulaAcurate Mass[M + H] +Retention Time (min)
Methyl (+/−)-jasmonateC13H20O3224.1412224.14128.13
Gibberellic acidC19H22O6346.1416347.14954.65
AdenineC5H5N5135.0545136.06231.32
AdenosineC10H13N5O4267.0968268.10461.33
Indole-3-acetonitrileC10H8N2156.0687157.07666.65
Indole-3-acetic acidC10H9NO2175.0633176.07125.68
2-oxindole-3-acetic acidC10H9NO3191.0582192.06614.42
TryptophanC11H12N2O2204.0899205.09773.51
(+/−)-Abscisic acidC15H20O4264.1362265.14405.95
SAC7H6O3138.0317139.03955.87
Methyl-indole-3-acetic acidC11H11NO2189.0790190.08687.10
TryptamineC10H12N2160.1000161.10783.82
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Reddy, P.; Plozza, T.; Ezernieks, V.; Stefanelli, D.; Scalisi, A.; Goodwin, I.; Rochfort, S. Metabolic Pathways for Observed Impacts of Crop Load on Floral Induction in Apple. Int. J. Mol. Sci. 2022, 23, 6019. https://doi.org/10.3390/ijms23116019

AMA Style

Reddy P, Plozza T, Ezernieks V, Stefanelli D, Scalisi A, Goodwin I, Rochfort S. Metabolic Pathways for Observed Impacts of Crop Load on Floral Induction in Apple. International Journal of Molecular Sciences. 2022; 23(11):6019. https://doi.org/10.3390/ijms23116019

Chicago/Turabian Style

Reddy, Priyanka, Tim Plozza, Vilnis Ezernieks, Dario Stefanelli, Alessio Scalisi, Ian Goodwin, and Simone Rochfort. 2022. "Metabolic Pathways for Observed Impacts of Crop Load on Floral Induction in Apple" International Journal of Molecular Sciences 23, no. 11: 6019. https://doi.org/10.3390/ijms23116019

APA Style

Reddy, P., Plozza, T., Ezernieks, V., Stefanelli, D., Scalisi, A., Goodwin, I., & Rochfort, S. (2022). Metabolic Pathways for Observed Impacts of Crop Load on Floral Induction in Apple. International Journal of Molecular Sciences, 23(11), 6019. https://doi.org/10.3390/ijms23116019

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop