Soybean Seed Vigor: Uniformity and Growth as Key Factors to Improve Yield

Emergence uniformity and the time required for emergence are essential factors for obtaining highly productive potential in plants. Factors such as sowing depth and soil moisture affect uniformity and emergence, but little is known about the impacts of seed vigor. Thus, we determined the impacts of seed vigor on uniformity and growth as well as development and yield in soybean (Glycine max [L.] Merrill) plants. The treatments consisted of four vigor levels (89%, 57%, 47%, and 43%) obtained by accelerated aging, in a randomized block design, with five replicates at three sites. Seeds with the highest vigor level showed higher uniformity and faster emergence. Dominated plants through their phenotypic plasticity modify stem diameter and internode length. The plants that emerged earlier had a larger leaf area in the three stages (V1, V4, and R2), allowing them to accumulate more photoassimilates in the initial stages. Consequently, these plants exhibited an increase in yield components, especially in the number of fertile nodes and the number of pods per plant, making them more productive.


Introduction
The importance of soybean (Glycine max [L.] Merrill) is undoubted given its role as a major source of protein in food-producing animals and its increasing significance in human consumption [1]. Recently, these considerations, along with changes to human diets (especially in developing countries), have increased the demand for this product. In order to increase the soybean supply to meet this demand, a larger area of land must be cultivated, leading to deforestation in the tropical forests of South America [2].
Such an environmental liability makes it necessary for researchers to develop strategies to fill the narrowing yield gap; through an improvement in yield, the pressure to open new fields for cultivation would be reduced. In this sense, the use of better-quality seeds would be an efficient way to address the challenge of increased food demand [3,4].
Like all products, seeds are also subject to quality standards for their commercialization, with purity (physical quality) and germination (biological quality) being the main standards [4]. However, germination data cannot be extrapolated to the field; this is mainly because germination tests assume ideal conditions, which is contrary to the sub-optimal conditions actually present in the field [5]. Therefore, biological quality should also be measured by another concept, namely vigor, which can be more assertive in the measurement of germination rate, uniformity, and growth [6].
Among the various vigor tests, the mean germination test evaluates the speed of emergence [5]. There are differences among seed lots with different vigor levels in terms of the time required for 47% and 43% of vigor.
Additionally, the date of emergence of the plants in the field was used to establish another factor, ranging from five to 14 days. For these data, we used a completely randomized design, with 200 repetitions in each of the places. We took this measure because we assumed that the day of plant emergence would be a determinant for the results of the variables related to plant development and productivity, supporting the hypothesis that the vigor level decreases the emergence time and that the shorter emergence time is reflected in the most productive plants.

Procedure
The experiment was performed in a no-tillage system with oat (Avena strigosa Schreb) as the previous crop at the three sites; all sites are located in a humid, subtropical climate region with humic dystrophic Red Latosol; soil chemistry properties (Supplementary Table S1) and water soil content (Supplementary Figure S1) had distinct properties at each site. Sowing was accompanied by fertilization with 6 kg·ha −1 N, 69 kg·ha −1 P2O5, and 69 kg·ha −1 K2O. Seeds were inoculated with Bradyrhizobium japonicum (250 mL of inoculant with 7 × 10 9 CFU/mL for 100 kg of seeds) and treated with insecticide (50 g of Fipronil for 100 kg of seeds) and fungicides (5 g of Piraclostrobin and 45 g of thiophanate-methyl for 100 kg of seeds), according to recommendations for soybean crops. Site I was sown on November 20, and sites II and III were sown on November 21 ( Figure  1). Phytosanitary management was made to preventively control diseases, pests, and weeds equally at all sites.  Figure S2). For the plants within these lines, we daily marked the emergence in days after sowing (EDAS) for 15 days. The emergence values of each plot were the means of the emergences of each of the five lines. The plant stand was measured in each plot by counting 6.75 m of the sowing row at 20 days after sowing.
The uniformity index was established by the sum of the percentage of plants emerged on the emergence peak day, plus the day before and after the emergence peak. The uniformity index (Equation (1)) was calculated in each plot by the average uniformity index of all the plants that emerged in five lines of 1.35 m: where NPEP is the number of plants, which emerge in the day with higher plant emergence; NPEBP is the number of plants, which emerge in the day before NPEP; and NPEAP is the number of plants, which emerge in the day after NPEP.

