Next Article in Journal
Effect of Stand Reduction at Different Growth Stages on Yield of Paprika-Type Chile Pepper
Next Article in Special Issue
Use of In Situ Soil Solution Electric Conductivity to Evaluate Mineral N in Commercial Orchards: Preliminary Results
Previous Article in Journal
Quantitative Trait Loci (QTL) Analysis of Fruit and Agronomic Traits of Tropical Pumpkin (Cucurbita moschata) in an Organic Production System
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Genetic Variation in Response to N, P, or K Deprivation in Baby Leaf Lettuce

Crop Improvement and Protection Research Unit, United States Department of Agriculture, Agricultural Research Service, Salinas, CA 93905, USA
Horticulturae 2020, 6(1), 15; https://doi.org/10.3390/horticulturae6010015
Submission received: 31 December 2019 / Revised: 2 February 2020 / Accepted: 13 February 2020 / Published: 3 March 2020
(This article belongs to the Special Issue Fertilization Management of Horticultural Crops)

Abstract

:
Lettuce harvested at the baby leaf stage is a popular component of mixed salads in ready-to-use packages. Little is known, however, about response of baby leaf lettuce to nitrogen (N), phosphorus (P), and potassium (K) fertilization. Eight lettuce accessions were subjected to five fertilization treatments to investigate genetic differences in reaction to N, P, and K fertilization. The control treatment provided optimal levels of macronutrients for plant growth, while other treatments deprived plants of either N, P, or K. Potassium deprivation had no obvious effect on plant weight or composition, apart from substantially decreased potassium content. Nitrogen and phosphorus deprivations, however, extensively decreased fresh weight and affected plant composition. Phosphorus and nitrogen deprivation considerably increased anthocyanin content in red-colored accessions, but anthocyanin was decreased in dark green-colored accessions, indicating differences in regulation of anthocyanin biosynthesis. Correlations between fresh weight, chlorophyll, anthocyanin, nitrogen, phosphorus, and potassium content were substantially affected by selection of datasets used for analyses; some relationships were revealed when analyzed separately by individual treatments, while others were more likely to be detected when analyzed by individual accessions. Absolute (ΔABS) and relative (2ΔREL) parameters described in this study were suitable for detecting over- and underperforming accessions. The ΔABS identified the absolute Lb-fold (logarithm to the base of 2, binary logarithm) change in performance of an accession in a treatment as compared to its performance in control conditions. The 2ΔREL parameter showed relative Lb-fold change for an accession as compared to the overall mean of ΔABS values of all accessions tested in control and treatment conditions.

Graphical Abstract

1. Introduction

Lettuce (Lactuca sativa L.) for baby leaf production is harvested approximately 30 days after planting when the first four true leaves reach a length of about 5–13 cm [1]. Baby leaf lettuce is the primary component of spring mix salads that are popular in ready-to-use packages [2]. Raw lettuce leaves contain dietary fiber, several important dietary minerals, vitamins (e.g., vitamins A, B9, and C), and bioactive compounds (e.g., chlorophyll, carotenoids, and phenolic compounds) [3,4,5] that contribute to human nutritional benefits [6]. Fertilization of lettuce with different rates of nitrogen (N), phosphorus (P), and potassium (K) affects plant yield, nutritional quality [7,8], postharvest quality [9,10], resistance to pathogens [11], microbiota in phyllosphere [12], and possibly also survival of human enteric pathogens on leaves [13]. In baby spinach, biomass yield, chlorophyll content, and leaf area increased significantly with increasing rates of N and P; while K had no effect on the yield, chlorophyll content, stomatal conductance, and leaf area [14]. In lettuce, chlorophyll concentration increased with increasing levels of N [7,15]. Anthocyanin accumulation in plants is most commonly associated with P deficiency [16], though deficiencies in N and K have also been noted to induce anthocyanin formation [17,18].
Lettuce requires a sufficient supply of essential macronutrients for optimal production, which are frequently provided in the form of mineral fertilizers. To minimize the environmental impact and financial cost associated with using fertilizers, it is important to improve the efficiency with which macronutrients from fertilizers are used by plants [19]. Many of the modern crops have low fertilizer use efficiency, in particular nitrogen use efficiency, because they have been bred in the presence of non-limiting mineral fertilization conditions [20]. To develop breeding lines with high fertilizer use efficiency, it will be necessary to screen a range of genotypes for their responses to limited supply of N, P, and K, and to assess relationships between supplied macronutrients and plant composition. Genotypes that can produce high yield and maintain good quality at low levels of macronutrients may be selected for future studies and integrated into breeding programs focused on developing lettuces with improved fertilizer use efficiency. Although several studies were previously performed to investigate the effect of N, P, and K fertilization on yield and quality of leafy vegetables, data that describe genetic variation in response to mineral nutrition in baby leaf lettuce are lacking. Therefore, the objectives of this study were (1) to assess genetic differences among eight lettuce accessions in their response to N, P, or K deprivation at the baby leaf stage; and (2) to determine the relationship between plant fresh weight and content of chlorophyll, anthocyanin, N, P, and K in leaves at both optimal and suboptimal growing conditions. Four N, P, and K, deprivation treatments used in this study were aimed to markedly change plant composition, thus allowing a determination of how altered plant composition affected correlations between traits.

2. Materials and Methods

Five different treatments were used in the present study to investigate genetic differences in reaction to nutrients: one treatment provided optimal levels of macronutrients for plant growth; and four treatments where plants were deprived of either N, P, or K. Differences in the performance of accessions at four suboptimal conditions were analyzed using several methodologies, including an absolute (ΔABS) parameter and a novel, relative (2ΔREL) parameter that compares relative performances of all accessions under the particular treatment after adjusting for their performance under control conditions (optimal levels of macronutrients). In this study, the term “performance” is used to indicate the change in values of an accession with the particular fertilizer treatment as compared to the values measured on plants of the same accession grown under control conditions (i.e., changes in yield, and in anthocyanin, chlorophyll, N, P, and K contents).

2.1. Plant Materials and Growth Conditions

A diverse set of eight lettuce accessions evaluated in this study included seven cultivars, namely ‘Annapolis’ (romaine, dark red), ‘Balady Barrage’ (B. Barrage; stem type, green), ‘Darkland’ (romaine, dark green), ‘Eruption’ (Latin, red), ‘La Brillante’ (Batavia, light green), ‘Merlot’ (leaf type, dark red), ‘Salinas’ (iceberg, green), and a single accession of prickly lettuce UC96US23 (Lactuca serriola L., dark green), which is a species from which desirable traits and alleles can be introgressed into cultivated lettuce. Lettuce seeds were planted in potting soil (Premium Growers Mix, Sun Land Garden Products, Watsonville, CA, USA), covered with sand, and then watered. Trays with seeds were kept for 48 h at 10 °C in the dark to improve uniformity of germination. Afterwards, trays were transferred to a growth room with a temperature of 20 °C and 16 h/8 h light/dark photoperiod for germination and initial growth. Established plants of uniform development were transplanted to 7.6 cm pots (~514 mL volume) containing 1:2 mix of Espoma VM8 8-Quart Organic Vermiculite (Espoma, Millville, NJ, USA) and sand. Plants were fertilized with Non-Coated Ammonium Nitrate 34-0-0 Prill Form Fertilizer (Intermountain Farmers Association, Salt Lake City, UT, USA), Triple Super Phosphate 0-46-0 Easy Peasy Plants 99% pure (Easy Peasy Plants, Alvin, IL, USA), and All-Natural Muriate of Potash-Easy Peasy 0-0-60 Potassium (Easy Peasy Plants, Alvin, IL, USA) to achieve concentrations of N, P, and K shown in Table 1. All plants were sprayed with an identical amount of micronutrient mix solution, which was prepared by dissolving 1 g of Axilo Mix 5 (0-0-0) Micronutrient Mix (Valagro USA, Houston, TX, USA) in 1 L of distilled water. Pots were then watered to saturation. Plants were grown in a greenhouse in a randomized complete block design with four replications (8 × 5 × 4 = 160 plants in total) and watered daily as needed. The average daily temperature in the greenhouse (April and May 2018) ranged from 20 to 24 °C, daylength ranged from approximately 13 to 14 h, and outdoor average daily light integrals were between 40 and 55 mol m−2 d−1. When four true leaves on the majority of control plants (NPK treatments) reached about 10 cm, all plants were harvested. Each of the four individual plants collected per accession was used for quantification of compounds.

2.2. Quantification of Fresh Weigh, Pigments, and N, P, K

