Original Research ARTICLE
A Vegetal Biopolymer-Based Biostimulant Promoted Root Growth in Melon While Triggering Brassinosteroids and Stress-Related Compounds
- 1Department for Sustainable Food Process, Università Cattolica del Sacro Cuore, Piacenza, Italy
- 2Department of Agricultural Sciences, University of Naples Federico II, Portici, Italy
- 3Consiglio per la Ricerca in Agricoltura e l’Analisi dell’Economia Agraria, Centro di Ricerca Orticoltura e Florovivaismo, Pontecagnano Faiano, Italy
- 4NGA Laboratory, Tarragona, Spain
- 5Department of Agricultural and Forestry Sciences, University of Tuscia, Viterbo, Italy
Plant biostimulants are receiving great interest for boosting root growth during the first phenological stages of vegetable crops. The present study aimed at elucidating the morphological, physiological, and metabolomic changes occurring in greenhouse melon treated with the biopolymer-based biostimulant Quik-link, containing lateral root promoting peptides, and lignosulphonates. The vegetal-based biopolymer was applied at five rates (0, 0.06, 0.12, 0.24, or 0.48 mL plant-1) as substrate drench. The application of biopolymer-based biostimulant at 0.12 and 0.24 mL plant-1 enhanced dry weight of melon leaves and total biomass by 30.5 and 27.7%, respectively, compared to biopolymer applications at 0.06 mL plant-1 and untreated plants. The root dry biomass, total root length, and surface in biostimulant-treated plants were significantly higher at 0.24 mL plant-1 and to a lesser extent at 0.12 and 0.48 mL plant-1, in comparison to 0.06 mL plant-1 and untreated melon plants. A convoluted biochemical response to the biostimulant treatment was highlighted through UHPLC/QTOF-MS metabolomics, in which brassinosteroids and their interaction with other hormones appeared to play a pivotal role. Root metabolic profile was more markedly altered than leaves, following application of the biopolymer-based biostimulant. Brassinosteroids triggered in roots could have been involved in changes of root development observed after biostimulant application. These hormones, once transported to shoots, could have caused an hormonal imbalance. Indeed, the involvement of abscisic acid, cytokinins, and gibberellin related compounds was observed in leaves following root application of the biopolymer-based biostimulant. Nonetheless, the treatment triggered an accumulation of several metabolites involved in defense mechanisms against biotic and abiotic stresses, such as flavonoids, carotenoids, and glucosinolates, thus potentially improving resistance toward plant stresses.
A well-developed root system is an important agronomic trait of horticultural crops, with implications on crop productivity, abiotic stress tolerance as well as nutrient uptake and assimilation (Koevoets et al., 2016). In fruiting vegetables (i.e., tomato, eggplant, pepper, melon, and watermelon), the most intense root growth development is usually concentrated during the first phenological stages and decreases over time due to the limited translocation of photosynthates (i.e., soluble sugars) resulting from the strong sink demand for reproductive growth (Heuvelink, 2005). Therefore, a vigorous and well-developed root system during early phenological stages able to explore the soil volume is a primary requisite for securing yield stability under both optimal and suboptimal conditions (Jung and McCouch, 2013).
Breeding and genetic engineering efforts to improve crop productivity are practically focused on shoot traits, whereas the root traits are still not fully exploited source of crop improvement (Wachsman et al., 2015; Koevoets et al., 2016). More recently, a variety of biostimulant substances (i.e., humic and fulvic acids, protein hydrolysates, and seaweed extracts) and microbial inoculants (i.e., mycorrhizal fungi and plant growth promoting rhizobacteria-PGPR) has been introduced as an efficient, safe, and sustainable tool to optimize root system thus boosting crop performance, and nutrient use efficiency as well as enhancing tolerance to environmental stressors (du Jardin, 2015; Colla et al., 2015a, 2017a; Rouphael et al., 2015, 2017b,c). The stimulation of biomass production in response to biostimulant application has been usually associated to the action of specific signaling molecules on plant metabolism and physiology (Colla et al., 2015b, 2017a).
Over the past decade, vegetal-based biopolymers like peptides and lignosulphonates have gained prominence worldwide as biostimulant molecules in vegetable cropping systems (Ertani et al., 2011; Matsumiya and Kubo, 2011; Colla et al., 2013, 2014, 2017b). Matsumiya and Kubo (2011) have shown that bioactive peptides isolated from soybean seeds had phytohormone-like activities since they able to promote the root hair formation (root hair length and numbers) of annual baby’s-breath, cabbage, Italian clover, and lettuce. Furthermore, Colla et al. (2014) also demonstrated a root stimulation effect induced by a commercial legume seeds-derived protein hydrolysate Trainer. The mode of action through which Trainer improves the absorption and transport of nutrients may be at least partially associated to changes in root growth with lateral root formation through an auxin-signaling mediated pathway (Colla et al., 2014, 2015a, 2017a). In addition, also lignosulphonates have been shown to elicit both auxin and gibberellin-like activities incurring a significant increase in maize biomass (Ertani et al., 2011). In their study Ertani et al. (2011) demonstrated that lignosulphonates treatment improved N assimilation in maize plants through the stimulation of key enzymes (glutamate-synthase and glutamine-synthetase) and also by promoting photosynthetic activity through the stimulation of both rubisco enzyme activity as well as chlorophyll biosynthesis. Besides the biostimulant-mediated enhancement of root growth, vegetal-based biopolymers can also upregulate the expression of genes encoding for plasma membrane H+-ATPases and nitrate transporters, resulting in greater efficiency in nitrate uptake and assimilation (Tavares et al., 2017).
Despite the above described effects of biopolymers in plants, conclusive evidences regarding their molecular targets are far from being elucidated. In this perspective, untargeted profiling approaches by metabolomics have proven to be a powerful tool for shedding light on the mode of action of vegetal-based biopolymers (Lucini et al., 2015, 2016; Ntatsi et al., 2017). Indeed, metabolomics allows identifying those compounds altered by treatment, thus driving optimal target uses of biopolymers-based products in agriculture. In turn, this information on the mechanisms through which biostimulants act in plant, can drive a better use of these products by highlighting those scenarios where they provide a real and essential contribution to crop production.
The elucidation of fundamental plant physiological and biochemical responses to vegetal-based biopolymer application can be instrumental to understand the mechanisms behind the well-developed root system; hence the objective of the current study was to assess the morphological, physiological, and metabolomic changes in response to vegetal-based biopolymer application on greenhouse melon. Melon plants treated with vegetal-based biopolymer by substrate drench were compared to untreated plants in terms of biomass production and partitioning, root morphology, SPAD index, chlorophyll fluorescence as well as metabolic profiling.
Materials and Methods
Plant Material, Growth Conditions, Biostimulant Treatments, and Experimental Design
The greenhouse trial was carried out in the spring 2016 growing season in a 300 m2 polyethylene greenhouse at the Experimental Farm of Tuscia University, central Italy (latitude 42°25′N, longitude 12°08′E, altitude 310 m). Transplants of melon (Cucumis melo L. - cv. Giorillo commercialized by Seminis, Milan, Italy) grown in polystyrene plug trays (84 holes) were planted at the second true-leaf stage on 4th May in black plastic pots (15 cm diameter; one plant per pot) containing 1.8 L of sandy soil. Fertilization was performed prior planting by mixing for each kilogram of soil 2 g of a slow mineral fertilizer containing (g kg-1) 150 N, 39.2 P, 124.5 K, 12.1 Mg, 80.0 S, 3.0 Fe, 0.1 B, 0.1 Mn, 0.02 Cu, and 0.02 Zn. Plants were grown under natural light conditions. The greenhouse was maintained at daily temperatures between 18 and 28°C, and day/night relative humidity of 60/85%. Just after transplanting, plants were irrigated with 200 mL of water per pot.
The five biostimulant application treatments were four rates of the commercial Quik-Link product: 0, 0.06, 0.12, 0.24, or 0.48 mL plant-1 (= 0, 0.3, 0.6, 1.2, or 2.4 L ha-1, respectively; the rates per hectare were calculated considering a localized biostimulant placement near the transplants with a plant density of 5,000 plants ha-1). Quik-link is a biopolymer-based biostimulant manufactured by Italpollina S.p.A., Rivoli Veronese, Italy; the product contains Lateral root promoting peptides (LRPP) and other biopolymers (e.g., lignosulphonates) with high biological activity on plants, and micronutrients (10 g kg-1 Fe, 7 g kg-1 Mn, 3 g kg-1 Zn, 1 g kg-1 Cu, 0.2 g kg-1 Mo). Quik-link is allowed in organic agriculture according to the Council Regulation (EC) No. 834/2007 of 28 June 2007.
Treatments were arranged in a randomized block design with four replicates. Each experimental unit consisted of 10 plants. Fifty milliliters of a water solution containing Quik-link product was applied manually after 2 days from transplanting around the collar level of each plant. During the experiment, plants were irrigated as needed with a drip irrigation system.
Growth Measurements and Root Characteristics
At the end of the greenhouse experiment (16th May; 12 days after transplanting-DAT) melon plants were separated into leaves, stems, and roots. All plant tissues were dried at 80°C for 72 h until they reached a constant weight which corresponded to their dry biomasses. Shoot dry weight was equal to the sum of the aerial vegetative parts (leaves + stems), and the root-to-shoot ratio was also calculated. Two plants per experimental plot were used for the root morphology determination. Root system collection and sample preparation were performed following the protocol described previously by Rouphael et al. (2017a). Briefly, the melon root were gently washed with fresh water, until the roots were free from any sandy particles. The determination of the root system architecture components was done using a WinRHIZO Pro (Regent Instruments Inc., Canada), connected to a STD4800 scanner. The following root morphology characteristics were recorded: total root length, mean root diameter, and total root surface area (Colla et al., 2014).
Soil Plant Analysis Development Index and Chlorophyll Fluorescence
Soil Plant Analysis Development (SPAD) index and the fluorescence measurements of melon leaves were also recorded at the end of the experiment. The relative leaf chlorophyll concentration expressed as a rational unit was measured using a portable chlorophyll meter SPAD-502 (Minolta Corporation, Ltd., Osaka, Japan). The SPAD measurements were made on the fully expanded leaves which correspond to the third leaf starting from the apical shoot, as described by Colla et al. (2011). Twenty random readings per experimental unit were taken and averaged to a single SPAD value for each biostimulant application treatment.
On the same date, the chlorophyll fluorescence was measured every 20 min on dark-adapted leaves (two measurements per plant) by means of a chlorophyll fluorometer Handy PEA (Hansatech Instruments Ltd., King’s Lynn, United Kingdom) with excitation source intensity higher than 3,000 μmol m-2 s-1 at the sample surface as described previously by Borgognone et al. (2016). Briefly, the minimal fluorescence intensity (F0) in a dark-adapted state was measured in the presence of a background weak light signal (about 2–3 μmol photons m-2 s-1). The maximal fluorescence level in the dark-adapted state (Fm) was induced by 0.8 s saturating light pulse (3,000 μmol photons m-2 s-1). The maximum quantum yield of open photosystem II (PSII) (Fv/Fm) was calculated as (Fm–F0)/Fm, as also described by Borgognone et al. (2016).
Collection of Samples and Metabolomic Analysis
The first fully expanded leaf (the second leaf from the growing tip) and root (terminal roots) samples were collected from two plants per experimental plot treated with 0 or 0.24 mL plant-1 of the biostimulant Quik-link. Leaf and root tissues were then quenched in liquid nitrogen and grounded into a fine powder using mortar and pestle, thereafter leaf and root samples were stored at -80°C for successive metabolomic analysis.
Samples (1.0 g) of five replicates per treatment were extracted by homogenization in 15 mL of 80% methanol added with 0.1% HCOOH, using an Ultra-Turrax (Ika T-25, Staufen, Germany). Extracts were centrifuged (9000 × g), filtered through a 0.22 μm cellulose syringe filter and then transferred to an amber vial for analysis. Metabolomic analysis was done using a 1290 UHPLC liquid chromatography system coupled to a G6550 quadrupole-time-of-flight mass spectrometer (UHPLC-ESI/QTOF-MS – Agilent Technologies, Santa Clara, CA, United States). The mass spectrometer was equipped with a JetStream dual Electrospray ionization source.
Instrumental conditions were taken from previous experiments (Pretali et al., 2016). Briefly, UHPLC separation was achieved on an Agilent Zorbax Eclipse-plus column (75 mm × 2.1 mm i.d., 1.8 μm) with a mobile phase consisting of water (A) and methanol (B), and a flow of 220 μL min-1 at 35°C. The gradient was operated from 5 to 90% B in 35 min, whereas the QTOF mass spectrometer was set in positive polarity, scan acquisition (100–1200 m/z+) and extended dynamic range mode. Nebulizer pressure was 60 psig, sheath gas was nitrogen at 10 L min-1 (350°C), drying gas was nitrogen at 10 L min-1 (280°C) and capillary voltage was 4 kV.
Features deconvolution, as well as the following mass and retention time alignment, were done in Profinder B.05 (from Agilent Technologies). Annotation was achieved using accurate mass, isotope accurate spacing and isotope ratio against the database PlantCyc 9.5 (Plant Metabolic Network1; released November 2014). Hence, a Level 2 of identification (i.e., putatively annotated compounds) was achieved, with reference to COSMOS Metabolomics Standards Initiative2. Compounds were finally filtered by frequency, retaining only those compounds that were present in 100% of replications within at least one treatment, then exported for statistics.
Statistical Analysis of Experimental Data
Analysis of variance (ANOVA) of the experimental data was made using the SPSS software package (IBM SPSS Statistics version 20.0.0). Orthogonal contrasts (Gomez and Gomez, 1983) were used to compare the biostimulant concentration effects on morphological and physiological parameters. Duncan test was also performed at P = 0.05 on each of the significant variables measured.
Regarding metabolomics, the dataset was interpreted in Agilent Mass Profiler Professional B.12.06 (from Agilent Technologies) as previously reported (Lucini et al., 2017). Compounds abundance was normalized at the 75th percentile and baselined to the median of control following the adoption of a threshold of 10000 counts. Pairwise comparisons were done in Volcano Plot analysis, by combining analysis of variance (P < 0.05, Bonferroni multiple testing correction) and fold-change analysis (cut-off = 5). The dataset was next exported into SIMCA 13 (Umetrics, Malmo, Sweden), pareto-scaled and elaborated for Principal Component Analysis. Thereafter, Orthogonal Projections to Latent Structures Discriminant Analysis (OPLS-DA) modeling was carried out for leaves and roots, separately. OPLS-DA supervised multivariate analysis targeted separating variation between the groups into predictive and orthogonal (i.e., ascribable to technical and biological variation) components. Outliers were excluded according to Hotelling’s T2 and, adopting 95 and 99% confidence limits for suspect and strong outliers, respectively. The OPLS-DA model was validated through cross validation CV-ANOVA (p < 0.01) and overfitting was excluded by permutation testing (n = 100). OPLS-DA goodness-of-fit R2Y and goodness-of-prediction Q2Y were recorded and finally variables importance in projection (VIP analysis) was used to select those having the highest discrimination potential (VIP score > 1.46).
Morphological and Physiological Parameters
In the current study, no significant differences on dry weight of stem, SPAD index and the maximum quantum use efficiency of PSII in dark-adapted state (Fv/Fm) were observed between biostimulant application rates. The mean dry weight of stem, SPAD index and Fv/Fm were 32.9 mg plant-1, 51.2 and 0.82, respectively (Table 1).
TABLE 1. Effect of biopolymer-based biostimulant applications on dry weight of leaves, stems, total biomass, SPAD index, and maximum quantum use efficiency of photosystem II in dark-adapted state (Fv/Fm) of melon plants grown under greenhouse conditions.
The growth parameters in particular the dry weight of leaves as well as the total dry biomass were influenced by the vegetal-based biopolymer applications (Table 1). For instance, increasing the concentration of the biopolymer-based biostimulant from 0.0 to 0.48 mL plant-1 incurred a quadratic behavior of both leaves and total shoot biomass (Figure 1 and Table 1). Our experimental data indicated that biopolymer-based biostimulant at 0.12 and 0.24 mL plant-1 enhanced dry weight of melon leaves and total biomass by 30.5 and 27.7%, respectively, compared to biopolymer applications at 0.06 mL plant-1 and untreated plants, with no significant difference found between 0.12 and 0.24 mL plant-1 rates (Table 1).
FIGURE 1. Twelve-day-old plants of melon treated with 0 (A), 0.06 (B), 0.12 (C), 0.24 (D), or 0.48 (E) mL plant-1 of biopolymer-based biostimulant.
Root Morphological Characteristics
The root-to-shoot ratio as well as the average root diameter were not affected by the vegetal-based biopolymer applications. The mean root-to-shoot ratio and root diameter of melon plants recorded in the current study were 0.25 and 0.36 mm, respectively (Table 2). Similarly, to the effects on plant growth parameters (dry weight of leaves and total biomass), the root dry biomass and the root morphology characteristics (total root length and surface) in biostimulant-treated plants were significantly higher at 0.24 mL plant-1 and to a lesser extent at 0.12 and 0.48 mL plant-1, in comparison to untreated melon plants (Table 2). In fact, the root application of the commercial biopolymer biostimulant Quik-link at a rate of 0.24 mL plant-1 elicited significant increase in 12-days melon transplants root biomass, total root length, and surface amounting to 34.1, 32.0, and 32.6%, respectively, compared to biostimulant applications at 0.06 mL plant-1 and untreated plants (Figure 1 and Table 2).
TABLE 2. Effect of biopolymer-based biostimulant applications on dry weight of roots, root-to-shoot ratio, total root length, average root diameter, and total root surface of melon plants grown under greenhouse conditions.
Metabolic Profiling of Melon Leaves and Roots
The changes in metabolic profile of roots and leaves of melon, following treatment with the biostimulant Quik-link at the rate of 0.24 mL plant-1, was investigated through UHPLC-ESI/QTOF-MS untargeted metabolomics. Overall, more than 2,300 compounds were annotated and passed the filter thresholds adopted. The following unsupervised Principal Component Analysis carried out from metabolomic profiles allowed to identify two main clusters, one per matrix (Figure 2). Even if matrix type was, as expected, the main classification factor, two clearly distinct sub-clusters could be evidenced for root profiles, whereas leaf profiles were partially overlapped. This indicated that the effect of Quik-link was more evident in roots, and suggested the need of more focused multivariate analyses to better point out those compounds being responsible of differences across treatments. The supervised OPLS-DA model, carried out separately for roots and leaves (Figure 3), provided with an excellent separation of control and Quik-link treated samples. In more detail, goodness-of-fit R2Y was 0.62 and 0.65 and prediction ability Q2Y was 0.68 and 0.66, for leaves and roots respectively. Therefore, prediction ability was above the acceptability threshold of 0.5 (Rombouts et al., 2017) and was considered as acceptable. Coherently, cross-validation ANOVA resulted in in a Fischer’s probability of 0.004 for both OPLS-DA models. Furthermore, overfitting was excluded by permutation testing and no outliers could be pointed out by Hotelling’s T2 (Supplementary Material). Provided that differences between treatments were represented in the dataset, Volcano analysis (P < 0.05, fold-change > 5) allowed exporting differential compounds for roots and leaves separately (Tables 3 and 4, respectively). These compounds were ascribed into classes by taking into count their biochemical function and physiological roles.
FIGURE 2. Principal Component Analysis (PCA)-based unsupervised analysis of melon roots and leaves metabolomic profile following root application of 0 (control) or 0.24 mL plant-1 of the biopolymer-based biostimulant.
FIGURE 3. Orthogonal Projections to Latent Structures Discriminant Analysis (OPLS-DA) of melon leaf (A) and root (B) metabolomic profile following root application of 0 (control) or 0.24 mL plant-1 of the biopolymer-based biostimulant. Ellipses denote Hotelling’s T2 (95%).
TABLE 3. Differential metabolites changing in melon root after applying the biopolymer-based biostimulant, as identified through Volcano Plot analysis (P < 0.05, Bonferroni multiple testing correction, and fold-change cut-off = 5).
TABLE 4. Differential metabolites changing in melon leaves after applying the biopolymer-based biostimulant, as identified through Volcano Plot analysis (P < 0.05, Bonferroni multiple testing correction, and fold-change cut-off = 5).
Regarding roots, most of the compounds were related to stress response (i.e., flavonoids, carotenoids, and glucosinolates) rather than hormone network (mainly brassinosteroids and cytokinins), followed by membrane lipids and other compounds. Although a well-defined trend cannot be observed for membrane lipids, flavonoids and carotenoids, glucosinolates were down-accumulated (excepting for 2/3-(5′-methylthio)pentylmalate). Cytokinins were also down-accumulated whereas brassinosteroids were generally up-accumulated (Table 3). VIP analysis from OPLS-DA supervised modeling (Supplementary Material) confirmed the involvement of brassinosteroids, being 6-alpha-hydroxycampestanol (isobaric with 3-epi-6-deoxocathasterone and 6-deoxocathasterone), campest-4-en-3beta-ol and campesterol among the compounds with the highest VIP score. Notably, other top-ranking OPLS-DA VIP score compounds such as the carotenoids zeaxanthin and 15,15′-dihydroxy-beta-carotene, were in common with Volcano plot analysis.
However, Volcano Plot analysis highlighted that the changes in leaves metabolic profile following application of Quik-link was distinct from response in roots (Table 4). A complex and wide alteration of hormone profile, involving among other brassinosteroids, auxin, abscisic acid (ABA), cytokinins, and gibberellins, was observed. The most clear trends were related to the increase of ABA intermediates (excepting the ABA glucose ester), brassinosteroids (excepting the precursor campesterol), and cytokinins. Furthermore, several stress-related compounds were selected by Volcano Plot, namely the oxidative stress-related 3-hexenal and hydroxycaprate, few alkaloids, carotenoids (all of them up-accumulated excepting for bixin) and flavonoids (also found to be up-accumulated). Membrane lipids and sterols were also altered by the treatment, together with compounds from other classes. Among these latter, pyrrole-related (delta1-pyrroline-2-carboxylate and 4-(1-methyl-2-pyrrolidinyl)-3-oxobutanoate methyl ester) and photosynthesis-related (2-methoxy-6-all trans-octaprenyl-2-methoxy-1,4-benzoquinol and a primary fluorescent chlorophyll catabolite) metabolites were pointed out. Discriminating compounds gained by VIP analysis following OPLS-DA confirmed the involvement of brassinosteroids (26-hydroxybrassinolide), membrane lipids (a lysophosphatidylcholine and a digalactosylglycerol), and photosynthesis-related compounds (haematoporphyrin IX, menaquinol-8, and dimethyl phylloquinone) in response to Quik-link treatment (Supplementary Material).
The use of biopolymer-based biostimulant in vegetable cropping systems can promote root growth leading to a better transplant establishment, and higher crop productivity. The current study demonstrated that substrate drench of a biopolymer-based biostimulant elicit dose-dependent (especially at 0.12 and 0.24 mL plant-1) increases of biomass production of melon transplants. A presumed mechanism behind the stimulation of biomass production in response to root application of vegetal-based biopolymers could be the presence of signaling compounds in particular bioactive peptides (e.g., LRPP) as well as lignosulfonates. The former signaling molecules in the Quik-link formulation which is easily perceived by root organ may have triggered a signal transduction pathway through modulation of various endogenous phytohormone biosynthesis (Matsumiya and Kubo, 2011; Colla et al., 2017a). Our findings are in agreement with the results of several research teams (Ertani et al., 2011; Matsumiya and Kubo, 2011; Colla et al., 2013, 2014) who showed that foliar of root applications of small peptides and lignosulfonate-humate elicited hormone-like activities on a wide range of agronomic and horticultural crops (baby’s-breath, cabbage, Italian clover, and lettuce, maize and tomato), thus boosting plant growth and yield.
Another putative mechanism behind the enhancement of plant growth parameters (dry weight of leaves and total biomass) induced by Quik-link application could be attributed to the stimulation of root morphology characteristics in particular the total root length and surface area (at 0.24 mL plant-1 and to a lesser extent at 0.12 and 0.48 mL plant-1), which may improve nutrient uptake and utilization efficiency with beneficial effect on biomass production (Matsumiya and Kubo, 2011). Our results are in line with previous studies assessing the stimulation action of a plant-derived protein hydrolysate containing bioactive peptides and commercial lignosulfonate-humate on root biomass and morphology (Ertani et al., 2011; Colla et al., 2014). For instance, Ertani et al. (2011) demonstrated that the application of lignosulfonate-humates to the nutrient solution (applied at 0.0, 0.5, or 1.0 mg L-1) elicited a dose-dependent increase of root dry weight (from 17 to 24%) in hydroponically-maize compared to the untreated control. Similarly, Colla et al. (2014) showed that the root application of 6 mL L-1 of a plant-derived protein hydrolysate increased root dry weight, length, and surface area of tomato cuttings by 24–35% in comparison to untreated plants.
Quik-link application resulted in a metabolic reprogramming, with roots undergoing a stronger change in biochemical profile. Indeed, despite Volcano Plot analysis highlighted a higher number of differential compounds in leaves, PCA unsupervised multivariate statistics showed that the metabolic profile of roots was markedly different in treated samples, as compared to control. Metabolomics was used to better understand the changes induced by the treatment and to highlight the processes underlying the morphological changes observed. Roots response to Quik-link application involved compounds falling into two main processes, namely stress response and hormone profile.
Among the latter, brassinosteroids appeared to play a pivotal role, with several compounds pointed out both by Volcano Plot analysis (most of them being up-accumulated) and by OPLS-DA VIP analysis. Brassinosteroids are polyhydroxy steroid lactone hormones that are essential for normal plant growth, having effect on various developmental processes of plants (Li and Chory, 1999; Zhu et al., 2013). They are produced from plant sterols, via teasterone, typhasterol, and castasterone, by an isoprenoid biosynthetic pathway that includes acetyl CoA (Ali, 2017). Brassinosteroids signaling pathway leaded to the identification of a putative brassinosteroid receptor as well as brassinosteroid-responsive genes (Li and Chory, 1999; Lin et al., 2017). These hormones are implicated in a wide range of physiological and biochemical responses in plants, including seed germination, cell division and elongation, vascular differentiation, photomorphogenesis, photosynthesis, and senescence (Ali, 2017). They have also been found to protect plants from both abiotic and biotic stress factors, such as salinity, drought, heavy metals, and pathogens (Ali, 2017). Interestingly, beside promoting crop yield, they are reported to regulate plant architecture (Lin et al., 2017), to promote root growth in a dose-dependent manner (Mussig, 2003) and to have an important role in directing epidermal cell fate in roots, regulating differentiation into hair or non-hair cells (Zhu et al., 2013). Therefore, the changes of brassinosteroids in plant tissues can help to explain the increased root development following Quik-link application. Indeed, these hormones display direct interaction with gibberellins signaling pathway, promoting mitotic activity and expression of cyclins in the root meristem (Zhu et al., 2013). In general, brassinosteroids are reported to have a wide interplay with several hormones including an antagonistic role toward ABA (Clouse, 2016) and ethylene biosynthesis (Zhu et al., 2015).
It is described that exogenous brassinosteroids are taken up through the roots to be translocated, unchanged, to the shoot (Symons et al., 2008). However, the same authors reported that, when applied directly to shoot tissues, exogenous brassinosteroids are relatively immobile. Coherently, brassinosteroids were found among differential compounds also in leaves both following Volcano Plot analysis (with 6-alpha-hydroxy-castasterone, cathasterone, and 22alpha-hydroxy-campest-4-en-3-one being up-accumulated) and OPLS-DA VIP analysis (26-hydroxybrassinolide). Several other hormones were found among differential compounds in leaves; the ethylene precursor 1-aminocyclopropane-1-carboxylate and the jasmonate precursor 13-hydroperoxylinoleate were both down-accumulated. Moreover, ABA precursors were up-accumulated, and its storage form beta-D-glucopyranosyl abscisate was down-accumulated. This might suggest that ABA biosynthesis was hampered by brassinosteroids. Furthermore, cytokinins and the most of gibberellins were both up-accumulated, whereas a clear trend could not be observed for auxins. Generally, brassinosteroids induced response appeared to affect the complex interplay and cross-talk between plant hormones, in agreement with previous literature (Li and Chory, 1999; Zhu et al., 2015; Clouse, 2016; Lin et al., 2017).
A correlation between net photosynthesis and brassinosteroids action has been also reported (Ali, 2017); indeed, a benzoquinol, a chlorophyll catabolite, and two pyrrole-related intermediates were up-accumulated in leaves following application of Quik-link. Similarly, hematoporphyrin IX, menaquinol-8, and demethylphylloquinone were among the compounds with the highest VIP score in OPLS-DA.
In a previous study focused on the effect of brassinosteroids on chill injury of fruits and vegetables during post-harvest, Soleimani Aghdam et al. (2016) reported: (i) enhanced membrane integrity, (ii) enhancing antioxidant capacity, and (iii) alteration in phenylalanine ammonia-lyase (PAL) and polyphenol oxidase (i.e., two key enzymes in phenolics biosynthesis and degradation) activities. Similarly, carotenoids were induced in tomato fruits following application of brassinosteroids (Zhu et al., 2015). This might explain the involvement of flavonoids and carotenoids in roots and, more markedly, in leaves (where both were up-accumulated). Nonetheless, these responses related to flavonoids and carotenoids, together with the changes in alkaloids and glucosinolates (two additional classes of metabolites displaying a pivotal role in stress response) are likely connected, at least in part, into the brassinosteroids-dependent response to biotic and abiotic stresses.
The application of the commercial vegetal-based biopolymer Quik-link as substrate drench improved the melon growth parameters in a dose-dependent manner. The improve in plant biomass production was associated to the stimulation of the root growth and morphology traits, thus inducing a ‘nutrient acquisition response’ that favors nutrient uptake and utilization efficiency. Metabolomics allowed to point out a wide and complex biochemical response to the treatment, in which hormone profile, and brassinosteroids in particular, appeared to play a pivotal role. Root metabolic profile was markedly altered following application of the test product. In roots, brassinosteroids could have been responsible of changes in root development. These hormones are known to be transported to shoots, where they might have interfered with other hormones through a complex and wide cross talk that resulted in the hormonal imbalance. Besides potentially shaping plant structure, brassinosteroids were likely responsible of the changes in other secondary metabolites such as phenolics and carotenoids, as well as the modulation of photosynthetic activity. Regardless the eventual direct role of brassinosteroids, the treatment triggered an up-accumulation of several metabolites involved in defense mechanisms against biotic and abiotic stresses, such as flavonoids, carotenoids, and glucosinolates. Although more focused and dedicated experiments are needed in this direction, this could result in an improved resistance toward plant stress.
LL performed the metabolomic study, analyzed the data, and wrote many parts of the manuscript. YR performed the agronomic study, analyzed the data, and wrote many parts of the manuscript. MC performed the root growth measurements and implemented the manuscript. PB gave support in the data analysis and interpretation. CB collaborated in the metabolomic study. GC provided intellectual inputs for defining the experimental design, data analysis and interpretation, and improved the manuscript.
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2018.00472/full#supplementary-material
FIGURE S1 | Hotelling’s T2 investigation on outliers as carried out in UHPLC/QTOF-MS metabolomic investigations on root and leaf metabolic profile following root application of the biopolymer-based biostimulant.
FIGURE S2 | Visual representation of the VIP-selected discriminating compounds from the OPLS-DA modeling on roots (A) and leaves (B) metabolic profile, following root application of the biopolymer-based biostimulant.
TABLE S1 | VIP-selected discriminating compounds from OPLS-DA modeling on root and leaf metabolic profile, following root application of the biopolymer-based biostimulant.
TABLE S2 | Whole dataset gained from UHPLC/QTOF-MS metabolomic investigations on melon metabolic profile following root application of the biopolymer-based biostimulant. Compounds abundance, together with retention time and composite spectra (mass-abundance pairwise) are provided.
Borgognone, D., Rouphael, Y., Cardarelli, M., Lucini, L., and Colla, G. (2016). Changes in biomass, mineral composition, and quality of cardoon in response to NO-3:Cl- ratio and nitrate deprivation from the nutrient solution. Front. Plant Sci. 7:978. doi: 10.3389/fpls.2016.00978
Colla, G., Hoagland, L., Ruzzi, M., Cardarelli, M., Bonini, P., Canaguier, R., et al. (2017a). Biostimulant action of protein hydrolysates: unraveling their effects on plant physiology and microbiome. Front. Plant Sci. 8:2202. doi: 10.3389/fpls.2017.02202
Colla, G., Cardarelli, M., Bonini, P., and Rouphael, Y. (2017b). Foliar applications of protein hydrolysate, plant and seaweed extracts increase yield but differentially modulate fruit quality of greenhouse tomato. HortScience 52, 1214–1220. doi: 10.21273/hortsci12200-17
Colla, G., Nardi, S., Cardarelli, M., Ertani, A., Lucini, L., Canaguier, R., et al. (2015a). Protein hydrolysates as biostimulants in horticulture. Sci. Hortic. 196, 28–38. doi: 10.1016/j.scienta.2015.08.037
Colla, G., Rouphael, Y., Di Mattia, E., El-Nakhel, C., and Cardarelli, M. (2015b). Co-inoculation of Glomus intraradices and Trichoderma atroviride acts as a biostimulant to promote growth, yield and nutrient uptake of vegetable crops. J. Sci. Food. Agric. 95, 1706–1715. doi: 10.1002/jsfa.6875
Colla, G., Rouphael, Y., Canaguier, R., Svecova, E., and Cardarelli, M. (2014). Biostimulant action of a plant-derived protein hydrolysate produced through enzymatic hydrolysis. Front. Plant Sci. 5:448. doi: 10.3389/fpls.2014.00448
Colla, G., Rouphael, Y., Mirabelli, C., and Cardarelli, M. (2011). Nitrogen-use efficiency traits of mini-watermelon in response to grafting and nitrogen-fertilization doses. J. Plant Nutr. Soil Sci. 174, 933–994. doi: 10.1002/jpln.201000325
Colla, G., Svecova, E., Rouphael, Y., Cardarelli, M., Reynaud, H., Canaguier, R., et al. (2013). Effectiveness of a plant –derived protein hydrolysate to improve crop performances under different growing conditions. Acta Hortic. 1009, 175–179. doi: 10.17660/actahortic.2013.1009.21
Ertani, A., Francioso, O., Tugnoli, V., Righi, V., and Nardi, S. (2011). Effect of commercial lignosulfonate-humate on Zea mays L. metabolism. J. Agric. Food Chem. 59, 11940–11948. doi: 10.1021/jf202473e
Gomez, K. A., and Gomez, A. A. (1983). “Comparison between treatment means,” in Statistical Procedures for Agricultural Research, 2nd Edn, eds K. A. Gomez and A. A. Gomez (Hoboken NJ: John Wiley & Sons), 187–240.
Koevoets, I. T., Venema, J. H., Elzenga, J. T. M., and Testerink, C. (2016). Roots withstanding their environment: exploiting root system architecture responses to abiotic stress to improve crop tolerance. Front. Plant Sci. 7:1335. doi: 10.3389/fpls.2016.01335
Lin, Y., Zhao, Z., Zhou, S., Liu, L., Kong, W., Chen, H., et al. (2017). Top Bending Panicle1 is involved in brassinosteroid signaling and regulates the plant architecture in rice. Plant Physiol. Biochem. 121, 1–13. doi: 10.1016/j.plaphy.2017.10.001
Lucini, L., Borgognone, D., Rouphael, Y., Cardarelli, M., Bernardi, J., and Colla, G. (2016). Mild potassium chloride stress alters the mineral composition, hormone network, and phenolic profile in artichoke leaves. Front. Plant Sci. 7:948. doi: 10.3389/fpls.2016.00948
Lucini, L., Rocchetti, G., Kane, D., and Trevisan, M. (2017). Phenolic fingerprint allows discriminating processed tomato products and tracing different processing sites. Food Control 73, 696–703. doi: 10.1016/j.foodcont.2016.09.020
Lucini, L., Rouphael, Y., Cardarelli, M., Canguier, R., Kumar, P., and Colla, G. (2015). The effect of a plant-derived biostimulant on metabolic profiling and crop performance of lettuce grown under saline conditions. Sci. Hortic. 182, 124–133. doi: 10.1016/j.scienta.2014.11.022
Ntatsi, G., Aliferis, K. A., Rouphael, Y., Napolitano, F., Makris, K., Kalala, G., et al. (2017). Salinity source alters mineral composition and metabolism of Cichorium spinosum. Environ. Exp. Bot. 141, 113–123. doi: 10.1016/j.envexpbot.2017.07.002
Pretali, L., Bernardo, L., Butterfield, T. S., Trevisan, M., and Lucini, L. (2016). Botanical and biological pesticides elicit a similar induced systemic response in tomato (Solanum lycopersicum) secondary metabolism. Phytochemistry 130, 56–63. doi: 10.1016/j.phytochem.2016.04.002
Rombouts, C., Hemeryck, L. Y., Van Hecke, T., De Smet, S., De Vos, W. H., and Vanhaecke, L. (2017). Untargeted metabolomics of colonic digests reveals kynurenine pathway metabolites, dityrosine and 3-dehydroxycarnitine as red versus white meat discriminating metabolites. Sci. Rep. 7:42514. doi: 10.1038/srep42514
Rouphael, Y., Cardarelli, M., Bonini, P., and Colla, G. (2017a). Synergistic action of a microbial-based biostimulant and a plant derived-protein hydrolysate enhances lettuce tolerance to alkalinity and salinity. Front. Plant Sci. 8:131. doi: 10.3389/fpls.2017.00131
Rouphael, Y., Colla, G., Giordano, M., El-Nakhel, C., Kyriacou, M. C., and De Pascale, S. (2017b). Foliar applications of a legume-derived protein hydrolysate elicit dose dependent increases of growth, leaf mineral composition, yield and fruit quality in two greenhouse tomato cultivars. Sci. Hortic. 22, 353–360. doi: 10.1016/j.scienta.2017.09.007
Rouphael, Y., De Micco, V., Arena, C., Raimondi, G., Colla, G., and De Pascale, S. (2017c). Effect of Ecklonia maxima seaweed extract on yield, mineral composition, gas exchange and leaf anatomy of zucchini squash grown under saline conditions. J. Appl. Phycol. 29, 459–470. doi: 10.1007/s10811-016-0937-x
Rouphael, Y., Franken, P., Schneider, C., Schwarz, D., Giovannetti, M., Agnolucci, M., et al. (2015). Arbuscular mycorrhizal fungi act as biostimulants in horticultural crops. Sci. Hortic. 196, 91–108. doi: 10.1016/j.scienta.2015.09.002
Soleimani Aghdam, M., Babalar, M., and Askari Sarcheshmeh, M. A. (2016). “Impact of brassinosteroids on postharvest physiology of fruits and vegetables,” in Eco-Friendly Technology for Postharvest Produce Quality, ed. M. W. Siddiqui (Amsterdam: Elsevier Science), 203–218. doi: 10.1016/B978-0-12-804313-4.00006-2
Tavares, O. C. H., Santos, L. A., Ferreira, L. M., Sperandio, M. V. L., da Rocha, J. G., Garcia, A. C., et al. (2017). You have full text access to this content Humic acid differentially improves nitrate kinetics under low- and high-affinity systems and alters the expression of plasma membrane H+-ATPases and nitrate transporters in rice. Ann. Appl. Biol. 170, 89–103. doi: 10.1111/aab.12317
Keywords: biostimulants, Cucumis melo L., hormone-like activity, lignosulfonates, metabolomics, peptides
Citation: Lucini L, Rouphael Y, Cardarelli M, Bonini P, Baffi C and Colla G (2018) A Vegetal Biopolymer-Based Biostimulant Promoted Root Growth in Melon While Triggering Brassinosteroids and Stress-Related Compounds. Front. Plant Sci. 9:472. doi: 10.3389/fpls.2018.00472
Received: 03 March 2018; Accepted: 26 March 2018;
Published: 10 April 2018.
Edited by:Barbara De Lucia, Università degli Studi di Bari Aldo Moro, Italy
Reviewed by:Georgios Soteriou, Agricultural Research Institute, Cyprus
Vincenzo Candido, University of Basilicata, Italy
Copyright © 2018 Lucini, Rouphael, Cardarelli, Bonini, Baffi and Colla. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Giuseppe Colla, firstname.lastname@example.org