Cold Stress during Flowering Alters Plant Structure, Yield and Seed Quality of Different Soybean Genotypes

Abstract

The objective of the study was to evaluate the effect of cold stress at flowering stage on plant structure, yield and chemical composition of seeds of 15 soybean cultivars. The study was conducted in 2019–2020, using the complete randomization method, in three replications. Fifteen soybean cultivars belonging to three maturity groups: early (EC), middle-early (MC) and late (LC) cultivars were included. Weekly cold stress (17/13 °C day/night) was applied at plant flowering stage. In the control treatment, plants were kept under natural conditions (24/17 °C day/night). Our research showed that cold stress negatively influenced the elements of plant structure: height, number of nodes, stem dry mass, number and weight of pods, number and weight of seeds per plant, as well as yield of soybean seeds, reducing it on average by 24%, as compared to the control treatment. The highest yield decrease was found in LC cultivars (31.2%), while a smaller and similar one in EC and MC cultivars (by 19.7 and 20.1%, respectively). Significant varietal differences were found for plant biometric traits and seed yield. EC cultivars had the lowest set first pod, as well as the lowest number of nodes, number of pods and seeds, pod and seed weight per plant, and seed yield. MC cultivars had the highest set first pod, and the smallest stem DM and seed yield average by 29.2% higher compared to EC cultivars. LC cultivars had the highest stem DM, number of pods and seeds, and pod and seed weight per plant compared to the other cultivar groups, and yield by 22.8% higher compared to EC cultivars. The experimental factors significantly affected crude protein, crude fat, and crude ash content, while they did not differentiate water-soluble carbohydrate and crude fiber content. Cold stress at the flowering stage caused a significant increase in protein content (by 4.1% on average) and ash content (by 3.8%) and a decrease in fat concentration (by 6.9%) in soybean seeds. Differences in nutrient content among cultivars were a genetic trait not related to cultivar maturity.

1. Introduction

Population growth and strong demand for oil and feed protein contributed to a global increase in soybean acreage in the late 20th and early 21st centuries. Compared to 1970, soybean acreage in 2019 increased from 29.5 to 120.5 million hectares and production increased more than four times from 81 to nearly 334 million tons. There has also been an increase in soybean yields from 1.48 to 2.77 t·ha⁻¹ [1]. Soybean seeds contain around 330 to 450 g kg⁻¹ of protein with a favorable amino acid composition and a biological value similar to beef protein, about 180 to 240 g kg⁻¹ of crude fat, of which more than half are unsaturated fatty acids, and around 55 to 80 g kg⁻¹ of crude fiber [2]. Moreover, they are a source of lecithin, vitamins, mineral salts and biologically active compounds—isoflavones [3]. In addition to the consumption and fodder purposes of soybean, the cultivation of soybean as a legume is additionally economically and ecologically beneficial, resulting from the nitrogen fixation by nodule bacteria Bradyrhizobium japonicum living in symbiosis with the roots of these plants. Under favorable conditions, soybean can bind from 50 to 100 kg N ha⁻¹, which allows to cover the nitrogen demand by 30–60% [4]. Thus, mineral nitrogen fertilization can be significantly reduced, which is of great importance for the protection of the natural environment and brings material benefits.

Climatic conditions are not favorable for soybean cultivation in Poland. The limiting factors are high thermal requirements of this species and sensitivity to day length [5,6,7]. According to Câmara et al. [8], photoperiod and temperature are the most important factors affecting the growth and development of soybean, while at the same time severely limiting the range of its cultivation, with numerous studies indicating that early-maturing cultivars grown in higher latitudes (cooler climate) are less sensitive to day length than late-maturing grown in warmer climates [9,10]. Soybean is temperature-sensitive from emergence to maturity, but the critical period associated with particular sensitivity to low temperatures is the flowering stage, where a biological minimum is considered to be a temperature range of 17–18 °C, while the optimum is 22–25 °C [5]. A drop in average air temperature below 15 °C causes growth retardation and inhibition of new leaf and shoot production, while a drop below 10 °C may even cause plants to fail to enter the flowering stage [8]. Ohnishi et al. [11] showed that low temperatures (15/10 °C day/night) 3–4 days before flowering, had a negative effect on the fertilization process, and consequently on pod setting and seed yield. In a study by Nawracała [12], the number of flowers on the plant and on the main shoot depended on the temperature during the inflorescence formation period. Soybean, on the other hand, responds favorably to higher temperatures throughout the growing season, but only up to 32–38 °C. Very high temperatures during the flowering also result in a reduction in the number and weight of seeds per plant and consequently reduce yield [13,14,15]. According to Mourtzinis et al. [16], adverse weather events from 1994 to 2013 in the United States, associated with high variability in precipitation during the growing season and high temperatures resulted in a 30% reduction in soybean seed production and USD 11 billion in economic losses.

