Characterization and Evaluation of Salvia hispanica L. and Salvia columbariae Benth. Varieties for Their Cultivation in Southwestern Germany

: Rising consumer attraction towards superfoods and the steadily increasing demand for healthy, environmentally sustainable, and regionally produced food products has sharpened the demand for chia. Over the course of 4 years, two early ﬂowering chia varieties belonging to Salvia hispanica L., and Salvia columbariae Benth. Species were identiﬁed to complete their phenological development and, therefore, able to reach maturity under a photoperiod > 12 h, thus enabling the cultivation of chia in central Europe—more speciﬁcally, in southwestern Germany—consistently for the ﬁrst time. Results obtained by the conducted ﬁeld trial in 2018 showed that chia seed yields and thousand-seed mass ranged from 284.13 to 643.99 kg ha − 1 and 0.92 to 1.36 g, respectively. Further, the statistical analyses showed that the protein content of the cultivated chia varieties ranged from 22.14 to 27.78%, the mucilage content varied from 10.35 to 20.66%, and the crude oil content amounted up to 28.00 and 31.73%. Fatty acid proﬁles were similar to previously reported data with α -Linolenic acid being the most prominent one, ranging from 60.40 to 65.87%, and we obtained ω 6: ω 3 ratios between 0.2 and 0.3. In conclusion, chia could represent a promising raw material from a nutritional point of view, while being able to diversify the local food basis of southwestern Germany.


Introduction
The human food supply relies on an increasingly dwindling number of species [1]. This dependence is especially critical in the context of climate change and the related environmental problems [2]. Diversifying agricultural production and therefore broadening the food base will inevitably become an essential task to adjust to future climate change [2]. In this context, the cultivation of chia species (Salvia hispanica L. and Salvia columbariae Benth.) could play a key role in meeting the steadily increasing demand for healthy, environmentally sustainable, and regionally produced food products in Europe [3,4]. Being able to provide food, feed, and industrial products/raw materials for the pharmaceutical and cosmetic industry would make this low-input, multi-purpose crop a valuable and beneficial addition to the present agricultural landscape [5].
Chia (Salvia hispanica L.), also known as the gold of the Aztecs, a summer annual herbaceous plant belonging to the Lamiacea family, is a crop of high economical value and has gained popularity in Europe in recent years due to its exceptional nutritional composition and the ever increasing demand

Preliminary Trial Part 1: Climatic Chamber
In order to select new chia varieties which could possibly be cultivated under the long day conditions of southwestern Germany, a climate chamber trial was conducted at the University of Hohenheim in 2015/2016 where, in total, 13 Salvia hispanica L. varieties and one Salvia columbariae Benth. variety were cultivated (Table 1). The pot experiment was arranged in a completely randomized block design with six replicates. In total, 360 pots (Ø 11 cm) were filled with two-thirds fertile loam topsoil (Filderlehm) and one-third sand. No additional fertilizer was applied. Sowing took place manually by hand on 4 December 2015. Two seeds were sown per pot. Sowing depth was about 1 cm. Plants were thinned to one plant per pot on 19 December 2015. Twelve trays of 15 pots were placed on top of the tables in each of the two climatic chambers.
One chamber was set at 16 h of illumination to mimic long day conditions whereas the other chamber was set at 8 h of illumination in order to mimic short day conditions. During the whole experiment, which lasted from 100 to 191 days from sowing to harvest (depending on the development of each variety), the plants were illuminated by six 400-watt high-pressure sodium lamps Son-T-Agro (Philips, Amsterdam, The Netherlands), which were mounted at a height of about 80 cm above the two tables. The temperature was first set to 30 • C during the illuminated hours and 25 • C during the non-illuminated hours. On 16th December 2015, the temperature settings were changed to 25 • C during the illuminated hours and 18 • C during the non-illuminated hours. The temperature settings were switched off on 14th March 2016. Temperatures cycled sinusoidally during the day, ranging from 27 to 35 • C [24]. The humidity was kept constant at about 58%.
Through the above-outlined climate chamber trial, two previously unknown chia varieties, able to produce inflorescences and mature under the mimicked long day conditions, could be selected to be cultivated under field conditions in subsequent years, namely Golden Chia and SALV66.

Preliminary Trial Part 2: Field Trials
In 2016 and 2017, single-row field trials took place at the experimental station Ihinger Hof and solely served the purpose of seed multiplication of SALV66 and Golden Chia, which were identified through the previously conducted climate chamber trial. In 2018, a field trial was set up to test the overall performance of the selected varieties shown in Figure 1 in comparison to G8, which has already been cultivated for several years at Ihinger Hof [22,23].  field trial was set up to test the overall performance of the selected varieties shown in Figure 1 in comparison to G8, which has already been cultivated for several years at Ihinger Hof [22,23].

Site Description
The field trial was conducted in 2018 at the experimental station Ihinger Hof (University of Hohenheim, Upper Neckarland, 48°44′40,70" N 8°55′26,36" E). Precipitation during the experimental period amounted to 208.5 mm. Mean temperature during the respective growing season was 17.6 °C. Meteorological data were obtained by the weather station at Ihinger Hof ( Figure 2). Growing degree days were calculated using the growing degree-day Equation (1) of McMaster and Wilhelm [25], where TMAX is defined as the daily maximum air temperature, TMIN is defined as the daily minimum air temperature, and TBASE as the crop base

