Investigating the Impact of Biostimulants on the Row Crops Corn and Soybean Using High-E ﬃ ciency Phenotyping and Next Generation Sequencing

: Row crops represent the most important crops in terms of global cultivated area. Such crops include soybean, corn, wheat, rice, rapeseed, sunﬂower, and cotton. Row crops agriculture is generally an intensive system of farming used to obtain high yields by employing elevated quantities of organic and mineral fertilizers. Considering this, and the decrease in area of arable land, it becomes crucial to ensure high yield and quality using alternative strategies, such as the use of plant biostimulants. These compounds are increasingly recognized as sustainable solution to optimize nutrient uptake, crop yield, quality, and tolerance to abiotic stresses. In this work, by means of high-throughput plant phenotyping, we evaluated the e ﬀ ectiveness of a set of three new foliar biostimulant prototypes (coded as 52096, 52097, 52113) applied on corn and soybean at application rates 2.5 and 5 mL / L (corresponding to 1 and 2 L / ha respectively). This allowed us to select the most e ﬀ ective prototype (52097, commercial name “YieldOn ® ”) in increasing digital biovolume (DB) and greener area (GGA) either in soybean (both application rates) or corn (rate 5 mL / L) and decreasing Stress Index (SI) in soybean (both application rates). Molecular mechanism of action of selected prototype 52097 was subsequently characterized through Next Generation Sequencing (NGS). In corn, genes involved in hormone (cytokinin and auxin) metabolism / catabolism, maltose biosynthesis, sugar transport and phloem loading were upregulated after application of prototype 52097. In soybean, genes involved in nitrogen metabolism, metal ion transport (mainly zinc and iron), sulfate reduction, and amino acid biosynthesis were induced. The proposed approach supports the integration of multiple omics to open new perspectives in the discovery, evaluation, and development of innovative and sustainable solutions to meet the increasing needs of row-crops agriculture. and Formal Analysis; N.V. Methodology, Next Generation Sequencing Data Curation and G.P. Conceptualization, writing, review, editing and supervision.


Introduction
The increase in global population and the uncertainty produced by climate change represent big challenges for current and future agriculture [1,2]. Agricultural activity should ensure crop production plant's ability to intercept light. The use of high-throughput imaging analysis system allows to successfully integrate the experiments involving many variables, a large number of samples, and multiple comparisons [36]. Moreover, the high-throughput image analysis system is a non-invasive method that has the potential to determine the plant phenotypic response to experimental variable(s) (e.g., abiotic stress conditions), throughout the growing (or part of it) of experimental crops [24,25,37].
Hence, very recent papers showed that the use of a "multi-omics" approach, in particular metabolomics and plant phenotyping, represents an effective tool to examine plant performances under different experimental conditions [24,34]. The application of such integrated approach could offer a better explanation of the mechanism(s) of action of different PBS molecules or compounds on crops. This can be obtained by the identification of several biomarkers of PBS action as reported in Paul et al. (2019b) and Ugena et al. (2018) [34,36].
The aim of this study was the selection and characterization of a novel biostimulant formulation conceived to increase the yield of row crops. To achieve such objective, using a phenomic approach we investigated-on corn and soybean-the effect and physiological mechanism(s) of action of three different foliar biostimulant formulations/prototypes. This allowed the selection of the most effective one, subsequently characterized at transcriptomic level to understand its molecular action. This study, based on the integration of phenomic and genomic tools, opens new perspectives to release effective formulations for row-crops agriculture.