The basis of breeding new cultivars is the consolidation of favorable traits affecting the technological and utilitarian value of plants, such as yield potential, chemical composition, or resistance to stress factors. With new breeding techniques, it is possible to accelerate varietal progress to produce genotypes adapted to specific environmental conditions [17]. The most important adaptive traits in soybean include length of the growing season and flowering time [18]. Under southwestern Polish conditions, as the growing season lengthened, soybean yields increased by, an average of 30 kg ha⁻¹ per day [2]. However, cultivation of late and very late cultivars in the northern part of Poland is risky due to a shorter growing season (by about 30–35 days) [19] and problems with plant maturity. The warming of climate favors the cultivation of thermophilic species in higher latitudes, but breeding work on soybean in Europe and in Poland should lead, first of all, to the production of cultivars with a shorter growing period, increased productivity, and a higher embedded first pod [20,21,22].

Soybean has a wide range of variation in morphological and physiological traits that are modified by habitat conditions. Sensitivity to flowering transition is a characteristic of this species, but there is also considerable variation in temperature response among soybean genotypes [23,24]. Gass et al. [5] showed variation among 10 soybean genotypes in cold tolerance of 3 °C. Other authors have also demonstrated genotypic differences in pod setting ability after low temperature exposure correlated with cold tolerance [25,26]. Cold stress disrupts bioenergetic processes, causes changes in metabolism, and contributes to damage to cellular structures, leading to stunted plant growth and development, resulting in reduced yield and poorer quality [27]. Changes in plant metabolism depend on the temperature range, duration of stress, and rate of temperature decline [28].

The chemical composition of soybean seeds depends on the genetic characteristics of the cultivar and environmental factors such as length of the growing season, soil type, weather conditions, and biotic and abiotic stresses. Protein and fat contents are largely shaped by the pattern of moisture and thermal conditions during vegetation [29,30,31], with a negative correlation shown between protein content and fat content and seed yield [32,33]. Mourtzinis et al. [34] showed an increase in crude fat content in soybean seeds under higher temperature conditions, while Kumar et al. [35] proved a positive relationship between temperature and protein content and a negative relationship between temperature and fat concentration.

The aim of this study was to evaluate the effect of 7-day cold stress at flowering stage on plant structure as well as yield and chemical composition of seeds of 15 soybean cultivars belonging to three maturity groups.

2. Materials and Methods

2.1. Plant Material

The pot experiment was conducted in 2019–2020 in the greenhouse of Institute of Soil Science and Plant Cultivation-State Research Institute in Puławy, Poland [51°24′59″ N, 21°58′09″ E] under partially controlled conditions. The research treatments were 15 cultivars of soybean Glycine max (L.) Merr., listed in the EU Common Catalogue of Varieties of Agricultural Plant Species [36]. The tested cultivars belonged to three maturity groups [37]: very early and early cultivars (EC): Aldana, Annushka, Augusta, Erica, Oressa (length of growing season 124–129 days), medium–early cultivars (MC): Abelina, Maja, Mavka, Merlin, Sculptor (136–139 days) and late and very late cultivars (LC): Aligator, GL Melanie, Lissabon, Madlen, Petrina (140–147 days).

2.2. Experimental Conditions

The soil from the arable layer of the field, was mixed with sand in a ratio of 5:2. The soil was characterized by the following content of available mineral components (g·100 g⁻¹ of soil): P₂O₅ 32.7, K₂O 12.2, Mg 11.0, S 1.34 and microelements (mg·kg⁻¹ of soil): Cu 7.6, Zn 24.1, Mn 213, B 7.2, Mb ˂ 0.2. Organic C content was 0.87%, soil reaction pH_KCl 6.9. Pre-sowing fertilization was performed with macronutrients (g pot⁻¹): nitrogen 0.5 (NH₄NO₃), phosphorus 1.0 (KH₂PO₄), potassium 1.50 (K₂SO₄), magnesium 0.5 (MgSO₄), and microelement nutrient solution. The plants were watered daily automatically (drip system) to ensure optimum moisture conditions. The thermal conditions were natural (mean daily temp. VI–VIII 19.3–22.9 °C). During the day, the pots stood outdoors, and at night and during the rain, they stood in an open, covered greenhouse.