Site Description
The field trial was conducted in 2018 at the experimental station Ihinger Hof (University of Hohenheim, Upper Neckarland, 48 • 44 40,70 N 8 • 55 26,36 E). Precipitation during the experimental period amounted to 208.5 mm. Mean temperature during the respective growing season was 17.6 • C. Meteorological data were obtained by the weather station at Ihinger Hof ( Figure 2). field trial was set up to test the overall performance of the selected varieties shown in Figure 1 in comparison to G8, which has already been cultivated for several years at Ihinger Hof [22,23].

Site Description
The field trial was conducted in 2018 at the experimental station Ihinger Hof (University of Hohenheim, Upper Neckarland, 48°44′40,70" N 8°55′26,36" E). Precipitation during the experimental period amounted to 208.5 mm. Mean temperature during the respective growing season was 17.6 °C . Meteorological data were obtained by the weather station at Ihinger Hof ( Figure 2). Growing degree days were calculated using the growing degree-day Equation (1) of McMaster and Wilhelm [25], where T MAX is defined as the daily maximum air temperature, T MIN is defined as the daily minimum air temperature, and T BASE as the crop base temperature. The crop base temperature of 10 • C according to Baginsky et al. [6] was subtracted from the daily average air temperature.

of 21
If the daily maximum air temperature was less than the crop base temperature, The soil on the trial site was characterized as Pelosol brown soil according to the IIUSS Working Group WRB [26]. The soil characteristics are shown in Table 2.

Trial Setup
The field trial was set up as a randomized complete block design with three replications for each variety with a single plot size of 4 × 15 m. Triticale (Triticosecale Tscherm.-Seys. ex Müntzing) was cultivated as a previous crop. Harrowing for seedbed preparation took place shortly before sowing with the help of a mechanical combination seed drill (Rau Landtechnik GmbH; Brigachtal, Germany) in order to loosen the topsoil layer and achieve a good crumbling. Prior to sowing, one pass of seedbed preparation was conducted with a roller (Güttler GmbH, Kirchheim/Teck, Germany) representing the device equipped with a 724 rotary harrow (Kuhn Machinen-Vertrieb GmbH, Schopsdorf, Germany) to a depth of 5 cm. Mechanical sowing took place on 28th May 2018 using a Fendt 207 (AGCO GmbH, Duluth, GA, USA) equipped with a 121 Deppe D82 Sower (Agrar Markt DEPPE GmbH, Rosdorf, Germany) and a roller (Güttler GmbH, Kirchheim/Teck, Germany). Chia was sown at a row spacing of 50 cm to a depth of 1 cm at a rate of 1.5 kg ha −1 .
A universal fertilizer spreader "UKS 230" (RAUCH Landmaschinenfabrik GmbH, Sinzheim, Germany), which was attached to a Massey Ferguson (AGCO GmbH, Duluth, GA, USA), was used to fertilize the plots with 20 kg N ha −1 , broadcasted as calcium ammonium nitrate 39 days after sowing. Manual weed control was performed 37 days after sowing on 3rd July 2018. To eliminate flea beetle infestation, 15 days after sowing, Bulldock (Nufarm Deutschland GmbH, Köln, Germany) was applied at a rate of 300 mL ha −1 with β-Cyfluthrin as an active ingredient dissolved in 300 L H 2 O ha −1 with the help of a spraying system (Dammann, Buxtehude-Hedendorf, Germany) mounted on a Fendt equipment device (AGCO GmbH, Duluth, GA, USA). In order to control biting and sucking insects, Biscaya (Bayer CropScience Deutschland GmbH, Langenfeld, Germany) was applied 54 and 68 days after sowing at a rate of 300 mL ha -1 with thiacloprid dissolved in 300 L H 2 O ha −1 .

Phenological and Agronomic Traits
In 2018, three plants per plot were recorded and assessed in detail every week starting nine days after sowing. Days until flowering and maturity were monitored. Single plant height of the matured plants was measured. Additionally, the number of non-senescent first-degree branches at the main stem, the number of inflorescences (>5 cm for Salvia hispanica L.), and the length and the width of the central axis inflorescence of three plants per plot were recorded. Finally, the number of seeds produced per single plant was counted, and the respective seed weight per single plant was assessed.

Yield Traits
Manual harvest of Golden Chia took place on 29th August 2018. Mechanical harvest of the G8 and SALV66 took place on 5th October 2018 with a plot combine Classic (Wintersteiger AG, Ried, Austria). The plot seed yield was determined on the basis of absolute dry matter (0% grain moisture) and thousand Agronomy 2020, 10, 2012 6 of 21 seed mass (TSM) was determined by weighing and counting (Contador, Pfeuffer GmbH, Kitzingen, Germany) one thousand seeds per plot in order to calculate the mean thousand seed mass.

