Avaliação nutricional e aceitação pelo consumidor de extrato de soja de cultivares especiais e convencionais

bjft Brazilian Journal of Food Technology Braz. J. Food Technol. 1981-6723 Instituto de Tecnologia de Alimentos - ITAL Resumo O extrato de soja é um produto pronto para o consumo, nutritivo e naturalmente livre de colesterol. Diversas cultivares especiais de soja foram desenvolvidas no Brasil com o objetivo de aumentar o consumo humano de produtos derivados da soja. O objetivo deste trabalho foi avaliar a composição nutricional e a aceitação pelo consumidor do extrato de soja de três cultivares de soja com características especiais, como alto teor de proteínas e isoflavonas, sabor suave e livre de lipoxigenases, em relação a três cultivares convencionais. O extrato de soja foi obtido após branqueamento com solução de bicarbonato, trituração, centrifugação e pasteurização. Os resultados foram analisados por meio do teste de Anova e Tukey, para verificar diferenças entre as médias, e os dados do consumidor, por meio de análise de agrupamento e mapa de preferência interno. As lipoxigenases foram inativadas após o processamento térmico. Houve diferenças significativas entre as cultivares de soja e os extratos obtidos em relação aos teores de proteína, isoflavonas, lipídios e açúcares, e para o rendimento de extrato (p < 0,05). As cultivares BRS 133 e BRS 284 apresentaram o maior teor de isoflavonas totais (p < 0,05). Em relação à aceitação geral pelo consumidor, houve diferenças significativas entre as cultivares (p < 0,05), sendo que a menor média foi obtida por BRS 267. No entanto, a análise de cluster identificou três segmentos de consumidores de acordo com a similaridade de preferência. As cultivares BRS 284 e BRS 267 alcançaram a maior média entre os consumidores do segmento 1 (6,8 e 6,7, respectivamente). Os consumidores do segmento 2 preferiram a cultivar sem lipoxigenases, BRS 213 (média de 7,5), enquanto o segmento 3 preferiu a BRS 133 (média de 7,4). O processo de branqueamento da soja com solução de bicarbonato, trituração com água quente e pasteurização atenuou as diferenças entre as cultivares convencionais e especiais, e as seis cultivares avaliadas foram adequadas para a produção de extrato de soja. 1 Introduction Efforts to increase soybean consumption are significant as scientific studies continue to demonstrate that soybean intake contributes to maintaining health and, possibly, is related to the reduction of chronic diseases such as breast and prostate cancer, osteoporosis, and heart diseases (Messina, 2016; Nachvak et al., 2019; Li et al., 2020; Sansai et al., 2020; Dietary Guidelines for Americans, 2020; Cai et al., 2021). The soymilk stands out among soybean products as a ready-to-eat, nutritious, and naturally cholesterol-free product. However, sensory barriers must yet be overcome. In the traditional processing to obtain the soymilk, the cold-water grain soaking step favors the development of a beany flavor due to reactions catalyzed by the lipoxygenase enzymes, which result in different compounds with undesirable sensory characteristics (Nelson et al., 1976; Han et al., 2021). Such compounds limited the soymilk acceptance among western consumers for many years (Al Mahfuz et al., 2004). Lipoxygenases (LOX) are present in soybean as three isozymes: LOX1, LOX2, and LOX3 with different characteristics (optimal pH and formation of compounds) (Kumar et al., 2003). In that regard, technological alternatives have been employed for soybean processing resulting in products with a more pleasant aroma and flavor (Nelson et al., 1976; Felberg et al., 2009). The research focused on specialty soybean cultivars for human consumption in Brazil has been developed by the breeding program at Embrapa (Empresa Brasileira de Pesquisa Agropecuária – in English Brazilian Agricultural Research Corporation) (Carrão-Panizzi et al., 2001). The cultivars BRS 213 and BRS 257 are genetically free of the lipoxygenase enzymes (Oliveira et al., 2010). Specialty cultivars also include BRS 216 which shows a high protein content (Silva et al., 2009). The cultivar BRS 267 is not a lipoxygenase free but has a pleasant flavor, possibly due to its particular sugar profile (Oliveira et al., 2010). Some cultivars such as BRS 232 are considered conventional and have no specific characteristics for human consumption, but show a light hilum color. Some specialty and conventional soybean cultivars were evaluated regarding the nutritional composition (proteins, fatty acids, sugars) and presence of bioactive compounds such as isoflavones in which Silva et al. (2009) and Ciabotti et al. (2019) showed the differences among them. However, consumer acceptance of soymilk was not performed. On the other hand, the acceptability of the beverage obtained with powdered soymilk from BRS 213 (lipoxygenase free) was higher than the commercial products for consumers of different cities in Brazil (Silva et al., 2007), showing the relevance of consumer acceptance evaluation. There are few results regarding the sensory evaluation of soymilk using these specialty cultivars in Brazil. This study aimed to investigate the performance of different soybean cultivars, specialty and conventional to produce soymilk regarding the soymilk yield, nutritional characteristics, presence of isoflavones and consumer acceptance. 2 Material and methods 2.1 Material The soybean cultivars were selected based on previous studies that investigated their nutritional characteristics and absence of LOX. The cultivars were provided by Embrapa Soybean (Londrina, PR, Brazil), and their characteristics are described in Table 1. Table 1 Soybean cultivars from the genetic breeding program at Embrapa Soybean. Cultivars Characteristics BRS 216* Specialty type for soybean sprouts, small grains, high protein and isoflavone contents BRS 267* Specialty type with pleasant flavor, large grains BRS 133* Conventional type with high isoflavone contents BRS 213* Specialty type with absence of lipoxygenases (LOX1, LOX2, LOX3) BRS 232 Conventional type and mild flavor BRS 284 Conventional type and mild flavor * Silva et al. (2009). 2.2 Processing Soymilk processing was based on Felberg et al. (2009) and Antoniassi et al. (2008) and consisted of dehulling, cooking for 10 minutes in sodium bicarbonate solution (0.25%) at a ratio of 1:3 (dehulled soybean: solution), draining, grinding with boiling water using a Waring® heavy-duty blender at a ratio of 1:8 (blanched soybeans: water) for two minutes, and centrifugation using an IEC® K7165 centrifuge with a 150 μm nylon filter at 4000 rpm (~2500 g). After the centrifugation step, the soymilk was pasteurized at 95-98 °C for 10 minutes. The process was carried out five times for each cultivar. The soymilk yield was measured as liters of soymilk resulting from 1 kg of dehulled soybean grains. The soymilk had about 7% of total solids and was formulated by adding 3% sugar and 0.2% salt (w/w). 2.3 Physicochemical analysis of soybean and soymilk 2.3.1 Proximate analysis The proximate analysis of soybean and soymilk was performed in six replicates according to the Official Methods of Analysis of AOAC International (Association of Official Analytical Chemists, 2010). Moisture analysis was performed by oven-drying to constant weight at 105 °C (AOAC 925.45B). The ash content was measured at 550 °C in a muffle furnace (AOAC 923.03/32.1.05). The protein content was calculated based on the nitrogen content obtained by the Kjeldahl method (AOAC 4.2.11) using a 6.25 factor. Oil extraction was performed by acid hydrolysis followed by ethyl ether and petroleum ether extraction according to AOAC method 922.06. 2.3.2 Lipoxygenase analysis Lipoxygenase activity in soymilk (after pasteurization) was determined by spectrophotometry according to Axelrod et al. (1981) and Kumar et al. (2003), using six replicates. The soymilk was lyophilized (LIOBRAS K120) and defatted with petroleum ether (30 ºC to 60 ºC) in a Soxhlet apparatus for 16 h before extraction. Enzyme extraction was performed by stirring for one hour at 0 ºC to 4 ºC with 0.2M phosphate buffer (pH 6.8), followed by centrifugation at 10,000 rpm for 10 minutes (4 ºC) (Thermo Scientific Sorvall Legend XTR). LOX1 activity was measured at pH 9 (borate buffer), whereas LOX2+3 activity was measured at pH 6.8 using sodium linoleate as substrate. The absorbance of LOX1 and LOX2+3 was read at 234 and 238 nm, respectively, and was measured for 3 min at 30-second intervals using an Agilent 8453 UV-visible spectrophotometer. The slope of the linear part of the curve corresponded to the enzymatic activity (Ludikhuyze et al., 1998; Wang et al., 2008). 2.3.3 Chromatographic analysis Sugar analysis was performed with six replicates according to the method proposed by Macrae (1998), based on the chromatographic separation of the sample into an amino column using a High Performance Liquid Chromatography (HPLC) system (Waters™, Waltham, MA, USA), with quantification by external standardization. The chromatographic conditions used in the determinations were a Zorbax carbohydrate column (4.6 x 250 mm; 5mm – Agilent™) at 30 °C, isocratic elution mode with acetonitrile: water (75:25, v/v), refractive index detector mode W2410 (Waters™) at an internal temperature of 45 ºC, and a 20 µL injection volume. Sample extraction was conducted with 1 g of sample weighed in a 25 mL volumetric flask containing 10 mL of ultrapure water, followed by an ultrasonic bath for 10 min. Subsequently, 5 mL of acetonitrile was added, and the volume was completed with ultrapure water. Finally, the extract obtained after this process was filtered directly into the autosampler vial. The isoflavone analysis of soybean and soymilk was carried out according to AOAC method 2001.10 (Association of Official Analytical Chemists, 2010) using six and twelve replicates, respectively. The extractions were performed with methanol/water (80:20, v/v), followed by hydrolysis with NaOH solution, acidification, and filtration with Whatman filter paper No. 1. The analyses were carried out by HPLC-Photodiode Array detector (Alliance™ 2695 and 2996, Waters, Waltham, MA, USA) using a YMC-Pack Pro C18 column (5 μm, 4.6 mm × 250 mm, YMC, Kyoto, Japan) and a gradient with acetic acid solution and methanol running at 1.3 mL/min. The data were acquired at 260 nm. The identification and quantification of isoflavones were performed by comparing the peak retention time under investigation with those of the respective standards injected as a pool. The peak identities were confirmed by UV spectra to avoid coelutions. The conversion of the isoflavone concentrations (genistin, glycitin, and daidzin) into aglycon equivalents was calculated by multiplying the mass of each isoflavone form by the ratio of its aglycone molecular weight to the molecular weight of the individual form (Song et al., 1998). The total isoflavones as aglycon equivalents were determined by summing the concentrations of daidzein, glycitein, and genistein to the aglycon equivalent concentrations of daidzin, glycitin, and genistin. 2.4 Evaluation of soymilk consumer acceptance The soymilk from the specialty and conventional cultivars (Table 1) were evaluated for the consumer. In total, 102 consumers (64 women and 38 men) aged 18 to 66 years took part in the study. The subjects were recruited at a supermarket in the West Zone of Rio de Janeiro, in the state of Rio de Janeiro (RJ). The samples were presented individually in white 50 mL disposable plastic cups codified with three-digit numbers and served at 8 ± 2 ºC. The samples were evaluated using a 9-point hedonic scale (1: dislike it extremely and 9: like it extremely). The presentation order was balanced according to MacFie et al. (1989). The data were analyzed by Analysis of Variance (ANOVA) (p ≤ 0.05), Tukey’s test to assess differences between means, Internal Preference Mapping, and cluster analysis. In addition to the acceptance, participants also recorded what they liked most and least about each product. The study was submitted to and approved by the Research Ethics Committee of the Federal University of Rio de Janeiro (UNIRIO) - (TTDD:232/2011). 2.5 Statistical analysis The one-way ANOVA was performed using SAS/Stat 9.2 (Statistical Analysis System Institute, 2008) and Base SAS 9.2 (Statistical Analysis System Institute, 2009). The mean values of each parameter were further compared by Tukey’s test to check the difference among means. The significance level for all tests was set at 5%. The soybean and soymilk results were compared on a dry and wet basis. However, the comparison shown in the tables was performed on a wet basis. 3 Results and discussion 3.1 Physicochemical analysis of soybean and soymilk The cultivars showed significant differences with regard to the proximate composition and sugar content (p < 0.05) (Table 2). The cultivar BRS 216 showed the highest protein content (46.8 g 100 g-1) and the lowest oil content (22.2 g 100 g-1). On the other hand, BRS 284 showed the lowest protein content (38.3 g 100 g-1) and the highest lipid content (28.0 g 100 g-1). These results confirmed the inverse relationship between the protein and oil contents in soybean (Carrão-Panizzi et al., 2021). The ash content ranged from 5.2 to 6.4 g 100 g-1, and the highest and lowest ash contents found in BRS 267 and BRS 216 were attributed to large and small grains, respectively. The levels of protein, oil and ash contents were higher than those reported by Silva et al. (2009), Ciabotti et al. (2019) and Carrão-Panizzi et al. (2021). The sucrose content ranged from 3.1 to 5.3 g 100 g-1 and was higher in BRS 216 and BRS 133 (p < 0.05), whereas the raffinose content ranged from 0.9 to 1.8 g 100 g-1, with BRS 267 showing the highest values (p < 0.05). No significant differences were observed among cultivars for stachyose (2.9 to 4.8 g 100 g-1), except for the lower content observed in BRS 267 (p < 0.05). Silva et al. (2009) evaluated five Brazilian soybean cultivars as well as BRS 213, BRS 216, BRS 267 and BRS 133 and the ranges observed for sucrose, raffinose and stachyose were 3.4 to 4.3; 0.4 to 1.0 and 2 to 3.5 g 100 g-1, respectively. However, the results of this work were higher for raffinose and stachyose than those reported by Silva et al. (2009) while the sucrose content was higher or lower depending on the cultivar evaluated. In another study, Oliveira et al. (2010) evaluated 28 Brazilian and foreign genotypes/cultivars and observed very wide ranges of values for sucrose (2.4 to 5.9 g 100 g-1), stachyose (from 2.7 to 4.4 g 100 g-1) and raffinose (from 0.4 to 1.2 g 100 g-1). These features pointed out the influence of edaphoclimatic conditions and different years on the composition of soybeans. Table 2 Proximate composition and sugar content of six soybean varieties (g 100 g-1) (wet basis)*. Cultivar Chemical composition Sugars Protein Oil Ash Moisture Sucrose Raffinose Stachyose BRS 216 46.8a 22.2d 5.2c 9.2cd 5.28a 1.33bc 4.54a BRS 267 44.5b 23.0c 6.4a 10.0b 3.70b 1.80a 2.91b BRS 133 42.4c 24.0c 5.6b 11.0a 5.33a 0.90d 4.24a BRS 213 41.8c 24.3b 5.8b 9.4cd 3.82b 1.49b 4.47a BRS 232 41.6c 24.4b 5.9b 10.1b 3.42b 0.92d 4.37a BRS 284 38.3d 28.0a 5.7b 9.0d 3.11b 1.22c 4.79a * Results expressed as the average of six replicates. Means with different lowercase letters in the same column are significantly different (P < 0.05) by Tukey’s test. The soymilk samples showed significant differences with regard to nutritional characteristics (p < 0.05) (Table 3). The total solids content of the soymilk samples differed significantly, and the statistical evaluation was performed on a wet and dry basis, with similar results. The results in the tables are shown on a wet basis. The soymilk from the cultivar BRS 232 showed the highest protein content (3.37 g 100 g-1) and the lowest lipid content (2.05 g 100 g-1) (p < 0.05). However, its grain protein content was one of the lowest. BRS 284 soymilk showed the lowest protein content (2.85 g 100g-1) and the highest lipid content (3.15 g 100 g-1). On the other hand, the extracts from cultivars BRS 133, 232, and 267 showed no significant differences with regard to the protein content. It should be noted that cultivar BRS 267 showed the highest ash contents agreeing with the result obtained for the soymilk. There was a significant difference in the yield of soymilk (p < 0.05) and BRS 133 showed the lowest figure (2.8 L/kg). Although no significant difference was observed for the others cultivars evaluated, there is a trend toward higher yield for the conventional cultivars BRS 232 and BRS 284, indicating a more favorable extraction of proteins from soybean grains. This characteristic was related to grain softening during blanching, enhancing the grinding efficiency. Table 3 Proximate composition (g 100 g-1)* and yield** of soymilk obtained from six soybean cultivars (wet basis). Cultivar Protein Oil Ash Moisture Soymilk yield ** BRS 216 3.20b 2.26b 0.66b 89.1bc 3.1ab BRS 267 3.13bc 2.20bc 0.73a 89.7a 3.2ab BRS 133 3.11c 2.94a 0.55c 88.8cd 2.8b BRS 213 3.11c 2.41b 0.57c 88.3d 3.3ab BRS 232 3.37a 2.05c 0.63b 89.9a 3.7a BRS 284 2.85d 3.15a 0.52d 89.3b 3.45a * Results expressed as the average of six replicates. **Average soymilk yield obtained from five processes (Liters of soymilk without formulation/kg of dehulled soybeans). Means with different lowercase letters in the same column are significantly different (p < 0.05) by Tukey’s test. The soymilk samples showed no residual LOX1 activity. However, the cultivar BRS 284 showed residual LOX2+3 activity, indicating that the blanching conditions (boiling sodium bicarbonate 0.25% solution for 10 minutes), grinding with boiling water and pasteurization inactivated the lipoxygenases in most soybean cultivars. These blanching and grinding conditions were similar to those reported by Felberg et al. (2004, 2009). The effect of the blanching conditions on lipoxygenase activity was evaluated by Antoniassi et al. (2008), with different bicarbonate solution concentrations (0.25 to 0.5%) and 10 to 12 minutes of blanching resulting in the highest lipoxygenase inactivation. The total average content of isoflavones of the samples evaluated expressed as aglycone equivalents (wet basis) are shown in Table 4. The isoflavone contents showed significant differences among the soybean and soymilk samples evaluated (p < 0.05). The main isoflavones were genistin, daidzin, and genistein, followed by lower amounts of daidzein. Glycitein was not detected in any of the soymilk and soybean samples evaluated. In contrast, glycitin was quantified in the soybean samples of cultivars BRS 267, BRS 133, and BRS 284 (up to 6 mg 100 g-1), and in the soymilk of BRS 133 (up to 0.19 mg 100 g-1)(data not shown). Cultivar BRS 133 showed the highest total isoflavone content (183 mg 100 g-1) agreeing with Silva et al. (2009). On the other hand, BRS267 showed the lowest content of this component (44.70 mg 100 g-1). The USDA isoflavone database (Bhagwat et al., 2015) showed significant differences in the isoflavone content of soybean samples around the world. Mature and raw soybeans (from different countries and cultivars) showed daidzein contents ranging from 2.64 to 191.43 mg.100 g-1, genistein contents ranging from 5.56 to 276.21 mg 100 g-1, and glycitein contents ranging from 0 to 121.69 mg.100 g-1. The total isoflavones ranged from 10.04 to 440.72 mg 100 g-1. The isoflavone contents in soybean vary due to genetic differences, geographic location, soil, year, and environmental conditions during growth (Wang & Murphy, 1994; Carrão-Panizzi et al., 2009). The isoflavone profile obtained was within the range observed for other soybean cultivars grown in Brazil and around the world. Table 4 Isoflavone content (mg 100 g-1, wet basis) of soybean and soymilk from different Brazilian soybean cultivars. Cultivars Soybean Soymilk Daidzein Genistein Daidzin* Genistin* Total isoflavone* Daidzein Genistein Daidzin* Genistin* Total isoflavone* BRS 216 1.78d 4.20e 33.96c 65.44d 105.38d 0.13d 0.22d 2.08c 5.64e 8.07d BRS 267 2.14c 7.27c 8.67e 25.30e 44.69f 0.12de 0.32b 0.55e 2.25f 3.26e BRS 133 5.66b 8.78b 69.10a 94.01b 183.10a 0.21b 0.30bc 4.94a 10.21b 15.84a BRS 213 8.49a 20.25a 38.71b 76.67c 144.12c 0.47a 1.09a 2.54b 7.23c 11.33c BRS 232 0.55f 2.15f 27.52d 68.69d 98.92e 0.18c 0.16e 1.84d 6.47d 8.64d BRS 284 0.72e 5.15d 26.95d 112.58a 151.01b 0.10e 0.27cd 1.93cd 11.02a 13.32b * Individual isoflavone glycosides were normalized for their molecular weight differences into aglycone forms and summed to obtain the total content. Means with different lowercase letters in the same column are significantly different (p < 0.05) by Tukey’s test. Results expressed as the average of six soybean replicates and twelve soymilk replicates . The total isoflavone content of soymilk ranged from 8.07 to 15.84 mg 100 g-1 (wet basis) (Table 4), with the isoflavone profile showing a similar pattern to that of soybean. However, from the total isoflavone contents of soybean and soymilk (on a dry basis), it was possible to observe (Figure 1) different transference ratios of isoflavone from soybean to soymilk. The highest ratio between soymilk and soybean isoflavones was observed for conventional cultivars BRS 232 and BRS 284, which showed the highest soymilk yield (Table 3). Grain softening during processing favored the extraction of both proteins and isoflavones into soymilk. Figure 1 Total isoflavone content of soybean and soymilk (mg.100 g-1, dry basis). 3.2 Consumer acceptance of soymilk There were significant differences (p < 0.05) with regard to the overall consumer acceptance (n = 102) among soybean cultivars (Table 5), which ranged from 5.9 to 6.7, and, therefore, in the “like” portion of the 9-point hedonic scale (1:disliked it extremely and 9:liked it extremely). The cultivar BRS 232 showed the highest acceptance, which did not differ from the BRS 284. The conventional cultivars BRS 284 and BRS 133 did not differ from the special BRS 216 and BRS 213, the latter of which had no lipoxygenases. This result suggests that the process used to obtain the soymilk, with blanching and grinding, inactivated the lipoxygenase enzymes and contributed to mitigating the differences between cultivars. Table 5 Average consumer acceptance#: overall and the identified segments' results. Cultivars Acceptance Overall (n = 102) Cluster 1 (n = 28) Cluster 2 (n = 32) Cluster 3 (n = 35) BRS 213 6.3b 5.4dC 7.5aA 5.8dB BRS 232 6.7a 6.3bB 6.6bA 7.0bA BRS 284 6.5ab 6.8aA 6.1cC 6.4cB BRS 133 6.4b 4.9cC 6.4bB 7.4aA BRS 216 6.3b 5.6dB 6.6bA 6.3cA BRS 267 5.9c 6.7aA 6.2cB 4.7eC # Evaluated on a 9-point hedonic scale (1: disliked it extremely and 9: liked it extremely). Means followed by lowercase letters in the columns and uppercase letters in the rows do not differ by Tukey’s test (p ≤ 0.05). The cluster analysis identified three consumer segments according to preference similarity. The acceptance means are shown in Table 5. The results of the internal preference mapping, a tool that considers the individual perception of participants, are shown in Figures 22B. The first two dimensions explained 52.3% of the variance, which is expected taking into account that the study was carried out by untrained people (consumers). Dimension 1 separated the cultivar BRS 267 from the others. In contrast, dimension 2 separated the soymilk prepared with cultivars BRS 284 and BRS 133. The soymilk prepared with BRS 267 was the least appreciated by most participants of this study (segments 2 and 3), differing from the conclusion reported by Silva et al. (2009) in which this cultivar was the best cultivar for soymilk preparation among the same cultivars evaluated in this work based on its chemical composition. However, this cultivar reached the highest mean among the consumers of segment 1 as well as BRS 284. BRS 267 is a vegetable-type cultivar with a more appropriate flavor to be consumed as edamame (Carrão-Panizzi et al., 2018; Felberg et al., 2020). Figure 2 Internal preference map showing (A): the position of the samples and (B): the position of the consumers in the three clusters. The consumers of segment 2 preferred the cultivar without lipoxygenases, BRS 213 (acceptance mean of 7.5). On the other hand, the consumers of segment 3 preferred the cultivar BRS 133 (acceptance mean of 7.4) which showed the highest total isoflavone content. These results highlight the effect of processing over the intrinsic characteristics of the cultivars, reducing the effect of the specialty cultivars in the soymilk. On the other hand, the internal preference mapping and cluster analysis showed the segmentation of different consumer groups which are useful tools to select among cultivars and products. 4 Conclusions The studied cultivars differed with regard to nutrients and bioactive compounds, remaining within the ranges expected for soybean. There was a difference among the cultivars related to the transference rates of protein and isoflavones from soybean to soymilk as well as the soymilk yield. Thermal treatment applied to soybean grains and soymilk inactivated the enzymes when these were present, and the soymilks did not differ from the one obtained of the cultivar without lipoxygenase. The results indicated that the conventional cultivars and specialty cultivars are promising for preparing the soymilk from the perspective of consumer acceptance since the technological process and thermal treatment attenuated the differences among cultivars. Acknowledgements The authors thank Embrapa for funding, and Joana de Novais Pereira and David Régis de Oliveira for their technical assistance and Andre L. N. Gomes and Marcos O. Moulin for assistance with the figures. Cite as: Felberg, I., Carrão-Panizzi, M. C., Deliza, R., Freitas, S. C., Gonçalves, E. B., Pacheco, S., Santiago, M. C. P. A., & Antoniassi, R. (2023). Nutritional and consumer acceptance evaluation of soymilk from specialty and conventional soybean cultivars. Brazilian Journal of Food Technology, 26, e2022075. https://doi.org/10.1590/1981-6723.07522 Funding: Empresa Brasileira de Pesquisa Agropecuária (EMBRAPA). References Al-Mahfuz A. Tsukamoto C. Kudou S. Ono T. 2004 Chances of adstringent sensation of soy milk during tofu curd formation Journal of Agricultural and Food Chemistry 52 23 7070 7074 http://dx.doi.org/10.1021/jf0491557 Al-Mahfuz, A., Tsukamoto, C., Kudou, S., & Ono, T. (2004). Chances of adstringent sensation of soy milk during tofu curd formation. Journal of Agricultural and Food Chemistry, 52(23), 7070-7074. http://dx.doi.org/10.1021/jf0491557 Antoniassi R. Felberg I. de Aguiar P. F. de Freitas S. C. de Aguiar A. C. Pereira L. M. Mesquita D. L. 2008 Efeito do cozimento na inativação de lipoxigenases da variedade de soja IAS-5. Congresso Brasileiro de Ciência e Tecnologia de Alimentos. Belo Horizonte SBCTA Antoniassi, R., Felberg, I., de Aguiar, P. F., de Freitas, S. C., de Aguiar, A. C., Pereira, L. M., & Mesquita, D. L. (2008). Efeito do cozimento na inativação de lipoxigenases da variedade de soja IAS-5. In: Congresso Brasileiro de Ciência e Tecnologia de Alimentos. Belo Horizonte: SBCTA. Association of Official Analytical Chemists AOAC International 2010 Official methods of analysis of AOAC International 18th Gaithersburg AOAC International Association of Official Analytical Chemists – AOAC International (2010). Official methods of analysis of AOAC International18thGaithersburgAOAC International Axelrod B. Cheesbrough T. M. Laakso S. 1981 Lipoxygenases in soybean Methods in Enzymology 71 441 451 http://dx.doi.org/10.1016/0076-6879(81)71055-3 Axelrod, B., Cheesbrough, T. M., & Laakso, S. (1981). Lipoxygenases in soybean. Methods in Enzymology, 71, 441-451. http://dx.doi.org/10.1016/0076-6879(81)71055-3 Bhagwat S. Haytowitz D. B. Holden J. M. 2015 USDA Database for the Isoflavone Content of Selected Foods Release 2.0 Retrieved in 2022, July 7 from http://www.ars.usda.gov/SP2UserFiles/Place/80400525/ Data/isoflav/Isoflav_R2.pdf Bhagwat, S., Haytowitz, D. B., & Holden, J. M. (2015). USDA Database for the Isoflavone Content of Selected Foods Release 2.0. Retrieved in 2022, July 7, from http://www.ars.usda.gov/SP2UserFiles/Place/80400525/ Data/isoflav/Isoflav_R2.pdf Cai J. S. Feng J. Y. Ni Z. J. Ma R. H. Thakur K. Wang S. Hu F. Zhang J. F. Wei Z. J. 2021 An update on the nutritional functional, sensory characteristics of soy products and applications of new processing strategies Trends in Food Science & Technology 112 676 689 http://dx.doi.org/10.1016/j.tifs.2021.04.039 Cai, J. S., Feng, J. Y., Ni, Z. J., Ma, R. H., Thakur, K., Wang, S., Hu, F., Zhang, J. F., & Wei, Z. J. (2021). An update on the nutritional functional, sensory characteristics of soy products and applications of new processing strategies. Trends in Food Science & Technology, 112, 676-689. http://dx.doi.org/10.1016/j.tifs.2021.04.039 Carrão-Panizzi M. C. Almeida L. A. Souza Kiihl R. A. Miranda L. C. Kikuchi A. 2001 Developing sustainable agricultural systems: determinants, future approaches and roles of different partners, as viewed from the soybean breeding program for human nutrition, at the National Soybean Research Center of EMBRAPA. JIRCAS International Symposium Series Japan Retrieved in 2022, July 7 from https://www.jircas.go.jp/en/publication/intlsymp/9/107 Carrão-Panizzi, M. C., Almeida, L. A., Souza Kiihl, R. A., Miranda, L. C., & Kikuchi, A. (2001). Developing sustainable agricultural systems: determinants, future approaches and roles of different partners, as viewed from the soybean breeding program for human nutrition, at the National Soybean Research Center of EMBRAPA. In JIRCAS International Symposium Series, Japan. Retrieved in 2022, July 7, from https://www.jircas.go.jp/en/publication/intlsymp/9/107 Carrão-Panizzi M. C. Berhow M. Mandarino J. M. G. Oliveira M. C. N. 2009 Environmental and genetic variation of isoflavone content of soybean seeds grown in Brazil Pesquisa Agropecuária Brasileira 44 11 1444 1451 http://dx.doi.org/10.1590/S0100-204X2009001100011 Carrão-Panizzi, M. C., Berhow, M., Mandarino, J. M. G., & Oliveira, M. C. N. (2009). Environmental and genetic variation of isoflavone content of soybean seeds grown in Brazil. Pesquisa Agropecuária Brasileira, 44(11), 1444-1451. http://dx.doi.org/10.1590/S0100-204X2009001100011 Carrão-Panizzi M. C. Oliveira M. A. Santos H. P. Mandarino J. M. G. Oliveira M. C. N. 2018 Boletim de Pesquisa e Desenvolvimento: Características de vagens e grãos de cultivares de soja para utilização como edamame Londrina Embrapa Soja Retrieved in 2022, July 7 from https://ainfo.cnptia.embrapa.br/digital/bitstream/item/190841/1/BOLETIM-20.pdf Carrão-Panizzi, M. C., Oliveira, M. A., Santos, H. P., Mandarino, J. M. G., & Oliveira, M. C. N. (2018). Boletim de Pesquisa e Desenvolvimento: Características de vagens e grãos de cultivares de soja para utilização como edamame. Londrina: Embrapa Soja. Retrieved in 2022, July 7, from https://ainfo.cnptia.embrapa.br/digital/bitstream/item/190841/1/BOLETIM-20.pdf Carrão-Panizzi M. C. Pires J. L. F. Strieder Leite R. S. Oliveira M. C. N. Corassa G. M. 2021 Teores de óleo e proteína em genótipos de soja em diferentes situações de manejo Passo Fundo Embrapa Trigo Retrieved in 2022, July 7 from http://ainfo.cnptia.embrapa.br/digital/bitstream/item/222440/1/CircTec-60-online-2021.pdf Carrão-Panizzi, M. C., Pires, J. L. F., Strieder, M. L., Leite, R. S., Oliveira, M. C. N., & Corassa, G. M. (2021). Teores de óleo e proteína em genótipos de soja em diferentes situações de manejo. Passo Fundo: Embrapa Trigo. Retrieved in 2022, July 7, from http://ainfo.cnptia.embrapa.br/digital/bitstream/item/222440/1/CircTec-60-online-2021.pdf Ciabotti S. Juhász A. C. P. Mandarino J. M. G. Costa L. L. Corrêa A. D. Simão A. A. Santos E. N. F. 2019 Chemical composition and lipoxygenase activity of soybean (Glycine max L. Merrill.) genotypes, specific for human consumption, with different tegument colours Brazilian Journal of Food Technology 22 22 http://dx.doi.org/10.1590/1981-6723.00318 Ciabotti, S., Juhász, A. C. P., Mandarino, J. M. G., Costa, L. L., Corrêa, A. D., Simão, A. A., & Santos, E. N. F. (2019). Chemical composition and lipoxygenase activity of soybean (Glycine max L. Merrill.) genotypes, specific for human consumption, with different tegument colours. Brazilian Journal of Food Technology, 22, 22. http://dx.doi.org/10.1590/1981-6723.00318 Dietary Guidelines for Americans DGA 2020 DGA: Make every bite count with the dietary guidelines. 9th ed. Washington, DC USDA Retrieved in 2022, July 7 from https://www.dietaryguidelines.gov/sites/default/files/2020-12/Dietary_Guidelines_for_Americans_2020-2025.pdf Dietary Guidelines for Americans – DGA. 2020. DGA: Make every bite count with the dietary guidelines. (9th ed.). Washington, DC: USDA. Retrieved in 2022, July 7, from https://www.dietaryguidelines.gov/sites/default/files/2020-12/Dietary_Guidelines_for_Americans_2020-2025.pdf Felberg I. Antoniassi R. Deliza Freitas S. C. Modesta R. C. D. 2009 Soy and Brazil nut beverage: Processing, composition, sensory and color evaluation Ciência e Tecnologia de Alimentos 29 3 609 617 http://dx.doi.org/10.1590/S0101-20612009000300024 Felberg, I., Antoniassi, R., Deliza, R., Freitas, S. C., & Modesta, R. C. D. (2009). Soy and Brazil nut beverage: Processing, composition, sensory and color evaluation. Ciência e Tecnologia de Alimentos, 29(3), 609-617. http://dx.doi.org/10.1590/S0101-20612009000300024 Felberg I. Deliza R. Gonçalves E. B. Antoniassi R. de Freitas S. C. 2004 Bebida mista de extrato de soja integral e castanha-do-brasil: Caracterização físico-química, nutricional e aceitabilidade do consumidor Alimentos e Nutrição 15 2 163 174 Felberg, I., Deliza, R., Gonçalves, E. B., Antoniassi, R., & de Freitas, S. C. (2004). Bebida mista de extrato de soja integral e castanha-do-brasil: Caracterização físico-química, nutricional e aceitabilidade do consumidor. Alimentos e Nutrição, 15(2), 163-174. Felberg I. Torrezan R. Pinto M. S. V. Nascimento F. Neto 2020 Processamento de edamame em vagens para agroindústria de pequeno porte Rio de Janeiro Embrapa Agroindústria de Alimentos Felberg, I., Torrezan, R., Pinto, M. S. V., & Nascimento Neto, F. (2020). Processamento de edamame em vagens para agroindústria de pequeno porte. Rio de Janeiro: Embrapa Agroindústria de Alimentos. Han H. Choi J. K. Park J. Im H. C. Han J. H. Huh M. H. Lee Y. B. 2021 Recent innovations in processing technologies for improvement of nutritional quality of soymilk CYTA: Journal of Food 19 1 287 303 http://dx.doi.org/10.1080/19476337.2021.1893824 Han, H., Choi, J. K., Park, J., Im, H. C., Han, J. H., Huh, M. H., & Lee, Y. B. (2021). Recent innovations in processing technologies for improvement of nutritional quality of soymilk. CYTA: Journal of Food, 19(1), 287-303. http://dx.doi.org/10.1080/19476337.2021.1893824 Kumar V. Rani A. Tindwani C. Jain M. 2003 Lipoxygenase isozymes and trypsin inhibitor activities in soybean as influenced by growing location Food Chemistry. London. 83 1 79 83 http://dx.doi.org/10.1016/S0308-8146(03)00052-9 Kumar, V., Rani, A., Tindwani, C., & Jain, M. (2003). Lipoxygenase isozymes and trypsin inhibitor activities in soybean as influenced by growing location. Food Chemistry. London., 83(1), 79-83. http://dx.doi.org/10.1016/S0308-8146(03)00052-9 Li N. Wu X. Zhuang W. Xia L. Chen Y. Zhao R. Yi M. Wan Q. Du L. Zhou Y. 2020 Soy and isoflavone consumption and multiple health outcomes: Umbrella review of systematic reviews and meta-analyses of observational studies and randomized trials in humans Molecular Nutrition & Food Research 64 4 e1900751 31584249. http://dx.doi.org/10.1002/mnfr.201900751 Li, N., Wu, X., Zhuang, W., Xia, L., Chen, Y., Zhao, R., Yi, M., Wan, Q., Du, L., & Zhou, Y. (2020). Soy and isoflavone consumption and multiple health outcomes: Umbrella review of systematic reviews and meta-analyses of observational studies and randomized trials in humans. Molecular Nutrition & Food Research, 64(4), e1900751. PMid:31584249. http://dx.doi.org/10.1002/mnfr.201900751 Ludikhuyze L. Indrawati Van den Broeck I. Weemaes C. Hendrickx M. 1998 Effect of combined pressure and temperature of soybean lipoxygenase 2 Modeling inactivation kinetics under static and dynamic conditions Journal of Agricultural and Food Chemistry 46 10 4081 4086 http://dx.doi.org/10.1021/jf9802575 Ludikhuyze, L., Indrawati, Van den Broeck, I., Weemaes, C., & Hendrickx, M. (1998). Effect of combined pressure and temperature of soybean lipoxygenase 2 Modeling inactivation kinetics under static and dynamic conditions. Journal of Agricultural and Food Chemistry, 46(10), 4081-4086. http://dx.doi.org/10.1021/jf9802575 MacFie H. J. H. Bratchell N. Greenhoff K. Vallis L. V. 1989 Designs to balance the effect of order of presentation and first-order carry-over effects in hall tests Journal of Sensory Studies 4 2 129 148 http://dx.doi.org/10.1111/j.1745-459X.1989.tb00463.x MacFie, H. J. H., Bratchell, N., Greenhoff, K., & Vallis, L. V. (1989). Designs to balance the effect of order of presentation and first-order carry-over effects in hall tests. Journal of Sensory Studies, 4(2), 129-148. http://dx.doi.org/10.1111/j.1745-459X.1989.tb00463.x Macrae R. 1998 HPLC in food analysis. 2nd New York Academic Press Macrae, R. (1998). HPLC in food analysis. (2nd ed.). New York: Academic Press. Messina M. 2016 Soy and health update: Evaluation of the clinical and epidemiologic literature review Nutrients 8 12 754 http://dx.doi.org/10.3390/nu8120754 Messina, M. (2016). Soy and health update: Evaluation of the clinical and epidemiologic literature review. Nutrients, 8(12), 754. http://dx.doi.org/10.3390/nu8120754 Nachvak S. M. Moradi S. Anjom-Shoae J. Rahmani J. Nasiri M. Maleki V. Sadeghi O. 2019 Soy, soy isoflavones and protein intake in relation to mortality from all causes cancers and cardiovascular diseases: A systematic review and dose-response meta-analysis of prospective cohort studies Journal of the Academy of Nutrition and Dietetics 119 9 1483 1500.e17 31278047. http://dx.doi.org/10.1016/j.jand.2019.04.011 Nachvak, S. M., Moradi, S., Anjom-Shoae, J., Rahmani, J., Nasiri, M., Maleki, V., & Sadeghi, O. (2019). Soy, soy isoflavones and protein intake in relation to mortality from all causes cancers and cardiovascular diseases: A systematic review and dose-response meta-analysis of prospective cohort studies. Journal of the Academy of Nutrition and Dietetics, 119(9), 1483-1500.e17. PMid:31278047. http://dx.doi.org/10.1016/j.jand.2019.04.011 Nelson A. I. Steinberg M. P. Wei L. S. 1976 Illinois process for preparation of soymilk Journal of Food Science. Chicago. Ill. 41 1 57 61 http://dx.doi.org/10.1111/j.1365-2621.1976.tb01100.x Nelson, A. I., Steinberg, M. P., & Wei, L. S. (1976). Illinois process for preparation of soymilk. Journal of Food Science. Chicago. Ill., 41(1), 57-61. http://dx.doi.org/10.1111/j.1365-2621.1976.tb01100.x Oliveira M. Carrão-Panizzi M. Mandarino J. Leite R. Campos P. Filho Vicentini M. 2010 Quantificação dos teores de açúcares, oligossacarídeos e amido em genótipos/ cultivares de soja (Glycine Max (L) Merril) especiais utilizados para alimentação humana Brazilian Journal of Food Technology 13 1 23 29 http://dx.doi.org/10.4260/BJFT2010130100004 Oliveira, M., Carrão-Panizzi, M., Mandarino, J., Leite, R., Campos Filho, P., & Vicentini, M. (2010). Quantificação dos teores de açúcares, oligossacarídeos e amido em genótipos/ cultivares de soja (Glycine Max (L) Merril) especiais utilizados para alimentação humana. Brazilian Journal of Food Technology, 13(1), 23-29. http://dx.doi.org/10.4260/BJFT2010130100004 Sansai K. Na Takuathung M. Khatsri R. Teekachunhatean S. Hanprasertpong N. Koonrungsesomboon N. 2020 Effects of isoflavone interventions on bone mineral density in postmenopausal women: A systematic review and meta-analysis of randomized controlled trials Osteoporosis International 31 10 1853 1864 32524173. http://dx.doi.org/10.1007/s00198-020-05476-z Sansai, K., Na Takuathung, M., Khatsri, R., Teekachunhatean, S., Hanprasertpong, N., & Koonrungsesomboon, N. (2020). Effects of isoflavone interventions on bone mineral density in postmenopausal women: A systematic review and meta-analysis of randomized controlled trials. Osteoporosis International, 31(10), 1853-1864. PMid:32524173. http://dx.doi.org/10.1007/s00198-020-05476-z Silva J. B. D. Prudêncio S. H. Felberg I. Deliza R. Carrão-Panizz M. C. 2007 Aceitabilidade de bebidas preparadas a partir de diferentes extratos hidrossolúveis de soja Pesquisa Agropecuária Brasileira 42 12 1779 1784 http://dx.doi.org/10.1590/S0100-204X2007001200016 Silva, J. B. D., Prudêncio, S. H., Felberg, I., Deliza, R., & Carrão-Panizz, M. C. (2007). Aceitabilidade de bebidas preparadas a partir de diferentes extratos hidrossolúveis de soja. Pesquisa Agropecuária Brasileira, 42(12), 1779-1784. http://dx.doi.org/10.1590/S0100-204X2007001200016 Silva J. B. Carrão-Panizzi M. C. Prudêncio S. H. 2009 Chemical and physical composition of grain-type and food-type soybean for food processing Pesquisa Agropecuária Brasileira 44 7 777 784 http://dx.doi.org/10.1590/S0100-204X2009000700019 Silva, J. B., Carrão-Panizzi, M. C., & Prudêncio, S. H. (2009). Chemical and physical composition of grain-type and food-type soybean for food processing. Pesquisa Agropecuária Brasileira, 44(7), 777-784. http://dx.doi.org/10.1590/S0100-204X2009000700019 Song T. Barua K. Buseman G. Murphy P. A. 1998 Soy isoflavone analysis: Quality control and a new internal standard The American Journal of Clinical Nutrition 68 6 1474S 1479S http://dx.doi.org/10.1093/ajcn/68.6.1474S Song, T., Barua, K., Buseman, G., & Murphy, P. A. (1998). Soy isoflavone analysis: Quality control and a new internal standard. The American Journal of Clinical Nutrition, 68(6), 1474S-1479S. http://dx.doi.org/10.1093/ajcn/68.6.1474S Statistical Analysis System Institute SAS 2008 SAS /STAT 9.2 users’s guide. Cary SAS Statistical Analysis System Institute – SAS. (2008). SAS /STAT 9.2 users’s guide. Cary: SAS. Statistical Analysis System Institute SAS 2009 Base SAS 9.2 procedures guide. Cary SAS Statistical Analysis System Institute – SAS. (2009). Base SAS 9.2 procedures guide. Cary: SAS. Wang H. J. Murphy P. A. 1994 Isoflavone composition of American and Japanese soybeans in Iowa: Effects of variety, crop year and location Journal of Agricultural and Food Chemistry 42 8 1674 1677 http://dx.doi.org/10.1021/jf00044a017 Wang, H. J., & Murphy, P. A. (1994). Isoflavone composition of American and Japanese soybeans in Iowa: Effects of variety, crop year and location. Journal of Agricultural and Food Chemistry, 42(8), 1674-1677. http://dx.doi.org/10.1021/jf00044a017 Wang R. Zhou X. Chen Z. 2008 High pressure inactivation of lipoxygenase in soy milk and crude soybean extract Food Chemistry 106 2 603 611 http://dx.doi.org/10.1016/j.foodchem.2007.06.056 Wang, R., Zhou, X., & Chen, Z. (2008). High pressure inactivation of lipoxygenase in soy milk and crude soybean extract. Food Chemistry, 106(2), 603-611. http://dx.doi.org/10.1016/j.foodchem.2007.06.056