Fresh weight of above ground parts was quantified immediately after harvest. The relative content of chlorophyll and anthocyanin was determined two days before harvest using SPAD-502 (Spectrum Technologies, Plainfield, IL, USA) and ACM-200 plus (Opti-Sciences, Hudson, NH, USA) hand-held meters, respectively. These devices use light transmittance to provide good in situ estimates of relative contents of the two pigments [21,22]. Chlorophyll and anthocyanin were measured on three leaves of similar age (avoiding youngest and oldest leaves) and size per plant. The measuring clip was positioned about 1 cm from the edge of the leaf, ensuring that major veins were avoided. The content of chlorophyll is expressed in SPAD units; the content of anthocyanin is expressed in anthocyanin content index (ACI) units. For each plant, the averages of three measurements of chlorophyll and anthocyanin were recorded and used in statistical analyses.
Leaf tissue of harvested plants was used for quantification of N, P, and K contents. All laboratory analyses were performed by UC Davis Analytical Laboratory (https://anlab.ucdavis.edu) using common quantification procedures developed for plant tissue samples. Prior to NPK analyses, samples were dried and ground to pass through a 40-mesh screen (<0.4 mm) on the Wiley Mini-Mill (Thomas Scientific, Swedesboro, NJ, USA). N was quantified using a method that measures soluble nitrate (NO3-N) in plant materials based on an extraction with a solution of 2% acetic acid [23]. Nitrate was determined by a Quikchem 8500 Flow Injection Analysis (FIA) system (Lachat Instruments, Loveland, CO, USA) using the reduction to nitrite via a copperized cadmium column, diazotization with sulfanilamide, followed by coupling with N-(1-naphthyl) ethlyenediaminie dihydrochloride. The absorbance of the product was measured at 520 nm. P was quantified using a method that determines the amount of phosphorus (PO4-P) in plant materials by extraction with a solution of 2% acetic acid [24]. Phosphorus concentration in the extract was determined spectrophotometrically by reacting with ammonium molybdate and antimony potassium tartrate under acidic conditions to form a complex. This complex was reduced with ascorbic acid to form a blue complex that absorbs light at 880 nm. The absorbance is proportional to the concentration of phosphorus in the sample. Samples were analyzed using an automated Quikchem 8500 FIA system. K was quantified using a method that measures soluble potassium (K) present in plant materials by extraction with a solution of 2% acetic acid [25]. Potassium was quantitatively determined in the extract using inductively coupled plasma atomic emission spectrometry (ICP-AES) with a Thermo Scientific iCAP 6000 (Waltham, MA, USA). N, P, and K content are reported in g per 1 kg of fresh weight (g kg−1 FW).

2.3. Data Analyses

Three statistical approaches were used to compare the performance of eight accessions tested in five environments. The distinctions of those statistical approaches can be explained by a simplified example that considers only three accessions, namely A, B, and C, grown in control (AC = 10, BC = 15, CC = 30) and treatment (AT = 20, BT = 25, CT = 34) environments, where numerical values indicate the averages of their respective performances. Statistical analyses were performed on data from all accessions grown in all environments (e.g., AC, BC, CC, AT, BT, CT) using analysis of variance (ANOVA) followed by a multiple comparison procedure (Tukey’s honest significant difference [HSD] test) to identify significant differences. However, because this kind of post hoc test compares all possible pairs of means, when eight accessions are tested in five environments, the approach leads to a very large number or pairwise comparisons (8 × 5 × (8 × 5 − 1)/2 = 780) that affect the experiment-wise error rate. Moreover, testing of all paired comparisons is not always necessary or relevant (e.g., comparing performance of a cultivar grown in one environment with another cultivar grown in a different environment).
Two other statistical tests were performed on absolute (ΔABS) and relative (2ΔREL) differences in values detected for an accession grown using a specific treatment and control environment. Here, ΔABS denotes an absolute difference in performance of an accession in two environments. Using the above example, the ΔABS for accession A is calculated as: AT − AC, 20 − 10 = 10. Because only performance of a single accession in two environments is compared using this approach, the ΔABS parameter is independent from performances of other accessions tested in the same experiments. Statistical analyses of ΔABS significance in this study were performed using Student’s t-test. Here, 2ΔREL denotes a relative difference in performance of an accession and is calculated by comparing the ΔABS of all accessions tested in the same two environments (control and a particular treatment). Therefore, this parameter shows a relative performance that is affected by the performance of other tested accessions. When the goal is to compare the ΔABS of accession A to the overall mean of all ΔABS, the 2ΔREL for this accession is calculated as ΔABS of A minus the average of ΔABS of A, ΔABS of B, and ΔABS of C. Using the above example, 2ΔREL for A equals (20 − 10) – ((20 − 10) + (25 − 15) + (34 − 30))/3 = 2. The sum of all 2ΔREL values equals zero. The 2ΔREL values detected in this study were statistically evaluated using analysis of means (ANOM). Alternative statistical analyses of ΔABS values could be performed using multiple comparisons procedures.
Pairwise correlations between values of six tested traits were performed on 14 datasets to identify how treatments and genotypes affect relationships among traits. Eight of the datasets comprised data from a single accession each, tested for all five treatments. These datasets show how correlations differ among accessions tested across all environments. Five of the datasets comprised data from a single treatment each, but using all eight accessions. These datasets show how correlations differ among treatments when data from all accessions are considered. Analyses were also performed on the complete dataset, which included all combinations of tested treatments and accessions, to identify correlations that could be detected in combined data.
Principal component analysis (PCA) of traits was done using the correlation matrix as an input dataset. Effect size (ω2), which estimates how much variance in the response variable is accounted for by individual factors (accessions and treatments) and their interaction, was calculated from ANOVA output [26]. All statistical analyses were performed using JMP software v. 14.2.0 (SAS Institute, Cary, NC, USA) and Microsoft Excel for Mac v. 16.16.5 (Microsoft, Redmond, WA, USA).

2.4. Data Transformation

It has been previously observed that lognormal distribution in plant weight develops naturally [27], and therefore assumption of lognormality of this trait is probably appropriate in most applications [28]. Logarithmic transformation of weight and other traits is not performed only because of statistical convenience, but it is appropriate due to the fundamentally multiplicative nature of biological phenomena [29]. In the present study, data for all six traits were statistically analyzed using original (measured) values and also logarithmically transformed values (logarithm to the base of 2, binary logarithm, Lb). Because statistical power for all but one trait was greater when using Lb-transformed values, only these results are presented in this paper (original values are shown in Supplementary Table S1). The only exception was the SPAD index, which showed practically identical statistical power when using both the original and the transformed values. However, to keep uniformity of data across all traits, only results obtained using Lb-transformed SPAD values are shown. “Lb” was added in front of trait names to indicate that data were logarithmically transformed (Lb FW, Lb SPAD, Lb ACI, Lb N, Lb P, and Lb K). Use of logarithmic transformation means that the difference between two “Lb” values is the logarithm of the ratio of two original values (e.g., Lb (AT) − Lb (AC) = Lb (AT/AC)). Therefore, the difference between two transformed values equals the 2difference-fold change between the original values (e.g., for AT = 20 and AC = 10, Lb 20 − Lb 10 = 1, which in turn is 2 1 = 2-fold change from AC = 10 to AT = 20, or 20/10 = 2). The use of transformed values also means that back transformation of averages calculated from logarithms equals geometric means calculated from original values.

3. Results

3.1. Effect of Treatments

Plants of eight accessions were grown at control (NPK) conditions, with all nutrients provided at an optimal level and four treatments having a limited amount of nitrogen (nPK and PK), phosphorus (NK), and potassium (NP), with treatments having a profound effect on six analyzed traits. Results of two-way ANOVA (Table 2) indicate that compared to control, partial N restriction (nPK) led to lower fresh weight (Lb FW) and lower content of N (Lb N). Total N deprivation (PK) decreased values of all six tested traits: plant fresh weight (Lb FW), content of chlorophyll (Lb SPAD), anthocyanin (Lb ACI), nitrogen (Lb N), phosphorus (Lb P), and potassium (Lb K). P deprivation (NK) reduced the fresh weight (Lb FW), content of chlorophyll (Lb SPAD), phosphorus (Lb P), and potassium (Lb K), but increased the content of anthocyanin (Lb ACI) and nitrogen (Lb N). K deprivation (NP) did not show any significant effect on change of tested traits compared to control, with the exception of potassium content (Lb K), which was decreased. PCA revealed that data from NK and PK treatments were generally well separated from data of the other three treatments and from each other (Figure 1a). Overall, there was a close relationship detected between Lb SPAD and Lb N traits (Figure 1b).

3.2. Fresh Weight

Fresh weight was influenced by both accessions and nutrient treatment. Fresh weight, expressed as a binary logarithm of FW in g (Lb FW), ranged from 0.80 for UC96US23 under PK conditions to 5.93 for ‘Darkland’ under NPK conditions (Table 3). The overall Lb FW mean for accessions was lowest for UC96US23 (2.63) and highest for ‘Darkland’ (4.16), while the overall mean for treatments ranged from 1.60 under PK conditions to 5.03 under NPK conditions (Table 2). A large effect size was found for treatments (ω2 = 0.852), while small effect size was found for accessions (ω2 = 0.076) and treatment by accession interaction (ω2 = 0.009). All accessions had significant negative ΔABS under PK and NK treatments. B. Barrage and ‘Darkland’ also had significant negative ΔABS values with nPK. No significant ΔABS for Lb FW was detected with NP (Figure 2). The 2ΔREL analysis revealed that fresh weight of ‘Darkland’ declined more than in other tested accessions with PK treatment.

3.3. Chlorophyll Content Index (SPAD)

The chlorophyll content index, expressed as the binary logarithm of SPAD (Lb SPAD), ranged from 3.82 for ‘Merlot’ with PK to 5.58 for UC96US23 with NP (Table 3). The overall Lb SPAD mean for accessions was lowest for ‘Merlot’ (4.32) and highest for UC96US23 (5.28), while the overall mean for treatments ranged from 4.32 with PK to 5.11 with NP (Table 2). Effect size was similar for accessions (ω2 = 0.377) and treatments (ω2 = 0.316), but substantially smaller for interaction (ω2 = 0.032). Significant negative ΔABS was detected for all accessions with PK, with the exception of ‘Salinas’, and also for B. Barrage with NK. Significant positive ΔABS, indicating an increase in chlorophyll index, was detected for ‘Merlot’ with NP. No 2ΔREL showed significant difference from zero (Figure 3).

3.4. Anthocyanin Content Index (ACI)

The anthocyanin content index, expressed as the binary logarithm of ACI (Lb ACI), ranged from 1.38 for ‘La Brillante’ with PK to 7.08 for ‘Annapolis’ with NK (Table 3). The overall Lb ACI mean for accessions was lowest for ‘La Brillante’ (1.75) and highest for ‘Annapolis’ (5.48), while the overall mean for treatments ranged from 2.94 with PK to 3.70 with NK (Table 2). The effect size was high for accessions (ω2 = 0.725) but low for treatments (ω2 = 0.037). Lb ACI had the largest effect size for interaction (ω2 = 0.146) from all tested traits. Significant positive ΔABS was detected for ‘Annapolis’, ‘Eruption’, and ‘Merlot’ with both PK and NK, indicating that anthocyanin content in these accessions increased as compared to the control (NPK). In comparison, significant negative ΔABS was detected for B. Barrage, ‘Darkland’, and UC96US23 with PK, and for B. Barrage with NK (Figure 4). All accession–treatment combinations with significant ΔABS values also had significant 2ΔREL values. In addition, 2ΔREL values for ‘Darkland’ and UC96US23 with NK were also significantly lower than the overall mean of ΔABS values.

3.5. Nitrogen Content

Nitrogen content, expressed as a binary logarithm of N content in g kg−1 FW (Lb N), ranged from 0.28 for ‘Darkland’ with PK to 2.78 for UC96US23 with NK (Table 3). The overall Lb N mean for accessions was lowest for ‘Darkland’ (1.50) and highest for UC96US23 (2.18), while the overall mean for treatments ranged from 0.78 with PK to 2.38 with NK (Table 2). Effect sizes for treatments, accessions, and treatment by accession interaction were ω2 = 0.759, ω2 = 0.108, and ω2 = 0.027, respectively. The ΔABS values were negative and significant for all accessions with PK, for all but two accessions (B. Barrage and UC96US23) at nPK, and ‘La Brillante’ at NP. In comparison, significant positive ΔABS values, indicating increase in N content, were observed for B. Barrage with NK and NP, ‘Darkland’ with NK, UC96US23 with NK, and Salinas with NP (Figure 5). Only two 2ΔREL values were significant, namely the increase for B. Barrage and decrease for ‘Merlot’, both with NK.

3.6. Phosphorus Content

Phosphorus content, expressed as the binary logarithm of P content in g kg−1 FW (Lb P), ranged from −3.65 for B. Barrage with NK to −0.21 for ‘Merlot’ with NPK (Table 3). The overall Lb P mean for accessions was lowest for B. Barrage (−1.65) and highest for ‘Merlot’ (−1.03), while the overall mean for treatments ranged from −3.18 with NK to −0.53 with NP (Table 2). The effect size was very high for treatments (ω2 = 0.892) and low for accessions (ω2 = 0.031). No effect size was detected for treatment by accession interaction (ω2 = 0.000). Results for effect size matched with those observed for ΔABS and 2ΔREL values. All accessions had significant negative ΔABS values with PK and NK, while none of the 2ΔREL values were significant, which indicated that all accessions responded very similarly to treatments (Figure 6).

3.7. Potassium Content

Potassium content, expressed as the binary logarithm of K content in g kg−1 FW (Lb K), ranged from 0.84 for B. Barrage with NP to 2.60 for ‘Eruption’ with PK (Table 3). The overall Lb K mean for accessions was lowest for B. Barrage (1.47) and highest for ‘Eruption’ (2.24), while the overall mean for treatments ranged from 1.36 with NP to 2.28 with nPK (Table 2). The effect size was ω2 = 0.435 for treatments, ω2 = 0.320 for accessions, and ω2 = 0.092 for treatment by accession interaction. Significant negative ΔABS values were observed for all accessions with NP, all but two accessions (‘Eruption’ and UC96US23) with NK, and four accessions (B. Barrage, ‘Darkland’, ‘La Brillante’, and ‘Salinas’) with PK (Figure 7). Additionally, a significant positive ΔABS was detected, indicating increases in K content compared to control (NPK) treatment were found for ‘Darkland’ and ‘Eruption’ with nPK, and for ‘Eruption’ with PK. Four 2ΔREL values were significant at p < 0.05, two of them were lower than the overall mean (for ‘Darkland’ and ‘La Brillante’ with PK), and another two were higher than the overall mean (for ‘Eruption’ with PK and NK).

3.8. Correlations Between Traits

Seven of the 15 pairwise, linear correlations between six traits were significant at p < 0.05 when evaluated on data from all accessions and treatments (Figure 8, Table 4). These positive correlations were found between (listed in descending order of r size) Lb FW and Lb P (r = 0.56), Lb FW and Lb SPAD (r = 0.41), Lb SPAD and Lb N (r = 0.41), Lb FW and Lb N (r = 0.40), Lb SPAD and Lb ACI (r = 0.30), Lb N and Lb K (r = 0.19), and Lb P and Lb K (r = 0.18). When correlations were performed using data from each treatment separately (thus showing correlations across accessions at each treatment), all five correlation coefficients were significant between Lb N and Lb P (r from 0.74 to 0.90, Figure 9), Lb P and Lb K (r from 0.72 to 0.88), and Lb N and Lb K (r from 0.62 to 0.97). Four of the correlation coefficients were significant between Lb FW and Lb N (r from −0.43 to −0.54), while the fifth correlation with nPK treatment (r = −0.31) was not significant at P < 0.05. When correlations were performed using data from each accession separately (thus showing correlations across treatments for each accession), all eight correlation coefficients were significant between Lb FW and Lb P (r from 0.54 to 0.77), seven between Lb FW and Lb SPAD (except ‘Salinas’) and Lb SPAD and Lb N (except ‘Salinas’), and six for Lb FW and Lb ACI (three negative for ‘Annapolis’, ‘Eruption’, and ‘Merlot’, and three positive for B. Barrage, ‘Darkland’, and UC96US23). Frequently, correlations that were significant when calculated from individual treatments were not significant when calculated from individual accessions, and vice versa.

4. Discussion

Commercial lettuce production requires adequate levels of NPK to provide high yield and good product quality. The composition of lettuce accessions tested in this study depended on macronutrients used in each treatment (Table 2 and Table 3). Generally, untransformed values (provided in Supplementary Table S1) for SPAD, ACI, N content, P content, and K content matched well with those previously reported for lettuce [3,4,30,31,32]. Fresh weight and SPAD values significantly decreased at low rates of N (PK) and P (NK), but not K (NP) (Table 2). These results agree with those observed for baby spinach, where both biomass yield and chlorophyll content increased significantly with increasing N and P rates, although not with increasing rates of K [14]. Similarly, chlorophyll levels in lettuce increased linearly with N rate from 0 to 960 mg per kg of soil [15].
The overall increase in anthocyanin content (Lb ACI) was associated with phosphorus deprivation with NK treatment (Figure 4). This observation corroborates previous results obtained with diverse plant species, such as petunia (Petunia × hybrida hort. ex E. Vilm.) [33], maize (Zea mays L.) [17], and cabbage (Brassica oleracea L. var. capitata) [18], where higher anthocyanin production was related to phosphorus deficiency. Besides phosphorus deficiency, nitrogen deficiency can also lead to higher anthocyanin accumulation [16,17,18]. Results from this study validate the effect of nitrogen deficiency on increased accumulation of anthocyanin. However, current results also point out differences among tested accessions in their reaction to P and N deprivation (Figure 4). While accessions that appeared as dark-red (‘Annapolis’ and ‘Merlot’) and red (‘Eruption’) colored in control treatment conditions (NPK) all responded to both P and N deprivation by increasing their anthocyanin content, accessions that appeared dark green in control conditions (B. Barrage, ‘Darkland’, and UC96US23) had lower accumulation of anthocyanin than other tested accessions when deprived of P or N. Accessions that were light green (‘La Brillante’) or green (‘Salinas’) in control conditions only had minimal changes in their content of anthocyanin in NK and PK treatments. Differences among accessions in their accumulation of anthocyanin under stress conditions were previously observed under both high and low temperature conditions [34]. At low temperature (3 °C), accumulation of anthocyanin increased as compared to optimal temperature (21 °C), but only in red-colored accessions and not in green ones. In comparison, under high temperature (39 °C), anthocyanin accumulation decreased in red-colored accessions and increased in all that were green-colored [34]. Because several genes are involved in controlling red leaf color in lettuce [35,36,37], it is possible that these genes respond differently to different environmental stimuli. To better understand the regulatory mechanism of anthocyanin biosynthesis in lettuce and investigate why plants with high levels of anthocyanin content responded differently to nutrient stresses, further studies using N and P deficiency treatments are needed.
Correlation coefficients calculated with the complete dataset (all accessions and all fertilizer treatments) frequently differed from those determined from either individual accessions tested across all treatments or from individual treatments that included all accessions (Figure 9). Such differences in results emphasize the importance of using appropriate experimental design when evaluating correlations between traits. For example, a significant positive correlation detected between Lb N and Lb SPAD corroborated the previously reported strong relationship between concentration of N and chlorophyll (as determined by SPAD) [38]. It is possible that such a strong relationship exists because up to 80% of leaf nitrogen is allocated in chloroplasts [39]. However, significant correlation between the two traits was detected only when performed on data from individual accessions tested in multiple environments (treatments), but not when data from all accessions tested in a single treatment were used for analysis. Comparable results were previously observed for rice (Oryza sativa L.) plants tested at different developmental stages [40]. Coefficients of determination between SPAD values and leaf nitrogen concentration at each of the four developmental stages ranged from R2 = 0.61 to 0.88, but decreased to only R2 = 0.11 when data from all developmental stages were pooled [40].
Strong, positive, linear correlations were detected between values of N, P, and K (Lb N, Lb P, and Lb K) when tested in individual treatments, but not when tested using individual accessions. Previously, significant positive correlations were found between N, P, and K contents in grains of rice (r from 0.89 to 0.97) [41] and in plantlets of stevia (Stevia rebaudiana Bertoni) (r from 0.59 to 0.85) [42], indicating that uptake and accumulation of these compounds may be related. It needs to be emphasized, however, that the strong correlations detected between contents of N, P, and K in the present study were determined on a fresh weight basis. It is possible that calculations performed on a dry weight basis would provide somewhat different results.
Use of ΔABS and 2ΔREL parameters allowed the identification of accessions with significantly changed values (performance) with different treatments. While the ΔABS parameter indicated an absolute change in performance of a particular accession, the 2ΔREL parameter indicated a relative change compared to other tested accessions. The difference between the two approaches can be seen, for example, by the decrease of fresh weight (Lb FW) with PK treatment (Figure 2). All accessions had a significant decrease in their ΔABS values; however, only ‘Darkland’ had a significant decrease in its 2ΔREL value, indicating that the fresh weight of this accession decreased more than the overall average of other accessions. Another example is the change in potassium content (Lb K) with NK treatment for ‘Eruption’ (Figure 7). Although a change was not detected when ΔABS was tested, the 2ΔREL value was higher than the overall mean. Using the 2ΔREL parameter as an indicator, a decrease was detected for ‘Darkland’ fresh weight with PK, ‘Merlot’ nitrogen content with NK, ‘Darkland’ and ‘La Brillante’ potassium content with PK, and B. Barrage, ‘Darkland’, and UC96US23 anthocyanin content with both PK and NK. On the other hand, higher than overall values of 2ΔREL were found for B. Barrage nitrogen content with NK, ‘Eruption’ potassium content with both PK and NK, and ‘Annapolis’, ‘Eruption’, and ‘Merlot’ anthocyanin content with both PK and NK (Figure 2, Figure 3, Figure 4, Figure 5, Figure 6 and Figure 7). Values of ΔABS and 2ΔREL can, thus, be used to identify over- or underperforming accessions. These two approaches complement each other by helping to find lettuce accessions with absolute or relative (Lb-fold) changes in their performance. Results can then be used to select accessions with the most desirable traits. Additional accessions of both cultivated and prickly lettuce should be tested using different fertilizer treatments to identify genotypes for which the composition is minimally affected by decreased fertilization with N, P, and K. Such genotypes could then be used to study nitrogen, phosphorus, and potassium accumulation in lettuce and to develop novel cultivars with improved efficiency of absorption and use of these macronutrients. Improving the efficiency of N, P, and K utilization would lead to more economic production of lettuce, while enhancing environmental sustainability.

5. Conclusions

Results of the present study were as follows.
Potassium deprivation had no effect (or no consistent effect across lettuce accessions) on any of the tested traits, with the exception of potassium content in fresh tissue, which substantially and significantly decreased.
Nitrogen and phosphorus deprivations extensively decreased fresh weight and affected plant composition.
The largest effects of treatments were observed on content of phosphorus and fresh weight; the largest effects of accessions were detected on anthocyanin content; and the largest effects of interactions were also found for anthocyanin content.
In phosphorus and nitrogen deprivation treatments, anthocyanin substantially increased in red-colored accessions, but decreased in those that were dark green-colored, indicating differences in the regulatory mechanism(s) of anthocyanin biosynthesis.
Potassium content in ‘Eruption’ lettuce increased with nitrogen and phosphorus deprivation treatments.
Correlations between traits were substantially affected by selection of datasets used for analyses. Some of the correlations were very high (e.g., between content of nitrogen, phosphorus, and potassium) when individual treatments were analyzed separately but not with combined data, while others were more likely to be detected when analyzed with individual accessions (e.g., between fresh weight and anthocyanin content, although in opposite directions with different accessions).
Parameters ΔABS and 2ΔREL, as described in this study, allow the detection of over- or underperforming accessions. These two parameters indicated Lb-fold changes in absolute and relative performance of an accession, were complementary to each other, and should be considered jointly.

Supplementary Materials

The following are available online at https://www.mdpi.com/2311-7524/6/1/15/s1. Table S1: Original measured values of fresh weight, and content of chlorophyll, anthocyanin, nitrogen, phosphorus, and potassium for eight lettuce accessions tested with five fertilizer treatments and harvested at baby leaf stage.

Funding

This research received no external funding.

Acknowledgments

The author would like to thank L. Rosental for excellent technical assistance in various phases of this research.

Conflicts of Interest

The author declares no conflict of interest.

Disclaimer

The mentioning of trade names or commercial products in this publication is solely for the purpose of providing specific information and does not imply recommendation or endorsement by the U.S. Department of Agriculture (USDA).

References

  1. Simko, I.; Hayes, R.J.; Mou, B.; McCreight, J.D. Lettuce and Spinach. In Yield Gains in Major U.S. Field Crops; CSSA Special Publications; Smith, S., Diers, B., Specht, J., Carver, B., Eds.; American Society of Agronomy, Inc., Crop Science Society of America, Inc., Soil Science Society of America, Inc.: Madison, WI, USA, 2014; pp. 53–86. [Google Scholar]
  2. Saini, R.K.; Ko, E.Y.; Keum, Y.-S. Minimally processed ready-to-eat baby-leaf vegetables: Production, processing, storage, microbial safety, and nutritional potential. Food Rev. Int. 2017, 33, 644–663. [Google Scholar] [CrossRef]
  3. Mou, B. Nutrient content of lettuce and its improvement. Curr. Nutr. Food Sci. 2009, 5, 242–248. [Google Scholar] [CrossRef]
  4. Simko, I. Genetic variation and relationship among content of vitamins, pigments, and sugars in baby leaf lettuce. Food Sci. Nutr. 2019, 7, 3317–3326. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Sytar, O.; Zivcak, M.; Bruckova, K.; Brestic, M.; Hemmerich, I.; Rauh, C.; Simko, I. Shift in accumulation of flavonoids and phenolic acids in lettuce attributable to changes in ultraviolet radiation and temperature. Sci. Hortic. 2018, 239, 193–204. [Google Scholar] [CrossRef]
  6. Dillard, C.J.; German, J.B. Phytochemicals: Nutraceuticals and human health. J. Sci. Food Agric. 2000, 80, 1744–1756. [Google Scholar] [CrossRef]
  7. Becker, C.; Urlić, B.; Špika, M.J.; Kläring, H.-P.; Krumbein, A.; Baldermann, S.; Ban, S.G.; Perica, S.; Schwarz, D. Nitrogen limited red and green leaf lettuce accumulate flavonoid glycosides, caffeic acid derivatives, and sucrose while losing chlorophylls, β-carotene and xanthophylls. PLoS ONE 2015, 10, e0142867. [Google Scholar] [CrossRef]
  8. Gruda, N.; Savvas, D.; Colla, G.; Rouphael, Y. Impacts of genetic material and current technologies on product quality of selected greenhouse vegetables–A review. Eur. J. Hortic. Sci. 2018, 83. [Google Scholar] [CrossRef]
  9. Hoque, M.M.; Ajwa, H.; Othman, M.; Smith, R.; Cahn, M. Yield and postharvest quality of lettuce in response to nitrogen, phosphorus, and potassium fertilizers. HortScience 2010, 45, 1539–1544. [Google Scholar] [CrossRef] [Green Version]
  10. Mampholo, B.M.; Maboko, M.; Soundy, P.; Sivakumar, D. Postharvest responses of hydroponically grown lettuce varieties to nitrogen application rate. J. Integr. Agric. 2019, 18, 2272–2283. [Google Scholar] [CrossRef]
  11. Lecompte, F.; Abro, M.; Nicot, P. Can plant sugars mediate the effect of nitrogen fertilization on lettuce susceptibility to two necrotrophic pathogens: Botrytis cinerea and Sclerotinia sclerotiorum? Plant Soil 2013, 369, 387–401. [Google Scholar] [CrossRef]
  12. Vorholt, J.A. Microbial life in the phyllosphere. Nat. Rev. Microbiol. 2012, 10, 828. [Google Scholar] [CrossRef] [PubMed]
  13. Brandl, M.; Amundson, R. Leaf age as a risk factor in contamination of lettuce with Escherichia coli O157: H7 and Salmonella enterica. Appl. Environ. Microbiol. 2008, 74, 2298–2306. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Nemadodzi, L.E.; Araya, H.; Nkomo, M.; Ngezimana, W.; Mudau, N.F. Nitrogen, phosphorus, and potassium effects on the physiology and biomass yield of baby spinach (Spinacia oleracea L.). J. Plant Nutr. 2017, 40, 2033–2044. [Google Scholar] [CrossRef]
  15. Fontes, P.; Pereira, P.; Conde, R. Critical chlorophyll, total nitrogen, and nitrate-nitrogen in leaves associated to maximum lettuce yield. J. Plant Nutr. 1997, 20, 1061–1068. [Google Scholar] [CrossRef]
  16. Boldt, J.K.; Meyer, M.H.; Erwin, J.E. Foliar anthocyanins: A horticultural review. Hortic. Rev. (Am. Soc. Hortic. Sci.) 2014, 42, 209–252. [Google Scholar]
  17. Lawanson, A.; Akindele, B.; Fasalojo, P.; Akpe, B. Time-course of anthocyanin formation during deficiencies of nitrogen, phosphorus and potassium in seedlings of Zea mays Linn. var. ES 1. Z. Pflanzenphysiol. 1972, 66, 251–253. [Google Scholar] [CrossRef]
  18. Yuan, Y.; Chiu, L.-W.; Li, L. Transcriptional regulation of anthocyanin biosynthesis in red cabbage. Planta 2009, 230, 1141–1153. [Google Scholar] [CrossRef]
  19. White, P.J.; Brown, P.H. Plant nutrition for sustainable development and global health. Ann. Bot. 2010, 105, 1073–1080. [Google Scholar] [CrossRef] [Green Version]
  20. Verzeaux, J.; Hirel, B.; Dubois, F.; Lea, P.J.; Tétu, T. Agricultural practices to improve nitrogen use efficiency through the use of arbuscular mycorrhizae: Basic and agronomic aspects. Plant Sci. 2017, 264, 48–56. [Google Scholar] [CrossRef] [Green Version]
  21. Parry, C.; Blonquist, J.M., Jr.; Bugbee, B. In situ measurement of leaf chlorophyll concentration: Analysis of the optical/absolute relationship. Plant Cell Environ. 2014, 37, 2508–2520. [Google Scholar] [CrossRef]
  22. Van den Berg, A.K.; Perkins, T.D. Nondestructive estimation of anthocyanin content in autumn sugar maple leaves. HortScience 2005, 40, 685–686. [Google Scholar] [CrossRef] [Green Version]
  23. Miller, R.O. Extractable Chloride, Nitrate, Orthophosphate, Potassium, and Sulfate-Sulfur in Plant Tissue: 2% Acetic Acid Extraction. In Handbook of Reference Methods for Plant Analysis; Kalra, Y.P., Ed.; CRC Press: Boca Raton, FL, USA, 1998; pp. 115–118. [Google Scholar]
  24. Prokopy, W.R. Phosphorus in Acetic Acid Extracts; Lachat Instruments: Milwaukee, WI, USA, 1995. [Google Scholar]
  25. Jones, B.J.J. Extraction of Chloride (Cl), Nitrate (NO3), Orthophosphate (PO4), Potassium (K), and Sulfate (SO4) from Plant Tissue Using 2% Acetic Acid. In Laboratory Guide for Conducting Soil Tests and Plant Analysis; CRC Press: Boca Raton, FL, USA, 2001; pp. 228–229. [Google Scholar]
  26. Nakagawa, S.; Cuthill, I.C. Effect size, confidence interval and statistical significance: A practical guide for biologists. Biol. Rev. 2007, 82, 591–605. [Google Scholar] [CrossRef]
  27. Poorter, H.; Garnier, E. Plant growth analysis: An evaluation of experimental design and computational methods. J. Exp. Bot. 1996, 47, 1343–1351. [Google Scholar] [CrossRef] [Green Version]
  28. Hoffmann, W.A.; Poorter, H. Avoiding bias in calculations of relative growth rate. Ann. Bot. 2002, 90, 37–42. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  29. Kerkhoff, A.J.; Enquist, B. Multiplicative by nature: Why logarithmic transformation is necessary in allometry. J. Theor. Biol. 2009, 257, 519–521. [Google Scholar] [CrossRef]
  30. Johnstone, P.; Hartz, T.; Cahn, M.; Johnstone, M. Lettuce response to phosphorus fertilization in high phosphorus soils. HortScience 2005, 40, 1499–1503. [Google Scholar] [CrossRef] [Green Version]
  31. Mampholo, B.M.; Maboko, M.M.; Soundy, P.; Sivakumar, D. Phytochemicals and overall quality of leafy lettuce (Lactuca sativa L.) varieties grown in closed hydroponic system. J. Food Qual. 2016, 39, 805–815. [Google Scholar] [CrossRef]
  32. Zhang, G.; Johkan, M.; Hohjo, M.; Tsukagoshi, S.; Maruo, T. Plant growth and photosynthesis response to low potassium conditions in three lettuce (Lactuca sativa) types. Hortic. J. 2017. [Google Scholar] [CrossRef] [Green Version]
  33. Hosokawa, M.; Yamauchi, T.; Takahama, M.; Goto, M.; Mikano, S.; Yamaguchi, Y.; Tanaka, Y.; Ohno, S.; Koeda, S.; Doi, M. Phosphorus starvation induces post-transcriptional CHS gene silencing in Petunia corolla. Plant Cell Rep. 2013, 32, 601–609. [Google Scholar] [CrossRef]
  34. Simko, I.; Hayes, R.J.; Furbank, R.T. Non-destructive phenotyping of lettuce plants in early stages of development with optical sensors. Front. Plant Sci. 2016, 7, 1985. [Google Scholar] [CrossRef] [Green Version]
  35. Gurdon, C.; Poulev, A.; Armas, I.; Satorov, S.; Tsai, M.; Raskin, I. Genetic and phytochemical characterization of lettuce flavonoid biosynthesis mutants. Sci. Rep. 2019, 9, 3305. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Mamo, B.E.; Hayes, R.J.; Truco, M.J.; Puri, K.D.; Michelmore, R.W.; Subbarao, K.V.; Simko, I. The genetics of resistance to lettuce drop (Sclerotinia spp.) in lettuce in a recombinant inbred line population from Reine des Glaces × Eruption. Theor. Appl. Genet. 2019, 132, 2439–2460. [Google Scholar] [CrossRef] [PubMed]
  37. Su, W.; Tao, R.; Liu, W.; Yu, C.; Yue, Z.; He, S.; Lavelle, D.; Zhang, W.; Zhang, L.; An, G. Characterization of four polymorphic genes controlling red leaf colour in lettuce that have undergone disruptive selection since domestication. Plant Biotechnol. J. 2020, 18, 479–490. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  38. Ercoli, L.; Mariotti, M.; Masoni, A.; Massantini, F. Relationship between nitrogen and chlorophyll content and spectral properties in maize leaves. Eur. J. Agron. 1993, 2, 113–117. [Google Scholar] [CrossRef]
  39. Makino, A.; Osmond, B. Effects of nitrogen nutrition on nitrogen partitioning between chloroplasts and mitochondria in pea and wheat. Plant Physiol. 1991, 96, 355–362. [Google Scholar] [CrossRef] [Green Version]
  40. Wang, Y.; Wang, D.; Shi, P.; Omasa, K. Estimating rice chlorophyll content and leaf nitrogen concentration with a digital still color camera under natural light. Plant Methods 2014, 10, 36. [Google Scholar] [CrossRef] [Green Version]
  41. Lavakush; Yadav, J.; Verma, J.P.; Jaiswal, D.K.; Kumar, A. Evaluation of PGPR and different concentration of phosphorus level on plant growth, yield and nutrient content of rice (Oryza sativa). Ecol. Eng. 2014, 62, 123–128. [Google Scholar] [CrossRef]
  42. Vafadar, F.; Amooaghaie, R.; Otroshy, M. Effects of plant-growth-promoting rhizobacteria and arbuscular mycorrhizal fungus on plant growth, stevioside, NPK, and chlorophyll content of Stevia rebaudiana. J. Plant Interact. 2014, 9, 128–136. [Google Scholar] [CrossRef]
Figure 1. Principal component analysis score plot (a) and loading plot (b). Data from five treatments are indicated by icons of different shapes and colors. NPK treatment (control) supplied optimal concentrations of all three macronutrients, nPK treatment supplied only 1/3 of N compared to control, PK treatment did not supply N, NK treatment did not supply P, and NP treatment did not supply K. NK and PK treatment data are generally well separated from those of the other three treatments and from each other. The loading plot results indicate a close relationship between Lb SPAD (chlorophyll content) and Lb N (nitrogen content) traits.
Figure 1. Principal component analysis score plot (a) and loading plot (b). Data from five treatments are indicated by icons of different shapes and colors. NPK treatment (control) supplied optimal concentrations of all three macronutrients, nPK treatment supplied only 1/3 of N compared to control, PK treatment did not supply N, NK treatment did not supply P, and NP treatment did not supply K. NK and PK treatment data are generally well separated from those of the other three treatments and from each other. The loading plot results indicate a close relationship between Lb SPAD (chlorophyll content) and Lb N (nitrogen content) traits.
Horticulturae 06 00015 g001
Figure 2. Values of ΔABS and 2ΔREL parameters for fresh weight (Lb FW). The ΔABS bars illustrate (Lb-fold) changes from values observed under control conditions (NPK treatment) for each accession and treatment. NPK treatment (control) supplied optimal concentrations of all three macronutrients, nPK treatment supplied only 1/3 of N compared to control, PK treatment did not supply N, NK treatment did not supply P, and NP treatment did not supply K. Significant (P ≤ 0.05) ΔABS values detected by Student’s t-test are indicated by asterisks. Solid horizontal red lines show overall mean of ΔABS at each treatment. Dashed horizontal red lines show 95% upper and lower decision limits for analysis of means (ANOM). The 2ΔREL parameter is calculated as the difference between ΔABS of an accession and the overall mean of ΔABS for all accessions tested with the particular treatment. The significance of the 2ΔREL parameter is indicated by capital letters L (significantly lower than overall mean) or H (significantly higher than overall mean).
Figure 2. Values of ΔABS and 2ΔREL parameters for fresh weight (Lb FW). The ΔABS bars illustrate (Lb-fold) changes from values observed under control conditions (NPK treatment) for each accession and treatment. NPK treatment (control) supplied optimal concentrations of all three macronutrients, nPK treatment supplied only 1/3 of N compared to control, PK treatment did not supply N, NK treatment did not supply P, and NP treatment did not supply K. Significant (P ≤ 0.05) ΔABS values detected by Student’s t-test are indicated by asterisks. Solid horizontal red lines show overall mean of ΔABS at each treatment. Dashed horizontal red lines show 95% upper and lower decision limits for analysis of means (ANOM). The 2ΔREL parameter is calculated as the difference between ΔABS of an accession and the overall mean of ΔABS for all accessions tested with the particular treatment. The significance of the 2ΔREL parameter is indicated by capital letters L (significantly lower than overall mean) or H (significantly higher than overall mean).
Horticulturae 06 00015 g002
Figure 3. Values of ΔABS and 2ΔREL parameters for chlorophyll index (Lb SPAD). The ΔABS bars illustrate (Lb-fold) changes from values observed under control conditions (NPK treatment) for each accession and treatment. NPK treatment (control) supplied optimal concentrations of all three macronutrients, nPK treatment supplied only 1/3 of N compared to control, PK treatment did not supply N, NK treatment did not supply P, and NP treatment did not supply K. Significant (p ≤ 0.05) ΔABS values detected by Student’s t-test are indicated by asterisks. Solid horizontal red lines show overall mean of ΔABS at each treatment. Dashed horizontal red lines show 95% upper and lower decision limits for ANOM. The 2ΔREL parameter is calculated as the difference between ΔABS of an accession and the overall mean of ΔABS for all accessions tested for the particular treatment. The significance of the 2ΔREL parameter is indicated by capital letters L (significantly lower than overall mean) or H (significantly higher than overall mean).
Figure 3. Values of ΔABS and 2ΔREL parameters for chlorophyll index (Lb SPAD). The ΔABS bars illustrate (Lb-fold) changes from values observed under control conditions (NPK treatment) for each accession and treatment. NPK treatment (control) supplied optimal concentrations of all three macronutrients, nPK treatment supplied only 1/3 of N compared to control, PK treatment did not supply N, NK treatment did not supply P, and NP treatment did not supply K. Significant (p ≤ 0.05) ΔABS values detected by Student’s t-test are indicated by asterisks. Solid horizontal red lines show overall mean of ΔABS at each treatment. Dashed horizontal red lines show 95% upper and lower decision limits for ANOM. The 2ΔREL parameter is calculated as the difference between ΔABS of an accession and the overall mean of ΔABS for all accessions tested for the particular treatment. The significance of the 2ΔREL parameter is indicated by capital letters L (significantly lower than overall mean) or H (significantly higher than overall mean).
Horticulturae 06 00015 g003
Figure 4. Values of ΔABS and 2ΔREL parameters for anthocyanin index (Lb ACI). The ΔABS bars illustrate (Lb-fold) changes from values observed under control conditions (NPK treatment) for each accession and treatment. NPK treatment (control) supplied optimal concentrations of all three macronutrients, nPK treatment supplied only 1/3 of N compared to control, PK treatment did not supply N, NK treatment did not supply P, and NP treatment did not supply K. Significant (p ≤ 0.05) ΔABS values detected by Student’s t-test are indicated by asterisks. Solid horizontal red lines show overall mean of ΔABS at each treatment. Dashed horizontal red lines show 95% upper and lower decision limits for ANOM. The 2ΔREL parameter is calculated as the difference between ΔABS of an accession and the overall mean of ΔABS for all accessions tested at the particular treatment. The significance of the 2ΔREL parameter is indicated by capital letters L (significantly lower than overall mean) or H (significantly higher than overall mean).
Figure 4. Values of ΔABS and 2ΔREL parameters for anthocyanin index (Lb ACI). The ΔABS bars illustrate (Lb-fold) changes from values observed under control conditions (NPK treatment) for each accession and treatment. NPK treatment (control) supplied optimal concentrations of all three macronutrients, nPK treatment supplied only 1/3 of N compared to control, PK treatment did not supply N, NK treatment did not supply P, and NP treatment did not supply K. Significant (p ≤ 0.05) ΔABS values detected by Student’s t-test are indicated by asterisks. Solid horizontal red lines show overall mean of ΔABS at each treatment. Dashed horizontal red lines show 95% upper and lower decision limits for ANOM. The 2ΔREL parameter is calculated as the difference between ΔABS of an accession and the overall mean of ΔABS for all accessions tested at the particular treatment. The significance of the 2ΔREL parameter is indicated by capital letters L (significantly lower than overall mean) or H (significantly higher than overall mean).
Horticulturae 06 00015 g004
Figure 5. Values of ΔABS and 2ΔREL parameters for nitrogen content (Lb N). The ΔABS bars illustrate (Lb-fold) changes from values observed under control conditions (NPK treatment) for each accession and treatment. NPK treatment (control) supplied optimal concentrations of all three macronutrients, nPK treatment supplied only 1/3 of N compared to control, PK treatment did not supply N, NK treatment did not supply P, and NP treatment did not supply K. Significant (p ≤ 0.05) ΔABS values detected by Student’s t-test are indicated by asterisks. Solid horizontal red lines show overall mean of ΔABS at each treatment. Dashed horizontal red lines show 95% upper and lower decision limits for ANOM. The 2ΔREL parameter is calculated as the difference between ΔABS of an accession and the overall mean of ΔABS for all accessions tested for the particular treatment. The significance of the 2ΔREL parameter is indicated by capital letters L (significantly lower than overall mean) or H (significantly higher than overall mean).
Figure 5. Values of ΔABS and 2ΔREL parameters for nitrogen content (Lb N). The ΔABS bars illustrate (Lb-fold) changes from values observed under control conditions (NPK treatment) for each accession and treatment. NPK treatment (control) supplied optimal concentrations of all three macronutrients, nPK treatment supplied only 1/3 of N compared to control, PK treatment did not supply N, NK treatment did not supply P, and NP treatment did not supply K. Significant (p ≤ 0.05) ΔABS values detected by Student’s t-test are indicated by asterisks. Solid horizontal red lines show overall mean of ΔABS at each treatment. Dashed horizontal red lines show 95% upper and lower decision limits for ANOM. The 2ΔREL parameter is calculated as the difference between ΔABS of an accession and the overall mean of ΔABS for all accessions tested for the particular treatment. The significance of the 2ΔREL parameter is indicated by capital letters L (significantly lower than overall mean) or H (significantly higher than overall mean).
Horticulturae 06 00015 g005
Figure 6. Values of ΔABS and 2ΔREL parameters for phosphorus content (Lb P). The ΔABS bars illustrate (Lb-fold) changes from values observed with control (NPK treatment) for each accession and treatment. NPK treatment (control) supplied optimal concentrations of all three macronutrients, nPK treatment supplied only 1/3 of N compared to control, PK treatment did not supply N, NK treatment did not supply P, and NP treatment did not supply K. Significant (p ≤ 0.05) ΔABS values detected by Student’s t-test are indicated by asterisks. Solid horizontal red lines show the overall mean of ΔABS with each treatment. Dashed horizontal red lines show 95% upper and lower decision limits for ANOM. The 2ΔREL parameter is calculated as the difference between ΔABS of an accession and the overall mean of ΔABS for all accessions tested with the particular treatment. The significance of the 2ΔREL parameter is indicated by capital letters L (significantly lower than overall mean) or H (significantly higher than overall mean).
Figure 6. Values of ΔABS and 2ΔREL parameters for phosphorus content (Lb P). The ΔABS bars illustrate (Lb-fold) changes from values observed with control (NPK treatment) for each accession and treatment. NPK treatment (control) supplied optimal concentrations of all three macronutrients, nPK treatment supplied only 1/3 of N compared to control, PK treatment did not supply N, NK treatment did not supply P, and NP treatment did not supply K. Significant (p ≤ 0.05) ΔABS values detected by Student’s t-test are indicated by asterisks. Solid horizontal red lines show the overall mean of ΔABS with each treatment. Dashed horizontal red lines show 95% upper and lower decision limits for ANOM. The 2ΔREL parameter is calculated as the difference between ΔABS of an accession and the overall mean of ΔABS for all accessions tested with the particular treatment. The significance of the 2ΔREL parameter is indicated by capital letters L (significantly lower than overall mean) or H (significantly higher than overall mean).
Horticulturae 06 00015 g006
Figure 7. Values of ΔABS and 2ΔREL parameters for potassium content (Lb K). The ΔABS bars illustrate (Lb-fold) changes from values observed under control conditions (NPK treatment) for each accession and treatment. NPK treatment (control) supplied optimal concentrations of all three macronutrients, nPK treatment supplied only 1/3 of N compared to control, PK treatment did not supply N, NK treatment did not supply P, and NP treatment did not supply K. Significant (p ≤ 0.05) ΔABS values detected by Student’s t-test are indicated by asterisks. Solid horizontal red lines show overall mean values of ΔABS with each treatment. Dashed horizontal red lines show 95% upper and lower decision limits for ANOM. The 2ΔREL parameter is calculated as the difference between ΔABS of an accession and the overall mean of ΔABS for all accessions tested with the particular treatment. The significance of the 2ΔREL parameter is indicated by capital letters L (significantly lower than overall mean) or H (significantly higher than overall mean). The missing dashed line for NP is outside of the range indicated on the y-axis.
Figure 7. Values of ΔABS and 2ΔREL parameters for potassium content (Lb K). The ΔABS bars illustrate (Lb-fold) changes from values observed under control conditions (NPK treatment) for each accession and treatment. NPK treatment (control) supplied optimal concentrations of all three macronutrients, nPK treatment supplied only 1/3 of N compared to control, PK treatment did not supply N, NK treatment did not supply P, and NP treatment did not supply K. Significant (p ≤ 0.05) ΔABS values detected by Student’s t-test are indicated by asterisks. Solid horizontal red lines show overall mean values of ΔABS with each treatment. Dashed horizontal red lines show 95% upper and lower decision limits for ANOM. The 2ΔREL parameter is calculated as the difference between ΔABS of an accession and the overall mean of ΔABS for all accessions tested with the particular treatment. The significance of the 2ΔREL parameter is indicated by capital letters L (significantly lower than overall mean) or H (significantly higher than overall mean). The missing dashed line for NP is outside of the range indicated on the y-axis.
Horticulturae 06 00015 g007
Figure 8. Scatterplot matrix of six traits evaluated in eight accessions in five treatments. Relationships between fresh weight (Lb FW) and content of chlorophyll (Lb SPAD), anthocyanin (Lb ACI), nitrogen (Lb N), phosphorus (Lb P), and potassium (Lb K) are shown in individual windows. Data from five treatments are indicated by icons with different shapes and colors. NPK treatment (control) supplied optimal concentrations of all three macronutrients, nPK treatment supplied only 1/3 of N compared to control, PK treatment did not supply N, NK treatment did not supply P, and NP treatment did not supply K. Note that regressions between traits calculated on data from different fertilizer treatments may have dissimilar slopes or intercepts but similar correlation coefficients (Table 4).
Figure 8. Scatterplot matrix of six traits evaluated in eight accessions in five treatments. Relationships between fresh weight (Lb FW) and content of chlorophyll (Lb SPAD), anthocyanin (Lb ACI), nitrogen (Lb N), phosphorus (Lb P), and potassium (Lb K) are shown in individual windows. Data from five treatments are indicated by icons with different shapes and colors. NPK treatment (control) supplied optimal concentrations of all three macronutrients, nPK treatment supplied only 1/3 of N compared to control, PK treatment did not supply N, NK treatment did not supply P, and NP treatment did not supply K. Note that regressions between traits calculated on data from different fertilizer treatments may have dissimilar slopes or intercepts but similar correlation coefficients (Table 4).
Horticulturae 06 00015 g008
Figure 9. Correlation between the content of nitrogen (Lb N) and phosphorus (Lb P) in baby leaf lettuce tested with different fertilizer treatments. Strong, positive, linear correlations ranging from r = 0.74 to r = 0.90 (p < 0.01, Table 4) were detected between Lb N and Lb P with each of the five treatments when analyzed separately. The correlation coefficient dropped to r = −0.02 (black, dashed line) when data from all treatments were considered together. Five treatments are indicated by icons with different shapes and colors. NPK treatment (control) supplied optimal concentrations of all three macronutrients, nPK treatment supplied only 1/3 of N compared to control, PK treatment did not supply N, NK treatment did not supply P, and NP treatment did not supply K.
Figure 9. Correlation between the content of nitrogen (Lb N) and phosphorus (Lb P) in baby leaf lettuce tested with different fertilizer treatments. Strong, positive, linear correlations ranging from r = 0.74 to r = 0.90 (p < 0.01, Table 4) were detected between Lb N and Lb P with each of the five treatments when analyzed separately. The correlation coefficient dropped to r = −0.02 (black, dashed line) when data from all treatments were considered together. Five treatments are indicated by icons with different shapes and colors. NPK treatment (control) supplied optimal concentrations of all three macronutrients, nPK treatment supplied only 1/3 of N compared to control, PK treatment did not supply N, NK treatment did not supply P, and NP treatment did not supply K.
Horticulturae 06 00015 g009
Table 1. Treatments and combinations of fertilizers.
Table 1. Treatments and combinations of fertilizers.
TreatmentAmmonium Nitrate z (g)Phosphate y (g)Potash x (g)Micronutrients w (mL)
NPK (control)632.33
nPK232.33
PK032.33
NK602.33
NP6303
z Non-Coated Ammonium Nitrate 34-0-0 Prill Form; y Triple Super Phosphate 0-46-0 Easy Peasy Plants 99% pure; x All-Natural Muriate of Potash-Easy Peasy 0-0-60 Potassium; w 1 g of Axilo Mix 5 (0-0-0) Micronutrient Mix in 1L of distilled water.
Table 2. Main effect of factors, calculated from binary logarithm (Lb) values, on fresh weight (Lb FW) and contents of chlorophyll (Lb SPAD), anthocyanin (Lb ACI), nitrogen (Lb N), phosphorus (Lb P), and potassium (Lb K).
Table 2. Main effect of factors, calculated from binary logarithm (Lb) values, on fresh weight (Lb FW) and contents of chlorophyll (Lb SPAD), anthocyanin (Lb ACI), nitrogen (Lb N), phosphorus (Lb P), and potassium (Lb K).
DescriptorLb FWLb SPADLb ACILb NLb PLb K
Effect Sizeω 2 z
Model0.937 ***0.726 ***0.908 ***0.894 ***0.922 ***0.847 ***
Accession0.076 ***0.377 ***0.725 ***0.108 ***0.031 ***0.320 ***
Treatment0.852 ***0.316 ***0.037 ***0.759 ***0.892 ***0.435 ***
Interaction0.009 *0.032 *0.146 ***0.027 ***0.0000.092 ***
AccessionsMean y
Annapolis3.68 b5.08 ab5.48 a1.88 bc−1.21 ab1.79 b
B. Barrage3.88 ab4.83 b2.30 f1.54 d−1.65 c1.47 c
Darkland4.16 a5.01 b2.73 de1.50 d−1.49 bc1.57 c
Eruption4.00 ab4.92 b4.26 b1.84 c−1.16 a2.24 a
La Brillante3.86 ab4.40 c1.75 g2.07 ab−1.09 a2.14 a
Merlot3.82 ab4.32 c3.73 c1.95 bc−1.03 a2.11 a
Salinas3.75 b4.93 b2.39 ef1.94 bc−1.26 ab2.12 a
UC96US232.63 c5.28 a2.96 d2.18 a−1.17 a2.16 a
TreatmentsMean
NPK5.03 a5.02 a3.06 bc2.16 b−0.60 a2.22 a
nPK4.53 b4.94 ab3.07 bc1.80 c−0.71 a2.28 a
PK1.60 d4.32 c2.94 c0.78 d−1.26 b1.95 b
NK2.42 c4.85 b3.70 a2.36 a−3.18 c1.94 b
NP5.02 a5.11 a3.23 b2.21 b−0.53 a1.36 c
z Effect size (ω2) and significance of factors were calculated from two-way ANOVA. Asterisks indicate value significance at p ≤ 0.05 (*), p ≤ 0.01 (**), and p ≤ 0.001 (***). y Values for accessions are means of 20 plants, while values for treatments are means of 32 plants. Mean values within a column and a factor followed by different letters are significantly different at p ≤ 0.05 using the Tukey’s HSD test.
Table 3. Combined effect of factors on fresh weight (Lb FW), content of chlorophyll (Lb SPAD), anthocyanin (Lb ACI), nitrogen (Lb N), phosphorus (Lb P), and potassium (Lb K).
Table 3. Combined effect of factors on fresh weight (Lb FW), content of chlorophyll (Lb SPAD), anthocyanin (Lb ACI), nitrogen (Lb N), phosphorus (Lb P), and potassium (Lb K).
AccessionTreatment zLb FWLb SPADLb ACILb NLb PLb K
AnnapolisNPK4.97 a–ey5.28 a–f4.69 b–e2.17 b–g−0.62 a–d2.12 a–i
B. BarrageNPK5.37 a–d5.15 a–g2.62 j–o1.59 g–j−1.11 b–f1.76 g–m
DarklandNPK5.93 a5.31 a–e3.07 h–l1.77 e–i−0.74 a–f1.87 e–m
EruptionNPK5.13 a–d5.06 a–g3.63 e–j2.21 a–f−0.52 a–c2.21 a–h
La BrillanteNPK5.02 a–d4.72 c–j1.93 m–p2.54 a–c−0.33 ab2.57 a
MerlotNPK4.86 a–f4.28 h–l2.80 i–n2.37 a–e−0.21 a2.43 a–d
SalinasNPK5.09 a–d5.01 a–h2.45 k–p2.22 a–f−0.63 a–d2.41 a–e
UC96US23NPK3.90 e–h5.39 a–d3.29 g–k2.38 a–d−0.61 a–d2.36 a–f
AnnapolisnPK4.43 c–g5.13 a–g4.96 b–d1.86 d–i−0.60 a–d1.97 c–l
B. BarragenPK4.27 d–g5.01 a–h2.44 k–p1.55 h–j−1.02 a–f1.84 f–m
DarklandnPK4.69 b–g5.16 a–g2.61 j–o1.36 i–k−0.90 a–f2.20 a–h
EruptionnPK4.92 a–f5.05 a–g4.00 d–h1.72 f–i−0.69 a–e2.55 ab
La BrillantenPK4.68 b–g4.46 g–l1.69 n–p1.96 c–h−0.59 a–d2.50 a–c
MerlotnPK4.66 b–g4.22 i–l3.13 h–l1.98 c–h−0.53 a–c2.42 a–d
SalinasnPK4.72 b–g5.09 a–g2.61 j–o1.82 d–i−0.79 a–f2.42 a–d
UC96US23nPK3.83 f–i5.44 a–c3.16 h–l2.15 b–h−0.58 a–d2.33 a–f
AnnapolisPK1.43 m–n4.55 f–k5.78 b1.11 j–l−1.06 a–f2.14 a–i
B. BarragePK2.10 k–m4.15 j–l1.50 o–p0.39 m–n−1.57 f1.41 m–o
DarklandPK1.85 k–n4.13 j–l1.72 n–p0.28 n−1.49 ef1.17 n–p
EruptionPK1.78 l–n4.46 g–l4.59 c–f0.88 k–m−1.06 a–f2.60 a
La BrillantePK1.40 m–n3.95 k–l1.38 p0.88 k–n−1.21 c–f1.92 d–m
MerlotPK1.71 l–n3.82 l4.42 c–g0.91 k–m−1.08 b–f2.30 a–g
SalinasPK1.74 l–n4.66 d–k2.04 l–p0.75 l–n−1.39 d–f1.90 d–m
UC96US23PK0.80 n4.81 b–j2.11 l–p1.07 j–l−1.26 c–f2.19 a–h
AnnapolisNK2.31 k–m5.09 a–g7.08 a2.21 a–f−3.22 g1.74 h–m
B. BarrageNK2.54 k–l4.60 e–k2.04 l–p2.22 a–f−3.65 g1.48 l–o
DarklandNK2.77 i–l4.98 a–h2.76 j–n2.26 a–f−3.56 g1.51 k–n
EruptionNK2.88 h–k4.95 a–h5.26 bc2.21 a–f−3.03 g2.33 a–f
La BrillanteNK2.73 j–l4.45 g–l1.92 m–p2.73 ab−2.89 g2.06 a–j
MerlotNK2.91 h–k4.58 e–k5.22 bc2.12 c–h−3.00 g1.94 d–m
SalinasNK2.37 k–m4.98 a–h2.48 j–p2.36 a–e−3.26 g2.03 b–k
UC96US23NK0.89 n5.16 a–g2.81 i–n2.78 a−2.85 g2.40 a–e
AnnapolisNP5.27 a–d5.37 a–d4.88 b–d2.04 c–h−0.55 a–d0.96 op
B. BarrageNP5.12 a–d5.27 a–f2.89 h–m1.96 c–i−0.91 a–f0.84 p
DarklandNP5.53 ab5.49 ab3.49 f–k1.86 d–i−0.74 a–f1.10 n–p
EruptionNP5.31 a–d5.10 a–g3.84 d–i2.19 a–g−0.50 a–c1.50 k–o
La BrillanteNP5.44 a–c4.44 g–l1.83 m–p2.23 a–f−0.44 a–c1.63 i–n
MerlotNP4.94 a–e4.70 d–j3.10 h–l2.38 a–d−0.31 ab1.47 l–o
SalinasNP4.82 b–f4.91 a–i2.38 k–p2.55 a–c−0.23 a1.86 f–m
UC96US23NP3.73 g–j5.58 a3.45 f–k2.50 a–c−0.56 a–d1.54 j–n
z NPK treatment (control) supplied optimal concentrations of all three macronutrients; nPK treatment supplied only 1/3 of N compared to control; PK treatment did not supply N; NK treatment did not supply P; and NP treatment did not supply K. y Values are means of four plants (replicates). Mean values within a column followed by different letters are significantly different at p ≤ 0.05 using the Tukey’s HSD test.
Table 4. Pearson correlation coefficients between fresh weight (Lb FW) and content of chlorophyll (Lb SPAD), anthocyanin (Lb ACI), nitrogen (Lb N), phosphorus (Lb P), and potassium (Lb K), evaluated from diverse data sets
Table 4. Pearson correlation coefficients between fresh weight (Lb FW) and content of chlorophyll (Lb SPAD), anthocyanin (Lb ACI), nitrogen (Lb N), phosphorus (Lb P), and potassium (Lb K), evaluated from diverse data sets
Treatmentz Accessiony
Trait 1Trait 2AllNPKnPKPKNKNPAnnapolisB. BarrageDarklandEruptionLa BrillanteMerlotSalinasUC96US23
Lb FWLb SPAD0.41 **x0.08−0.25−0.15−0.32−0.290.81 **0.82 **0.78 **0.58 **0.57 **0.45 *0.440.79 **
Lb FWLb ACI−0.03−0.05−0.09−0.06−0.020.14−0.71 **0.81 **0.75 **−0.72 **0.44−0.62 **0.370.79 **
Lb FW Lb N0.40 **−0.51 **−0.31−0.48 **−0.54 **−0.43*0.51 *0.390.440.62 **0.55 *0.78 **0.56 *0.35
Lb FWLb P0.56 **−0.23−0.19−0.17−0.17−0.200.55 *0.62 **0.58 **0.54 *0.61 **0.56 *0.68 **0.77 **
Lb FWLb K−0.11−0.46 **0.09−0.32−0.30−0.33−0.43−0.050.29−0.52 *0.21−0.140.43−0.23
Lb SPADLb ACI0.30 **0.54 **0.340.100.69 **0.31−0.220.92 **0.82 **−0.140.80 **0.300.69 **0.85 **
Lb SPADLb N0.41 **−0.28−0.180.17−0.200.050.65 **0.55 *0.71 **0.61 **0.67 **0.54 *0.380.54 *
Lb SPADLb P0.09−0.36 *−0.21−0.13−0.290.040.110.420.240.140.21−0.16−0.030.38
Lb SPADLb K−0.14−0.35 *−0.270.12−0.41 *0.20−0.50 *−0.080.25−0.260.42−0.55 *0.26−0.40
Lb ACILb N0.11−0.010.040.50 **−0.23−0.47 **−0.020.59 **0.75 **−0.120.75 **−0.390.320.69 **
Lb ACILb P−0.110.010.110.49 **−0.110.09−0.86 **0.420.21−0.75 **0.01−0.76 **−0.100.40
Lb ACILb K0.04−0.18−0.200.62 **−0.42 *0.040.08−0.170.080.370.26−0.150.12−0.30
Lb NLb P−0.020.90 **0.87 **0.83 **0.87 **0.74 **−0.11−0.25−0.26−0.05−0.120.210.08−0.13
Lb NLb K0.19 *0.97 **0.62 **0.86 **0.87 **0.65 **−0.16−0.010.26−0.390.33−0.180.270.09
Lb PLb K0.18 *0.88 **0.72 **0.77 **0.81 **0.88 **0.050.050.24−0.130.300.230.29−0.22
z NPK treatment (control) supplied optimal concentrations of all three macronutrients, nPK treatment supplied only 1/3 of N compared to control, PK treatment did not supply N, NK treatment did not supply P, and NP treatment did not supply K. Correlation coefficients were calculated on data from individual treatments across all accessions (n = 32), except “all”, where correlation coefficients were calculated on complete datasets that included all data (n = 160). y Correlation coefficients calculated on data from individual accessions across all treatments (n = 20). x Asterisks indicate correlation coefficients significant at p ≤ 0.05 (*) or p ≤ 0.01 (**).2.

Share and Cite

MDPI and ACS Style

Simko, I. Genetic Variation in Response to N, P, or K Deprivation in Baby Leaf Lettuce. Horticulturae 2020, 6, 15. https://doi.org/10.3390/horticulturae6010015

AMA Style

Simko I. Genetic Variation in Response to N, P, or K Deprivation in Baby Leaf Lettuce. Horticulturae. 2020; 6(1):15. https://doi.org/10.3390/horticulturae6010015

Chicago/Turabian Style

Simko, Ivan. 2020. "Genetic Variation in Response to N, P, or K Deprivation in Baby Leaf Lettuce" Horticulturae 6, no. 1: 15. https://doi.org/10.3390/horticulturae6010015

APA Style

Simko, I. (2020). Genetic Variation in Response to N, P, or K Deprivation in Baby Leaf Lettuce. Horticulturae, 6(1), 15. https://doi.org/10.3390/horticulturae6010015

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

Article Metrics

Back to TopTop