Application of Biostimulants Containing Amino Acids to Tomatoes Could Favor Sustainable Cultivation: Implications for Tyrosine, Lysine, and Methionine

: Tomato ( Solanum lycopersicum L.) is one of the most important crops worldwide as per its production and the surface cultivated. The use of biostimulant products plays a fundamental role in mitigating the negative e ﬀ ects of climate change and reducing the use of conventional fertilizers. Many of these products are formulated with amino acids (AAs). This study was conducted to elucidate the e ﬀ ects of the foliar application of tyrosine (Tyr) (15 mM), lysine (Lys) (15 mM), methionine (Met) (15 mM), and a Tyr + Lys + Met (15 mM + 15 mM + 15 mM) mixture on the physiological and metabolic processes, vegetative growth, and nutritional state of Optima variety tomato plants. The results showed that application of the AAs, individually and combined, was beneﬁcial for the growth of the aerial part, net assimilation of CO 2 , and water use e ﬃ ciency (WUE). Application of Tyr resulted in the best WUE. The metabolomics study revealed that AA treatments increased the concentration of proline, fructose, and glucose, whose role was to stimulate glycolysis and the Krebs cycle. Thus, the plants could have greater reduction power and energy, as well as more carbon molecules for their growth processes.


Introduction
The tomato (Solanum lycopersicum L.) is one of the most important crops worldwide on the economic and social level [1]. Its fruit is of great interest to today's society, as it is a source of compounds with great human health benefits [2]. According to the statistical data of the Food and Agriculture Organization (FAO) (2020) [3], Spain is the second-largest producer of tomatoes in Europe, with approximately 56 thousand hectares cultivated and a yield of 85 t/ha.
In the last few years, agricultural production has been threatened by climate change, which brings with it high temperatures and periods of drought that are long-lasting and more intense. Along with this, the world population has increased considerably. This is said to be the general trend in the coming and ethylene, aside from favoring resistance to abiotic stresses through the S-adenosyl methionine (SAM) pathway [23,24]. Rai [25] showed that exogenous application of this amino acid could help plants tolerate or resist various abiotic stresses. Nevertheless, limited scientific information has been found on the specific effect that individual application of these compounds has on the physiological and metabolic processes of plants, and even less information is available on the simultaneous application of these AAs as a mixture. Knowing all of these aspects could be very useful when formulating the newest generation of biostimulants, either through the addition of free amino acids or through the selection of raw materials that contain AAs in the ratio and concentration desired. Therefore, the main objective of the present study was to understand the effects that exogenous application of the amino acids tyrosine, lysine, and methionine, either individually or as a mixture, has on the physiological, metabolic, and nutritional processes of tomato plants of the 'Optima' variety, and to determine if these changes have an effect on the growth of the plants.

Growing Conditions and Plant Material
The study took place in a multi-tunnel greenhouse at the experimental field of the CEBAS "La Matanza", located in the municipality of Santomera (Murcia). To monitor the temperature, a cooling system and aluminum shade cloths (30%) were utilized, which maintain the optimum environmental conditions for the crops. The environmental characteristics of the greenhouse during the experiment were as follows: photosynthetically active radiation (PAR) of 1000 nmol m −2 s −1 , day/night temperatures of 32/19 ± 3 • C, relative humidity (RH) day/night of 52/80 ± 5%, and a natural photoperiod of 15 h.
For this study, tomato plants (Solanum lycopersicum L.) var. 'Optima' were utilized, obtained from the nursery Babyplant S.L. (Santomera, Murcia, Spain). The seedlings were transplanted to 3.5 L pots containing a "Universal" substrate (Projar professional, 2018, Quart de Poblet, Spain), after which they were placed in the greenhouse described above. The substrate was composed of 50% yellow and black peat, 40% coconut fiber, and 10% perlite. Once the plants were transplanted, water was applied at pH 8.3 and an electric conductivity (EC) of 0.9 dS m −1 in an irrigation system composed of lines containing 4 L h −1 self-compensating drippers. The volume of water applied was sufficient for drainage to be observed after each irrigation event. Two weeks after irrigating with only water, the plants were irrigated with 50% Hoagland nutrient solution for one minute a day, composed of KNO 3 (27 g 100 L −1 ), Ca(NO 3 ) 2 (42 g 100 L −1 ), KH 2 PO 4 (7 g 100 L −1 ), MgSO 4 (13 g 100 L −1 ), Fe-EDTA (1 g 100 L −1 ), and micronutrients (1 g 100 L −1 ; Hidromix S, Valagro). After two weeks, the fertilization was increased to 100% Hoagland solution, twice a day, five minutes in the morning and five minutes in the afternoon. Additionally, there was 15% drainage after each irrigation. One week after irrigating with the 100% solution, the different AA treatments were applied.