Two days before sowing, soybean seeds were treated with SAROX 75 WS antifungal dressing, and on the day of sowing, they were inoculated with bacterial inoculum (nitragine) containing strains of Bradyrhizobium japonicum. Seeds were sown in early May (5 June 2019, 5 April 2020) into Mitcherlich pots contain 7 kg of soil (20 seeds per pot) and placed in the greenhouse. After 21 days, some seedlings were removed and five plants per pot were left.

At the flowering stage (BBCH 62–67), the pots were placed in the MICRO-CLIMA 1750 phytotron (SNIJDERS LABS) for 7 days, where cold stress was inflicted (17/13 °C day/night), and then moved to the greenhouse where further plant growth took place. The control treatments were under natural conditions at all times (mean temperature 24/17 °C day/night).

2.3. Scope of Tests and Measurements

Plants were harvested at full seed maturity (BBCH 97) with 10–12% moisture content of seeds. At harvest, biometric measurements were taken for 10 plants: height, height to the first pod, number of nodes, stem dry mass, number and weight of pods per plant, and number and weight of seeds per plant. The dry mass of the stem was determined after drying at a temperature of 55 °C. Yield (seed weight per pot) and 1000 grain weight were analyzed from three pots (pot as a replication). In the certified Chemical Laboratory in Puławy, the content of total nitrogen in seeds was determined by the Kjeldahl method [38], crude fat by the Soxhlet method [39], water-soluble carbohydrates by the Bertrand method [40], crude fiber by the enzyme-weighing method [41], and crude ash by conventional method after wet mineralization [41], based on the average samples from objects. The protein content (CP) of the seeds was converted according to the formula: CP = N × 6.25 [42].

2.4. Statistical Analyses

The normal distribution of variables was tested using the Shapiro–Wilk test. The experiment was conducted using the complete randomization method in three replications. Seed yield results were statistically processed by analysis of variance, determining confidence intervals with Tukey’s HSD test at a significance level of p = 0.05, using STATGRAPH Plus for Windows software. The two-way analysis of variance (ANOVA) for plant structure elements, yield, and 1000 seed weight were performed as the mean for 2019 and 2020, where cold stress (two regimes: stress and control) and soybean cultivar (15 cultivars) were the experimental factors. Seed nutrient content was determined only in 2020 and developed by analysis of variance for a two-factor, no-replication experiment, determining half confidence intervals by Tukey’s test at a significance level of p = 0.05, where the interaction was treated as an error, using STATGRAPH Plus for Windows software.

3. Results

3.1. Plant Morphological Features and Elements of the Yield Structure

Cold stress inflicted on plants at flowering as well as cultivar generally significantly affected morphological traits of plants and elements of the yield structure. Cold stress significantly differentiated almost all the studied traits, i.e., plant height (PH), number of nodes per shoot (NN), stem dry mass (SDM), number of pods per plant (NP), pod weight per plant (WP), number of seeds per plant (NS) and seed weight per plant (WS), except for height to first pod (H1P). On the other hand, genetic factor (cultivar) significantly affected all studied traits (PH, NN, SDM, NP, WP, NS, WS, H1P). Statistical analysis of the experimental results showed significant interaction of the experimental factors and its effect only for the number of nodes (

Table 1. Results (p-values) of a two-way ANOVA on the effects of cold stress and cultivar treatments, and their interactions on the plant morphological features and elements of the yield structure.

Treatment	PH	H1P	NN	SDM	NP	WP	NS	WS	SY	TSW
Stress	˂0.0001	0.1946	˂0.0001	˂0.0001	˂0.0001	˂0.0001	˂0.0001	˂0.0001	˂0.0001	0.4856
Cultivar	˂0.0001	˂0.0001	˂0.0001	˂0.0001	˂0.0001	˂0.0001	˂0.0001	˂0.0001	˂0.0001	0.4172
Stress × Cultivar	0.1052	0.2130	0.0430	0.1461	0.4762	0.7115	0.7802	0.5672	0.8385	0.9999

PH plant height, H1P height to the first pod, NN number of nodes, SDM stem dry mass, NP number of pods, WP weight of pods, NS number of seeds per plant, WS weight of seeds per plant, SY seed yield, TSW 1000 seed weight.

).

