Structural Characterization and Antioxidant Capacity of Quinoa Cultivars Using Techniques of FT-MIR and UHPLC/ESI-Orbitrap MS Spectroscopy

The existence of more of 16,000 varieties of quinoa accessions around the world has caused a disregard on their structural and phytochemical characteristics. Most of such accessions belong to cultivars settled in Colombia. The goal of this research was to evaluate the structural attributes and antioxidant capacities from six quinoa cultivars with high productive potential from central regions in Colombia. This study used middle-range infrared spectroscopy (IR-MIR) to determine the proteins, starch and lipids distinctive to quinoa grains. Ultra-high-performance liquid chromatography electrospray ionization Orbitrap, along with high-resolution mass spectrometry (UHPLC/ESI-Orbitrap MS), were also used to identify the existence of polyphenols in cultivars. The antioxidant capacity was determined through DPPH, ABTS and FRAP. The spectrums exhibited significant variances on the transmittance bands associated with 2922 cm−1, 1016 cm−1 and 1633 cm−1. Moreover, the intensity variations on the peaks from the secondary protein structure were identified, mainly on the bands associated with β-Sheet-1 and -2, random coil α elice and β-turns-2 and -3. Changes found in the ratios 996 cm−1/1014 cm−1 and 1041 cm−1/1014 cm−1 were associated with the crystalline/amorphous affinity. Regarding the antioxidant capacity, great differences were identified (p < 0.001) mainly through FRAP methods, while the phenolic acids and flavonoids were determined by UHPLC/ESI-Orbitrap MS techniques. The presence of apigenin and pinocembrin on grains was reported for the first time. Titicaca and Nariño were the most phytochemically diverse quinoa seeds.


Introduction
According to different scientific reports, quinoa had its origins in South America between Colombia and Chile. However, its greater abundance in terms of morphologic, ecophysiology and nutritional diversity happens in countries like Ecuador, Bolivia, Perú and Chile. The distribution of cultivars happens according to their source of origin as quinoa from: the altiplano Inter-Andean Valley, from the coast and from Salares and Yungas [1]. Quinoa species show a great adaptability to harsh edaphoclimatic conditions, mainly related to dramatic changes in the temperature, salinity stress and lower availability of water and nutrients [2].

Quinoa Grain Characteristics
Color has recently become a classification factor for quinoa grains, as well as a quality attribute ( Figure 1). Colombian cultivars show a cream-white color associated with genomic quinoas from the inter-Andean Valley region, with typical short days with a high sensitivity to luminosity [8,19]. However, the genetic diversity and increment of its consumption has brought an ease to its spreading to new cultivars, showing attributes linked to pericarpic coloring, dimensions and constitutive molecules.
The evaluated cultivar grains showed significant differences between the L*a*b* coordinates linked to external integument coloration. Such a condition happened to all the cultivars except Titicaca and Puno (on L*) and Salcedo and Puno on a*. The proteins, carbohydrates and fats were studied on all the cultivars. Titicaca and Pasankalla had greater protein values; Puno and Nariño were rich in carbohydrates, while the fats were greater in Puno region cultivars (Table 1).

Quinoa Grain Characteristics
Color has recently become a classification factor for quinoa grains, as well as a quality attribute ( Figure 1). Colombian cultivars show a cream-white color associated with genomic quinoas from the inter-Andean Valley region, with typical short days with a high sensitivity to luminosity [8,19]. However, the genetic diversity and increment of its consumption has brought an ease to its spreading to new cultivars, showing attributes linked to pericarpic coloring, dimensions and constitutive molecules.

Titicaca
Salcedo Soracá Faro Pasankalla Nariño The evaluated cultivar grains showed significant differences between the L*a*b* coordinates linked to external integument coloration. Such a condition happened to all the cultivars except Titicaca and Puno (on L*) and Salcedo and Puno on a*. The proteins, carbohydrates and fats were studied on all the cultivars. Titicaca and Pasankalla had greater protein values; Puno and Nariño were rich in carbohydrates, while the fats were greater in Puno region cultivars (Table 1).