Amino Acid Preparation
In this study, the amino acids tyrosine (Tyr), lysine (Lys), and methionine (Met) were assayed and were acquired from Calid Iberica (Barcelona, Spain). For the preparation of each AA (each was prepared independently), a stock solution for each was prepared by weighing 7 g and diluting with 100 mL of milliQ water. Once the stock solutions were prepared, the treatment solutions were prepared for foliar application, including (i) control treatment (milliQ water), (ii) treatment with Tyr (15 mM), (iii) treatment with Lys (15 mM), (iv) treatment with Met (15 mM), and (v) treatment with Tyr + Lys + Met (15 mM + 15 mM + 15 mM; a mix of the three AAs). The pH of each solution was adjusted to 5.5-7.5, and Tween-20 was added at a final concentration of 0.1%. Tween-20 is a wetting agent that improves and facilitates the adherence of solutions to the leaves of plants.

Application of the Treatments
Diverse treatments were applied so that most of the foliar area was covered by the formulations. This was done one week after starting the irrigation with 100% Hoagland nutrient solution. The day after application, the plants were sprayed with distilled water to ensure greater absorption of the product, as the leaf cuticle swells and its permeability increases. Treatments were applied to 12 plants (four experimental units with three plants per experimental unit).

Gas Exchange Parameters
Gas exchange parameters were measured in fully expanded and developed leaves between 8:30 and 11:00 a.m. in two plants per experimental unit one week after foliar application. For this, a portable photosynthesis analyzer was utilized (PP System Ciras2, Hitchin, UK), whose measurement conditions were adjusted to a PAR = 1200 mmol m −2 s −1 and 400 ppm of CO 2 in the measurement chamber. This equipment measured the net assimilation rate of CO 2 (A CO2 ) and stomatal conductance (g s ), and it calculated the water use efficiency (WUE = A CO2 /E leaf , where E leaf corresponds to the value of leaf transpiration obtained in each measurement) [26].

Chlorophyll Fluorescence and Concentration Parameters
An FMS-2 portable modulated pulse fluorimeter (Hansatech Instruments Ltd., Pentney, UK) was utilized. Measurements were taken with the same plants used to determine the gas exchange parameters. The measurement conditions are described in Simón-Grao et al. [26]. The CL-01 portable measurement device (SPAD units, Hansatech) was utilized to measure the chlorophyll concentration, expressed as SPAD units.

Growth Parameters
Height was recorded after measuring the chlorophyll fluorescence and gas exchange parameters before the final harvest, and the diameter of the plants was measured at three different heights. In the final sampling, the shoots (leaves and stem) were weighed with a Sartorius precision scale (Acculab). Afterwards, four leaflets were taken, weighted, and washed with deionized water. They were then dried in an oven at 60 • C for at least 48 h. This material was used to determine the concentration of mineral nutrients. Leaflets were also collected, frozen with liquid nitrogen, and stored at −80 • C for posterior metabolic analysis. Leaflet samples for the ionomic and metabolomic study were chosen randomly at a middle height of the plant, ensuring that these were fully developed.

Measuring the Concentration of Mineral Nutrients in Leaf Tissues
The concentrations of Na, K, Mg, Ca, P, S, Fe, Cu, Mn, Zn, and B were determined from oven-dried leaf samples. These were ground to a fine powder and digested with HNO 3 :H 2 O 2 (5:3 v/v) using a microwave (CERM Mars Xpress, Matthews, NC, USA) with a temperature ramp that reached 200 • C. These were posteriorly analyzed using inductively coupled plasma mass spectrometry (ICP-MS, Iris Intrepid II, Thermo Electron Corporation, Franklin, TN, USA) [23]. Total C and N were analyzed with a Thermo Finnigan C/N elemental analyzer (Milan, Italy).