Plant Growth and Development
For the destructive evaluations, 10 plants in one of the marked emergence (EDAS) lines (see Section 2.3.1 and Supplementary Figure S2) were collected from each of the 60 plots. The area of the first trifoliate leaf was measured at stage V 2 (stage 12 102 in the Biologische Bundesanstalt, Bundessortenamt, und Chemische Industrie (BBCH) scale) [14], and the leaf area was determined at stages V 4 (stage 14 104 in BBCH scale) and R 2 (stage 65 605 in BBCH scale). All analyses were performed destructively using a leaf area meter (LI-3100C, LI-COR Biosciences, Lincoln, Nebraska).
Dry root and shoot masses were measured at stage V 4 [14]. For this, the plants were divided into shoots and roots and placed in a drying oven (60 • C) until reaching constant mass. After drying, the dry masses of roots and shoots were determined per plant (g·plant −1 ). Plant height was measured in V 4 and R 8 from the soil surface to the apical meristem of the main stem. The diameter of the stem was measured at ground level with the use of a caliper in R 2 [14].

Nodulation
Nodulation was evaluated destructively in 10 plants from lines with day of emergence marked (EDAS) (see Section 2.3.1 and Supplementary Figure S2) per plot sequentially along the sowing row by counting the number of nodules at stage V 4 . The nodules were oven-dried at 60 • C until reaching constant mass [14]; after drying, the nodule dry mass per plant was determined (g·plant −1 ).

Grain Yield Components
The yield components were evaluated for each plant based in the EDAS in the R 8 stage (stage 99 909 in the BBCH scale). The total number of nodes as well as the number of fertile nodes, pods on the main stem, and seeds on the whole plant (main stem plus branches) were determined [13]. Average node length was determined by the division of the plant height by the total number of nodes. Harvesting was made in the three central rows and the five central meters of each plot, using a plot harvester (WINTERSTEIGER Classic, AT). The samples were weighed and corrected for 13% moisture, and the grain yield (kg·ha −1 ) and thousand-seed weight (TSW) were defined.

Statistical Analysis
The data (seedling emergence, plant stand, uniformity index, and yield) were subjected to analysis of variance as a randomized block design with two factors (vigor levels versus sites), comprised of four vigor levels with five repetitions in three sites, in a total of 60 experimental units. The means were compared with a Tukey test at 0.05 error probability.
The other variables were subjected to analysis of variance as a completely randomized design with EDAS as the only factor, with each of the 200 plants analyzed at each site as an experimental unit. The means were compared with a Tukey test at 0.05 error probability.

Seedling Emergence, Plant Stand, and Uniformity
The differences between the average time of emergence of the highest (89%) and the lowest vigor levels (43%) were found on day 1 for site I, and from 0.5 days for sites II and III (Table 1). We observed no significant differences in plant stands among the sites and treatments, with the treatments showing plant stands between 222.768 and 226.060 plants·ha −1 . Means followed by the same uppercase letter within a row and lowercase letter within a column were not significantly different by Tukey's test (p = 0.05). CV = coefficient of variation.
The uniformity index reduced as the vigor reduced, showing that the decrease in vigor enhanced variability among the plants in the field. There were differences among the sites, with site I being distinct from sites II and III. This was due to the distinct flows of emergence (Figure 2), where site I had beyond its peak a small peak in day 8, causing the low uniformity index.

