Stevia rebaudiana (Bertoni) as a Multifunctional and Sustainable Crop for the Mediterranean Climate

: Stevia rebaudiana (Bertoni) is a promising medicinal and aromatic plant for Mediterranean agroecosystems given its positive agronomic attributes and interesting quality features. It has both food and pharmaceutical applications, since its leaves contain sweet-tasting steviol glycosides (SVglys) and bioactive compounds, such as phenolics, ﬂavonoids, and vitamins. We evaluated the agronomic and qualitative performances of nine stevia genotypes cultivated, in open ﬁeld conditions, for two consecutive years under the Mediterranean climate of central Italy. Growth, biomass production, and accumulation of bioactive compounds (SVglys, polyphenols, and their related antioxidant activities) were evaluated, considering the effect of harvest time and crop age (ﬁrst and second year of cultivation). The results showed high variability among genotypes in terms of both morpho-productive and phytochemical characteristics. In general, greater leaf dry yields, polyphenol accumulation, and antioxidant activities were found in the second year of cultivation, harvesting the plants in full vegetative growth. On the other hand, total SVglys leaf content reached the highest values in the ﬁrst year when plants were at the beginning of the reproductive phase. On the other hand, although the SVglys proﬁle (Rubusoside, Dulcoside A, Stevioside, Rebaudioside A, C, D, E, and M) remained stable over harvest times, it differed signiﬁcantly depending on the crop age and genotype. Our ﬁndings provide useful information on the inﬂuence of crop age and harvest time in deﬁning quanti-qualitative traits in stevia, with PL, SL, BR5, and SW30 being the best performing genotypes and thus suitable for breeding programs. Our study highlighted that stevia, in the tested environment, represents a promising semi-perennial crop which offers new solutions in terms of cropping system diversiﬁcation and marketing opportunities.


Introduction
In many Mediterranean areas, high evapotranspiration rates, increased precipitation variability, and intense summer drought are the main environmental constraints to agricultural management and crop production. Under such conditions, sustainable agricultural practices and the introduction of new crops that are able to diversify cropping systems and to mitigate climate changes are very promising strategies for addressing cropping system sustainability.
Compared to continuous cropping systems, crop diversification involves better utilization of land resources, lower risks from pests and diseases, and greater yield stability. In this context, medicinal and aromatic plants (MAPs) can be included in low-input productive systems: they can enhance the multifunctionality traits of the agricultural sector, produce safe final products, respect the environment, and promote rural areas. Various MAPs find Weather parameters (daily minimum, T min , maximum, T max , and mean temperatures, T mean , and cumulative rainfall) were recorded by an automated weather station nearby the experimental site, from the beginning of vegetative plant development (June) to the full plant flowering (October) in each year of cultivation. Throughout the experiment, total rainfalls were 164.2 and 233.2 mm, in 2018 and 2019, respectively, with a mean average temperature of 19.0 • C and 16.4 • C in 2018 and 2019, respectively. Mean average T max and T min temperatures did not differ notably between 2018 and 2019 ( Table 2). Before beginning the experiment, soil samples were collected at 0 to 30 cm depth in order to evaluate the physical and chemical characteristics of the soil ( Table 3). The soil was loam with a sub-alkaline reaction, a good content of total nitrogen, organic matter, and exchangeable potassium, and with a low level of available phosphorus and salinity.