Infrared Spectroscopy
The results allowed to identify a 1633-cm −1 transmittance linked to C=O stretching, which is also related to its protein composition [20]. Table 2 shows meaningful differences of the transmittance values from the cultivars. Since there is a direct relationship between the transmittance and protein, it is possible to say that Pasankalla is the one with a lesser protein value in contrast to Nariño, Salcedo and Puno, with greater amounts.
This way, band differences associated with C-H stretching and C-O-C and C-O folding were observed. About starch, a major proportion was observed on the Titicaca and Puno cultivars, while, in Soracá and Pasankalla, the rates were lower. The stretching of the C-H path group linked to lipids made it possible to identify the major differences among the cultivars.  Table 3 shows the bands after performing a deconvolution on amide I in the region between 1600 and 1700 cm −1 . This is the most spectral-sensitive area to secondary protein structures but also responsible for the vibration and strain of the peptic C=O bond [21]. The amide I band spectrum allowed to determine the secondary structure protein for all six cultivar quinoa grains, as in the previously described parameter analysis protocols. Table 3. Secondary protein structures in six Colombian quinoa cultivars. There were noticeable peaks observed on 800 cm −1 and 1300 cm −1 attributed to starch and corresponding to a flexing COH vibration. According to the data, the evaluated cultivars showed distinguished results (p < 0.05) on each region, the most intense at 996 cm −1 and 1076 cm −1 . As for the 1041-cm −1 crystalline structure band, it had a greater intensity in Pasankalla and Nariño, while the 1014-cm −1 region had more shapeless structural intensity in Nariño (Table 4). Results are expressed as the mean ± standard deviation (n = 3). Different letters had significant differences (p < 0.05) using Tukey's test.

Polyphenols and Antioxidant Capacity by ABTS, DPPH and FRAP
Quinoa crops showed polyphenol values of 0.6681-1.7737-mg AG/g/sample, which were quite different among them (p < 0.001). A richer concentration was found in the Nariño cultivar. About the antioxidant capacity, according to the DPPH technique, Nariño had the highest value. Based on the FRAP technique, Titicaca and, on the ABTS method, both cultivars had the greatest values (Table 5). Results are expressed as the mean ± standard deviation (n = 3). Different letters had significant differences (p < 0.001) using Tukey's test. AG: gallic acid; T: Trolox; AA: ascorbic acid.

Phenolic Compounds under UHPLC-ESI-Orbitrap MS
The results showed phenolic acid, xanthine, flavonoids and anthocyanin varieties among farms. Contrasting values were observed (p < 0.001), highlighting the Nariño cultivar for its apigenin (0.08 ± 0.01) presence and the highest value (42%) of identified elements. Salcedo had an elevated value of caffeic and vanillic acid and Soracá of pcoumaric acid and pinocembrin (Table 6).

Phytochemical Multivariate Structure Analysis
The results showed organization in relation to the phenolic components. As noticed in the heatmap, the groupings showed a high diversity in bioactive compounds ( Figure 2A). The bootstrap exhibited two main groups: the first includes one xanthine and one anthocyanin and an antioxidant capacity method, while the second combined phenolic acid, a flavonoid, total polyphenols and an antioxidant technique ( Figure 2B).
An analysis of the quinoa cultivars' main components allowed to evaluate the differences in their bioactive composition ( Figure 3). The cumulative variance of the two first components was 67.2% (CP); however, when including a third component, the CP variance was 90.8%. It is possible to notice a discrimination of the first two CP from the variance of all six cultivars. The first group's Nariño cultivar showed the highest amounts of pinocembrin, p-coumaric, total polyphenol and high antioxidant capacity with DPPH. Secondly, Titicaca quinoa had greater values of quercetin-3-glucoside, caffeine and antioxidant capacity with ABTS and FRAP. Pasankalla and Soracá presented an observable similarity when adding PC3 as an average cultivar. Finally, Salcedo and Puno were similar because of the presence of some phenolic acids. Results are expressed as the mean ± standard deviation (n = 3). Different letters had significant differences (p < 0.001) using Tukey's test.