Soybean Growth and Development
The leaf area of the first trifoliate leaf (V 1 stage) showed that the plants that emerged earlier had a larger leaf area ( Table 2). When comparing plants that emerged within six days with those that emerged within eight days, there was a percentage difference of 11% for site I, 37% for site II, and 21% for site III. The differences between the plants of six and eight EDAS at the V 4 stage were 18% for site I and 29% for site III (site II percentage differences for plants with six EDAS at stage V 4 were not determined due to the absence of an evaluation of these plants). For stage R 2 , this difference was 2, 36, and 18% for locations I, II, and III, respectively.

Nodulation
There was no clear relationship between these factors when assessing nodule dry mass ( Table 2). Sites II and III showed a decreasing pattern with increasing EDAS values, although at site I, the response to the EDAS increase was random. A similar effect was observed for nodule number, where only site II showed a response pattern.  The emergence flow of the soybean seeds. Emergence (%) at each day after sowing (5, 6, 7, 8, 9, 10, 11, 12, 13, and 14) of soybean seeds with different vigor levels (89%, 57%, 47%, and 43%) in three sites (I, II, and III). Table 2. First trifoliate leaf area (stage V 1 ), leaf area (V 4 and R 2 stage), dry nodule mass (stage V 4 ), and number of nodules (stage V 4 ) in soybean plants with different emergence days after sowing (EDAS)* at three sites (I, II, and III). EDAS  I  II  III  I  II  III  I  II  III  I  II  III  I  II  III   5

Shoot, Root, and Stem
The dry mass of shoot and root in V 4 (Table 3) showed a reduction when the EDAS was increased. The weighted average of the dry mass of shoot (Table 3) and the percentage of plant emergence ( Figure 2) multiplied by 220,000 (plant population per hectare) resulted in a biomass accumulation of 1633.32 kg·ha −1 in the vigor of 89% and 1447.21 kg·ha −1 in the vigor of 43% at site I. Therefore, due to the difference of emergence caused by the distinct vigor, we saw a reduction of 11.4% of the total dry mass accumulate in the V 4 stage. Upon analysis, site III had the smallest difference of biomass accumulation of 1370.07 kg·ha −1 in the vigor of 89% and 1301.62 kg·ha −1 in the vigor of 43%, meaning a decrease of 5% of the biomass. We did not estimate the biomass of site II due to the absence of an evaluation of plants with six and 12 EDAS in the V 4 stage.
The ratio of dry mass of shoot by total dry mass remained stable despite the differences of EDAS, with about 80% for shoot and 20% for roots ( Table 3). The diameter of the stem in the R 2 stage showed an inverse proportional response to EDAS (Table 4). When comparing EDAS, we could see a reduction of 5, 12, and 11% (in sites I, II, and III, respectively) in the diameter of plants emerging within six days from plants emerging within eight days.
The shoot length was not very responsive for EDAS in stage V 4 (Table 3), as only the plants with a much higher EDAS (11, 12, and 14) showed a minor decrease in shoot length when compared with earlier emerging plants. This result remains when analyzing shoot length in R 2 (Table 4), where, with a few exceptions, even late emerging plants showed the same shoot length as that of earlier emerging plants. However, the number of nodes of plants decreased as the EDAS increased, resulting in an increase of the average length of internodes.

Yield Components
The yield components were affected by the EDAS in terms of the number of fertile nodes (Table 4), number of pods, number of grains, and grain mass (Table 5). However, the number of grains per plant and grain mass were only a result of the distinct number of pods, once the grains for pods was around 2.40 to 2.55 for all of the EDAS; in addition, the thousand-seed weight did not show a cohesive response to the EDAS, with only some late emerging plants showing an increase in this variable.
When analyzing grain mass, EDAS can be divided into three groups: rapid, with EDAS values between 5 and 7; medium, with values of 8 and 9; or slow, with EDAS values of 10 and more. The rapid group had a grain mass (Table 5)  The final yield ( Figure 3) showed that growth and uniformity variations in seedlings with different vigor levels need to be considered. In the present study, the differences between the vigors of 43% and 89% were 599 (+15.2%), 410 (+10.3%), and 476 (+10.4%) kg·ha −1 for sites I, II, and III, respectively. Table 3. Shoot dry mass, root dry mass, dry mass of shoot ratio by total dry mass (SDM/TDM), shoot length, and root length in stage V 4 in soybean plants from distinct emergence days after sowing (EDAS)* at the three sites (I, II, and III).