Metabolomic Analysis of Leaf Tissues
A "non-directed" metabolomic analysis was performed with leaflet samples frozen in liquid nitrogen. These samples were lyophilized and ground to obtain a fine powder. The extraction protocol utilized to prepare the samples for nuclear magnetic resonance (NMR) was the same as that utilized by Van der Sar et al. [27], and the equipment utilized for the analysis was a 500 MHz Bruker spectrometer (Bruker Biospin, Rheinstetten, Germany) with a 5 mm broadband Prodigy BBO N 2 cryoprobe.
The NMR equipment provided a frequency graph of the detected signals. The resulting spectra were sent to the Chenomx NMR suite version 8.3 processor (Chenomx, Edmonton, Alberta, Canada) to quantify the metabolites, where the samples were standardized with Trimethylsilylpropionic acid sodium salt (TSP) internal standard, and the pH was maintained at around 6. This processor contains a large database that is used to detect metabolites found at concentrations greater than parts-per-million. Among the metabolites found and/or quantified, the following are underlined: aspartate, glutamate, alanine, glutamine, isoleucine, valine, tyrosine, proline, phenylalanine, citrate, formate, fumarate, malate, fructose, glucose, sucrose, 4-Aminobutyrate (GABA), and trigonelline [28,29].

Statistical Analysis
This study utilized a single-factor experimental design, in which the application of different AA treatments was studied. Statistical analysis consisted of one-way analysis of variance (ANOVA) (treatments) performed with the statistical package SPSS version 24. The values presented for each treatment were from a total of four biological units (n = 4). When ANOVA was significant (p < 0.05), Duncan's multiple range test was utilized to separate the means. Pearson's correlation coefficient was utilized with the same statistical package, which determined the linear relationship between the diverse variables. Student's t test was utilized to verify the existence of significant differences detected between the metabolites and the control. For this, the data were transformed through the square root of the value divided by 100. Likewise, a principal component analysis (PCA) and a cluster analysis (CA) were performed.

Vegetative Growth Parameters
The vegetative growth parameters that were significantly affected by the treatments were the stem diameter and the dry weight of the shoot (Sh; Table 1). The control plants obtained values of 13 g plant −1 and 11 mm plant −1 for the shoot biomass and diameter of the stem, respectively. Exogenous foliar application of Tyr, Lys, Met, and Tyr + Lys + Met showed shoot (Sh) dry weight (dw) values greater than those of the control plants (between 19 and 21 g dry weight (dw)). Control plants had a greater stem diameter, although significant differences were only found for the Tyr and Tyr + Lys + Met treatments. As for the height, significant differences between the different treatments were not found (Table 1).

Physiological Study
For the gas exchange parameter, significant differences were observed between the treatments. The control plants obtained values of A CO2 = 16.8 µmol m −2 s −1 , g s = 1196.3 mmol m −2 s −1 , and WUE = 2.38 µmol CO 2 mmol −1 H 2 O (Table 2). These values increased for A CO2 and WUE and decreased for g s with the AA treatments. Individual AA treatments increased A CO2 up to 23% as compared to the control, while the Tyr + Lys + Met treatment increased A CO2 up to 16% without significant difference from the control. As for the g s parameter, the control plants obtained the highest values as compared to plants treated with the AAs. Among the latter, the plants treated with Tyr obtained the lowest values, followed by Met, Lys, and Tyr + Lys + Met. As for WUE, the lowest values were found in the control treatment, followed by Lys, Met, and Tyr + Lys + Met, and Tyr obtained the lowest values among AA-treated plants.  In the ANOVA, 'ns' indicates non-significant differences with a confidence interval of 95%; * and *** indicate significant differences at p < 0.05 and 0.001, respectively. Different lowercase letters indicate significant differences (p < 0.05) between the means, as established by Duncan's test (n = 4).