Phytochemical Multivariate Structure Analysis
The results showed organization in relation to the phenolic components. As noticed in the heatmap, the groupings showed a high diversity in bioactive compounds ( Figure  2A). The bootstrap exhibited two main groups: the first includes one xanthine and one anthocyanin and an antioxidant capacity method, while the second combined phenolic acid, a flavonoid, total polyphenols and an antioxidant technique ( Figure 2B An analysis of the quinoa cultivars' main components allowed to evaluate the differ ences in their bioactive composition ( Figure 3). The cumulative variance of the two firs components was 67.2% (CP); however, when including a third component, the CP vari ance was 90.8%. It is possible to notice a discrimination of the first two CP from the vari ance of all six cultivars. The first group's Nariño cultivar showed the highest amounts o pinocembrin, p-coumaric, total polyphenol and high antioxidant capacity with DPPH Secondly, Titicaca quinoa had greater values of quercetin-3-glucoside, caffeine and anti oxidant capacity with ABTS and FRAP. Pasankalla and Soracá presented an observable similarity when adding PC3 as an average cultivar. Finally, Salcedo and Puno were simi lar because of the presence of some phenolic acids.

Discussion
The visual aspects of quinoa seeds are determinant on its variety and even compositional differentiation [22]. Over recent years, research evaluating colorimetric coordinates allowed compound identification to boost quinoa cultivars with black, pink, purple and red seeds. Furthermore, it brings about the relationship between tannins, saponins and phytic acid compounds [23,24].
In this sense, protein, lipidic and starchy attributes on grains are fundamental to consider agro-industrial processing or even more, when aiming towards hyper-protein starch and fat enriched for highly functional quality food [14,25].
To elaborate a structural analysis with the goal to deeply learn about the secondary protein structures of each cultivar, spectrum deconvolution techniques were performed on 1633 cm −1 [26]. Deconvolution of the 1016 cm −1 band can relate to the crystallinity of a shapeless starch granule [27]. This turned out to be an important fact, because the functional flour does not only depend on the protein and starch concentrations but, also, on its organization and molecular conformation.
The intensity from the amide I band from quinoa flour resulted in 1632 cm −1 and 1645 cm −1 . This condition is relevant to crops' characterization and differentiation. Nariño and Pasankalla offer a band intensity at 1645 cm −1 , Puno and Soracá at 1643 cm −1 , Salcedo at 1632 cm −1 and Titicaca at 1633 cm −1 .

Discussion
The visual aspects of quinoa seeds are determinant on its variety and even compositional differentiation [22]. Over recent years, research evaluating colorimetric coordinates allowed compound identification to boost quinoa cultivars with black, pink, purple and red seeds. Furthermore, it brings about the relationship between tannins, saponins and phytic acid compounds [23,24].
In this sense, protein, lipidic and starchy attributes on grains are fundamental to consider agro-industrial processing or even more, when aiming towards hyper-protein, starch and fat enriched for highly functional quality food [14,25].
To elaborate a structural analysis with the goal to deeply learn about the secondary protein structures of each cultivar, spectrum deconvolution techniques were performed on 1633 cm −1 [26]. Deconvolution of the 1016 cm −1 band can relate to the crystallinity of a shapeless starch granule [27]. This turned out to be an important fact, because the functional flour does not only depend on the protein and starch concentrations but, also, on its organization and molecular conformation.
The intensity from the amide I band from quinoa flour resulted in 1632 cm −1 and 1645 cm −1 . This condition is relevant to crops' characterization and differentiation. Nariño and Pasankalla offer a band intensity at 1645 cm −1 , Puno and Soracá at 1643 cm −1 , Salcedo at 1632 cm −1 and Titicaca at 1633 cm −1 .
According to Yang et al. (2015) [26], the absorbance peak at the amide I band withstood a deconvolution process to determine the structures like β-sheet (1624-1642 cm −1 ), Random coil (1645-1648 cm −1 ), α-elice (1653-1656 cm −1 ) and β-turns (1667-1694 cm −1 ). All the cultivars displayed contrasting intensity peaks on the secondary structures ( Figure 4). Titicaca had a greater β-sheet-1 and -2 proportions, Soracá β-sheet-3, Random coil and α-elice. Titicaca had β-turns-1, while there was β-turns-2 and -3 in Nariño.  In this sense, the discovered results were dramatically different according to other research because of its major manifestation, the peak intensity was found using random coil and α-elice. The results obtained agree with what was stated by Wang, Zhao, and Yuan (2020) [28], who found secondary structures of proteins in the same peak position in the original spectra in this research. Wolkers et al. (1998) [29] explained that such variability could be associated with effects on its production regions because of the edaphoclimatic conditions.
It is also important to outline that secondary structure spectrums could be also used like a fingerprint distinct to each cultivar. Improving the spectroscopy methods could contribute to identify not only the genetic variability but to offer quality control towards adulterated flour [30].
Then, when looking at each cultivar protein fingerprint, it could be possible to assert that the structural differences directly affect the technical and functional capacity to produce gel, foams or emulsions [31].
Ratios 1041/1014 and 996/1014 cm −1 were used with greater frequency to measure the crystallinity range in starch. Such a fact made it possible to establish that the range 1041/1014 cm −1 varied among all six evaluated cultivars ( Figure 5). The organized structure of granules indicates that starch changes mostly because of the two conditions determining photoassimilate transports from organs to grains. The first was related to K + movement on apoplast, determining the flow pressure. The second, the invertase acid activity influencing the transforming saccharose in hexoses and restructuring dispose organs for the nutritional quality of seeds [32,33].
Such results are quite important, because to obtain gel or translucent films, is important to use starch with a high crystallinity index, to which the Soracá cultivar is the most suitable (1041/1014 = 4.02). Retrograde phenomenon must be avoided when using flour as the complementing nutritional value. Then, starch with a lower crystallinity index should be used, in which 1.3 and 1.48 indexes from Pasankalla and Nariño are the most appropriate. In this sense, the discovered results were dramatically different according to other research because of its major manifestation, the peak intensity was found using random coil and α-elice. The results obtained agree with what was stated by Wang, Zhao, and Yuan (2020) [28], who found secondary structures of proteins in the same peak position in the original spectra in this research. Wolkers et al. (1998) [29] explained that such variability could be associated with effects on its production regions because of the edaphoclimatic conditions. It is also important to outline that secondary structure spectrums could be also used like a fingerprint distinct to each cultivar. Improving the spectroscopy methods could contribute to identify not only the genetic variability but to offer quality control towards adulterated flour [30].
Then, when looking at each cultivar protein fingerprint, it could be possible to assert that the structural differences directly affect the technical and functional capacity to produce gel, foams or emulsions [31].
Ratios 1041/1014 and 996/1014 cm −1 were used with greater frequency to measure the crystallinity range in starch. Such a fact made it possible to establish that the range 1041/1014 cm −1 varied among all six evaluated cultivars ( Figure 5). The organized structure of granules indicates that starch changes mostly because of the two conditions determining photoassimilate transports from organs to grains. The first was related to K + movement on apoplast, determining the flow pressure. The second, the invertase acid activity influencing the transforming saccharose in hexoses and restructuring dispose organs for the nutritional quality of seeds [32,33].
Such results are quite important, because to obtain gel or translucent films, is important to use starch with a high crystallinity index, to which the Soracá cultivar is the most suitable (1041/1014 = 4.02). Retrograde phenomenon must be avoided when using flour as the complementing nutritional value. Then, starch with a lower crystallinity index should be used, in which 1.3 and 1.48 indexes from Pasankalla and Nariño are the most appropriate. In this sense, the high diversity of quinoa makes it possible to find variations in the foliar area of the plants, which is related to the CO2 net assimilation and to obtain triose phosphate in the CO2 fixation phase, which allows the obtaining of starches for their translocation to sink organs [34,35]. Additionally, Reguera and coworkers found changes in the synthesis of secondary metabolites of three quinoa cultivars and highlighted that these compounds seem to be genotype-specific and could vary significantly under stressful conditions [9] The polyphenol values vary in all quinoa cultivars and are linked to climatic conditions, agronomical practices and genetic interaction [17]. Li, Lietz & Seal (2021) [36] reported 1.44-mg AG/g of total phenols and antioxidant activity between 8.61 and 522 µ m of Trolox/g under ABTS and DPPH. Such results are different from the findings in this paper. The Pasankalla and Salcedo farms reported an antioxidant activity of 5.6-and 3.18µ m Trolox/g under DPPH [37] and a 27.9-µ m Trolox/g value with the ABTS technique on a Bolivian farm [38]. Such results are meaningfully high compared to the ones observed here.
This parameter turns to be important to perform a study on how to enrich nutritional models with a low antioxidant capacity, as well as to identify cultivars that might bring about this quality [39]. Likewise, the cultivars with the highest antioxidant capacity were Nariño and Titicaca, because their measurements were highly recognized over other samples during at least two different techniques.
Polyphenols have been recently studied in quinoa, because their diversity is equal or greater than other cereals. Findings of 0.14-0.33 mg/kg of caffeic acid and 0.4-and 0.94mg/kg vanillic acid were in ten different accessions in Chilean quinoa [18].
The results observed by the heatmap and bootstrap determined the bioactive differences from each evaluated cultivar. It is relevant to mention that many polyphenols are not synthetized on a regular basis in plants, but this could be a result from external conditions such as biotic and abiotic factors [41]. In this sense, the high diversity of quinoa makes it possible to find variations in the foliar area of the plants, which is related to the CO 2 net assimilation and to obtain triose phosphate in the CO 2 fixation phase, which allows the obtaining of starches for their translocation to sink organs [34,35]. Additionally, Reguera and coworkers found changes in the synthesis of secondary metabolites of three quinoa cultivars and highlighted that these compounds seem to be genotype-specific and could vary significantly under stressful conditions [9].
The polyphenol values vary in all quinoa cultivars and are linked to climatic conditions, agronomical practices and genetic interaction [17]. Li, Lietz & Seal (2021) [36] reported 1.44-mg AG/g of total phenols and antioxidant activity between 8.61 and 522 µm of Trolox/g under ABTS and DPPH. Such results are different from the findings in this paper. The Pasankalla and Salcedo farms reported an antioxidant activity of 5.6-and 3.18-µm Trolox/g under DPPH [37] and a 27.9-µm Trolox/g value with the ABTS technique on a Bolivian farm [38]. Such results are meaningfully high compared to the ones observed here.
This parameter turns to be important to perform a study on how to enrich nutritional models with a low antioxidant capacity, as well as to identify cultivars that might bring about this quality [39]. Likewise, the cultivars with the highest antioxidant capacity were Nariño and Titicaca, because their measurements were highly recognized over other samples during at least two different techniques.
Polyphenols have been recently studied in quinoa, because their diversity is equal or greater than other cereals. Findings of 0.14-0.33 mg/kg of caffeic acid and 0.4-and 0.94-mg/kg vanillic acid were in ten different accessions in Chilean quinoa [18].
The results observed by the heatmap and bootstrap determined the bioactive differences from each evaluated cultivar. It is relevant to mention that many polyphenols are not synthetized on a regular basis in plants, but this could be a result from external conditions such as biotic and abiotic factors [41].
The similarities between quercetin-3-glucoside, FRAP and ABTS happen because of the high capacity from anthocyanins to capture free radicals. In the same way, the concentration resemblance in some phenolic acids occurs because of methylation from a few compounds like caffeic acid, the forefather of vanillic acid [42]. The correlation between ferulic acid and p-coumaric acid is a consequence, because they both belong to the hydroxycinnamic group.
It is accurate to mention that this paper reports, for the first time, the presence of pinocembrin and apigenin in quinoa compounds. For other species, its synthesis is attributed to pathogen agents affecting the plant during primary production. Additionally, the influence of edaphoclimatic conditions like salinity, water deficit or very cold weather in amaranth plants was found [43,44].
High contributions on the PCA of ferulic acid and vanillic acid, among other acids, coincide with the ones reported by Antognoni et al. (2021) [45]. They identified that Titicaca has the greatest polyphenol values expressed by these two acids and that its concentration could be greater than 100 mg/kg in quinoa flour.
Carciochi, Manrique & Dimitrov (2014) [46] found a relationship between ferulic acid and vanillic acid concentrations in seeds. Outlining that the two elements are the ones supporting the seed antioxidant activity, in parallel, its content value rises along with the plant growing towards germination.

Materials
We used Salcedo, Puno, Pasankalla, Soracá, Titicaca and Nariño cultivars collected from crops established in the municipalities of Oicatá, Ventaquemada and Tuta in the Department of Boyacá, and others were taken from the seed collection of the "Laboratorio de Biotecnología Vegetal de la Gobernación de Boyacá". The cultivars were grown in a greenhouse located at Victoria Granja Agroecologica Company, Ventaqumada, Colombia During the study, the irrigation and preventive sanitary management were carried out in order to maintain plant health. The seeds were manually cleaned to remove foreign matter and broken and immature seeds.
The total polyphenol reagents were Folin Ciocalteau, sodium acid carbonate (99%) and gallic acid. The

Proximal Characteristics
To determinate the protein content, the Kjeldahl method was used. AOAC 960.52 was multiplied by a conversion factor of 6.25 to measure the total raw protein from the seeds. Fats were defined by the AOAC 922.06 method using Soxhlet (Soxtec 2050). Finally, the total carbohydrates were calculated by the difference (i.e., protein + ash + fat + moisture − 100).

CIEL*a*b* Coordinates Determination
The spectrophotometer used was Konica Minolta (Model CM-5; Tokyo, Japan) to determine the seed color hue from the quinoa cultivars. The results were recorded as L* a* b* values.

Mid-Range Infrared Spectroscopy
IR spectrum was obtained using IRAFFINITY-1S equipment (Shimadzu Corp., Kyoto, Japan). The spectrum was performed by the reflection mode between 400 cm −1 and 4000 cm −1 and proportional over 32 scans with 4 cm −1 and a temperature of 25 • C. In the previous measurements, each cultivar had a blank background spectrum.
A spectrum analysis was made by using OriginPro, 7th version. Measurements were fixed and standardized over a baseline between 0 and 1 (represented in the figures) to show the highest transmittance peaks. Deconvolution was performed on a spectrum using Fourier transform to determine the protein secondary structure changes (band 1600-1700 cm −1 ). Later, a modeling of the Gaussian function spectrum was applied. The starch concentration was studied from its formation on a short-range band between 875 cm −1 and 1175 cm −1 . Lipids were present on the 2800-2900 cm −1 range. All samples were analyzed in triplicate.

Polyphenols
From each cultivar, 0.5 g of quinoa flour were mixed. Then, to each sample was added 30 mL of sodium hydroxide 2M and stirring for 10 s with vortex (Fisher Scientific, G560, Bohemia, Waltham, MA, USA). It was left for extraction over 24 h on amber glass at room temperature (25 • C) and mixed in darkness. After that, 5.7 mL of hydrochloride acid were added to the samples upon reaching a pH value of 2. Finally, the samples were centrifugated for 15 min with a rotational speed of 10,000 rpm (HERMLE model Z306, La-Bortechnik GmbH, Wasserburg, Germany), then vacuum-filtrated with Whatman paper (No 1).
Phenolic compounds were determined through the Folin Ciocalteau (FC) method. Triple samples were disposed with 40 µL on a separate 1800 µL with FC reagent. For 15 min, a vortex mix was performed to later rest for 5 min. Afterwards, 1200 µL of sodium acid carbonate (NaHCO 3 ) at 7.5% (m/v) was added and let to settle for 60 min. Absorbance at 765 nm was measured on a V-630 UV-VIS model spectrophotometer (JASCO Inc., Easton, MD, USA). The phenol amounts were measured by sample grams expressed in gallic acid milligrams (mg AG/g).

FRAP
In a 50-mL falcon tube, 0.4 g of the sample were placed. Then, 10 mL of ethanol were added to begin the extraction process. This was mixed by vortex for 15 s and left to later extract the liquid-solid in the mixer over 30 min at 37 • C. After, the samples were centrifugated (15 min. at 5000 rpm and 4 • C) and later vacuum-filtered with Whatman paper (No 1).
Test tubes with a lid were prepared for 3 different samples: blank, experiment and calibration. For neutral (blank) reagents, the sample had 60 µL of methanol. To begin the reaction, 180 µL of distill water were put in each tube. Then, there was 15 s of vortexing to later add 1800 µL of the FRAP reagent. Next, it was incubated at 37 • C for 30 min.
Finally, 595-nm absorbance was measured on a V-630 UV-VIS model spectrophotometer (JASCO Inc., Easton, MD, USA). The antioxidant activity was expressed with µmol of ascorbic acid/g of the sample (µmol de AA/g).

ABTS
In a 50-mL falcon test tube, 0.4 g from the sample was disposed. To begin the extraction, 10 mL of ethanol were added, vortexed for 15 s and left for solid-liquid extraction over 16 h at 4 • C. Subsequently, the samples were centrifugated over 15 min at 10,000 rpm, then vacuum-filtered with Whatman paper (No 1).
In the test tube, 4 mL of ABTS solution was placed and totally covered with aluminum foil. To begin the process, 135 µL of standard solution were added and then mixed by vortexing for 5 s. Blank reagent consisted of 4 mL of buffer acetate and 135 µL of ethanol. The zero point was mixed with 4.5-mL ABTS solution and 135 µL of ethanol. The test tube was closed and left to react for30 min to finally measure the absorbance at a 729.7-nm wavelength.

DPPH
In a 50-mL falcon test tube, 0.4 g were measured before adding 10 mL of ethanol. Vortexing for 15 s was done and then, the solid-liquid extract was left to settle for 16 h at 4 • C. After such a time, the samples were centrifugated for 15 min at 10,000 rpm, then vacuum-filtered with Whatman paper (No 1).
In a test tube, 3.9 mL of DPPH solution and 100 µL of standard solution were applied to launch a reaction by vortex stirring for 5 s. Blank reagents (control) were carried out with ethanol. The zero point was adjusted by adding 3.9 mL of DPPH solution and 100-µL ethanol. The sample was covered for 30 min to initiate the reaction to later measure the absorbance at 517 nm.

Phenolic Compounds Determination by UHPLC/ESI-Orbitrap MS
The seeds were macerated for four minutes with a blade grinder (Krups; Solingen, Germany). Later, 100 mL of 80% ethanol at 60 • C was added to 5.0 g of quinoa flour. The resultant extract was centrifugated at 10,000 rpm over 10 min. The leftover sample was recovered and stored at 4 • C for later use.
The extract was analyzed on a liquid high-efficiency chromatograph (UHPLC) Dionex Ultimate 3000 (Thermo Scientific, Sunnyvale, CA, USA, EE.UU.). This equipment consists of a binary gradient pump (HP G3400RS) and automatic sample injector (WPS 300TRS) and a thermostatic unit for a TCC 3000 column. The interface of LC-MS was electrospray ionization (ESI). The spectrometer Orbitrap had high mass resolution with ion current detecting system. Chromatography separation was performed by Hypersil GOLD Aq 100 × 2.1 mm, 1.9 µm (Thermo Scientific, Sunnyvale, CA, USA, EE.UU.) column at 30 • C.
Mobile phase A consisted of a 0.2% aqueous solution of ammonium formate. Mobile phase B: 0.2% on acetonitrile ammonia formate. The initial gradient condition was 100%. A moved linearly towards 100% of B in 8 min. Such a phenomena was held for 4 min, then returned to the initial conditions in one minute. The total running time lasted 13 min, plus 3 post-running additional minutes.
All compound identification was achieved by using the full scan mode and ion extraction (EIC) corresponding to the [M+H] + of the target compounds. A mass assessment was performed with the accuracy precision of ∆ ppm < 0.001 and using a standard mix solution of phenolic compounds. To quantify the samples, the external standard calibration was applied using the results factor (Rf) stablished by the standard solutions analysis under different concentrations.

Statistical Analysis
The values presented are the means ± standard deviation of three replicates. The data were analyzed using the Shapiro-Wilk normality test and Bartlett's homogeneity test, and the variance analysis (ANOVA) by Tukey's test was done to observe significant differences (p < 0.05). A principal component analysis (PCA) was collected to stablish the relationship of the antioxidant and polyphenol changes among the cultivars using the Factoextra and FactoMinerR libraries. A bootstrap analysis was elaborated using the pvclust function [47]. A double-grouping analysis was performed by Manhattan distance with the gplots and RColorBrewer libraries [48]. All the analyses were conducted using R version 3.6.1 software.

Conclusions
A spectroscopy analysis of different quinoa cultivars showed significant variations attributed to their genetic character. This showed that not only external factors but, also, morphological and physiological diversity are influential on the structural configuration of proteins. The same conditions happen to the crystalline/amorphous properties of starch granules.
It is relevant to mention that some alloys like caffeic acid, quercetin 3-glucosid, caffeine and pinocembrin are not very well-identified in quinoa grains; therefore, the influence over its antioxidant activity is still uncertain.
The teamwork between the spectrum deconvolution technique, chromatographic analysis and multivariate analysis are strong tools used for the seed characterization of cultivars from the Colombian Inter-Andean Valley.