Quality Traits
Quality traits and the following mentioned procedures were determined and amended according to Grimes et al. [22]. Mucilage extraction took place according to Muñoz et al. [27] and total seed nitrogen content was determined according to Dumas [28]. In order to estimate the crude protein content, the total seed nitrogen content values were multiplied by the conversion factor 5.71 [29]. The determination of the crude oil content took place according to the European Commission's Regulation 152/III H procedure B [30]. Fatty acid methyl esters (FAMEs) were generated by rapid saponification and esterification in order to identify the fatty acid profile [31][32][33].

Statistical Analysis
For agronomical traits including plant height, plant weight, and number of inflorescences, branches, and seeds, as well as the lengths of the main inflorescences, single-plant data of three plants per plot were available. The following linear mixed model was used to detect variety differences among these traits. Each trait was used as a response variable in the subsequent model: where y ijk is the observation taken on plant k in block j with variety i; µ is the constant term; b j is the fixed effect of the j-th block; τ i is the fixed effect of the i-th variety; p ij is the random effect of the ij-th plot with p ij ∼ N 0, σ 2 p ; and e ijk are the residual errors with e ijk ∼ N 0, σ 2 e . In the case of some variables, separate error variances for each variety were estimated if residuals plots suggested strong, variety-specific heterogeneity of variance. This was the case with the following traits: length of main inflorescence and number of branches. Variance components were estimated by the method of restricted maximum likelihood through the MIXED procedure in SAS v9.4 (SAS Institute, Cary, NC, USA). The model assumption normal distribution of residuals and homogeneity of variance were graphically assessed by the inspection of quantile-quantile plots for the former assumption and scatter plots of residuals versus predicted values for the latter assumption. Plots of studentized residuals were used. If necessary, response variables were transformed to meet assumptions. Number of inflorescences and number of seeds were log-transformed. Seed weight was transformed by the fourth root. The null hypothesis of equal variety effects was tested using sequential Wald-type F-tests. Denominator degrees of freedom were adjusted using the method of Kenward and Roger. This method additionally corrects the standard errors of means and differences. If the null hypothesis was rejected, the three variety means were compared by pairwise t-tests. Throughout all statistical analyses, a significance level of α = 5% was used. Estimates of traits, which were transformed for analysis, were back-transformed for presentation. Standard errors were back-transformed using the delta method.
For the quality traits mucilage content, protein content, oil content, fatty acids (palmitic, stearic, oleic, vaccenic, linoleic, and α-linolenic acids), polyunsaturated fatty acid to saturated fatty acid (PUFA:SFA) and omega-6 to omega-3 ratios, as well as for the agronomical traits yield and thousand seed mass, only one observation per plot was available. However, the number of plants per plot varied strongly. To adjust for varying number of plants, (a) the number of plants per m 2 were included as a weighting variable into statistical analysis and (b) the number of plants was, in addition, included as a covariate into the model. The following model was used: where y ij is the observation in block j with variety i, i.e., the plot observation; k ij is the corresponding number of plants per plot, which is used as a covariate and weighting variable; µ is the constant term; b j is the fixed effect of the j-th block; τ i is the fixed effect of the i-th variety; β 1 is the common slope in a regression on the covariate k ij ; β 2i is the variety-specific slope for the regression on the covariate; e ij are the residual errors with e ij ∼ N 0, k −1 ij σ 2 e . Parameter estimation, adjustment of denominator degrees of freedom, and the F-test followed the same method as described for agronomical traits. First of all, in the interaction of covariate and treatment, the factor (β 2i k ij ) was tested for significance. Only without such interaction, the covariate adjustment of the treatment variable was performed. Non-significant effects were removed from the model, and in case of a rejection of the null hypothesis of equal variety effects, the effects were compared by pairwise t-tests.

Results and Discussion
The general feasibility of chia cultivation in southwestern Germany was proven by Grimes et al. as varieties, either adapted or insensitive to photoperiod, have been made available [22]. On this basis, a suitable cropping system was developed for southwestern Germany [23]. Within the genus Salvia, two varieties, namely SALV66 (Salvia hispanica L.) and Golden Chia (Salvia columbariae Benth.), were characterized and examined. Both varieties were compared to G8 (Salvia hispanica L.) for their performance under the given climatic conditions since G8 has already been cultivated for several years and was therefore considered as the reference. As reported in more detail in the subsequent sections, the two Salvia hispanica L. varieties G8 and SALV66 differed strongly in (i) growth habitus, (ii) photosensitivity reaction, (iii) flower induction/harvest maturity, and (iv) seed yield compared to the Salvia columbariae Benth. variety Golden Chia. To our knowledge, this is the first time the agronomical and quality traits of SALV66 and Golden Chia have been reported.

