ORIGINAL RESEARCH article
Studying the Impact of Different Field Environmental Conditions on Seed Quality of Quinoa: The Case of Three Different Years Changing Seed Nutritional Traits in Southern Europe
- 1Department of Biology, Universidad Autónoma de Madrid, Madrid, Spain
- 2Castile-Leon Agriculture Technology Institute (ITACyL), Valladolid, Spain
- 3Agrarian Research Institute “La Orden-Valdesequera” of Extremadura (CICYTEX), Badajoz, Spain
- 4Department of Agroforestry Sciences, Universidad de Huelva, Huelva, Spain
- 5Department of Food & Health, Instituto de la Grasa, CSIC, Seville, Spain
- 6Cereal Group, Institute of Agrochemistry and Food Technology (IATA-CSIC), Valencia, Spain
- 7Department of Abiotic Stress and Plant Pathology, Centro de Edafología y Biología Aplicada del Segura (CSIC), Murcia, Spain
- 8Departamento de Producción y Caracterización de Nuevos Alimentos, Instituto de Investigación enCiencias de la Alimentación (CIAL) (CSIC-UAM), Madrid, Spain
- 9Sección Departamental de Ciencias de la Alimentación, Facultad de Ciencias, Universidad Autónoma de Madrid, Madrid, Spain
Chenopodium quinoa Willd (quinoa) has acquired an increased agronomical and nutritional relevance due to the capacity of adaptation to different environments and the exceptional nutritional properties of their seeds. These include high mineral and protein contents, a balanced amino acid composition, an elevated antioxidant capacity related to the high phenol content, and the absence of gluten. Although it is known that these properties can be determined by the environment, limited efforts have been made to determine the exact changes occurring at a nutritional level under changing environmental conditions in this crop. To shed light on this, this study aimed at characterizing variations in nutritional-related parameters associated with the year of cultivation and different genotypes. Various nutritional and physiological traits were analyzed in seeds of different quinoa cultivars grown in the field during three consecutive years. We found differences among cultivars for most of the nutritional parameters analyzed. It was observed that the year of cultivation was a determinant factor in every parameter studied, being 2018 the year with lower yields, germination rates, and antioxidant capacity, but higher seed weights and seed protein contents. Overall, this work will greatly contribute to increase our knowledge of the impact of the environment and genotype on the nutritional properties of quinoa seeds, especially in areas that share climatic conditions to Southern Europe.
Chenopodium quinoa Willd, commonly known as quinoa, belongs to the Amaranthaceae family native to the Andean region (Wilson, 1990; Alandia et al., 2020). Its natural distribution is extended from northern Colombia to the southern region of Chile, and it can be cultivated in a wide range of altitudes, from sea level up to 4,000 m above sea level (Zurita-Silva et al., 2014). In the last decades, the cultivation of this crop has expanded worldwide, although the main producers in the world are still Bolivia and Peru (Bazile et al., 2016). One of the reasons for the increased interest in cultivating quinoa is the capacity of adaptation and its resilience to extreme conditions (Jacobsen et al., 2003). Quinoa can tolerate drought, high soil salinity, frost, and low temperatures (Jacobsen et al., 2005, 2012; Pulvento et al., 2010; Adolf et al., 2012), which makes it an ideal crop to be exploited and introduced in marginal environments (Choukr-Allah et al., 2016). On the other hand, the remarkable nutritional traits of quinoa seeds are key for its recent risen popularity for human consumption (Abugoch James, 2009). Quinoa is a very valued food for its high protein content, which is higher than that of cereal crops like barley, wheat, maize, and rice (Koziol, 1992). This protein is also of higher quality since it contains all amino acids essential for human consumption (including lysine, methionine, and cysteine). Furthermore, their proportions are well balanced and close to the amino acidic profile recommended by the Food and Agriculture Organization (FAO) (Filho et al., 2017). Quinoa seeds also contain fiber, vitamins, and minerals like calcium, zinc, magnesium, iron, potassium, phosphorus, manganese, copper, and sodium, and stand out for their unsaturated fatty acid contents and their large antioxidant capacity (Abugoch James, 2009; Paśko et al., 2009). The high content of vitamins (A, B, and E) and polyphenols, like phenolic acids or flavonoids, contribute to this high antioxidant capacity, and make of quinoa seeds an excellent example of “functional food,” since these antioxidants may prevent cancer, cardiovascular and other chronic diseases (Paśko et al., 2009; Tang and Tsao, 2017). Another interesting aspect of quinoa seed composition is their lack of gluten, which makes this food suitable for people with coeliac disease (Peñas et al., 2014). However, it should be noted that quinoa seeds also present important concentrations of antinutrient components such as saponins, which cause the characteristic bitter taste in non-desaponificated seeds (Ruales and Nair, 1993; Mastebroek et al., 2000). Nevertheless, saponins are also of current popularity as phytochemicals with bioactive properties for health (Marrelli et al., 2016; Singh et al., 2017; Navarro del Hierro et al., 2018; El Hazzam et al., 2020). Given quinoa’s nutritional quality and its ability to grow in a wide range of climatological conditions, the FAO considers that this crop has the potential for playing an important role in worldwide food security (Ruiz et al., 2014; FAO and CIRAD, 2015).
Given the high demand for quinoa, numerous breeding programs have been developed aiming to obtain new varieties better adapted to new agricultural areas. These breeding programs have been mainly focused on the generation of varieties less sensitive to the photoperiod, more resistant to downy mildew, with low saponin contents (sweet varieties), and with increased yields, all these aiming to satisfy the high demand for quinoa (Zurita-Silva et al., 2014). Grain size has also been used as a selection criterion in response to commercial appeal (Zurita-Silva et al., 2014). However, little attention has been paid to the improvement of nutritional properties in new varieties, even though food quality is as important as food quantity for food security (Tester and Langridge, 2010).
Nutritional quality in quinoa seeds is variable (Craine and Murphy, 2020), and this variability results from the interaction of genetic and environmental factors (Wimalasekera, 2015). Several studies have reported effects of the agrological conditions (Gonzalez et al., 2012; Prado et al., 2014; Präger et al., 2018; Reguera et al., 2018) and environmental factors such as water availability, soil salinity, climatic conditions, and temperature (Pulvento et al., 2012; Miranda et al., 2013; Bascuñán-Godoy et al., 2015; Aloisi et al., 2016; Lesjak and Calderini, 2017; Curti et al., 2018) on nutritional traits of quinoa seeds. However, there is still limited knowledge about the mechanisms that trigger these changes in nutrient content and whether the seed nutritional traits are stable depending on the genotype, environment, and GXE interaction.
Thus, the aim of this study was to identify and evaluate changes in the nutritional properties of quinoa seeds linked to the genotype and those related to changes in the environmental conditions. For this purpose, six bred cultivars were grown during three consecutive years (2017, 2018, and 2019) and different physiological and nutritional parameters were determined. The results presented here suggest that environmental factors heavily influenced most of the nutritional parameters analyzed.
Materials and Methods
Plant Material, Experimental Design, and Location
Quinoa cultivars Regalona (registered variety of BAER, Chile), Puno, Titicaca, and Vikinga (Quinoa Quality, Denmark), Q3 and Q5 (International Center for Biosaline Agriculture (ICBA), Dubai, United Arab Emirates) were cultivated in the Experimental Station of Zamadueñas, which belongs to the Instituto Tecnológico Agrario de Castilla y León (ITACyL), in Valladolid (Spain) (41°42′N and 4°42′W, 690 m.a.n.l) under irrigation in three consecutive growing seasons (2017, 2018, and 2019).
Climatological data (as a monthly average of daily data), including total precipitation and total irrigation, is presented in Supplementary Table 1. The data was obtained from a local climatological station located at the field experimental station of Zamadueñas. In 2017 and 2019, quinoa plants were sown in the spring (April and May, respectively). In 2018, due to heavy precipitations, the sowing date was delayed until June 18th.
The soil, containing 42% silt, 17% clay, and 41% sand, was a clay-silty-loam type presenting a pH ranging from 8.2 to 8.55. It contained 2.03–1.19% organic matter, and 0.053–0.083 dS.m–1 of electrical conductivity (EC) of the saturated paste. Phosphorous content (as ppm of P2O5) ranged from 48 to105, potassium content (as ppm of K2O) ranged from 0.33 to 0.472 and total nitrogen (%) ranged from 0.067 to 0.075 (Supplementary Table 2). The basal fertilization in the plots consisted of 300 kg ha–1 of 8-15-15 (NPK) in addition to top dressing of 500 kg ha–1 of NSA (26%) at six-leaves stage. The trials were kept free of weeds and pests.
The experimental design consisted of randomized blocks with 3 replications. Each block was 8.0 m long and 3.0 m broad (24 m2) with 6 rows separated 0.50 m. Sowing was carried out mechanically using a sowing density of 10 kg/ha between 1 cm and 2 cm depth. Harvesting was also carried out mechanically.
Seed Weight and Seed Area
Seeds were manually counted and weighed in an analytical balance. Seed area was analyzed using the open-source software ImageJ1. Images were taken using an Olympus SZ61 stereomicroscope (Olympus Corporation, Shinjuku, Tokyo, Japan) and processed with the AnalySIS GetIT image software (analysis getIT 5.1, Olympus Corporation).
Color parameters were determined as described by Guiotto et al. (2020).
Seed Germination Rate
Quinoa seeds were sterilized first in ethanol 70% (2 min), followed by a wash in bleach 50% with a droplet of Tween-20 (2 min) and then rinsed several times in distilled water (H20). Sterilized seeds were sown on a double layer of paper filter wet with distilled water on Petri dishes and then transferred to a growth chamber under darkness and a controlled temperature of 23°C. Germinated seeds (considered as germinated when the radicle protrusion was longer than 2 mm) were counted daily for the first week after sowing.
Seed viability was performed using the tetrazolium method (2,3,5-triphenyl-2H-tetrazolium chloride). First, seeds were imbibed in distilled water at 30°C for an hour in order to facilitate longitudinal and superficial cuts of the embryo and to ensure a homogeneous dying of the seed tissues. After cutting, seeds were submerged in 1% tetrazolium chloride at 30°C for 2 h. Seeds with more than 50% staining in the embryonic tissue were considered viable.
Saponins were quantified following the protocol described by Navarro del Hierro et al. (2020) based on a previous extraction with methanol assisted by ultrasound and a subsequent analysis by HPLC-DAD.
Amino Acid Quantification
Amino acid analysis was performed following the protocol described by Villanueva et al. (1999).
The mineral content was analyzed following the official methods of analysis of the Spanish Ministry of Agriculture (MAPA, 1995). The phosphorus content was determined using a spectrophotometer UV-VIS (Hitachi U-2810) (yellow coloration, 430 nm). Potassium was determined using flame atomic emission spectroscopy. Calcium, magnesium, sodium, copper, manganese, zinc, and iron content were assessed using flame atomic absorption spectroscopy (AAS) (SpectrAA 110, Agilent) after mineralizing the samples with H20 and HCl (35%).
Evaluation of Antioxidant Capacity: Ferric Reducing Antioxidant Power (FRAP) Assay, Total Phenolic Content (TPC), and Total Flavonoid Content (TFC)
Total extracts were obtained from 100 mg of ground seeds, that were homogenized in 1 ml of an extraction buffer consisting of methanol (50%), acetic acid (1%), and distilled water (49%). The samples were then vortexed for 2 min and centrifuged for 15 min at 10,000 rpm. The supernatants were stored at −20°C until their use in the FRAP, phenol, and flavonoid assays.
Ferric Reducing Antioxidant Power (FRAP) Assay
The antioxidant capacity of seed samples was determined following the procedure described by Benzie and Strain (1996). The FRAP reagent consisted of a mix of 300 mM acetate buffer (pH 3.6), with 10 mM TPTZ in 40 mM HCl and 20 mM FeCl3⋅6H2O at a ratio of 10:1:1 (v/v/v). Twenty μl of sample extract and 180 μl of FRAP reagent were added into a 96-well microplate and, after 4 min, absorbance was read at 593 nm using a microplate reader (Lector Multi-ModalSynergy HTX, BioTek Instruments, Inc., United States). The antioxidant capacity was calculated from a calibration curve obtained with iron (II) sulfate (FeSO4). FRAP value was expressed as μmol of Fe2+/g of seed.
Total Phenol Content (TPC)
The content of polyphenols was measured following the protocol described by Tang et al. (2015). Briefly, the mixture of 50 μl of sample extract, 50 μl of the Folin-Ciocalteu reagent 10%, and 100 μl of sodium carbonate 13% was incubated for 60 min. Absorbance was read at 750 nm using a microplate reader (Lector Multi-ModalSynergy HTX, BioTek Instruments, Inc., United States). The TPC was expressed as mg of gallic acid equivalents per g of quinoa seed (mg GAE/g).
Total Flavonoid Content (TFC)
Flavonoid content was determined following the procedure described by Valenzuela-Bustamante (2015). Briefly, 30 μl of sample extract, 10 μl of aluminum chloride (AlCl3) 10%, 10 μl of sodium acetate (NaC2H3O2) 1M, and 250 μl of dH2O were mixed and incubated for 30 min. The absorbance was read at 415 nm using a microplate reader (Lector Multi-ModalSynergy HTX, BioTek Instruments, Inc., United States). The results were expressed in mg of quercetin equivalents per g of quinoa seed (mg QE/g).
To analyze the Genotype x Year interaction, two-way ANOVA was performed. Normality and equality of variances of the data were tested through a Kolmogorov-Smirnov’s test and a Levene’s, respectively. For variables where normality and equal variances could be assumed, a One-way ANOVA test was performed, followed by a Tukey post-hoc test, to perform multiple comparisons at a probability level of 5% (p < 0.05). A Krustal-Wallis test by ranks was performed when data did not present a normal distribution and a Welch’s ANOVA test followed by a Games-Howell post-hoc test was performed when variances were not equal, both at a probability level of 5% (p < 0.05).
Correlations amongst variables were evaluated with a Pearson’s correlation coefficient test and simple linear regressions were performed to analyze the relation between yield and qualitative seed variables. A sequential path analysis was performed to evaluate the specific contribution of nutritional seed traits to germination rate. This analysis allows ordering different variables as predictors of seed germination rate of first, second, or third-order (Mohammadi et al., 2003). For this purpose, a stepwise multiple linear regression procedure was used where variables that showed weak contribution (p > 0.05) to the dependent variable (germination rate) or high multicollinearity, were automatically dropped from the model. The variables entered into the model were considered as first-order predictors and the procedure was repeated using these variables as the response variable to identify traits that function as second-order predictors of germination rate. Tolerance and variance inflation factor (VIF) were used to measure the level of multicollinearity for each predictor trait, considering tolerance lower than 0.1 or VIF values higher than 10 as high levels of collinearity. Tolerance (1- R2i, where R2i is the coefficient of determination for the prediction of variable i by the predictor variables) is the amount of variance of the selected independent variable not explained by other independent variables. VIF (1/Tolerance) indicates the extent of effects of other independent variables on the variability of the selected independent variable. Principal component analysis was performed for viability and germination rates, yield, 1000 seeds’ weight, seed area, all color parameters, protein content, amino acid contents, FRAP value, phenols and flavonoids contents, and mineral contents. The SPSS Statistics 23.0 (SPSS Inc.) package was used for the statistical analyses.
Plant Performance During the Three Consecutive Years
Field trials were performed in three consecutive growing seasons (2017, 2018, and 2019) using 6 different cultivars: Titicaca, Vikinga, Regalona, Puno, Q3, and Q5. In seasons 2017 and 2019, quinoa seeds were sown in April and May, respectively. However, in 2018, due to heavy rainfalls during the spring, sowing was delayed till June, which explains why higher minimum temperatures (accompanied with higher precipitations) were registered during the first month of cultivation (Supplementary Table 1 and Figure 1). Quinoa life cycle spans about 5 months, therefore, plants that were sown in 2017 and 2019 were harvested in the summer (August 31st and September 20th, respectively), facing higher maximum temperatures during their last months of growth, which include grain maturation (Supplementary Table 1). On the contrary, plants cultivated in 2018 were exposed to lower temperatures (5–10°C) in the last month before harvesting on October 25th (Supplementary Table 1). Furthermore, differences appeared among varieties when evaluating their lifespan (Supplementary Figure 1). Thus, while Puno cv., Titicaca cv., Q5 cv., and Regalona cv. showed similar phenological stages, Q3 cv. and Vikinga cv. presented longer life cycles.
Figure 1. Environmental conditions (A) Total water supply and mean temperature for each month of development. The 2017 data is shown in black bars and lines, 2018 data in light gray bars and lines and 2019 data in dark gray bars and lines. (B) Mean temperatures in 2017 (black line), 2018 (light gray line), and 2019 (dark gray line) and day-length (dashed black line) in each month.
Precipitations along the different growing seasons were very variable. In 2019, precipitations were low the first two months of cultivation (6.0–7.0 mm) increasing during the last three months. In 2017, the range of precipitations among months was wider, with rainfalls concentrated during the second and fourth months of cultivation (May and July). The second growing season, 2018, presented more extreme conditions with very high precipitations during the first month of cultivation (June, 70.6 mm) and no rainfall in the third month (August) (Supplementary Table 1 and Figure 1). Nonetheless, the irrigated conditions minimized the differences caused by the low precipitations in terms of water supply, being 2019 the year that showed a reduced water supply especially during the first two months of growth.
Total seed yield ranged from 0.70 t/ha to 3.25 t/ha, being Vikinga cv. 2019 an outlier, with a seed yield of 0,23 t/ha (Supplementary Figure 1). There were important yield differences among years: 2017 was the growing season that presented the highest yields (from 2.2 t/ha the Q3 cultivar to 3.25 t/ha Puno cv.). In contrast, 2018 showed the lowest yields except for Puno cv. and Q3 cv., which produced almost three times as much seed as the other cultivars.
Seed area and seed weight showed an effect related to the cultivar, year of sowing, and cultivar x year interaction (p < 0.05) (Figure 2 and Supplementary Table 3). One thousand seeds’ weight ranged from 1.7 to 3.4 g and was significantly higher in seeds harvested in 2018 and 2017 compared to 2019 seeds (Figure 2A). Seeds from the 2018 harvest were also larger compared to the seeds from other years, while seeds from 2019 presented the smallest areas (Figure 2B). When comparing cultivars, Titicaca cv. and Q5 cv. seeds were both heavier and larger while Puno cv. seeds were the smallest and lightest. Interestingly, Vikinga cv. seeds showed lower weights than those from the Q3 cultivar, but larger areas. Q5 cv. 2017 seeds presented the highest seed weight, followed by Titicaca cv. 2018. Both also showed the largest areas, being Titicaca cv. 2018 seeds twice as wide as the lightest and smallest seeds harvested from Puno cv. in 2019 (4.4 and 2.2 mm2, respectively). No correlation between these two parameters was found (Supplementary Figure 4), which indicates that seeds may differ in density and shapes depending on the cultivar and the environmental conditions.
Figure 2. Seed weight and area. Black bars represent 2017 values, light gray bars show 2018 values, and 2019 values are represented by dark gray bars. (A) Thousand seed weight An ANOVA test followed by a post-hoc test Tukey was performed. (B) Seed area. A Welch’s ANOVA test followed by a Games-Howell post-hoc test was performed. Bars that do not share the same letters show statistically significant differences.
Germination Rates and Seed Viability
To evaluate the germination capacity of the seeds, germination rates were determined for all cultivars harvested in the three consecutive years (Figure 3). The year of cultivation, the cultivar, and the interaction of these two factors had a significant influence on the germination rates (p < 0.05) (Supplementary Table 3). The results showed that 2018 was the year in which seeds showed lower germination rates, being Regalona the cultivar with the highest rate and Q3 the only variety that did not germinate. Most of the varieties showed similar germination trends in 2017 and 2019, except for Vikinga cv. and Titicaca cv. In 2019, these two cultivars showed similar rates to those obtained for the 2018 growing season.
Figure 3. Germination rates and seed viability. (A) Germination rates. Evolving germination rate is shown over the 7 first days after sowing seeds harvested in three consecutive years. (B) Seeds viability was evaluated using the tetrazolium method. Black bars, light gray bars, and dark gray bars indicate 2017 seeds, 2018 seeds, and 2019 seeds, respectively. A Welch’s ANOVA test followed by a Games-Howell post-hoc test was used to compare the data of germination rate at seven days after sowing, and a Krustal-Wallis test by ranks was performed for viability rate. Bars that do not share letters are significantly different.
Seed viability also presented high influence of year, cultivar, and the year x cultivar interaction, and was well correlated with the germination data (Supplementary Figure 4). Thus, 2018 seeds showed a steeped decrease in seed viability and Vikinga cv. seeds’ viability from the 2019 harvest was hardly measurable (Figure 3B).
Saponins are secondary metabolites whose content may vary when changing the environmental conditions (Szakiel et al., 2011). To analyze the effects of environmental conditions on saponin content in quinoa seeds, we studied the saponin content of two different cultivars: Titicaca cv., considered a bitter cultivar due to the higher concentration of saponin in its seeds, and Vikinga cv., a sweet cultivar with lower saponin content. Indeed, significant differences were found between the two cultivars, Titicaca cv. seeds containing more saponins than Vikinga cv. seeds (an average of 1.92 g/100 g and 1.26 g/100 g, respectively). Furthermore, there was an important increase in saponins in 2019, especially in Vikinga cv. seeds (2.04 g/100 g), that increased the saponin levels till reaching Titicaca cv. levels (Figure 4). The statistical analysis showed that both the year and the cultivar, played a significant role in determining saponin content (p < 0.05), but not the interaction between the two factors (p = 0.075) (Supplementary Table 3).
Figure 4. Saponin content in quinoa seeds. Saponin content of Titicaca cv. and Vikinga cv. seeds sown and harvested in different years is presented. Different lower-case letters indicate significant differences between cultivars within a year and upper-case letter show differences between years within a cultivar. T-student comparisons were performed between pairs of samples.
Protein Content and Amino Acids
Total protein content ranged from 13.8 to 19.1% of seed weight (Figure 5). The year, cultivar, and year x cultivar interaction factored in the determination of the protein content (p < 0.05) (Supplementary Table 3). On average, protein content was significantly higher in Vikinga cv. seeds compared to Q5 cv., Regalona cv., or Puno cv. Seeds harvested in 2018 showed higher protein content compared to the other two years of cultivation (Figure 5). Besides, differences were found within the same variety among growing seasons (e.g., Titicaca cv., Puno cv., and Q5 cv. harvested in 2017 showed lower content than the harvesting of 2018) and within the same year among varieties (e.g., in 2017, Vikinga cv. protein content was significantly higher than the content found in Puno cv.).
Figure 5. Protein contents. Protein content presented as percentage of protein per seed dry weight. Black bars represent 2017 values, light gray bars show 2018 values, and 2019 values are represented by dark gray bars. Bars that do not share the same letters show statistically significant differences, following the ANOVA test and post-hoc test Tukey.
Regarding the amino acid profile, in all samples, glutamic acid was the most abundant amino acid (21.6 mg/g of seed, on average), followed by arginine and aspartic acid, while tryptophan showed the lowest concentrations (0.7 mg/g of seed) followed by cysteine and methionine (Supplementary Figure 2). The year of cultivation was a determining factor for all amino acid contents (p < 0.05), but the cultivar was not a significant factor in the case of threonine, proline, tyrosine, methionine, tryptophan, phenylalanine, and isoleucine, and the interaction between these factors only influenced the aspartic acid, cysteine, and arginine contents. Generally (except for methionine, cysteine, and tryptophan), 2018 yielded higher amino acid concentrations and 2019 the lowest. When comparing cultivars, Vikinga cv. and Titicaca cv. showed higher concentrations; on the contrary, lower amino acid concentrations were found in Puno cv., Q3 cv., and Q5 cv. However, cysteine and tryptophan concentrations were higher in 2017 seeds compared to 2018 seeds, and cysteine and methionine contents were higher in Puno cv. and Q3 cv. seeds compared to Vikinga cv. or Titicaca cv. (Table 1). In the case of proline, a higher content was found in Titicaca cv., Vikinga cv., Regalona cv., and Puno cv. seeds harvested in 2017 (6.2–8.2 mg/g of seed), and lower in 2019 and in Q3 cv. and Q5 cv. seeds harvested in 2017 (2.5–4.0 mg/g of seed) (Table 1).
The total content of phosphorous (P), potassium (K), calcium (Ca), magnesium (Mg), sodium (Na), copper (Cu), iron (Fe), manganese (Mn), and zinc (Zn) in quinoa seeds were determined in order to analyze the effect of genotype and environment on the mineral content (Table 2).
Overall, it was observed that mineral content was greatly influenced by the year of cultivation and by the cultivar (p < 0.05), except for Na (p = 0.419 and p = 0.063, respectively), and by the interaction between the two factors (except for Mg, Fe, and Zn content) (Supplementary Table 3). In the case of P, Cu, Mn, and Zn (Table 2), contents were higher in seeds harvested in 2018, with a difference especially remarkable in Zn. On the other hand, this trend was inverted in the case of K, and the contents of Ca, Mg, and Fe were significantly lower in 2017 seeds (Table 2). For P, Mg, Fe, Mn, and Zn contents, Puno cv. seeds showed the lowest levels while Vikinga cv. showed the highest contents, while for K, the content was higher in Titicaca cv. and Regalona cv.
We evaluated the antioxidant capacity by performing the FRAP assay, together with the quantification of total polyphenols (TPC) and flavonoids (TFC) contents (Figure 6). These three variables showed strong correlations among them (p = 0,000, Supplementary Figure 4) and a significant influence of the year of cultivation, the cultivar (p < 0.05), and their interaction in the case of FRAP value and TFC.
Figure 6. Antioxidant capacity of quinoa seeds. Black bars represent 2017 values, light gray bars show 2018 values, and 2019 values are represented by dark gray bars. (A) Antioxidant power of quinoa seeds was measured using the ferric reducing antioxidant power (FRAP) assay and is expressed as μmol of Fe2+per gram of seed. Statistical differences were analyzed through a Welch’s ANOVA test followed by a Games-Howell post-hoc test. (B) Total polyphenol content (TPC) is expressed as milligrams of gallic acid equivalents (GAE) per gram of seeds. The statistical analysis performed was a One-way ANOVA test followed by a post-hoc Tukey test. (C) Total flavonoid content (TFC) is expressed as milligrams of quercetin equivalents (QE) per gram of seeds. A Krustal-Wallis test by ranks was performed for multiple comparisons. Bars that do not share the same letters show statistically significant differences.
Significant differences were found when comparing cultivars. Titicaca cv. showed the highest antioxidant capacity, TPC, and TFC, followed by Q5 cv., while Regalona cv. presented the lowest values (Figure 6). The lowest FRAP, TPC, and TFC values appeared in the seed samples harvested in the second year (2018), while samples from the 2017 season showed the highest. When analyzing changes within the same variety, it was observed that Q3 cv. seeds harvested in 2017 presented the highest phenolic content and those harvested in 2019 had the largest flavonoid content, although in this last year, Q3 cv. seeds showed an especially low phenolic content.
The color was determined as parameters related to the technological and functional quality of the seeds (Guiotto et al., 2020). In this analysis, the color values indicated that there were significant differences when comparing cultivars or among year of sowing. Both factors and their interaction had a significant influence on luminosity (L∗) and in hue value (h°), but only the cultivar factored in chroma (C∗, color saturation) and color distance (ΔE). In general, Regalona cv. and Puno cv. were the clearest (high luminosity (L∗), Table 3). This characteristic was significantly affected by the year of cultivation, being 2018 the year associated with darker seeds. Regarding cultivars, the a∗ parameter, that corresponds to red (+a∗) or green (-a∗) component of color, Titicaca cv. and Q5 cv. seeds were the reddest (Table 3 and Supplementary Figure 3). Although there were slight differences, 2018 was the year that yielded higher a∗ values. The b∗ component of the color (+b∗ yellow, -b∗ blue) ranged between 18.23 and 27.43, being Q5 cv., followed by Titicaca cv., Regalona cv., and Puno cv., the cultivars that yielded seeds associated with a highest yellow component.
Correlation, Linear Regression, Path Analysis, and Principal Components Analysis (PCA)
A Pearson’s correlation coefficient test was performed to analyze the correlation between variables (Supplementary Figure 4). All amino acids contents (except for cysteine and tryptophan contents) showed high correlation coefficients with total protein content and among them, besides showing high positive correlation with P and Zn contents and negative correlations with flavonoid and K contents, yield, seed viability, and germination rates. Flavonoid content also correlated with the phenolic content and the antioxidant capacity (r = 0.87), and with yield, germination, and seed viability rate, as well as with color parameters and some mineral contents, especially Zn (r = −0.6). Both, germination rate and seed viability presented high correlation coefficients with each other (r = 0.81) and with yield (r = 0.69). Seed viability and germination rate correlated negatively with protein (r = −0.65 and r = −0.58, respectively) and amino acids contents and with all the minerals except for Na. The correlation coefficients of Zn and K contents with seed viability (r = −0.71 and r = 0.75, respectively) and germination rate (r = −0.59 and r = 0.62, respectively) were also remarkable.
This analysis was followed up by simple linear regressions to evaluate the influence of yield on different seed nutritional quality-related traits (Table 4). Seed viability showed the largest linear correlation (R2 = 0.455) with yield, where yield explains 67.5% of the viability rate’s variance. K content also showed a large relation with yield (R2 = 0.293), being this trait responsible for 54.2% of the K content variance. The model obtained relating yield to total protein content (R2 = 0.289) predicted a decrease in seed protein content of 0.12% when the cultivar yield increases in 1 t/ha. Other qualitative seed traits such as germination rate, antioxidant capacity and flavonoid, phenols, Mg, Fe, and Zn contents proved to be influenced by yield, producing adequate linear models (p < 0.05). However, seed weight and P, Ca, Na, and Cu contents could not be explained by yield performance (p > 0.05).
Table 4. Regression analysis summary of predictive models for yield predicting qualitative seed variables.
Going further, path analysis was performed to define the direct and indirect contributions of each trait on seed germination rate. First, a predictive multiple linear regression model was performed following the stepwise method in order to find seed traits with a direct effect on germination rates. Individual amino acids contents were not included in the analysis because of their high multicollinearity. Yield, antioxidant capacity, phenols, flavonoids, protein, Zn, Fe, Ca, and Zn contents were eliminated from the model, while Mg, K, and Mn contents, entered the model as first-order predictors, explaining 55.8% of the germination rate variability (Table 5) and generating the following model:
Table 5. Direct effects of predictor variables of first-, second-, third-, and fourth-order on germination rate, tolerance and variance inflation factor of the path analysis.
(being the germination rate expressed as %, K and Mg content as % of seed weight, and Mn content as ppm of dry seed). K content showed a strong positive effect in the germination rate, accounting for 43% of the germination rate, while Mg, and Mn contents were in turn responsible for 30 and 28.2% of the germination rate’s variance, respectively, presenting a negative effect. Going further, TFC and Ca contents acted as second-order predictors of the germination rate through Mg content, having the earlier a negative effect on Mg and the latter a strong positive effect (61.5% of Mg’s variance). Ca was also a second-order predictor through Mn content, together with Zn, Cu, and phenols contents. Zn showed the largest direct effect on Mn (60.1% of Mn’s variance, Table 5). Through the K content path, two second-order predictors were found, phenols and Zn contents, showing the latter a strong negative contribution to K’s variance (68.2%, Table 5). The analysis also showed third-order predictors like Fe content (through Ca content), TFC (through a positive effect on TPC and a negative one on Zn), yield (through TPC), total protein content (through Zn), and antioxidant activity (via negative influence on Cu), as well as two fourth-order predictors, Fe content and Yield, which explained protein content with a positive effect of 30.9% and a negative effect of 43.2%, respectively, on protein’s variance (Table 5). No collinearity was found in the analysis.
Furthermore, a principal component analysis (PCA) was performed to reduce the number of variables (Figure 7). This analysis identified five principal components that were able to explain 69% of the variance. Component 1, which contributed to 37% of the variance, was mainly explained by the amino acid contents (except for cysteine, methionine, and tryptophan), yield, seed viability, and by the content of some minerals such as P, K, and Zn. In line with this, a strong positive correlation was found between the content of amino acids, P and Zn contents, which have high component 1 values, and a negative correlation of these variables with yield, seed viability, or K, with low component 1 values (Supplementary Figure 4 and Figure 7). Component 2 (explaining 11% of the variance) was accounted mainly for the antioxidant parameters (FRAP value, TPC, and TFC), yield, seed viability, the germination rate, and K content (high component 2 values), and for protein and Zn content, which correlated negatively with the other variables and showed low component 2 values (Supplementary Figure 4 and Figure 7). Component 3 (explaining 9% of the variance) considered the color parameters L∗, a∗, h°, the FRAP value, and the protein content. The color parameters C∗ and b∗ together with the Ca, Mg, and Fe contents and the germination rate, contributed to component 4 (explaining 6% of the variance). Component 5 (explaining 5% of the variance) included the content of P, Cu, Mn, and Zn, cysteine, and seed weight.
Figure 7. Principal components analysis. Biplot of main components 1 and 2 for the cultivars sown in each year of the experiment and for the variables tested. Component 1 (X axis) is contributed mainly by seed viability rate (VR) and yield, amino acid content (except methionine, cysteine, and tryptophan), P, K, and Zn contents. Component 2 (Y axis) includes yield, viability and germination rate (GR), FRAP value, and protein, tryptophan, polyphenols (TPC), flavonoids (TFC), K, and Zn contents.
By reducing the number of variables to main components it was possible to classify the quinoa genotypes for each year of harvesting in different groups (Figure 7). As consistently found in previous analyses, the year of cultivation was the most important factor when grouping the different variables studied. Seeds that belong to the 2017 harvest (except Puno cv. and Q5 cv.) formed the first group, presenting positive values for component 1 (meaning higher amino acids, P, and Zn contents, and lower K content, yield, and viability) and positive values for component 2 (related to higher antioxidants contents, viability, germination rate, tryptophan, and K content, and lower protein and Zn content). All 2018 seeds belonged to the second group and showed negative values for component 2, which explains the higher protein and Zn content and lower antioxidant capacity, germination, viability, and yield found that year. This group also shows positive values (except Q3 cv. and Puno cv.) for component 1, explaining the higher amino acids contents. Meanwhile, all 2019 seeds, comprising the third group, presented negative values for component 1, which coincides with the lower amino acids’ concentrations found that year. However, 2017 seeds from Puno cv. and Q5 cv. appeared closer in the PCA analysis to the seeds from the 2019 harvest than other 2017 seeds, showing also lower amino acid contents. For component 1, 2018 samples showed higher values than 2019, and Vikinga cv. seeds showed higher levels than Regalona cv., Puno cv., Q3 cv., or Q5 cv. For component 2, 2018 samples were significantly lower than 2019, and 2017, and Vikinga cv., Puno cv., and Q3 cv. had lower values than Titicaca cv., Regalona cv., or Q5 cv. Vikinga cv. seeds from 2019 showed the lowest component 2 level for 2019 samples and one of the lowest overall, which coincides with its poor yield performance, and its low viability and germination rate, and flavonoid and polyphenol contents.
Differences among cultivars were also observed and were reflected in Figure 7 through the PCA. Titicaca cv. seeds presented the highest values for component 2, and, in 2018, for component 1, although with larger variations. These seeds showed higher yields, germination rates, and antioxidant contents, although they were negatively affected by the environment in 2018. In contrast, Vikinga cv. seeds showed higher values for component 1 in 2017 and 2018 but very low values for component 2 in 2018 and 2019, being the only cultivar with lower values in 2019 seeds than in 2018 seeds. This can be seen in their consistently amino acids contents and the steeped decrease in yield, seed germination rates and viability, and antioxidants contents in both 2018 and 2019. On the other hand, Puno cv. seeds values are low for both components and show low germination rates, seed weights, and protein, amino acids, and antioxidants contents.
In the last decade, quinoa has acquired an increased agronomical and nutritional relevance related to the capacity of adaptation of this crop to different environments together with the exceptional nutritional properties of their seeds, which include high protein contents, an optimal amino acid balance, and an excellent antioxidant capacity, this later largely related to the high phenol content (Bazile et al., 2016; Jacobsen, 2017; Angeli et al., 2020). However, the establishment of this crop in many agronomical areas outside South America is still limited. It could be considered that quinoa cultivar selection process remains unfinished for new cultivation areas, including those located in Southern Europe. Furthermore, although the potential of this crop has been comprehensively analyzed in nearby areas (Jacobsen, 2017) there is still very limited information regarding the stability of seed nutritional characteristics under changing environments.
Multiple environmental factors, such as temperature, water status, photoperiod and light quality, and soil nutrient content, together with genetic features, are responsible for determining the quality of seeds. Other factors like the physiological status of the plant during growth or postharvest parameters such as moisture and temperature during the storage of seeds may play pivotal roles in determining quality as well. Ultimately, genotype, environment, and their interaction are the main factors determining the status of seeds (Hakeem, 2015). In this study, six different cultivars were used in order to examine the effect of the genotype on different physiological and nutritional traits of quinoa seeds, which were sown in three consecutive years aiming to analyze the environmental effect on those parameters.
Saponins are antinutrients that can diminish the nutritional value of quinoa seeds. These compounds can alter the absorption of minerals such as Fe and Zn (Ruales and Nair, 1993) and they give bitterness to the seeds decreasing their palatability. Therefore, extensive efforts have been made through breeding programs toward reducing their concentrations in seeds (Mastebroek et al., 2000; Zurita-Silva et al., 2014). Nevertheless, the increasing evidence regarding their multiple bioactivities for health may be worth considering, such as anti-inflammatory, antitumor, hypocholesterolemic, or immunomodulatory (Marrelli et al., 2016; Singh et al., 2017; Navarro del Hierro et al., 2018; El Hazzam et al., 2020). The limit established to classify quinoa varieties as sweet or bitter is 0.11% of saponin per seed fresh weight (Koziol, 1991). Thus, Titicaca cv. (used in this study) can be classified as bitter, as had been previously reported (Medina-Meza et al., 2016), and, on the other hand, Vikinga cv. could be considered a “low saponin” variety (Figure 4; Medina-Meza et al., 2016). The results showed that the main difference in saponin contents was determined by the cultivar, although higher contents in 2019 Vikinga cv. seeds could be observed compared to previous years (Figure 4). As previously described, 2019 was the year that showed a reduced water supply (Supplementary Table 1 and Figure 1). This finding is in agreement with that of Reguera et al. (2018), who suggested that saponin content is mainly a genotypic-dependent trait, although it can change as well under stress conditions including drought or soil salinity (Gómez-Caravaca et al., 2012; Pulvento et al., 2012). In fact, according to Figure 4, it seems that specific environmental conditions might even lead to the production of bitter quinoa seeds but from sweet quinoa varieties. Furthermore, it should not be ruled out an influence of the genetic purity of the seeds used in this study in the variations observed in the saponin content, as it may reflect genetic segregation or genetic cross-contamination, which might as well explain the reduced Vikinga cv. seed yield observed in 2019. Therefore, further research to understand the mechanisms or environmental conditions that impact the content of saponins of quinoa seeds would be of interest. This would benefit future breeding efforts if they were to target, either to reach low-saponin quinoa seeds, or even saponin-enriched seeds as a selection trait, considering the current popularity of these phytochemicals as bioactive compounds.
Quinoa seeds are well known for possessing high protein contents. This study shows a protein range that went from 14.1 to 18.6% of seed weight, with an average of 16.5% protein per seed weight (Figure 5). These results were similar to those obtained in other works on quinoa (Koziol, 1992; Miranda et al., 2012; Reguera et al., 2018) and reflect higher protein contents compared to important staple cereal crops such as barley (11%), wheat (10%), maize (14%), or rice (8%) (Koziol, 1992; Nowak et al., 2016; Filho et al., 2017). However, not only protein quantity but also protein quality may impact human diet. Most proteins of plant origin have very low levels of essential amino acids, especially tryptophan, methionine, and lysine compared to those of animal origin (Friedman, 1996). Particularly, lysine and tryptophan are present in low levels in grains of cereal crops like barley and wheat, while sulfur amino acids (methionine and cysteine) show lower levels in legume seeds including soybean or beans (Friedman, 1996). In line with this, quinoa seeds are considered a “complete protein” source since they provide all amino acids essential for human consumption. They have higher levels of lysine, methionine, and cysteine compared to cereals or legumes which makes quinoa a great food complement in healthy diets (Koziol, 1992; Repo-Carrasco et al., 2003; Abugoch James, 2009; Filho et al., 2017). However, different studies have reported limiting essential amino acid contents in quinoa according to the daily requirements established by the Food Agriculture Organization (FAO) (World Health Organization and United Nations University, 2007; Gonzalez et al., 2012; Miranda et al., 2012; Präger et al., 2018; Craine and Murphy, 2020). In this work, all samples analyzed met the daily lysine and leucine requirements established by the FAO for all age groups (World Health Organization and United Nations University, 2007; Table 1). All samples from 2017 harvest and 2019 Puno cv. and Q5 cv. met the daily requirements of tryptophan for adults, but only 2017 Titicaca cv. and Q5 cv. seeds contained enough tryptophan to meet the requirements for children (World Health Organization and United Nations University, 2007; Table 1). The other samples did not meet the requirements of either group. In the case of sulfur amino acids (combining cysteine and methionine contents), only 2017 seeds of all cultivars and those from Puno cv. and Q3 cv. harvested in 2018 met the daily children and adults’ recommendations, especially due to the low levels of methionine (Table 1). Both methionine and tryptophan may be limiting in quinoa (Mahoney et al., 1975; Gonzalez et al., 2012; Craine and Murphy, 2020). Therefore, overall, and based on these and previous analysis (Craine and Murphy, 2020), increased lysine and sulfur amino acids levels should be breeding targets in quinoa toward enhancing quality to reach Food Security.
It has been hypothesized that protein and amino acids contents can be determined by nitrogen availability in the soil, the environmental and agroecological conditions, and the genotype (Thanapornpoonpong et al., 2008; Craine and Murphy, 2020). In this study, a significant variability throughout samples was found in both protein and amino acid contents (Figure 5 and Table 1). While the evidences presented by Miranda et al. (2012) suggest a strong influence of the genotype in the amino acids contents, the studies performed by Reguera et al. (2018) and Prieto et al. (2021) support the hypothesis that environmental factors can influence the protein content of quinoa seeds. Reguera et al. (2018) found differences in protein quantity among seeds harvested in different countries but not among cultivars in a certain location, while Prieto et al. (2021) showed an important increase in quinoa protein content when the crop underwent heat stress. In the present study, there is a strong influence in both protein and amino acids contents by the year of cultivation, being 2018 the growing season with higher contents except for cysteine and tryptophan (Figure 5 and Table 1). Some amino acids did not show significant changes among cultivars, but in most cases Vikinga cv. seeds showed higher contents and Puno cv. seeds lower. It is worth mentioning that Titicaca cv. seeds showed larger differences for most amino acids, being generally 2018 seeds richer, while Regalona cv. and Q3 cv. seeds had a more stable amino acid content among years of harvest. This differential response from different genotypes to changing environmental conditions had already been observed by Präger et al. (2018), who described higher essential amino acid contents in Jessie cv. and Zeno cv. in 2016 compared to 2015, but not in Puno cv. or Titicaca cv. Thus, the modulation of protein and amino acid contents by environmental conditions is largely dependent on the genotype. Understanding the mechanisms responsible for these changes in the amino acidic profile will be of importance in the study of stress tolerance in different quinoa cultivars since some essential amino acids also play a role as osmolytes. For instance, the accumulation of branched-chain amino acids (BCAAs) (valine, leucine, isoleucine) is induced by osmotic stress (Joshi et al., 2010).
Mineral content is an important determinant of seed quality. Quinoa seeds are characterized by presenting a high content of Ca, Mg, Fe, Cu, and Zn; moreover, Ca, Mg, and K are found in sufficient quantities in quinoa seeds to meet a balanced human diet, since they are in bioavailable forms (Repo-Carrasco et al., 2003). The seeds here analyzed met Fe daily requirements for all age groups, even for women at a menstruating age (World Health Organization, 2004), but K contents did not meet the requirements for pregnant or lactating women (Turck et al., 2016). The mineral contents here presented are similar to previous reports in quinoa (Table 2; Aranda et al., 2013; Prado et al., 2014; Reguera et al., 2018). In fact, Regalona cv. and Titicaca cv. levels matched in range with those reported by Reguera et al. (2018), who used these same varieties, but not for Na and Ca contents (higher in the present study). In the preceding study, stark differences were found between sowing locations, which were attributed to differences in soil composition (Reguera et al., 2018). In this study, both genotype and environment influenced mineral contents (except for Na), even though all seeds were sown at the same location in the three consecutive years. For instance, contents of P, Mg, Fe, Mn, and Zn were significantly higher in Vikinga cv. seeds and in 2018 seeds (Table 2). It should be noted that further analysis should be made to determine if the increased mineral content is correlated with mineral bioavailability and to determine the role of the genotype and environment regulating this aspect, considering the differences previously observed in quinoa (Vidueiros et al., 2015) and the effect of components like saponins or phytic acid on the mineral’s bioavailability in quinoa (Ruales and Nair, 1993).
Also, antioxidants may condition the nutritional quality and shelf life of seeds. From a nutritional point of view, these compounds add health benefits as they can reduce the risk of cancer and cardiovascular diseases (Vega-Gálvez et al., 2010; Tang and Tsao, 2017). In line with this, quinoa seeds are an excellent source of antioxidants, especially because of their high contents of phytochemicals like polyphenols, flavonoids, and vitamin E (Tang and Tsao, 2017), which also exceed the levels found in cereals (Gorinstein et al., 2007). The contents of these phytochemicals are variable in quinoa and are genotype-dependent. Besides, there is evidence of changes related to environmental constraints such as salt or water stress reflecting the environmental control in the synthesis of these compounds (Fischer et al., 2013; Aloisi et al., 2016; Ismail et al., 2016). In the present study, the amounts of polyphenols, flavonoids, and FRAP capacity were found at comparable levels to those reported by Ismail et al. (2016), Paśko et al., 2009, Abderrahim et al. (2015), and Fischer et al. (2013). Both genotype and environmental conditions were determinant factors of these parameters, being lower in seeds harvested in 2018 and in Regalona cv. and Puno cv. seeds, while higher in Titicaca cv. (Figure 6). Nonetheless, it should be highlighted that the response to different environmental conditions seems to differ between cultivars, and Titicaca cv. seeds have a steeper decrease in FRAP value and TFC compared to Vikinga cv. seeds (Figure 6). This supports the hypothesis that the antioxidant capacity depends on the genotype, the environment, and the interaction of these factors. Both, Ismail et al. (2016) and Aloisi et al. (2016) propose that those varieties more tolerant to stress require a lower production of antioxidants. For this reason, more research needs to be performed investigating the mechanisms responsible for the changes in the antioxidants of quinoa seeds. This would benefit future breeding efforts if they were to target antioxidant capacity as a selection trait.
The content of phenols and other secondary metabolites can affect seed parameters such as seed color (Ballester-Sánchez et al., 2019). Similarly, seed color might be indicative of the content of these compounds. In this study color differences were related to the cultivar, being Titicaca cv. and Q5 cv. seeds darker and redder, and Regalona cv. and Puno cv. seeds lighter (Table 3 and Supplementary Figure 3). The year of cultivation only affected the a∗ component of color (redness), being 2018 seeds the reddest. The L∗ component of color negatively correlated with FRAP and TPC levels (Supplementary Figure 4), i.e., darker seeds provided more antioxidants. This characteristic had been previously reported in quinoa seeds and in other plant species (Szydłowska-Czerniak et al., 2011; Abderrahim et al., 2015).
Yield is a common selection criterion for quinoa breeding programs, aiming to increase productivity (Zurita-Silva et al., 2014). The varieties cultivated in this study yielded between 0.7 and 3.25 t/ha, and fell within the ranges previously reported for quinoa (Supplementary Figure 1; Bertero et al., 2004; Zurita-Silva et al., 2014; Choukr-Allah et al., 2016; Lesjak and Calderini, 2017). In this study, the year of harvest was an important determinant of yield, with yields plummeting in 2018 in most cultivars. Q3 cv. and Puno cv. were less affected by environmental changes among years compared to the rest of cultivars. Regalona cv. experienced larger yield penalties in the present study than in previous works assessing drought stress (Fischer et al., 2013). Thus, besides the genetic background, other factors might be responsible for lowering yields of quinoa under field conditions, which might include high night temperatures or high temperatures combined with low precipitations or long photoperiods during sensitive phases that go from flowering to the end of grain filling stage (Bertero et al., 1999; Stikic et al., 2012; Lesjak and Calderini, 2017). Furthermore, on average, seed size showed larger values in 2018 (Figure 2B). Seed size might impact yields, and the environmental conditions seem to play an important role in controlling it. For instance, in quinoa, temperatures are shown to influence seed size during grain filling, especially in certain genotypes, which fits well with the results here presented (Bertero, 2021).
Yield can also be an important trait for breeding programs when pursuing improvements in quinoa seed quality. In fact, yield is positively correlated with fat and fiber contents and lower protein contents (De Santis et al., 2016; Curti et al., 2018). Besides, the relationship between yield and nutritional quality-related traits such as the protein content, is influenced by the genotype and the environment, in quinoa and other crops (Caballero et al., 2015; Halford et al., 2015; Curti et al., 2018). For this reason, a deeper analysis of yield influence in seed quality traits was performed (Table 4) and an important negative effect of yield on protein content was found, agreeing with the trade-off described by Curti et al. (2018) for winter sowing and with other similar findings in quinoa (De Santis et al., 2016; Präger et al., 2018; Reguera et al., 2018; Prieto et al., 2021) and other crops, including cereals (Simmonds, 1995; Rondanini et al., 2019). Nonetheless, both, yield and protein content, have also been shown to increase simultaneously when fertilizing with N (Caballero et al., 2015). Likewise, yield also correlated with most seed quality traits, but not with morphologic ones like 1000 seed weight (Table 4 and Supplementary Figure 4), indicating that seed weight is not the main contributor to yield performance (Curti et al., 2014, 2018).
Nonetheless, a correlation analysis provides a limited view of the complex interrelation that can occur between different seed traits, not showing directionality or indirect effects on other nutrient levels (Dewey and Lu, 1959). A sequential path analysis allows the classification of different variables as first-order, second-order, third-order (and so on) predictors of a response variable, so the different intercorrelations can be unraveled. Other studies have performed path analysis to study the contributions of physiological crop traits to yield performance in quinoa and other crops (Mohammadi et al., 2003; Bhargava et al., 2007; Mubai et al., 2020). However, to our knowledge, this is the first study attempting to explain quinoa seed germination through other seed characteristics. The model postulated in this work suggests that different mineral contents play an important role in determining germination rates of quinoa seeds, mainly K content, positively, and Mn and Mg contents, negatively. Different studies have linked mineral nutrition with seed germination capacity, including fertilization treatments with K increasing the germination rates of cotton (Sawan et al., 2011), or seedling growth inhibition in horse gram, after applying high concentrations of Mn (Kumari et al., 2016). However, none of these studies paid attention to the seed content of these elements. Noteworthy, it is important to highlight that the path model here presented is designed within the context of a particular environment and may work differently in other areas of cultivation. For instance, when growing quinoa in high salinity soil, the antioxidant contents may play a major role (not a secondary one) in germination capacity, where a lower germination rate may be achieved in seeds that accumulate K and Mg (Koyro and Eisa, 2008; Panuccio et al., 2014). Besides, environmental factors such as temperature, photoperiod lengths, and precipitations (throughout the crop life cycle) could influence germination rates and seed viability. Thus, the lower germination rates and seed viability observed in 2018 could be related to the shorter photoperiods after anthesis compared to 2017 and 2019 and/or with the lower temperatures at seed filling stage. Nonetheless, the specific role of temperatures, photoperiod lengths, and/or high precipitations regulating germination and viability should be further analyzed.
Overall, differences depending on genotype and on the environmental factors together with the genotype x environment interaction were found in most parameters measured (Supplementary Table 3). The cultivars analyzed had different genetic backgrounds and were grown under different environmental contexts. These variations can be seen in Figure 7, where variables have been reduced to principal components through a PCA. The factor that made the difference between seeds was the year of cultivation, as Figure 7 shows, being 2019 seeds low on component 1, 2017 seeds high on component 2, and 2018 seeds low on component 2 and high on component 1. This means that 2018 seeds show low yields and antioxidants, but high protein, amino acids, P, Cu, Mn, and Zn contents, while 2019 seeds have lower amino acids contents but higher yields and germination rates.
In 2018, the sowing date was delayed because of heavy rainfalls, so in the first two months of development, plants were exposed to higher temperatures than those sown in 2017 or 2019 (Supplementary Table 1 and Figure 1). All three growing seasons showed long photoperiods, which peaked in June at 15 h long days. Plants in 2017, 2018, and 2019 were sown with daylengths of 13 h, 15 h, and 14 h, respectively; they reached the highest photoperiod in their second month, first week, and first month, and were harvested at daylengths of 13 h, 10.75 h, and 12.25 h, respectively (Figure 1). Quinoa is a facultative short-day plant where photoperiods longer than 12 h can disrupt seed filling and maturation, although day-length neutral varieties have been developed in order to introduce the crop to higher latitudes (Bendevis et al., 2014; Zurita-Silva et al., 2014). The combination of long photoperiods and high temperatures can cause yield penalties in quinoa, especially during flowering (Bertero et al., 1999; Bertero, 2003; Lesjak and Calderini, 2017). The combination of these factors was observed in 2018 (Supplementary Table 1 and Figure 1) which could explain the poorer yield performance. Furthermore, according to the path model presented in this work, lower yields would have caused higher protein contents in seeds (Figure 5) but, ultimately, lower germinations (Figure 3).
It is important to point out that an important nutritional seed trait such as the protein content negatively correlated with yield, which is an important criterion for breeding programs (Zurita-Silva et al., 2014), and with the antioxidant capacity, which is another very interesting trait for human nutrition due to the accompanying health benefits (Figure 7 and Supplementary Figure 4; Vega-Gálvez et al., 2010). Thus, from an agronomical point of view, Puno cv. could be the best choice for cultivation, since it was the cultivar that better performed in 2018 in comparison with the rest of cultivars. However, that would have meant having smaller seeds with low values for most of the nutritional parameters. Therefore, from a nutritional perspective, a more appropriate cultivar for this area of study would be Titicaca cv., which showed higher protein and amino acids contents, especially in 2018, and also a greater antioxidant capacity.
These differences between cultivars and their performance according to the environment, together with the trade-offs among important crop characteristics, pose a challenge for breeders, who will need to study closely different cultivars for each location of cultivation. Therefore, the information presented in this work will greatly help the efforts of quinoa establishment in Northwestern Spain. This also highlights the need for in-depth research to unravel the mechanisms that cause the variations observed in the nutritional traits due to changes in the environmental conditions and agroecological contexts, with the ultimate goal of obtaining better adapted and more nutritious seeds toward contributing to food security worldwide.
The results here presented highlight a great influence of the environmental conditions on the nutritional and physiological characteristics of quinoa seeds, which affects overall seed quality. Particularly, this work has shown the important effect of the environment on the amino acid balance and content, the impact of the genotype on the antioxidant capacity, protein amount, the negative correlation between protein and antioxidant contents and the existence of stable nutritional components such as the Na content. Indeed, plants grown during the second year (2018) showed lower yields and heavier seeds presenting worse germination powers. These seeds also presented higher amino acids, phosphorous, copper, manganese, and zinc contents, and lower potassium and antioxidants. According to the analyses performed, yields were associated with seed viability and protein, phenol, K, and Zn contents. Furthermore, germination rates were found to be directly influenced by K, Mg, and Mn seed contents. It is expected that these findings will help to maximize quinoa productivity and/or nutritional quality, especially for comparable climatic areas of analysis. Moreover, although the main goal of this study was to evaluate the impact of the environment on quality-related traits, it also highlights that there are still important limitations in the agronomical adaption of quinoa to these areas of cultivation which are characterized by having intense precipitations at early growth stages and high temperatures at later stages of the crop development. In line with this, we found that agronomically, varieties such as Puno cv. or Q3 cv. might be better adapted to these conditions, but Titicaca cv. and Vikinga cv. showed better nutritional properties, as they possessed higher protein contents. Altogether, this study supports the huge potential of this crop by choosing the appropriate variety according to the area of interest. Therefore, the selection of the cultivar must be well informed, paying careful attention to how the seeds respond to the climatological characteristics in a particular location.
Data Availability Statement
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author/s.
MR, NA, LFP-R, NF-G, and JM conceived and planned the experiments. SG-R, LFP-R, IG, JM, NA JP, and CH carried out the experiments. SG-R, MR, IM, JM, NA, JP, CH, NF-G, and LB contributed to the interpretation of the results. MR and SG-R took the lead in writing the manuscript. All authors provided critical feedback and helped shape the research, analysis, and manuscript.
This work was supported by the Ministerio de Ciencia e Innovación (MICINN, Spain) (PID2019-105748RA-I00), the Madrid Government (Comunidad de Madrid-Spain) under the Multiannual Agreement with Universidad Autónoma de Madrid in the line of action encouraging youth research doctors, in the context of the V PRICIT (Regional Programme of Research and Technological Innovation) (SI1/PJI/2019-00124), the CYTED (ValSe-Food 119RT0567), the FPI UAM Fellowship Programme 2019 (to SG-R), and the Ramón y Cajal Programme 2019 (to MR).
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
The authors greatly thank Sven Jacobsen (Quinoa Quality), the International Center for Biosaline Agriculture (ICBA), and Semillas Baer (Chile), for providing the quinoa seeds used in this study.
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2021.649132/full#supplementary-material
Abderrahim, F., Huanatico, E., Segura, R., Arribas, S., Gonzalez, M. C., and Condezo-Hoyos, L. (2015). Physical features, phenolic compounds, betalains and total antioxidant capacity of coloured quinoa seeds (Chenopodium quinoa Willd.) from Peruvian Altiplano. Food Chem. 183, 83–90. doi: 10.1016/j.foodchem.2015.03.029
Adolf, V. I., Shabala, S., Andersen, M. N., Razzaghi, F., and Jacobsen, S. E. (2012). Varietal differences of quinoa’s tolerance to saline conditions. Plant Soil 357, 117–129. doi: 10.1007/s11104-012-1133-7
Alandia, G., Rodriguez, J. P., Jacobsen, S.-E., Bazile, D., and Condori, B. (2020). Global expansion of quinoa and challenges for the Andean region. Glob. Food Sec. 26:100429. doi: 10.1016/j.gfs.2020.100429
Aloisi, I., Parrotta, L., Ruiz, K. B., Landi, C., Bini, L., Cai, G., et al. (2016). New Insight into Quinoa Seed Quality under Salinity: Changes in Proteomic and Amino Acid Profiles, Phenolic Content, and Antioxidant Activity of Protein Extracts. Front. Plant Sci. 7:656. doi: 10.3389/fpls.2016.00656
Angeli, V., Silva, P. M., Massuela, D. C., Khan, M. W., Hamar, A., Khajedi, F., et al. (2020). Quinoa (Chenopodium quinoa Willd.): An Overview of the Potentials of the “Golden Grain” and Socio-Economic and Environmental Aspects of Its Cultivation and Marketization. Foods 9, 1–31. doi: 10.3390/foods9020216
Aranda, M., Fuentes, F., Miranda, M., Vega-gálvez, A., Martínez, E. A., López, J., et al. (2013). Influence of contrasting environments on seed composition of two quinoa genotypes: nutritional and functional properties. Chil. J. Agric. Res. 73, 108–116. doi: 10.4067/S0718-58392013000200004.1
Ballester-Sánchez, J., Gil, J. V., Fernández-Espinar, M. T., and Haros, C. M. (2019). Quinoa wet-milling: Effect of steeping conditions on starch recovery and quality. Food Hydrocoll. 89, 837–843. doi: 10.1016/j.foodhyd.2018.11.053
Bascuñán-Godoy, L., Reguera, M., Abdel-tawab, Y. M., and Blumwald, E. (2015). Water deficit stress-induced changes in carbon and nitrogen partitioning in Chenopodium quinoa Willd. Planta 243, 591–603. doi: 10.1007/s00425-015-2424-z
Bendevis, M. A., Sun, Y., Shabala, S., Rosenqvist, E., Liu, F., and Jacobsen, S. E. (2014). Differentiation of Photoperiod-Induced ABA and Soluble Sugar Responses of Two Quinoa (Chenopodium quinoa Willd.) Cultivars. J. Plant Growth Regul. 33, 562–570. doi: 10.1007/s00344-013-9406-9
Bertero, H. D., De La Vega, A. J., Correa, G., Jacobsen, S. E., and Mujica, Á (2004). Genotype and genotype-by-environment interaction effects for grain yield and grain size of quinoa (Chenopodium quinoa Willd.) as revealed by pattern analysis of international multi-environment trials. F. Crop. Res. 89, 299–318. doi: 10.1016/j.fcr.2004.02.006
Bhargava, A., Shukla, S., and Ohri, D. (2007). Genetic variability and interrelationship among various morphological and quality traits in quinoa (Chenopodium quinoa Willd.). F. Crop. Res. 101, 104–116. doi: 10.1016/j.fcr.2006.10.001
Caballero, A., Maceda, W., Miranda, R., and Bosque, H. (2015). Rendimiento y contenido de proteína de la quinua (Chenopodium quinoa Willd), en cinco fases fenológicas, bajo cuatro niveles de incorporación de estiércol. Rev. Investig. Innov. Agrop. Recurs. Natur. 2015, 68–75.
Choukr-Allah, R., Rao, N. K., Hirich, A., Shahid, M., Alshankiti, A., Toderich, K., et al. (2016). Quinoa for Marginal Environments: Toward Future Food and Nutritional Security in MENA and Central Asia Regions. Front. Plant Sci. 7:346. doi: 10.3389/fpls.2016.00346
Curti, R. N., de la Vega, A. J., Andrade, A. J., Bramardi, S. J., and Bertero, H. D. (2014). Multi-environmental evaluation for grain yield and its physiological determinants of quinoa genotypes across Northwest Argentina. F. Crop. Res. 166, 46–57. doi: 10.1016/j.fcr.2014.06.011
Curti, R. N., Sanahuja, M., del, C., Vidueiros, S. M., Pallaro, A. N., and Bertero, H. D. (2018). Trade-off between seed yield components and seed composition traits in sea level quinoa in response to sowing dates. Cereal Chem. 95, 734–741. doi: 10.1002/cche.10088
De Santis, G., Maddaluno, C., D’Ambrosio, T., Rascio, A., Rinaldi, M., and Troisi, J. (2016). Characterisation of quinoa (Chenopodium quinoa willd.) accessions for the saponin content in Mediterranean environment. Ital. J. Agron. 11, 277–281. doi: 10.4081/ija.2016.774
Dewey, D. R., and Lu, K. H. (1959). A Correlation and Path−Coefficient Analysis of Components of Crested Wheatgrass Seed Production. Agron. J. 51, 515–518. doi: 10.2134/agronj1959.00021962005100090002x
El Hazzam, K., Hafsa, J., Sobeh, M., Mhada, M., Taourirte, M., Kacimi, K. E. L., et al. (2020). An insight into saponins from Quinoa (Chenopodium quinoa Willd): A review. Molecules 25:1059. doi: 10.3390/molecules25051059
Filho, A. M. M., Pirozi, M. R., Borges, J. T. D. S., Pinheiro Sant’Ana, H. M., Chaves, J. B. P., and Coimbra, J. S. D. R. (2017). Quinoa: Nutritional, functional, and antinutritional aspects. Crit. Rev. Food Sci. Nutr. 57, 1618–1630. doi: 10.1080/10408398.2014.1001811
Fischer, S., Wilckens, R., Jara, J., and Aranda, M. (2013). Variation in antioxidant capacity of quinoa (Chenopodium quinoa Will) subjected to drought stress. Ind. Crops Prod. 46, 341–349. doi: 10.1016/j.indcrop.2013.01.037
Gómez-Caravaca, A. M., Iafelice, G., Lavini, A., Pulvento, C., Caboni, M. F., and Marconi, E. (2012). Phenolic compounds and saponins in quinoa samples (Chenopodium quinoa Willd.) grown under different saline and nonsaline irrigation regimens. J. Agric. Food Chem. 60, 4620–4627. doi: 10.1021/jf3002125
Gonzalez, J. A., Konishi, Y., Bruno, M., Valoy, M., and Prado, F. E. (2012). Interrelationships among seed yield, total protein and amino acid composition of ten quinoa (Chenopodium quinoa) cultivars from two different agroecological regions. J. Sci. Food Agric. 92, 1222–1229. doi: 10.1002/jsfa.4686
Gorinstein, S., Vargas, O. J. M., Jaramillo, N. O., Salas, I. A., Ayala, A. L. M., Arancibia-Avila, P., et al. (2007). The total polyphenols and the antioxidant potentials of some selected cereals and pseudocereals. Eur. Food Res. Technol. 225, 321–328. doi: 10.1007/s00217-006-0417-7
Halford, N. G., Curtis, T. Y., Chen, Z., and Huang, J. (2015). Effects of abiotic stress and crop management on cereal grain composition: Implications for food quality and safety. J. Exp. Bot. 66, 1145–1156. doi: 10.1093/jxb/eru473
Ismail, H., Dragišic Maksimovic, J., Maksimovic, V., Shabala, L., Branka, D. Ž, Tian, Y., et al. (2016). Rutin, a fl avonoid with antioxidant activity, improves plant salinity tolerance by regulating K+ retention and Na+ exclusion from leaf mesophyll in quinoa and broad beans. Function 43, 75–86. doi: 10.1071/FP15312
Jacobsen, S. E., Monteros, C., Christiansen, J. L., Bravo, L. A., Corcuera, L. J., and Mujica, A. (2005). Plant responses of quinoa (Chenopodium quinoa Willd.) to frost at various phenological stages. Eur. J. Agron. 22, 131–139. doi: 10.1016/j.eja.2004.01.003
Joshi, V., Joung, J. G., Fei, Z., and Jander, G. (2010). Interdependence of threonine, methionine and isoleucine metabolism in plants: Accumulation and transcriptional regulation under abiotic stress. Amino Acids 39, 933–947. doi: 10.1007/s00726-010-0505-7
Kumari, M. K., Varaprasad, D., Narasimham, D., Paramesh, K., and Chandrasekhar, T. (2016). Impacts of Cadmium and Manganese on in Vitro Seed Germination and Seedling Growth of Horsegram. Ind. J. Plant Sci. 5, 119–125.
Lesjak, J., and Calderini, D. F. (2017). Increased night temperature negatively affects grain yield, biomass and grain number in Chilean quinoa. Front. Plant Sci. 8, 1–11. doi: 10.3389/fpls.2017.00352
Marrelli, M., Conforti, F., Araniti, F., and Statti, G. A. (2016). Effects of saponins on lipid metabolism: A review of potential health benefits in the treatment of obesity. Molecules 21:1404. doi: 10.3390/molecules21101404
Mastebroek, H. D., Limburg, H., Gilles, T., and Marvin, H. J. P. (2000). Occurrence of sapogenins in leaves and seeds of quinoa (Chenopodium quinoa Willd). J. Sci. Food Agric. 80, 152–156. doi: 10.1002/(SICI)1097-0010(20000101)80
Medina-Meza, I. G., Aluwi, N. A., Saunders, S. R., and Ganjyal, G. M. (2016). GC-MS Profiling of Triterpenoid Saponins from 28 Quinoa Varieties (Chenopodium quinoa Willd.) Grown in Washington State. J. Agric. Food Chem. 64, 8583–8591. doi: 10.1021/acs.jafc.6b02156
Miranda, M., Vega-Gálvez, A., Martinez, E., López, J., Rodríguez, M. J., Henríquez, K., et al. (2012). Genetic diversity and comparison of physicochemical and nutritional characteristics of six quinoa (Chenopodium quinoa willd.) genotypes cultivated in Chile. Food Sci. Technol. 32, 835–843. doi: 10.1590/s0101-20612012005000114
Miranda, M., Vega-gálvez, A., Martínez, E. A., López, J., Marín, R., Aranda, M., et al. (2013). Influence of contrasting environments on seed composition of two quinoa genotypes: nutritional and functional properties. Chil. J. Agric. Res. 73, 108–116. doi: 10.4067/S0718-58392013000200004.1
Mohammadi, S. A., Prasanna, B. M., and Singh, N. N. (2003). Sequential path model for determining interrelationships among grain yield and related characters in maize. Crop Sci. 43, 1690–1697. doi: 10.2135/cropsci2003.1690
Mubai, N., Sibiya, J., Mwololo, J., Musvosvi, C., Charlie, H., Munthali, W., et al. (2020). Phenotypic correlation, path coefficient and multivariate analysis for yield and yield-associated traits in groundnut accessions. Cogent Food Agric. 6:1823591. doi: 10.1080/23311932.2020.1823591
Nascimento, A. C., Mota, C., Coelho, I., Gueifão, S., Santos, M., Matos, A. S., et al. (2014). Characterisation of nutrient profile of quinoa (Chenopodium quinoa), amaranth (Amaranthus caudatus), and purple corn (Zea mays L.) consumed in the North of Argentina: Proximates, minerals and trace elements. Food Chem. 148, 420–426. doi: 10.1016/j.foodchem.2013.09.155
Navarro del Hierro, J., Herrera, T., García-Risco, M. R., Fornari, T., Reglero, G., et al. (2018). Ultrasound-assisted extraction and bioaccessibility of saponins from edible seeds: quinoa, lentil, fenugreek, soybean and lupin. Food Res. Int. 109, 440–447. doi: 10.1016/j.foodres.2018.04.058
Navarro del Hierro, J., Reglero, G., and Martin, D. (2020). Chemical characterization and bioaccessibility of bioactive compounds from saponin−rich extracts and their acid−hydrolysates obtained from fenugreek and quinoa. Foods 9:1159. doi: 10.3390/foods9091159
Panuccio, M. R., Jacobsen, S. E., Akhtar, S. S., and Muscolo, A. (2014). Effect of saline water on seed germination and early seedling growth of the halophyte quinoa. AoB Plants 6, 1–18. doi: 10.1093/aobpla/plu047
Paśko, P., Bartoń, H., Zagrodzki, P., Gorinstein, S., Fołta, M., and Zachwieja, Z. (2009). Anthocyanins, total polyphenols and antioxidant activity in amaranth and quinoa seeds and sprouts during their growth. Food Chem. 115, 994–998. doi: 10.1016/j.foodchem.2009.01.037
Peñas, E., Uberti, F., di Lorenzo, C., Ballabio, C., Brandolini, A., and Restani, P. (2014). Biochemical and Immunochemical Evidences Supporting the Inclusion of Quinoa (Chenopodium quinoa Willd.) as a Gluten-free Ingredient. Plant Foods Hum. Nutr. 69, 297–303. doi: 10.1007/s11130-014-0449-2
Prado, F. E., Fernández-Turiel, J. L., Tsarouchi, M., Psaras, G. K., and González, J. A. (2014). Variation of seed mineral concentrations in seven quinoa cultivars grown in two agroecological sites. Cereal Chem. 91, 453–459. doi: 10.1094/CCHEM-08-13-0157-R
Präger, A., Munz, S., Nkebiwe, P. M., Mast, B., and Graeff-Hönninger, S. (2018). Yield and quality characteristics of different quinoa (Chenopodium quinoa Willd.) cultivars grown under field conditions in southwestern Germany. Agronomy 8:197. doi: 10.3390/agronomy8100197
Prieto, J., Rodríguez, M. J., Cruz, V., Calvo, P., and Reguera, M. (2021). Heat stress lowers yields, alters nutrient uptake, and changes seed quality in quinoa grown under Mediterranean field conditions. J. Agron. Crop Sci. 1–11. doi: 10.1111/jac.12495
Pulvento, C., Riccardi, M., Lavini, A., D’Andria, R., Iafelice, G., and Marconi, E. (2010). Field Trial Evaluation of Two Chenopodium quinoa Genotypes Grown Under Rain-Fed Conditions in a Typical Mediterranean Environment in South Italy. J. Agron. Crop Sci. 196, 407–411. doi: 10.1111/j.1439-037X.2010.00431.x
Pulvento, C., Riccardi, M., Lavini, A., Iafelice, G., Marconi, E., and D’Andria, R. (2012). Yield and Quality Characteristics of Quinoa Grown in Open Field Under Different Saline and Non-Saline Irrigation Regimes. J. Agron. Crop Sci. 198, 254–263. doi: 10.1111/j.1439-037X.2012.00509.x
Reguera, M., Conesa, C. M., Gil-Gómez, A., Haros, C. M., Pérez-Casas, M. Á, Briones-Labarca, V., et al. (2018). The impact of different agroecological conditions on the nutritional composition of quinoa seeds. PeerJ 6:e4442. doi: 10.7717/peerj.4442
Repo-Carrasco, R., Espinoza, C., and Jacobsen, S. E. (2003). Nutritional value and use of the andean crops quinoa (Chenopodium quinoa) and kañiwa (Chenopodium pallidicaule). Food Rev. Int. 19, 179–189. doi: 10.1081/FRI-120018884
Rondanini, D. P., Borrás, L., and Savin, R. (2019). “Improving Grain Quality in Oil and Cereal Crops,” in Crop Science, eds R. Savin and A. Slafer (Netherland: Springer Science+Business Media), 269–285. doi: 10.1007/978-1-4939-8621-7_108
Ruiz, K. B., Biondi, S., Oses, R., Acuña-rodríguez, I. S., Antognoni, F., Martinez-Mosqueira, E. A., et al. (2014). Quinoa biodiversity and sustainability for food security under climate change. A review. Agron. Sustain. Dev. 34, 349–359. doi: 10.1007/s13593-013-0195-0
Sawan, Z. M., Fahmy, A. H., and Yousef, S. E. (2011). Effect of potassium, zinc and phosphorus on seed yield, seed viability and seedling vigor of cotton (Gossypium barbadense L.). Arch. Agron. Soil Sci. 57, 75–90. doi: 10.1080/03650340903222328
Stikic, R., Glamoclija, D., Demin, M., Vucelic-Radovic, B., Jovanovic, Z., Milojkovic-Opsenica, D., et al. (2012). Agronomical and nutritional evaluation of quinoa seeds (Chenopodium quinoa Willd.) as an ingredient in bread formulations. J. Cereal Sci. 55, 132–138. doi: 10.1016/j.jcs.2011.10.010
Szydłowska-Czerniak, A., Bartkowiak-Broda, I., Karlović, I., Karlovits, G., and Szłyk, E. (2011). Antioxidant capacity, total phenolics, glucosinolates and colour parameters of rapeseed cultivars. Food Chem. 127, 556–563. doi: 10.1016/j.foodchem.2011.01.040
Tang, Y., Li, X., Chen, P. X., Zhang, B., Hernandez, M., Zhang, H., et al. (2015). Characterisation of fatty acid, carotenoid, tocopherol/tocotrienol compositions and antioxidant activities in seeds of three Chenopodium quinoa Willd. genotypes. Food Chem. 174, 502–508. doi: 10.1016/j.foodchem.2014.11.040
Tang, Y., and Tsao, R. (2017). Phytochemicals in quinoa and amaranth grains and their antioxidant, anti-inflammatory, and potential health beneficial effects: a review. Mol. Nutr. Food Res. 61:1600767. doi: 10.1002/mnfr.201600767
Thanapornpoonpong, S.-N., Vearasilp, S., Pawelzik, E., and Gorinstein, S. (2008). Influence of various nitrogen applications on protein and amino acid profiles of amaranth and quinoa. J. Agric. Food Chem. 56, 11464–11470. doi: 10.1021/jf802673x
Valenzuela-Bustamante, P. (2015). Evaluación de la actividad antioxidante y determinación del contenido de fenoles totales y flavonoides de hojas de diferentes genotipos de Ugni molinae Turcz. Ph D. Thesis, (Netherland: Springer Science+Business Media).
Vega-Gálvez, A., Miranda, M., Vergara, J., Uribe, E., Puente, L., and Martínez, E. A. (2010). Nutrition facts and functional potential of quinoa (Chenopodium quinoa Willd.), an ancient Andean grain: A review. J. Sci. Food Agric. 90, 2541–2547. doi: 10.1002/jsfa.4158
Vidueiros, S. M., Curti, R. N., Dyner, L. M., Binaghi, M. J., Peterson, G., Bertero, H. D., et al. (2015). Diversity and interrelationships in nutritional traits in cultivated quinoa (Chenopodium quinoa Willd.) from Northwest Argentina. J. Cereal Sci. 62, 87–93. doi: 10.1016/j.jcs.2015.01.001
Villanueva, A., Vioque, J., Sânchez-Vioque, R., Clémente, A., Pedroche, J., Bautista, J., et al. (1999). Peptide characteristics of sunflower protein hydrolysates. JAOCS 76, 1455–1460. doi: 10.1007/s11746-999-0184-2
Wimalasekera, R. (2015). “Role of Seed Quality in Improving Crop Yields,” in Crop Production and Global Environmental Issues, ed. K. R. Hakeem (Cham: Springer), 153–168. doi: 10.1007/978-3-319-23162-4
Keywords: quinoa, nutritional traits, seeds (grains), environmental adaptability, emerging crops, environment × genotype
Citation: Granado-Rodríguez S, Aparicio N, Matías J, Pérez-Romero LF, Maestro I, Gracés I, Pedroche JJ, Haros CM, Fernández-Garcia N, Navarro del Hierro J, Martin D, Bolaños L and Reguera M (2021) Studying the Impact of Different Field Environmental Conditions on Seed Quality of Quinoa: The Case of Three Different Years Changing Seed Nutritional Traits in Southern Europe. Front. Plant Sci. 12:649132. doi: 10.3389/fpls.2021.649132
Received: 03 January 2021; Accepted: 16 April 2021;
Published: 12 May 2021.
Edited by:Heng Zhang, Shanghai Center for Plant Stress Biology, Shanghai Institute for Biological Sciences, Chinese Academy of Sciences (CAS), China
Reviewed by:Jorge Antonio Rojas-Beltrán, University of San Simón, Bolivia
Hector Daniel Bertero, University of Buenos Aires, Argentina
Copyright © 2021 Granado-Rodríguez, Aparicio, Matías, Pérez-Romero, Maestro, Gracés, Pedroche, Haros, Fernandez-Garcia, Navarro del Hierro, Martin, Bolaños and Reguera. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: María Reguera, firstname.lastname@example.org