Crop Sampling and Agronomic Measurements
Phenological observations, from the beginning of vegetative plant development until full plant flowering, were performed in the ratooning crop, starting from the regrowth of new shoots, per each genotype and in each plot. For every observation, six plants per genotype were randomly selected in the field, with 54 total plants observed. The length of both the vegetative and reproductive phases was evaluated using accumulated thermal time ( • C/day) and the accumulated growing degree days (GDDs) were calculated daily according to the equation (1) presented by McMaster and Wilhelm [24], where T max is the daily maximum air temperature, T min is the daily minimum air temperature, and T b is the base temperature, for which 10 • C was used [22]. For each growing season, two destructive samplings (namely HT1 and HT2) were carried out at the following phenological stages (according to the BBCH scale by Le Bihan et al. [25]): HT1-Stage 48 = about 80% of final leaf biomass is developed; and HT2-Stage 55 = 50% of apex leaves are differentiated and present inflorescence, but flower buds are still closed.
HT1 corresponded to the full vegetative growth which was reached between 27 July-2 August in the 1st year and between 15-17 July in the 2nd year, depending on the genotype. HT2 corresponded to the beginning of flowering, occurred in early September in the 1st year (10-13 September 2018) and between the end of August and the beginning of September in the 2nd year (from 20 August to 2-3 September 2019). Samplings were manually performed, collecting six plants per genotype (1 sample = 1 plant). After each harvest, plant height (cm), basal stem diameter (mm), branching (n . stem/plant), fresh and dry weight of leaves and stems (g plant −1 ), harvest index (HI), and specific leaf weight (LSW) were measured. Leaves of each sample were air-dried in a ventilated oven from 30 to 40 • C until constant weight. Dry leaves were ground to a fine powder by a laboratory mill (Grindomix GM 200, Retsch, Pedrengo (BG), Italy) and stored until the subsequent analyses. The harvest index (HI), which represents the plant's efficiency at producing leaves, was calculated as the ratio between leaf dry yield (g plant −1 ) and total aboveground biomass yield (leaves plus stems). Specific leaf weight (SLW) was calculated as the ratio between dry leaf weight and leaf area (mg cm −2 ). The leaf area was measured by collecting and subjecting to a color scan, two apical leaves, two middle leaves, and two basal leaves. The leaf images were then analyzed with ImageJ (Fiji Particle Analysis plug-in). hydrate (3,4,5-trihydroxybenzoic acid), 2,4,6-tri(2-pyridyl)-triazine (TPTZ), trizma acetate, Folin-Ciocalteu reagent, Trolox (6-hydroxy-2,5,7,8-tetramethylchroman-2-carboxylic acid), sodium carbonate, and ferric chloride were obtained from Sigma-Aldrich Chemical Co. (Milan, Italy). All chemicals used in the present study, including solvents, were of analytical grade.

Sample Extraction
A total of 0.1 g of leaf powder per sample was extracted with 10 mL of 70% (v/v) EtOH, and sonicated for 30 min at 60 • C. At the end of the sonication, the extracts were centrifuged (3500 rpm for 10 min) and filtered with a syringe filter (Ø 0.45 µm) to remove any suspended material. The extracts obtained were stored at 4 • C until subsequent analyses.

Steviol Glycosides Determination
The extraction procedure and the determination of steviol glycosides (SVglys) were carried out following Zimmerman et al. [26]. Steviol glycosides were analyzed using a Jasco PU980 HPLC system (JASCO Benelux B.V., Utrecht, Netherlands) coupled with a UV-visible wavelength detector. A hydrophilic column (Luna HILIC 200A, Phenomenex Inc., Torrance, CA, USA), 5 µm, 250 mm × 4.6 mm (Phenomenex Inc., Torrance, CA, USA) in combination with the corresponding pre-column (4 × 3.0 mm) was used. UV detection was carried out at 205 nm at room temperature with a flow rate of 0.68 mL/min and a run time of 20 min. Separation was achieved in acetonitrile/water (80:20) as isocratic mobile phase at pH 3.6 regulated with acetic acid. Chromatograms were acquired online, and data were collected using a Jasco interface (Hercules 2000 Interface Chromatography). Steviol glycosides were quantified using authentic standards, through calibration curves (0.005-1.00 g L −1 ), obtained from standard mixtures containing Rubusoside (Rub), Dulcoside A (Dulc A), Stevioside (Stev), and Rebaudiosides A, C, D, E, and M (Reb A, C, D, E, and M).

Analysis of Total Phenols and Flavonoids
Total phenols were determined using the Folin-Ciocalteu method according to Dewanto et al. [27] and expressed as gallic acid equivalents (mg GA g −1 dry leaf). This method involved the reduction of Folin-Ciocalteu reagent by phenolic compounds, with a blue complex formation determined at 765 nm by UV-Vis spectrophotometer (Varian Cary 1E, Palo Alto, CA, USA). Total flavonoids were determined using the aluminum trichloride method according to Jia et al. [28] and expressed as catechin equivalents (mg CE g −1 dry leaf). The flavonoids-aluminum reaction created a pink complex formation measured at 510 nm using a UV-Vis spectrophotometer (Varian Cary 1E, Palo Alto, CA, USA).