Plant Development
Both Salvia hispanica L. varieties (G8, SALV66) accumulated almost twice as many growing degree days (451.85, 58 days after sowing (DAS)), flowering simultaneously, compared to Golden Chia (247.15, 34 DAS) in order to induce (Table 3). Radiation, average day length, and average mean temperature until flower induction amounted to 10,399 Wh/m 2 , 15.74 h, and 18.1 • C, respectively, for G8 and SALV66 and were found to be similar to Golden Chia at 10,563 Wh/m 2 , 15.84 h, and 17.7 • C, respectively (Table 3).
Previously conducted field trials by Grimes et al. [22,23] showed that cultivated Salvia hispanica L. varieties required 65 to 74 days in order to induce flowering, while the radiation and day length obtained were similar (~10,200 Wh/m 2 ; 15.8 h) to the present study. In addition, it could be shown that during those previously conducted field trials, cultivated Salvia hispanica L. varieties accumulated more growing degree days (463-545), and the average mean temperature until flower induction tended to be lower (16.8-18.0 • C) compared to the present study [22,23]. The results of the present study therefore verify the findings of Grimes et al. [22] and Baginsky et al. [6], which identified the threshold of 15.8 h of daylength for flower induction to be adequate for the two newly cultivated varieties (SALV66 and Golden Chia), being of higher importance for their development than the accumulated growing degree days. Both Salvia hispanica L. varieties, G8 and SALV66, reached harvest maturity 131 DAS, after having accumulated 1006.85 growing degree days (GDD), slightly more compared to previously conducted trials in southwestern Germany [22,23]. Golden Chia, on the other hand, reached harvest maturity after 93 DAS, while having accumulated 802.1 GDD. Average mean temperature from flower induction to harvest maturity was 17.4 • C (G8, SALV66) and 20.5 • C (Golden Chia). From the results obtained in this present study, it is reasonable to assume that temperature led to an accelerated plant development of G8 and SALV66 accompanied with an earlier beginning of flowering, as already stated by Grimes et al. [22]. The relatively short vegetation period of Golden Chia compared to the results of the trial conducted by Ayerza and Coats [34] in the Tucuman/Bolivian Forest implies that the dry and hot climatic conditions in southwestern Germany (2018) were more favorable as Golden Chia is, in general, more adapted to desert climates, preferring dry and sandy soils ( Table 3).
As Grimes et al. [22] discussed previously, sowing date is of great significance in regard to plant development, as it affects day length, GDD, temperature, and vegetation period. Sowing dates in later May limit the risk of frost damage to the emerging Salvia hispanica L. plants as the last spring frost often occurs at the end of April/beginning of May (ice saints), while simultaneously exposing plants to the risk of early frost damage occurring in October/November in case of late ripening under the given climatic conditions of southwestern Germany. Therefore, the relatively short vegetation period of Salvia columbariae Benth. compared to Salvia hispanica L. constitutes an advantage, as maturity is most certainly reached prior to partially unfavorable climatic conditions in late autumn. Furthermore, the short vegetation period of Salvia columbariae Benth. enables a greater range of subsequent crops as the necessary fieldwork (soil tillage, etc.) can be conducted earlier during the year. In general, it can be stated that the cultivated varieties (G8, SALV66, Golden Chia) completed their phenological development and were able to reach maturity, leading to the assumption that the thermal and photoperiod requirements were fulfilled [22,23,35,36].
In terms of germination, we could not detect any problems regarding Golden Chia compared to former studies conducted in California (greenhouse, laboratory, and field experiments), which indicated that the ambient field temperatures seemed to be adequate to alleviate the seed dormancy barrier for better germination to achieve a uniform number of plants per m 2 [37].

Plant Growth
The statistical analysis of the data indicated a significant difference between the Salvia hispanica L. and Salvia columbariae Benth. varieties concerning the obtained plant height, number of branches and inflorescences, and the length of the central axis main inflorescence at harvest time (p = 0.0005, p = 0.0002, p = 0.0145, p = 0.0020) ( Table 4).  (Figure 1). Means which share a common lowercase letter are not significant at α = 0.05.
The morphological differences shown in Table 4 can be attributed to the different growth habitus of the cultivated varieties as shown in Figure 1. To state only some observed differences, which are in accordance with literature, Salvia columbariae Benth., for example, grows 10 to 50 cm tall. It develops a tap rooted, rosetted, and unbranched or sporadically branched stem. The stems emerge singly or in several branches at the base and form terminal inflorescences, subtended by a few leaves; one or more spherical clusters, each containing one dozen to several dozen flowers reaching 8 to 20 mm in diameter, also appear [38,39]. Salvia hispanica L., on the other hand, reaches growth heights of up to 1.75 m, has shallow roots, and its stem is sparsely branched and forms terminal pseudo whorl with blue or white flowers which can be more than 20 cm long [1].
Golden Chia produced significantly shorter plants (37.48 cm) in comparison to the other two Salvia hispanica L. varieties, G8 (74.23 cm) and SALV66 (78.60 cm), due to the difference in growth habitus as already indicated (Table 4). With a mean of 2.08 branches per plant, Golden Chia produced significantly less branches compared to G8 (7.77) and SALV66 (8.27), which were in the range (7.9-9.3) of what Souza and Chavez [40] observed under greenhouse conditions in Brazil (Table 4). [22]. In regard to the number of branches developed, Sosa Baldivia and Ibarra [41] pointed to the fact that Salvia hispanica L. is subject to high plasticity as it is directly influenced by row spacing and seeding density, which has to be matched to the specific, given environmental conditions. G8 (7.98) produced a significantly lower number of inflorescences compared to Golden Chia (29.42) and SALV66 (15.97), decreasing from Golden Chia > SALV66 > G8. The number of inflorescences produced by Salvia hispanica L. varieties in 2016 cultivated at Ihinger Hof amounted to 11.8-13.3 [22]. As for the number of branches, the number of developed inflorescences is most probably also subject to the plasticity of chia [41]. G8 and SALV66 formed pseudo whorl inflorescences whereas Golden Chia developed compact, round inflorescences as shown in Figure 1.
Therefore, the length of the inflorescence of the central axis of SALV66 (12.26) and G8 (10.40) was measured whereas the width of the central axis inflorescence of Golden Chia (2.76) was recorded. The length of the central axis inflorescence of Salvia hispanica L. was in the range (7.92-14.82 cm) of what Hüller Goergen et al. [42] observed during a field trial conducted in Brazil but was considerably shorter compared to the results (16.8-19.2 cm) obtained at Ihinger Hof in 2016 by Grimes et al. [22]. The difference in length is most probably due to environmental influences and weather conditions, crop management, or a combination of those influences and specific species responses [43]. To our knowledge, the above-mentioned data of Golden Chia were presented for the first time; therefore, it would be of great interest to examine different environmental and crop management effects in regard to the cultivation of Golden Chia under the given conditions in southwestern Germany.