On average, the plants subjected to the cold stress at the flowering stage had: 6.0% lower height (Figure 1a), 7.3% fewer nodes (Figure 1c), 14.9% lower stem dry mass (Figure 1d), 17.4% lower number of pods per plant (Figure 2a), 22.0% lower pod weight per plant (Figure 2b), 22.1% lower number of seeds per plant (Figure 2c), and 19.1% lower seed weight per plant (Figure 2d) compared to the control. Irrespective of the thermal conditions, significant variation in cultivars was shown for plant morphological features and elements of the yield structure. Early and very early cultivars (EC) had the lowest H1P (Figure 1b), the lowest NN, NP, NS, WP, and WS. Among this group of cultivars, Erica showed the best traits in terms of SDM, WP, NS, and WS, while Oressa showed the best traits in terms of PH, H1P, NN, and NP. On the other hand, the least favorable traits were found in cultivars: Aldana in terms of NP and NS, Augusta—PH, H1P, NN, and SDM and Annushka—WP and WS. Medium-term cultivars (MC) had the highest H1P and the lowest SDM, while the other traits were at an average level, compared to the EC and LC groups. Among the cultivars of MC group, Abelina in terms of PH, H1P, SDM, WP, WS, and Merlin in terms of NN, NP, NS stood out favorably. In contrast, the weakest traits were exhibited by: Sculptor (PH, NN, NP, WP, NS), Maja (SDM, WS) and Mavka (H1P). The late and very late cultivars (LC) showed the highest SDM, NP NS, WP and WS compared to the EC and MC groups. Compared to the EC group, these traits were larger by, respectively, 14.6, 23.4, 29.9, 29.4, and 24.2%. Among LC cultivars, the most favorable traits were found for GL Melanie in terms of SDM, WP, NP, NS, and WS and Petrina in terms of PH, NP. However, Madlen was the weakest cultivar in terms of almost all plant morphological traits and elements of the yield structure (PH, NN, SDM, NP, WP, NS, WS).

3.2. Seed Yield and 1000 Seed Weight (TSW)

Cold stress and cultivar significantly differentiated soybean seed yield, but did not affect thousand seed weight (Figure 3a,b). Cold stress at the flowering stage significantly reduced soybean seed yield (on average by 24.0%) compared to the control. In EC and MC cultivars the average yield decrease was 19.7 and 20.1%, respectively, while in LC cultivars—31.2%.

Irrespective of thermal conditions, significant differences in seed yield were found among soybean cultivars. EC cultivars yielded the lowest, while seed yield of MC and LC cultivars was higher by 29.2 and 22.8%, respectively, compared to EC. Among the tested cultivars, the lowest yield was obtained from the early cultivar Annushka and the late cultivar Madlen, while the highest seed yield was obtained from the late cultivar GL Melanie and medium–early Abelina. Among the cultivars in the EC group, Erica gave significantly higher seed yield compared to the other EC cultivars (with the exception of Oressa), in the MC group Abelina gave significantly higher yield compared to other MC cultivars (with the exception of Merlin). In the LC group Petrina, GL Melanie, Alligator, and Lissabon gave significantly higher yield than Madlen. There was no significant interaction between the studied factors and its effect on yield and TSW of soybean.

3.3. Chemical Composition of Seeds

Regardless of the experimental factors, soybean seeds contained an average of 400 g kg⁻¹ of crude protein, 211 g kg⁻¹ of crude fat, 110 g kg⁻¹ of water-soluble carbohydrates, 62.7 g kg⁻¹ of crude fiber, and 53.3 g kg⁻¹ of crude ash. The experimental factors significantly affected protein, fat, and ash content, while they did not differentiate sugar and fiber content (

Table 2. Results (p-values) of a two-way ANOVA on the effects of cold stress and cultivar treatments on the chemical composition of soybean seeds.

Treatment	CP	CFa	WSC	CF	CA
Stress	0.0006	˂0.0001	0.2586	0.0973	0.0050
Cultivar	0.0002	0.0029	0.0860	0.7155	0.0185

CP—crude protein, CFa—crude fat, WSC—water soluble carbohydrates, CF—crude fiber, Ca—crude ash.

). Cold stress at the flowering stage caused a significant increase in protein content (by 4.1% on average) and ash content (by 3.8%) and a decrease in fat content (by 6.9%) in soybean seeds (Figure 4a–c). There was also a tendency for crude fiber and water-soluble carbohydrate contents to decrease under this factor, but the differences were not statistically proven (Figure 4d,e).

