Environmental and Genetic Regulation of Plant Height in Soybean

doi:10.21203/rs.2.20549/v3

PDF

Research article

Environmental and Genetic Regulation of Plant Height in Soybean

https://doi.org/10.21203/rs.2.20549/v3

This work is licensed under a CC BY 4.0 License

Journal Publication

published 25 Jan, 2021

Read the published version in BMC Plant Biology →

You are reading this older preprint version

Read the latest preprint version →

Background:

Shoot architecture is fundamentally crucial to crop growth and productivity. As a key component of shoot architecture, plant height is known to be controlled by both genetic and environmental factors, though specific details remain scarce.

Results:

In this study, 308 representative soybean lines from a core collection and 168 F₉ soybean progeny were planted at distinct field sites. The results demonstrated the presence of significant genotype × environment interaction (G × E) effects on traits associated with plant height in a natural soybean population. In total, 19 loci containing 51 QTLs (quantitative trait locus) for plant height were identified across four environments, with 23, 13 and 15 being QTLs for SH (shoot height), SNN (stem node number) and AIL (average internode length), respectively. Significant LOD ranging from 2.50 to 16.46 explained 2.80% - 26.10% of phenotypic variation. Intriguingly, only two loci, Loc11 and Loc19-1, containing 20 QTLs, were simultaneously detected across all environments. Results from Pearson correlation analysis and PCA (principal component analysis) revealed that each of the five agro-meteorological factors and four soil properties significantly affected soybean plant height traits, and that the corresponding QTLs had additive effects. Among significant environmental factors, AD (average day-length), AMaT (average maximum temperature), pH, and AN (available nitrogen) had the largest impacts on soybean plant height. Therefore, in spite of uncontrollable agro-meteorological factors, soybean shoot architecture might be remolded through combined efforts to produce superior soybean genetic materials while also optimizing soil properties.

Conclusions:

Overall, the comprehensive set of relationships outlined herein among environment factors, soybean genotypes and QTLs in effects on plant height opens new avenues to explore in work aiming to increase soybean yield through improvements in shoot architecture.

Plant Physiology and Morphology

Plant Molecular Biology and Genetics

Shoot architecture

Plants height

Genotype

QTLs

Agro-meteorological factors

Soil properties

With the world population continually increasing, the demands placed upon agriculture to supply enough food will remain a great challenge for the foreseeable future [1]. Increasing crop yields has often been highlighted as a potential solution for meeting the challenge of feeding our growing population [2, 3]. Field-scale plant traits, such as plant density and lodging resistance, are critical determinants of grain yield for many crops. As such, ideal shoot architecture is considered one of the most important breeding targets for many crops [4, 5]. A key component of ideal shoot architecture is plant height. In rice, wheat and maize, shorter stem lengths contribute to higher yields through improved resistance to lodging [6-10]. For instance, the wide-spread incorporation of semi-dwarf cultivars into wheat and rice breeding programs throughout Asia in the 1960s and 1970s was an important factor in ushering the Green Revolution [11-13]. In soybean, advantageous shoot architectures are considered important components of numerous high yielding semi-dwarf cultivars, such as Hobbit87, Charleston and Apex [8, 9]. At present, the optimal height for current commercial soybean cultivars is typically 70 - 90 cm, with shorter or taller stands leading to yield reductions [14-17]. In short, existing evidence strongly suggests that ideal shoot architectures often depend on suitable plant heights, which, as yet, remains to be fully exploited for developing new high yielding cultivars.