Yield Traits
Seed weight per single plant was not significantly influenced by variety (p > 0.05) ( Table 5). It ranged from 3.11 ± 0.6 (G8) to 3.63 ± 0.8 (Golden Chia) and 4.46 ± 0.7 (SALV66) respectively ( Figure 3). The same applied for the number of seeds produced per single plant, which varied from 2996 ± 540 (G8) to 3260 ± 695 (Golden Chia) and 3162 ± 525 (SALV66) (Figure 3). Table 5. Results of the statistical evaluation of seed weight (g) per single plant, number of seeds per single plant, yield (kg ha −1 ) ‡ , and thousand seed mass (TSM) (g) ‡ of three chia varieties (G8, Golden Chia, SALV66) cultivated at Ihinger Hof in 2018 (n = 3, α = 0.05).  Reported manually obtained seed yields and amounts of seeds produced per single plant were in the range of 6.3 to 6.4 g and 4051 to 4289 for Salvia hispanica L., respectively, being distinctly higher compared to the values obtained by the present study [22]. It is, in general, stated that obtained seed yields per single plant and seeds produced per single plant in different plant species could be influenced by various factors, such as plant density, genotypic differences, climatic conditions, and mineral nitrogen content [44][45][46][47]. Reported manually obtained seed yields and amounts of seeds produced per single plant were in the range of 6.3 to 6.4 g and 4051 to 4289 for Salvia hispanica L., respectively, being distinctly higher compared to the values obtained by the present study [22]. It is, in general, stated that obtained seed yields per single plant and seeds produced per single plant in different plant species could be influenced by various factors, such as plant density, genotypic differences, climatic conditions, and mineral nitrogen content [44][45][46][47]. For both SALV66 and Golden Chia, it can be assumed that single plant yield and seeds produced were most likely predominantly influenced by the environmental conditions of the respective growing season. The assumption which was made in regard to the high plasticity of Salvia hispanica L. enabling the compensation of lower plant densities with higher seed numbers per single plant could not be verified by the present study [22,41]. To our knowledge, this is the first time these specific yield traits of Salvia columbariae Benth. were published. Further field trials would have to be carried out in order to make well-founded statements about the possible specific effects influencing seed yields and numbers of seeds produced per single plant.