The genetic factor significantly differentiated the protein, fat and ash content. The least protein in seeds was accumulated by cultivars Lissabon, Petrina, and Annushka (respectively, 371, 373, and 376 g kg⁻¹), while considerably more by Madlen, Sculptor, Erica and Maja (respectively, 419, 423, 426, and 437 g kg⁻¹). The remaining cultivars had an average content of this component. The content of crude fat in the seeds of the examined cultivars ranged from 194 to 226 g kg⁻¹. Madlen and Erica accumulated the least fat, while Lissabon and Annushka accumulated significantly more of it. The genetic factor did not significantly affect the content of water-soluble and structural carbohydrates, but it differentiated the content of crude ash. The lowest content of this component in seeds was found in Abelina cultivar (49.0 g kg⁻¹), while considerably higher in Annushka (55 g kg⁻¹).

4. Discussion

Weather conditions are of great importance for the growth and development, and thus the yield of crops. The basis for determining natural and agrotechnical factors affecting the yield of crops are phenological data and dates of particular developmental stages. Linking them to the yield and elements of its structure, and weather conditions occurring at that time allows for finding the weather–yield relationship. This is especially important for crops from other climatic zones, such as soybeans. Soybean is sensitive to thermal conditions throughout its life cycle, but it is particularly sensitive to cold stress during the emergence and flowering [5,11,43]. Analysis of thermal conditions during the critical period, which is the flowering stage, and their effect on final yield and elements of plant structure and yield was useful to determine the effect of subsequent cold stress on plant structure and soybean seed yield quantity and quality.

The results of our own study showed that the weekly cold stress (17/13 °C day/night) negatively affected plant height, number of nodes, stem dry mass, number and weight of pods and seeds per plant, and soybean seed yield, compared to natural conditions (averagely 24/17 °C day/night). In a study by Kurosaki and Yumoto [44], 2-week cold stress (18/13 °C day/night) at the flowering stage, negatively influenced mainly generative organs of soybean (number of pods per plant), while to a small extent vegetative parts such as number of nodes or plant height. According to the authors, plant growth was almost complete at the flowering, hence the differences were not statistically confirmed. Kumagai and Sameshima [10] showed significantly higher pod and seed number and seed yield under higher temperature conditions (daily mean 20–27 °C) compared to cultivation under cooler conditions (19.4–22.6 °C). This was due to prolonged flowering period of soybean, increased photosynthetic rate and increased leaf area. Additionally, Sionit et al. [45] showed that an increase in daily temperature from 18 to 26 °C increased seed number and yield. According to Gass et al. [5], cold stress at the flowering stage causes insufficient pollination of flowers and consequently a low number of pods set. Slight overcooling of plants may result in non-opening of flowers and lack of pollination, resulting in the appearance of small, seedless, or deformed pods. Moderate cold stress leads to flower drooping and reduced pod set or the appearance of barren pods. Severe stress can result in complete flower drop from the plant and a complete lack of yield [11]. In the United States, a temperature increase of 0.8 °C during the postanthesis phase resulted in a 2.4% decrease in soybean seed yield in the southern states of USA (average temperature 26.7 °C), but the same temperature increase in the middle eastern states of USA (average temperature 22.5 °C) resulted in a yield increase of 1.7% [46].