Shoot Dry Mass (g)
Root Dry Mass (g) SDM/TDM Shoot Length (cm) Root Length (cm)   EDAS  I  II  III  I  II  III  I  II  III  I  II  III  I  II I  II  III  I  II  III  I  II  III  I  II  III  I  II  III   5 I  II  III  I  II  III  I  II  III  I  II  III  I  II  III   5  The final yield (Figure 3) showed that growth and uniformity variations in seedlings with different vigor levels need to be considered. In the present study, the differences between the vigors of 43% and 89% were 599 (+15.2%), 410 (+10.3%), and 476 (+10.4%) kg·ha −1 for sites I, II, and III, respectively. Our results showed that soybean seed vigor interfered with the time required for seedling emergence. In addition, we observed that the emergence date was a determining factor in the development of soybean plants throughout the cycle. Thus, crops from seeds with low vigor levels are less productive because they have higher percentages of late emergence plants. These plants have lower yield per plant because they presented a smaller number of pods due to the lower accumulation of photoassimilates and the smaller leaf area.
The same situation observed in the laboratory was found in the field when the emergence speed was analyzed, as seeds with a higher vigor level produced a quick emergence [7]. Nevertheless, the differences among the values in the sites showed that the magnitude of the dates were impacted by the environment and that better condition allowed for more distinction in emergence speed: once at site I, where seeds started to emerge at 5 EDAS, we had a higher difference among the speeds of plant emergence.
The greater uniformity established by the uniformity index of the lots with higher vigor levels was already expected; although the environmental conditions interfered with the magnitude of the uniformity data, the lots with greater vigor presented a higher uniformity over a wide range of environmental conditions (Figure 4) [5,15,16]. Again, the best conditions that allowed for quick emergence provided more distinction of the uniformity caused by seed quality, as was found at site I.

