- Department of Plant and Environmental Sciences, Clemson University, Clemson, SC, United States
Unlike quantitative changes, the compositional changes of plant phenolics and changes in their tissue association as influenced by the nutrient supply are less well understood. We evaluated the quantity, composition, and tissue association of phenolics in leaves of two Fragaria ananassa cultivars in response to different levels of nitrogen (N) fertilization using global metabolomic approaches. Influence of N supply on phenolic content in both cultivars was similar, but the magnitude of this response was compound specific. Ellagitannins, the most abundant class of phenolic oligomers, were less responsive to the applied N treatments, whereas proanthocyanidins, the less abundant class of phenolic oligomers, exhibited higher fold change. Within mono-phenolics, the hydroxycinnamates were more abundant but showed lower fold change than the hydroxybenzoates. Among flavonoids, the hydroxylated flavonols showed higher abundances than the flavones, with a preferential accumulation of dihydroxylated flavonol at lower N levels. Furthermore, glycosylated flavonols were higher than the acylated forms. The extractable fraction of phenolics was more influenced by the N treatment than the fiber-bound fraction. The extensive compositional modification of phenolics and a greater response of non-bound fractions in response to N rates highlight the potential to use precise management of N supply as an effective strategy to enhance the bioactive compounds in crops.
Introduction
Secondary metabolites play a myriad of essential roles in plant–organismal and plant–environmental interactions, including functioning as signaling molecules, structural polymers, attractants and deterrents, regulators of phytohormones, and stress mitigators (Dixon et al., 2005). Secondary metabolite production in plants is influenced by the plant genotype, the environment, and their interaction; hence, the composition of secondary metabolites in plants is highly diverse (Wurtzel and Kutchan, 2016). This compositional diversity confers a broader array of biological functions to the secondary metabolites. The molecular identity also influences the association of metabolites with different cellular compartments, which in turn guides the overall function of these compounds in plants. Although the changes in the overall quantity of secondary metabolites in relation to environmental stressors have been widely studied, we lack an in-depth understanding of the stress-induced changes in the composition and tissue association of secondary metabolites. This knowledge could be critical to forecasting the biological role of these compounds in facilitating overall plant fitness.
Phenolic compounds (phenolics) constitute one key class of carbon-based secondary metabolites that are of ecological and nutritive importance. Phenolic metabolites are biosynthesized from the derivatives of the shikimate pathway and are broadly categorized as phenolic acids, flavonoids, and tannins (Vogt, 2010). The phenolic acids and flavonoids together constitute the phenylpropanoids, whereas the tannins encompass two classes, namely, hydrolyzable tannins (HTs) and proanthocyanidins (condensed tannins; Suseela and Tharayil, 2018). The shikimate pathway intermediate viz. dehydroshikimate is utilized in the biosynthesis of gallic acid, the monomer unit for HT (Figure 1). Gallic acid monomers ester-link to glucose form simple galloyl glucoses, which further condense to form gallotannins (GTs) or ellagitannins (ETs; Salminen and Maarit, 2011; Top et al., 2017). The terminal product of the shikimate pathway, chorismic acid, is the precursor of phenylalanine, the precursor of phenylpropanoids. Deamination of phenylalanine generates phenolic acids, hydroxybenzoates (C6–C1), and hydroxycinnamates (C6–C3), which are then utilized for the synthesis of diverse salicinoids and flavonoid compounds. The flavonoids are further categorized into six broad groups based on the chemical structure derived from the enzyme-catalyzed modification of the basic flavan scaffold (C6–C3–C6) (Winkel-Shirley, 2002; Vogt, 2010): flavanones, flavones, flavonols, flavan-3-ols, anthocyanins, and proanthocyanidins (Ferreyra et al., 2012; Alseekh et al., 2020). Diversity within the polymeric phenolic compounds that include the proanthocyanidins and HTs is attributed to the variation in the type, number, and linkage between the monomers (Xie and Dixon, 2005; Salminen and Maarit, 2011). Although previous studies have tracked the changes in content of phenolic compounds as a function of environmental stress, the environment-induced changes in the composition and localization of these metabolites remain less understood.