Number of Seeds
Golden Chia obtained a significantly lower mean yield (284.13 ± 74.44 kg ha −1 ) compared to SALV66 (643.99 ± 42.92 kg ha −1 ), whereas the obtained yield of G8 (491.91 ± 67.45 kg ha −1 ) did not differ significantly from Golden Chia and SALV66, as shown in Figure 4. As presented in Table 5 (41)(42)(43)(44)(45)(46)(47) at Ihinger Hof. In this regard, it must be stated that three days after sowing, a storm with downpour (35 mm) took place, followed by a dry period (Figure 2), leading to soil crusting and a partially irregular germination, which probably significantly influenced the obtained number of plants per m 2 and yields of the cultivated varieties [48]. The myxocarpic characteristic of the Salvia genus generally constitutes an advantage to seeds coping with restricted and irregular water availability during germination and early seedling development, as its mucilage may aid in supplying water to the seed [49,50]. However, the potential advantage of the mucilaginous chia seeds might have been overruled by the susceptibility of the seedlings, which could not penetrate the soil crust [51].
The myxocarpic characteristic of the Salvia genus generally constitutes an advantage to seeds coping with restricted and irregular water availability during germination and early seedling development, as its mucilage may aid in supplying water to the seed [49,50]. However, the potential advantage of the mucilaginous chia seeds might have been overruled by the susceptibility of the seedlings, which could not penetrate the soil crust [51]. In this present study, G8 obtained a noticeably lower seed yield compared to previously obtained seed yields of G8 at Ihinger Hof in 2016 and 2017 (752-1170 kg ha −1 ) cultivated at the same sowing rate, row spacing, and nitrogen fertilization rate. The relatively low seed yields obtained by both of the Salvia hispanica L. varieties in this study could be explained first and foremost due to unfavorable field conditions and the In this present study, G8 obtained a noticeably lower seed yield compared to previously obtained seed yields of G8 at Ihinger Hof in 2016 and 2017 (752-1170 kg ha −1 ) cultivated at the same sowing rate, row spacing, and nitrogen fertilization rate. The relatively low seed yields obtained by both of the Salvia hispanica L. varieties in this study could be explained first and foremost due to unfavorable field conditions and the accelerated plant development resulting from them, as addressed in Section 3.1 and shown in Table 3. In this regard, Mathan et al. stated that flower induction, which marks the transition from vegetative to generative plant growth, is, inter alia, decisive for optimal yield in crop plants [52]. The average mean temperature until flower induction was slightly higher (0.1-1.3 • C) compared to the previous studies conducted at Ihinger Hof in 2015 to 2017; furthermore, the average mean temperature from flower induction to harvest maturity was considerably higher (3-5.5 • C). Those high temperatures during the grain filling phase probably led to an increase in the photosynthesis rate above its optimum, leading to reduced metabolic activity during grain filling and maturation and, thus, to yield losses [6,16,22]. Alongside the relatively high temperatures, precipitation was partly substantially lower (between 38-193 mm) compared to the studies conducted at Ihinger Hof between 2015 and 2017, which certainly contributed to the comparably low seed yields obtained. Even though it is stated that Salvia hispanica L. is able to grow in dry conditions and is semi-tolerant to drought, this only applies to already established plant stocks [35,[53][54][55]. It is indicated that Salvia hispanica L. profits from precipitation and general water availability ranging from 300 to 1000 mm during the entire vegetation period, having a higher water demand during its vegetative phase and while establishing the reproductive organs, whereas drier conditions are required during the generative phase and seed maturity [56][57][58]. In this regard, Steduto et al. stated two possible negative effects on yield in general: (1) the inhibition of water stress on pollination and successful formation of the embryo, which could have led to a reduced number of set fruits (or grains) and, therefore, limited harvest index and seed yield, respectively, and (2) the under-filling and abortion of younger fruits resulting from a lack of photosynthetic assimilates available to be used for fruit developing [58].
In contrast to that, the seed yield of Golden Chia could, in comparison to the obtained seed yield reported by Ayerza and Coats (133 kg ha −1 ), be considered as quite high. The disparity between the seed yield (g) produced per single plant and the seed yield per kg ha −1 of Golden Chia and the significantly lower seed yield (kg ha −1 ) compared to the other Salvia hispanica L. varieties can be explained by the degree of domestication. Domestication of Salvia hispanica L. began, according to Ayerza and Coats, more than 2000 years ago, whereas Golden Chia remains highly undomesticated to date [1]. The high seed shattering level is a distinct attribute of the wild stage (lack of domestication), leading to enormous seed yield losses due to threshing [34,59]. Further breeding approaches towards reduced seed shattering should, therefore, be developed in order to be able to cultivate Golden Chia on an economically profitable basis. Based on our results, it can be stated that Golden Chia holds great potential, especially with regard to its photosensitive reaction, as flower induction and harvest maturity were reached significantly earlier, reducing the risk of frost damage due to late harvests. One general factor possibly leading to lower obtained seed yields is the choice of harvesting time and harvest intensity (manually/mechanically). It has been reported that sub-optimal combine operations and harsh weather conditions could lead to yield losses of up to 37% as chia is highly susceptible to seed shattering [60,61]. In this context, the non-uniform (top-down) maturation of Salvia hispanica L. should be mentioned, posing an issue in regard to optimal harvest time as the central inflorescence matures and ripens while inflorescences on side branches are still immature [23,62]. The same applies to Salvia columbariae Benth., as it could be observed that the flowers generally ripen from the inside out in addition to a delayed ripening of the side branches. In addition, the short vegetative period could, in this regard, negatively influence the potential seed yield of Salvia columbariae Benth.
Thousand seed mass as another yield determining trait seemed to be influenced exclusively by variety (p < 0.0001). As illustrated in Figure 4, the TSM of G8 was 1.17 ± 0.01 g, while the TSM of Golden Chia was 0.92 ± 0.01 g, and SALV66 obtained a TSM of 1.36 ± 0.01 g. TSM differed significantly between each of the varieties. The TSMs for G8 and SALV66 were in the range of current literature (1.1-1.4 g) [19,22,63]. In regard to obtained Salvia hispanica L. TSMs, it could be shown that a high number of plants per m 2 leads to interplant competition and therefore to a decrease in TSM [23]. Other studies related to wheat (Triticum aestivum L.) showed that plant density and nitrogen fertilization influenced TSM significantly [64][65][66]. Identifying factors which influence the TSM of the cultivated Salvia varieties could be advantageous in order to further improve agronomic traits of the cultivated varieties in general and especially under the given climatic conditions in southwestern Germany.
With regard to Golden Chia, it should be mentioned that there might be a need to develop a specific cultivation system which differs in row distance, seeding density, and nitrogen fertilizer rates, as well as a specific harvest process from that established for Salvia hispanica L. in order to exploit its full potential, leading to higher seed yields and TSM [23].

Protein, Mucilage and Oil content
Variety, numbers of plants per m 2 , and the interaction of variety and numbers of plants per m 2 had no significant effect on crude protein content (p > 0.05). Crude protein content was 24.78 ± 0.9% (G8), 22.14 ± 1.1% (Golden Chia), and 24.46 ± 1.0% (SALV66) ( Table 6). The results of this present study are in line with data obtained in previous field trials conducted in southwestern Germany and current literature [22,23]. The decrease in protein content due to high temperatures and due to an increase in altitude and elevation could not be verified, as obtained protein contents can be classified as high even though the average mean temperatures until flower induction (0.1-1.3 C) and harvest maturity (3-5.5 C) were higher compared to previously conducted studies at Ihinger Hof between 2015 and 2017 as mentioned in Section 3.2.2 [22,23,63,67]. Salvia hispanica L. and Salvia columbariae Benth. species do not only exceed the protein content quantitatively but also qualitatively, compared to other starch crops such as wheat, corn, rice, oats, and barley. It was shown that chia can be incorporated into human diets as a plant-based protein enrichment in order to produce more balanced protein food sources [1,16].
The statistical analysis indicated significant differences in the mucilage content between the examined varieties (p ≤ 0.001). Additionally, it seemed that the numbers of plants per m 2 did influence the mucilage content (p ≤ 0.01). As presented, the total mucilage content of the seeds was 10.78 ± 0.3% (G8), 20.66 ± 0.3% (Golden Chia), and 10.35 ± 0.2% (SALV66) ( Table 6). The mucilage content of Golden Chia was displayed and statistically evaluated although the values are not representative since it could be demonstrated that the applied mucilage extraction method was only suitable for G8 and SALV66, being in line with the values obtained by current literature [22,23,68]. Seeds of Golden Chia crumbled during the extraction process, which led to a major surface enlargement. The statistically significant difference which was found regarding the mucilage content between the Salvia hispanica L. varieties G8 and SALV66 and the Salvia columbariae Benth. variety Golden Chia can be explained by this fact. The mucilage extraction method for Golden Chia could be adapted by decreasing the temperature and hydration time. According to Tavares et al., first attempts examining the applicability of cold extraction method for Salvia hispanica L. seemed to be promising [69]. Further analysis should be conducted in order to make well-founded statements about the possible content of mucilage extracted from Golden Chia. In regard to the mucilage content, chia could be used inter alia as a vegan thickener, foam stabilizer, suspending agent, emulsifier, adhesive, or binder due to its water-holding capacity and viscosity [5,70]. Another field of application that is of particular interest in terms of environmental protection is the constant and recurring theme around packaging waste. In the packaging industry, mucilage obtained from chia seeds may be used for functional coating and edible films, representing a promising alternative to synthetic packaging [5,71].
The crude oil content of the seeds was not significantly affected by variety, plants per m 2 , or the interaction between variety and plants per m 2 (p > 0.05, Table 6). Crude oil content was within the range stated in recent literature, varying from 30.52 ± 1.0% (G8) to 28.00 ± 0.6% (Golden Chia) and 31.73 ± 1.1% (SALV66), though still being lower compared to the crude oil contents obtained by G8 in previous trials conducted at Ihinger Hof in southwestern Germany (30.9-33.8%) [6,12,22,23,34,72]. In this regard, it must be pointed out that the non-uniform maturation of Salvia hispanica L. (top-down) and Salvia Columbariae Benth. (inside-out) makes it difficult to determine the optimal harvest time and, hence, may have contributed to the low oil content due to a possible high number of immature seeds. Generally, the variation in chia seed oil yields depend on factors such as climatic conditions [67,73].
As summarized by Hrnčič et al., chia oil obtained from Salvia hispanica L. is one of the most valuable oils on the market today and is used in the pharmaceutical and food industries [5]. In this regard, Hrnčič et al. draw particular attention to nano-emulsion-based delivery systems based on chia seed oil from Salvia hispanica L. for prospective applications to encapsulate lipophilic bioactive components in food, cosmetics, and pharmaceuticals [5]. In order to meet the given EU requirements of the Novel Food guideline, crude oil contents of Salvia hispanica L. must range from 30 to 34%, which is also in accordance to current literature [22,63,74,75]. As the crude oil content of Golden Chia obtained by this present study was below 30%, it would not meet those requirements yet. Subsequent trials would be necessary in this respect in order to validate or falsify the obtained results, as seeds of Golden Chia would also constitute a promising raw material for the production of edible, high-quality oil and for the previously mentioned applications [12].

Fatty Acids
Results of the ANOVA showed that the fatty acid composition of the investigated Salvia varieties was significantly affected by variety (p ≤ 0.01, Table 7). As shown in Table 2, contents of all individual fatty acids, except for vaccenic acid, and respective ratios of Golden Chia differed significantly from the Salvia hispanica L. varieties G8 and SALV66. All samples were characterized by a high content of ω-3 and ω-6 polyunsaturated fatty acids, namely α-linolenic acid as the main fatty acid (60.4 to 65.9 g/100 g) and linoleic acid (15.1 to 20.4 g/100 g), being in line with values obtained by previously conducted field trials in southwestern Germany [12,22,23,[75][76][77][78][79]. We also verified the findings of Matthäus and Özcan showing that within the genus Salvia, the amount of α-linolenic acid is highest in Salvia columbariae Benth., making it highly interesting from a nutritional point of view [12].
The high seed lipid concentration (~30%), especially of omega-3 and omega-6 fatty acids, is due to proteins connected to the production and storage of plant lipids [80]. In accordance with Ayerza and Coats [34], the study could demonstrate that the amount of α-linolenic acid obtained by Golden Chia was higher compared to that of the Salvia hispanica L. varieties G8 and SALV66, whereas it was the opposite for linoleic acid, as presented in Table 7. The amounts of monounsaturated fatty acids obtained from G8 and SALV66, namely oleic (6.9 to 9.1 g/100 g) and vaccenic (0.8 to 0.9 g/100 g), were in line with current literature [7,22,23,81,82]. The same applied to the amounts of saturated fatty acids, which varied between 5.7 and 7.7 g/100 g for palmitic acid and between 3.2 and 3.8 g/100 g for stearic acid. The only published article to date which examines the lipid fraction of seeds from Salvia columbariae Benth. grown in Arizona obtained slightly lower mean contents of vaccenic (1.0 g/100 g), palmitic (4.8 g/100 g), and stearic (2.1 g/100 g) fatty acids and a slightly higher amount of oleic acid (9.8 g/100 g) [12].
In the context of fatty acid composition, it should be pointed out that a high amount of polyunsaturated fatty acids determines the low oxidative stability of chia oil, which represents a solvable challenge as respective production processes are available [83].
According to current literature, chia fatty acid compositions are highly dependent on environmental conditions and, therefore, are affected by location [67,73,84,85]. In addition, it was shown that the relationship between elevation, fatty acid composition, and oil saturation for chia is most likely related to a temperature-elevation interaction since elevation is strongly negatively related to temperature, which, as such, affects oil content variation during seed development [67,73,84,85].This applies to individual fatty acids, which are influenced by irrigation levels, climate conditions including temperature, elevation, and proteins [67,73,80,86,87].
In this context, various approaches are discussed in current literature. Herman et al. [86], for example, mentioned that water stress contributes to an increased production of α-linolenic acid, whereas Amato et al. showed inter alia that mineral fertilization increased free acidity in seeds while it reduced oxidative stability.
The importance of qualitative analysis regarding fatty acids becomes clear when looking at the changes in Western diets during the last decades. Total fat intake in general has continuously decreased, shifting towards a highly carbohydrate-focused diet [88]. The enormous decrease in omega-3 fatty acid intake in addition to the tremendous increase in omega-6 fatty acid intake has resulted in an increased omega-6 to omega-3 ratio up to 20:1 or even higher as Simopoulos stated [88]. This shift in the composition of fatty acids has led to a rise in overweight and obesity as well as pathogenesis of other diseases such as cardiovascular illnesses, cancer, and inflammatory and autoimmune diseases, whereas increased levels of omega-3 PUFA (a low omega-6 to omega-3 ratio) exert suppressive effects [88,89]. Therefore, a balanced omega-6 to omega-3 ratio must be emphasized in regard to general health and in the prevention and management of obesity [88].
The results of the present study indicate that due to the composition of fatty acids and vitamin E active compounds of Salvia hispanica L. and Salvia columbariae Benth., an increased intake of chia seeds and seed oil could contribute to a more balanced omega-6 to omega-3 diet and high-value food products, possibly leading to an overall improved nutrition [12,22]. Table 7. Fatty acid composition (% of total fatty acid) of three chia varieties (G8, Golden Chia, SALV66) cultivated at IHO in 2018 (n = 3, α = 0.05, mean ± standard error).

Conclusions
It can be concluded that within the given biodiversity of Salvia species, there are existing Salvia hispanica L. and Salvia columbariae Benth. varieties adapted to day lengths greater than 12 h. Besides intensive breeding approaches, the given range in biodiversity enables the cultivation of chia under the existing day length conditions of southwestern Germany. With seed yields obtained being more or less in line with those of their countries of origin, the selected chia varieties SALV66 and Golden Chia represent very promising raw materials from a nutritional point of view. Chia cultivation in Germany would fulfill the steadily increasing demand for healthy, environmentally sustainable, and regionally produced food products while simultaneously representing a profitable economic source of income for local farmers, potentially leading to an advantageous local, economically and environmentally profitable way of broadening the food chain and food base under the given conditions in southwestern Germany. However, seed shattering of Golden Chia is a significant drawback in regard to its commercial production. Further breeding efforts and scientific studies on agronomic management practices are therefore necessary.

Data Availability Statement:
The data presented in this study are available on request from the corresponding author. The data are not publicly available due to project restrictions.