Plant Material and Growing Conditions
Experiments were performed on corn (Zea mays L., hybrid P0423, Pioneer) from November to December 2016, and soybean (Glycine max L. Merr.), from May to June 2017.
Both species were sowed directly into white pots (16 cm diameter, 20 cm height), containing 3.5 L of substrate consisting of a 50:50 mixture of peat and river sand. The day before sowing, 20 units of nitrogen, 40 units of phosphate, and 20 units of potassium oxide were added to the substrate mixture. Both soybean and corn plants were irrigated daily with 100 mL water, for 12 days. Afterwards, water was increased to 200 mL until the end of experiment to ensure an adequate water supply.
Three different biostimulant prototypes based on different combinations of seaweed and plant extracts formulated with selected micronutrients such as Mn, Zn, Mo (prototypes coded as 52096, 52097, 52113; proprietary composition Valagro SpA) were sprayed (using a portable atomizer sprayer) on both species at the third true leaf stage. Two rates were applied: 2.5 (lower dosage) and 5 mL/L (higher dosage) (corresponding to 1 and 2 L/ha respectively) during the experiment on soybean, while 5 mL/L was the unique rate applied during the experiment on corn. Untreated control plants were sprayed with distilled water. The experimental setup was composed by 5 biological replicates (plants) for each experimental condition using a completely randomized experimental design.
For NGS and qRT-PCR analysis, samples were collected just before (t = 0) and after the foliar application of PBS 52097 (5 ml/L) at 8, 24, 48, and 168 h. Both corn and soybean 3 leaves were removed from the plants and immediately submerged in liquid nitrogen. For each experimental condition (treatment and time-point), three biological replicates were collected from different plants at the third fully expanded leaf stage. Each biological replicate consisted of three entire leaves (central position) collected from three individual plants and pooled.
Plants used for sampling leaves were excluded from subsequent imaging acquisition or additional leaf sampling. All samples were collected at around 9:00 a.m.

Non-Destructive Measurements
The morphological and physiological characterization of the plants was carried out, non-destructively, by plant imaging. Images of plants were acquired, throughout the experiment, with the plant phenotyping platform Scanalyzer 3D (LemnaTec GmbH). Detailed information on the platform is reported in Briglia et al. (2019) [38]. Briefly, it is composed of 2 imaging chambers visible light (RGB) and fluorescence (FLUO), respectively. FLUO images, recorded into the fluorescence imaging chamber, were used to evaluate the photosynthetic performance through the "Stress Index" (see below). For each imaging chamber three images per plant were taken, one from top view of the plant and two from side view (0 • and 90 • ).

Digital Biovolume Assessment
Plant growth was assessed through the digital biovolume (DB) [39] as follow: where pixel sideview 0 • , 90 • and top view are the plant pixel areas from all sides and top view images.

Color Classification
During the experiment the resulting RGB images were then analyzed by categorizing the pixels according to their color.
After the color segmentation process, that allow to separate the plant from the background, the RGB images were converted to HSI color space (Hue, Saturation, Intensity) and then the hue histogram was calculated. According to Casadesús et al. (2007) [40] the relative greener area (GGA) of each image was calculated as the sum of frequencies of the histogram classes included in the hue angle ranging from 80 • to 180 • . The GGA were used to evaluate the health status of the plant via colour classification (e.g., green: healthy and active leaf surface; yellow: degree of the plant senescence) [40].

Stress Index
The performance of photosynthetic system is not constant and depends mainly on the health and stress condition of a plant. When a plant is placed under stress, more fluorescent light of higher energy is released and this change in the pixel distribution can be measured using the fluorescence imaging chamber.
The Stress Index was calculated according to Petrozza et al. (2014) [25] as (F x − F y )/(F x + F y ) where F is the number of pixels in the x, y color classes, under the assumption that any impairments of the photosynthesis result in a change of pixel number at the x and y color class [25].
The x, y color classes were determined experimentally, by examining the hue histogram. Values of photosynthetic Stress Index vary from +1, poor photosynthetic efficiency, to −1, greater photosynthetic efficiency and should be considered only as a relative level when compared to other plants in the same experiment.

RNA-Seq Analysis
Single samples from leaves from untreated control plants (UTC) and from those treated with PBS 52097 (24 h after application) were used for RNA-sequencing. For each sample, total RNA was isolated using a CTAB-based protocol as described by Chang et al. (1993) [41]. RNA-seq libraries were prepared according to the so-called "dUTP method" to generate mRNA-seq libraries [42,43]. In short, mRNA was purified from 4 µg total RNA using oligo-dT beads, fragmented, and converted to cDNA. Libraries were subsequently made using the Illumina mRNA-Seq Sample Preparation Kit according to the manufacturer's instructions. An amount of 4 pmol of each library was sequenced by BaseClear B.V. (The Netherlands) using the Illumina HiSeq2500 system, with a read length of 50 nucleotides.
Single-end sequencing reads were filtered using the Illumina Casava pipeline version 1.8.3 and Illumina Chastity filtering. Additional filtering on the remaining reads was performed using the FASTQC quality control tool version 0.10.0. For RNA-seq analysis, sequence reads were mapped (per sample) to the reference using CLC Bio Genomics Workbench software (version 5.1.5). As a reference for corn, the publicly available B73 reference sequence (AGPv3.22, downloaded from the ZmGDB genome browser) consisting of 63,241 sequences was used. For soybean, the Glycine max_275_Wm82.a2.v1 primary transcripts [44], consisting of 88,647 sequences were used as a reference. To determine gene expression levels and differential gene expression, RPKM values (read counts corrected for library size and transcript length) were calculated using the CLC Bio software. Differentially expressed genes (DEGs) were selected by calculating the ratio of the RPKM value of treated samples over the RPKM value of untreated samples. Only genes with at least 50 reads and an RPKM value of over 5 in at least one sample were considered.
To functionally categorize the corn and soybean gene sequences, gene ontology (GO) terms were assigned to each assembled contig using Blast2GO software (version 3.1). GO terms provide a controlled vocabulary to describe the functions of genes across species. Blast2GO is an automated tool for the assignment of GO terms based on sequence similarity [45]. Statistical assessment of GO term enrichment in groups of DEGs were done using Fisher s Exact Test in combination with false discovery rate (FDR) correction for multiple testing.