As a characteristically quantitative trait, plant height displays significant variation among genetic backgrounds [18, 19]. Traditional breeding processes are time and labor consuming efforts, while, in modern breeding programs, marker-assisted selection (MAS), which has been successfully applied for many crops, allows for rapid selection of desirable traits [20]. Exploring more genetic resources in MAS efforts might, therefore, accelerate the process of breeding soybeans with suitable plant heights for optimal yields. To date, numerous QTLs associated with plant height have been identified in many crops, and several corresponding genes also have been identified through map-based cloning. Moreover, a portion of these genes have been proven to play critical roles in multiple breeding programs [7, 21]. For example, the Green Revolution gene, sd1 (semi-dwarf), is associated with shortened rice plants, where it also improves lodging resistance [22]. Meanwhile, the well-known wheat dwarf gene, Rht (reduced height), which confers lodging resistance and increased allocation of assimilates to the grain, has been found in 70% of commercial wheat cultivars worldwide [21, 23]. Moreover, overexpressing ZmPIN1a (PIN-FORMED proteins) in maize significantly reduces plant height, internode length and ear height, all of which leads to shoot architectures that thrive in high-density cultivation conditions [24]. Finally, similar functions have been well documented for other QTLs or genes responsible for plant height, including qDH1, qDTH8, D18 (DWARF 18), D61 (DWARF 61), brd1 (brassinosteroid-dependent 1), HTD2 (high tillering and dwarf 2), IPA1 (ideal plant architecture1), MPH1 (MYB-like gene of plant height 1), SLR1 (slender rice1), and Sdd(t) (dominant semi-dwarf) [7, 25-31]. In soybean, 239, 37 and 28 QTLs distributed across most of the 20 soybean chromosomes have been associated with plant height, stem node number and internode length, respectively, according to publicly available data (https://www.soybase.org). Among these QTLs, two loci, Dt1 (indeterminate growth 1) and Dt2 (semideterminate growth 2), have been well documented as associated with soybean shoot architecture and final grain yield. The existing evidence suggests that dt1(determinate growth 1) and Dt2 act synergistically in regulating stem development through termination of apical growth, which leads to decreases in plant height and stem node number [20, 32-35].

Beyond genetic effects, crop shoot architecture is also influenced considerably by environmental conditions [36]. For example, elongation of the main stem may be promoted while the outgrowth of lateral buds is inhibited under low light intensity conditions [5]. Additionally, most plant physiological processes remain active only within the 0 - 40°C temperature range. However, optimal temperatures vary among different physiological processes, with, for example, 18°C being reported as the optimal night temperature for tomato stem elongation [37]. As a facultative short-day plant, soybean is influenced by both day-length and temperature, both of which also play critical roles in the formation of shoot architecture [38-40]. For example, soybean exposed to short photoperiods and high temperatures in low latitude regions typically exhibit early flowering, short periods of vegetative growth, short plant heights, and great reductions in yield [41, 42]. Meanwhile, several sensory loci or genes conferring sensitivity of shoot architecture to photoperiod have been identified and cloned from soybean, most notably E1 - E9 (early flowering and maturity) and J (long-juvenile locus) [43-49]. However, despite these numerous reports on the impacts of photoperiod on soybean shoot architecture, only few experiments have addressed temperature effects. Beyond meteorological factors, soil properties are also known to affect crop shoot architecture. As is widely known, crops require a suitable range of nutrients to meet the demands of growth and development. However, most agriculture soils cannot supply adequate quantities of all nutrients necessary to meet high yielding crop demands, so farmers continue to rely on fertilization. Evidence gathered to date suggests that the main fertilizer elements, N, P and K, might function in shoot architecture development. For instance, in cotton, wheat and rice, supplying N fertilizer leads to increases in plant height through formation of longer internode segments, but not increases in node number [7, 50, 51]. In contrast, K fertilization significantly reduces internode length [52, 53].

Soybean is a major source of oil and protein for food and feed [54, 55], though average yields globally are lower than obtained for other grain crops, such as rice, maize and wheat [2]. According to published data, in 2016, 81.3% of global soybean production was occurring in three countries in North and South America, including the United States, Brazil and Argentina. On the other hand, China is the largest consumer of soybeans, despite the fact that China only accounts for 3.57% of the global soybean production [56]. In China, the major soybean production areas include six disparate regions [57], each with soil properties that are distinct from the other regions. Most soils in the South China region belong to acidic soil types with low pH values and poor nutrient conditions, which is similar to soil conditions in Brazil and Argentina [41, 42, 56, 58]. Soil from the Huanghuaihai region and the lower-middle reaches of the Yangtze River basin tends to have higher pH values and more available nutrients than South China counterparts, which makes them similar to many soil types found across the USA [56, 59]. Despite these similarities between Chinese soils and soils found elsewhere, and in spite of Chinese farmers applying plentiful and, at times, excessive fertilizers in the field, average soybean yields in China (1.8 t ha^-1) are far lower than the average yields obtained in the USA (3.51 t ha^-1), Argentina (3.02 t ha^-1), or Brazil (2.91 t ha^-1) [56]. These situations imply that neither soil properties nor the amount of fertilizers applied are limiting factors for soybean yield in China. This suggests that fertilizer management, which is typically neglected by Chinese farmers and breeders, might be the critical factor for increasing soybean yields to levels in line with the yields reported from leading soybean producing countries. Moreover, excessive fertilization of soybean often leads to significant yield reductions, possibly due to development of poor shoot architectures producing taller and, more massive plants yielding less grain and growing less resistant to lodging [60]. Elucidating the effects of soil properties, particularly nutrient supplies, on shoot architecture development might, therefore, help breeders and farmers to realize higher yields with lower input costs. Previously, soybean genetic resources have been extensively studied for variation in soybean shoot architecture traits (e.g. plant height, node number and internode length). However, information on genetic and environmental impacts on soybean shoot architecture remain largely unknown. In order to address this issue, a recombinant inbred line (RIL) population containing 168 F₉ lines was investigated for the presence of QTLs associated with three shoot architecture traits under four distinct environmental conditions. Further analysis was also conducted to identify correlations between QTLs and important environmental factors. The results presented here may contribute to efforts to breed soybean cultivars optimized for both shoot architectures and adaptation to diverse ranges in soil health properties.

Effects of genotype × environment interactions on plant height traits in soybean

In order to evaluate whether G × E impact plant height in a natural soybean population, 308 representative cultivars from a core soybean germplasm collection [61] were selected and planted at two distinct experimental sites, Boluo (BL, 114.29°E, 23.17°N) and Hainan (HN, 109.48°E, 18.31°N). Three traits related to plant height (SH, shoot height; SNN, stem node number; and AIL, average internode length) were determined from field samples. In these tests, the mean values of SH, SNN and AIL were 81.46%, 34.05% and 36.24% higher, respectively, in BL than that in HN (P value < 0.001) (Fig. 1a, b and c). This demonstrated that the plant height in soybean significantly varied between two distinct environments. Furthermore, genetic analysis suggested that the distributions for the three tested traits measured in two environments were approximately normal according to Kurtosis and Skewness values calculated over three replicates (Table 2). Broad-sense heritability (h²_b) for all the traits under the tested environments varied from 0.74 to 0.92, with generally higher values being observed for SH than for the other two traits (Table 2). Regardless of these relatively small differences among traits, the results herein clearly suggest that variation in SH, SNN and AIL depend mainly on genotypic effects in a single environment. Across locations, however, values of h²_b for SH, SNN and AIL ranged between 0.38 and 0.40, all of which were significantly lower than in individual environments. Taken together, these results strongly suggest that SH, SNN and AIL are all greatly affected by both genotype and environment. In order to further determine G × E, two-way ANOVA was performed. As expected, the results showed that SH, SNN and AIL were significantly all affected (P value < 0.001) by environment, genotype and G × E (Table 3). However, the environment itself consists of many factors, including temperature, day-length, precipitation, soil properties and so on. To sort through these myriad environmental influences, we further evaluated the effects of several primary environmental factors, along with QTLs and QTL × environmental (QTL × E) on the tested traits. Analyzing specific environmental factors in this way might contribute to breeding soybean with shoot architectures optimized for specific sets of environmental conditions.

Phenotypic variation among recombinant inbred lines

Given the prevalence of G × E identified for soybean in the plant height experiments above, two representative soybean accessions were, therefore, selected for developing a RIL population designed to explore QTL × E more fully in soybean. In addition, field characterizations were performed in an expanded set of four geographically distinct growth environments. In these trials, plant height traits of the parental lines, BX10 with the genotype of E1E2E3E4E9dt1dt2tof11Tof12J and BD2 with the genotype of E1E2E3E4E9Dt1dt2Tof11tof12J, significantly varied across the four tested environments, with observed ranges falling between 33.56 and 122.00 for SH, 9.63 and 23.00 for SNN, and 3.43 and 5.27 for AIL (Table 4). Although there were no significant differences observed between parental lines within individual environments, data from the RIL population exhibited maximum and minimum values beyond the parental extremes, and most of the distributions for traits tested across four environments were approximately normal according to Kurtosis and Skewness values calculated over three replicates (Fig. 2). These results suggest that soybean plant height traits are typical quantitative traits and both parents contain one or more genes contributing additively towards the tested traits. When sites were observed separately, the mean values of SH, SNN and AIL significantly varied in the ranges of 33.20 - 112.39, 10.07 - 22.70 and 3.36 - 5.06, respectively (Fig. 1d, e, f and Table 4), implying large impacts of environmental factors on the tested traits. Furthermore, ANOVA results revealed that the variation observed for SH, SNN and AIL among RILs was significantly affected by environment and genotype, individually or in interaction terms (P value < 0.001) (Table 5). This was consistent with the results obtained from using the core collection germplasm cultivars (Table 3). Overall, the results herein demonstrate that the observed RIL population was suitable for further analysis.

Identification of QTLs contributing to plant height traits

A high-density genetic linkage map consisting of 3319 recombinant bin markers had been constructed using the RIL population developed in a previous study [62]. In order to identify significant QTLs, trait mean values were calculated for each RIL line. Subsequent QTL analysis identified a total of 19 significant loci containing 51 QTLs for the three tested traits, with 23, 13 and 15 QTLs being associated with SH, SNN and AIL, respectively. The LOD values of these QTLs ranged from 2.50 to 16.46, and explained 2.80% - 26.10% of phenotypic variation (Additional file 1: Table S1). Within environments, 13, 16, 13 and 9 QTLs were identified at the Zhao County (ZC, 114.48°E, 37.50°N), Hangzhou (HZ, 120.69°E, 30.51°N), Yangzhong (YZ, 118.20°E, 26.17°N) and BL field sites, respectively. However, only two loci, Loc11 and Loc19-1, containing a total of 20 QTLs, were identified in each of the four distinct environments. Interestingly, the additive effect of Loc11 was derived from BX10 and BD2 as determined in the two southern (including YZ and BL) and two northern (ZC and HZ) experimental stations, respectively. In addition, seven loci (QTLs) were significant only for single trait observed within one of the four tested environments. Other loci contributed to variation in two or more traits and/or at least two environments (Additional file 1: Table S1). The variation in significant QTL numbers and the extent of the additive effects of these QTLs suggests that soybean height QTLs might depend in part on specific environmental conditions present within individual sites, resulting in plant height influenced by genotype, environment, and G × E.

QTL contributions to soybean plant height traits under varied environmental conditions

In order to explore the stability of detected QTL contributions to plant height traits, QTL and plant height data from the four tested environments were subjected to principal components analysis (PCA). In this case, the first two principal components accounted for 44.3 and 25.7% of the total trait variation and QTL additive effects, respectively (Fig. 3a). Traits associated with plant height (SH, SNN and AIL) tended to group together, indicating a high correlation among them. In contrast, the total additive QTL effects for plant height traits (i.e. qSHt, qSNNt and qAILt) tended to group separately, to the extent that nearly 90° angles were observed among the directional vectors (Fig. 3a), which is indicative of these effects acting independently. These results suggest that the detected QTLs do not fully explain the extent of variation in plant height traits observed across varied environments, with the fact that most of these 51 QTLs were not significant in one or more tests reinforcing the conclusion that site specific conditions significantly influenced soybean height outcomes. To test this hypothesis, qSHt, qSNNt and qAILt were replaced by total additive QTL effects (qSHs, qSNNs and qAILs) from the corresponding environments in further PCA. Consistent with the previous PCA results, the first two principal components in this test accounted for 59.2% and 16.8% of the total variation, respectively (Fig. 3b). Besides the vector for qSNNs, the other 5 vectors grouped closely together (Fig. 3b), which suggests, consistent with our hypothesis, that the studied traits are highly correlated. On the other hand, the unexpected PCA results for qSNNs, the vector of which deviated considerably from the vector for SNN, strongly implied that environment differences greatly affected the QTLs for SNN. To minimize environment effects, plant height trait data (SH, SNN and AIL) were replaced by corrected data (SHc, SNNc and AILc) and subjected to PCA again. As expected, the first two principal components accounted for most of the variation, in this case, 42.9% and 24.5% of total variation, respectively (Fig. 3c). Additionally, all three vectors of additive effects (qSHs, qSNNs and qAILs) were relatively close to their corresponding traits (SHc, SNNc and AILc). Taken together, all of the results above strongly indicate that both G × E and QTL × E contribute to plant height phenotypes in the tested soybean population.

Genotype × environmental factor interaction effects on plant height traits expressed in RILs

In order to further evaluate the effects of the main environmental factors on soybean plant height traits, correlation analysis and PCA were conducted with data collected for the tested traits, agro-meteorological factors and basic soil chemical properties. Results from PCA clearly showed that the first two principal components accounted for more than 88% of the total variation, and the vectors of AD and AMaT grouped closely with the vectors of SH, AIL and SNN (Fig. 4a). This suggests that both AD and AMaT contribute to enhance SH, SNN and AIL. Although, AMiT, EAT and AT grouped separately from most of the other vectors, their placement below 90°, implies that these three environmental factors might also enhance SH, SNN and AIL (Fig. 4a). This was further supported by the results from Pearson correlation analysis, in which significant correlations were identified among tested traits and agro-meteorological factors and correlation coefficients varied between 0.220 - 0.827 (P value <0.01) (Table 6). Contrasting results were obtained when no vectors for soil factors grouped closely with SH, SNN or AIL (Fig. 4b). Except for the angle between pH and AN, all other angles between the AP and AK vectors and plant height traits were larger than 90°, which suggests that there were positive or negative interaction effects of pH and AN, or AP and AK on plant height traits (Fig. 4b). This was further confirmed in Pearson correlation analysis, in which significant positive correlations were established for pH and AN, and negative correlations for AP and AK with SH, SNN and AIL (Table 6). These results strongly demonstrate that both agro-meteorological and soil properties influence plant height traits, but the agro-meteorological factors largely predominate.

QTL × environmental factor interactions in RILs

In order to further explore the main factors imparting QTL additive effects, Pearson correlation analysis and PCA were also performed for agro-meteorological factors, soil properties and QTLs additive effects. Here, AD and AMaT closely grouped with qSHs and qAILs, while, AMiT, EAT and AT distributed separately (Fig. 5a), which is consistent with the relationships obtained in PCA of environmental factors and plant height traits (Fig. 4a). Interestingly, qSNNs aligned very closely with AMiT, yet were far from AMaT, suggesting that the additive effects of qSNNs increased with either increases in AMiT or reductions in AMaT. The positive relationship between qSNNs and AMiT, as well as, the negative relationship between qSNNs and AMaT were further confirmed by correlation analysis, in which the Pearson correlation coefficient was 0.491 between qSNNs and AMiT, or -0.263 between qSNNs and AMaT (P value < 0.01) (Table 6). Further evaluation of soil properties and plant height traits showed that qSHs were significantly negatively correlated with AP, but positively correlated with pH. Meanwhile, qSNNs exhibited significant negative correlations with AN, and positive correlations with AK, while qAILs had significant positive correlations with two soil factors (pH and AN), but was negatively correlated with AK (Fig. 5b, Table 6). Taken together, these results demonstrate that both agro-meteorological factors and soil properties can significantly affect the additive effects of QTLs in regulating soybean plant height.

In contrast to environmental factors, genetic factors can be easily predicted and manually designed through traditional or modern techniques, such as cross-breeding or genetic modification. Furthermore, once genetic factors have been established, further monitoring of markers is unnecessary. Therefore, mining favorable alleles of QTLs conferring development of ideal plant heights became one of the most economic strategies employed to promote crop yield. Over recent decades, many researchers have attempted to identify stable QTLs regulating soybean plant height under varied environments, with a subset of these efforts seeking to clone the underlying genes [63-71]. To date, more than 304 QTLs have been documented in Soybase (https://www.soybase.org), however, many of the reported effects could not be confirmed in different environments, or their additive effects declined considerably in different conditions [16, 63, 67]. This reinforces the point suggested herein that QTLs effects depend on the specific environment conditions present where the soybeans are being grown. Therefore, it is unsurprising that only 2 loci (20 QTLs) out of the identified 19 loci (51 QTLs) were detected across all of the four distinct environments (Additional file 1: Table S1), and that the 51 detected QTLs could not explain a majority of the phenotype variation observed among RILs grown in the 4 diverse environments (Fig. 3a). Unfortunately, these “environmental QTLs”, which might play critical roles under specific environmental conditions, have been typically neglected in previous studies, possible due to more attention being devoted to detecting QTLs that remain stable under varied environmental conditions. Meanwhile, every advantage have its disadvantage, some QTLs possibly be omitted by using only one method to detect, especially for minor QTLs. In order to detect more genetic loci regulating soybean plant height, two algorithms, MQM and ICIM, were employed in this study. Among the 51 QTLs, more than half QTLs could be simultaneously detected by two methods, whereas 8 and 15 QTLs could be only detected by MQM and ICIM, respectively. We speculated that these inconsistent QTLs mainly derived from different algorithms between MQM and ICIM. However, most of the major QTLs could simultaneously detected by two methods, especially for QTLs which clustered in Loc11 and Loc19-1 and these inconsistent QTLs could explained more genetic variation under specific environmental conditions.

Under a given environment, shoot architecture were considerably regulated by flowering, maturity and growth habit of the soybean plants, and some genes underlined well-known genetic locus, such as E1 - E9, J, Tof11(time of flowering 11), Tof12(time of flowering 12), Dt1and Dt2, were cloned[32-34, 43-49, 72]. Therefore, to further evaluate the affection of these flowering, maturity and growth habit-relate gene on soybean shoot architecture in our RIL population, the genotype of BX10 and BD2 were analysis basing on our recently published re-sequence data [62]. As expected, the sequence of three flowering, maturity and growth habit-related genes, Tof11, Tof12 and Dt1, displayed significant variation between BX10 (tof11Tof12dt1) and BD2 (Tof11tof12Dt1) which possible could cause phenotype variation. Moreover, Tof11 and Dt1 were just located in the two environment stable locus (Loc11 and Loc19-1) which strongly that Tof11 and Dt1 underlying Loc11 and Loc19-1, respectively. Interestingly, the additive effect of Loc11 was derived from BX10 and BD2 as determined in the two southern (including YZ and BL) and two northern (ZC and HZ) experimental stations, respectively. We assumed that this possibly due to the genetic roles of Tof11 relied on photoperiod central gene, E1. Under short-day environments, the expression of E1 was greatly suppressed [47, 72] which significantly impaired the function of Tof11, while the expression of two key FT homologs, FT2a (FLOWERING LOCUS T) and FT5a (FLOWERING LOCUS T), were significantly increased that leading to an earlier time of flowering and maturity and a relatively lower soybean plant height. Contrastingly, under long-day environments, dominant Tof11 gene could significant enhance the expression of E1 whereas significantly impair FT2a and FT5a expression which resulting in a later time of flowering and maturity and a relatively higher soybean plant height. Therefore, in our study, the contrasted functions of Loc11 in low and high latitudes were largely dependent on the expression of its central gene, E1 [72]. In addition, Dt1, which was just located in Loc19-1, is the most well-known gene of growth habit and plant height in soybean [33] and high expression of FT5a could accelerate terminating apical stem growth through inhabiting Dt1 expression in post-flowering stage [73], which strongly implied that Dt1 functions also relied on the expression of E1. Therefore, it is not surprised that the additive effect of Loc19-1 was higher in long-day environments than that in short-day environments (Additional file 1: Table S1). Whatever, exploring and incorporating environmental factors that can regulate effective QTLs into breeding efforts should facilitate the development of new cultivars selected through marker assistant selection (MAS) that are adapted to produce grains in wide ranges of environmental conditions.

In order to facilitate the development of such breeding programs, various ecological environments have been classified and characterized throughout the main soybean producing countries [74-78]. For example, photoperiod and temperature are critical environmental factors that influence soybean shoot architecture development [38-40, 79-82]. In soybean, the effect of photoperiod on a variety of developmental processes has been well described, and more than 10 genetic loci sensitive to photoperiod changes have been cloned [40-42, 44, 83]. The sensitive alleles of these loci may enhance the duration of the soybean juvenile phase under long-day conditions, which leads to taller plants. Moreover, these photoperiod sensitive alleles have also been shown to play critical roles in the process of domestication and improvement, due of their ability to alter shoot architecture and enhance grain yields [44].

In contrast to the number of genes known to be photoperiod sensitive, temperature effects, though well documented, have not yet been adequately explained, and genetic loci sensitive to temperature remain rare. In this study, in order to explain the effects of temperature on soybean plant height, four temperature factors and three plant height traits were observed along with day-length. Interestingly, AMaT appeared to exert influence over the three tested plant height traits, whereas, AMiT, EAT and AT exhibited relatively small impacts (Fig. 4a). In addition, AMaT also affected AIL more than AD (Fig. 4a), which led to considerable impacts of AMaT on the QTLs of AIL (Fig. 5a). On the other hand, while the vector of SNN in PCA grouped with the vectors of SH and AIL (Fig. 4a), the vector of qSNNs was very distinct from those of qAILs and qSHs (Fig. 5a) which seemed that qSNNs did not significantly affect stem node number across environment. In this RIL population, both genotype of E1, the central gene of photoperiod [44, 72], were consistence in two parents. Recently, it was reported that GmFT5a and GmAP1s (APETALA1) could effectively terminate post-flowering stem node number [73, 84]. Then, under long day condition, high expression of E1 could considerably inhibit GmFT5a and GmAP1s expression and significantly increase stem node number [47, 73, 84]. Therefore, the major variation of SNN across environments might be aroused by E1 or E1-depended gene × environment interactions. Whatever, these results strongly indicate that variation of SNN across the tested environments is mainly regulated by G × E, but not QTL × E.

Higher temperatures are known to facilitate soybean node development. For instance, soybean node numbers increased from 18 to 29 and to 40 per plant when the temperature was increased from 30/22 °C to 38/30 °C and to 42/34 °C day/night regimes, respectively [85]. It has also been reported that the number of main stem nodes, plant height and mean internode length of crops increases with increasing temperature [86, 87]. However, no research has yet been conducted to determine the effects of diurnal temperature changes on soybean. For soybean, regions with large diurnal variations in temperature, such as Xinjiang Province in China, typically produce higher soybean yields [88]. In this study, we found that AMiT had a positive impact, and AMaT had a negative impact on enhancing the additive effects of QTLs for SNN. This might help to explain why large fluctuations in diurnal temperature can be beneficial for increasing soybean yield, though further work is needed to reveal the underlying molecular and genetic mechanisms.

Based on the present results, soil pH values appear to exert extensive influence over plant height (Fig. 5b and Table 6), possibly due to the fact that soils with low pH values offer limited bioavailability of N and P. The significant positive correlations were all established for pH and SH, SNN or AIL, as well as pH and qSH or qAIL, except for qSNN according to the PCA and Pearson correlation analysis (Fig. 4b, 5b and Table 6). Indicating that soybean plant height were increased by appropriately increasing the soil pH. So the additive effect of QTL related to soybean plant height traits might were promoted in ZC and HZ experimental stations and suppressed in YZ and BL experimental stations, which caused the additive effect of QTL were derived from BD2 detected in ZC and HZ, and the additive effect of QTL were derived from BX10 detected in YZ and BL. Therefore, on acid soils, fertilizers that can increase soil pH values should be first considered. In contrast, alkaline soils tend to have better nutrient availability conditions, and higher biological nitrogen fixation (BNF) capacities for soybean than their acidic counterparts. Over 70% of the N required for soybean growth can be derived from BNF [89], and excess N fertilizer input not only impairs the BNF capacity for soybean [62], but also leads to taller plants (Fig. 4b), which leads to poor lodging resistance. In addition, long-term fertilization with excessive amounts of N causes soil acidification [90-92], which often leads to deteriorating soil conditions. Therefore, in regions harboring alkaline soils, the amount of N fertilizers should be strictly controlled. Contrasted to K, fertilizers rich in P possible enhance AIL and decline SNN which unfavorable for final yield. However, P, which is critical for flower number, poding and filling grain [93, 94], is easily fixed by soil particles [95] or quickly leached out into water supplies, especially in acid soil of southern field and inadequate P may cause more serious yield lost. Therefore, fertilizers both rich in P and K should be considered for more extensive application.

On the whole, the present study provides comprehensive results that contribute to understanding the relationships among environment, genotype, QTLs and soybean shoot architecture. Most importantly, these results also suggest that shoot architecture can be regulated not only by genetic modulators, but also by management strategies designed to optimize soil properties for soybean production. As such, this research opens new avenues for formulating strategies to breed soybean cultivars with improved shoot architectures geared towards sustainable production of high soybean yields in diverse environments.

Plant materials

A total of 308 representative soybean cultivars selected from an applied core germplasm collection [61] were included along with 168 F₉ RIL progeny in tests for interactions between genotype and environment in effects on plant height traits, including shoot height (SH), stem node number (SNN) and average internode length (AIL). Experiments were conducted in geographically distinct field environments. Two cultivars, BX10 and BD2 with contrasted phenotype in flowering, photoperiod sensitivity, shoot architecture, and adaption ability in acidic soils, were selected to construct the RIL population using the single seed descent (SSD) method [96]. This RIL population was used to construct a genetic linkage map of QTLs for soybean plant height traits, as well as, to explore the genetic mechanisms underlying QTL × E. In addition, basing on analysis of re-sequence data [61], the genotype of some well-known flowering, maturity, and growth habit-related genes in BX10 and BD2 were E1E2E3E4E9dt1dt2tof11Tof12J and E1E2E3E4E9Dt1dt2Tof11tof12J, respectively.

Field trials

The 308 soybean germplasm selections were planted in Boluo (BL, Guangdong province 114.29°E, 23.17°N), in 2018, and at the Hainan (HN, Hainan province 109.48°E, 18.31°N) experimental station in 2019. The 168 RIL progeny were grown at four experimental stations differing in agro-meteorological conditions and basic soil properties (Table 1). Specifically, these sites included the Zhao County (ZC, Hebei province 114.48°E, 37.50°N) experimental farm of the Institute of Genetics and Developmental Biology, Chinese Academy of Sciences, the Hangzhou (HZ, Zhejiang province 120.69°E, 30.51°N) experimental farm of the Institute of Crop and Nuclear Technology Utilization, Zhejiang Academy of Agricultural Sciences, BL and the Yangzhong (YZ, Fujian province 118.20°E, 26.17°N) experimental station of Fujian Agriculture and Forestry University. All of the trials were laid out as randomized complete block designs with three replications. Thirty seeds of each genotype were sown per plot in single 3 m rows spaced 0.5 m apart. None of the experiments were fertilized during soybean growth, and all of them incorporated consistent field management practices.

Plant sampling and genetic analysis

SH and SNN of three representative plants from each line were directly measured in the field at the R6 stage, and AIL was calculated as the ratio of SH to SNN. All the data were used to determine the effect of genotype and environment on the tested traits through Two-Way ANOVA in SPSS 19 [97], and to estimate the broad sense heritability of each trait in each or all environments using the formula h²b = VG/(VG + VE), with VG and VE as the respective variance between and within RILs.

Measurements of agro-meteorological and basic soil chemical properties

The agro-meteorological data from each field site are listed in Table 1. Among measured conditions, maximum temperature (MaT, °C), minimum temperature (MiT, °C) and average day-length (AD, hours·day^-1) were obtained from weather data deposited at http://tianqi.2345.com/ and https://www.51240.com/. Average maximum temperature (AMaT, °C), average minimum temperature (AMiT, °C) and accumulated temperature (AT, °C·d) were calculated as:

Meanwhile, due to effective accumulated temperature (EAT, °C·d), which mean the sum of the difference value between the daily average temperature and biological zero point, could significantly affects plant growth and development [98], the EAT for soybean during seeding to R6 stage were also evaluated. For soybean, the biological zero point is 10°C and EAT for soybean was calculated as follows:

The basic soil chemical properties of the top 0 - 20 cm of soil at each location were determined using 10 randomly collected soil samples from each experimental site. The soil pH, available nitrogen (AN), available phosphorus (AP) and available potassium (AK) as measured according to soil and agricultural chemistry analysis protocols [99] are listed in Table 1.

Genetic linkage map and QTL Mapping

Based on a previously constructed genetic linkage map [62], the mean value of each trait from three plants in each plot was used to identify significant quantitative trait loci (QTL) using QTL IciMapping version 4.1 running the inclusive composite interval mapping (ICIM) method [100], as well as in MapQTL6.0 running interval mapping (IM) and Multiple-QTL model (MQM) algorithms [101]. For QTL IciMapping version 4.1, the mapping method of ICIM-ADD was selected to identify QTLs with the step width, probability in stepwise regression and LOD threshold being set to 1 cM, 0.001 and 2.5, respectively. First, IM analysis was conducted for QTL analysis, and the markers with the highest LOD scores were selected as cofactors to carry out MQM analysis. QTLs with LOD score exceeding 2.5 were considered as high confidence QTLs in MQM mapping. To precisely evaluate extensive affection of the candidate QTLs, separated QTLs which clustered together on linkage group were combined as a genetic locus.

Evaluation of genetic effects on plant height parameters measured at four distinct locations

The sum of additive effects for QTLs identified in each single environment (qSHs, qSNNs and qAILs) or all four environments (qSHt, qSNNt and qAILt) were evaluated by principal component analysis (PCA). The qSHs, qSNNs, qAILs, qSHt, qSNNt and qAILt were calculated as follows:

PCA analysis

QTL × E effects on soybean plant height were assessed by PCA performed with the tested traits, agro-meteorological data, basic soil properties, additive eﬀects of QTLs and corrected values for each trait using R with the packages ‘factoextra ’, ‘factoMineR’ and ‘ggplot2’, along with the function ‘fviz_pca_var’ [102-104] (version 3.6.1, https://www.r-project.org/).

SH: shoot height; SNN: stem node number; AIL: average internode length; SSD: single seed descent; ZC: Zhao County; HZ: Hangzhou; YZ: Yangzhong; BL: Boluo; HN: Hainan; MaT: maximum temperature; MiT: minimum temperature; AD: average day-length; AMaT: Average maximum temperature; AMiT: average minimum temperature; AT: accumulated temperature; EAT: effective accumulated temperature; AN: available nitrogen; AP: available phosphorus; AK: available potassium; QTL: quantitative trait locus; h²b: broad-sense heritability; qSHs: sum of additive eﬀects of QTLs for SH in single environments; qSNNs: sum of additive eﬀects of QTLs for SNN in single environments; qAILs: sum of additive eﬀects of QTLs for AIL in single environments; qSHt: sum of additive eﬀects of QTLs for SH under all environments; qSNNt: sum of additive eﬀects of QTLs for SNN under all environments; qAILt: sum of additive eﬀects of QTLs for AIL under all environments; SHc: corrected values for soybean SH; SNNc: corrected values for soybean SNN; AILc: corrected values for soybean AIL; Env: environment; SD: standard deviation; MIN: minimum values; MAX: maximum values; CV%: coefficient of variation; Kurt: kurtosis; Skew: skewness; PVE(%): percentage of phenotypic variance explained by the QTL; Add: additive eﬀects

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Availability of data and materials

The datasets generated during and/or analysed during the current study are available from the corresponding author on reasonable request ([email protected]).

Competing interests

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.

Funding

This work was supported by National Natural Science Foundation of China (31830083) and Innovative Research Groups of the Natural Science Foundation of Hebei province (C2020301020). The funding bodies had no role in the design of the study, the collection, analysis, and interpretation of data and in writing the manuscript.

Authors’ contributions

HL and YY designed the experiments and critically revised the manuscript. QY and YY analyzed the data and wrote the manuscript. QY, YY, GL, HLv and CW, carried out the experiments. All authors have read and approved the manuscript.

Acknowledgments

We would like to acknowledge staff members of the Institute of Genetics and Developmental Biology, Chinese Academy of Sciences, Institute of Crop and Nuclear Technology Utilization, Zhejiang Academy of Agricultural Sciences and Institute of Cereal and Oil Crops, Hebei Academy of Agricultural and Forestry Sciences. We also would like to thank the graduate students of the Root Biology Center, Fujian Agriculture and Forestry University for soybean harvesting, and Dr. Thomas Walk of Golden Fidelity LLC for critical reviewing.

Khan MS, Zaidi A, Wani PA. Role of phosphate-solubilizing microorganisms in sustainable agriculture - a review. Agron Sustain Dev. 2007;27(1):29-43.
Ray DK, Mueller ND, West PC, Foley JA. Yield trends are insufficient to double global crop production by 2050. PloS One. 2013;8(6):e66428.
Jez JM, Lee SG, Sherp AM. The next green movement: plant biology for the environment and sustainability. Science. 2016;353(6305):1241-4.
Reinhardt D, Kuhlemeier C. Plant architecture. EMBO Rep. 2002;3(9):846-51.
Wang B, Steven MS, Li JY. Genetic regulation of shoot architecture. Annu Rev Plant Biol. 2018;69(1):437-68.
Chairi F, Sanchez-Bragado R, Serret MD, Aparicio N, Nieto-Taladriz MT, Luis Araus J. Agronomic and physiological traits related to the genetic advance of semi-dwarf durum wheat: The case of Spain. Plant Sci. 2019;110210.
Chen X, Xu P, Zhou J, Tao D, Yu D. Mapping and breeding value evaluation of a semi-dominant semi-dwarf gene in upland rice. Plant Diversity. 2018;40:238-44.
Cooper RL, Martin RJ, St. Martin SK, Calip-DuBois A, Fioritto RJ, Schmitthenner AF. Registration of ‘Charleston’ soybean. Crop Sci. 1995;35(2):593.
Cooper RL, Mendiola T, St. Martin SK, Fioritto RJ, Dorrance AE. Registration of ‘Apex’ soybean. Crop Sci. 2003;43(4):1563-4.
Su C. QTL mapping, validation and candidate genes analysis for plant height in maize. Indian J Genet Pl Br. 2018;78(4):443-53.
Hedden P. The genes of the green revolution. Trends Genet. 2003;19(1):5-9.
Khush GS. Green revolution: the way forward. Nat Rev Genet. 2001;2(10):815-22.
Peng J, Richards DE, Hartley NM, Murphy GP, Devos KM, Flintham JE, et al. ‘Green revolution’ genes encode mutant gibberellins response modulators. Nature. 1999;400(6741):256-61.
Chen Y, Nelson RL. Variation in early plant height in wild soybean. Crop Sci. 2006;46(2):865-9.
Huang ZW, Wang W, Xu XJ, Wen ZX, Li HC, Li JY, et al. Relationship of dynamic plant height and its relative growth rate with yield using recombinant inbred lines of soybean. Acta Agronomica Sinica. 2011;37(3):559-62.
Josie J, Alcivar A, Rainho J, Kassem MA. Genomic regions containing QTL for plant height, internodes length, and flower color in soybean [Glycine max (L.) Merr.]. Bios. 2007;78(4):119-26.
Thompson J, Bernard RL, Nelson RL. A third allele at the soybean dtl Crop Sci. 1997;37(3):757-62.
Wang Y, Cheng LR, Sun Y, Zhou Z, Zhu LH, Xu ZJ, et al. Effect of genetic background on QTLs for heading date and plant height and interactions between QTL and environment using reciprocal introgression lines in rice. Acta Agronomica Sinica. 2009;35(8):1386-94.
Zhang J, Song Q, Cregan PB, Nelson RL, Wang X, Wu J, et al. Genome-wide association study for flowering time, maturity dates and plant height in early maturing soybean (Glycine max) germplasm. BMC Genomics. 2015;16:217.
Cao Y, Li SG, He XH, Chang FG, Kong JJ, et al. Mapping QTLs for plant height and flowering time in a Chinese summer planting soybean RIL population. Euphytica. 2017;213:39.
Miedaner T, Herter CP, Ebmeyer E, Kollers S, Korzun V, Buerstmayr H. Use of non-adapted quantitative trait loci for increasing Fusarium head blight resistance for breeding semi-dwarf wheat. Plant Breeding. 2019;138:140-7.
Srivastava D, Shamim M, Mishra A, Yadav P, Kumar D, Pandey P, et al. Introgression of semi-dwarf gene in Kalanamak rice using marker-assisted selection breeding. Curr Sci India. 2019;116(4):597-603.
Grover G, Sharma A, Gill HS, Srivastava P, Bains NS. Rht8 gene as an alternate dwarfing gene in elite Indian spring wheat cultivars. PloS One. 2018;13(6):e0199330.
Li Z, Zhang X, Zhao Y, Li Y, Zhang G, Peng Z, et al. Enhancing auxin accumulation in maize root tips improves root growth and dwarfs plant height. Plant Biotechnol J. 2018;16:86-99.
Wei X, Xu J, Guo H, Jiang L, Chen S, Yu C, et al. DTH8 suppresses flowering in rice, influencing plant height and yield potential simultaneously. Plant Physiol. 2010;153(4):1747-58.
Mori M, Nomura T, Ooka H, Ishizaka M, Yokota T, Sugimoto K, et al. Isolation and characterization of a rice dwarf mutant with a defect in brassinosteroid biosynthesis. Plant Physiol. 2002;130(3):1152-61.
Liu W, Wu C, Fu Y, Hu G, Si H, Zhu L, et al. Identification and characterization of HTD2: a novel gene negatively regulating tiller bud outgrowth in rice. Planta. 2009;230(4):649-58.
Lu Z, Yu H, Xiong G, Wang J, Jiao Y, Liu G, et al. Genome-Wide Binding Analysis of the Transcription Activator IDEAL PLANT ARCHITECTURE1 Reveals a Complex Network Regulating Rice Plant Architecture. Plant Cell. 2013;25(10):3743-3759.
Zhang YX, Yu CS, Lin JZ, Liu J, Liu B, Wang J, et al. OsMPH1 regulates plant height and improves grain yield in rice. PloS One. 2017;12(7):1-17.
Zhang YH, Bian XF, Zhang SB, Ling J, Wang YJ, Wei XY, et al. Identification of a novel gain-of-function mutant allele, slr1-d5, of rice DELLA protein. J Integr Agr. 2015;15(7):1441-1448.
Liu BM, Wu YJ, Fu XD, Qian Q. Characterizations and molecular mapping of a novel dominant semi-dwarf gene Sdd(t) in rice (Oryza sativa). Plant Breeding, 2008;127(2):125-130.
Bernard RL. Two genes affecting stem termination in soybeans. Crop Sci. 1972;12(2):235-9.
Liu B, Watanabe S, Uchiyama T, Kong F, Kanazawa A, Xia Z, et al. The soybean stem growth habit gene Dt1 is an ortholog of Arabidopsis TERMINAL FLOWER1. Plant Physiol. 2010;153(1):198-210.
Ping J, Liu Y, Sun L, Zhao M, Li Y, She M, et al. Dt2 is a gain-of-function MADS-domain factor gene that specifies semideterminacy in soybean. Plant Cell. 2014;26(7):2831-42.
Ting CL. Genetic Studies on the Wild and Cultivated Soybeans. Journal of The American Society of Agronomy. 1946;38(5):381-393.
Osiru MO, Olanya OM, Adipala E, Kapinga R, Lemaga B. Yield stability analysis of Ipomoea batatus cultivars in diverse environments. Aust J Crop Sci. 2009;3(4):213-20.
Went FW. The effect of temperature on plant growth. Annu Rev Plant Physiol Plant Mol Bioi. 1953;4(1):347-62.
Zhang SR, Wang H, Wang Z, Ren Y, Niu L, Liu J, et al. Photoperiodism dynamics during the domestication and improvement of soybean. Sci China Life Sci. 2017;60(12):1416-27.
Alliprandini LF, Abatti C, Bertagnolli PF, Cavassim JE, Gabe HL, Kurek A, et al. Understanding soybean maturity groups in Brazil: environment, cultivar classification, and stability. Crop Sci. 2009;49(3):801-8.
Gupta S, Bhatia VS, Kumawat G, Thakur D, Singh G, Tripathi R, et al. Genetic analyses for deciphering the status and role of photoperiodic and maturity genes in major Indian soybean cultivars. J Genet. 2017;96(1):147-54.
Abrahão GM, Costa MH. Evolution of rain and photoperiod limitations on the soybean growing season in Brazil: The rise (and possible fall) of double-cropping systems. Agr Forest Meteorol. 2018;256-257:32-45.
Cober ER, Curtis DF, Stewart DW, Morrison MJ. Quantifying the effects of photoperiod, temperature and daily irradiance on flowering time of soybean isolines. Plants. 2014;3(4):476-97.
Li XM, Fang C, Xu ML, Zhang FG, Lu SJ, Nan HY, et al. Quantitative trait locus mapping of soybean maturity gene E6. Crop Sci. 2017;57(5):2547-54.
Lu S, Zhao X, Hu Y, Liu S, Nan H, Li X, et al. Natural variation at the soybean J locus improves adaptation to the tropics and enhances yield. Nat Genet. 2017;49(5):773-9.
Watanabe S, Xia Z, Hideshima R, Tsubokura Y, Sato S, Yamanaka N, et al. A map-based cloning strategy employing a residual heterozygous line reveals that the GIGANTEA gene is involved in soybean maturity and flowering. Genetics. 2011;188(2):395-407.
Xu M, Xu Z, Liu B, Kong F, Tsubokura Y, Watanabe S, et al. Genetic variation in four maturity genes affects photoperiod insensitivity and PHYA-regulated post-flowering responses of soybean. BMC Plant Biol. 2013;13:91.
Xu M, Yamagishi N, Zhao C, Takeshima R, Kasai M, Watanabe S, et al. The soybean-specific maturity gene E1 family of floral repressors controls night-break responses through down-regulation of FLOWERING LOCUS T Plant Physiol Bioch. 2015;168(4):1735-46.
Zhao C, Takeshima R, Zhu J, Xu M, Sato M, Watanabe S, et al. A recessive allele for delayed flowering at the soybean maturity locus E9 is a leaky allele of FT2a, a FLOWERING LOCUS T BMC Plant Biol. 2016;16:20.
Zhao L, Li M, Xu C, Yang X, Li D, Zhao X, et al. Natural variation in GmGBP1 promoter affects photoperiod control of flowering time and maturity in soybean. Plant J. 2018;96(1):147-62.
Malik MNA, Edwards DG, Evenson JP. Effects of flower bud removal and nitrogen supply on growth and development of cotton (Gossypium hirsutum). Aust J Plant Physiol. 1981;8(3):285-91.
Xiong WB, Xu FY, Wang XY. Effect of different nitrogen application rate on rice stem characteristics. Agricultural Biotechnology. 2018;7(5):204-7.
Moura WM, Soares YJ, Amaral Junior AT, Gravina GA, Barili LD, Vieira HD. Biometric analysis of arabica coffee grown in low potassium nutrient solution under greenhouse conditions. Genet Mol Res. 2016;15(3):gmr.15038753.
Sun JW, Li N, Wang CY, Zhao JH, Zhang SW, Jiang MJ, et al. Effects of transplanting methods and potassium rates on lodging resistance of hybrid rice. Journal of Nuclear Agricultural Sciences. 2017;31(12):2408-17.
Lee S, Jun TH, Michel AP, Rouf Mian MA. SNP markers linked to QTL conditioning plant height, lodging, and maturity in soybean. Euphytica. 2014;203(3):521-32.
Zhang X, Wang W, Guo N, Zhang Y, Bu Y, Zhao J, et al. Combining QTL-seq and linkage mapping to fine map a wild soybean allele characteristic of greater plant height. BMC genomics. 2018;19(1):1-12.
Terzić D, Popović V, Tatić M, Vasileva V, Đekić V, Ugrenović, et al. Soybean area, yield and production in world. Eco-Conference. 2018;10:135-44.
Wang YS, Gai JY. Study on the ecological regions of soybean in China Ⅱ. Ecological environment and representative varieties. Chinese Journal of Applied Ecology. 2002;13(1):71-5.
Maki T, Nomachi M, Yoshida S, Ezawa T. Plant symbiotic microorganisms in acid sulfate soil: significance in the growth of pioneer plants. Plant Soil. 2008;310(1-2):55-65.
Xu X, He P, Pampolino MF, Li Y, Liu S, Xie J, et al. Narrowing yield gaps and increasing nutrient use efficiencies using the nutrient expert system for maize in northeast china. Field Crop Res. 2016;194:75-82.
Zhang XC，Chen H，Huang SL，Yin XW，Du CZ，Zhang JJ. Optimal combination of nitrogen fertilizer and spring soybean varieties in Chongqing. Soybean Science. 2012;31(2):255-8.
Zhao J, Fu JB, Liao H, He Y, Nian H, Hu YM, et al. Characterization of root architecture in an applied core collection for phosphorus efficiency of soybean germplasm. Chinese Sci Bull. 2004;49(15):1611-20.
Yang Q, Yang YQ, Xu RN, Lv HY, Liao H. Genetic analysis and mapping of QTLs for soybean biological nitrogen fixation traits under varied field conditions. Front Plant Sci. 2019;10:75.
Eskandari M, Cober ER, Rajcan I. Genetic control of soybean seed oil: II. QTL and genes that increase oil concentration without decreasing protein or with increased seed yield. Theor Appl Genet. 2013;126(6):1677-87.
Gai JY, Wang YJ, Wu XL, Chen SY. A comparative study on segregation analysis and QTL mapping of quantitative traits in plants-with a case in soybean. Front Agr China. 2007;1(1):1-7.
Lark KG, Chase K, Adler F, Mansur LM, Orf JH. Interactions between quantitative trait loci in soybean in which trait variation at one locus is conditional upon a specific allele at another. P Natl Acad Sci USA. 1995;92(10):4656-60.
Lee SH, Bailey MA, Mian MAR, Shipe ER, Ashley DA, Parrott WA, et al. Identification of quantitative trait loci for plant height, lodging, and maturity in a soybean population segregating for growth habit. Theor Appl Genet. 1996;92(5):516-23.
Orf JH, Chase K, Jarvik T, Mansur LM, Cregan PB, Adler FR, et al. Genetics of soybean agronomic traits: I. comparison of three related recombinant inbred populations. Crop Sci. 1999;39(6):1642-51.
Pathan SM, Vuong T, Clark K, Lee JD, Shannon JG. Genetic mapping and confirmation of quantitative trait loci for seed protein and oil contents and seed weight in soybean. Crop Sci. 2013;53(3):765-74.
Sun DS, Li WB, Zhang ZC, Chen QS, Ning HL, Qiu LJ, et al. Quantitative trait loci analysis for the developmental behavior of soybean (Glycine max Merr.). Theor Appl Genet. 2006;112(4):665-73.
Wang D, Graef GL, Procopiuk AM, Diers BW. Identification of putative QTL that underlie yield in interspecific soybean backcross populations. Theor Appl Genet. 2004;108(3):458-67.
Yao D, Liu ZZ, Zhang J, Liu SY, Qu J, Guan SY, et al. Analysis of quantitative trait loci for main plant traits in soybean. Genet Mol Res: GMR. 2015;14(2):6101-9.
Lu S, Dong L, Fang C, Liu S, Cheng Q, Kong L, et al. Stepwise selection on homeologous PRR genes controlling flowering and maturity during soybean domestication. Nat Genet. 2020;52(4):1-9.
Takeshima R, Nan HY, Harigai K, Dong LD, Zhu JH, Lu SJ, et al. Functional divergence between soybean FLOWERING LOCUS T orthologues, FT2a and FT5a, in post-flowering stem growth. J Exp Bot. 2019;70(15): 3941-3953.
Cober ER, Voldeng HD. A new soybean maturity and photoperiod-sensitivity locus linked to E1 and T. Crop Sci. 2001;41:698-701.
Hao G, Chen X, Pu M. Classification of the Chinese soybean cultivars into maturity groups. Acta Agronomica Sinica. 1992;18(4):275-81.
Hartwig E. Growth and reproductive characteristics of soybeans [Glycine max (L.) Merr.] grown under short-day conditions. Trop Sci. 1970;12:47-53.
Wang G. The research about the ecological division of soybean cultivars in China. Sci Agr Sin. 1981;39-46.
Zhang LX, Kyei-Boahen S, Zhang J, Zhang MH, Freeland TB, Watson CE, et al. Modifications of optimum adaptation zones for soybean maturity groups in the USA. Crop Manage. 2007;6:1.
Hatfield JL, Prueger JH. Temperature extremes: effect on plant growth and development. Weather and Climate Extremes. 2015;10:4-10.
Kiss T, Dixon LE, Soltesz A, Banyai J, Mayer M, Balla K, et al. Effects of ambient temperature in association with photoperiod on phenology and on the expressions of major plant developmental genes in wheat (Triticum aestivum). Plant Cell Environ. 2017;40(8):1629-42.
Luan WJ, Chen HZ, Fu YP, Si HM, Peng W, Song SS, et al. The effect of the crosstalk between photoperiod and temperature on the heading-date in rice. PloS One. 2009;4(6):e5891.
Tian L, Wang SX, Song XH, Zhang J, Liu P, Chen Z, et al. Long photoperiod affects the maize transition from vegetative to reproductive stages: a proteomic comparison between photoperiod-sensitive inbred line and its recurrent parent. Amino Acids. 2018;50(1):149-61.
Xia Z, Watanabe S, Yamada T, Tsubokura Y, Nakashima H, Zhai H, et al. Positional cloning and characterization reveal the molecular basis for soybean maturity locus E1 that regulates photoperiodic flowering. P Natl Acad Sci Usa. 2012;109(32):E2155-64.
Chen LY, Nan HY, Kong LP, Yue L, Yang H, Zhao QS, et al. Soybean AP1 homologs control flowering time and plant height. J Integr Plant Biol. 2020;00(00):1-12.
Allen LH, Zhang L, Boote KJ, Hauser BA. Elevated temperature intensity, timing, and duration of exposure affect soybean internode elongation, mainstem node number, and pod number per plant. Crop J. 2018;6(2):148-61.
Thomas JF, Raper CD, Jr. Morphological response of soybeans as governed by photoperiod, temperature, and age at treatment. Botanical Gazette. 1977;138(3):321-8.
Van Schaik PH, Probst AH. Effects of some environmental factors on flower production and reproductive efficiency in soybeans. Agron J. 1958;50(4):192-7.
Wang XG, Zhao NL, Wei JJ, Dong Z. Case analysis of super-high-yielding soybean variety，Zhonghuang 35. Soybean science. 2011;30(6):1051-3.
Peoples MB, Brockwell J, Herridge DF, Rochester IJ, Alves BJR, Urquiaga S, et al. The contributions of nitrogen-fixing crop legumes to the productivity of agricultural systems. Symbiosis. 2009;48(1/3):1-17.
Hao T, Zhu Q, Zeng M, Shen J, Shi X, Liu X, et al. Quantification of the contribution of nitrogen fertilization and crop harvesting to soil acidification in a wheat-maize double cropping system. Plant Soil. 2019;434(1-2):167-84.
Lin H, Jing CM, Wang JH. The Influence of long-term fertilization on soil acidification. Advanced Materials Research. 2014;955-9:3552-5.
Schroder JL, Zhang H, Girma K, Raun WR, Penn CJ, Payton ME. Soil acidification from long-term use of nitrogen fertilizers on winter wheat. Soil Sci Soc Am J. 2011;75(3):957-64.
Yang YQ, Tong Y, Li XX, He Y, Xu RN, Liu D, et al. Genetic analysis and fine mapping of phosphorus efficiency locus 1 (PE1) in soybean. Theor Appl Genet. 2019;132: 2847-58.
Chen LY, Qin L, Zhou LL, Chen ZC, Sun LL, Wang WF, et al. A nodule-localized phosphate transporter GmPT7 plays an important role in enhancing symbiotic N2 fixation and yield in soybean. New Phytol. 2019;221(4):2013-25.
Schröder JJ, Smit AL, Cordell D, Rosemarin A. Improved phosphorus use efficiency in agriculture: a key requirement for its sustainable use. Chemosphere. 2011;84(6):822-31.
Knott DR, Kumar J. Comparison of early generation yield testing and a single seed descent procedure in wheat breeding. Crop Sci. 1975;15(3):295-9.
Gray CD, Kinnear PR. IBM SPSS statistics 19 made simple. Am Stat. 2012;66(2):143.
Li ZF, Xu MG, Zhang HM, Zhang SX, Zhang WJ. Sustainability of crop yields in China under long-term fertilization and different ecological conditions. Chinese Journal of Applied Ecology. 2010;21(5):1246-69.
Bao SD. Soil and agricultural chemistry analysis. In: Li GZ, Yang GD, editors. Beijing: China Agriculture Press; 2000. p. 56-108.
Meng L, Li H, Zhang L, Wang J. QTL IciMapping: Integrated software for genetic linkage map construction and quantitative trait locus mapping in biparental populations. Crop J. 2015;3(3):269-83.
Van Ooijen J W. MapQTL 6, software for the mapping of quantitative trait loci in experimental populations of diploid species. Wageningen: Kyazma B.V. 2009.
Brzosko E, Bajguz A. Nectar composition in moth-pollinated Platanthera bifolia and chlorantha and its importance for reproductive success. Planta. 2019;250(1):263-79.
Kassambara A, Mundt F. Factoextra: extract and visualize the results of multivariate data analyses. R Package Version. 2017;1.0. 5. http://spout.ussg.indiana.edu/CRAN/web/packages/factoextra/index.html.
Lê S, Josse J, Husson F. FactoMineR: an R package for multivariate analysis. J Stat Softw. 2008;25, doi:10.18637/jss.v025.i01.

Table 1 Agro-meteorological data and basic soil chemical properties of the experimental locations

Experimental Locations		Agro-meteorological data					Soil basic chemical properties
Experimental Locations		AMaT	AMiT	AT	EAT	AD	pH	AN	AP	AK
ZC	114.48°E, 37.50°N	31.16	20.03	2353.50	1473.50	13.80	8.12	90.25	16.17	89.53
HZ	120.69°E, 30.51°N	29.31	21.56	2515.50	1635.50	12.98	7.85	65.56	14.74	132.42
YZ	118.20°E, 26.17°N	30.28	21.18	2701.50	1651.50	13.50	5.47	129.91	153.17	52.37
BL	114.29°E, 23.17°N	28.13	19.39	1913.50	1133.50	12.65	5.77	78.02	88.89	93.66
HN	109.48°E, 18.31°N	28.15	21.41	2329.50	1389.50	11.11	6.99	69.62	24.74	234.98

Note: ZC, Zhao county; HZ, Hangzhou; YZ, Yangzhong; BL, Boluo; HN, Hainan; AMaT, average maximum temperature, °C; AMiT, average minimum temperature, °C; AT, accumulated temperature, °C·d; EAT, effective accumulated temperature, °C·d; AD, average day-length, hours·day^-1; AN, available nitrogen, mg·kg^-1; AP, available phosphorus, mg·kg^-1; AK, available potassium, mg·kg^-1

Table 2 Phenotypic variation and genetic analysis of plant height traits among 308 soybean germplasm varieties grown in two distinct environments

Trait ^a	Env ^b	Mean ± SD	MIN	MAX	CV%	Kurt	Skew	*h²b*
SH	BL	48.08 ± 14.68	12.00	110.00	30.52	2.42	1.14	0.92
	HN	26.50 ± 6.35	5.00	57.00	23.96	1.34	0.41	0.90
	Total	36.50 ± 15.4	5.00	110.00	42.20	2.37	1.29	0.39
SNN	BL	12.77 ± 3.01	3.00	24.00	23.60	0.27	0.27	0.89
	HN	9.53 ± 1.67	4.00	17.00	17.58	0.81	0.38	0.74
	Total	11.02 ± 2.88	3.00	24.00	26.13	0.81	0.85	0.40
AIL	BL	3.81 ± 0.94	1.23	10.00	24.77	3.90	1.41	0.81
	HN	2.80 ± 0.56	0.56	4.89	19.94	0.87	0.17	0.83
	Total	3.27 ± 0.91	0.56	10.00	27.96	4.06	1.34	0.38

Note: ^a SH, shoot height, cm·plant^-1; SNN, stem node number, #·plant^-1; AIL, average internode length, cm·per node^-1; ^bBL and HN represented the experimental sites of Boluo and Hainan, respectively; Total, includes data from both environments; MIN and MAX, the minimum and maximum values, respectively, of plant height traits; CV%, coefficient of variation; Kurt and Skew represent the Kurtosis and Skewness of plant height traits, respectively; h²b, broad-sense heritability

Table 3 ANOVA for variation of plant height traits among 308 soybean germplasm varieties grown in two distinct environments

Trait	F values
Trait	Environment	Genotype	Environment × Genotype
SH	3420.19***	17.70***	17.67***
SNN	948.83***	15.80***	14.97***
AIL	942.26***	14.93***	12.05***

Note: SH, shoot height; SNN, stem node number; AIL, average internode length; *** indicates significant differences at the 1‰ level (P<0.001)

Table 4 Phenotypic variation and genetic analysis of plant height traits among 168 F₉ soybean RILs grown in four distinct environments

Trait ^a	Env ^b	Parents		RILs ^c
Trait ^a	Env ^b	BX10 ± SD	BD2 ± SD	Mean ± SD	MIN	MAX	CV%	*h²b*
SH	ZC	85.00 ± 5.29	122.00 ± 12.70	112.39 ± 25.51	63.33	191.67	22.70	0.82
	HZ	54.33 ± 4.62	54.67 ± 0.58	55.07 ± 17.15	29.33	102.00	31.14	0.91
	YZ	65.33 ± 4.62	61.00 ± 7.37	72.70 ± 15.32	35.00	141.33	21.08	0.84
	BL	33.56 ± 4.48	36.56 ± 6.67	33.20 ± 6.95	19.11	53.33	20.94	0.93
	Total	50.94 ± 20.53	57.83 ± 32.67	68.35 ± 33.94	19.11	191.67	49.65	0.76
SNN	ZC	22.33 ± 1.53	23.00 ± 1.41	22.70 ± 3.03	13.33	31.33	13.34	0.62
	HZ	12.33 ± 1.53	12.33 ± 1.53	13.73 ± 2.97	7.67	22.00	21.65	0.89
	YZ	13.00 ± 1.00	14.00 ± 2.65	14.22 ± 1.87	8.33	19.33	13.12	0.67
	BL	9.63 ± 1.30	10.11 ± 2.32	10.07 ± 1.35	6.78	14.89	13.37	0.77
	Total	12.94 ± 4.85	12.71 ± 4.62	15.18 ± 5.23	6.78	31.33	34.41	0.72
AIL	ZC	3.81 ± 0.14	5.27 ± 1.03	5.06 ± 1.20	2.69	9.51	23.76	0.72
	HZ	4.45 ± 0.24	4.45 ± 0.49	4.00 ± 0.70	2.75	6.74	17.59	0.70
	YZ	5.05 ± 0.55	4.37 ± 0.44	5.16 ± 0.93	3.50	8.81	18.00	0.67
	BL	3.43 ± 0.34	3.73 ± 0.89	3.36 ± 0.48	2.38	5.00	14.40	0.70
	Total	3.96 ± 0.71	4.15 ± 0.89	4.39 ± 1.15	2.38	9.51	26.16	0.75

Note: ^a SH, shoot height, cm·plant^-1; SNN, stem node number, #·plant^-1; AIL, average internode length, cm·per node^-1; ^bZC, HZ, YZ, and BL represent the experimental sites of Zhao county, Hangzhou, Yangzhong and Boluo, respectively; Total, all the four environments combined; ^cRILs, recombinant inbred lines; MIN and MAX, the minimum and maximum values of plant height traits; CV%, coefficient of variation; h²b, broad-sense heritability

Table 5 ANOVA for variation in plant height traits among 168 F₉soybean RILs grown in four distinct environments

Trait	F values
Trait	Environment	Genotype	Environment × Genotype
SH	5100.59***	15.15***	5.59***
SNN	2223.34***	4.46***	2.62***
AIL	585.38***	5.35***	2.23***

Note: SH, shoot height; SNN, stem node number; AIL, average internode length; *** indicates significant differences at the 1‰ level (P<0.001)

Table 6 Pearson correlation coefficients (r) for relationships among soybean plant height traits, agro-meteorological data, basic soil chemical properties, and additive eﬀects of QTLs in individual environments

	SH	SNN	AIL	*qSHs*	*qSNNs*	*qAILs*
AMaT	0.827**	0.798**	0.602**	0.613**	-0.263**	0.483**
AMiT	ns	ns	0.220**	ns	0.491**	0.252**
AT	0.369**	0.267**	0.472**	0.207**	0.237**	0.411**
EAT	0.368**	0.297**	0.424**	0.244**	0.295**	0.399**
AD	0.821**	0.780**	0.613**	0.591**	-0.318**	0.475**
pH	0.461**	0.595**	0. 092*	0.517**	ns	0.176*
AN	0.280**	0.115**	0.463**	ns	-0.203**	0.279**
AP	-0.250**	-0.415**	0.109**	-0.382**	ns	ns
AK	-0.226**	-0.079*	-0.374**	ns	0.349**	-0.182**

Note: SH, shoot height; SNN, stem node number; AIL, average internode length; AMaT, average maximum temperature; AMiT, average minimum temperature; AT, accumulated temperature; EAT, effective accumulated temperature; AD, average day-length; AN, available nitrogen; AP, available phosphorus; AK, available potassium; qSHs, qSNNs and qAILs represent the sum of additive eﬀects of QTL for soybean on shoot height, stem node number and average internode length in single environments; * and ** indicate significant correlations at the 5% (P<0.05) and 1% (P<0.01) levels, respectively

Additionalfile1Revision20201223.xlsx
Additional file 1: Table S1 Putative QTLs detected for plant height traits by MapQTL 6.0 and QTL IciMaping 4.1 using 168 F9 soybean RILs under different environments. (XLSX 18.8 KB)

PDF

Journal Publication

published 25 Jan, 2021

Read the published version in BMC Plant Biology →

Submission checks completed at journal
31 Dec, 2020
Editorial decision: Minor revision
28 Dec, 2020

You are reading this older preprint version

Read the latest preprint version →

Environmental and Genetic Regulation of Plant Height in Soybean

Status:

Journal Publication

Version 3

Abstract

Figures

Background

Results

Discussion

Conclusions

Methods

Abbreviations

Declarations

References

Tables

Supplementary Files

Status:

Journal Publication

Version 3