Ferric Reducing Antioxidant Power (FRAP) Assay and Free Radical-Scavenging Assay
The determination of ferric reducing antioxidant power and the free radical-scavenging activity (FRAP and DPPH assay) of stevia leaf extracts followed Tavarini et al. [29]. The FRAP method is based on the ability of the antioxidant compounds to reduce Fe 3+ to Fe 2+ which, in the presence of TPTZ (2,4,6-tris(2-pyridyl)-s-triazine) led to a blue complex formation (Fe 2+ -TPTZ), measured at 593 nm using a UV-Vis spectrophotometer (Varian Cary 1E, Palo Alto, CA, USA). The DPPH assay is based on the reducing activity of the antioxidant molecules against the 1,1-diphenyl-2-picryl-hydrazil (DPPH) radical which was characterized by a purple red color. The extent of the disappearance of DPPH is directly proportional to the amount of antioxidant present in the reaction measured at 517 nm using a UV-Vis spectrophotometer (Varian Cary 1E, Palo Alto, CA, USA). Total antioxidant activity and free radical-scavenging capacity were expressed as trolox equivalents (mmol TE g −1 dry leaf).

Statistical Analyses
All data were subjected to analysis of variance (ANOVA) using GraphPad Prism v. 8.0.2 (GraphPad Software, Inc., La Jolla, CA, USA). A three-way ANOVA analysis Agriculture 2021, 11, 123 6 of 17 was conducted to assess the effect of genotype (G), crop age (CA), harvest time (HT), and their reciprocal interactions, on the agronomic characteristics and on total phenols and flavonoids, total SVglys and antioxidant activities (FRAP and DPPH). Means were separated on the basis of the least significant difference (LSD) only when the ANOVA F test showed significance at 0.05 or 0.01 probability level.
Hierarchical cluster (HC) and principal component (PC) analyses were performed on (i) the total content of phenols, flavonoids, and steviol glycosides, and (ii) on the individual steviol glycoside (Dulcoside A; Rebaudiosides A, C, D, E and M; Rubusoside, Stevioside) concentrations with JMP ® Pro 13.2.1 (SAS Institute Inc., Cary, NC, USA). As unsupervised methods, the groups of samples obtained with hierarchical cluster analysis (HCA) and PCA can be observed even when there are no reference samples that can be used as a training set to establish the model. For both observation groups, the hierarchical cluster analysis (HCA) was conducted on the normalized average values, with Ward's algorithm, using Euclidean distances as a measure of similarity among the samples. In addition, principal component analyses (PCA) were carried out in order to reduce the dimensionality of the multivariate data of the matrix, whilst preserving most of the variance [30]. For the PCA of the total content of phenols, flavonoids, and steviol glycosides, a 35 × 3 (35 samples, 3 metabolite contents, 105 total data) dimensional matrix was used. The score plot obtained was defined by a PC1 and a PC2 covering 78.3 and 21.3% of the variance, respectively, for a total explained variance of 99.6%. For the PCA of the total content of individual steviol glycoside (Dulcoside A; Rebaudiosides A, C, E, M, and D; Rubusoside, Stevioside) concentrations, a 35 × 8 (35 samples, 8 steviol glycosides, 280 total data) dimensional matrix was used. The score plot obtained was defined by a PC1 and a PC2 covering 68.2 and 19.3% of the variance, respectively, for a total explained variance of 87.5%. Genotype (G), crop age (CA), harvest time (HT) and their reciprocal interactions significantly affected biometric and productive traits, except for (i) branching, which did not vary depending on CA, (ii) basal stem diameter, in relation to G × HT interaction and (iii) plant height and SLW which were not affected by CA×HT (Table 4).