Discussion
Our results showed that soybean seed vigor interfered with the time required for seedling emergence. In addition, we observed that the emergence date was a determining factor in the development of soybean plants throughout the cycle. Thus, crops from seeds with low vigor levels are less productive because they have higher percentages of late emergence plants. These plants have lower yield per plant because they presented a smaller number of pods due to the lower accumulation of photoassimilates and the smaller leaf area.
The same situation observed in the laboratory was found in the field when the emergence speed was analyzed, as seeds with a higher vigor level produced a quick emergence [7]. Nevertheless, the differences among the values in the sites showed that the magnitude of the dates were impacted by the environment and that better condition allowed for more distinction in emergence speed: once at site I, where seeds started to emerge at 5 EDAS, we had a higher difference among the speeds of plant emergence. The greater uniformity established by the uniformity index of the lots with higher vigor levels was already expected; although the environmental conditions interfered with the magnitude of the uniformity data, the lots with greater vigor presented a higher uniformity over a wide range of environmental conditions (Figure 4) [5,15,16]. Again, the best conditions that allowed for quick emergence provided more distinction of the uniformity caused by seed quality, as was found at site I.
Thinking in terms of uniformity is necessary because it is a quality parameter of a crop [8] as the crop is a community of plants that influence each other via intraspecific competition [17]. The uniformity of plant emergence reduces this competition, enabling a better use of available resources.
Nonuniformity in seedling emergence generates dominant and dominated plants in the crops, with substantial differences in leaf area and number of nodes [18]. These differences could be mediated by a decrease in the red-far red (R:FR) ratio, impinging on the dominated plants [19]. However, it is important to note that even under non-restrictive conditions such a change in light quality, reduced biomass, increased internodes, and reduced branching in soybean can be observed, albeit without any effects on the harvest index [19]. Thinking in terms of uniformity is necessary because it is a quality parameter of a crop [8] as the crop is a community of plants that influence each other via intraspecific competition [17]. The uniformity of plant emergence reduces this competition, enabling a better use of available resources.
Nonuniformity in seedling emergence generates dominant and dominated plants in the crops, with substantial differences in leaf area and number of nodes [18]. These differences could be mediated by a decrease in the red-far red (R:FR) ratio, impinging on the dominated plants [19]. However, it is important to note that even under non-restrictive conditions such a change in light quality, reduced biomass, increased internodes, and reduced branching in soybean can be observed, albeit without any effects on the harvest index [19].
Some studies have found that the effects of seed vigor would only be relevant in the early stages of the crop [20,21]. However, in this study, the variable leaf area in the V1, V4, and R2 stages maintained the differences observed due to the different EDAS values, so the effects did not dissipate.
Maintenance of the differences in terms of leaf area can be explained by the greater leaf area since the initial stages (V1) of the crop result in greater ground coverage and resource use [13,22]. Thus, the plant is provided with more carbohydrates for development and a growth cycle is initiated in which the larger leaf area provides more carbohydrates and, subsequently, more energy for a positive cycle of plant development [23].
Due to the greater leaf area and carbon assimilation obtained by the dry shoot mass in the V4 stage in plants with lower EDAS, it was expected that the results for the dry mass of nodules and number of nodules would point in the same direction. Since the increased availability of carbohydrates would provide more sugars for nitrogen-fixing bacteria, this would promote an increase in nodule number and dry mass [24].
However, this does not occur. The fact that this variable presents high unevenness among plants may be the core of this problem. To mitigate such problems, it is recommended that evaluations of these variables be made at later stages of the culture, when the nodules are larger and the variations tend to decrease. Moreover, a previous study states that carbohydrate partition between shoot and roots is responsible for modulating the development of nodules [25], showing that the mere increase in dry shoot mass or the larger leaf area is not, in isolation, a good parameter for inferences with nodulation.
It was expected that with increasing EDAS values, the root dry mass would decrease more significantly than the shoot dry mass. In soybean plants, the embryonic axis is highly susceptible to deterioration because the radicle of the embryonic axis is in close contact with the end of the funicle and ensures the interaction with the environment, resulting in enhanced lipid peroxidation [26,27].
Thus, plants with high EDAS values must be more deteriorated, and the radicle must be the structure more affected. However, this does not occur, the differences are generally small, and the partition of dry mass between roots and aerial parts tends to remain stable, independent of the EDAS value. The distribution of carbohydrates, despite being affected by a number of factors (nutrients, water potential, luminosity, temperature, and CO2 concentration), tends to maintain a proportion that favors the balance of shoot and root activities [28].
The shoot length results (Table 3) support differences in the initial developmental capacity of plants with different EDAS values. In this sense, the larger leaf area provides more resources for Some studies have found that the effects of seed vigor would only be relevant in the early stages of the crop [20,21]. However, in this study, the variable leaf area in the V 1 , V 4 , and R 2 stages maintained the differences observed due to the different EDAS values, so the effects did not dissipate.
Maintenance of the differences in terms of leaf area can be explained by the greater leaf area since the initial stages (V 1 ) of the crop result in greater ground coverage and resource use [13,22]. Thus, the plant is provided with more carbohydrates for development and a growth cycle is initiated in which the larger leaf area provides more carbohydrates and, subsequently, more energy for a positive cycle of plant development [23].
Due to the greater leaf area and carbon assimilation obtained by the dry shoot mass in the V 4 stage in plants with lower EDAS, it was expected that the results for the dry mass of nodules and number of nodules would point in the same direction. Since the increased availability of carbohydrates would provide more sugars for nitrogen-fixing bacteria, this would promote an increase in nodule number and dry mass [24].
However, this does not occur. The fact that this variable presents high unevenness among plants may be the core of this problem. To mitigate such problems, it is recommended that evaluations of these variables be made at later stages of the culture, when the nodules are larger and the variations tend to decrease. Moreover, a previous study states that carbohydrate partition between shoot and roots is responsible for modulating the development of nodules [25], showing that the mere increase in dry shoot mass or the larger leaf area is not, in isolation, a good parameter for inferences with nodulation.
It was expected that with increasing EDAS values, the root dry mass would decrease more significantly than the shoot dry mass. In soybean plants, the embryonic axis is highly susceptible to deterioration because the radicle of the embryonic axis is in close contact with the end of the funicle and ensures the interaction with the environment, resulting in enhanced lipid peroxidation [26,27].
Thus, plants with high EDAS values must be more deteriorated, and the radicle must be the structure more affected. However, this does not occur, the differences are generally small, and the partition of dry mass between roots and aerial parts tends to remain stable, independent of the EDAS value. The distribution of carbohydrates, despite being affected by a number of factors (nutrients, water potential, luminosity, temperature, and CO 2 concentration), tends to maintain a proportion that favors the balance of shoot and root activities [28].
The shoot length results (Table 3) support differences in the initial developmental capacity of plants with different EDAS values. In this sense, the larger leaf area provides more resources for faster development [29]. The capacity of the phenotypic plasticity of plants undergoing competition will only be able to match plant height at more advanced stages of development [19].
The reduction of stem diameter with increasing EDAS values is one of the phenotypic responses of soybean plants to shading because, when exposed to low photosynthetically active conditions and low R:FR ratios, plants respond by decreasing stem diameter and increasing stem length [30].
Shoot length at the R 8 stage was more uniform among the different EDAS values than at the V 4 stage (Tables 3 and 4), which is a result of etiolation, in which the changes in light quality stimulate an increase in the average internode length [19]. The increase in the mean length of the internodes led to a reduction in the number of total nodes (Table 4). However, the number of fertile nodes presented a greater relevance, since, in addition to the effects of etiolation, this variable is still impacted by the plant's ability to flower adhesion in the formation of a fertile node [31]. The maintenance of a reproductive node depends on two characteristics: the amount of photoassimilates and the distribution of flowering [32]. The differences in fertile nodes can partly be attributed to the larger leaf area of the plants with lower EDAS values, providing a greater amount of photoassimilates and thus reducing the abscission of flowers and pods.
At the end of the cycle at stage R 8 , the main yield components affected by delays in EDAS were pod and grain number, which decreased (Table 5). Reductions in grain number are generally observed when intraspecific competition occurs at the field level [10]. Two hypotheses can be constructed about this phenomenon. The first is that plants with lower EDAS values have greater flowering and fruiting periods by altering intra-nodal pod interactions, similar to the effects that occur when the photoperiod increases, generating more pods [33]. Another hypothesis is that the greater accumulation of photoassimilates leads to greater fruiting [34]. The number of grains per pod is a characteristic with high uniformity and is not impacted by EDAS (Table 5), since it does not respond to intraspecific competition [35]. The absence of significant variations in the number of grains per pod implies that the number of grains per plant only corroborates the differences caused by the EDAS in the number of pods per plant.
The differences in yield among the distinct groups are a result of the higher trifoliate leaf (Table 2) area of the rapid group, allowing the plants to use more resources (light) for their development [13,22]. This difference generates a higher growth [29] and dry mass accumulation in the V 4 stage (Tables 2  and 3) for the rapid plants [10], making them more competitive for light. These plants therefore do not need to increase their average internode length [30] and have a higher number of total nodes and a higher pod number [34], which explains the differences in grain mass ( Figure 5).
will only be able to match plant height at more advanced stages of development [19].
The reduction of stem diameter with increasing EDAS values is one of the phenotypic responses of soybean plants to shading because, when exposed to low photosynthetically active conditions and low R:FR ratios, plants respond by decreasing stem diameter and increasing stem length [30].
Shoot length at the R8 stage was more uniform among the different EDAS values than at the V4 stage (Tables 3 and 4), which is a result of etiolation, in which the changes in light quality stimulate an increase in the average internode length [19]. The increase in the mean length of the internodes led to a reduction in the number of total nodes (Table 4). However, the number of fertile nodes presented a greater relevance, since, in addition to the effects of etiolation, this variable is still impacted by the plant's ability to flower adhesion in the formation of a fertile node [31]. The maintenance of a reproductive node depends on two characteristics: the amount of photoassimilates and the distribution of flowering [32]. The differences in fertile nodes can partly be attributed to the larger leaf area of the plants with lower EDAS values, providing a greater amount of photoassimilates and thus reducing the abscission of flowers and pods.
At the end of the cycle at stage R8, the main yield components affected by delays in EDAS were pod and grain number, which decreased (Table 5). Reductions in grain number are generally observed when intraspecific competition occurs at the field level [10]. Two hypotheses can be constructed about this phenomenon. The first is that plants with lower EDAS values have greater flowering and fruiting periods by altering intra-nodal pod interactions, similar to the effects that occur when the photoperiod increases, generating more pods [33]. Another hypothesis is that the greater accumulation of photoassimilates leads to greater fruiting [34]. The number of grains per pod is a characteristic with high uniformity and is not impacted by EDAS (Table 5), since it does not respond to intraspecific competition [35]. The absence of significant variations in the number of grains per pod implies that the number of grains per plant only corroborates the differences caused by the EDAS in the number of pods per plant.
The differences in yield among the distinct groups are a result of the higher trifoliate leaf ( Table  2) area of the rapid group, allowing the plants to use more resources (light) for their development [13,22]. This difference generates a higher growth [29] and dry mass accumulation in the V4 stage (Tables 2 and 3) for the rapid plants [10], making them more competitive for light. These plants therefore do not need to increase their average internode length [30] and have a higher number of total nodes and a higher pod number [34], which explains the differences in grain mass ( Figure 5). However, when comparing these results with the results of studies without any adjustment of the plant population, vigor variations from 48% to 90% generated greater productivity differences of 1180; 1376; and 1307 kg·ha −1 for different locations [13]. Thus, it is suggested that adjusting the sowing rate slows down some of the effects of vigor, but does not eliminate them, since in addition to the final plant stand, uniformity and growth will also reverberate in the final yield.
The differences in yield were due to the higher photosynthetic activity in the initial stages, where shading by the crop was still low, and thus the larger leaf area at these times allowed for a greater assimilation of carbon [23]. Consequently, fertile node number and pod number were higher [34], resulting in increased yield.