Figure 1. Scheme of the biosynthesis of phenolics in plants focusing on the flavonoid classes and the potential trade-off between the HTs and flavonoids. Discontinuous lines indicate multiple steps that are not all included in the figure. Phenolic acids and flavonoids together constitute the phenylpropanoids. The key enzymes that regulate this biosynthesis are indicated: PAL, phenylalanine ammonia lyase; CHS, chalcone synthase; CHI, chalcone isomerase; FS, flavone synthase; F3H, flavanone3-hydroxylase; FLS, flavonol synthase; DFR, dihydroflavonol 4-reductase; LDOX, leucoanthocyanidin dioxygenase; LCR, leucoanthocyanidin reductase; ANR, anthocyanidin reductase; AT, acyltransferase; UFGT, UDP flavonoid glucosyltransferase.
The stress mitigation efficiency of phenolic compounds is greatly influenced by their molecular structure (Heim et al., 2002; Alseekh et al., 2020). For example, compared to mono-phenolics, flavonoids have a higher antioxidant capacity due to the enhanced electron-donating properties of the catechol group in the B-ring. Since the radical scavenging activity of flavonoids is affected by the hydrogen donation of its phenolic hydroxyl group, the antioxidant capacity increases with the number of hydroxyl groups in the B-ring (Heim et al., 2002). Along with the number of hydroxyl groups, their position on the B-ring also contributes to the antioxidant potential. For example, compared to an ortho-diphenolic structure of quercitin, the meta-arrangement of dihydroxyl groups reduces the antioxidant potential of the morin by half. Within the flavonoids, the flavonols and flavones have conjugations between the A- and B-ring and thus a more efficient electron delocalization across the molecule. Hence, they have a higher antioxidant capacity than flavanones, which have a saturated C-ring. Furthermore, saturation of the C-ring and the hydroxylation of the A- and B-ring interact to modify the overall antioxidant potential of flavonoids (Heim et al., 2002). Similarly, the glycosylation of flavonoids reduces their antioxidant potential, as it affects the delocalization of electrons but this modification facilitates their transport within plants (Kähkönen and Heinonen, 2003; Cyboran-Mikołajczyk et al., 2018).
Along with the above compositional changes, under stress, the structural changes determine the association (localization) of the phenolics with cell wall. Phenolics exist either as extractable (soluble) forms present mainly in vacuoles or as non-extractable forms that are attached to cell wall proteins and carbohydrates (fiber-bound). The localization of phenolics significantly influences their function, with polyphenols conjugated to cell walls exhibiting a higher photoprotection efficiency than those sequestered in cell vacuoles (Clarke and Robinson, 2008). For example, the flavonoid molecules sequestered in vacuoles serve as antioxidants by quenching reactive oxygen species (ROS), and the same molecule when bound to the cell wall serves as a photoprotectant by absorbing UV-A radiation. The localization of phenolics also affects their transport between sources (e.g., leaf) and sinks (e.g., fruits) (Cadot et al., 2011; Allegro et al., 2016; Shahidi and Yeo, 2016) and thus can influence the nutritive quality of produce. Despite the influence that the chemical composition and localization have on their biological activity, phenolic compounds have less frequently been subjected to finer-level compositional studies.
Nitrogen (N) is a key macronutrient, and its availability dictates the partitioning and transportation of photosynthetic assimilates between the primary and secondary metabolism in plants (Bot et al., 2009; Tegeder and Masclaux-Daubresse, 2018; Fernie et al., 2020). The major theories proposed to explain the variation in the biosynthesis of secondary metabolites in response to N availability include the carbon nutrient balance theory (Coley et al., 1985), the growth–differentiation balance theory (Herms and Mattson, 1992), and the protein competition model (Jones and Hartley, 1999). According to these theories, under N limitation, plants direct the accumulation of carbon-based assimilates toward the synthesis of secondary compounds. Along these lines, many crop species have shown increased phenolic content under N limitation (Fritz et al., 2006; Larbat et al., 2012; Becker et al., 2015). However, the compositional changes of phenolics as influenced by plant ontogeny and the amount and form of N fertilization are less known. Foliar phytochemistry influences plant vigor through stress mitigation (Bobeica et al., 2015), susceptibility to pest and pathogens through deterrence (Kim et al., 2016), and nutritive quality of the yield through source–sink interactions (Fernie et al., 2020). Phytochemicals produced in the leaves are transported to fruits (Petrussa et al., 2013; Gutierrez et al., 2017), making the foliar phytochemistry a key determinant of the content and quality of phytochemicals of the fruits. This transport of the phytochemicals out of leaves is primarily determined by the association of the compounds with the cellular matrix (Top et al., 2017). Compared to phytochemicals that are ester/ether-linked to cell walls (cell wall-bound), compounds that are sequestered in vacuoles are readily transported across the cells and organs. Thus, the association of phytochemicals within the foliage could influence the chemical composition of the fruits. However, little is known about the influence of N limitation on the reprogramming of the tissue association of phenolics.
The metabolome, which is the end product of gene expression, could effectively capture the physiology of stress mitigation in plants (Alseekh and Fernie, 2018). Thus, metabolomics could be an ideal tool to understand plant ecophysiology, which could pave the way for developing new biotechnological solutions benefiting both plants and the environment (Alseekh et al., 2018). The main objectives of this study were to evaluate how different levels of N fertilization affect the content, composition, and association of different classes of foliar phenolics. We hypothesized that (i) the significant reprogramming of phenylpropanoid pathway under N stress would not only result in quantitative changes in the total phenolic content but also impart changes in their composition, whereby the proportional abundance of different classes of phenolics and the glycosylation within a class would be influenced by N; (ii) the altered N supply would lead to the differential association of the biosynthesized phenolics with the tissue fraction (soluble vs. bound), and the proportion of soluble-to-bound phenolics would be higher with increased N limitation; (iii) above changes in composition and localization of phenolic compounds could exhibit a non-linear relationship with N availability.
Materials and Methods
Plant Source
Strawberry plants (Fragaria ananassa) cv. Camarosa (June-bearing) and cv. Albion (day-neutral) were procured as plugs (3–5 leaves) from Lassen Canyon Nursery, Inc., Redding, CA, United States. The plants were grown at Clemson University, SC, United States, in a greenhouse in 7-L plastic pots with ∼2 kg potting mixture, which had a base N level of 7.5 μg/g (150 mg N/pot). The treatment application was initiated 2 months after the plugs were transplanted to the pots. Four levels of N treatment were evaluated in the study, consisting of the control (0 N), 8 mM N, 16 mM N, and 30 mM N. Nitrogen fertilizer was supplied in the form of ammonium nitrate, and the treatment solution was applied as 200 ml/pot at an interval of 10 days. This fertilizer dosage corresponds to an N addition of 0, 43, 86, and 160 mg/pot every 10 days. The treatment time and rate of application were selected based on the preliminary experiments on N addition, with the 30-mM N treatment being closer to the recommend rate of N application for strawberries in the Southeastern United States under annual plasticulture production systems (Rysin et al., 2015). The treatment consisted of a completely randomized design with five biological replicates. Twelve weeks after treatment, mature actively growing young leaves were harvested from the third whorl of each plant and immediately frozen using dry ice. The leaves of similar leaf area were selected across the treatment, and at the time of harvest, none of the plants, including the control treatment, exhibited visible chlorosis. The tissues were ground into a fine powder with dry ice and immediately transferred to −80°C for storage until further analysis.
Estimation of Primary Metabolites by Gas Chromatography–Mass Spectrometry
Polar metabolites were extracted from the ground tissue as per Suseela et al. (2015), with slight modification. Briefly, 300 mg of powdered plant material was extracted with 3 ml of methanol using a bead-beater homogenizer (ceramic beads, 6,000 rpm, 2 min) and sonication (25% amplitude, 1 min) and centrifuged at 1,500 × g. The samples were kept at 4°C during the processing. One milliliter of supernatant was transferred to a glass tube, and the extract was further partitioned to polar and non-polar phase by successive addition and mixing 1 ml of ice-cold CHCl3 and water. The tubes were centrifuged at 1,500 × g to complete the phase separation between the lower CHCl3 phase and the top aqueous methanol phase. Forty microliters of the aqueous methanolic supernatant was transferred to silanized glass inserts with 20 μl of 200 μg/ml ribitol (internal standard) and 500 μg/ml d27-myristic acid (retention time lock). The samples and standards were completely dried in a vacuum evaporator at room temperature. Twenty microliters of a 40-mg/ml solution of methoxylamine hydrochloride in pyridine was added to the glass insets and incubated at 40°C for 60 min with intermittent mixing. The second derivatization consisted of the addition of 90 μl of N-methyl-N-(trimethylsilyl) trifluoroacetamide (MSTFA) with 1% trimethylchlorosilane (TMCS) and incubating the samples for 30 min at 55°C. The derivatized samples were analyzed on a gas chromatograph coupled to a mass spectrometer (GC-MS; Agilent 7980, Agilent Technologies, Santa Clara, CA, United States). Compounds were separated on a DB-5 column. The GC-MS instrument settings were as per Maroli et al. (2015). Mass spectra were processed with an automatic mass spectral deconvolution and identification system (AMDIS v2.71, NIST, Gaithersburg, MD, United States) with a minimum match factor of 70. The identification of the compounds was based on their retention time index and mass and fragmentation pattern comparison with the Fiehn library (G1676AA; Agilent Technologies, Wilmington, DE, United States; Suseela et al., 2015). The peak intensities were normalized to ribitol.
Estimation of Proanthocyanidins by Acid Butanol Assay
The total proanthocyanidins in the leaf sample were estimated by acid butanol assay modified from Top et al. (2017). The aqueous methanol phase was dried at room temperature in a vacuum evaporator, followed by the addition of 1 ml of acid butanol reagent and incubation at 95°C for 60 min. A spectrophotometric measurement of the color that developed was taken at 550 nm, and the amount of proanthocyanidins was calculated using a standard curve derived from the cyanidin standard.
Estimation of Extractable Phenolics by Liquid Chromatography–Mass Spectrometry
Briefly, 600 mg of ground leaf tissue was extracted with 4 ml of methanol and homogenized at 4°C as described above and centrifuged 3,000 rpm. The sediment was then sequentially extracted twice with methanol and thrice with 80% aqueous acetone to completely extract the vacuolar metabolites. The supernatant from all six extractions was pooled and was used for the determination of soluble (extractable) phenolics. Isotope-labeled resveratrol (13C6; 0.5 ng/μl) was added as the internal standard. The samples were analyzed on an Ultimate 3000 HPLC (ultrahigh-pressure liquid chromatography) coupled to an Orbitrap Fusion Tribrid mass spectrometer equipped with electrospray ionization (ESI; Thermo Scientific, Waltham, MA, United States) as per Bowers et al. (2018). Separation of the metabolites was carried out using an HSS T3 column (Waters Corp., Milford, MA, United States; 150 mm × 2.1 mm, 1.8 μm) at 30°C. The following gradient program utilizing water with 0.1% formic acid as mobile phase A and acetonitrile as mobile phase B was employed: 0 min, 5% B; 2 min, 5% B; 24 min, 60% B; 27 min, 90% B; followed by a 3-min washing step at 90% B and a subsequent re-equilibration for 6 min at 5% B. The flow rate was set to 0.22 ml/min, and the injection volume was set as 2 μl. The mass spectrometer was operated as per Bowers et al. (2018) in negative ionization mode with a data-dependent MS2 HCD-CID method for compound confirmation. The interface conditions were as follows: emitter voltage, −2,600 V; vaporizer temperature, 325°C; ion transfer tube, 325°C; sheath gas, 55 (arb); aux gas, 10 (arb); and sweep gas, 1 (arb).
Estimation of Bound Phenolics by Liquid Chromatography–Mass Spectrometry
The bound (non-extractable) phenolics were determined on the residue left after the sequential extraction with methanol and acetone. The residue was dried at 60°C and acid hydrolyzed using 2 M methanolic HCl with 0.04% ascorbic acid at 85°C for 4 h. The supernatant after centrifugation was used for the determination of the bound phenolics (Top et al., 2017). The samples were analyzed using an Orbitrap Fusion Tribrid mass spectrometer as described above.
Determination of Antioxidant Capacity
The extract from the extractable phenolics for liquid chromatography–mass spectrometry (LC-MS) analysis was used for the determination of the antioxidant capacity of the strawberry leaves from the different N treatments using a DPPH assay (Xie and Schaich, 2014). Briefly, 100 μl of the extract was added to 900 μl of freshly prepared 1 mM methanolic DPPH reagent and incubated at room temperature in the dark. After 30 min, the absorbance was measured at 515 nm against the reagent blank. The percentage of DPPH radical scavenging was calculated as given in Eq. 1, where A0 is the absorbance of the reagent without the sample, and A is the absorbance of the reaction mixture with the sample solution.
Data Analysis and Statistics
Statistical significance was determined by one-way ANOVA followed by Tukey’s multiple comparisons for normally distributed data and by Kruskal–Wallis ANOVA for non-parametric data using JMP Pro13.1.0 (SAS Institute Inc., Cary, NC, United States) with a statistical significance level of P < 0.05. To evaluate the effect of the applied N treatments on the overall foliar metabolite profile, the metabolomics data were subjected to multivariate analysis using Metaboanalyst 4.0 (Chong et al., 2018). The metabolomic data of the two cultivars were autoscaled or log-transformed to satisfy the assumptions of normality. The metabolites were further analyzed with partial least squares discriminant analysis (PLS-DA) along with built-in cross-validation and permutation testing. The tentative identification of the phenolic metabolites was based on the literature reporting the accurate masses of the precursor ion and fragment ions using Xcalibur 4.1 (Thermo Fisher Scientific, MA, United States). The data were processed by the node-based workflow of the Compound Discoverer 2.2 software (Thermo Fisher Scientific, MA, United States). The 13C resveratrol normalized peak area was utilized to evaluate the phenolic metabolites quantitatively. Graphs were plotted using GraphPad Prism 7.05 (GraphPad software, La Jolla, CA, United States).
Results
Variation Among the Primary Metabolites as a Response to the Nitrogen Supply
The composition of primary metabolites in strawberry leaves, which included sugars, sugar alcohols, organic acids, and amino acids, was quantified using a GC-MS-based metabolomic approach. The analysis identified 40 and 49 metabolites from the leaf extract of cv. Albion and cv. Camarosa, respectively. PLS-DA of these metabolites along the first two-component axis (Figure 2) explained ca. 73–78% of the total variation in the data. The PLS-DA analysis showed good model fitness with R2 and Q2 value of >0.6. The robustness of the class discrimination was verified through permutation testing of separation distance based on the ratio of the between-group sum of squares and the within-group sum of squares and had P < 0.05 over 2,000 iterations. The metabolites with variable importance in projection (VIP) scores >1 are considered relevant for influencing the grouping (Supplementary Figure 1). In both cultivars, the component axis 1 separated control treatment from the other N treatments, and all other N treatments grouped into distinct clusters along the component axis 2. Although the metabolites possessing the highest discrimination potential varied between the two cultivars, they predominantly constituted tagatose, trehalose, citric acid, fumaric acid, and galactitol. All metabolites that discriminated between the treatments did not exhibit a linear response to the level of N additions.

