Experimental design and field trial setup.
This 2-year study was conducted on certified organic land at the University of Florida Plant Science Research and Education Unit in Citra, FL, during the 2018–19 and 2019–20 production seasons. The soil type at the research site is Gainesville loamy sand (hyperthermic, coated Typic Quartzipsamments) with 97.0% sand, 2.2% clay, and 0.8% silt, and an average soil organic matter content at 0.8%. In both seasons, granular organic fertilizer Nature Safe 10N–0.9P–6.6K (Darling Ingredients Inc., Irving, TX, USA) was applied preplant at an N rate of 84.1 kg⋅ha−1. False raised beds were formed first and then followed by granular organic fertilizer application on the bed top, and tillage was used to mix organic fertilizer with soil at a depth of 15 cm on 3 Oct 2018 and 30 Sep 2019. On the same day, the final raised beds were made and covered with 0.03 mm black totally impermeable film plastic mulch (Intergro, Inc., Clearwater, FL, USA) and a single drip tape (30.5 cm emitter spacing, 3.4 L⋅h−1⋅m−1 flow rate; Jain Irrigation, Inc., Haines City, FL, USA) was put on the soil surface in the middle of each bed.
In both years, a split-plot design with four replications of each biostimulant treatment was used for the field trials. In the 2018–19 season, biostimulant treatments were arranged in the whole plots following a randomized complete block design. Three biostimulant treatments were assessed along with a no-biostimulant (water) control, including seaweed extract biostimulant Stimplex® (4.7 L⋅ha−1 biweekly application rate; Acadian Seaplants Limited, Dartmouth, NS, Canada), microbial biostimulant TerraGrow® (1.1 kg⋅ha−1 biweekly application rate; BioSafe Systems, LLC, East Hartford, CT, USA), and a combination of Stimplex® plus TerraGrow®. According to the product label, the seaweed extract biostimulant is derived from Ascophyllum nodosum; the microbial biostimulant contained Bacillus subtilis, Bacillus pumilus, Bacillus licheniformis, Bacillus megaterium, Bacillus amyloliquefaciens, and Trichoderma harzianum. The subplots consisted of two strawberry cultivars: Sweet Sensation® Florida127 and Florida Brilliance. In the 2019–20 season, the Sweet Sensation® Florida127 and Florida Brilliance strawberry cultivar treatments were arranged in the whole plots according to a randomized complete block design. The subplots consisted of the no-biostimulant control and the same biostimulant treatments evaluated in the 2018–19 season.
In both years, the experimental unit consisted of 40 plants in each subplot. Strawberry plants were grown in a staggered pattern with two rows (30.5-cm spacing between two rows) in each bed at 30.5 cm spacing between plants, and bed centers were spaced at 1.5 m. Containerized strawberry plants were transplanted on 4 Oct 2018 and 9 Oct 2019. Before strawberry transplanting, all seedling roots were dipped in water, TerraGrow® solution (7.5 g⋅L−1), Stimplex® solution (1 mL⋅L−1), or a combination solution of TerraGrow® and Stimplex® for 2 min, respectively. After transplanting, the strawberry plants in each plot receiving the microbial biostimulant treatment were immediately treated again with TerraGrow® solution at a rate of 1.7 kg⋅ha−1 (dissolved in 7014 L⋅ha−1 water to thoroughly drench the root zone with the product) through drip irrigation lines. The biostimulant combination treatment was drip applied on a biweekly basis from 18 Oct 2018 to 19 Apr 2019 and from 24 Oct 2019 to 13 Mar 2020. Liquid fish fertilizer Aqua PowerTM 5N–0.4P–0.8K (JH Biotech, Inc., Ventura, CA, USA) in combination with Big-KTM sulfate of potash 0N–0P–41.5K (JH Biotech, Inc.) were used for in-season fertigation through drip irrigation lines for the strawberry production season. In-season fertigation was applied weekly at 4.71 kg⋅ha−1 N and 4.71 kg⋅ha−1 K for 1 to 3 weeks after transplanting (WAT) in both seasons and then increased to 7.85 kg⋅ha−1 N and 6.28 kg⋅ha−1 K for 4 to 29 WAT in the 2018–19 season and for 4 to 22 WAT in the 2019–20 season. Total organic fertilizer application rates were 302.3 and 247.4 kg⋅ha−1 N, 232.9 and 189.0 kg⋅ha−1 K, 25.0 and 20.6 kg⋅ha−1 P in 2018–19 and 2019–20 seasons, respectively.
All strawberry plants were covered using the AgroFabric® System frost protection fabric (Agrifabrics LLC, Alpharetta, GA, USA) for 7 nights in the 2018–19 season and 14 nights in the 2019–20 season when the night air temperature was projected to be below 4 °C.
Biochemical analyses of seaweed extract biostimulant Stimplex®.
The presence and abundance of proteins, lipids, and metabolites in the seaweed extract were qualitatively analyzed using liquid chromatography–tandem mass spectrometry (LC-MS/MS). Biostimulant samples (100 µL each) were extracted for proteins, lipids, and metabolites. The proteins were extracted by acetone/methanol precipitation, lipids were extracted with a modified Folch method (Folch et al. 1957), and metabolites by ice-cold methanol. The details of the biochemical analysis protocols were described by Li et al. (2021).
Photosynthetic parameters and vegetative growth.
Photosynthetic parameters including photosynthetic rate (A, μmol⋅m−2⋅s−1 CO2), transpiration rate (E, mmol⋅m−2⋅s−1 H2O), intercellular CO2 concentration (ci, µmol⋅mol−1), and stomatal conductance (gsw, mol⋅m−2⋅s−1) were measured twice (early harvest and peak harvest) in each season. Three representative plants were randomly chosen from each subplot and measurements were conducted on the most recent, fully developed, mature leaf from the top of the canopy of each plant with an LI-6800 portable photosynthesis system (Software Version 1.2.; LI-COR Inc., Lincoln, NE, USA). The light intensity and the CO2 reference concentration in the leaf chamber (6 cm2 aperture) were set at 1000 μmol⋅m−2⋅s−1 and 400 μmol⋅mol−1, respectively.
The leaf relative chlorophyll content was measured at 1, 30, 60, 90, 120, 150, and 180 (first season only) d after transplanting (DAT) on fully developed and most recent mature leaves on the top of plants using a SPAD-502 Plus chlorophyl meter (Konica Minolta, Inc., Tokyo, Japan) in both seasons. Nine SPAD readings were randomly taken from the trifoliate leaves of three different plants in each subplot and then nine readings were averaged and used for the SPAD value of each subplot. Crown diameter at the base of the plant (including all crowns) was measured at 1, 30, 60, 90, 120, 150, and 180 (first season only) DAT using a caliper in both seasons. Destructive sampling was conducted at 50 and 200 DAT in the 2018–19 season and at 60 and 155 DAT during the 2019–20 season for measuring crown number, leaf number, leaf area, and aboveground biomass. The aboveground tissues of two representative plants were removed from each subplot and plants were separated into two parts: crowns (with petioles included) and leaves (fully expanded leaf blades). The numbers of crowns and fully expanded trifoliate leaves were counted and recorded, respectively. The total area of fully expanded trifoliate leaves of two strawberry plants were measured using an LI-3100 Area Meter (LI-COR Inc.). The average leaf area per leaf was determined based on the total leaf area and the total number of leaves.
Fruit yield evaluation.
Strawberry harvests started on 4 Dec 2018 (61 DAT) and 6 Dec 2019 (58 DAT) and ended on 25 Apr 2019 (203 DAT) and 16 Mar 2020 (159 DAT). Ripe fruit with 75% to 100% red color from all plants were harvested twice per week in each experimental unit (subplot) and then classified into marketable and unmarketable categories following the US Department of Agriculture standards for grades of strawberries (Santos et al. 2009; USDA-AMS 2006). The small (less than 10 g) and deformed, diseased, and pest damaged fruit were considered to be unmarketable fruit. Fresh marketable and unmarketable fruit were counted and weighed, and data were reported on a per-plant basis. Total fruit yield included marketable and unmarketable fruit yield.
N accumulation in aboveground biomass and fruit mineral content.
After measuring growth parameters as aforementioned, samples of crown and leaf parts collected from destructive sampling at 200 DAT in the 2018–19 season and at 155 DAT during the 2019–20 season were dried at 65 °C for 7 d to determine dry weight. Dried plant biomass was then ground through a 1-mm sieve using a Wiley laboratory mill (Model 4; Thomas Scientific, Swedesboro, NJ, USA), before plant tissue N content analysis (by Waters Agricultural Laboratories, Inc., Camilla, GA, USA). Ten ripe fruit were collected in the early (20 Dec 2018 and 19 Dec 2019), peak (21 Feb 2019 and 13 Feb 2020), and late (4 Apr 2019 and 2 Mar 2020) harvest seasons for determination of dry matter and N contents. The drying procedure and N content analysis for fruit followed the same protocol of plant biomass measurement. In addition, the fruit collected in the late season (4 Apr 2019 and 2 Mar 2020) were analyzed for the contents of macronutrients P, K, Ca, Mg, and S and micronutrients B, Fe, Mn, and Zn using the ICAP-open vessel wet digestion Digi Block 3000 (by Waters Agricultural Laboratories, Inc.). The N accumulation in crown and leaf parts was estimated by multiplying the plant tissue N content with the tissue dry weight per plant (g/plant). The N accumulation in fruit during early, peak, and late harvest seasons was estimated by multiplying the corresponding fruit N content of the fruit sample with the total fruit dry weight for the harvest period (g/plant). The total N accumulation (kg/ha) was then calculated by multiplying the amount of N accumulated in crown, leaf, and fruit tissues per plant by the plant population (43,056 plants/ha).
Root growth parameters.
Root architecture was also compared between the water control and the biostimulant combination treatment. Two representative strawberry plants from each of the corresponding subplots were sampled for root assessment at 201 DAT on 23 Apr 2019 and at 155 DAT on 12 Mar 2020. The whole root system was collected from the field using a shovel. The samples were submerged in plastic bags filled with deionized water for 1 h. The root system was carefully washed to remove the soil particles and then scanned using a root scanning apparatus (EPSON color image scanner LA1600+; EPSON, Toronto, ON, Canada). The determination of the root system characteristics, including total root length (cm/plant), average root diameter (mm), total root surface area (cm2/plant), and numbers of tips, forks, and crossings was performed using the WinRhizo image analysis system (Version 2012b; Regent Instruments, Quebec, QC, Canada). Then the strawberry roots were dried at 65 °C for 7 d to determine the dry weight.
Soil respiration.
Soil respiration measurements were conducted at 60, 90, and 130 DAT during the 2019–20 strawberry growing season. One cylindrical PVC collar (20 cm in diameter and 10 cm in height) was inserted into the soil to a depth of 2 cm between the strawberry plants in each subplot. Soil collars were installed 1 week before the CO2 measurements to get stabilized CO2 flux inside collars. Soil CO2 efflux was measured using an LI-6800 portable photosynthesis system with an LI-6800 soil CO2 flux chamber (LI-COR Inc.). Measurements were consecutively replicated three times at each collar and the duration of each measurement was 160 s. All weeds inside the collar were removed immediately once the collar was placed in the field and the inside of the collars was kept weed free during the strawberry growing season. Soil moisture content and temperature as useful ancillary data were measured with a Stevens HydraProbe (Stevens Water Monitoring Systems, Inc., Portland, OR, USA) inserted into the soil in the vicinity of the collar. The soil respiration rate was expressed as μmol⋅m−2⋅s−1 CO2.
Fruit quality attributes.
Fruit quality attributes including titratable acidity (TA), soluble solids content (SSC), and total anthocyanin content (TAC) were measured during the early, peak, and late harvest periods in each production season. After removal of the calyxes, 20 fully ripe strawberries per biostimulant treatment replicate with marketable fruit quality were homogenized in a blender (Model: HBB908; Hamilton Beach Brands, Inc., Glen Allen, VA, USA). A portion (∼45 g) of the homogenized mixture was centrifuged at 4 °C for 20 min at 19,319 gn using a Sorvall LYNX 4000 refrigerated centrifuge (Thermo Fisher Scientific, Waltham, MA, USA). The supernatant from the centrifuged sample was filtered through double-layer cheesecloth before further analysis.
The filtered supernatant clear juice was analyzed to measure TA and SSC. The determination of TA was conducted using an automated titrator (905 Titrando; Metrohm, Riverview, FL, USA) by titrating 6 mL strawberry juice (filtered supernatant) with a solution of 0.1 mol⋅L−1 sodium hydroxide (NaOH) to an endpoint of pH 8.2. The TA was expressed as a percentage based on citric acid equivalents. The SSC was assessed using an automatic, temperature compensated refractometer (r2i300; AMETEK Reichert Technology, Munich, Germany) and expressed as °Brix.
TAC was measured based on the pH differentiation method (Giusti and Wrolstad 2001). Briefly, anthocyanins were extracted using acidified methanol solution. A 2.0-g sample of the blended mixture was added into 5 mL of HCl (1%)-methanol solution and shaken by a vortex for 20 s. The samples covered with aluminum foil were equilibrated at 4 °C for 15 min and then centrifuged at 4 °C for 20 min at 19,319 gn using a Sorvall LYNX 4000 refrigerated centrifuge to collect the supernatant. The extractions were conducted in duplicate. Two dilutions of the same extractant sample were prepared by adding 600 μL of extract to 2.4 mL of potassium chloride (0.025 M, pH = 1.0) and to 2.4 mL of sodium acetate (0.4 M, pH = 4.5), respectively. All samples were vortexed for 10 s and then equilibrated in the dark at room temperature for 15 min. Their absorbance was measured in triplicate (200 μL for each) at 515 nm and 700 nm for each solution mixture at pH = 1.0 and pH = 4.5, respectively, using a microplate reader (Epoch 2; BioTek Instruments, Inc., Winooski, VT, USA). Results were expressed as milligrams of pelargonidin 3-glucoside equivalents 100 g−1 fresh strawberry. TAC was calculated as follows:where A is the final total absorbance calculated by (A515 – A700) (pH = 1.0) – (A515 – A700) (pH = 4.5); MW is the molecular weight (433.2 g⋅mol−1) of the reference anthocyanin compound (pelargonidin-3-glucoside); DF is the dilution factor; 1000 is the factor for conversion from g to mg; L is the pathlength (1 cm); and ε is molar absorptivity coefficient (36,000 M⋅cm−1).
Statistical analyses.
Data analyses were conducted separately for the two production seasons due to significant seasonal effects shown in a preliminary analysis (Linear Mixed Model). A linear mixed model with the GLIMMIX procedure in SAS (Version 9.4; SAS Institute, Cary, NC, USA) was used for data analysis. Data transformation was not needed after checking normality, homogeneity of variance, and linearity for each dataset. In both seasons, a two-way Linear Mixed Model was performed with strawberry cultivar, biostimulant, and their interactions as the fixed effects and block as the random effect in the model. Multiple comparisons of different measurements among treatments were performed using Tukey’s test at P ≤ 0.05.