Conclusions
In this work, we showed that the vigor of soybean seeds impacts the time required for seedling emergence that is, a higher vigor level leads to a faster emergence. In addition, we have shown that the emergence date is a determining factor in the development of soybean plants throughout their growth cycle, as faster emergent plants have higher leaf area indexes from the early stages (V 1 ) to the reproductive stages (R 2 ).
Late emergence plants, through their phenotypic plasticity, adapt to intraspecific competition. By reducing the stem diameter and increasing the average internode length, these late plants are equal in height to early emergence plants. However, despite their phenotypic plasticity, late emergence plants have lower dry grain mass due to the reduction in the number of fertile nodes and the number of pods per plant. Thus, crops from seeds with low vigor levels are less productive because they have higher percentages of late emergence plants, which have a lower yield per plant.
Supplementary Materials: The following are available online at http://www.mdpi.com/2073-4395/10/4/545/s1, Figure S1: Water soil content. Volume moisture (m3/m3) at Field Capacity (Vm FC) and Permanent Wilting Point (Vm PWP) in different soil depth (0.00-0.20 m) at three sites (I, II and III), Figure S2: Plot design. Schematic representation of each plot (7m × 14 rows), with 2 rows and 1 m of each side as borders control to plants with emergence days after sowing (EDAS) evaluated in the five rows of 1.35 m located in the third and twelfth sowing row. Table S1: Chemical analysis of the soils of the three sites (I, II and III).