qRT-PCR Analysis
qRT-PCR analysis was performed on single samples collected 8, 24, 48, and 168 h after PBS(s) application. Total RNA was isolated as described above. First-strand cDNA was prepared using 80 ng total RNA and qScript™ cDNA Supermix (Quanta Biosciences). Two and a half µL of 1:5 diluted first-strand cDNA was used as a template in the subsequent PCR, performed on a Bio-Rad CFX using 5 pmol of both primers (sequence of primers reported in Supplemental Table S2) and PerfeCta SYBR Green SuperMix UNG (Quanta Biosciences) in a final volume of 12.5 µL per reaction, according to the manufacturer s instructions. All transcript levels were normalized using a eukaryotic translation initiation factor gene (corn) or actin gene (soybean) as a control.

Statistical Analysis of Data
The statistical analysis was performed using R software (3.3.2 version; R foundation for Statistical Computing, Vienna, Austria). Phenotyping results were analyzed using one-way analysis of variance (ANOVA), and the means were compared with Duncan's New Multiple Range Test (p < 0.05).

Digital Biovolume
Plant development was assessed through the DB which is a morphometric measurement previously employed in high-throughput (HTP) studies to monitor the influence of abiotic stresses, mainly drought, on plant growth [34,38].
The DB of soybean plants was significantly improved after 10 days from the application of each of the three prototypes tested both at lower (2.5 mL/L) and higher (5 ml/L) concentration, in comparison with UTC ( Figure 1A,B). However, the most consistent results (higher DB increase) were obtained using formulation 52097 at the lower dosage (2.5 mL/L; Figure 1A) reaching-at the end of experiment-the mean DB value of 49.58 k units (+72% compared to UTC plants). Such improvement in DB was clearly observed after 10 days from prototype application and maintained during time until the end of DB measurements in soybean, which was 26 days after treatment ( Figure 1A).
Parallel measurements on corn confirmed that the 52097 PBS prototype exerted the higher increasing effect on DB in comparison with the other treatments and UTC plants ( Figure 1C). In this case, beside an early plant response to the treatment observed already 7 days after treatment, the plants treated with 52097 maintained a greater DB throughout the experiment. No statistically significant differences between the 52113 and control plants were noted. At day 10 after treatment the plants treated with the 52097 reached a DB level of 27.92 k units, 81% higher than that of the UTC plants.
For both soybean and corn, plants treated with prototype 52097 showed a constant and consistent increase in DB accumulation compared to plants treated with 52096 or 52113.