Growth, Biometric, and Productive Measurements
Taking into account the genetic effect (G), plant height ranged from 47.31 to 58.19 cm, reaching the highest values in RG and NU genotypes and the lowest values in SL, PL, and SW30. On the other hand, SW30 was characterized by the greatest stem basal diam- not vary depending on CA, (ii) basal stem diameter, in relation to G × HT interaction and (iii) plant height and SLW which were not affected by CA×HT (Table 4).
Taking into account the genetic effect (G), plant height ranged from 47.31 to 58.19 cm, reaching the highest values in RG and NU genotypes and the lowest values in SL, PL, and SW30. On the other hand, SW30 was characterized by the greatest stem basal diameter (8.80 mm), while BR16 showed the lowest value (5.06 mm), although along with RG, BR16 had the greatest branching than the other genotypes (PL, SL, NU, CO, BR5, SW30, and BR1) ( Table 4).
Regarding leaf dry yield per plant, NU and BR1 were the most productive, while BR16 appeared to be the least productive genotype. The other genotypes showed medium leaf dry yields. Harvest index ranged from 73.51% (SL) to 52.79% (BR16), reflecting the plant architecture, in terms of height and branching, and leaf yield. PL and SL exhibited the highest specific leaf weight (SLW), while RG and NU showed the lowest weights.
In terms of crop age, all biometric and productive measurements showed significant increases in the second year of cultivation, except for the stem basal diameter, which was significantly smaller in 2-year-old plants, and branching, which did not vary in the two years of cultivation (Table 4).
Additionally, harvest time played a key role in defining the agronomic responses of stevia genotypes. In general, maximum plant height, stem basal diameter, branching and leaf dry yield increased from the point of full vegetative growth to the beginning of flowering (Table 4). On the other hand, the highest values of specific leaf weight and harvest index were recorded during the full vegetative growth, with a significant decrease in the subsequent harvest.