Gas Exchange Parameters
As for chlorophyll fluorescence, ANOVA did not reveal significant differences between the different treatments ( Table 2). The chlorophyll content measurements were significantly different in the fully developed leaves (DL) and the buds (LB). In this case, the control plants had chlorophyll values of 24.8 and 20.1 SPAD units for DL and LB, respectively ( Table 2). The plants treated with Met had the highest chlorophyll (Chl) DL values, with 30.8 SPAD units (although only significant with the Tyr + Lys + Met treatment). On the other hand, plants treated with Lys had the highest LB values as compared with the other treatments (although this was only significant with the control and Tyr + Lys + Met treatments).

Ionomic Study
The nutritional states of the plants one week after foliar application of the AAs were significantly different between the Ca, K, Mg, Na, P, B, Cu, Fe, and Mn treatments ( Table 3). The control plants had the highest foliar concentrations, reaching values of 4.22, 4.25, 1.03, 0.34, and 0.46 g 100 g −1 dw for Ca, K, Mg, Na, and P, respectively. AA treatments, independently of their nature, produced significant decreases in the foliar concentration, although a significant difference between AA treatments was found only for Na. Thus, the plants treated with Met and Lys obtained values that were lower than those treated with Tyr + Lys + Met. Table 3. Concentration of macro (g 100 g −1 dw) and micronutrients (ppm) quantified in leaves from the tomato plant of var. Optima one week after exogenous application of AA treatments: Control (without AAs), L-Tyrosine (L-Tyr), L-Lysine (L-Lys), L-Methionine (L-Met), and the Tyr + Lys + Met combination.
In the ANOVA, 'ns' indicates non-significant differences with a confidence interval of 95%; ** and *** indicate significant differences at p < 0.01 and 0.001, respectively. The different lowercase letters indicate significant differences (p < 0.05) between the means, as established by Duncan's test (n = 4).
As for the micronutrients analyzed, the control plants obtained values of 49.1, 9.1, 172.6, 78.3, and 28.0 ppm for B, Cu, Mn, and Zn, respectively. As for B and Cu, the same was found for the micronutrients, in that the greatest concentration was found for plants not treated with AAs. In the case of Fe concentration, the greatest was found in the plants treated with Lys (210.8 g 100 g −1 dw), while for the remaining treatments, the concentration ranged between 145 and 158 ppm (Table 3). Zn was the only microelement that was not significantly different between treatments.

Amino Acids
Nuclear Magnetic Resonance (NMR) was utilized to detect the following AAs: glutamate, proline, glutamine, aspartate, alanine, phenylalanine, valine, tyrosine, and isoleucine. These amino acids can be divided into two groups according to their concentration. The first is composed of glutamate, proline, glutamine, and aspartate, whose concentrations were higher than 2 mg g −1 dw, with glutamate having the highest concentration of 7.7 mg g −1 dw. These amino acids represented 81% of the total; of this, 37% was glutamate, 25% proline, and 19% glutamine (Table 4, Figure S1). The second group was composed of alanine, phenylalanine, valine, tyrosine, and isoleucine, with concentrations lower than 1 mg g −1 dw. Table 4. Concentration of the amino acids (mg g −1 dw) quantified with NMR in leaves of the tomato plant of var. Optima one week after exogenous application of the AA treatments: Control (without AAs), L-Tyrosine (L-Tyr), L-Lysine (L-Lys), L-Methionine (L-Met), and the Tyr + Lys + Met combination.
Foliar application of AAs changed the plant profile in the leaves with respect to those not treated with AAs. Applications of Tyr, Lys, Met, and their combination increased the concentration of proline from 5 to 9.18 mg g −1 dw, with no differences observed between AA treatments. The concentration of glutamine decreased with the application of the four AA treatments with respect to the control treatment. This decrease was similar to the individual treatments with Tyr, Lys, and Met, and was lower than for the application of the Tyr + Lys + Met combination. The concentration of aspartate only increased with applications of Tyr and the Tyr + Lys + Met mix, and the concentration of isoleucine increased with the application of all four AA treatments.