Greener Area
The value of greener area (GGA) in soybean showed a general decrease during time ( Figure  2A,B), as expected, due to the progression of phenological phases, which lead to senescence and yellowing. This was observed in untreated soybean plants, but also in treated with prototypes 52096 and 52113. Interestingly, soybean plants treated with prototype 52097 at both rates showed a consistent persistence of optimal (ranging from 0.7 to 0.8) GGA values during time (Figure 2A,B), even during the latest phenomic measurements (19)(20)(21)(22)(23)(24)(25)(26).
For the first 12 days after treatment (DAT), no statistically significant differences were observed between the treatments, since all plants showed the same mean GGA value around 0.8 (Figure 2A,B). Starting from 14 days after treatment the first yellowing/sign of senescence were recorded on the untreated control plants and the plants treated with 52113 and 52096 prototypes. At the end of the experiment it was possible to see how, at both rates tested, the application of 52097 on soybean allowed a higher level of GGA, in particular around 0.77 (2.5 ml\L dosage) and 0.75 (5 mL\L dosage). On the other hand, control plants and plants treated with 52113 and 52096 prototypes reached values between 0.40 and 0.28 respectively (Figure 2A, B). It can be concluded that only prototype 52097 was able to preserve and improve GGA in comparison with the other experimental conditions. This can be attributed to a positive effect of prototype 52097 on the "stay-green" condition, that is known to allow plants to maintain high photosynthetic activity [46], and gain benefits on biomass accumulation, as confirmed by the data previously shown on DB (Figure 1).
The positive results on GGA observed after application of prototype 52097 were visibly clear by looking at the set of pictures taken during the cycle, by mean of the visible camera ( Figure 2C).
Considering the same test on corn, all the formulations exerted a slight increase in GGA in comparison with untreated control (UTC) plants ( Figure 2D), statistically significant at 7, 10, and 13 days after treatment with prototypes 52096 and 52097.

Greener Area
The value of greener area (GGA) in soybean showed a general decrease during time (Figure 2A,B), as expected, due to the progression of phenological phases, which lead to senescence and yellowing. This was observed in untreated soybean plants, but also in treated with prototypes 52096 and 52113. Interestingly, soybean plants treated with prototype 52097 at both rates showed a consistent persistence of optimal (ranging from 0.7 to 0.8) GGA values during time (Figure 2A,B), even during the latest phenomic measurements (19)(20)(21)(22)(23)(24)(25)(26).
For the first 12 days after treatment (DAT), no statistically significant differences were observed between the treatments, since all plants showed the same mean GGA value around 0.8 (Figure 2A,B). Starting from 14 days after treatment the first yellowing/sign of senescence were recorded on the untreated control plants and the plants treated with 52113 and 52096 prototypes. At the end of the experiment it was possible to see how, at both rates tested, the application of 52097 on soybean allowed a higher level of GGA, in particular around 0.77 (2.5 ml/L dosage) and 0.75 (5 mL/L dosage). On the other hand, control plants and plants treated with 52113 and 52096 prototypes reached values between 0.40 and 0.28 respectively (Figure 2A,B). It can be concluded that only prototype 52097 was able to preserve and improve GGA in comparison with the other experimental conditions. This can be attributed to a positive effect of prototype 52097 on the "stay-green" condition, that is known to allow plants to maintain high photosynthetic activity [46], and gain benefits on biomass accumulation, as confirmed by the data previously shown on DB (Figure 1).
The positive results on GGA observed after application of prototype 52097 were visibly clear by looking at the set of pictures taken during the cycle, by mean of the visible camera ( Figure 2C). Considering the same test on corn, all the formulations exerted a slight increase in GGA in comparison with untreated control (UTC) plants ( Figure 2D), statistically significant at 7, 10, and 13 days after treatment with prototypes 52096 and 52097.

Stress Index
The measurement of Stress Index of treated and untreated soybean plants did not show statistically significant differences during the first four data acquisitions time ( Figure 3). As expected, due to the plant cycle progression and senescence, during the last three data acquisitions a higher Stress Index-ranging from 0.4 to 0.7-was observed for UTC. Treating plants with 52096 or 52113 did not affect the Stress Index. Interestingly, soybean plants treated with 52097 showed a lower level of Stress Index than the other treatments, with values stable around 0.2 ( Figure 3).   Table S3).

Molecular Analyses
Based on the obtained results, although some differences were observed between soybean and corn, we selected compound 52097 as the best candidate for further NGS analyses. Leaf tissue from untreated soybean and corn plants was compared to its 52097-treated counterpart by RNA-seq analysis. Per sample, over 25 million single reads were generated and mapped to the relevant reference transcriptome (see Material and Methods). For corn, around 77% of the available sequencing reads could be mapped to this assembly, for soybean this was 89%. By comparing 52097treated samples to untreated controls, differentially expressed genes (DEGs) were identified. Naturally, the number of DEGs depended on the fold-change threshold applied (Table 1). In general, the number of DEGs was higher in corn than in soybean. Lists of the 20 most upregulated genes for both crops are provided in Supplemental Tables S4 (soybean) and Table S5 (corn).

Functional Annotation Using Gene Ontology
Of the 63,241 sequences present in the reference transcriptome from corn, 49,426 sequences (78%)  Table S3).