Figure 2. Partial least squares discriminant analysis (PLS-DA) score plot of the first two principal components of the analysis of polar primary metabolites of strawberry cv. Albion (A) and cv. Camarosa (B) exposed to four levels of nitrogen (Green: 30 mM N, Red: 16 mM N, Blue: 8 mM N, and Turquoise: Control). The permutation cross-validation of the PLS-DA model had P = 0.0005 over 2,000 iterations. The ellipse represents a 95% confidence interval, and data points represent biological replicates (n = 5).
The relative abundances of the total non-structural sugars (Figures 3A,B) exhibited the least variation between N-replete (30 mM N) and -deficient (control) treatment for both cultivars. Sucrose was the most abundant sugar in both cultivars, contributing to ca. 69–80% of the total sugar and did not vary with N treatment. Other soluble sugars like glucose, trehalose, and tagatose exhibited non-linear responses (Supplementary Figure 1) with accumulation at N-limited treatment for the two cultivars. The tricarboxylic acid (TCA) cycle intermediates showed differential response across the treatments in the two cultivars (Figure 3C). Dehydroascorbate and quinate constituted the non-TCA organic acids and exhibited much lower percentage variation across the N treatments.

Figure 3. Abundance of sugars in cv. Albion (A) and cv. Camarosa (B) and intermediaries of tricarboxylic acid (TCA) cycle (C) of cv. Albion (open bars) and cv. Camarosa (closed bars), as influenced by the four levels of N treatments. The data represent the mean ± SD (n = 5) ribitol normalized abundance of the respective samples and were tested by one-way ANOVA with Tukey’s honestly significant difference (HSD). Bars with the same letters are not significantly different at the 95% confidence interval.
Variation in the Content and Composition of Polyphenols in Response to Varied Nitrogen Supply
After accounting for multiple adducts, potential source fragmentation, UHPLC-ESI-Orbitrap-MS data analysis of the strawberry leaf tissue using Compound Discoverer 2.2, positively annotated 789 potential compounds across all samples. Among them, 738 and 680 compounds in cv. Albion and cv. Camarosa, respectively, were classified as significant based on the log fold change (2), P value, and false discovery rate (<0.05). Score plot from PLS-DA of the 789 metabolites plotted for two discriminatory components explained ca. 74–83% of the total variation and separated the different N treatments for both cultivars (Figure 4). The PLS-DA model showed a good model fit with R2 and Q2 value >0.6 and P value <0.01 based on the permutation testing. The tentative molecular identity of 149 metabolites (Supplementary Table 1) was established based on manual curation based on their accurate molecular masses (<5 ppm), MS2 fragmentation pattern from the literature, and online mass spectral libraries (mzCloud, Massbank). The identified phenolic compounds belonged to non-anthocyanin categories of phenolics, which included simple galloyl glucoses, ellagic acid derivatives, ETs, hydroxybenzoates, hydroxycinnamates, flavanones, flavones, flavonols, flavan-3-ols, and proanthocyanidins. Compounds having the same molecular masses (mass error <2 ppm of theoretical mass) but having different retention times within the same sample (>0.5 min) were treated as a unique compound. The cumulative peak area of the identified metabolites constituted ca. 54% of the total peak area of the 789 compounds. The PLS-DA of the identified 149 metabolites also differentiated the N treatments along the principal component axis 1 that explained >61% variation, with the control treatments significantly different from the N-addition treatments (Supplementary Figure 2).

Figure 4. Partial least squares discriminant analysis (PLS-DA) score plot of the first two principal components of the analysis of 789 secondary metabolites of strawberry cv. Albion (A) and cv. Camarosa (B) exposed to four levels of nitrogen (Green: 30 mM N, Red: 16 mM N, Blue: 8 mM N, and Turquoise: Control). The permutation cross-validation of the PLS-DA model had P = 0.0005 over 2,000 iterations. The ellipse represents a 95% confidence interval, and data points represent biological replicates (n = 5).
Ellagitannins were the most abundant compound class in the leaf metabolites, comprising 48 out of the 149 tentatively identified phenolic compounds and accounts for 23% of the total peak area of the 789 metabolites. For both cultivars, the total ET content increased (P = 0.01) with a decrease in N application (Figure 5A). For both cultivars, the content of hexahydroxydiphenoyl (HHDP) ET and dehydro-HHDP esters (geraniin) doubled under lower N treatments (control and 8 mM N) as against N-replete treatment (Figures 5B,C). Galloyl glucoses and tetragalloyl glucoses, the precursors of ETs, were higher in control treatment for both cultivars (Supplementary Figure 3A). Simple ellagic acid derivatives also showed ca. 1.5-fold increase (P < 0.002) in control treatment as compared to N-replete treatment for both cultivars (Supplementary Figure 3B). Catechin was the most abundant metabolite in both cultivars, accounting for ca. 7% of the total peak area of the 789 metabolites. Oligomeric proanthocyanidins and its precursor, flavan-3-ols, increased with a decrease in N availability in both cultivars (Figures 6A,B). In both cultivars, the total flavan-3-ols and oligomeric proanthocyanidins increased linearly with an increase in N addition. The dimers constituted 68–71% of the total proanthocyanidins identified in both cultivars; however, in cv. Camarosa, the trimer and tetramers showed a higher fold increase under N limitation as compared to the dimers (Figures 6C,D). In general, the hydroxycinnamates were 4- to 10-fold higher than hydroxybenzoates. In both cultivars, the hydroxybenzoates increased with decrease in N supply (Supplementary Figure 4).