Organic Acids
The organic acids identified and quantified with H-NMR in the tomato plant leaves were malate, citrate, fumarate, and formate (Figure 1), whose concentrations decreased from malate (39.2 mg g −1 dw) to formate (0.07 mg g −1 dw), so that the concentrations of malate and citrate represented 79% and 21%, respectively, of the total organic acids quantified ( Figure S2). AA treatments changed the concentrations of organic acids in the leaves with respect to the control treatment. All AA treatments increased the concentration of malate significantly as compared to the control, with the Tyr + Lys + Met treatment reaching the highest value of 57.2 mg g −1 dw. No significant differences were found between the control plants and those treated with lysine. As for citrate and fumarate, all AA-treated plants, independently from single AAs or their mixture, showed lower values as compared to the control. Significant differences were also observed between the AA treatments.
Thus, the plants treated with the Tyr + Lys + Met mixture had the lowest concentration of citrate as compared to the rest of the treatments, and the plants treated with Lys had the greatest concentration of fumarate (Figure 1). An increase in formate concentration was observed in the following order: Tyr < Lys = Met < Tyr + Lys + Met. However, despite the changes in plants after application of the diverse treatments, all AA concentrations relatively followed the same distribution pattern with respect to the control plants ( Figure S2).

Sugars
NMR analysis identified and quantified the sugars fructose, glucose, and sucrose, with the highest being fructose (44%), followed by sucrose (31%) and glucose (25%) (Figures 2 and S3). The application of AA treatments affected the concentrations of these sugars. NMR analysis identified and quantified the sugars fructose, glucose, and sucrose, with the highest being fructose (44%), followed by sucrose (31%) and glucose (25%) (Figure 2 and S3). The application of AA treatments affected the concentrations of these sugars.
For fructose, the AA treatments increased the concentration significantly with respect to the control. Among the AA treatments, the greatest value was found for Met and the lowest for Lys. For sucrose, two different responses were observed: an increase with the application of Tyr, reaching a value of 16.2 mg g −1 dw, and a decrease with the application of Lys, Met, and the AA mixture, with the lowest concentration obtained for Lys (4.1 mg g −1 dw). The concentration of glucose increased significantly in all AA treatments with respect to the control treatment. Between the AA treatments, the lowest value was found for Lys.

Other Metabolites
Other metabolites were also quantified, such as 4-Aminobutyrate (GABA) and trigonelline ( Figure 3); in control conditions (without AAs), their concentrations were 0.56 and 0.79 mg g −1 dw, For fructose, the AA treatments increased the concentration significantly with respect to the control. Among the AA treatments, the greatest value was found for Met and the lowest for Lys. For sucrose, two different responses were observed: an increase with the application of Tyr, reaching a value of 16.2 mg g −1 dw, and a decrease with the application of Lys, Met, and the AA mixture, with the lowest concentration obtained for Lys (4.1 mg g −1 dw). The concentration of glucose increased significantly in all AA treatments with respect to the control treatment. Between the AA treatments, the lowest value was found for Lys.

Other Metabolites
Other metabolites were also quantified, such as 4-Aminobutyrate (GABA) and trigonelline ( Figure 3); in control conditions (without AAs), their concentrations were 0.56 and 0.79 mg g −1 dw, respectively. Both compounds showed significant differences in treatment application. The concentration of GABA increased with the application of AAs as compared to the control, with the Tyr + Lys + Met treatment greatly increasing its concentration (3.8 mg g −1 dw). Trigolline, on the other hand, suffered a decrease when the AAs were applied to the plants. This decrease was similar for the four AA treatments, with a range of values between 0.55 and 0.62 mg g −1 dw. respectively. Both compounds showed significant differences in treatment application. The concentration of GABA increased with the application of AAs as compared to the control, with the Tyr + Lys + Met treatment greatly increasing its concentration (3.8 mg g −1 dw). Trigolline, on the other hand, suffered a decrease when the AAs were applied to the plants. This decrease was similar for the four AA treatments, with a range of values between 0.55 and 0.62 mg g −1 dw.