Molecular Analyses
Based on the obtained results, although some differences were observed between soybean and corn, we selected compound 52097 as the best candidate for further NGS analyses. Leaf tissue from untreated soybean and corn plants was compared to its 52097-treated counterpart by RNA-seq analysis. Per sample, over 25 million single reads were generated and mapped to the relevant reference transcriptome (see Material and Methods). For corn, around 77% of the available sequencing reads could be mapped to this assembly, for soybean this was 89%. By comparing 52097-treated samples to untreated controls, differentially expressed genes (DEGs) were identified. Naturally, the number of DEGs depended on the fold-change threshold applied (Table 1). In general, the number of DEGs was higher in corn than in soybean. Lists of the 20 most upregulated genes for both crops are provided in Supplemental

Functional Annotation Using Gene Ontology
Of the 63,241 sequences present in the reference transcriptome from corn, 49,426 sequences (78%) could be functionally annotated using GO, meaning that one or more biological processes, molecular functions, or cellular localizations could be linked to these sequences based on sequence homology. For soybean, from the 88,647 genes present in the soybean transcriptome, 70,776 sequences (80%) could be functionally annotated using GO.
For corn, all DEGs up-or downregulated more than 3-fold were used for enrichment analysis. Several biological processes, including nitrogen assimilation, maltose biosynthesis, and cytokinin metabolism were enriched among the 331 upregulated genes from corn ( Figure 4A). Analysis of the 593 downregulated corn transcripts resulted in 55 enriched GO-terms (biological process). By filtering out the most reduced GO-terms, that is, removing parent terms of already present statistically significant child GO terms, a list of 13 significantly enriched biological processes remained ( Figure 4B). These terms included divalent metal ion transport, response to carbohydrate, phenylalanine degradation, and flavonoid biosynthesis.
For soybean, first, all DEGs up-or downregulated more than 3-fold were used for enrichment analysis. Analysis of the 65 upregulated soybean transcripts resulted in 25 enriched GO-terms (biological process), including metal ion transport, sulfate reduction, asparagine biosynthesis, and serine metabolism (not shown). However, there were no significantly enriched GO terms among the 59 downregulated soybean transcripts. For this reason, it was decided to include all soybean contigs with a FC greater than 2. This resulted in 16 significantly enriched (reduced) GO terms for the 278 upregulated contigs (Figure 4C), and 13 for the 321 downregulated contigs, including auxin-activated signaling, sulfur amino acid metabolism, iron transport, and sulfur compound biosynthesis ( Figure 4D).
From both crops, two individual DEGs were selected using the results from the GO enrichment analysis. For corn, these were a cytokinin dehydrogenase (CKX; Figure 5A) and a glutamine synthetase (GS; Figure 5B). CKX catabolizes the plant hormone cytokinin and plays an important role in cytokinin regulated processes [47]. GS is required for nitrogen assimilation and allocation within the plant, and for nitrogen remobilization in both source and sink tissues [48]. GS is important for ammonium assimilation in roots, during senescence, and during photorespiration. Several studies have indicated that GS plays an essential role in plant development and yield formation in cereals. For example, in corn, leaf-localized GS are of specific importance for the development of the cob with respect to kernel number and kernel size [49]. A putative GS gene is induced by treatment with product 52097 ( Figure 5B). In addition, 3 more genes from corn were selected on the basis of their functionality: An SPX domain-containing protein (named after SYG1/Pho81/XPR1 proteins; Figure 5C), a NRT1/PTR FAMILY (NPF) protein ( Figure 5D, and a polyol/monosaccharide transporter (PMT; Figure 5E). It is well described that plant growth and development are highly dependent on the availability of inorganic phosphate (Pi). Among the many proteins involved in the plant response to Pi starvation, proteins containing the SPX domain are key players involved in the maintenance of internal levels of Pi. Indeed, SPX genes have been reported to be induced upon Pi starvation in roots and shoots, and proteins harboring the SPX domain have been shown to be involved in P use efficiency [50]. Members of the plant NITRATE TRANSPORTER 1/PEPTIDE TRANSPORTER (NRT1/PTR) family display protein sequence homology with peptide transporters in animals. In comparison to their animal and bacterial counterparts, the plant NRT1/PTR family proteins transport a wide variety of substrates: nitrate, peptides, amino acids, dicarboxylates, glucosinolates, IAA, and ABA [51]. The transcript identified here shows the highest similarity to the first identified member of the NRT1/PTR family: NRT1.1. NRT1.1 is an Arabidopsis nitrate transporter that also functions as a nitrate sensor and can transport auxins. As such, it links nutrient and hormone signaling [51]. PMTs are proteins capable of transporting a range of sugar alcohols and monosaccharides including glucose, fructose, sorbitol, mannitol, xylitol, xylose, and galactose [52]. PMTs are believed to be involved in phloem loading [52]. Hence, the induction of a PMT gene identified in this study could point towards increased phloem loading.
For all these genes, expression was determined in leaf samples collected on several timepoints after application of prototype 52097. It was observed that data obtained by qRT-PCR corroborated the NGS results. Twenty-four hours after application, the differences in gene expression between treated and untreated leaves were very comparable. The additional timepoints showed that GS and SPX expression already peaked 8 h after application ( Figure 5B,C), whereas the other 5 genes reached their maximum expression after 24 h (Figure 5A,D-G).  represented by the bars in the graphs on the right (numbers of contigs in the test/reference sets are reproduced next to each bar).

Conclusions
This study highlights the use of high-throughput/efficiency plant phenotyping (phenomics) together with Next Generation Sequencing to investigate the effectiveness and mechanism of action of new biostimulant formulations. Such formulations were specifically conceived as foliar applications to increase yield of different row/industrial crops, such as corn and soybean.
Phenomic-based measurements of digital biovolume, Greener Area, and Stress index allowed us to select 52,097 (commercial name "YieldOn ® ") as the most effective prototype among the ones tested. Subsequently, through NGS a deep characterization of the molecular mechanisms by which the biostimulant under investigation exerts its positive effect was performed. This analysis explained the mechanism of action of the biostimulant under investigation, which in corn upregulated specific processes like nitrogen and phosphate assimilation and metabolism, maltose biosynthesis, sugar For soybean, we selected a ZIP (ZRT, IRT-like protein; Figure 5F) transporter and an asparagine synthetase (AS; Figure 5G). ZIP transporters are important during uptake and transport of zinc and iron and other divalent metal cations [53]. AS, like GS, functions in nitrogen metabolism [54].
For all these genes, expression was determined in leaf samples collected on several timepoints after application of prototype 52097. It was observed that data obtained by qRT-PCR corroborated the NGS results. Twenty-four hours after application, the differences in gene expression between treated and untreated leaves were very comparable. The additional timepoints showed that GS and SPX expression already peaked 8 h after application ( Figure 5B,C), whereas the other 5 genes reached their maximum expression after 24 h (Figure 5A,D-G).

Conclusions
This study highlights the use of high-throughput/efficiency plant phenotyping (phenomics) together with Next Generation Sequencing to investigate the effectiveness and mechanism of action of new biostimulant formulations. Such formulations were specifically conceived as foliar applications to increase yield of different row/industrial crops, such as corn and soybean.
Phenomic-based measurements of digital biovolume, Greener Area, and Stress index allowed us to select 52,097 (commercial name "YieldOn ® ") as the most effective prototype among the ones tested. Subsequently, through NGS a deep characterization of the molecular mechanisms by which the biostimulant under investigation exerts its positive effect was performed. This analysis explained the mechanism of action of the biostimulant under investigation, which in corn upregulated specific processes like nitrogen and phosphate assimilation and metabolism, maltose biosynthesis, sugar transport and phloem loading, hormone (cytokinin) metabolism. In soybean nitrogen metabolism, metal ion transport (mainly zinc and iron), sulfate reduction, and amino acid biosynthesis were upregulated.
In conclusion, the results showed in this work support the integration of multiple "omics" as robust and objective tools in the discovery, evaluation, and development of innovative, sustainable, and targeted solutions to meet the emerging needs of row-crops agriculture.
Supplementary Materials: The following are available online at http://www.mdpi.com/2073-4395/9/11/761/s1, Table S1: Temperatures, humidity (RH%), and radiation (PAR µmol m −2 s −1 ) parameters measured during the trial. Table S2: Sequence of primers used for qRT-PCR analysis; Table S3: Stress Index measurements in corn; Table  S4: lists of the 20 most upregulated genes for soybean; Table S5: lists of the 20 most upregulated genes for corn. Funding: The research conducted was funded by Valagro SpA, which provided biostimulant prototypes that were blindly identified by a numbered code to avoid any conflict of interest in the data production. The prototypes were named only at the end of the experiments.