Genetic variation in soybean response to high temperatures (above 30 °C) had been proven [47], but studies of varietal differences in response to low temperature at flowering and the effect of stress on seed yield, are not numerous. Our own study showed that LC cultivars responded more strongly to cold stress as well as exhibited higher yield loss under such conditions than EC and MC cultivars. Kumagai and Sameshima [10] showed that a few degrees increase in temperature under cooler climate conditions (mean daily temperature 19.4–22.6 °C) caused an increase in soybean seed yield in late-maturing cultivars that was not observed in early-maturing cultivars. The authors explain this differential response to temperature during the flowering period by the different growing season length and day length requirements for each cultivar. Under higher temperature conditions, the flowering period of late soybean cultivars was prolonged and, besides, the total number of open flowers increased, which resulted in the setting of more pods and seeds and, consequently, in a higher yield than in early cultivars. This was also confirmed by other authors [34,48]. The differential response to weather conditions of two soybean genotypes is also reported by Kołodziej and Pisulewska [49]. The authors show that the meteorological factors had a stronger influence on the yield of the small-seeded cultivar Nawiko than the large-seeded cultivar Aldana, with seed and fat yields showing the highest susceptibility to weather conditions (coefficient of variation 13.6 and 19.3% in Nawiko, 6.1 and 3.5% in Aldana, respectively). Additionally, in the study of Kurosaki and Yumoto [44], they showed varietal variation in response to low temperature at flowering stage. The more stress-sensitive cultivar (Toyomusume) had significantly lower number of pods per plant (by 64%), number of seeds per pod (by 28%), and 100 seed weight (by 44%) after chilling compared to the control, while the more tolerant cultivar (Hayahikari) had significantly reduced only the number of seeds per pod (by 13%), while the number of pods per plant and 100 seed weight decreased slightly (by 6 and 5%, respectively). Rahman et al. [9] indicate that in the cool climate of New Zealand, soybean development from flowering to maturity was mainly controlled by temperature, as indicated by the high regression coefficients obtained in the regression analysis for two cultivars: Northern Conquest (r² = 0.95) and March (r² = 0.88). Obtaining cultivars with greater tolerance to low temperatures, especially shortly before and during the flowering stage, is one of the most important objectives of soybean breeding in cold climates [11].

The results show that under the influence of cold stress at the flowering stage, protein and ash contents increased significantly, crude fat contents decreased, while the concentration of water-soluble and structural carbohydrates in soybean seeds did not change. In a meta-analysis of environmental studies, Rotundo and Westgate [50] showed a negative relationship between protein and fat concentrations in soybean seeds. This is a difficult task especially for growers who want to increase seed fat content while maintaining high protein levels. Oil is important to the soybean industry because of its high use in edible oil production, but also as a source of a major renewable feedstock in biodiesel production. Under low temperature conditions, various authors have shown an increase in protein content and a decrease in fat content in soybean seeds [51,52,53,54,55,56]. In a study by Kołodziej and Pisulewska [49], an increase by 1 °C of minimal air temperature (13.4 °C) at the soybean flowering stage, resulted in an increase in fat yield by 7.7 kg·ha⁻¹, while the same increase in average daily temperature (19.9 °C) contributed to an increase in fat yield by 23.3 kg·ha⁻¹, with similar values of correlation coefficients (R = 0.916 and 0.922 respectively) and their level of significance (p = 0.01).

Soluble sugars are the main source of energy in the fermentation of sweets in food products such as soy milk [57]. However, certain sugars (raffinose, stachyose) have negative effects on soybean seed quality, contributing to reduced animal performance [58]. Wolf et al. [51] investigated the content of soluble sugars in soybeans after planting in the temperature range 18/13 °C and 33/28 °C day/night. According to the authors, among the analyzed sugars, the sucrose concentration decreased by 56% with the temperature increase by 15 °C, and the stachyose showed a slight decrease, while the remaining sugars remained unchanged.

In this study it was shown that protein, fat and ash contents in soybean seeds significantly depended on the genetic factor (cultivar), which confirms the results of other authors. Piper and Boote [59] showed variation in protein and fat concentration in 20 soybean cultivars, with protein content significantly dependent only on cultivar, while fat content on cultivar and temperature related to latitude. Further, other authors proved significant effect of genotype on protein [54], fat [55,60,61], ash [61], and carbohydrate (raffinose and stachyose) contents [62]. The results presented in this study showed little anecdotal varietal variation in sugar content, which is consistent with the study of Alsajri et al. [55]. On the other hand, Kozak et al. [63] proved that the chemical composition of soybean seeds was more dependent on climatic conditions than on cultivar.

5. Conclusions

Cold stress at the flowering stage negatively affected the elements of plant structure and seed yield of soybean, with the greatest yield decrease shown in late cultivars, while a smaller and similar decrease in early and medium–early cultivars. Significant genetic differentiation was shown in morphological structure and yielding of cultivars. Early cultivars had the smallest parameters of plant structure and yield traits, while late cultivars by the largest ones. Cold stress at the flowering stage caused a significant increase in protein and ash content and a decrease in fat concentration in soybean seeds, while it did not differentiate soluble and structural carbohydrate content. Differences in nutrient content among cultivars were a genetic trait not related to cultivar maturity. The high genetic variability in soybean indicates that it has significant adaptive potential to different climatic conditions, so further research is needed to capture environmental variability in different regions to expand the range of producers in cooler climates.