Discussion
The application of amino acids is an agricultural practice that is becoming more widespread, as it supports sustainable agriculture thanks to the role that amino acids play in the physiological, metabolic, and biochemical processes of plants (e.g., water relations, photosynthesis, antioxidant capacity, nutrient absorption, and stress signaling, to name a few), which have an effect on the plant's vegetative and reproductive development as well as the quality of the crops [12]. Thus, amino acids are some of the most utilized components when formulating biostimulants. These can be added as free amino acids or through sources of botanical or animal origin that contain a great number and concentration of amino acids [11]. However, to achieve greater efficiency in the use of biostimulant products, it is important to first understand the effects that individual AAs have on the physiological, biochemical, and metabolic processes for each crop of interest, and also to determine the possible antagonistic, synergistic, or neutral effects when these are applied as mixtures.
In our assays, exogenous application of amino acids, either individually (Tyr, Lys, and Met) or as a mixture (Tyr + Lys + Met), significantly increased the growth rate of tomato plants as compared to the control plants. Therefore, foliar application of these AAs has beneficial effects on tomato plants (Table 1). In other studies, it was also observed that application of these AAs had positive effects on the crops. In a study conducted by El-Sherbeny and da Silva [30] in beet plants, foliar application of tyrosine at a concentration of 0.55 mM increased the vegetative growth of the plants. Additionally, an assay performed with lettuce plants grown in a hydroponic system to which methionine was applied through the roots showed that this application notably improved the performance of the crop [11]. Likewise, in an experiment with tomato plants grown in a hydroponic system and under salt stress, Met treatment mitigated the negative effects of salinity [23]. Lastly, these positive effects of Met application were also previously verified in caupi pea plants under drought conditions [31], and in bitter melon plants as well [32].
The combined application of Tyr + Lys + Met did not alter the plant response with respect to individual AA application. Therefore, the combination of these three AAs had a neutral effect, as it did not improve or worsen the outcome. However, synergistic or antagonistic interactions between the AAs have been observed in other studies. For example, in a study conducted with tomato plants grown under saline treatments, an antagonistic effect was observed when methionine and arginine were applied together, as the combined application of these two amino acids had negative effects on the plants. In the same study, a synergistic effect was found for the simultaneous application of proline and glutamic acid, as the response of the plants was better than when these AAs were applied separately [23]. However, in a study conducted with lettuce plants, where glycine and glutamine were applied individually or combined, no important differences were found in the results [33]. In a treatment with soy seeds, a synergistic effect with the AAs cysteine, phenylalanine, glycine, and glutamate was observed when these were applied as a mixture [12].
Individual application of Tyr, Lys, and Met stimulated physiological processes of plants, as observed in the gas exchange data (Table 2), in which A CO2 increased by 30% with respect to the control plants. AA application also downregulated stomatal conductance and leaf transpiration (data not shown). The ratio between A CO2 and leaf transpiration (a ratio known as water use efficiency, WUE) increased with the application of AA treatments with respect to the control plants, and this parameter was higher in plants from the Tyrosine treatment. High WUE values are beneficial to crops, as they indicate that for the same quantity of water transpired, they are able to produce a greater biomass, allowing the plants to regulate and save water in a sustainable manner during their development cycle. Thus, with these data, it can be concluded that the applications of Tyr, Lys, and Met induced an increase in WUE in plants by regulating non-stomatal factors (even when the concentration of substomatal CO 2 decreased and A CO2 increased) and stomatal factors (regulation of stomata opening and closing).
In photosynthesis, there are various non-stomatal factors, such as the photochemical efficiency of photosystem II, the carboxylation efficiency, the activity of photosynthetic enzymes, and the chlorophyll concentration [34,35]. The results of chlorophyll fluorescence in our assays allowed us to conclude that application of these AAs induced an increase in A CO2 through the biochemical reactions of C metabolism and not through the light-dependent processes of photosynthesis, as changes were not observed in the qP, Fv'/Fm', and Φ PSII parameters ( Figure 5). However, these AAs increased the concentration of chlorophylls, which could be the reason why A CO2 increased. Some evidence of A CO2 and g s regulation has already been published by Kong et al. [36], but the possible mechanisms involved in these effects are still unknown. Bakry et al. [37] observed that exogenous application of Tyr in concentrations of 0.27, 0.55, and 0.83 mM improved growth by increasing the photosynthetic pigments of the plants. Likewise, in an assay with lettuce plants to which Met was applied through the roots in a hydroponic system, better results were obtained with the lowest doses (0.0013 and 1.3·10 −7 mM), highlighting that these doses improved the photosynthetic activity and also showed an increase in performance. It was also demonstrated that higher doses damaged the photosynthesis machinery [11]. On the other hand, in an experiment conducted with bean plants irrigated with seawater, to which a commercial product based on amino acids was applied at different concentrations, a relationship between the AAs and photosynthetic production was established [38]. However, foliar application of two commercial products, 'Kadoustim' ® and 'Humi Fort' ® , composed of AAs, did not show an increase in photosynthetic parameters [39]. In conclusion, these studies confirm that the response mechanisms involved are still unknown, and this could be due to the great variability existing between amino acids, doses, crops to which these are applied, timing, and manner of application.

