<?xml version="1.0" encoding="utf-8"?>
<!DOCTYPE article PUBLIC "-//NLM//DTD Journal Publishing DTD v2.3 20070202//EN" "journalpublishing.dtd">
<article xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" article-type="research-article" dtd-version="2.3" xml:lang="EN">
<front>
<journal-meta>
<journal-id journal-id-type="publisher-id">Front. Sustain. Food Syst.</journal-id>
<journal-title>Frontiers in Sustainable Food Systems</journal-title>
<abbrev-journal-title abbrev-type="pubmed">Front. Sustain. Food Syst.</abbrev-journal-title>
<issn pub-type="epub">2571-581X</issn>
<publisher>
<publisher-name>Frontiers Media S.A.</publisher-name>
</publisher>
</journal-meta>
<article-meta>
<article-id pub-id-type="doi">10.3389/fsufs.2023.1223921</article-id>
<article-categories>
<subj-group subj-group-type="heading">
<subject>Sustainable Food Systems</subject>
<subj-group>
<subject>Original Research</subject>
</subj-group>
</subj-group>
</article-categories>
<title-group>
<article-title>Western Hemisphere quality and production capacity of soybean protein</article-title>
</title-group>
<contrib-group>
<contrib contrib-type="author" corresp="yes"><name><surname>Cerrudo</surname> <given-names>Anibal</given-names></name><xref rid="aff1" ref-type="aff"><sup>1</sup></xref><xref rid="aff2" ref-type="aff"><sup>2</sup></xref><xref rid="c001" ref-type="corresp"><sup>&#x002A;</sup></xref>
<uri xlink:href="https://loop.frontiersin.org/people/2313504/overview"/>
</contrib>
<contrib contrib-type="author"><name><surname>Miller-Garvin</surname> <given-names>Jill</given-names></name><xref rid="aff1" ref-type="aff"><sup>1</sup></xref>
</contrib>
<contrib contrib-type="author"><name><surname>Naeve</surname> <given-names>Seth L.</given-names></name><xref rid="aff1" ref-type="aff"><sup>1</sup></xref>
<uri xlink:href="https://loop.frontiersin.org/people/1470587/overview"/>
</contrib>
</contrib-group>
<aff id="aff1"><sup>1</sup><institution>Department of Agronomy and Plant Genetics, University of Minnesota</institution>, <addr-line>Saint Paul, MN</addr-line>, <country>United States</country></aff>
<aff id="aff2"><sup>2</sup><institution>Ecofisiolog&#x00ED;a de cultivos, Unidad Integrada Balcarce (INTA-FCA)</institution>, <addr-line>Balcarce</addr-line>, <country>Argentina</country></aff>
<author-notes>
<fn fn-type="edited-by" id="fn0005">
<p>Edited by: Krishnanand P. Kulkarni, Delaware State University, United States</p>
</fn>
<fn fn-type="edited-by" id="fn0006">
<p>Reviewed by: Sapna Langyan, Indian Council of Agricultural Research (ICAR), India; Somak Dutta, Iowa State University, United States</p>
</fn>
<corresp id="c001">&#x002A;Correspondence: Anibal Cerrudo, <email>cerru006@umn.edu</email></corresp>
</author-notes>
<pub-date pub-type="epub">
<day>07</day>
<month>11</month>
<year>2023</year>
</pub-date>
<pub-date pub-type="collection">
<year>2023</year>
</pub-date>
<volume>7</volume>
<elocation-id>1223921</elocation-id>
<history>
<date date-type="received">
<day>17</day>
<month>05</month>
<year>2023</year>
</date>
<date date-type="accepted">
<day>12</day>
<month>10</month>
<year>2023</year>
</date>
</history>
<permissions>
<copyright-statement>Copyright &#x00A9; 2023 Cerrudo, Miller-Garvin and Naeve.</copyright-statement>
<copyright-year>2023</copyright-year>
<copyright-holder>Cerrudo, Miller-Garvin and Naeve</copyright-holder>
<license xlink:href="http://creativecommons.org/licenses/by/4.0/">
<p>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.</p>
</license>
</permissions>
<abstract>
<p>Soybean&#x2019;s protein concentration and amino acid profile vary geographically due to genotype and environment. We used 1817 on-farm grain samples from the United States, Brazil, and the Southern Cone to assess soybean protein variability in the Western Hemisphere, which accounts for 85% of global production. Our analysis found consistent spatial patterns for protein concentration and amino acid composition, with differences in protein concentration of up to 5.4% among regions. While there was a negative correlation between protein concentration and the abundance of critical amino acids (lysine, cysteine, threonine, methionine, and tryptophane), their amount per ton of grain increased with protein concentration. Our findings can guide market and industry in benchmarking soybean protein quality across the Western Hemisphere, which according to our projection could supply close to 120 million metric tons of protein and 17 million metric tons of the most critical amino acids annually by 2030.</p>
</abstract>
<abstract abstract-type="graphical">
<title>Graphical abstract</title>
<p>
<graphic xlink:href="fsufs-07-1223921gr0001.tif" xmlns:xlink="http://www.w3.org/1999/xlink"/>
</p>
</abstract>
<kwd-group>
<kwd>grain composition</kwd>
<kwd>amino acid profile</kwd>
<kwd>essentials amino acids</kwd>
<kwd>crude protein</kwd>
<kwd>spatial pattern</kwd>
<kwd>Western Hemisphere</kwd>
</kwd-group>
<counts>
<fig-count count="7"/>
<table-count count="2"/>
<equation-count count="0"/>
<ref-count count="49"/>
<page-count count="10"/>
<word-count count="6979"/>
</counts>
<custom-meta-wrap>
<custom-meta>
<meta-name>section-at-acceptance</meta-name>
<meta-value>Crop Biology and Sustainability</meta-value>
</custom-meta>
</custom-meta-wrap>
</article-meta>
</front>
<body>
<sec sec-type="intro" id="sec1">
<label>1.</label>
<title>Introduction</title>
<p>Soybean [<italic>Glycine max</italic> (L.) Merr.] is a major protein source for monogastric livestock and comprises 70% of the world&#x2019;s protein-meal consumption (<xref ref-type="bibr" rid="ref11">FAO Faostat, 2022</xref>). Surprisingly, only 2% of soybean protein is consumed directly as human food, with the remaining 98% used for animal feed (<xref ref-type="bibr" rid="ref13">Goldsmith et al., 2004</xref>). Driven by the growing demand for oil and animal protein (<xref ref-type="bibr" rid="ref12">Godfray and Garnett, 2014</xref>), soybean production in the Americas has doubled in the last two decades, mainly through an expansion in the cropping area (<xref ref-type="bibr" rid="ref51">Zalles et al., 2019</xref>; <xref ref-type="bibr" rid="ref43">Song et al., 2021</xref>; <xref ref-type="bibr" rid="ref23">Marin et al., 2022</xref>). Currently, the United States, Brazil, and the Southern Cone (Argentina, Paraguay, and Uruguay) account for 95% of soybean production in the continent. It is worth noting that almost one-third of the arable land in the Western Hemisphere (95 million ha) is used to produce over 85% of the world&#x2019;s soybean grain (<xref ref-type="bibr" rid="ref26">OECD and FAO, 2022</xref>).</p>
<p>According to the relevance of the feed end-use of soybean and because modern monogastric rations are oriented to maximize feed conversion efficiency, the market demands high-protein concentration soybeans (<xref ref-type="bibr" rid="ref15">Hertsgaard et al., 2019</xref>; <xref ref-type="bibr" rid="ref32">Pope et al., 2023</xref>). In addition, the ultimate nutritional value of soybean protein is determined by the proportion of each amino acid in total protein, the amino acid profile (<xref ref-type="bibr" rid="ref27">Osborne et al., 1914</xref>; <xref ref-type="bibr" rid="ref44">Tessari et al., 2016</xref>). Optimizing the amino acid profile can lead to less excess amino acid deamination and excretion of nitrogen, resulting in increased conversion efficiency of vegetable into animal protein (<xref ref-type="bibr" rid="ref47">Waguespack et al., 2009</xref>; <xref ref-type="bibr" rid="ref4">Attia et al., 2020</xref>). Soybean protein analysis typically reports 18 amino acids (<xref ref-type="bibr" rid="ref19">Kovalenko et al., 2006</xref>), with lysine, threonine, tryptophan, cysteine, and methionine considered the most critical amino acids as they frequently limit the conversion efficiency of soybean meal in monogastric animal products (<xref ref-type="bibr" rid="ref45">Thakur and Hurburgh, 2007</xref>; <xref ref-type="bibr" rid="ref34">Ravindran et al., 2014</xref>).</p>
<p>Soybean grain composition is determined by a combination of environmental resources such as radiation, water, and nutrients, as well as regulators like temperature, all of which interact with the genotype (<xref ref-type="bibr" rid="ref30">Piper and Boote, 1999</xref>; <xref ref-type="bibr" rid="ref9">Carrera et al., 2009</xref>, <xref ref-type="bibr" rid="ref10">2022</xref>; <xref ref-type="bibr" rid="ref36">Rotundo and Westgate, 2009</xref>; <xref ref-type="bibr" rid="ref1">Anthony et al., 2012</xref>; <xref ref-type="bibr" rid="ref29">Pfarr et al., 2018</xref>; <xref ref-type="bibr" rid="ref7">Bosaz et al., 2019</xref>). The spatial structure of these factors varies due to differences in soil fertility, precipitation, temperature, cultivar adoption, and other factors. As a result of these combined effects, soybean grain protein concentration exhibits geographical variability. Moreover, year-to-year variability also exists due to weather-determined environmental factors (<xref ref-type="bibr" rid="ref38">Rotundo et al., 2016</xref>).</p>
<p>There is some limited evidence that suggest that changes in the amino acid profile are linked to variations in protein concentration, and thus may be also influenced by the environment (<xref ref-type="bibr" rid="ref28">Park and Hurburgh, 2002</xref>; <xref ref-type="bibr" rid="ref42">Song et al., 2016</xref>; <xref ref-type="bibr" rid="ref2">Assefa et al., 2018</xref>). Controlled experiments have shown that as grain protein concentration increases, the abundance of critical amino acids decreases (<xref ref-type="bibr" rid="ref29">Pfarr et al., 2018</xref>). This is particularly relevant because regions with lower protein concentration may produce soybean meal of better quality (<xref ref-type="bibr" rid="ref49">William et al., 2020</xref>).</p>
<p>Knowledge about the occurrence of spatial patterns in grain composition across the Western Hemisphere would be relevant for the market, trading, and the industry (<xref ref-type="bibr" rid="ref15">Hertsgaard et al., 2019</xref>). Defining spatial patterns would help in the segregation of production according to protein concentration or amino acid profile, which would allow for the optimization of quality for different end-uses (<xref ref-type="bibr" rid="ref32">Pope et al., 2023</xref>). Different soybeans could even be combined in a way that the amino acid profile results in a higher protein score (<xref ref-type="bibr" rid="ref25">Messina, 2022</xref>). Additionally, it is essential to consider spatial variability in the estimation of current and future protein and amino acids production capacities, especially since the future increase in production may not be evenly distributed among regions (<xref ref-type="bibr" rid="ref26">OECD and FAO, 2022</xref>; <xref ref-type="bibr" rid="ref35">Ricard et al., 2022</xref>).</p>
<p>Previous studies have demonstrated regional variations in soybean protein concentrations in the U.S., Brazil, and the Southern Cone (<xref ref-type="bibr" rid="ref8">Breene et al., 1988</xref>; <xref ref-type="bibr" rid="ref37">Rotundo et al., 2014</xref>; <xref ref-type="bibr" rid="ref3">Assefa et al., 2019</xref>; <xref ref-type="bibr" rid="ref7">Bosaz et al., 2019</xref>). While attempts have been made to assess the protein profile of some amino acids in the U.S. (<xref ref-type="bibr" rid="ref2">Assefa et al., 2018</xref>), there is no comprehensive analysis of the geographical patterns for protein concentration and amino acid profiles across the entire Western Hemisphere. Therefore, in this study, we aim to (i) investigate the spatial patterns of soybean protein concentration and amino acid profiles, (ii) explore the interrelationships between amino acids and protein content in soybean grain, and (iii) project the soybean protein and amino acid production capacity of the Western Hemisphere.</p>
</sec>
<sec id="sec2">
<label>2.</label>
<title>Materials and methods</title>
<sec id="sec3">
<label>2.1.</label>
<title>Sites and years</title>
<p>We analyzed a soybean compositional dataset comprising 1817 grain samples collected during the 2012&#x2013;2014 growing seasons across the United States, Brazil, and the Southern Cone (Argentina, Paraguay, and Uruguay). To capture the agro-climatic and geographic diversity within these regions, we further divided them into subregions (<xref rid="fig1" ref-type="fig">Figure 1</xref>). In the United States, the subregions were the Northern Corn Belt, Western Corn Belt, Eastern Corn Belt, East Coast, and South. Brazil&#x2019;s subregions were North Brazil, Cerrado, Atlantic Forest, and South Brazil. In the Southern Cone, we defined the subregions as Paraguay, Uruguay, Northern Pampas, Central Pampas, and Southern Pampas.</p>
<fig position="float" id="fig1">
<label>Figure 1</label>
<caption>
<p>Regions and subregions across the Western Hemisphere where the soybean grain sampling was conducted between 2012 and 2014 growing seasons.</p>
</caption>
<graphic xlink:href="fsufs-07-1223921-g001.tif"/>
</fig>
</sec>
<sec id="sec4">
<label>2.2.</label>
<title>Sampling, compositional analysis, and derived variables</title>
<p>Soybean grain samples were collected from each region and subregion (<xref rid="fig1" ref-type="fig">Figure 1</xref>). Each year, samples were collected by state, province, or department trying to match actual production statistics, so that the resulting data would reflect the regional soybean compositional profile. Sampling varied slightly by country. United States samples were obtained directly from farmers; grain sample kits were mailed to farmers. South American soybean crops were sampled by independent contractors working in each country. Brazilian soybeans were sampled at the farm during harvest, while in the Southern Cone, samples were recovered at discharge locations including processors and ports. Production location was provided by the truck driver directly when samples were pulled. All grain samples were shipped to the University of Minnesota, St. Paul, Minnesota. The resulting dataset comprised 636, 645, and 536 samples from the United States, Brazil, and the Southern Cone, respectively. All samples were sent to the Experiment Station Chemical Laboratories at the University of Missouri-Columbia for analysis. Crude protein concentration was measured by combustion (LECO), crude oil concentration by ether extraction, total carbohydrates and the complete amino acid profile using AOAC Official Methods 990.30, 920.39, and 982.30 and expressed in dry basis. We estimated net protein concentration from the sum of all 18 amino acids and calculated the apparent non-protein N fraction for each sample as the difference between crude protein and net protein. The amount of the most critical amino acids was calculated as the sum of lysine, threonine, tryptophane, methionine, and cysteine (<xref ref-type="bibr" rid="ref45">Thakur and Hurburgh, 2007</xref>). Finally, we presented sulfur-containing amino acids (methionine and cysteine) together for some analyzes due to their sparing functionality in digestion (<xref ref-type="bibr" rid="ref6">Ball et al., 2006</xref>).</p>
</sec>
<sec id="sec5">
<label>2.3.</label>
<title>Statistical analysis</title>
<p>To assess the seed composition of different regions and subregions, average and percentile descriptors of crude protein, true protein, apparent non-protein N, and the amino acid profile were calculated across 3&#x2009;years of data. Normal distributions were observed for all seed composition variables. Variance component analysis was conducted using R (<xref ref-type="bibr" rid="ref33">R Core Team, 2021</xref>) to evaluate the contribution of the year, region, and subregion, as well as their interactions, to the variability registered for the different compositional variables. The statistical model included year as a random factor, region as a fixed factor, and subregion as a fixed factor that was nested within the region (nesting was necessary because subregions belong to only one region). Year by region and year by subregion interactions were considered random factors. The variance component analysis was performed using the VCA R-package (<xref ref-type="bibr" rid="ref39">Schuetzenmeister et al., 2022</xref>) that allows the extraction of variance components from both random and fixed effects in linear mixed-effects models. We assessed random effects&#x2019; significance in our mixed linear model through a likelihood ratio test, comparing the full model (with random effects) to a reduced model (without them) using the anova base R function (<xref ref-type="bibr" rid="ref5">Bakare et al., 2022</xref>).</p>
<p>Principal component analysis was used to identify any grouping between relevant traits focusing on crude protein, net protein, apparent non-protein N, oil concentration, carbohydrates concentration, and the most limiting amino acids (<xref ref-type="bibr" rid="ref17">Ju et al., 2019</xref>). Principal component analysis was performed using the FactoMineR package in R (<xref ref-type="bibr" rid="ref22">L&#x00EA; et al., 2008</xref>), only dimensions with eigenvalues &#x003E;1 were considered. Variables that made the greatest contribution to total variation captured by each dimension were identified. Linear regression and correlation analysis were generated using R software. We have thoroughly examined the assumptions of the linear regression model, including the independence of residuals, and found no evidence of violations.</p>
</sec>
<sec id="sec6">
<label>2.4.</label>
<title>Protein and amino acids production projection</title>
<p>To project protein and amino acid production by 2030, we combined regional estimations of grain composition from our analysis with the 2030 projection for soybean grain production at a country level (<xref ref-type="bibr" rid="ref26">OECD and FAO, 2022</xref>). Grain production at a subregion level was then estimated considering the relative contribution of each subregion to each country&#x2019;s production during the last 5-year period using data retrieved from each country&#x2019;s national agency (United States Department of Agriculture,<xref rid="fn0001" ref-type="fn">
<sup>1</sup></xref> Brazilian Institute of Geography and Statistics,<xref rid="fn0002" ref-type="fn">
<sup>2</sup></xref> the Argentinean Secretary of Agriculture,<xref rid="fn0003" ref-type="fn">
<sup>3</sup></xref> and the board of exporters and traders of cereals and oilseeds of Paraguay<xref rid="fn0004" ref-type="fn">
<sup>4</sup></xref>). Crude protein and most critical amino acids production projections for 2030 are presented in the results section, while the production estimations for all the individual 18 amino acids are presented in <xref ref-type="supplementary-material" rid="SM1">Supplementary Table S3</xref>.</p>
</sec>
</sec>
<sec sec-type="results" id="sec7">
<label>3.</label>
<title>Results</title>
<p>Our analysis of 1817 soybean seed samples revealed an average crude protein concentration of 39.6%, ranging from 36.1 to 42.6% defined by the fifth and 95th percentiles (P5-P95), indicating a 16% variation relative to the average. On average, the sum of the 18 primary amino acids (net protein) accounted for 97.2% of grain crude protein, while the remaining 2.8% was represented by apparent non-protein N. The fraction of the most critical amino acids in protein averaged 14.7%, ranging from 14.1 to 15.2% (P5-P95), indicating only a 7% variation relative to the average. The complete compositional information at a subregion level is presented in <xref ref-type="supplementary-material" rid="SM1">Supplementary Tables S1, S2</xref>.</p>
<sec id="sec8">
<label>3.1.</label>
<title>Crude protein, net protein, and non-protein nitrogen concentration</title>
<p>The results of the variance component analysis (<xref rid="tab1" ref-type="table">Table 1</xref>) indicated that for crude and net protein concentrations, the combined effects of region and subregion were the most significant sources of variation, collectively accounting for close to 30 and 33% of the total variation, respectively. Year did not affect crude or true protein concentrations. The interaction between the region and the year did not have a relevant effect on either protein fraction, but it was a significant source of variation for non-protein N concentration. In turn, the interaction between the subregion and the year was in general significant but explained only a small proportion of the total variation for crude protein, net protein, and non-protein N (less than 5% of the total sum of squares; <xref rid="tab1" ref-type="table">Table 1</xref>).</p>
<table-wrap position="float" id="tab1">
<label>Table 1</label>
<caption>
<p>Analysis of variance components for soybean grain composition across the Western Hemisphere.</p>
</caption>
<table frame="hsides" rules="groups">
<thead>
<tr>
<th align="left" valign="middle" rowspan="2">Source of variation</th>
<th align="center" valign="middle" colspan="2">Crude protein</th>
<th align="center" valign="middle" colspan="2">Net protein</th>
<th align="center" valign="middle" colspan="2">NpN</th>
<th align="center" valign="middle" colspan="2">CAA</th>
<th align="center" valign="middle" colspan="2">lys</th>
<th align="center" valign="middle" colspan="2">thr</th>
<th align="center" valign="middle" colspan="2">try</th>
<th align="center" valign="middle" colspan="2">cys&#x2009;+&#x2009;met</th>
</tr>
<tr>
<th align="center" valign="middle" colspan="16">%</th>
</tr>
</thead>
<tbody>
<tr>
<td align="left" valign="middle">Region</td>
<td align="char" valign="middle" char=".">20</td>
<td align="center" valign="middle">&#x002A;&#x002A;&#x002A;</td>
<td align="char" valign="middle" char=".">23</td>
<td align="center" valign="middle">&#x002A;&#x002A;&#x002A;</td>
<td align="char" valign="middle" char=".">1.1</td>
<td align="center" valign="middle">ns</td>
<td align="char" valign="middle" char=".">11</td>
<td align="center" valign="middle">&#x002A;</td>
<td align="char" valign="middle" char=".">7.4</td>
<td align="center" valign="middle">ns</td>
<td align="char" valign="middle" char=".">9.7</td>
<td align="center" valign="middle">ns</td>
<td align="char" valign="middle" char=".">0.1</td>
<td align="center" valign="middle">ns</td>
<td align="char" valign="middle" char=".">3.7</td>
<td align="center" valign="middle">&#x002A;&#x002A;</td>
</tr>
<tr>
<td align="left" valign="middle">Year</td>
<td align="char" valign="middle" char=".">0.9</td>
<td align="center" valign="middle">ns</td>
<td align="char" valign="middle" char=".">0.3</td>
<td align="center" valign="middle">ns</td>
<td align="char" valign="middle" char=".">11</td>
<td align="center" valign="middle">ns</td>
<td align="char" valign="middle" char=".">4.8</td>
<td align="center" valign="middle">ns</td>
<td align="char" valign="middle" char=".">3.3</td>
<td align="center" valign="middle">ns</td>
<td align="char" valign="middle" char=".">30</td>
<td align="center" valign="middle">ns</td>
<td align="char" valign="middle" char=".">36</td>
<td align="center" valign="middle">&#x002A;&#x002A;&#x002A;</td>
<td align="char" valign="middle" char=".">2.3</td>
<td align="center" valign="middle">ns</td>
</tr>
<tr>
<td align="left" valign="middle">Region&#x2009;&#x00D7;&#x2009;Year</td>
<td align="char" valign="middle" char=".">6.0</td>
<td align="center" valign="middle">ns</td>
<td align="char" valign="middle" char=".">12</td>
<td align="center" valign="middle">&#x002A;&#x002A;&#x002A;</td>
<td align="char" valign="middle" char=".">8.7</td>
<td align="center" valign="middle">&#x002A;&#x002A;&#x002A;</td>
<td align="char" valign="middle" char=".">3.4</td>
<td align="center" valign="middle">.</td>
<td align="char" valign="middle" char=".">35</td>
<td align="center" valign="middle">&#x002A;&#x002A;&#x002A;</td>
<td align="char" valign="middle" char=".">19</td>
<td align="center" valign="middle">&#x002A;&#x002A;&#x002A;</td>
<td align="char" valign="middle" char=".">2.3</td>
<td align="center" valign="middle">&#x002A;&#x002A;&#x002A;</td>
<td align="char" valign="middle" char=".">1.5</td>
<td align="center" valign="middle">&#x002A;&#x002A;</td>
</tr>
<tr>
<td align="left" valign="middle">Subregion</td>
<td align="char" valign="middle" char=".">9.3</td>
<td align="center" valign="middle">&#x002A;&#x002A;&#x002A;</td>
<td align="char" valign="middle" char=".">7.5</td>
<td align="center" valign="middle">&#x002A;&#x002A;&#x002A;</td>
<td align="char" valign="middle" char=".">1.0</td>
<td align="center" valign="middle">ns</td>
<td align="char" valign="middle" char=".">4.5</td>
<td align="center" valign="middle">ns</td>
<td align="char" valign="middle" char=".">3.2</td>
<td align="center" valign="middle">.</td>
<td align="char" valign="middle" char=".">2.7</td>
<td align="center" valign="middle">ns</td>
<td align="char" valign="middle" char=".">2.9</td>
<td align="center" valign="middle">&#x002A;&#x002A;&#x002A;</td>
<td align="char" valign="middle" char=".">2.8</td>
<td align="center" valign="middle">&#x002A;&#x002A;&#x002A;</td>
</tr>
<tr>
<td align="left" valign="middle">Subregion&#x2009;&#x00D7;&#x2009;Year</td>
<td align="char" valign="middle" char=".">3.5</td>
<td align="center" valign="middle">&#x002A;&#x002A;&#x002A;</td>
<td align="char" valign="middle" char=".">2.3</td>
<td align="center" valign="middle">&#x002A;&#x002A;&#x002A;</td>
<td align="char" valign="middle" char=".">3.4</td>
<td align="center" valign="middle">&#x002A;&#x002A;&#x002A;</td>
<td align="char" valign="middle" char=".">4.8</td>
<td align="center" valign="middle">&#x002A;&#x002A;&#x002A;</td>
<td align="char" valign="middle" char=".">4.0</td>
<td align="center" valign="middle">&#x002A;&#x002A;&#x002A;</td>
<td align="char" valign="middle" char=".">4.7</td>
<td align="center" valign="middle">&#x002A;&#x002A;&#x002A;</td>
<td align="char" valign="middle" char=".">1.4</td>
<td align="center" valign="middle">ns</td>
<td align="char" valign="middle" char=".">1.7</td>
<td align="center" valign="middle">ns</td>
</tr>
<tr>
<td align="left" valign="middle">Residuals</td>
<td align="char" valign="middle" char=".">60</td>
<td/>
<td align="char" valign="middle" char=".">56</td>
<td/>
<td align="char" valign="middle" char=".">75</td>
<td/>
<td align="char" valign="middle" char=".">72</td>
<td/>
<td align="char" valign="middle" char=".">47</td>
<td/>
<td align="char" valign="middle" char=".">34</td>
<td/>
<td align="char" valign="middle" char=".">57</td>
<td/>
<td align="char" valign="middle" char=".">88</td>
<td/>
</tr>
<tr>
<td align="left" valign="middle">CV (%)</td>
<td align="char" valign="middle" char=".">5.2</td>
<td/>
<td align="char" valign="middle" char=".">5.6</td>
<td/>
<td align="char" valign="middle" char=".">62</td>
<td/>
<td align="char" valign="middle" char=".">3.3</td>
<td/>
<td align="char" valign="middle" char=".">2.7</td>
<td/>
<td align="char" valign="middle" char=".">4.4</td>
<td/>
<td align="char" valign="middle" char=".">15</td>
<td/>
<td align="char" valign="middle" char=".">6.4</td>
<td/>
</tr>
</tbody>
</table>
<table-wrap-foot>
<p>Proportion of the total sum of squares explained by each source of variation expressed as percentage and significance level. Variables included were crude protein, net protein, and non-protein N (NpN) concentrations, as well as lysine (lys), threonine (thr), tryptophan (try), cysteine and methionine (cys&#x2009;+&#x2009;met), and the sum of these critical amino acids (CAA&#x2009;=&#x2009;lysine + threonine + tryptophane + cysteine + methionine) expressed as a fraction of net protein. The dataset includes 1,817 samples collected over 3&#x2009;years (2012&#x2013;2014) across 13 subregions in the United States, Brazil, and the Southern Cone regions.</p>
<p>&#x002A;Significant at the 0.05 probability level.</p>
<p>&#x002A;&#x002A;Significant at the 0.01 probability level.</p>
<p>&#x002A;&#x002A;&#x002A;Significant at the 0.001 probability level.</p>
<p>ns, not significant.</p>
</table-wrap-foot>
</table-wrap>
<p>Crude and net protein concentrations in grain, in average, were higher for Brazil, intermediate for the United States, and lower for the Southern Cone (<xref rid="fig2" ref-type="fig">Figure 2</xref>). The Southern Cone showed higher variability in both protein fractions compared to the other two regions. Within Brazil and the United States, the protein concentration varied in a narrow range among subregions, while within the Southern Cone, the Southern and Central Pampas subregions had relatively lower grain concentration for these grain constituents compared to Uruguay, the Northern Pampas, or Paraguay (<xref rid="fig3" ref-type="fig">Figure 3</xref>). There was no clear regional pattern for non-protein N concentration. We observed a weak positive correlation between non-protein N concentration and crude protein concentration (<italic>r</italic> =&#x2009;0.17, <italic>p</italic> &#x003C;&#x2009;0.01).</p>
<fig position="float" id="fig2">
<label>Figure 2</label>
<caption>
<p>Mean and distribution for soybean crude protein <bold>(A)</bold> and net protein <bold>(B)</bold> grain concentration, and the abundance of the most critical amino acids (CAA) in protein <bold>(C)</bold> for samples collected from the United States, Brazil, and the Southern Cone (Argentina, Paraguay, and Uruguay) regions across the Western Hemisphere. A total of 1817 samples were collected during the 2012, 2013, and 2014 cropping seasons. Abundance of CAA was calculated as [(lysine&#x2009;+&#x2009;threonine&#x2009;+&#x2009;tryptophane&#x2009;+&#x2009;cysteine&#x2009;+&#x2009;methionine)/net protein] x&#x2009;100.</p>
</caption>
<graphic xlink:href="fsufs-07-1223921-g002.tif"/>
</fig>
<fig position="float" id="fig3">
<label>Figure 3</label>
<caption>
<p>Mean and distribution for soybean crude protein concentration and the abundance of the most critical amino acids (CAA) in protein for different subregions within the United States, Brazil, and the Southern Cone based on on-farm samples collected between 2012 and 2014 cropping seasons. The subregions are as follows: United States (red)- Northern Corn Belt (nc), Eastern Corn Belt (ec), Western Corn Belt (wc), Northeast (ne), and South (su); Brazil (green)- Cerrado (ce), Southern Brazil (sb), Atlantic Forest (af), and Northern Brazil (nb); Southern Cone (blue)- Northern Pampas (np), Central Pampas (cp), Southern Pampas (sp), Paraguay (py), and Uruguay (uy). The gray strip contains the data between the 25th to the 75th quantile (50% of the data). Abundance of CAA was calculated as [(lysine&#x2009;+&#x2009;threonine&#x2009;+&#x2009;tryptophane&#x2009;+&#x2009;cysteine&#x2009;+&#x2009;methionine)/net protein] x&#x2009;100.</p>
</caption>
<graphic xlink:href="fsufs-07-1223921-g003.tif"/>
</fig>
</sec>
<sec id="sec9">
<label>3.2.</label>
<title>Amino acid profile &#x2013; spatial pattern</title>
<p>The relative abundance of the most limiting amino acids in protein was higher for the Southern Cone, intermediate for the United States, and lower for Brazil (<xref rid="fig2" ref-type="fig">Figure 2</xref>). The subregion and year interaction effect also affected this fraction, but as with protein concentrations, this effect explained a small proportion (less than 5%) of the total variation (<xref rid="tab1" ref-type="table">Table 1</xref> and <xref rid="fig3" ref-type="fig">Figure 3</xref>). For example, Uruguay in the Southern Cone during the 2012 and 2014 growing seasons showed a clear effect of the subregion and year interaction (<xref rid="fig3" ref-type="fig">Figure 3</xref>). When considering the most limiting amino acids individually, their abundance in protein was not consistently affected by the region or the subregion as the variance components for region and subregion generally captured a small proportion of the variation for each amino acid (<xref rid="tab1" ref-type="table">Table 1</xref>). Year was an important source of variation for the relative amount of some relevant amino acids, such as threonine and tryptophane. The interaction between region and year explained almost a third and a fifth of the variation in the abundance of lysine and tryptophane in protein, respectively (<xref rid="tab1" ref-type="table">Table 1</xref>). The effect of the subregion and year interaction on the abundance of individual amino acids in protein was mostly significant, but it explained a small proportion (less than 5%) of the total variation (<xref rid="tab1" ref-type="table">Table 1</xref>).</p>
</sec>
<sec id="sec10">
<label>3.3.</label>
<title>Amino acid profile &#x2013; relationships among grain constituents</title>
<p>Principal component analysis allowed us to explore relationships among amino acid abundance in protein and other grain constituents (<xref rid="fig4" ref-type="fig">Figure 4</xref>). From the principal component analysis performed, principal components (PC) 1 and 2 explained 52% of the total variation (<xref rid="fig4" ref-type="fig">Figure 4</xref>), and with the inclusion of PC3 and PC4, the explained variability increased to 72% (data not shown). PC1 indicated that the abundance of most critical amino acids in protein was negatively associated with protein concentration. To evaluate this relationship, we performed separate linear regression analyzes for each region, utilizing bootstrapping with reposition to compare slopes and intercepts. Our results indicated overlapping confidence intervals across all three regions (data not shown), thus justifying the use of a unique linear regression analysis to assess the relationship between these grain components (<xref rid="fig5" ref-type="fig">Figure 5</xref>). Abundance of most limiting amino acids decreased 0.7% for each percentage point increase in protein (<italic>R</italic><sup>2</sup>&#x2009;=&#x2009;0.26, <italic>p</italic>&#x2009;&#x003C;&#x2009;0.001). Principal component analysis also showed that when grouping amino acids by their biosynthetic pathway, alpha-ketoglutarate derived amino acids (arginine, glutamine, proline, histidine) abundance was negatively associated with that for oxaloacetate derived ones (lysine, methionine, threonine, isoleucine), and positively associated with protein concentration (<xref rid="fig4" ref-type="fig">Figure 4</xref>). The pool of alpha-ketoglutarate derived amino acids increased proportionally more than the pool of oxaloacetate derived amino acids as crude protein increased (<xref rid="fig6" ref-type="fig">Figure 6</xref>). In accordance, lysine and arginine were the most responsive to protein variation in the oxaloacetate and the alpha-ketoglutarate pathways, respectively. Lysine abundance in protein was negatively (<italic>r</italic>&#x2009;=&#x2009;&#x2212;0.45, <italic>p</italic>&#x2009;&#x003C;&#x2009;0.001) and arginine was positively (<italic>r</italic>&#x2009;=&#x2009;0.58, <italic>p</italic>&#x2009;&#x003C;&#x2009;0.001) correlated with crude protein concentration (data not shown). Accordingly, the ratio lysine/arginine decreased sharply with the increase in protein concentration (<xref ref-type="supplementary-material" rid="SM1">Supplementary Figure S1</xref>). Oil and carbohydrate concentrations did not help to explain the amino acid profile (<xref rid="fig4" ref-type="fig">Figure 4</xref>).</p>
<fig position="float" id="fig4">
<label>Figure 4</label>
<caption>
<p>Biplot derived from principal component analysis (PCA) for grain compositional variables for soybean, including crude protein, net protein, oil, carbohydrates (CHO), and non-protein N, expressed as a fraction of the grain. The amino acids lysine, threonine, tryptophane, cysteine + methionine, arginine, and the sum of the most critical amino acid (CAA&#x2009;=&#x2009;lysine + threonine + tryptophane + cysteine + methionine) are expressed as a fraction of net protein, were also included. The color of the arrows indicates the contribution of each variable to the first two dimensions of the PCA. The data used for this biplot analysis consisted of 1,817 samples collected across the Western Hemisphere between 2012 and 2014.</p>
</caption>
<graphic xlink:href="fsufs-07-1223921-g004.tif"/>
</fig>
<fig position="float" id="fig5">
<label>Figure 5</label>
<caption>
<p>Abundance of the most critical amino acids (CAA) in protein as a function of protein concentration in grain. Abundance of CAA was calculated as [(lysine&#x2009;+&#x2009;threonine&#x2009;+&#x2009;tryptophane&#x2009;+&#x2009;cysteine&#x2009;+&#x2009;methionine)/net protein] x&#x2009;100. The linear regression is presented by the red dotted line (<italic>R</italic><sup>2</sup>&#x2009;=&#x2009;0.27; <italic>y</italic>&#x2009;=&#x2009;0.18 &#x2013; x0.00096; <italic>p</italic>&#x2009;&#x003C;&#x2009;0.001). The black dotted lines connect points on the graph representing equal amounts of CAA per ton of grain.</p>
</caption>
<graphic xlink:href="fsufs-07-1223921-g005.tif"/>
</fig>
<fig position="float" id="fig6">
<label>Figure 6</label>
<caption>
<p>Abundance of oxaloacetate-derived amino acids (lysine, methionine, threonine, isoleucine, asparagine) as a function of the abundance of ketoglutarate-derived amino acids (arginine, glutamate, proline, histidine) in soybean grain protein. Color gradient indicates grain crude protein concentration (%). The linear regression is presented by the dotted line (<italic>R</italic><sup>2</sup>&#x2009;=&#x2009;0.55, <italic>y</italic>&#x2009;=&#x2009;44&#x2013;0.47x, <italic>p</italic>&#x2009;&#x003C;&#x2009;0.001), with the 95% confidence interval shaded in gray. The data is based on 1,817 samples collected across the Western Hemisphere between 2012&#x2013;2014. &#x002A;Glutamate is glutamate + glutamine.</p>
</caption>
<graphic xlink:href="fsufs-07-1223921-g006.tif"/>
</fig>
</sec>
<sec id="sec11">
<label>3.4.</label>
<title>Protein and amino acid production capacity</title>
<p>Combining protein concentration and amino acid profile data with the soybean production projected by 2030, we have estimated that the Western Hemisphere&#x2019;s protein production capacity could reach approximately 120 million metric tons annually. Furthermore, our analysis suggests that the combined production of the most critical amino acids could approach 17 million metric tons per year (<xref rid="tab2" ref-type="table">Table 2</xref>). Projected protein and critical amino acids production was also estimated at region (<xref rid="tab2" ref-type="table">Table 2</xref>) and subregion (<xref rid="fig7" ref-type="fig">Figure 7</xref> and <xref ref-type="supplementary-material" rid="SM1">Supplementary Table S3</xref>) levels.</p>
<table-wrap position="float" id="tab2">
<label>Table 2</label>
<caption>
<p>Projected annually production of soybean, grain crude protein, and the sum of the most critical amino acids (CAA&#x2009;=&#x2009;lysine + threonine + tryptophan + cysteine + methionine) in the Western Hemisphere by 2030 expressed in megatons (Mt).</p>
</caption>
<table frame="hsides" rules="groups">
<thead>
<tr>
<th align="left" valign="middle">Region</th>
<th align="center" valign="middle">Grain (Mt)</th>
<th align="center" valign="middle">Crude protein (Mt)</th>
<th align="center" valign="middle">CAA (Mt)</th>
</tr>
</thead>
<tbody>
<tr>
<td align="left" valign="middle">United States</td>
<td align="char" valign="middle" char=".">132.1</td>
<td align="char" valign="middle" char=".">43.8</td>
<td align="char" valign="middle" char=".">7.40</td>
</tr>
<tr>
<td align="left" valign="middle">Brazil</td>
<td align="char" valign="middle" char=".">149.3</td>
<td align="char" valign="middle" char=".">50.9</td>
<td align="char" valign="middle" char=".">3.37</td>
</tr>
<tr>
<td align="left" valign="middle">Southern Cone</td>
<td align="char" valign="middle" char=".">70.9</td>
<td align="char" valign="middle" char=".">22.8</td>
<td align="char" valign="middle" char=".">6.47</td>
</tr>
<tr>
<td align="left" valign="middle">Total</td>
<td align="char" valign="middle" char=".">352.4</td>
<td align="char" valign="middle" char=".">117.5</td>
<td align="char" valign="middle" char=".">17.24</td>
</tr>
</tbody>
</table>
<table-wrap-foot>
<p>Projection based on measured grain composition and soybean grain production projections derived from FAO (<xref ref-type="bibr" rid="ref26">OECD and FAO, 2022</xref>).</p>
</table-wrap-foot>
</table-wrap>
<fig position="float" id="fig7">
<label>Figure 7</label>
<caption>
<p>Projected crude protein <bold>(A)</bold> and most critical amino acids (lysine, threonine, tryptophan, cysteine, and methionine); <bold>(B)</bold> production for different cropping subregions across the Western Hemisphere by the year 2030. The projections are based on the presented grain composition analysis and projections for soybean grain production at a country level obtained from <xref ref-type="bibr" rid="ref26">OECD and FAO (2022)</xref>.</p>
</caption>
<graphic xlink:href="fsufs-07-1223921-g007.tif"/>
</fig>
</sec>
</sec>
<sec sec-type="discussions" id="sec12">
<label>4.</label>
<title>Discussion</title>
<p>Major soybean cropping regions and subregions were sampled according to their contribution to total grain production. This original method allowed us to create realistic comparisons of grain composition. In this way, the study stands out as the first comprehensive examination of protein and amino acid profiles in soybeans across the entire Americas, surpassing previous research that primarily focused on individual country-level analyzes. Furthermore, this research pioneered the exploration of novel relationships among the most relevant amino acids and protein for soybean production systems, shedding new light on their complex interactions. Additionally, our study breaks new ground by providing valuable insights into the continent&#x2019;s protein production capacity, an essential component of global soybean production, which has not been addressed in prior literature. This crucial perspective offers unique implications for the agricultural industry and food security on a global scale.</p>
<p>Our study revealed a consistent spatial pattern for crude protein concentration in soybeans across the United States, Brazil, and the Southern Cone soybean-producing regions. Brazil produced 2.0 and 5.4% more protein per ton of grain than the United States and the Southern Cone, respectively (<xref rid="fig2" ref-type="fig">Figure 2</xref>). Within the Southern Cone, we observed notable differences in protein concentrations among subregions, with a consistently higher crude protein concentration in the Northern Pampas, Paraguay, and Uruguay compared to the Southern and Central Pampas (<xref rid="fig3" ref-type="fig">Figure 3</xref>). In contrast, a relatively lower variability among subregions was found within the United States and Brazil.</p>
<p>It&#x2019;s worth noting that the market perceived soybeans from the U.S. Northern Cornbelt to have lower protein concentrations, leading to trading discounts for soybeans shipped at Pacific Northwest (<xref ref-type="bibr" rid="ref15">Hertsgaard et al., 2019</xref>; <xref ref-type="bibr" rid="ref49">William et al., 2020</xref>). This concept is also supported by an older study by <xref ref-type="bibr" rid="ref16">Hurburgh et al. (1990)</xref> which found that protein for soybean from the northern states were 1.5 to 2 percentage points lower compared to that from central states. However, from our dataset the effect was not consistent across years, and on average, soybeans produced in the Northern Corn Belt region had only 0.5 percentage points lower protein concentration than the other regions examined (<xref rid="tab1" ref-type="table">Table 1</xref> and <xref rid="fig3" ref-type="fig">Figure 3</xref>). This last difference aligns with more recent studies on soybean protein patterns across the United States (<xref ref-type="bibr" rid="ref38">Rotundo et al., 2016</xref>; <xref ref-type="bibr" rid="ref3">Assefa et al., 2019</xref>).</p>
<p>Crude protein content of commercial soybeans is typically determined by nitrogen combustion and multiplying the result by a conversion factor of 6.25 (<xref ref-type="bibr" rid="ref21">Krul, 2019</xref>). Our analysis showed that this method overestimated the net protein content by an average of 3%, due to the presence of non-protein nitrogen fractions (<xref rid="fig2" ref-type="fig">Figure 2</xref> and <xref ref-type="supplementary-material" rid="SM1">Supplementary Table S1</xref>). This finding is consistent with previous studies that were based on a more limited number of samples (<xref ref-type="bibr" rid="ref20">Krober and Gibbons, 1962</xref>; <xref ref-type="bibr" rid="ref45">Thakur and Hurburgh, 2007</xref>). Relaying solely on crude protein content can be misleading when assessing soybean meal quality, especially for monogastric animals that cannot utilize non-protein nitrogen fractions. Additional methods for evaluating the quality of soybean meal should be considered.</p>
<p>The protein amino acid profile in soybean exhibited consistent differences across regions, as shown in <xref rid="tab1" ref-type="table">Table 1</xref> and <xref rid="fig3" ref-type="fig">Figure 3</xref>. Among the assessed samples, the abundance of the most limiting amino acids varied up to 8.5% (<xref rid="fig3" ref-type="fig">Figures 3</xref>, <xref rid="fig5" ref-type="fig">5</xref>), with the Southern Cone having the highest abundance. The ranking among subregions was season-dependent, but averaged across years, the Central and Southern Pampas subregions showed the highest limiting amino acids abundance (<xref rid="fig3" ref-type="fig">Figure 3</xref>). Limiting amino acids concentration variability was also highest for the Southern Cone (<xref rid="fig3" ref-type="fig">Figures 3</xref>, <xref rid="fig4" ref-type="fig">4</xref>), suggesting that soybeans from the United States and Brazil may be more consistent and predictable in their amino acid profile. Understanding and predicting this variability could be relevant to adequate production according to market or industry requirements.</p>
<p>While being aware of the potential correlations between subregions and regions, our initial analysis, as presented in <xref rid="fig5" ref-type="fig">Figure 5</xref>, provided a valuable starting point for understanding relevant relationships between grain components. Critical amino acids abundance in protein was reduced when protein content increased. This dilution effect has been previously reported (<xref ref-type="bibr" rid="ref45">Thakur and Hurburgh, 2007</xref>; <xref ref-type="bibr" rid="ref24">Medic et al., 2014</xref>; <xref ref-type="bibr" rid="ref29">Pfarr et al., 2018</xref>) but never using such a wide spectrum of samples. The implications of this dilution process on the nutritional value of the protein (<xref ref-type="bibr" rid="ref44">Tessari et al., 2016</xref>) are significant and warrant further investigation. A higher abundance of the most limiting amino acids in protein would lead to a higher conversion efficiency of protein in meat (kg meat/kg protein; <xref ref-type="bibr" rid="ref40">Selle et al., 2020</xref>). Soybean meal made from low protein soybeans produced in the Central and Southern Pampas would likely have a higher efficiency to convert protein in meat than a meal made from high protein soybeans produced in Brazil or the United States, which could be relevant from a sustainability standpoint (<xref ref-type="bibr" rid="ref31">Pomar et al., 2021</xref>). Conversely our data showed that increasing protein concentration, though more diluted, also increased the amount of the most critical amino acids per ton of grain (<xref rid="fig5" ref-type="fig">Figure 5</xref>). Thus, meals made from higher protein concentration soybeans, like those produced in Brazil or the United States, would likely present a higher conversion efficiency of meal in meat (kg meat/kg meal) than meals made from lower protein concentration soybeans like those from the Central and Southern Pampas.</p>
<p>Our data revealed a strong negative association (<italic>R</italic><sup>2</sup> =&#x2009;&#x2212;0.55, <italic>p</italic> &#x003C;&#x2009;0.001) between oxaloacetate path-derived amino acids (lysine, methionine, threonine, isoleucine, asparagine) and alpha-ketoglutarate path-derived amino acids (arginine, glutamine, proline, histidine) in protein (<xref rid="fig6" ref-type="fig">Figure 6</xref>). Increasing protein concentration resulted in a greater increase in alpha-ketoglutarate pathway-derived amino acids compared to amino acids derived from the aspartate pathway. This association was primarily driven by an increase in arginine abundance and a decrease in lysine (<xref ref-type="supplementary-material" rid="SM1">Supplementary Figure S1</xref>) and is consistent with results from controlled experiments on soybean seed protein concentration alteration (<xref ref-type="bibr" rid="ref29">Pfarr et al., 2018</xref>).</p>
<p>Arginine, with the highest N content per mass unit (32%), is a major nitrogen storage and transport molecule in plants (<xref ref-type="bibr" rid="ref50">Winter et al., 2015</xref>). Conditions of abundant reduced N may promote the synthesis of this amino acid, which can store large quantities of N per unit of protein to support future N requirements of seedling. What is relevant from a feed end-use standpoint is that using the reduced nitrogen to produce a gram of arginine would allow the plant to synthesize almost twice the amount of lysine and tryptophane, three times the threonine and cysteine, and four times more methionine (<xref ref-type="bibr" rid="ref46">Tillman, 2019</xref>). Understanding the physiological basis, limits, and control of these relationships among amino acids through breeding, genetic engineering, or crop management could provide a via to increase soybean protein value by adapting the amino acid profile according to its end-use.</p>
<p>This study represents a significant step towards a comprehensive understanding of the variation in soybean seed composition at a regional scale in the Western Hemisphere. However, to achieve a more complete and accurate prediction of soybean composition, further research is required to investigate the physiological processes involved. While recent studies have shed light on the factors that impact soybean seed protein concentration (<xref ref-type="bibr" rid="ref30">Piper and Boote, 1999</xref>; <xref ref-type="bibr" rid="ref9">Carrera et al., 2009</xref>; <xref ref-type="bibr" rid="ref36">Rotundo and Westgate, 2009</xref>; <xref ref-type="bibr" rid="ref1">Anthony et al., 2012</xref>; <xref ref-type="bibr" rid="ref29">Pfarr et al., 2018</xref>; <xref ref-type="bibr" rid="ref7">Bosaz et al., 2019</xref>), more research is needed to fully elucidate the control of the amino acid profile. Nonetheless, the extensive and high-quality dataset of grain composition presented in this study (<xref ref-type="supplementary-material" rid="SM1">Supplementary Table S2</xref>) can be a valuable resource for validating emerging conceptual frameworks and modeling tools in this field.</p>
<p>By 2030, soybean production in the Western Hemisphere is expected to reach 352 Mt, accounting for 86% of global soybean production (<xref ref-type="bibr" rid="ref26">OECD and FAO, 2022</xref>). Based on this forecast, we projected the protein production capacity of the region to be approximately 120 Mt, containing only 17 Mt. of the most limiting amino acids for monogastric diets. While our estimations are subject to uncertainties, we used the best available data to estimate the protein and amino acid production capacity of the Western Hemisphere. Brazil is projected to be the most important region for soybean protein production (<xref rid="tab2" ref-type="table">Table 2</xref>), driven by high soybean grain production and protein concentration (<xref rid="fig3" ref-type="fig">Figure 3</xref>). Notably, the Cerrado subregion alone could produce almost 20% of the soybean protein in the Western Hemisphere. These findings have important implications for addressing global challenges related to producing enough protein for feed and food protein markets (<xref ref-type="bibr" rid="ref18">Kim et al., 2019</xref>; <xref ref-type="bibr" rid="ref14">Guthman et al., 2022</xref>).</p>
<p>Finally, while our study focused on analyzing the protein concentration and amino acid profile of soybeans in the Western Hemisphere, it&#x2019;s important to note that the complete nutritional and economic value of soybean grain includes not only the protein fraction but also the carbohydrate and the oil fractions (<xref ref-type="bibr" rid="ref24">Medic et al., 2014</xref>). Therefore, to fully assess the value of soybean grain produced in different regions or its potential applications in various industries, it&#x2019;s essential to consider not only the protein but also the oil and carbohydrates fractions of the grain.</p>
</sec>
<sec sec-type="conclusions" id="sec13">
<label>5.</label>
<title>Conclusion</title>
<p>This study presented a novel method for assessing soybean grain composition across the Western Hemisphere. Spatial processes played a relevant role in determining protein concentration, which exhibit a consistent pattern across the major soybean-producing regions. The analysis also uncovered differences in the amino acid profile across regions and established relationships among relevant amino acids which could be relevant for markets and industries. The generated dataset should be relevant for researchers seeking to understand processes involved in grain composition determination. Lastly, the study projected global soybean protein and amino acid production capacities for this major region by 2030, which is highly relevant to analyzing the global challenges of producing enough food. Overall, this research provides valuable information for industry, markets, and researchers useful as a benchmark for soybean grain composition across the Western Hemisphere.</p>
</sec>
<sec sec-type="data-availability" id="sec14">
<title>Data availability statement</title>
<p>The original contributions presented in the study are included in the article/<xref ref-type="supplementary-material" rid="SM1">Supplementary material</xref>, further inquiries can be directed to the corresponding author.</p>
</sec>
<sec sec-type="author-contributions" id="sec15">
<title>Author contributions</title>
<p>AC: conceptualization, data curation, formal analysis, investigation, methodology, visualization, and writing &#x2013; original draft. JM-G: writing &#x2013; review &#x0026; editing. SN: conceptualization, investigation, methodology, funding acquisition, and writing &#x2013; review &#x0026; editing. All authors contributed to the article and approved the submitted version.</p>
</sec>
</body>
<back>
<sec sec-type="funding-information" id="sec16">
<title>Funding</title>
<p>The United Soybean Board (USB) provided funding to support the conducted analysis.</p>
</sec>
<sec sec-type="COI-statement" id="sec17">
<title>Conflict of interest</title>
<p>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.</p>
</sec>
<sec id="sec100" sec-type="disclaimer">
<title>Publisher&#x2019;s note</title>
<p>All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.</p>
</sec>
<sec sec-type="supplementary-material" id="sec18">
<title>Supplementary material</title>
<p>The Supplementary material for this article can be found online at: <ext-link xlink:href="https://www.frontiersin.org/articles/10.3389/fsufs.2023.1223921/full#supplementary-material" ext-link-type="uri">https://www.frontiersin.org/articles/10.3389/fsufs.2023.1223921/full#supplementary-material</ext-link></p>
<supplementary-material xlink:href="Table_1.docx" id="SM1" mimetype="application/vnd.openxmlformats-officedocument.wordprocessingml.document" xmlns:xlink="http://www.w3.org/1999/xlink"/>
</sec>
<fn-group>
<fn id="fn0001">
<p>
<sup>1</sup>
<ext-link xlink:href="https://www.nass.usda.gov/" ext-link-type="uri">https://www.nass.usda.gov/</ext-link>
</p>
</fn>
<fn id="fn0002">
<p>
<sup>2</sup>
<ext-link xlink:href="https://www.ibge.gov.br/" ext-link-type="uri">https://www.ibge.gov.br/</ext-link>
</p>
</fn>
<fn id="fn0003">
<p>
<sup>3</sup>
<ext-link xlink:href="https://www.magyp.gob.ar/sitio/areas/estimaciones/" ext-link-type="uri">https://www.magyp.gob.ar/sitio/areas/estimaciones/</ext-link>
</p>
</fn>
<fn id="fn0004">
<p>
<sup>4</sup>
<ext-link xlink:href="https://capeco.org.py/" ext-link-type="uri">https://capeco.org.py/</ext-link>
</p>
</fn>
</fn-group>
<ref-list>
<title>References</title>
<ref id="ref1">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Anthony</surname> <given-names>P.</given-names></name> <name><surname>Malzer</surname> <given-names>G.</given-names></name> <name><surname>Sparrow</surname> <given-names>S.</given-names></name> <name><surname>Zhang</surname> <given-names>M.</given-names></name></person-group> (<year>2012</year>). <article-title>Soybean yield and quality in relation to soil properties</article-title>. <source>Agron. J.</source> <volume>104</volume>, <fpage>1443</fpage>&#x2013;<lpage>1458</lpage>. doi: <pub-id pub-id-type="doi">10.2134/agronj2012.0095</pub-id></citation>
</ref>
<ref id="ref2">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Assefa</surname> <given-names>Y.</given-names></name> <name><surname>Bajjalieh</surname> <given-names>N.</given-names></name> <name><surname>Archontoulis</surname> <given-names>S.</given-names></name> <name><surname>Casteel</surname> <given-names>S.</given-names></name> <name><surname>Davidson</surname> <given-names>D.</given-names></name> <name><surname>Kov&#x00E1;cs</surname> <given-names>P.</given-names></name> <etal/></person-group>. (<year>2018</year>). <article-title>Spatial characterization of soybean yield and quality (amino acids, oil, and protein) for United States</article-title>. <source>Sci. Rep.</source> <volume>8</volume>:<fpage>14653</fpage>. doi: <pub-id pub-id-type="doi">10.1038/s41598-018-32895-0</pub-id>, PMID: <pub-id pub-id-type="pmid">30279447</pub-id></citation>
</ref>
<ref id="ref3">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Assefa</surname> <given-names>Y.</given-names></name> <name><surname>Purcell</surname> <given-names>L. C.</given-names></name> <name><surname>Salmeron</surname> <given-names>M.</given-names></name> <name><surname>Naeve</surname> <given-names>S.</given-names></name> <name><surname>Casteel</surname> <given-names>S. N.</given-names></name> <name><surname>Kov&#x00E1;cs</surname> <given-names>P.</given-names></name> <etal/></person-group>. (<year>2019</year>). <article-title>Assessing variation in US soybean seed composition (protein and oil)</article-title>. <source>Front. Plant Sci.</source> <volume>10</volume>:<fpage>298</fpage>. doi: <pub-id pub-id-type="doi">10.3389/fpls.2019.00298</pub-id>, PMID: <pub-id pub-id-type="pmid">30915097</pub-id></citation>
</ref>
<ref id="ref4">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Attia</surname> <given-names>Y. A.</given-names></name> <name><surname>Bovera</surname> <given-names>F.</given-names></name> <name><surname>Wang</surname> <given-names>J.</given-names></name> <name><surname>Al-Harthi</surname> <given-names>M. A.</given-names></name> <name><surname>Kim</surname> <given-names>W. K.</given-names></name></person-group> (<year>2020</year>). <article-title>Multiple amino acid supplementations to low-protein diets: effect on performance, carcass yield, meat quality and nitrogen excretion of finishing broilers under hot climate conditions</article-title>. <source>Animals</source> <volume>10</volume>:<fpage>973</fpage>. doi: <pub-id pub-id-type="doi">10.3390/ani10060973</pub-id>, PMID: <pub-id pub-id-type="pmid">32503244</pub-id></citation>
</ref>
<ref id="ref5">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Bakare</surname> <given-names>M. A.</given-names></name> <name><surname>Kayondo</surname> <given-names>S. I.</given-names></name> <name><surname>Aghogho</surname> <given-names>C. I.</given-names></name> <name><surname>Wolfe</surname> <given-names>M. D.</given-names></name> <name><surname>Parkes</surname> <given-names>E. Y.</given-names></name> <name><surname>Kulakow</surname> <given-names>P.</given-names></name> <etal/></person-group>. (<year>2022</year>). <article-title>Exploring genotype by environment interaction on cassava yield and yield related traits using classical statistical methods</article-title>. <source>PLoS One</source> <volume>17</volume>:<fpage>e0268189</fpage>. doi: <pub-id pub-id-type="doi">10.1371/journal.pone.0268189</pub-id>, PMID: <pub-id pub-id-type="pmid">35849556</pub-id></citation>
</ref>
<ref id="ref6">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Ball</surname> <given-names>R. O.</given-names></name> <name><surname>Courtney-Martin</surname> <given-names>G.</given-names></name> <name><surname>Pencharz</surname> <given-names>P. B.</given-names></name></person-group> (<year>2006</year>). <article-title>The in vivo sparing of methionine by lysteine in sulfur amino acid requirements in animal models and adult humans</article-title>. <source>J. Nutr.</source> <volume>136</volume>, <fpage>1682S</fpage>&#x2013;<lpage>1693S</lpage>. doi: <pub-id pub-id-type="doi">10.1093/jn/136.6.1682S</pub-id>, PMID: <pub-id pub-id-type="pmid">16702340</pub-id></citation>
</ref>
<ref id="ref7">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Bosaz</surname> <given-names>L. B.</given-names></name> <name><surname>Gerde</surname> <given-names>J. A.</given-names></name> <name><surname>Borr&#x00E1;s</surname> <given-names>L.</given-names></name> <name><surname>Cipriotti</surname> <given-names>P. A.</given-names></name> <name><surname>Ascheri</surname> <given-names>L.</given-names></name> <name><surname>Campos</surname> <given-names>M.</given-names></name> <etal/></person-group>. (<year>2019</year>). <article-title>Management and environmental factors explaining soybean seed protein variability in Central Argentina</article-title>. <source>Field Crop Res.</source> <volume>240</volume>, <fpage>34</fpage>&#x2013;<lpage>43</lpage>. doi: <pub-id pub-id-type="doi">10.1016/j.fcr.2019.05.007</pub-id></citation>
</ref>
<ref id="ref8">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Breene</surname> <given-names>W. M.</given-names></name> <name><surname>Lin</surname> <given-names>S.</given-names></name> <name><surname>Hardman</surname> <given-names>L.</given-names></name> <name><surname>Orf</surname> <given-names>J.</given-names></name></person-group> (<year>1988</year>). <article-title>Protein and oil content of soybeans from different geographic locations</article-title>. <source>J. Am. Oil Chem. Soc.</source> <volume>65</volume>, <fpage>1927</fpage>&#x2013;<lpage>1931</lpage>. doi: <pub-id pub-id-type="doi">10.1007/BF02546009</pub-id></citation>
</ref>
<ref id="ref9">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Carrera</surname> <given-names>C.</given-names></name> <name><surname>Mart&#x00ED;nez</surname> <given-names>M. J.</given-names></name> <name><surname>Dardanelli</surname> <given-names>J.</given-names></name> <name><surname>Balzarini</surname> <given-names>M.</given-names></name></person-group> (<year>2009</year>). <article-title>Water deficit effect on the relationship between temperature during the seed fill period and soybean seed oil and protein concentrations</article-title>. <source>Crop Sci.</source> <volume>49</volume>, <fpage>990</fpage>&#x2013;<lpage>998</lpage>. doi: <pub-id pub-id-type="doi">10.2135/cropsci2008.06.0361</pub-id></citation>
</ref>
<ref id="ref10">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Carrera</surname> <given-names>C. S.</given-names></name> <name><surname>Rosas</surname> <given-names>M. B.</given-names></name> <name><surname>Gontijo Mandarino</surname> <given-names>J. M.</given-names></name> <name><surname>Leite</surname> <given-names>R. S.</given-names></name> <name><surname>Raspa</surname> <given-names>F.</given-names></name> <name><surname>Fava</surname> <given-names>F.</given-names></name> <etal/></person-group>. (<year>2022</year>). <article-title>Partial and total defoliation during the filling period affected grain industrial and nutraceutical quality in soybean</article-title>. <source>J. Sci. Food Agric.</source> <volume>102</volume>, <fpage>4109</fpage>&#x2013;<lpage>4120</lpage>. doi: <pub-id pub-id-type="doi">10.1002/jsfa.11760</pub-id>, PMID: <pub-id pub-id-type="pmid">34997583</pub-id></citation>
</ref>
<ref id="ref11">
<citation citation-type="other"><person-group person-group-type="author">
<collab id="coll1">FAO Faostat</collab>
</person-group>. (<year>2022</year>). Statistics database of the food and agriculture organization of the United Nations. Available at: <ext-link xlink:href="http://www.fao.org/statistics/databases/en/" ext-link-type="uri">http://www.fao.org/statistics/databases/en/</ext-link></citation>
</ref>
<ref id="ref12">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Godfray</surname> <given-names>H. C. J.</given-names></name> <name><surname>Garnett</surname> <given-names>T.</given-names></name></person-group> (<year>2014</year>). <article-title>Food security and sustainable intensification</article-title>. <source>Philos. Trans. R. Soc. Lond. Ser. B Biol. Sci.</source> <volume>369</volume>:<fpage>20120273</fpage>. doi: <pub-id pub-id-type="doi">10.1098/rstb.2012.0273</pub-id>, PMID: <pub-id pub-id-type="pmid">24535385</pub-id></citation>
</ref>
<ref id="ref13">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Goldsmith</surname> <given-names>P. D.</given-names></name> <name><surname>Gunjal</surname> <given-names>K.</given-names></name> <name><surname>Ndarishikanye</surname> <given-names>B.</given-names></name></person-group> (<year>2004</year>). <article-title>Rural&#x2013;urban migration and agricultural productivity: the case of Senegal</article-title>. <source>Agric. Econ.</source> <volume>31</volume>, <fpage>33</fpage>&#x2013;<lpage>45</lpage>. doi: <pub-id pub-id-type="doi">10.1111/j.1574-0862.2004.tb00220.x</pub-id></citation>
</ref>
<ref id="ref14">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Guthman</surname> <given-names>J.</given-names></name> <name><surname>Butler</surname> <given-names>M.</given-names></name> <name><surname>Martin</surname> <given-names>S. J.</given-names></name></person-group> (<year>2022</year>). <article-title>In the name of protein</article-title>. <source>Nat. Food</source> <volume>3</volume>, <fpage>391</fpage>&#x2013;<lpage>393</lpage>. doi: <pub-id pub-id-type="doi">10.1038/s43016-022-00532-9</pub-id></citation>
</ref>
<ref id="ref15">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Hertsgaard</surname> <given-names>D. J.</given-names></name> <name><surname>Wilson</surname> <given-names>W. W.</given-names></name> <name><surname>Dahl</surname> <given-names>B.</given-names></name></person-group> (<year>2019</year>). <article-title>Costs and risks of testing and blending for essential amino acids in soybeans</article-title>. <source>Agribusiness</source> <volume>35</volume>, <fpage>265</fpage>&#x2013;<lpage>280</lpage>. doi: <pub-id pub-id-type="doi">10.1002/agr.21576</pub-id></citation>
</ref>
<ref id="ref16">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Hurburgh</surname> <given-names>C. R.</given-names> <suffix>Jr.</suffix></name> <name><surname>Brumm</surname> <given-names>T. J.</given-names></name> <name><surname>Guinn</surname> <given-names>J. M.</given-names></name> <name><surname>Hartwig</surname> <given-names>R. A.</given-names></name></person-group> (<year>1990</year>). <article-title>Protein and oil patterns in US and world soybean markets</article-title>. <source>J. Am. Oil Chem. Soc.</source> <volume>67</volume>, <fpage>966</fpage>&#x2013;<lpage>973</lpage>. doi: <pub-id pub-id-type="doi">10.1007/BF02541859</pub-id></citation>
</ref>
<ref id="ref17">
<citation citation-type="other"><person-group person-group-type="author"><name><surname>Ju</surname> <given-names>J.</given-names></name> <name><surname>Banfelder</surname> <given-names>J.</given-names></name> <name><surname>Skrabanek</surname> <given-names>L.</given-names></name></person-group> (<year>2019</year>). Quantitative understanding in biology principal component analysis. Available at: <ext-link xlink:href="https://physiology.med.cornell.edu/people/banfelder/qbio/lecture_notes/3.4_Principal_component_analysis.pdf" ext-link-type="uri">https://physiology.med.cornell.edu/people/banfelder/qbio/lecture_notes/3.4_Principal_component_analysis.pdf</ext-link></citation>
</ref>
<ref id="ref18">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Kim</surname> <given-names>S. W.</given-names></name> <name><surname>Less</surname> <given-names>J. F.</given-names></name> <name><surname>Wang</surname> <given-names>L.</given-names></name> <name><surname>Yan</surname> <given-names>T.</given-names></name> <name><surname>Kiron</surname> <given-names>V.</given-names></name> <name><surname>Kaushik</surname> <given-names>S. J.</given-names></name> <etal/></person-group>. (<year>2019</year>). <article-title>Meeting global feed protein demand: challenge, opportunity, and strategy</article-title>. <source>Annu. Rev. Anim. Biosci.</source> <volume>7</volume>, <fpage>221</fpage>&#x2013;<lpage>243</lpage>. doi: <pub-id pub-id-type="doi">10.1146/annurev-animal-030117-014838</pub-id>, PMID: <pub-id pub-id-type="pmid">30418803</pub-id></citation>
</ref>
<ref id="ref19">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Kovalenko</surname> <given-names>I. V.</given-names></name> <name><surname>Rippke</surname> <given-names>G. R.</given-names></name> <name><surname>Hurburgh</surname> <given-names>C. R.</given-names></name></person-group> (<year>2006</year>). <article-title>Determination of amino acid composition of soybeans (<italic>Glycine max</italic>) by near-infrared spectroscopy</article-title>. <source>J. Agric. Food Chem.</source> <volume>54</volume>, <fpage>3485</fpage>&#x2013;<lpage>3491</lpage>. doi: <pub-id pub-id-type="doi">10.1021/jf052570u</pub-id>, PMID: <pub-id pub-id-type="pmid">19127714</pub-id></citation>
</ref>
<ref id="ref20">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Krober</surname> <given-names>O. A.</given-names></name> <name><surname>Gibbons</surname> <given-names>S. J.</given-names></name></person-group> (<year>1962</year>). <article-title>Composition of feedstuffs, nonprotein nitrogen in soybeans</article-title>. <source>J. Agric. Food Chem.</source> <volume>10</volume>, <fpage>57</fpage>&#x2013;<lpage>59</lpage>. doi: <pub-id pub-id-type="doi">10.1021/jf60119a019</pub-id></citation>
</ref>
<ref id="ref21">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Krul</surname> <given-names>E. S.</given-names></name>
</person-group> (<year>2019</year>). <article-title>Calculation of nitrogen-to-protein conversion factors: a review with a focus on soy protein</article-title>. <source>J. Sci. Food Agric.</source> <volume>96</volume>, <fpage>339</fpage>&#x2013;<lpage>364</lpage>. doi: <pub-id pub-id-type="doi">10.1002/aocs.12196</pub-id></citation>
</ref>
<ref id="ref22">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>L&#x00EA;</surname> <given-names>S.</given-names></name> <name><surname>Josse</surname> <given-names>J.</given-names></name> <name><surname>Husson</surname> <given-names>F.</given-names></name></person-group> (<year>2008</year>). <article-title>FactoMineR: an R package for multivariate analysis</article-title>. <source>J. Stat. Softw.</source> <volume>25</volume>, <fpage>1</fpage>&#x2013;<lpage>18</lpage>. doi: <pub-id pub-id-type="doi">10.18637/jss.v025.i01</pub-id></citation>
</ref>
<ref id="ref23">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Marin</surname> <given-names>F. R.</given-names></name> <name><surname>Zanon</surname> <given-names>A. J.</given-names></name> <name><surname>Monzon</surname> <given-names>J. P.</given-names></name> <name><surname>Andrade</surname> <given-names>J. F.</given-names></name> <name><surname>Silva</surname> <given-names>E. H.</given-names></name> <name><surname>Richter</surname> <given-names>G. L.</given-names></name> <etal/></person-group>. (<year>2022</year>). <article-title>Protecting the Amazon forest and reducing global warming via agricultural intensification</article-title>. <source>Nat. Sust.</source> <volume>5</volume>, <fpage>1018</fpage>&#x2013;<lpage>1026</lpage>. doi: <pub-id pub-id-type="doi">10.1038/s41893-022-00968-8</pub-id></citation>
</ref>
<ref id="ref24">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Medic</surname> <given-names>J.</given-names></name> <name><surname>Atkinson</surname> <given-names>C.</given-names></name> <name><surname>Hurburgh</surname> <given-names>C. R.</given-names></name></person-group> (<year>2014</year>). <article-title>Current knowledge in soybean composition</article-title>. <source>J. Sci. Food Agric.</source> <volume>91</volume>, <fpage>363</fpage>&#x2013;<lpage>384</lpage>. doi: <pub-id pub-id-type="doi">10.1007/s11746-013-2407-9</pub-id></citation>
</ref>
<ref id="ref25">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Messina</surname> <given-names>M.</given-names></name>
</person-group> (<year>2022</year>). <article-title>Perspective: soybeans can help address the caloric and protein needs of a growing global population</article-title>. <source>Front. Nutr.</source> <volume>9</volume>:<fpage>909464</fpage>. doi: <pub-id pub-id-type="doi">10.3389/fnut.2022.909464</pub-id>, PMID: <pub-id pub-id-type="pmid">35600829</pub-id></citation>
</ref>
<ref id="ref26">
<citation citation-type="book"><person-group person-group-type="author">
<collab id="coll2">OECD and FAO</collab>
</person-group>. (<year>2022</year>). <source>OECD-FAO Agricultural Outlook 2022-2031</source>. <publisher-name>OECD</publisher-name>, <publisher-loc>Paris, France.</publisher-loc></citation>
</ref>
<ref id="ref27">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Osborne</surname> <given-names>T. B.</given-names></name> <name><surname>Mendel</surname> <given-names>L. B.</given-names></name> <name><surname>Ferry</surname> <given-names>E. L.</given-names></name> <name><surname>Wakeman</surname> <given-names>A. J.</given-names></name></person-group> (<year>1914</year>). <article-title>Amino-acids in nutrition and growth</article-title>. <source>J. Biol. Chem.</source> <volume>17</volume>, <fpage>325</fpage>&#x2013;<lpage>349</lpage>. doi: <pub-id pub-id-type="doi">10.1016/S0021-9258(18)88398-3</pub-id></citation>
</ref>
<ref id="ref28">
<citation citation-type="other"><person-group person-group-type="author"><name><surname>Park</surname> <given-names>H.</given-names></name> <name><surname>Hurburgh</surname> <given-names>C.</given-names></name></person-group> (<year>2002</year>). Improving the US position in world soybean meal trade. Center for Agricultural and Rural Development (CARD) at Iowa State University. Available at: <ext-link xlink:href="http://lib.dr.iastate.edu/cgi/viewcontent.cgi?article=1006&#x0026;context=matric_workingpapers" ext-link-type="uri">http://lib.dr.iastate.edu/cgi/viewcontent.cgi?article=1006&#x0026;context=matric_workingpapers</ext-link></citation>
</ref>
<ref id="ref29">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Pfarr</surname> <given-names>M. D.</given-names></name> <name><surname>Kazula</surname> <given-names>M. J.</given-names></name> <name><surname>Miller-Garvin</surname> <given-names>J. E.</given-names></name> <name><surname>Naeve</surname> <given-names>S. L.</given-names></name></person-group> (<year>2018</year>). <article-title>Amino acid balance is affected by protein concentration in soybean</article-title>. <source>Crop Sci.</source> <volume>58</volume>, <fpage>2050</fpage>&#x2013;<lpage>2062</lpage>. doi: <pub-id pub-id-type="doi">10.2135/cropsci2017.11.0703</pub-id></citation>
</ref>
<ref id="ref30">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Piper</surname> <given-names>E. L.</given-names></name> <name><surname>Boote</surname> <given-names>K. I.</given-names></name></person-group> (<year>1999</year>). <article-title>Temperature and cultivar effects on soybean seed oil and protein concentrations</article-title>. <source>J. Sci. Food Agric.</source> <volume>76</volume>, <fpage>1233</fpage>&#x2013;<lpage>1241</lpage>. doi: <pub-id pub-id-type="doi">10.1007/s11746-999-0099-y</pub-id></citation>
</ref>
<ref id="ref31">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Pomar</surname> <given-names>C.</given-names></name> <name><surname>Andretta</surname> <given-names>I.</given-names></name> <name><surname>Remus</surname> <given-names>A.</given-names></name></person-group> (<year>2021</year>). <article-title>Feeding strategies to reduce nutrient losses and improve the sustainability of growing pigs</article-title>. <source>Front. Vet. Sci.</source> <volume>8</volume>:<fpage>742220</fpage>. doi: <pub-id pub-id-type="doi">10.3389/fvets.2021.742220</pub-id>, PMID: <pub-id pub-id-type="pmid">34778430</pub-id></citation>
</ref>
<ref id="ref32">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Pope</surname> <given-names>M.</given-names></name> <name><surname>Borg</surname> <given-names>B.</given-names></name> <name><surname>Boyd</surname> <given-names>R. D.</given-names></name> <name><surname>Holzgraefe</surname> <given-names>D.</given-names></name> <name><surname>Rush</surname> <given-names>C.</given-names></name> <name><surname>Sifri</surname> <given-names>M.</given-names></name></person-group> (<year>2023</year>). <article-title>Quantifying the value of soybean meal in poultry and swine diets</article-title>. <source>J. Appl. Poult. Res.</source> <volume>32</volume>:<fpage>100337</fpage>. doi: <pub-id pub-id-type="doi">10.1016/j.japr.2023.100337</pub-id></citation>
</ref>
<ref id="ref33">
<citation citation-type="other"><person-group person-group-type="author">
<collab id="coll3">R Core Team</collab>
</person-group> (<year>2021</year>). R: A language and environment for statistical computing. Available at: <ext-link xlink:href="http://r.meteo.uni.wroc.pl/web/packages/dplR/vignettes/intro-dplR.pdf" ext-link-type="uri">http://r.meteo.uni.wroc.pl/web/packages/dplR/vignettes/intro-dplR.pdf</ext-link></citation>
</ref>
<ref id="ref34">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Ravindran</surname> <given-names>V.</given-names></name> <name><surname>Abdollahi</surname> <given-names>M.</given-names></name> <name><surname>Bootwalla</surname> <given-names>S.</given-names></name></person-group> (<year>2014</year>). <article-title>Nutrient analysis, apparent metabolisable energy and ileal amino acid digestibility of full fat soybean for broilers</article-title>. <source>Anim. Feed Sci. Technol.</source> <volume>197</volume>, <fpage>233</fpage>&#x2013;<lpage>240</lpage>. doi: <pub-id pub-id-type="doi">10.1016/j.anifeedsci.2014.08.014</pub-id></citation>
</ref>
<ref id="ref35">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Ricard</surname> <given-names>M. F.</given-names></name> <name><surname>Mayer</surname> <given-names>M. A.</given-names></name> <name><surname>Viglizzo</surname> <given-names>E. F.</given-names></name></person-group> (<year>2022</year>). <article-title>The impact of beef and soybean protein demand on carbon emissions in Argentina during the first two decades of the twenty-first century</article-title>. <source>Environ. Sci. Pollut. Res. Int.</source> <volume>29</volume>, <fpage>20939</fpage>&#x2013;<lpage>20946</lpage>. doi: <pub-id pub-id-type="doi">10.1007/s11356-021-16744-8</pub-id>, PMID: <pub-id pub-id-type="pmid">34750757</pub-id></citation>
</ref>
<ref id="ref36">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Rotundo</surname> <given-names>J. L.</given-names></name> <name><surname>Westgate</surname> <given-names>M. E.</given-names></name></person-group> (<year>2009</year>). <article-title>Meta-analysis of environmental effects on soybean seed composition</article-title>. <source>Field Crop Res.</source> <volume>110</volume>, <fpage>147</fpage>&#x2013;<lpage>156</lpage>. doi: <pub-id pub-id-type="doi">10.1016/j.fcr.2008.07.012</pub-id></citation>
</ref>
<ref id="ref37">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Rotundo</surname> <given-names>J. L.</given-names></name> <name><surname>Borr&#x00E1;s</surname> <given-names>L.</given-names></name> <name><surname>De Bruin</surname> <given-names>J.</given-names></name> <name><surname>Pedersen</surname> <given-names>P.</given-names></name></person-group> (<year>2014</year>). <article-title>Soybean nitrogen uptake and utilization in Argentina and United States cultivars</article-title>. <source>Crop Sci.</source> <volume>54</volume>, <fpage>1153</fpage>&#x2013;<lpage>1165</lpage>. doi: <pub-id pub-id-type="doi">10.2135/cropsci2013.09.0618</pub-id></citation>
</ref>
<ref id="ref38">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Rotundo</surname> <given-names>J. L.</given-names></name> <name><surname>Miller-Garvin</surname> <given-names>J. E.</given-names></name> <name><surname>Naeve</surname> <given-names>S. L.</given-names></name></person-group> (<year>2016</year>). <article-title>Regional and temporal variation in soybean seed protein and oil across the United States</article-title>. <source>Crop Sci.</source> <volume>56</volume>, <fpage>797</fpage>&#x2013;<lpage>808</lpage>. doi: <pub-id pub-id-type="doi">10.2135/cropsci2015.06.0394</pub-id></citation>
</ref>
<ref id="ref39">
<citation citation-type="other"><person-group person-group-type="author"><name><surname>Schuetzenmeister</surname> <given-names>A.</given-names></name> <name><surname>Dufey</surname> <given-names>F.</given-names></name> <name><surname>Schuetzenmeister</surname> <given-names>M. A.</given-names></name> <name><surname>Suggests</surname> <given-names>V. F. P.</given-names></name></person-group>, (<year>2022</year>). <source>Package &#x2018;VCA&#x2019;</source>.</citation>
</ref>
<ref id="ref40">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Selle</surname> <given-names>P. H.</given-names></name> <name><surname>Dorigam</surname> <given-names>J. C. D. P.</given-names></name> <name><surname>Lemme</surname> <given-names>A.</given-names></name> <name><surname>Chrystal</surname> <given-names>P. V.</given-names></name> <name><surname>Liu</surname> <given-names>S. Y.</given-names></name></person-group> (<year>2020</year>). <article-title>Synthetic and crystalline amino acids: alternatives to soybean meal in chicken-meat production</article-title>. <source>Animals</source> <volume>10</volume>:<fpage>729</fpage>. doi: <pub-id pub-id-type="doi">10.3390/ani10040729</pub-id>, PMID: <pub-id pub-id-type="pmid">32331461</pub-id></citation>
</ref>
<ref id="ref42">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Song</surname> <given-names>W.</given-names></name> <name><surname>Yang</surname> <given-names>R.</given-names></name> <name><surname>Wu</surname> <given-names>T.</given-names></name> <name><surname>Wu</surname> <given-names>C.</given-names></name> <name><surname>Sun</surname> <given-names>S.</given-names></name> <name><surname>Zhang</surname> <given-names>S.</given-names></name> <etal/></person-group>. (<year>2016</year>). <article-title>Analyzing the effects of climate factors on soybean protein, oil contents, and composition by extensive and high-density sampling in China</article-title>. <source>J. Agric. Food Chem.</source> <volume>64</volume>, <fpage>4121</fpage>&#x2013;<lpage>4130</lpage>. doi: <pub-id pub-id-type="doi">10.1021/acs.jafc.6b00008</pub-id>, PMID: <pub-id pub-id-type="pmid">27022763</pub-id></citation>
</ref>
<ref id="ref43">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Song</surname> <given-names>X. P.</given-names></name> <name><surname>Hansen</surname> <given-names>M. C.</given-names></name> <name><surname>Potapov</surname> <given-names>P.</given-names></name> <name><surname>Adusei</surname> <given-names>B.</given-names></name> <name><surname>Pickering</surname> <given-names>J.</given-names></name> <name><surname>Adami</surname> <given-names>M.</given-names></name> <etal/></person-group>. (<year>2021</year>). <article-title>Massive soybean expansion in South America since 2000 and implications for conservation</article-title>. <source>Nat. Sust.</source> <volume>4</volume>, <fpage>784</fpage>&#x2013;<lpage>792</lpage>. doi: <pub-id pub-id-type="doi">10.1038/s41893-021-00729-z</pub-id>, PMID: <pub-id pub-id-type="pmid">34377840</pub-id></citation>
</ref>
<ref id="ref44">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Tessari</surname> <given-names>P.</given-names></name> <name><surname>Lante</surname> <given-names>A.</given-names></name> <name><surname>Mosca</surname> <given-names>G.</given-names></name></person-group> (<year>2016</year>). <article-title>Essential amino acids: master regulators of nutrition and environmental footprint?</article-title> <source>Sci. Rep.</source> <volume>6</volume>, <fpage>1</fpage>&#x2013;<lpage>13</lpage>. doi: <pub-id pub-id-type="doi">10.1038/srep26074</pub-id>, PMID: <pub-id pub-id-type="pmid">27221394</pub-id></citation>
</ref>
<ref id="ref45">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Thakur</surname> <given-names>M.</given-names></name> <name><surname>Hurburgh</surname> <given-names>C. R.</given-names></name></person-group> (<year>2007</year>). <article-title>Quality of US soybean meal compared to the quality of soybean meal from other origins</article-title>. <source>J. Sci. Food Agric.</source> <volume>84</volume>, <fpage>835</fpage>&#x2013;<lpage>843</lpage>. doi: <pub-id pub-id-type="doi">10.1007/s11746-007-1107-8</pub-id></citation>
</ref>
<ref id="ref46">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Tillman</surname> <given-names>P. B.</given-names></name>
</person-group> (<year>2019</year>). <article-title>Determination of nutrient values for commercial amino acids</article-title>. <source>J. Appl. Poult. Res.</source> <volume>28</volume>, <fpage>526</fpage>&#x2013;<lpage>530</lpage>. doi: <pub-id pub-id-type="doi">10.3382/japr/pfz010</pub-id></citation>
</ref>
<ref id="ref47">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Waguespack</surname> <given-names>S.</given-names></name> <name><surname>Powell</surname> <given-names>T. D.</given-names></name> <name><surname>Bidner</surname> <given-names>R. L.</given-names></name> <name><surname>Payne</surname> <given-names>L. L.</given-names></name></person-group> (<year>2009</year>). <article-title>Southern, effect of incremental levels of L-lysine and determination of the limiting amino acids in low crude protein corn-soybean meal diets for broilers1</article-title>. <source>Poult. Sci.</source> <volume>88</volume>, <fpage>1216</fpage>&#x2013;<lpage>1226</lpage>. doi: <pub-id pub-id-type="doi">10.3382/ps.2008-00452</pub-id></citation>
</ref>
<ref id="ref49">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>William</surname> <given-names>W.</given-names></name> <name><surname>Dahl</surname> <given-names>B.</given-names></name> <name><surname>Hertsgaard</surname> <given-names>D.</given-names></name></person-group> (<year>2020</year>). <article-title>Soybean quality differentials, blending, testing and spatial arbitrage</article-title>. <source>J. Commod. Mark.</source> <volume>18</volume>:<fpage>100095</fpage>. doi: <pub-id pub-id-type="doi">10.1016/j.jcomm.2019.100095</pub-id></citation>
</ref>
<ref id="ref50">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Winter</surname> <given-names>G.</given-names></name> <name><surname>Todd</surname> <given-names>C. D.</given-names></name> <name><surname>Trovato</surname> <given-names>M.</given-names></name> <name><surname>Forlani</surname> <given-names>G.</given-names></name> <name><surname>Funck</surname> <given-names>D.</given-names></name></person-group> (<year>2015</year>). <article-title>Physiological implications of arginine metabolism in plants</article-title>. <source>Front. Plant Sci.</source> <volume>6</volume>:<fpage>534</fpage>. doi: <pub-id pub-id-type="doi">10.3389/fpls.2015.00534</pub-id>, PMID: <pub-id pub-id-type="pmid">26284079</pub-id></citation>
</ref>
<ref id="ref51">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Zalles</surname> <given-names>V.</given-names></name> <name><surname>Hansen</surname> <given-names>M. C.</given-names></name> <name><surname>Potapov</surname> <given-names>P. V.</given-names></name> <name><surname>Stehman</surname> <given-names>S. V.</given-names></name> <name><surname>Tyukavina</surname> <given-names>A.</given-names></name> <name><surname>Pickens</surname> <given-names>A.</given-names></name> <etal/></person-group>. (<year>2019</year>). <article-title>Near doubling of Brazil&#x2019;s intensive row crop area since 2000</article-title>. <source>Proc. Natl. Acad. Sci.</source> <volume>116</volume>, <fpage>428</fpage>&#x2013;<lpage>435</lpage>. doi: <pub-id pub-id-type="doi">10.1073/pnas.1810301115</pub-id>, PMID: <pub-id pub-id-type="pmid">30559198</pub-id></citation>
</ref>
</ref-list>
</back>
</article>