Figure 5. Abundances of total ellagitannin (A) and different subgroups of ellagitannins in Fragaria ananassa cv. Albion (B) and cv. Camarosa (C) exposed to four N treatments. The data represent the mean ± SD (n = 5) 13C resveratrol (1 ng) normalized abundance and were tested by one-way ANOVA with Tukey’s honestly significant difference (HSD). Bars with the same letters are not significantly different at the 95% confidence interval.

Figure 6. Abundances of flavan-3-ols (A) and proanthocyanidins (B) and different subgroups of proanthocyanidins in Fragaria ananassa cv. Albion (C), and cv. Camarosa (D) exposed to four N treatments. The data represent the mean ± SD (n = 5) 13C resveratrol (1 ng) normalized abundance and were tested by one-way ANOVA with Tukey’s honestly significant difference (HSD). Bars with the same letters are not significantly different at the 95% confidence interval.
Naringenin and eriodictyol glucoside that constituted the flavanone class that is derived from the chalcones showed ca. 60% increase (P < 0.05) in control as compared to N-replete treatment in both cultivars (Supplementary Figure 4). Flavonols represented the second abundant class of phenolic compounds in strawberry leaves. The flavonol content increased (P < 0.005) in both cultivars, with a decrease in N concentration from 16 mM N to control (Figure 7A). Free aglycones of flavonoids in the soluble fraction were not detected, and glycosides of quercetin and kaempferol were the two key flavonols identified in both cultivars. The kaempferol:quercetin ratio decreased by more than half under N-replete treatment in both cultivars at lower N (control and 8 mM) (Figure 7B). About 61–80% of the total tentatively identified flavonoids were glycosylated and 26–38% acylated (Figures 7C,D). Hexoses and pentoses were the only sugars identified to be associated with flavonols along with glucuronide. Only coumaroyl hexoside was an acylated flavonol form identified in both cultivars and constituted about 30% of the total identified flavonols. Coumaroyl quercetin increased (P < 0.02) by >5-fold in control treatment of both cultivars as compared to 30 mM N treatment, whereas coumaroyl kaempferol hexoside showed non-significant variation across the N treatments.

Figure 7. Abundances of flavonols (A) and kaempferol:quercetin ratio (B) and different structural modifications of kaempferol and quercetin in Fragaria ananassa cv. Albion (C) and cv. Camarosa (D) exposed to four N treatments. The data represent the mean ± SD (n = 5) 13C resveratrol (1 ng) normalized abundance and were tested by one-way ANOVA with Tukey’s honestly significant difference (HSD). Bars with the same letters are not significantly different at the 95% confidence interval.
Variation in the Localization of the Polyphenols in Response to Varied Nitrogen Supply
Based on accurate mass and fragmentation patterns of phenolics from literature, 17 phenolic aglycone metabolites (Supplementary Table 1) could be tentatively identified from the acid hydrolysis of residue. Total bound phenolics did not show any significant variation to the applied N treatments. Thus, change in the soluble phenolics was the major contributor to the observed change in the proportion of bound to soluble phenolics. The total soluble phenolics were >1.5-fold higher under N limitation compared to N-replete in both cultivars (Figure 8).