Conclusions
In this study, foliar applications of Tyr, Lys, and Met at concentrations of 15 mM had beneficial effects on tomato plants, as an increase in vegetative growth was observed compared to the control plants. This was due to the treatments stimulating ACO2 and increasing the efficient use of nutrients; therefore, despite the reduced concentration of minerals, this did not limit the physiological processes. Another interesting result from the experiments was that treatments with AAs increased the water use efficiency as a result of the increase in ACO2 as well as a reduction in gs. Thus, application of these AAs could be useful for agronomic strategies directed towards saving water and, therefore, will support the practice of sustainable agriculture. Our metabolic study indicated that greater vegetative growth of the plants treated with AAs could be due to (i) increased proline concentration in the plants, which is an AA related to the processes of osmotic adjustment and protection against adverse environmental factors; (ii) the greater increase in ACO2 rate, promoting additional carbon in the shape of fructose and glucose, which are destined towards the processes of glycolysis and the Krebs cycle to obtain a greater reducing power and energy, as well as more C molecules needed for The white color indicates that no significant difference between the control and the corresponding AA treatment was found.
In the ionomic study, all the plants, independently of the treatment, had leaf nutrient concentrations within normal ranges for tomato crops (Table 3; [40]). However, this concentration was lower (Ca, K, Mg, and P) in the plants treated with AAs. This result indicates that the AA treatments inhibited absorption of the nutrients; however, this result may be related with a dilution affect that may occur when shoot mass increases by 46%. The decrease in concentration of minerals, without being affected by the growth parameters or A CO2 , could also indicate that application of these AAs increases the nutrient efficiency, meaning that for a given lower concentration of nutrients in the leaves, these plants are able to increase their growth and stimulate their physiological processes, as compared to plants not treated with AAs. A different model from the one stated previously was only observed for N and Fe. In the case of N, this could be due to the extra supply of N given to the plants treated with the AAs. The relationship between the AAs and N assimilation is more related to the general levels found in the cell than to the accumulation of a specific AA [41]. On the other hand, more recent studies have shown that Ca is an essential element for plants, and it also plays an important role in N assimilation. Lastly, an appropriate supply of Ca can considerably increase nitrate transport and also increase the concentration of amino acids in leaves and roots [42].
Our metabolic study showed that application of AAs, independent of treatment, changed the aminogram of the leaves: the concentration of proline increased and that of glutamine decreased. This indicates that for our plants treated with exogenous AAs, proline was the main amino acid, and glutamine increased its concentration near aspartate values (Table 4, Figure S1). Proline is an amino acid that plays multiple roles in plants, including osmotic adjustment, stabilization of membranes and proteins, reserving carbon and nitrogen, inhibition of reactive oxygen species (ROS), etc. [43,44]. The increase in this AA due to exogenous applications of Tyr, Lys, and Met led to a better response of the plants to environmental conditions. Although the plants were not subjected to abiotic stress during the experiment, the temperatures commonly reached inside a greenhouse located in the Mediterranean region, especially at midday, could affect their development, so that an increase in endogenous proline could have palliated the negative effects caused by environmental factors. Along the same line, we found an increase in isoleucine in all AA treatments and an increase in tyrosine in the treatments with Tyr and Tyr + Lys + Met. Tyrosine and isoleucine have specific functions in crops. It is known that branched amino acids such as isoleucine play a fundamental role in a plant's tolerance to stresses, and this could be due to the adjustment of intracellular pH, accumulation of osmolytes, and the detoxification of reactive oxygen species (ROS), xenobiotics, and heavy metals. This AA could also be an alternative source of C under stress, given that branched amino acids promote their own degradation during specific times, such as seed germination and senescence, or at times when sugars are lacking [45]. This is also corroborated by Li et al. [46], attesting that accumulation of isoleucine could palliate the negative effects of salinity by improving glycolysis. Among the many functions of Tyr, pollinator attraction, defense, and a greater structural support are attributed to the greater assimilation of lignin [19]. This amino acid is synthesized in the shikimate route, more specifically starting with chorismate. This route is a starting point for the biosynthesis of carbocyclic aromatic compounds, and it is also a precursor for other AAs, such as phenylalanine and tryptophan. In addition, chorismate is involved in the synthesis of vitamins K and E, folate, plastoquinones, and co-enzyme Q [21].
NMR analysis identified the compounds malate, citrate, fumarate, and formate as the major organic acids. These compounds participate in the tricarboxylic acid cycle (citrate, fumarate, and malate) and the pentose phosphate pathway (formate). The organic acid profile, according to their concentrations, followed the order malate > citrate > fumarate = formate ( Figure S2). This order did not change with the application of different AAs, although the concentration of these compounds was affected by the AA treatments with respect to the control plants. These data indicate that application of AAs modifies the metabolic processes in the plants. A generalized response was observed with the application of Tyr, Lys, Met, and their mixture, based on an increase in the concentration of malate and a decrease in the concentration of citrate and fumarate in all AA treatments. The increase in malate in the Krebs cycle (major acid) could indicate that treatments with AAs induced a greater influx of C to the glycolysis cycle, which could generate greater oxidation power (redox potential) and energy (NADH, FADH2, and ATP) with respect to the plants that were not treated with the AAs [47].
The major sugars in tomato plant leaves are fructose, sucrose, and glucose ( Figure 2). In control plants, their concentrations decreased in the following order: fructose > sucrose > glucose. However, when the AAs were applied, the order changed to fructose > glucose > sucrose. Therefore, the AAs could have promoted a greater influx of C into the plant, thereby increasing the concentrations of fructose and glucose, so that the increase in glucose shuttled more C to the glycolysis reactions as compared to the synthesis of sucrose. This would result in the plant having greater redox potential and more C and energy available for its growth [47].
GABA is a non-proteinogenic amino acid that can temporarily store nitrogen, thereby reducing excess accumulation of amines produced during photorespiration and maintaining appropriate pH. The synthesis of GABA is catalyzed by the enzyme glutamate decarboxylase (GAD) and liberates CO 2 under drought stress conditions, which promotes correct functioning of the Calvin cycle with closed or semi-closed stomata. Thus, pressure on the photosynthetic electron chain is relieved, the concentration of ROS decreases, and photodamage to leaf tissues is curbed. GABA can act as an osmolyte with no toxic effects in equilibrating the water potential during cellular dehydration, and it works as an antioxidant to stabilize and protect thylakoids and macromolecules [48].

Conclusions
In this study, foliar applications of Tyr, Lys, and Met at concentrations of 15 mM had beneficial effects on tomato plants, as an increase in vegetative growth was observed compared to the control plants. This was due to the treatments stimulating A CO2 and increasing the efficient use of nutrients; therefore, despite the reduced concentration of minerals, this did not limit the physiological processes. Another interesting result from the experiments was that treatments with AAs increased the water use efficiency as a result of the increase in A CO2 as well as a reduction in g s . Thus, application of these AAs could be useful for agronomic strategies directed towards saving water and, therefore, will support the practice of sustainable agriculture. Our metabolic study indicated that greater vegetative growth of the plants treated with AAs could be due to (i) increased proline concentration in the plants, which is an AA related to the processes of osmotic adjustment and protection against adverse environmental factors; (ii) the greater increase in A CO2 rate, promoting additional carbon in the shape of fructose and glucose, which are destined towards the processes of glycolysis and the Krebs cycle to obtain a greater reducing power and energy, as well as more C molecules needed for the growth of the plants. To link metabolic profiles with vegetative growth, artificial intelligence will be needed in the future to estimate how plant growth is initiated with these metabolic profiles.

Conflicts of Interest:
The authors declare no conflict of interest.