Phytochemical Evaluation and Antioxidant Activities
Interestingly, compared to the other genotypes, BR16, RG, and NU exhibited the lowest values of all secondary metabolites (total phenols, total flavonoids and total SVglys), as well as the lowest in vitro antioxidant activities (Table 5). Crop age (CA) positively influenced the phytochemical characteristics of stevia leaf extracts, with an increase in their values in the plants in the second year after transplanting. Only the total SVglys showed a significantly decrease passing from the first to the second year of cultivation. Taking into account harvest time, total phenols, total flavonoids, and the antioxidant activities decreased from the vegetative phase (July) to the beginning of the reproductive one (September). On the other hand, an opposite trend, with a significant increase from the 1st to the 2nd harvest time, was observed for total SVglys, confirming that the development of the first flower buds corresponds to the time of major accumulation of SVglys in stevia plants. Table 6 reports the SVglys profile and related statistical significance. A significant effect of genotype was observed for almost all steviol glycosides, with the exception of Dulcoside A. Interestingly, the identified compounds were not present in all genotypes. This difference was evident for Reb M which was not detected in RG, SL, NU, and CO, and for Reb D which was not detected in NU and CO. The reduction in SVglys content observed in plants at the second year of age was also seen in the most represented compounds such as Stev, Reb C, Reb A, and Reb C. On the other hand, Rubusoside increased with crop age, while no effect was observed for Reb M and Reb D. With regard to harvest time, as already noted for the total content, each steviol glycoside significantly increased from the vegetative phase to the beginning of flowering, except for Rubusoside which remained stable in both harvest times.
HCA and PCA were carried out for total phenols, flavonoids, and SVglys content. The two-way dendrogram of the HCA is reported in Figure 2. The first macro-cluster (red) was grouped by itself, while the second comprised three sub-clusters (green, blue, and yellow). Among the metabolites analyzed, the total contents of phenols and flavonoids were clustered together, while the total SVglys concentration was grouped by itself, thus highlighting a higher degree of dissimilarity based on this parameter. Based on their total SVglys concentration, samples in the yellow cluster, all from 2018, shared the highest SVglys content; samples of the red cluster, all from 2019, were, instead, characterized by intermediate concentrations of SVglys; the lowest concentrations of these metabolites were, instead, common to samples of the blue and green clusters, especially those from 2019. As evidenced by their grouping in a common macro-cluster, the total phenols and total flavonoids showed a common quantitative distribution among the samples. In fact, higher concentrations of both these chemical classes were evidenced for samples of the red cluster, all from 2019; also, intermediate concentrations were evidenced for samples from 2019 in the green cluster; finally, samples belonging to the blue and yellow clusters exhibited the lowest contents of these compounds. With the exception of the green cluster, the samples were homogeneously distributed in the groups based on the harvest time.
Some genotypes showed very similar total abundances of the analyzed compounds between the two harvest times ( Among the other groups, the higher content of SVglys detected in the yellow samples meant that they were plotted in the right quadrants (PC1 > 0), whereas all the other samples were in the left quadrants (PC1 < 0). The same multivariate analyses were conducted on the SVglys profile of all the genotypes studied (Figures 4 and 5). The two-way dendrogram (Figure 4) of the HCA on the compounds showed two main groups driving the dissimilarities: Rubusoside, Stevioside, Rebaudiosides A, M, and D for macro-cluster 1; Rebaudioside E, Rebaudioside C, and Dulcoside A for macro-cluster 2. The samples in the dendrogram (Figure 4) were distributed in two macro-clusters: the first comprised the red and green sub-groups; the second comprised the blue and yellow sub-groups, of which the latter only contained one sample (SW30_2_2018). The second macro-cluster was mainly composed of 2018 samples, while the first macro-cluster comprised all the 2019 samples, as well as a few of the 2018 ones ( Figure 4).
Compared to the HCA of the total compounds, the dendrogram of the SVglys profile suggested that the differences between the samples were less due to the year and more to the steviol glycosides pool produced by each genotype (i.e., all the PL samples are in the red cluster; all the BR16 samples are in the green cluster). In addition, the SVglys pool of samples RG and NU appeared very similar, thus suggesting the proximity of these two genotypes; the same seemed true for BR5, BR1, and BR16 ( Figure 4). Regarding SVglys profile, the influence of the year on the sample distribution was also confirmed by the PCA. The score plot is reported in Figure 5

Discussion
We assessed the agronomic and phytochemical performances of nine S. rebaudiana (Bertoni) genotypes for two consecutive growing seasons. The aim was to select the best productive genetic resources to develop new and improved cultivars for the Mediterranean area. To improve the competitiveness of stevia production as a semi-perennial crop in this environment, it is important to produce higher-performing crops in terms of yield and quality. Therefore, the identification of new varieties/cultivars of S. rebaudiana with long stand duration, consistently high and stable leaf and SVglys yields, as well as a high level of other beneficial substances, is the top priority for the successful introduction of stevia into Mediterranean cropping systems.
Our results showed that some of the genotypes tested were characterized by a longer vegetative period before flowering, as demonstrated by the higher thermal requirements ( Figure 1). This behavior promotes both higher leaf yield and steviol glycoside concentration, which in turn determine higher SVglys yield. This is important given the time required for the plant to synthesize and accumulate the steviol glycosides in the leaves. In fact, the maximum accumulation of these compounds, which depends on several environmental, agronomic, and physiological factors, is reached when the plant moves from the budding phase to an initial flowering stage with less than 10% flowers [31]. Therefore, the greater the vegetative period, the greater the accumulation of steviol

Discussion
We assessed the agronomic and phytochemical performances of nine S. rebaudiana (Bertoni) genotypes for two consecutive growing seasons. The aim was to select the best productive genetic resources to develop new and improved cultivars for the Mediterranean area. To improve the competitiveness of stevia production as a semi-perennial crop in this environment, it is important to produce higher-performing crops in terms of yield and quality. Therefore, the identification of new varieties/cultivars of S. rebaudiana with long stand duration, consistently high and stable leaf and SVglys yields, as well as a high level of other beneficial substances, is the top priority for the successful introduction of stevia into Mediterranean cropping systems.
Our results showed that some of the genotypes tested were characterized by a longer vegetative period before flowering, as demonstrated by the higher thermal requirements (Figure 1). This behavior promotes both higher leaf yield and steviol glycoside concentration, which in turn determine higher SVglys yield. This is important given the time required for the plant to synthesize and accumulate the steviol glycosides in the leaves. In fact, the maximum accumulation of these compounds, which depends on several environmental, agronomic, and physiological factors, is reached when the plant moves from the budding phase to an initial flowering stage with less than 10% flowers [31]. Therefore, the greater the vegetative period, the greater the accumulation of steviol glycosides in the leaves. When the plant starts to flower, nutrients accumulate in the reproductive organs and, as a result, vegetative growth declines.
Given that the leaves are the commercially important part of stevia, a delay in flowering can enhance vegetative growth and economic yield. In this regard, Ceunen and Geuns [31] observed large amounts of SVglys within the upper leaves during the budding phase and, as more and more buds become flowers, leaves contained lower amounts of SVglys. This is important for the choice of best harvesting time, since the maximum SVglys concentration is reached in the physiological stage of flowering-bud formation [32,33]. Flowering and, in general, plant growth are also affected by the photoperiod and temperature. Stevia is, in fact, a short-day plant with a critical photoperiod of between 12-13 h of day length [34,35]. Regarding air temperature, in our environment vegetative growth was lower when the maximum day temperature was below 10 • C or over 35 • C.
Overall, our results showed that genotype, crop age and harvest time represented key pre-harvest factors for defining the morphological, phenological, and quanti-qualitative traits in stevia.
Regarding differences in plant morphology and canopy architecture, we confirmed previous findings [21,36,37], which reported a very high variability, due to genetic characteristics, has been reported for plant height, basal stem diameter, branching, specific leaf weight, and leaf yield. In order to investigate this great variability among stevia genotypes and populations, some authors [21,22] have observed that increased yields were correlated to a high leaf area index (LAI), which is responsible for a greater light interception and, in turn, for a higher leaf photosynthesis.
We assessed the specific leaf weight (SLW), a leaf thickness index that is positively correlated with leaf photosynthesis: leaves with high SLW values are thicker and generally have a higher chlorophyll density per area unit (µg cm −2 ) and, therefore, greater photosynthetic capacities than thinner leaves [38]. SLW significantly changed depending on genotype, crop age and harvest time. In particular, in 2-year-old plants, an increase in SLW was accompanied by increased leaf dry yields, suggesting that, thanks to a greater SLW leaf photosynthesis may have been enhanced with a consequent increase in crop yield [39].
Plants with a high SLW have leaves with a low surface/volume ratio, which is more efficient in terms of water use [40]. This ratio plays an important role in leaf functioning and is related to the strategies for acquiring and using the energy resources [41], as well as a tool to evaluate plant productivity [42] with a view to sustainable agriculture. For a better comprehension of the relationships between morphotype and light interception and photosynthetic activity, future studies are needed in order to improve leaf yield in S. rebaudiana.
Regarding stevia productivity, we found that the least productive genotypes are those with more branching, as reported by Tateo et al. [43]. In addition, our findings suggested a positive relation between plant height and leaf dry yield in agreement with previous studies [36,37]. Finally, the aerial biomass variability among genotypes exhibited high heritability, and the differences can also be partly attributed to crop age, pedo-climatic conditions of the cultivation site and plant development stage, as already observed in previous studies [19,44,45].
In our environment, stevia productivity significantly increased from the first to the second year of cultivation, with a very high winter survival rate. This suggests that stevia could grow as a semi-perennial crop in central Italy. This could bring positive agronomic advantages, such as, reductions in soil erosion and nutrient leaching, C soil sequestration, which in turn have important implications in multifunctional and sustainable agriculture.
Our findings confirmed previous results [5,6] in which in the temperate area of central Italy, the cultivation of stevia is long-term (5-7 years), with a vegetative period from April-May to September-October. In the tested climatic conditions, during winter, above-ground parts of the plant became dry, and there was regrowth through new shoots in the buried rhizome the next spring, thus producing a new crop without replanting. On the other hand, at higher latitudes as in central Europe (such as Belgium and Germany) and Canada [46][47][48], with cooler winters, stevia behaves like an annual crop with annual replanting.
We thus believe that stevia can be cultivated as a promising new crop in the Mediterranean climate of central Italy, though the SVgly content and composition still need to be optimized.
We found that the content and composition in bioactive compounds varied greatly depending on the genotype, crop age and harvest time. Two-year-old plants generally exhibited the highest content of total phenols and flavonoids, as well as the antioxidant activities; conversely, total SVglys content peaked in the establishment year. The sampling carried out at the beginning of the reproductive phase only improved total SVglys, again confirming that, in stevia, plant development is a primary factor in influencing the biosynthetic pathway.
Harvesting when plants were in full vegetative development (July harvest) maximized the content of polyphenol compounds and improved the antioxidant activities. Similar results have been obtained in a previous work [49] highlighting how the identification of the optimal harvest time was able to maximize the bioactive compounds of interest and, consequently, the health-promoting properties of stevia leaves. However, since the biosynthesis of secondary metabolites is a plant defense mechanism against biotic and abiotic stresses, through many physiological, biochemical, and molecular changes in plant metabolism, it is not always easy to identify the optimal time in which the different metabolites are maximally expressed.
The accumulation of polyphenols during the beginning of flowering stage could be related to the ecological roles of these compounds, such as intensifying antifungal defenses and attracting pollinators [32]. In addition, since different interactions among environmental and agronomic factors may occur, it is difficult to select individual stimuli that can influence a single metabolic pathway.
Of the secondary metabolites that have been synthesized and accumulated in stevia leaves, the presence of phenolic compounds is the subject of increasing interest because of their significant practical use for nutritional and medicinal applications. In fact, there are important implications for the growing market of natural stevia products [50], where they are employed as natural preservatives, thanks to their capacity to delay the oxidative degradation of lipids and to improve the shelf life of foods and beverages. They are also involved in the prevention of oxidative stress in humans, thanks to the hydroxyl groups in their molecules which have antioxidant anti-inflammatory properties [51,52].
Our findings provide new knowledge about the dissimilarity among genotypes, which was above all generated by the total SVglys content with respect to total phenols and flavonoids. Previous studies have investigated the effect of farming practices and genotype × environment interaction on stevia crop productivity, Reb A and Stev accumulation and on their reciprocal ratio (Reb A/Stev) [19,21,22,[53][54][55]. Our results showed that the differences in the steviol glycoside profile are less due to crop age and more to genetics, as suggested by the dendrogram of the SVglys profile. It is well-known that, while the SVglys content in the leaves can vary depending on environmental and agronomic factors [56,57], the SVglys qualitative profile remains quite stable, indicating the high genotypic determinism of this trait [19].
We thus believe that our findings are a good starting point to screen the best genotypes in terms of SVglys profile, taking into account that breeding programs for stevia are increasingly aimed at varieties with an optimal Reb A/Stev ratio, which is considered a good qualitative measure of sweetness, or with a high content of Reb M and Reb D, characterized by a very sweet taste with no liquorice-like taste. We found that genotypes, such as SL, BR5, PL, and SW30 accumulated significantly high amounts of Reb A, Reb M, and Reb D. Among these genotypes, SL and BR5 were characterized by relatively low content of Stevioside. Conversely, PL and SW30 exhibited high contents of Stevioside, which negatively affect the Reb A/Stev ratio, lowering it below 1.
In our study, we try to identify the stevia ideotype, deriving from the combination of morphological, productive, and phytochemical traits. Consequently, considering not only the SVglys profile, but also leaf yield, growth crop cycle and the composition of their whole phytocomplex (polyphenols, total SVglys and their related antioxidant activities), SL, BR5, PL, and SW30 seemed to be the best performing genotypes, in the given environment, and thus suitable as starting point for future breeding programs.

Conclusions
Our findings highlighted that stevia, in the tested environment, represents a promising semi-perennial crop, which can contribute to the diversification of traditional cropping systems, thus increasing their sustainability and, at the same time, generating functional products with high added value.
We revealed a strong dependence of growth, crop yield and quality on genotype, harvest time, and crop age. The high variability among genotypes highlighted the importance of morpho-productive and qualitative traits for the identification of the best stevia ideotype, for our environment and for targeted end uses. The best combination of morphological, productive, and phytochemical traits was observed for SL, BR5, PL, and SW30 genotypes.
In conclusion, identifying genotypes characterized by specific SVglys profile (high Reb A, Reb D, Reb M), together with high content of all bioactive compounds and satisfactory yields, is crucial to establish future breeding programs for this crop.