Figure 8. Proportion of soluble vs. fiber-bound phenolics in Fragaria ananassa cv. Albion (A) and cv. Camarosa (B) exposed to four N treatments. The data represent the mean ± SD (n = 5). Bars with the same letters are not significantly different at the 95% confidence interval. “*” Denotes statistical significance (P < 0.05) between bound and soluble phenolics with a treatment.
Derivatives of protocatechuic acid, coumaric acid, gallic acid, caffeic acid, and isoferulic acid were used for computing the variation between bound and soluble simple phenolic acids. In cv. Albion, bound and soluble simple phenolic acids increased by >50% (P < 0.01) in control treatment as compared to N-replete treatments (Supplementary Figure 5A). In cv. Camarosa, the bound-to-soluble ratio decreased by half with a 1.2-fold increase (P < 0.001) in the soluble phenolic acids in the control treatment (Supplementary Figure 5B). For both cultivars, the soluble flavonols showed >50% increase (P < 0.001) with a decrease in N (Supplementary Figures 5C,D). The sum of ellagic and gallic acid content, the hydrolysis product of bound HTs, did not show any variation to the applied N treatment, and total soluble HT increased by 35–47% and ca. 56–91% in cv. Albion and cv. Camarosa, respectively, with the decrease in N from 16 mM N to control (Supplementary Figures 6A,B). In cv. Albion, the total proanthocyanidin content was higher [2.2 mg cyanidin equivalent/g tissue fresh weight (FW)] in N-deficient plants as compared to cv. Camarosa (1.1 mg cyanidin equivalent/g tissue FW). In cv. Albion, 45% of the total proanthocyanidins were in an extractable form in the control treatment, which dropped (P > 0.05) to 33% under N-replete levels (30 mM N) (Supplementary Figure 6C). However, in cv. Camarosa, under N limitation, 70% of the total proanthocyanidins were in extractable form, which significantly dropped to 35% under N-replete levels (30 mM N) (Supplementary Figure 6D). The N-deficient plants showed approximately 90 and 88% radical scavenging of 0.1 mM DPPH solution for cv. Camarosa and cv. Albion, respectively (Supplementary Figure 10), whereas the radical scavenging activity decreased to 48 and 59% for cv. Camarosa and cv. Albion, respectively, under N-replete conditions.
Discussion
Variation in the Content of Primary Metabolites in Response to the Nitrogen Supply
By regulating the C/N balance, N determines the feedback mechanisms of the biochemical reactions (Paul and Pellny, 2003). As observed on this study, higher N can drive higher tissue carbohydrate content, similar to the positive curvilinear relationship between the soluble sugars and N levels reported in potato leaves, grapes, and corn (Chen and Cheng, 2003; Jin et al., 2015; Braun et al., 2016). This response is thought to stem from a plant strategy of maintaining a balance between photochemistry and physiology, where the production of sugars acts as a non-radiative energy dissipation mechanism (Reyes et al., 2016). In addition to sucrose, which was the most abundant sugar in leaves, other soluble sugars, and sugar alcohols exhibited a much lower and non-linear response across the N treatments (Supplementary Figure 1). These metabolites are known for their role as osmoprotectants, scavengers of ROS, and signaling molecules and are reported to be upregulated under cold, heat, and drought stress (Nishizawa et al., 2008; Koehler et al., 2015). Organic acids derived from the anaplerotic TCA cycle act as transient carbon substrates for amino acid synthesis and as regulators of redox and energy levels (Igamberdiev and Eprintsev, 2016). In the present study, the non-linear variations of the organic acids (Figure 3C) point toward the multifaceted role of the TCA cycle in generating reducing equivalents and α-ketoglutarate for efficient N assimilation (Igamberdiev and Eprintsev, 2016). The abundance of α-ketoglutarate and glutamate under reduced N conditions is in accordance with previous reports in maize and barley (Amiour et al., 2012; Quan et al., 2016).
Variation in the Content of Phenolic Classes in Response to Varied Nitrogen Supply
The present study captures the changes in the content of phenolics in strawberry leaves under different N levels; however, the magnitude of this change varied among the phenolic classes. The shikimate pathway, the entry point of carbohydrates for the biosynthesis of phenylpropanoids and HTs, represents the first branch point of phenolic biosynthesis. ETs are the major subclass of HT reported in strawberry leaves and consist of C–C-linked gallic acid monomers. Despite being the most abundant class of phenolics in strawberry leaves, the ET content showed a lower response to the variation in N treatments than the phenylpropanoids, which exhibited an increase in their content as the N supply decreased (Figure 5A and Supplementary Figures 6A,B). Weaker induction of ET biosynthesis could arise due to substrate competition for the synthesis of gallic acid and chorismic acid, the precursors for phenylpropanoids. Thus, regulation by N at the branch point of 3-dehydroshikimate is proposed to explain the variation in the abundance pattern of ETs and phenylpropanoids (Salminen and Maarit, 2011). Similar tradeoffs between gallic acid and phenylpropanoids have been reported in strawberry leaves (Kim et al., 2016). These observations highlight that the stress stimuli can be utilized to regulate carbon partitioning between the shikimate and the phenylpropanoid pathways, and N appears to be a potential candidate for imposing this modulation.
The second branch point of phenolic biosynthesis occurs at coumarate in the phenylpropanoid pathway, resulting in the biosynthesis of the derivatives of simple phenolic acids and flavonoids. Among simple phenolic acids, hydroxycinnamates conjugated with flavonols exhibit a higher free radical scavenging capacity than free hydroxycinnamates (Martinez et al., 2016). Thus, the similar ratio of hydroxycinnamates:flavonoids (Supplementary Figure 9) across the N treatments observed in the present study could represent the effective regulation of carbon partitioning for the biosynthesis of substrates, which in turn could be used to generate conjugated phenolics that are more reactive. The abundance of different flavonoid classes, including flavanones, flavones, flavonols, flavan-3-ols, and proanthocyanidins, was inversely related to N supply, which is in agreement with previous research across different species (Strissel et al., 2008; Becker et al., 2015). The reduced N treatment of 8 mM resulted in a maximum increase in the flavonol content over that of the flavones in both cultivars. Flavonols, because of their 3’-hydroxyl group and lower redox potential, are more potent free radical scavengers than flavones (Dwivedi et al., 2016).
The synthesis of flavan-3-ols and their oligomers (proanthocyanidins/condensed tannins) shares common steps in the flavonoid pathway with the synthesis of anthocyanins and flavonols. In the present study, both flavonol and flavan-3-ols increased with the decrease in N supply, and the anthocyanin abundance was negligible. However, the flavan-3-ols exhibited a much higher fold change with the decrease in N than the flavonols. Dihydroflavonols are utilized in the biosynthesis of not only flavonols but also leucoanthocyanidins, which are utilized for the biosynthesis of anthocyanins and flavan-3-ols. The anthocyanin content observed in our study was lower than the content previously reported under cold treatment in leaves and fruits (Koehler et al., 2015; Tanou et al., 2017). The reduction of anthocyanins to flavan-3-ols by anthocyanin reductases (ANRs) has previously been reported to contribute a major share of the flavan-3-ols in tobacco (Luo et al., 2016), cotton (Zhu et al., 2015), white clover (Abeynayake et al., 2012), grape (Bogs et al., 2007), and strawberry (Fischer et al., 2014). Thus, the observed lower abundance of anthocyanin and higher abundance of flavan-3-ols under N limitation could be attributed to upregulation of the ANR. A similar decrease in anthocyanin content and rechanneling of the metabolic flux to flavan-3-ols was reported in strawberry fruits when UDP-glucose:anthocyanidin glucosyl transferase (FaGT1) was downregulated (Griesser et al., 2008). However, more investigation is required to establish the relationship between N and the activity of the enzymes involved in biosynthesis.
Variation in the Composition of Phenolic Classes in Response to Varied Nitrogen Supply
The simple galloyl glucoses, which are the first branch-off products from the shikimate pathway and the precursors of HT, were higher at lower N (control and 8 mM N) levels. Glucogallin constituted a major (>70%) galloyl glucose ester in the leaves of both strawberry cultivars across all N treatments. A decrease in glucogallin followed by an increase in digalloylglucose has been reported in mature birch leaves and is proposed to be a result of the efficient utilization of glucogallin (Salminen et al., 2001). The higher content of glucogallin in proportion with tri- and tetra-galloyl glucose across all N treatments in both cultivars could indicate reduced kinetics of the galloyltransferases compared to glucosyltransferase in the applied N treatment. Across plant species, it has been reported that GT and ET are not simultaneously produced (Salminen and Maarit, 2011) and that strawberry leaves have a higher abundance of ET compared to polymeric GTs (Kim et al., 2016; Kårlund et al., 2016; Michalska et al., 2017). The preferential biosynthesis of ET in strawberry leaves over GT could be due to preferential oxidation of the branch point substrate pentagalloylglucose, as opposed to the formation of a depside bond between the gallic acid subunits.
The simple phenolic acids from the study were categorized as hydroxybenzoates and hydroxycinnamates. The presence of a p-hydroxy group and an additional carbon side chain structure on the phenol ring has been shown to increase radical scavenging activity (Natella et al., 1999). The observed higher content of hydroxycinnamates (coumaric and ferulic acids) compared to the hydroxybenzoates (protocatechuic acid) under N limitation could represent a structure-dependent preferential accumulation of more effective antioxidants. In addition, the preferential utilization of hydroxycinnamates in conjugation to flavonols could explain the increased biosynthesis of hydroxycinnamates over hydroxybenzoates.
Further along the flavonoid pathway, the flavone apigenin along with the flavonols quercetin and kaempferol (via dihydroflavonols) constituted the most abundant compounds derived from the two branches of the flavanone naringenin in both cultivars across all N treatments. Apigenin is often reported as the only flavone in strawberry fruits (Akhatou et al., 2017) and has been shown to have minimal antioxidant capacity compared to kaempferol and quercetin (Pannala et al., 2001; Sroka et al., 2015). In the present study, the decrease in the kaempferol:quercetin ratio with decreases in the N treatment was attributed to an increase in quercetin and not to a decrease in kaempferol. Regardless of the stress factor, quercetin is reported as a more abundant flavonol in strawberry fruits and leaves compared to kaempferol (Kårlund et al., 2016; Michalska et al., 2017). The dihydroxylated flavonol quercetin has higher superoxide radical scavenging activity than the monohydroxylated kaempferol (Marković et al., 2014; Zhang et al., 2015). The variation in the relative proportions of kaempferol and quercetin could be due to their differential stability under different N conditions. Previous studies in Arabidopsis have reported that glycosides of kaempferol are more stable than quercetin following 2 days of N starvation and high temperatures (Olsen et al., 2009). In accordance with previously reported studies, hexoses, pentoses, and glucuronides were the key glycosylated forms in the strawberry leaves (Supplementary Table 1; Hanhineva et al., 2008; Kim et al., 2016). Glycosylation and acylation not only improve the biological activity of flavonoids but also modify the physical parameters of solubility, thermal/light stability, and membrane permeability (Ogo et al., 2015). Furthermore, differences in the bioactivity of these flavonoids associated with the decorative modification on the flavan structure is well documented (Alseekh et al., 2020), the exact role of individual metabolites under specific stress cue needs more detailed evaluation.
In our study, catechin monomers and proanthocyanidin dimers had higher abundances under N limitation than the trimers and tetramers. The enzymatic (via polyphenol oxidases) and non-enzymatic (chemical oxidation) mechanisms of proanthocyanidin polymerization from monomers are widely debated (Xie and Dixon, 2005), making it difficult to ascertain the effect of N fertilization. A positive correlation has been demonstrated between the mean degree of polymerization and the antioxidant capacity (Zhou et al., 2014). Thus, despite the higher abundances of dimers across the treatments in both cultivars, the higher percentage increase in trimers with the reduction in N supply could be construed as the production of compounds with increased antioxidant capacity. The increase in radical scavenging activity with the increase in the polymerization from n = 1 to 4 and the reversal of that trend at higher degrees of polymerization was attributed to the steric hindrances contributing to the accessibility of hydroxyl groups (de Gaulejac et al., 1999). These findings support the close association of structural conformation with biological activity and require further studies to determine the relationship between N limitation and proanthocyanidins at a higher degree of polymerization.
Variation in the Localization of Phenolic Classes in Response to Varied Nitrogen Supply
In the present study, acid hydrolysis of the residue left after the extraction of the soluble phenolics was performed for the determination of the bound form of the phenolic classes. Although the extractable form of the different phenolic classes increased with a decrease in N supply, the bound form of phenolics was similar across the different N treatments in both cultivars. A similar increase in soluble phenolics with an increase in N and relative stability in bound phenolics was reported in wheat (Stumpf et al., 2015). It has been proposed that genetic determinants primarily regulate the bound phenolics, and the soluble phenolics are primarily regulated by environmental cues (Fernandez-Orozco et al., 2010). Although the exact sequestration mechanism of phenolics onto the cell wall is currently unknown, it is proposed that the phenolic precursors are mobilized from the site of biosynthesis by vesicle trafficking or transporter proteins and polymerized on the structural components of the cell wall usually through peroxidase activity (Zhao et al., 2010). These cell wall-bound phenolics, characterized mainly by phenolic acids, act as a barrier for effectively absorbing harmful UV radiation and limiting water loss by maintaining the integrity and hydrophobicity of the cell wall (Hura et al., 2017). In accordance with the above theory, in the present study, there is a higher proportion of the simple phenolic acids that are bound to the cell wall than flavonols. In addition, an increase in porosity in the cell wall is proposed to facilitate tannin adsorption to the cell walls and decrease tannin extractability (Bindon et al., 2014). Thus, the observed higher proportion of phenolics in the soluble form could be because of the increased porosity and reduced stability of the cell membrane under N limitation.
Conclusion
The structural diversity within the phenolics differentially regulates their biological functions. Our results show that, along with the variation in the quantity of different phenolic classes, N supply influences the compositional variation both within and across different phenolic classes and their association with cellular components. Overall, the phytochemical content increased with decreasing N availability. However, the proportional increase in the quantity of proanthocyanidins compared to the other phenolic classes highlights the non-linear response of phytochemicals to the applied N treatments. The increased preferential accumulation of metabolites with selective structural modifications with enhanced antioxidant capacity, such as hydroxylation, glycosylation, and ring conjugation, under a low N supply, suggests a potential adaptation of these plants to abiotic stress. The specific influence of N supply on the transportable phenolics with higher antioxidant capacity, but not on the bound phenolics, suggests the potential for precise management of moderate nutrient stress as a means to increase the phytochemical content and nutritive value of food crops.
Data Availability Statement
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author/s.
Author Contributions
NT designed the experiment. AN conducted the experiment and analyzed the data. Both authors contributed to the manuscript preparation.
Funding
This project was supported by Agriculture and Food Research Initiative Grant no. 2016-67014-24578 from the USDA National Institute of Food and Agriculture. This material is based upon work supported by the NIFA/USDA, under Project Number SC-1700597. Technical Contribution No. 6958 of the Clemson University Experiment Station.
Conflict of Interest
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.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2021.613839/full#supplementary-material
References
Abeynayake, S. W., Panter, S., Chapman, R., Webster, T., Rochfort, S., Mouradov, A., et al. (2012). Biosynthesis of proanthocyanidins in white clover flowers: cross talk within the flavonoid pathway. Plant Physiol. 158, 666–678. doi: 10.1104/pp.111.189258
Akhatou, I., Sayago, A., González-Domínguez, R., and Fernández-Recamales, A. (2017). Application of targeted metabolomics to investigate optimum growing conditions to enhance bioactive content of strawberry. J. Agric. Food Chem. 65, 9559–9567. doi: 10.1021/acs.jafc.7b03701
Allegro, G., Pastore, C., Valentini, G., Muzzi, E., and Filippetti, I. (2016). Influence of berry ripeness on accumulation, composition and extractability of skin and seed flavonoids in cv. sangiovese (Vitis vinifera L.). J. Sci. Food Agric. 96, 4553–4559. doi: 10.1002/jsfa.7672
Alseekh, S., Bermudez, L., de Haro, L. A., Fernie, A. R., and Carrari, F. (2018). Crop metabolomics: from diagnostics to assisted breeding. Metabolomics 14, 1–13.
Alseekh, S., and Fernie, A. R. (2018). Metabolomics 20 years on: what have we learned and what hurdles remain? Plant J. 94, 933–942. doi: 10.1111/tpj.13950
Alseekh, S., Perez de Souza, L., Benina, M., and Fernie, A. R. (2020). The style and substance of plant flavonoid decoration; towards defining both structure and function. Phytochemistry 174:112347. doi: 10.1016/j.phytochem.2020.112347
Amiour, N., Imbaud, S., Clément, G., Agier, N., Zivy, M., Valot, B., et al. (2012). The use of metabolomics integrated with transcriptomic and proteomic studies for identifying key steps involved in the control of nitrogen metabolism in crops such as maize. J. Exp. Bot. 63, 5017–5033. doi: 10.1093/jxb/ers186
Becker, C., Urlić, B., Špika, M. J., Kläring, H., Krumbein, A., Baldermann, S., et al. (2015). Nitrogen limited red and green leaf lettuce accumulate flavonoid glycosides, caffeic acid derivatives, and sucrose while losing chlorophylls, β-carotene and xanthophylls. PLoS One 10:e0142867. doi: 10.1371/journal.pone.0142867
Bindon, K. A., Madani, S. H., Pendleton, P., Smith, P. A., and Kennedy, J. A. (2014). Factors affecting skin tannin extractability in ripening grapes. J. Agric. Food Chem. 62, 1130–1141. doi: 10.1021/jf4050606
Bobeica, N., Poni, S., Hilbert, G., Renaud, C., Gomès, E., Delrot, S., et al. (2015). Differential responses of sugar, organic acids and anthocyanins to source-sink modulation in cabernet sauvignon and sangiovese grapevines. Front. Plant Sci. 6:382.
Bogs, J., Jaffe, F. W., Takos, A. M., Walker, A. R., and Robinson, S. P. (2007). The grapevine transcription factor VvMYBPA1 regulates proanthocyanidin synthesis during fruit development. Plant Physiol. 143, 1347–1361. doi: 10.1104/pp.106.093203
Bot, J. L., Bénard, C., Robin, C., Bourgaud, F., and Adamowicz, S. (2009). The ‘trade-off’ between synthesis of primary and secondary compounds in young tomato leaves is altered by nitrate nutrition: experimental evidence and model consistency. J. Exper. Bot. 60, 4301–4314. doi: 10.1093/jxb/erp271
Bowers, J. J., Gunawardena, H. P., Cornu, A., Narvekar, A. S., Richieu, A., Deffieux, D., et al. (2018). Rapid screening of ellagitannins in natural sources via targeted reporter ion triggered tandem mass spectrometry. Sci. Rep. 8:10399.
Braun, H., Fontes, P. C. R., Silva, T. P. D., Finger, F. L., Cecon, P. R., and Ferreira, A. P. S. (2016). Carbohydrates concentration in leaves of potato plants affected by nitrogen fertilization rates. Rev. Ceres 63, 241–248. doi: 10.1590/0034-737x201663020016
Cadot, Y., Chevalier, M., and Barbeau, G. (2011). Evolution of the localisation and composition of phenolics in grape skin between veraison and maturity in relation to water availability and some climatic conditions. J. Sci. Food Agric. 91, 1963–1976. doi: 10.1002/jsfa.4401
Chen, L., and Cheng, L. (2003). Carbon assimilation and carbohydrate metabolism of Concord’grape (Vitis labrusca L.) leaves in response to nitrogen supply. J. Am. Soc. Hortic. Sci. 128, 754–760. doi: 10.21273/jashs.128.5.0754
Chong, J., Soufan, O., Li, C., Caraus, I., Li, S., Bourque, G., et al. (2018). MetaboAnalyst 4.0: towards more transparent and integrative metabolomics analysis. Nucleic Acids Res. 46, W486–W494.
Clarke, L. J., and Robinson, S. A. (2008). Cell wall-bound ultraviolet-screening compounds explain the high ultraviolet tolerance of the Antarctic moss, Ceratodon purpureus. New Phytol. 179, 776–783. doi: 10.1111/j.1469-8137.2008.02499.x
Coley, P. D., Bryant, J. P., and Chapin, F. S. III (1985). Resource availability and plant antiherbivore defense. Science 230, 895–899. doi: 10.1126/science.230.4728.895
Cyboran-Mikołajczyk, S., Jurkiewicz, P., Hof, M., and Kleszczyńska, H. (2018). The impact of O-glycosylation on cyanidin interaction with POPC membranes: structure-activity relationship. Molecules 23:2771. doi: 10.3390/molecules23112771
de Gaulejac, N. S., Vivas, N., de Freitas, V., and Bourgeois, G. (1999). The influence of various phenolic compounds on scavenging activity assessed by an enzymatic method. J. Sci. Food Agric. 79, 1081–1090. doi: 10.1002/(sici)1097-0010(199906)79:8<1081::aid-jsfa330>3.0.co;2-g
Dixon, R. A., De-Yu, X., and Sharma, S. B. (2005). Proanthocyanidins-a final frontier in flavonoid research? New Phytol. 165, 9–28. doi: 10.1111/j.1469-8137.2004.01217.x
Dwivedi, S. L., Upadhyaya, H. D., Chung, I., De Vita, P., Garcia-Lara, S., Guajardo-Flores, D., et al. (2016). Exploiting phenylpropanoid derivatives to enhance the nutraceutical values of cereals and legumes. Front. Plant Sci. 7:763. doi: 10.3389/fpls.2016.00763
Fernandez-Orozco, R., Li, L., Harflett, C., Shewry, P. R., and Ward, J. L. (2010). Effects of environment and genotype on phenolic acids in wheat in the HEALTHGRAIN diversity screen. J. Agric. Food Chem. 58, 9341–9352. doi: 10.1021/jf102017s
Fernie, A. R., Bachem, C. W. B., Helariutta, Y., Neuhaus, H. E., Prat, S., Ruan, Y.-L., et al. (2020). Synchronization of developmental, molecular and metabolic aspects of source-sink interactions. Nat. Plants 6, 55–66. doi: 10.1038/s41477-020-0590-x
Ferreyra, M. F., Rius, S. P., and Casati, P. (2012). Flavonoids: biosynthesis, biological functions, and biotechnological applications. Front. Plant Sci. 3:222. doi: 10.3389/fpls.2012.00222
Fischer, T. C., Mirbeth, B., Rentsch, J., Sutter, C., Ring, L., Flachowsky, H., et al. (2014). Premature and ectopic anthocyanin formation by silencing of anthocyanidin reductase in strawberry (Fragaria× ananassa). New Phytol. 201, 440–451. doi: 10.1111/nph.12528
Fritz, C., Palacios-Rojas, N., Feil, R., and Stitt, M. (2006). Regulation of secondary metabolism by the carbon-nitrogen status in tobacco: nitrate inhibits large sectors of phenylpropanoid metabolism. Plant J. 46, 533–548. doi: 10.1111/j.1365-313x.2006.02715.x
Griesser, M., Hoffmann, T., Bellido, M. L., Rosati, C., Fink, B., Kurtzer, R., et al. (2008). Redirection of flavonoid biosynthesis through the down-regulation of an anthocyanidin glucosyltransferase in ripening strawberry fruit. Plant Physiol. 146, 1528–1539. doi: 10.1104/pp.107.114280
Gutierrez, E., García-Villaraco, A., Lucas, J. A., Gradillas, A., Gutierrez-Mañero, F. J., and Ramos-Solano, B. (2017). Transcriptomics, targeted metabolomics and gene expression of blackberry leaves and fruits indicate flavonoid metabolic flux from leaf to red fruit. Front. Plant Sci. 8:472. doi: 10.3389/fpls.2017.00472
Hanhineva, K., Rogachev, I., Kokko, H., Mintz-Oron, S., Venger, I., Kärenlampi, S., et al. (2008). Non-targeted analysis of spatial metabolite composition in strawberry (Fragaria×ananassa) flowers. Phytochemistry 69, 2463–2481. doi: 10.1016/j.phytochem.2008.07.009
Heim, K. E., Tagliaferro, A. R., and Bobilya, D. J. (2002). Flavonoid antioxidants: chemistry, metabolism and structure-activity relationships. J. Nutr. Biochem. 13, 572–584. doi: 10.1016/s0955-2863(02)00208-5
Herms, D. A., and Mattson, W. J. (1992). The dilemma of plants: to grow or defend. Q. Rev. Biol. 67, 283–335. doi: 10.1086/417659
Hura, T., Tyrka, M., Hura, K., Ostrowska, A., and Dziurka, K. (2017). QTLs for cell wall-bound phenolics in relation to the photosynthetic apparatus activity and leaf water status under drought stress at different growth stages of triticale. Mol. Genet. Genom. 292, 415–433. doi: 10.1007/s00438-016-1276-y
Igamberdiev, A. U., and Eprintsev, A. T. (2016). Organic acids: the pools of fixed carbon involved in redox regulation and energy balance in higher plants. Front. Plant Sci. 7:1042. doi: 10.3389/fpls.2016.01042
Jin, X., Yang, G., Tan, C., and Zhao, C. (2015). Effects of nitrogen stress on the photosynthetic CO2 assimilation, chlorophyll fluorescence, and sugar-nitrogen ratio in corn. Sci. Rep. 5:9311.
Jones, C. G., and Hartley, S. E. (1999). A protein competition model of phenolic allocation. Oikos 86, 27–44. doi: 10.2307/3546567
Kähkönen, M. P., and Heinonen, M. (2003). Antioxidant capacity of anthocyanins and their aglycons. J. Agric. Food Chem. 51, 628–633. doi: 10.1021/jf025551i
Kårlund, A., Kati, H., Marko, L., McDougall Gordon, J., Derek, S., and Karjalainen, R. O. (2016). Non-targeted metabolite profiling highlights the potential of strawberry leaves as a resource for specific bioactive compounds. J. Sci. Food Agric. 97, 2182–2190. doi: 10.1002/jsfa.8027
Kim, M., Jin, J. S., Kwak, Y., and Hwang, G. (2016). Metabolic response of strawberry (Fragaria x ananassa) leaves exposed to the angular leaf spot bacterium (Xanthomonas fragariae). J. Agric. Food Chem. 64, 1889–1898. doi: 10.1021/acs.jafc.5b05201
Koehler, G., Rohloff, J., Wilson, R. C., Kopka, J., Erban, A., Winge, P., et al. (2015). Integrative “omic” analysis reveals distinctive cold responses in leaves and roots of strawberry, Fragaria× ananassa ‘Korona’. Front. Plant Sci. 6:826. doi: 10.3389/fpls.2015.00826
Larbat, R., Olsen, K. M., Slimestad, R., Løvdal, T., Bénard, C., Verheul, M., et al. (2012). Influence of repeated short-term nitrogen limitations on leaf phenolics metabolism in tomato. Phytochemistry 77, 119–128. doi: 10.1016/j.phytochem.2012.02.004
Luo, P., Shen, Y., Jin, S., Huang, S., Cheng, X., Wang, Z., et al. (2016). Overexpression of Rosa rugosa anthocyanidin reductase enhances tobacco tolerance to abiotic stress through increased ROS scavenging and modulation of ABA signaling. Plant Sci. 245, 35–49. doi: 10.1016/j.plantsci.2016.01.007
Marković, J. M. D., Milenković, D., Amić, D., Popović-Bijelić, A., Mojović, M., Pašti, I. A., et al. (2014). Energy requirements of the reactions of kaempferol and selected radical species in different media: towards the prediction of the possible radical scavenging mechanisms. Struct. Chem. 25, 1795–1804. doi: 10.1007/s11224-014-0453-z
Maroli, A. S., Nandula, V. K., Dayan, F. E., Duke, S. O., Gerard, P., and Tharayil, N. (2015). Metabolic profiling and enzyme analyses indicate a potential role of antioxidant systems in complementing glyphosate resistance in an Amaranthus palmeri biotype. J. Agric. Food Chem. 63, 9199–9209. doi: 10.1021/acs.jafc.5b04223
Martinez, V., Mestre, T. C., Rubio, F., Girones-Vilaplana, A., Moreno, D. A., Mittler, R., et al. (2016). Accumulation of flavonols over hydroxycinnamic acids favors oxidative damage protection under abiotic stress. Front. Plant Sci. 7:838. doi: 10.3389/fpls.2016.00838
Michalska, A., Carlen, C., Heritier, J., and Andlauer, W. (2017). Profiles of bioactive compounds in fruits and leaves of strawberry cultivars. J. Berry Res. 7, 71–84. doi: 10.3233/jbr-160146
Natella, F., Nardini, M., Di Felice, M., and Scaccini, C. (1999). Benzoic and cinnamic acid derivatives as antioxidants: structure– activity relation. J. Agric. Food Chem. 47, 1453–1459. doi: 10.1021/jf980737w
Nishizawa, A., Yabuta, Y., and Shigeoka, S. (2008). Galactinol and raffinose constitute a novel function to protect plants from oxidative damage. Plant Physiol. 147, 1251–1263. doi: 10.1104/pp.108.122465
Ogo, Y., Mori, T., Nakabayashi, R., Saito, K., and Takaiwa, F. (2015). Transgenic rice seed expressing flavonoid biosynthetic genes accumulate glycosylated and/or acylated flavonoids in protein bodies. J. Exper. Bot. 67, 95–106. doi: 10.1093/jxb/erv429
Olsen, K. M., Slimestad, R., Lea, U. S., Brede, C., Løvdal, T., Ruoff, P., et al. (2009). Temperature and nitrogen effects on regulators and products of the flavonoid pathway: experimental and kinetic model studies. Plant Cell Environ. 32, 286–299. doi: 10.1111/j.1365-3040.2008.01920.x
Pannala, A. S., Chan, T. S., O’Brien, P. J., and Rice-Evans, C. A. (2001). Flavonoid B-ring chemistry and antioxidant capacity: fast reaction kinetics. Biochem. Biophys. Res. Commun. 282, 1161–1168. doi: 10.1006/bbrc.2001.4705
Paul, M. J., and Pellny, T. K. (2003). Carbon metabolite feedback regulation of leaf photosynthesis and development. J. Exper. Bot. 54, 539–547. doi: 10.1093/jxb/erg052
Petrussa, E., Braidot, E., Zancani, M., Peresson, C., Bertolini, A., Patui, S., et al. (2013). Plant flavonoids-biosynthesis, transport and involvement in stress responses. Int. J. Mol. Sci. 14, 14950–14973. doi: 10.3390/ijms140714950
Quan, X., Qian, Q., Ye, Z., Zeng, J., Han, Z., and Zhang, G. (2016). Metabolic analysis of two contrasting wild barley genotypes grown hydroponically reveals adaptive strategies in response to low nitrogen stress. J. Plant Physiol. 206, 59–67. doi: 10.1016/j.jplph.2016.07.020
Reyes, T. H., Scartazza, A., Lu, Y., Yamaguchi, J., and Guglielminetti, L. (2016). Effect of carbon/nitrogen ratio on carbohydrate metabolism and light energy dissipation mechanisms in Arabidopsis thaliana. Plant Physiol. Biochem. 105, 195–202. doi: 10.1016/j.plaphy.2016.04.030
Rysin, O., McWhirt, A., Fernandez, G., Louws, F. J., and Schroeder-Moreno, M. (2015). Economic viability and environmental impact assessment of three different strawberry production systems in the southeastern United States. Hort. Technol. 25, 585–594. doi: 10.21273/horttech.25.4.585
Salminen, J. P., and Maarit, K. (2011). Chemical ecology of tannins and other phenolics: we need a change in approach. Funct. Ecol. 25, 325–338. doi: 10.1111/j.1365-2435.2010.01826.x
Salminen, J. P., Ossipov, V., Haukioja, E., and Pihlaja, K. (2001). Seasonal variation in the content of hydrolysable tannins in leaves of Betula pubescens. Phytochemistry 57, 15–22. doi: 10.1016/s0031-9422(00)00502-1
Shahidi, F., and Yeo, J. (2016). Insoluble-bound phenolics in food. Molecules 21:1216. doi: 10.3390/molecules21091216
Sroka, Z., Żbikowska, B., and Hładyszowski, J. (2015). The antiradical activity of some selected flavones and flavonols: experimental and quantum mechanical study. J. Mol. Model. 21:307.
Strissel, T., Halbwirth, H., Hoyer, U., Zistler, C., Stich, K., and Treutter, D. (2008). Growth-Promoting nitrogen nutrition affects flavonoid biosynthesis in young apple (Malus domestica Borkh.) leaves. Plant Biol. 7, 677–685. doi: 10.1055/s-2005-872989
Stumpf, B., Yan, F., and Honermeier, B. (2015). Nitrogen fertilization and maturity influence the phenolic concentration of wheat grain (Triticum aestivum). J. Plant Nutr. Soil Sci. 178, 118–125. doi: 10.1002/jpln.201400139
Suseela, V., and Tharayil, N. (2018). Decoupling the direct and indirect effects of climate on plant litter decomposition: accounting for stress-induced modifications in plant chemistry. Glob. Chang. Biol. 24, 1428–1451. doi: 10.1111/gcb.13923
Suseela, V., Tharayil, N., Xing, B., and Dukes, J. S. (2015). Warming and drought differentially influence the production and resorption of elemental and metabolic nitrogen pools in quercus rubra. Glob. Chang. Biol. 21, 4177–4195. doi: 10.1111/gcb.13033
Tanou, G., Minas, I. S., Scossa, F., Belghazi, M., Xanthopoulou, A., Ganopoulos, I., et al. (2017). Exploring priming responses involved in peach fruit acclimation to cold stress. Sci. Rep. 7:11358.
Tegeder, M., and Masclaux-Daubresse, C. (2018). Source and sink mechanisms of nitrogen transport and use. New Phytol. 217, 35–53. doi: 10.1111/nph.14876
Top, S. M., Preston, C. M., Dukes, J. S., and Tharayil, N. (2017). Climate influences the content and chemical composition of foliar tannins in green and senesced tissues of Quercus rubra. Front. Plant Sci. 8:423. doi: 10.3389/fpls.2017.00423
Winkel-Shirley, B. (2002). Biosynthesis of flavonoids and effects of stress. Curr. Opin. Plant Biol. 5, 218–223. doi: 10.1016/s1369-5266(02)00256-x
Wurtzel, E. T., and Kutchan, T. M. (2016). Plant metabolism, the diverse chemistry set of the future. Science 353, 1232–1236. doi: 10.1126/science.aad2062
Xie, D., and Dixon, R. A. (2005). Proanthocyanidin biosynthesis - still more questions than answers? Phytochemistry. 66, 2127–2144. doi: 10.1016/j.phytochem.2005.01.008
Xie, J., and Schaich, K. M. (2014). Re-evaluation of the 2,2-diphenyl-1-picrylhydrazyl free radical (DPPH) assay for antioxidant capacity. J. Agric. Food Chem. 62, 4251–4260. doi: 10.1021/jf500180u
Zhang, Y., De Stefano, R., Robine, M., Butelli, E., Bulling, K., Hill, L., et al. (2015). Different reactive oxygen species scavenging properties of flavonoids determine their abilities to extend the shelf life of Tomato. Plant Physiol. 169, 1568–1583. doi: 10.1104/pp.15.00346
Zhao, J., Pang, Y., and Dixon, R. A. (2010). The mysteries of proanthocyanidin transport and polymerization. Plant Physiol. 153, 437–443. doi: 10.1104/pp.110.155432
Zhou, H., Tam, N. F., Lin, Y., Ding, Z., Chai, W., and Wei, S. (2014). Relationships between degree of polymerization and antioxidant activities: a study on proanthocyanidins from the leaves of a medicinal mangrove plant Ceriops tagal. PLoS One 9:e107606. doi: 10.1371/journal.pone.0107606
Keywords: nitrogen, polyphenols, flavonoids, tannins, metabolomics
Citation: Narvekar AS and Tharayil N (2021) Nitrogen Fertilization Influences the Quantity, Composition, and Tissue Association of Foliar Phenolics in Strawberries. Front. Plant Sci. 12:613839. doi: 10.3389/fpls.2021.613839
Received: 15 October 2020; Accepted: 08 January 2021;
Published: 20 April 2021.
Edited by:
Bruno Mezzetti, Marche Polytechnic University, ItalyReviewed by:
María Teresa Ariza-Fernández, Andalusian Institute for Research and Training in Agriculture, Fisheries, Food and Ecological Production (IFAPA), SpainJose G. Vallarino, University of Malaga, Spain
Copyright © 2021 Narvekar and Tharayil. 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(s) 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: Nishanth Tharayil, bnRoYXJheUBjbGVtc29uLmVkdQ==
†Present address: Ashwini Sushil Narvekar, Ginkgo Bioworks, Boston, MA, United States