# Path Coefficient and Principal Component Analyses for Biomass Allocation, Drought Tolerance and Carbon Sequestration Potential in Wheat

Kwame Shamuyarira, Hussein Shimelis, Sandiswa Figlan, Vincent Chaplot

► To cite this version:

Kwame Shamuyarira, Hussein Shimelis, Sandiswa Figlan, Vincent Chaplot. Path Coefficient and Principal Component Analyses for Biomass Allocation, Drought Tolerance and Carbon Sequestration Potential in Wheat. Plants, 2022, 11 (11), pp.1407. 10.3390/plants1111407. [https://hal.science/hal-03772853v1](https://hal.science/hal-03772853v1)

HAL Id: hal-03772853

Submitted on 9 Sep 2022

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons CC BY 4.0 - Attribution - International License
<figure><img src="image_2.png" /></figure>

Article

# Path Coefficient and Principal Component Analyses for Biomass Allocation, Drought Tolerance and Carbon Sequestration Potential in Wheat

Kwame W. Shamuyarira <fnref n="1" />*, Hussein Shimelis <fnref n="1" />,*, Sandiswa Figlan <fnref n="2" /> and Vincent Chaplot <fnref n="1" />,<fnref n="3" />

**Abstract:** Increased root biomass allocation could serve as a proxy trait for selecting crop ideotypes with drought tolerance and carbon sequestration potential in agricultural soils. The objective of this study was to assess the magnitude of the relationship between root biomass and yield components and to identify influential traits so as to optimise genotype selection for enhanced biomass allocation, drought tolerance and carbon sequestration potential in bread wheat (*Triticum aestivum* L.). One-hundred wheat genotypes consisting of 10 parents and 90 derived F₂ families were evaluated under drought-stressed and non-stressed conditions at two different sites. Data were collected for days to heading (DTH), days to maturity (DTM), plant height, productive tiller number (TN), spike length, spikelets per spike (SPS), kernels per spike (KPS), thousand kernel weight (TKW), shoot biomass, root biomass, total plant biomass (PB), root-to-shoot ratio (RS) and grain yield. There was significant (*p* < 0.05) genetic variation in most assessed traits, TN and RS being exceptions. Root biomass had significant positive correlations with grain yield under drought-stressed (*r* = 0.28) and non-stressed (*r* = 0.41) conditions, but a non-significant correlation was recorded for RS and grain yield. Notably, both root biomass and shoot biomass had significant positive correlations under both water regimes, revealing the potential of increasing both traits with minimal biomass trade-offs. The highest positive direct effects on grain yield were found for KPS and PB under both water regimes. The present study demonstrated that selection based on KPS and PB rather than RS will be more effective in ideotype selection of segregating populations for drought tolerance and carbon sequestration potential.

**Keywords:** biomass allocation; drought tolerance; path analysis; principal component analysis; selection efficiency; wheat

<figure><img src="image_3.png" /></figure>

check for
updates

**Citation:** Shamuyarira, K.W.; Shimelis, H.; Figlan, S.; Chaplot, V. Path Coefficient and Principal Component Analyses for Biomass Allocation, Drought Tolerance and Carbon Sequestration Potential in Wheat. *Plants* **2022**, *11*, 1407. https://doi.org/10.3390/plants11111407

Academic Editors: Amidou N'Diaye and Jemanesh Haile

Received: 9 May 2022

Accepted: 23 May 2022

Published: 26 May 2022

**Publisher's Note:** MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_4.png" /></figure>

**Copyright:** © 2022 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).

## 1. Introduction

Soil organic carbon (SOC) is an essential component of soil health, affecting the chemical, physical and biological properties of agricultural soils. In sub-Saharan Africa (SSA), soils of most agricultural lands are degraded due to high soil carbon loss driven by intensive production, poor soil management and crop residue removal [1]. Crop residues present the largest source of carbon (45%) that can be assimilated into the soil through carbon sequestration to improve SOC content [2]. The source of crop residues are root biomass or above-ground plant biomass. However, subsoil additions of carbon from root residues in the form of root litter and root exudates are more efficient in increasing SOC than above-ground plant residues [3]. Furthermore, the residence time in soil of root-derived carbon is 2.4 times that of carbon derived from shoot biomass [4,5]. Hence, soil carbon derived from root biomass is critical for the long-term stabilisation of soil carbon stocks.

<footnote n="1">¹ School of Agricultural, Earth and Environmental Sciences, University of KwaZulu-Natal, Pietermaritzburg 3201, South Africa; shimelish@ukzn.ac.za (H.S.); vincent.chaplot@ird.fr (V.C.)</footnote>

<footnote n="2">² Department of Agriculture and Animal Health, University of South Africa, Florida 1709, South Africa; figlas@unisa.ac.za</footnote>

<footnote n="3">³ Laboratory of Oceanography and Climate, Experiments and Numerical Approaches (LOCEAN), UMR 7159, IRD/C NRS/UPMC/MNHN, IPSL, 75005 Paris, France</footnote>

<footnote>* Correspondence: 217060194@stu.ukzn.ac.za</footnote>
Promoting below-ground additions of carbon via increased crop root biomass may improve the rate of carbon sequestration in agricultural soils [6].

Breeding wheat (*Triticum aestivum* L., 2n = 6x = 42, AABBDD) with high biomass allocation to roots promises to increase carbon sequestration potential and contribute to yield stability in dryland farming [7]. Crops grown under dryland or rainfed systems are exposed to high fluctuations of water in the soil profile due to rainfall variability and recurrent drought. Plants respond to these fluctuations by altering their biomass allocation patterns [8]. The most unique mechanism by which plants adjust biomass allocation under drought stress is by increasing root-to-shoot ratios, whereby plants allocate more carbohydrates to roots [9]. The accumulation of carbohydrates in the roots stimulates more profound root growth as plants scavenge for moisture to support growth and development [10]. Different studies have shown that root sizes and depths vary significantly within species and between crop cultivars of major cereals, including wheat [11], maize and rice [12]. The genetic variation allows for the identification of contrasting parental genotypes with optimal growth and biomass allocation and agronomic traits to generate new breeding populations for the development of a wheat ideotype with high root biomass and better yield gains.

There is a dearth of information and some debate about whether large root systems impact the increased productivity of wheat in different soil moisture conditions [13,14]. Root biomass has been regarded as an inefficient sink for photo-assimilates, as twice the amount of carbohydrates needed to produce one unit of shoot biomass is needed to produce one unit of root biomass [15]. In high-input crop production systems that use inorganic fertilisers and irrigation, deeper root systems are underestimated, as water and nutrients are readily available [16]. However, recent studies have shown that root system plasticity can be targeted to provide resilience and yield stability under both irrigated and rainfed conditions in arid and semi-arid environments [17–19]. High root biomass increases the capacity for and the efficiency of water capture from the soil profile to support plant growth, biomass production, anthesis, grain filling and final grain yield [20]. Ref. [14] reported a 25% increase in shoot biomass production under full irrigation in wheat cultivars with high root biomass relative to those with small root systems. Positive correlations have been observed between seedling root attributes and yield-promoting traits, such as maturity time, kernels per spike, shoot biomass and grain yield, in wheat [21,22]. Therefore, increasing rooting capacity, yield and yield-influencing agronomic traits could stabilise productivity and ensure sustainable grain production of wheat cultivars in different soil moisture conditions [21,23].

The extent of the relationship between root biomass and yield components can be explored through path coefficient and principal component analyses to identify influential traits for cultivar development. In an attempt to develop next-generation wheat cultivars for drought tolerance and carbon sequestration potential, the research group at the University of KwaZulu-Natal identified lines with high genetic variation for desirable traits, such as root biomass, root-to-shoot ratio and grain yield under drought conditions. The selected lines were crossed to produce new breeding populations which were advanced to the F₂ generation for individual plant selection. Selected plants will be advanced to subsequent generations via pure line or pedigree breeding using complementary technologies, such as doubled haploidy and speed breeding to fast-track selection gains [24,25]. In light of the above background, the objective of this study was to assess the magnitude of the relationship between root biomass and yield components and identify influential traits to optimise selection for enhanced biomass allocation, drought tolerance and carbon sequestration potential in bread wheat. Evaluations of the parental lines and F₂ families were conducted in field and greenhouse trials. It was hypothesized that there are direct or indirect relationships between and variable contributions made by root biomass and above-ground traits, this hypothesis being used as a guide in the selection of unique genotypes with high root biomass and better grain yield gains.
## 2. Results

### 2.1. *Analysis of Variance*

The combined analysis of variance revealed significant differences among the test genotypes for the recorded traits, except for TN and RS (Table 1). There were significant differences for RS, TN, DTH, SPS, TKW, spike length and grain yield due to the genotype-by-water regime interaction effect. Significant differences were recorded for DTH, SPS, KPS, TKW, PB, plant height, spike length, shoot biomass, root biomass and grain yield for the genotype-by-site interaction. The genotype-by-water regime-by-site interaction effect was significant only for spike length, SPS and TKW.

### 2.2. *Mean Performance*

The mean values of measured traits of the ten best performing genotypes and the lowest-performing genotypes based on grain yield under drought-stressed conditions are presented in Table S1. Genotype LM70 × BW141 had the fewest days to heading under drought-stressed and non-stressed conditions. Three of the top four genotypes (LM26 × BW140, BW140 × LM71 and LM47 × BW152) matured later when compared with the lowest-performing genotypes. The top ten genotypes were generally taller than the lowest-performing genotypes, while no clear trends were found for TN for the two groups. Spike length (6.13%) and SPS (5.82%) were reduced slightly by drought stress. Among the lowest-performing genotypes, BW152 (399.15 g·m⁻²) and BW162 × LM71 (334.12 g·m⁻²) had higher shoot biomass under non-stressed conditions than all the top 10 genotypes. However, these values were significantly reduced under drought-stressed conditions. LM75 × LM71 had notably high shoot biomass (>300 g·m⁻²), root biomass (>40 g·m⁻²) and PB (>400 g·m⁻²) values under both moisture conditions but had low yield gain (72.11 g·m⁻²) under drought. The average root-to-shoot ratio for all genotypes increased by 60% due to drought stress. The lowest-performing genotypes had relatively high grain yields (>300 g·m⁻²) under non-stressed conditions. Among the top-performing genotypes, only genotype BW162 × LM75 (436.16 g·m⁻²) had a higher grain yield under non-stressed conditions than all the poor-performing genotypes.

### 2.3. *Principal Component and Biplot Analyses*

The rotated component matrix showing the percentage variance of different principal components (PCs) and the respective loadings of recorded traits is shown in Table 2. The first four PCs under non-stressed conditions had a cumulative variance of 79.93%. The first PC had the highest variation of 37.04%, followed by PC2 with 15.27%. Shoot biomass, PB and grain yield made the highest contributions to PC1, followed by root biomass, RS and TKW with positive contributions to PC2. The highest positive loadings for PC3 and PC4 were for RS and DTM, respectively. Under drought-stressed conditions, the first PC had a percentage variance of 39.19% and was positively correlated with PB, shoot biomass, SPS, spike length and plant height. The second and third PCs were correlated with DTH and TKW, respectively. Principal components 4 and 5 had high positive loadings for root-related traits, such as RS and root biomass, respectively.

Biplots based on the principal component analysis were drawn for drought-stressed and non-stressed conditions (Figures 1 and 2, respectively). The genotypes were evenly scattered across both PC1 and PC2 under drought-stressed conditions. High-yielding genotypes under stress, such as the parental line LM75 and cross BW140 × LM71 had high KPS and PB, while LM26 × BW140 and LM47 × BW152 were also associated with TN (Figure 1 and Table S1). High RS and late flowering were reflected in poor-yielding families, including LM47 × LM71, LM75 × LM71 and BW152 (Figure 1 and Table S1). Under non-stressed conditions, grain yield had strong correlations with PB and TN and high-yielding genotypes under these moisture conditions, including crosses LM71 × LM75 and LM71 × LM26, which excelled in these traits (Figure 2 and Table S1).
<table><caption>**Table 1.** Combined analysis of variance and significance tests for yield components and root attributes of 10 parental lines and their 90 F₂ progenies evaluated under drought-stressed and non-stressed conditions across two sites.</caption><thead><tr><th>SOV</th><th>d.f</th><th>DTH</th><th>DTM</th><th>PH</th><th>TN</th><th>SL</th><th>SPS</th><th>KPS</th><th>TKW</th><th>SB</th><th>RB</th><th>PB</th><th>RS</th><th>GY</th></tr></thead><tbody><tr><td>Replication</td><td>1</td><td>30.06</td><td>228.44 ***</td><td>142.45</td><td>335.34 ***</td><td>5.95 *</td><td>18.01 *</td><td>212.15</td><td>421.16 ***</td><td>4484.00</td><td>10,390.90 ***</td><td>220,237.00 *</td><td>0.49 ***</td><td>43,869.00 *</td></tr><tr><td>Block</td><td>18</td><td>29.41 ***</td><td>27.11</td><td>176.89 *</td><td>21.30</td><td>2.06 *</td><td>6.352</td><td>115.63 *</td><td>56.25 **</td><td>25,982.00 ***</td><td>994.50 *</td><td>67,304.00</td><td>0.04</td><td>21,502.00 **</td></tr><tr><td>Genotype (Gen)</td><td>99</td><td>60.16 ***</td><td>41.60 ***</td><td>243.27 ***</td><td>16.10</td><td>4.09 ***</td><td>11.61 ***</td><td>142.25 ***</td><td>58.08 ***</td><td>18,128.00 ***</td><td>769.50 **</td><td>69,654.00 ***</td><td>0.02</td><td>14,638.00 **</td></tr><tr><td>Water Regime (WR)</td><td>1</td><td>36.23</td><td>11,794.04 ***</td><td>11,091.44 ***</td><td>205.42 **</td><td>59.49 ***</td><td>202.55 ***</td><td>7189.67 ***</td><td>9737.33 ***</td><td>938,657.00 ***</td><td>10,709.50 ***</td><td>12,251,041.00 ***</td><td>1.68 ***</td><td>5,465,271.00 ***</td></tr><tr><td>Site</td><td>1</td><td>46,292.11 ***</td><td>128,227.14 ***</td><td>89,648.69 ***</td><td>15,460.90 ***</td><td>1630.17 ***</td><td>5094.70 ***</td><td>67,179.66 ***</td><td>64,551.23 ***</td><td>2,054,036.00 ***</td><td>24,588.50 ***</td><td>1,132,252.00 ***</td><td>0.05</td><td>165,356.00 ***</td></tr><tr><td>Gen * WR</td><td>99</td><td>13.34 *</td><td>20.84</td><td>104.61</td><td>20.34 *</td><td>1.59 *</td><td>7.20 ***</td><td>71.70</td><td>32.42 *</td><td>11,642.00</td><td>631.90</td><td>53,980.00</td><td>0.03 *</td><td>13,481.00 *</td></tr><tr><td>Gen * Site</td><td>98</td><td>37.36 ***</td><td>17.86</td><td>149.27 ***</td><td>17.57</td><td>2.28 ***</td><td>10.32 ***</td><td>94.62 *</td><td>32.11 *</td><td>17,491.00 ***</td><td>756.10 *</td><td>67,155.00 **</td><td>0.03</td><td>14,496.00 *</td></tr><tr><td>Gen.WR Site</td><td>97</td><td>12.37</td><td>17.73</td><td>116.64</td><td>13.40</td><td>1.73 **</td><td>7.13 ***</td><td>72.21</td><td>37.09 ***</td><td>12,166.00</td><td>577.70</td><td>45,619.00</td><td>0.03</td><td>12,154.00</td></tr><tr><td>Residual</td><td>368</td><td>9.79</td><td>19.02</td><td>91.79</td><td>14.67</td><td>1.148</td><td>4.352</td><td>66.81</td><td>23.58</td><td>10,071.00</td><td>533.00</td><td>43,026.00</td><td>0.03</td><td>10,167.00</td></tr></tbody></table>

<table><caption>**Table 2.** Rotated component matrix for 13 yield components and root attributes in 100 wheat genotypes assessed under drought-stressed and non-stressed conditions.</caption><thead><tr><th rowspan="2">Traits</th><th colspan="4">Non-Stress</th><th colspan="4">Drought-Stress</th></tr><tr><th>PC1</th><th>PC2</th><th>PC3</th><th>PC4</th><th>PC1</th><th>PC2</th><th>PC3</th><th>PC4</th><th>PC5</th></tr></thead><tbody><tr><td>DTH</td><td>0.08</td><td>-0.38</td><td>0.17</td><td>-0.33</td><td>0.05</td><td>0.44</td><td>0.35</td><td>0.22</td><td>0.33</td></tr><tr><td>DTM</td><td>0.06</td><td>-0.29</td><td>0.01</td><td>0.45</td><td>-0.11</td><td>-0.49</td><td>0.18</td><td>0.34</td><td>0.27</td></tr><tr><td>PH</td><td>0.32</td><td>0.00</td><td>0.25</td><td>-0.29</td><td>0.35</td><td>0.07</td><td>0.33</td><td>0.07</td><td>0.10</td></tr><tr><td>TN</td><td>0.29</td><td>0.25</td><td>-0.27</td><td>0.04</td><td>0.26</td><td>-0.09</td><td>-0.30</td><td>-0.25</td><td>0.26</td></tr><tr><td>SL</td><td>0.32</td><td>-0.24</td><td>0.36</td><td>-0.12</td><td>0.35</td><td>0.13</td><td>0.26</td><td>0.24</td><td>-0.16</td></tr><tr><td>SPS</td><td>0.34</td><td>-0.29</td><td>0.24</td><td>0.07</td><td>0.35</td><td>0.11</td><td>0.14</td><td>0.25</td><td>-0.21</td></tr><tr><td>KPS</td><td>0.26</td><td>-0.37</td><td>0.03</td><td>0.22</td><td>0.27</td><td>-0.16</td><td>0.00</td><td>0.14</td><td>-0.60</td></tr><tr><td>TKW</td><td>0.04</td><td>0.31</td><td>0.24</td><td>-0.56</td><td>-0.18</td><td>-0.37</td><td>0.42</td><td>0.20</td><td>0.25</td></tr><tr><td>SB</td><td>0.39</td><td>0.07</td><td>-0.28</td><td>-0.05</td><td>0.39</td><td>0.06</td><td>0.05</td><td>-0.20</td><td>0.23</td></tr><tr><td>RB</td><td>0.22</td><td>0.43</td><td>0.28</td><td>0.28</td><td>0.27</td><td>0.04</td><td>-0.43</td><td>0.25</td><td>0.39</td></tr><tr><td>PB</td><td>0.41</td><td>0.16</td><td>-0.21</td><td>0.04</td><td>0.39</td><td>-0.22</td><td>-0.06</td><td>-0.09</td><td>0.19</td></tr><tr><td>RS</td><td>-0.04</td><td>0.31</td><td>0.60</td><td>0.37</td><td>-0.08</td><td>0.07</td><td>-0.44</td><td>0.69</td><td>-0.03</td></tr><tr><td>GY</td><td>0.39</td><td>0.14</td><td>-0.19</td><td>0.05</td><td>0.23</td><td>-0.55</td><td>-0.05</td><td>-0.02</td><td>-0.10</td></tr><tr><td>Eigenvalue</td><td>4.82</td><td>1.99</td><td>1.40</td><td>1.35</td><td>5.10</td><td>1.79</td><td>1.52</td><td>1.23</td><td>1.02</td></tr><tr><td>Per Var (%)</td><td>37.04</td><td>15.27</td><td>10.74</td><td>10.35</td><td>39.19</td><td>13.77</td><td>11.66</td><td>9.45</td><td>7.81</td></tr><tr><td>Cum Var (%)</td><td>37.04</td><td>52.31</td><td>63.05</td><td>73.40</td><td>39.19</td><td>52.96</td><td>64.62</td><td>74.07</td><td>81.88</td></tr></tbody></table>

<footnote>* Significant at *p* < 0.05, ** *p* < 0.01, *** *p* < 0.001, SOV = source of variation, df = degrees of freedom, DTH = days to 50% heading, DTM = days to 50% maturity, PH = plant height (cm), TN = tiller number, SL = spike length (cm), SPS = spikelets per spike, KPS = kernels per spike, TKW = thousand kernel weight (g), SB = shoot biomass (g·m⁻²), RB = root biomass (g·m⁻²), PB = total plant biomass (g·m⁻²), RS = root-to-shoot ratio and GY = grain yield (g·m⁻²).</footnote>

<footnote>PC = principal component, Per Var = percentage variance, Cum Var = cumulative variance, DTH = days to 50% heading, DTM = days to 50% maturity, PH = plant height, TN = tiller number, SL = spike length, SPS = spikelets per spike, KPS = kernels per spike, TKW = thousand kernel weight, SB = shoot biomass, RB = root biomass, PB = total plant biomass, RS = root-to-shoot ratio and GY = grain yield.</footnote>
<figure><img src="image_5.png" /><figcaption>**Figure 1.** Principal component biplot of 10 bread wheat parental lines and 90 F₂ families under drought-stressed conditions. Genotypes are coded with numbers, as recorded in Table S1. Smaller angles between vectors of recorded traits indicate a high correlation between the traits in discriminating genotypes. Genotypes plotted closer to and further along a vector line scored highly in that particular trait. PC = principal component, DTH = days to 50% heading, DTM = days to 50% maturity, PH = plant height (cm), TN = tiller number, SL = spike length (cm), SPS = spikelets per spike, KPS = kernels per spike, TKW = thousand kernel weight (g), SB = shoot biomass (g·m⁻²), RB = root biomass (g·m⁻²), PB = total plant biomass (g·m⁻²), RS = root-to-shoot ratio and GY = grain yield (g·m⁻²).</figcaption></figure>

## 2.4. Correlations of Root Biomass and Yield Components with Grain Yield under Drought-Stressed and Non-Stressed Conditions

Phenotypic correlation coefficients showing the relationships between root attributes and yield components under drought-stressed and non-stressed conditions are shown in Table 3. Grain yield showed higher positive correlations with PB (*r* = 0.67) and KPS (*r* = 0.55) under drought-stressed conditions. Productive tiller number and shoot biomass had moderate correlations of *r* = 0.33 and *r* = 0.39, respectively, with grain yield, whereas root biomass had low (*r* < 0.30) positive correlations with grain yield. Only DTH had significant negative correlations (*r* = −0.38) with grain yield. Root biomass was positively correlated with RS and all above-ground traits except DTH and DTM. Negative correlations were recorded for RS with plant height (*r* = −0.24), shoot biomass (*r* = −0.31) and PB (*r* = −0.19).
<figure><img src="image_6.png" /></figure>

**Figure 2.** Principal component biplot of 10 bread wheat parental lines and 90 F₂ families under non-stressed conditions. Genotypes are coded with numbers, as recorded in Table S1. Smaller angles between vectors of recorded traits indicate a high correlation between the traits in discriminating genotypes. Genotypes plotted closer to and further along a vector line scored highly in that particular trait. PC = principal component, DTH = days to 50% heading, DTM = days to 50% maturity, PH = plant height (cm), TN = tiller number, SL = spike length (cm), SPS = spikelets per spike, KPS = kernels per spike, TKW = thousand kernel weight (g), SB = shoot biomass (g·m⁻²), RB = root biomass (g·m⁻²), PB = total plant biomass (g·m⁻²), RS = root-to-shoot ratio and GY = grain yield (g·m⁻²).

## 2.5. Path Coefficient Analysis of Root Biomass and Yield Components on Grain Yield

Direct and indirect effects of yield components and root attributes on grain yield under drought-stressed and non-stressed conditions are summarised in Table S2 and Figure 3 and Table S3 and Figure 4, respectively. PB (1.001) had the greatest direct effect on grain yield under drought-stressed conditions, followed by KPS (0.34) and TKW (0.17). Conversely, shoot biomass and root biomass had negative direct effects on grain yield at −0.410 and −0.153, respectively. Shoot biomass (0.88) and root biomass (0.65) had the greatest positive indirect effects on grain yield through PB. Indirect effects on grain yield through PB were observed for plant height, TN, spike length, SPS and KPS. Root-to-shoot ratio had negative direct effects on grain yield through PB.
<table><caption>**Table 3.** Pearson's correlation coefficients and significance tests for yield components and root attributes of 10 bread wheat parental lines and 90 F₂ progenies under drought-stressed (lower diagonal values) and non-stressed (upper diagonal values) conditions at two sites.</caption><thead><tr><th>Traits</th><th>DTH</th><th>DTM</th><th>PH</th><th>TN</th><th>SL</th><th>SPS</th><th>KPS</th><th>TKW</th><th>SB</th><th>RB</th><th>PB</th><th>RS</th><th>GY</th></tr></thead><tbody><tr><td>DTH</td><td>1</td><td>0.03</td><td>0.23*</td><td>-0.05</td><td>0.21*</td><td>0.13</td><td>0.12</td><td>0.02</td><td>0.1</td><td>-0.16</td><td>0.03</td><td>-0.17</td><td>0.02</td></tr><tr><td>DTM</td><td>-0.29**</td><td>1</td><td>-0.01</td><td>-0.01</td><td>0.05</td><td>0.32**</td><td>0.05</td><td>-0.21*</td><td>0.11</td><td>0.03</td><td>0.03</td><td>-0.01</td><td>0.02</td></tr><tr><td>PH</td><td>0.28**</td><td>-0.25*</td><td>1</td><td>0.32**</td><td>0.55***</td><td>0.51***</td><td>0.28**</td><td>0.32**</td><td>0.53***</td><td>0.31**</td><td>0.51***</td><td>-0.05</td><td>0.48***</td></tr><tr><td>TN</td><td>-0.04</td><td>-0.1</td><td>0.30**</td><td>1</td><td>0.22*</td><td>0.24*</td><td>0.15</td><td>0.05</td><td>0.59***</td><td>0.42***</td><td>0.62***</td><td>0.03</td><td>0.60***</td></tr><tr><td>SL</td><td>0.27**</td><td>-0.27**</td><td>0.74***</td><td>0.31**</td><td>1</td><td>0.76***</td><td>0.49***</td><td>0.1</td><td>0.45***</td><td>0.22*</td><td>0.41***</td><td>-0.06</td><td>0.41***</td></tr><tr><td>SPS</td><td>0.13</td><td>-0.15</td><td>0.64***</td><td>0.32**</td><td>0.79***</td><td>1</td><td>0.62***</td><td>-0.04</td><td>0.49***</td><td>0.2</td><td>0.47***</td><td>-0.07</td><td>0.48***</td></tr><tr><td>KPS</td><td>-0.04</td><td>-0.06</td><td>0.39***</td><td>0.25*</td><td>0.49***</td><td>0.55***</td><td>1</td><td>-0.25*</td><td>0.29**</td><td>0.04</td><td>0.40***</td><td>-0.11</td><td>0.48***</td></tr><tr><td>TKW</td><td>-0.13</td><td>0.37***</td><td>-0.17**</td><td>-0.22*</td><td>-0.24*</td><td>-0.38***</td><td>-0.29**</td><td>1</td><td>0.05</td><td>0.09</td><td>0.1</td><td>0.08</td><td>0.13</td></tr><tr><td>SB</td><td>0.15</td><td>-0.23*</td><td>0.71***</td><td>0.55***</td><td>0.63***</td><td>0.59***</td><td>0.41***</td><td>-0.40***</td><td>1</td><td>0.37***</td><td>0.84***</td><td>-0.24*</td><td>0.72***</td></tr><tr><td>RB</td><td>0.06</td><td>-0.03</td><td>0.37***</td><td>0.51***</td><td>0.36***</td><td>0.40***</td><td>0.23*</td><td>-0.33***</td><td>0.58***</td><td>1</td><td>0.49***</td><td>0.65***</td><td>0.41***</td></tr><tr><td>PB</td><td>-0.06</td><td>-0.03</td><td>0.60***</td><td>0.57***</td><td>0.54***</td><td>0.58***</td><td>0.50***</td><td>-0.26**</td><td>0.88***</td><td>0.64***</td><td>1</td><td>-0.05</td><td>0.94***</td></tr><tr><td>RS</td><td>0.03</td><td>0.1</td><td>-0.24*</td><td>-0.06</td><td>-0.13</td><td>-0.13</td><td>-0.11</td><td>-0.05</td><td>-0.31**</td><td>0.39***</td><td>-0.19***</td><td>1</td><td>-0.03</td></tr><tr><td>GY</td><td>-0.38***</td><td>0.24*</td><td>0.26**</td><td>0.33***</td><td>0.24*</td><td>0.25*</td><td>0.55***</td><td>0.09</td><td>0.39***</td><td>0.28**</td><td>0.67***</td><td>-0.14</td><td>1</td></tr></tbody></table>

<figure><img src="image_7.png" /><figcaption>**Figure 3.** Path analysis model diagram displaying causal relationships of assessed traits on grain yield in 10 bread wheat parental lines and 90 derived F₂ families of wheat assessed under drought-stressed conditions. Values in parenthesis are direct path coefficients, while other values are correlation coefficients. R = residual effect; DTH = days to 50% heading, DTM = days to 50% maturity, PH = plant height (cm), TN = tiller number, SL = spike length (cm), SPS = spikelets per spike, KPS = kernels per spike, TKW = thousand kernel weight (g), SB = shoot biomass (g·m⁻²), RB = root biomass (g·m⁻²), PB = total plant biomass (g·m⁻²), RS = root-to-shoot ratio and GY = grain yield (g·m⁻²).</figcaption></figure>

<footnote>* Significant at *p* < 0.05, ** *p* < 0.01, *** *p* < 0.001, DTH = days to 50% heading, DTM = days to 50% maturity, PH = plant height (cm), TN = tiller number, SL = spike length (cm), SPS = spikelets per spike, KPS = kernels per spike, TKW = thousand kernel weight (g), SB = shoot biomass (g·m⁻²), RB = root biomass (g·m⁻²), PB = total plant biomass (g·m⁻²), RS = root-to-shoot ratio and GY = grain yield (g·m⁻²). The majority of the traits had significant positive correlations with grain yield under non-stressed conditions. Highly significant correlations were recorded between grain yield and PB (r = 0.94) and between shoot biomass (r = 0.72) and TN (0.60). Plant height (r = 0.48), SPS (r = 0.48), KPS (r = 0.48), root biomass (r = 0.41) and spike length (r = 0.41) showed moderate correlations with grain yield. Productive tiller number and SL had positive correlations with shoot biomass, PB and root biomass. However, spike characteristics, including SPS and KPS, were positively correlated with shoot biomass and PB but had non-significant correlations with root biomass. Root-to-shoot ratio had significant correlations with root biomass (r = 0.65) and shoot biomass (r = -0.24). Positive correlations were recorded for root biomass with PH, TN, spike length, shoot biomass and PB. Notably, under both water regimes, RS had non-significant and negative correlations with grain yield.</footnote>
<figure><img src="image_8.png" /></figure>

**Figure 4.** Path analysis model diagram displaying causal relationships of assessed traits on grain yield in 10 bread wheat parental lines and 90 derived F₂ families of wheat assessed under non-stressed conditions. Values in parenthesis are direct path coefficients, while other values are correlation coefficients. R = residual effect; DTH = days to 50% heading, DTM = days to 50% maturity, PH = plant height (cm), TN = tiller number, SL = spike length (cm), SPS = spikelets per spike, KPS = kernels per spike, TKW = thousand kernel weight (g), SB = shoot biomass (g·m⁻²), RB = root biomass (g·m⁻²), PB = total plant biomass (g·m⁻²), RS = root-to-shoot ratio and GY = grain yield (g·m⁻²).

Under non-stressed conditions, PB (1.06) and KPS (0.13) had positive direct effects on grain yield. Shoot biomass (–0.19) and root biomass (–0.13) had negative direct effects on grain yield; however, the same traits had large positive indirect effects of 0.89 and 0.52, respectively, on grain yield through PB. Except for RS, assessed traits had positive indirect effects on grain yield through PB, including plant height (0.55), TN (0.66), spike length (0.44), SPS (0.50) and KPS (0.43). The residual value for the path analysis model was 0.077 under drought-stressed conditions and 0.231 under non-stressed conditions (Figures 3 and 4).

**3. Discussion**

Genotypic differences observed for yield components and root biomass under differ-
ent environments indicate sufficient genetic variation for the development of new wheat
ideotypes for grain production and carbon sequestration. The rotated component ma-
trix showed that PB, shoot biomass and grain yield were the most discriminating traits
under non-stressed conditions, followed by root biomass, RS and TKW. Therefore, iden-
tifying genotypes based on biomass production and allocation will allow more effective
explanation of the differences among individual genotypes and families. Similarly, under
drought-stressed conditions, biomass production, as shown by high contributions of PB
and shoot biomass and spike-related traits, such as SPS and spike length, were the most
essential traits in differentiating the test genotypes.

Positive correlations were found between shoot biomass and grain yield under both drought-stressed and non-stressed conditions, indicating the importance of shoot biomass for improving grain yield. High shoot biomass contributes to grain yield by providing greater leaf area for light interception and carbon assimilation to support grain filling [26,27]. This agrees with [28], the authors of which reported that yield gains in wheat are influenced by biomass production and sink size, which is determined by the number of fertile spikes
per unit area. To increase the number of fertile spikes per unit area, breeders need to identify plant genotypes with high TN [29]. The importance of TN was notable under drought-stressed conditions, as high yielding genotypes such as LM75 and BW140 × LM71 were strongly correlated with TN. Therefore, increasing shoot biomass and TN will have a positive impact on grain yield.

The positive relationship of root biomass and grain yield shows the importance of root traits for increasing productivity under both drought-stressed and non-stressed conditions. Increase in root size improves the capacity and efficiency of nutrient and moisture acquisition by plants, increasing the resilience of agro-ecosystems. This is more important under drought conditions, where there is limited water in the soil profile and a deeper and more abundant root system can forage for water in the soil profile [30]. Though larger root systems in crops are important, especially in dry areas, they may, in some cases, be inefficient or even entail a yield penalty in wet seasons or in regions with sufficient water and capacity to supply supplementary irrigation [31,32]. However, evidence from this study reflects that root biomass has a positive influence on productivity, even under conditions of sufficient water availability. This is supported by the authors of [33], who reported that deeper and more profuse root growth can be attained without any adverse effect on grain yield. These increased root sizes can be harnessed to increase soil carbon in agricultural soils.

Root-to-shoot ratio was negatively correlated with plant height, shoot biomass and PB under drought, which agrees with [34], whose authors observed negative correlations between root-to-shoot ratio and plant height and above-ground biomass. Some studies debated this and have reported carbon tradeoffs between root and shoot biomass, leading to reduced wheat productivity [35–37]. Therefore, increasing RS alone in pursuit of higher carbon inputs may have a negative impact on grain yield, as reflected in families such as LM47 × LM71, LM75 × LM71 and parental line BW152, which had high RS values but low yields. On the other hand, PB had the highest positive correlations with grain yield under both water regimes. Hence, instead of altering RS, yield gains may be obtained by increasing PB without compromising the carbon sequestration potential from the roots. Increasing root biomass and shoot biomass simultaneously to achieve high PB is possible, as reflected by positive correlations between root biomass and shoot biomass under drought-stressed and non-stressed conditions. However, shoot biomass showed stronger indirect effects on grain yield through PB than root biomass under both water regimes. This may indicate the greater contribution to grain yield by shoot biomass than root biomass under different soil moisture environments. Increasing shoot biomass in such cases rather than root biomass may be more beneficial for yield increases but will lead to a considerable reduction in potential carbon inputs in the soil due to reduced root sizes [3,38].

KPS had high correlations with and direct effects on grain yield under both water regimes, corroborating the report of [39]. According to [40], KPS determines the level of grain yield that is achieved in wheat and understanding the fundamental mechanisms that influence KPS will be useful in increasing grain yield. Drought stress reduces KPS due to kernel abortion, which reduces the sink strength in the spikes and causes yield loss [41]. Thus, maintaining a high kernel number, as in the families LM47 × BW152 and LM70 × BW141, under drought may have a positive influence on grain yield. This explains the close relationship between KPS and dry matter accumulation in the spikes, which is reflected in the strong correlations of KPS with spike length and SPS under both drought-stressed and non-stressed conditions [40,42].

The path coefficient analysis models under NS conditions had lower residual values than under drought-stressed conditions. Therefore, the model with non-stressed conditions was more effective in explaining the total variation in grain yield than that with drought-stressed conditions. Both models showed that PB and KPS can be used for efficient selection of grain yield under different soil moisture conditions while pursuing the goal of increasing carbon inputs in croplands. According to [43], SOC levels in the soil are determined by land management practices, such as crop residue retention. Selecting plants with high PB will
be beneficial at crop farms that practice minimum tillage and other conservation practices; shoot biomass will provide more crop residues to incorporate into the soil profile after harvest. On the other hand, larger root systems in the same plants will make a high carbon contribution to croplands, especially in areas where conventional agriculture is practised and where there is little contribution from above-ground plant residues.

# 4. Materials and Methods

## 4.1. Plant Material and Population Development

The study used 10 contrasting wheat parental lines selected based on their genetic variation for grain yield and better shoot and root biomass production under drought conditions. The names and pedigrees of the parental genotypes are summarised in Table 4. Eight of the lines were candidate drought- and heat-tolerant lines acquired from the International Maize and Wheat Improvement Center's (CIMMYT) drought and heat nurseries, and the two lines were local checks adapted to dryland wheat production. The parents were crossed in controlled environment facilities (CEFs) at the University of KwaZulu-Natal (UKZN), South Africa. The parental lines were stagger-planted on three occasions with two-week intervals to allow for synchronised flowering for emasculation and pollination. A total of 90 F1 crosses were developed, including 45 direct crosses and 45 reciprocals using a full diallel mating design. The F1 seeds were harvested and bulked and advanced to the F2 generation. The 10 parental lines and their 90 F2 families were evaluated under field and greenhouse conditions, as described below.

<table><caption>**Table 4.** List of the parental lines used in the crosses.</caption><thead><tr><td>Genotype</td><td>Pedigree</td></tr></thead><tbody><tr><td>LM26</td><td>ATTILA * 2/PBW65//TAM200/TUI</td></tr><tr><td>LM47</td><td>FRET2/KUKUNA//FRET2/3/YANAC/4/FRET2/KIRITATI</td></tr><tr><td>LM48</td><td>FRET2/KUKUNA//FRET2/3/PASTOR//HXL7573/2 * BAU/5/FRET2 *2/4/SNI/TRAP#1/3/KAUZ * 2/TRAP//KAUZ</td></tr><tr><td>LM71</td><td>BABAX/3/PRL/SARA//TSI/VEE#5/4/CROC_1/AE.SQUARROSA (224)//2 * OPATA</td></tr><tr><td>LM75</td><td>BUC/MN72253//PASTOR</td></tr><tr><td>BW141</td><td>CGSS05B00243T-099TOPY-099M-099NJ-099NJ-1WGY-0B</td></tr><tr><td>BW152</td><td>CGSS05B00258T-099TOPY-099M-099NJ-1WGY-0B</td></tr><tr><td>BW162</td><td>CGSS05B00304T-099TOPY-099M-099NJ-099NJ-3WGY-0B</td></tr><tr><td>LM70</td><td>Local check</td></tr><tr><td>BW140</td><td>Local check</td></tr></tbody></table>

## 4.2. Phenotyping for Root Biomass and Yield Components

### 4.2.1. Field Evaluation

* Experimental design and planting

The test genotypes (10 parental lines and 90 F2 families) were planted in the field at Ukulinga Research Farm of the University of KwaZulu-Natal (29°40' S, 30°24' E; 806 m above sea level). The trial was laid out in a 10 × 10 alpha lattice design with two replications. Genotypes were planted on a 2 m-long plot spaced 30 cm apart, and plants were spaced 20 cm within the row. Three seeds were planted per hole in the experimental unit and thinned to two plants per hill after two weeks of germination. Border rows were used to reduce the risk of yield inflation in test plots in the outer rows. Water and fertiliser were applied through an automated drip irrigation system to ensure that all plants received equal amounts of water.

* Imposing drought stress

The trial was conducted under two water regimes (drought-stressed and non-stressed conditions). Drought-stress was imposed at heading growth stage [44] by withholding irrigation until soil moisture dropped to 35% field capacity, followed by re-watering to 80%

<footnote>* backcross.</footnote>
field capacity. The non-stressed treatment received water by an automated drip irrigation system to maintain the soil at 80% field capacity. The soil moisture content was monitored by digital moisture sensors (HOBO UX120, Onset, Bourne, MA, USA) inserted to a depth of 60 cm and placed in each replicate of the different water regimes. Standard agronomic practices were used in both water regimes for the duration of the trials according to wheat production guidelines in South Africa [45]. Weather data were recorded for the duration of growth during the trial and are summarised in Table 5.

<table><caption>**Table 5.** Monthly weather data during the field trial at Ukulinga Research Farm, Pietermaritzburg, South Africa (2020).</caption><thead><tr><td>Month</td><td>Max Temp (°C)</td><td>Min Temp (°C)</td><td>Humidity (%)</td><td>Rain (mm)</td></tr></thead><tbody><tr><td>July</td><td>19.9</td><td>5.4</td><td>60</td><td>31</td></tr><tr><td>August</td><td>21.8</td><td>7.7</td><td>60</td><td>37</td></tr><tr><td>September</td><td>23.1</td><td>10</td><td>67</td><td>59</td></tr><tr><td>October</td><td>23.3</td><td>12</td><td>75</td><td>100</td></tr><tr><td>November</td><td>23.7</td><td>13.5</td><td>79</td><td>121</td></tr><tr><td>December</td><td>24.8</td><td>15.3</td><td>81</td><td>137</td></tr></tbody></table>

Max Temp = average maximum temperature, Min Temp = average minimum temperature.

### 4.2.2. Greenhouse Evaluation

A greenhouse trial was established at the CEF at the UKZN in Pietermaritzburg, South Africa. The temperature inside the greenhouse was maintained at 25 °C during the day and 15 °C at night. The experimental design followed the field evaluation. Plants were grown in 10 L-capacity plastic pots filled with composted pine bark. Eight seeds were planted in each pot and thinned to five seedlings at two weeks after germination. Water and fertiliser were supplied as described above.

* Imposing drought stress

The trial was conducted under drought-stressed and non-stressed conditions. Irrigation was withheld at the heading stage in the drought-stressed treatment to impose water stress. Rewatering to 80% of field capacity was initiated when moisture dropped to 35% of field capacity to avoid total crop failure. The non-stressed treatment received irrigation four times a day until crop maturity to maintain growing media at 80% of field capacity. Aphids and powdery mildew were chemically controlled using Chess (active ingredient: pyridine azomethine) and Tilt (active ingredient: triazole), respectively. Any emerging weeds were removed manually.

### 4.3. Data Collection

During both field and greenhouse evaluations, data were collected on root biomass, biomass traits and yield components. Days to heading were recorded as the number of days from planting to when 50% of the genotypes in each test plot had fully exposed spikes; days to maturity were recorded as the number of days from planting to when 50% of the genotypes in each test plot had dried spikes; plant height was measured as the height of the plant from the soil to the tip of the spike in centimeters (cm); productive tiller number was recorded as the number of tillers with spikes containing grain at harvest; spike length was the length of the spike, excluding the awns, and expressed in cm; the number of spikelets per spike was counted as the number of spikelets per individual spike; the number of kernels per spike was counted and represented the number of kernels harvested from an individual spike; and thousand kernel weight was taken as the weight of 1000 kernels randomly selected from a genotype and weighed in grams (g). Shoot biomass was recorded as the total above-ground biomass cut from the base of the plant, excluding the grain. The shoots were oven-dried at 70 °C for 48 h, weighed and expressed in g·m⁻². Root biomass was recorded as the total root dry matter harvested per genotype per plot. Root samples for each plot were harvested to a depth of 50 cm. Large roots were separated from the soil by hand and washed under running water to remove all soil particles. The remaining soil was
mixed with water and the suspension was sieved through a 2 mm sieve. Fine roots were
collected from the sieve residue and added to the large roots. The roots were oven-dried
at 70 °C for 48 h. The dried roots were weighed on a balance to obtain the root biomass,
which was adjusted to g·m⁻². Total plant biomass was the sum of all dry plant material for
each genotype, including root biomass, shoot biomass and grain yield, harvested from the
test plots and recorded in g·m⁻². Root-to-shoot ratio was the ratio of the root biomass to
the shoot biomass, as recorded above. Grain yield was the weight of harvested grain at
12.5% moisture content per genotype per plot and expressed in g·m⁻².

4.4. *Data Analysis*

The data collected across the two sites were individually analysed and subjected
to Bartlet's homogeneity test of variance. A combined analysis of variance (ANOVA)
was computed to test genotype effects using a lattice procedure involving a three-way
interaction of genotype, water regime and site, using Genstat, 18th edition [46]. Mean values
of the test genotypes for the assessed traits were compared at the 5% significance level using
Fisher's least significance difference (LSD) procedure. Pearson's correlation coefficients (r)
were calculated separately for drought-stressed and non-stressed conditions using Genstat,
18th edition [46]. Correlation coefficients were partitioned into direct and indirect effects
using the phenotypic correlation matrix: A = B × C, where A is the phenotypic correlation
coefficients of grain yield with yield components and root attributes, B is the phenotypic
correlations among all the recorded traits in all combinations and C is the value of the
path coefficients. The inverse of matrix B was calculated in Microsoft Excel 2016, using
the Matrix Inverse function (MINVERSE). The direct and indirect path coefficients of yield
components and root biomass were calculated as the products of phenotypic correlations
and inverse matrix B, according to [47]. Path coefficient analysis diagrams were developed
in R software using structural equation modelling [48]. The rotated component matrix
and principal component analysis (PCA) biplots were generated for yield components and
root biomass based on the correlation matrix in Genstat, 18th edition. Both analyses were
separated for drought-stressed and non-stressed conditions.

**5. Conclusions**

Increasing RS was not effective with respect to improving wheat productivity. The positive relationship between root biomass and shoot biomass informs the possibility of increasing both traits with minimum biomass trade-offs to develop drought-tolerant cultivars with high carbon sequestration potential. The large positive direct effects of KPS and PB on grain yield under both water regimes point to plants with high root and shoot biomass production as the crop ideotype that will be most resilient and environmentally friendly under varying moisture conditions. This was observed for high-yielding genotypes, including LM75 and BW140 × LM71 under drought-stress and LM71 × LM75 and LM71 × LM26 under non-stressed conditions. Therefore, selection based on KPS and PB rather than RS will be more effective in ideotype selection of segregating F₂ populations with enhanced drought tolerance and carbon sequestration potential. Selected plants and families will be advanced using pedigree breeding to develop climate-smart cultivars for dryland farming in SSA. The study recommends further evaluations across multiple environments and test populations to develop genetic models and to guide the breeding of climate-smart wheat varieties with drought tolerance and high carbon sequestration potential.

**Supplementary Materials:** The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/plants11111407/s1, Table S1: Mean values of yield components and root attributes of the ten parental lines and 90 F₂ families of wheat based on grain yield under drought-stressed conditions. Table S2. Direct (bold face values) and indirect effects for yield components and root attributes on grain yield of 10 bread wheat parental lines and 90 F₂ progenies under drought-stressed conditions at two sites. Table S3. Direct (bold face values) and indirect effects for yield components and root attributes on grain yield of 10 bread wheat parental lines and 90 F₂ progenies under non-stressed conditions at two sites.
**Author Contributions:** Each author participated sufficiently in the completion of this work. Conceptualization, K.W.S.; methodology, K.W.S.; formal analysis, K.W.S.; writing—original draft preparation, K.W.S.; writing—review and editing, H.S., S.F. and V.C.; supervision, H.S., S.F. and V.C.; project administration, H.S., S.F. and V.C.; funding acquisition, H.S., S.F. and V.C. All authors have read and agreed to the published version of the manuscript.

**Funding:** This research was funded by the Water Research Commission of South Africa (WRC), grant number C2020/2021-00646 and the National Research Foundation of South Africa (NRF).

**Acknowledgments:** The authors acknowledge the African Centre for Crop Improvement (ACCI) for technical assistance and overall research support.

**Conflicts of Interest:** The authors declare no conflict of interest.

## References

1. Corbeels, M.; Cardinael, R.; Naudin, K.; Guibert, H.; Torquebiau, E. The 4 per 1000 goal and soil carbon storage under agroforestry and conservation agriculture systems in sub-Saharan Africa. *Soil Tillage Res.* **2019**, *188*, 16–26. [CrossRef]

2. Lal, R. Residue management, conservation tillage and soil restoration for mitigating greenhouse effect by CO₂-enrichment. *Soil Tillage Res.* **1997**, *43*, 81–107. [CrossRef]

3. Chenu, C.; Angers, D.A.; Barre, P.; Derrien, D.; Arrouays, D.; Balesdent, J. Increasing organic stocks in agricultural soils: Knowledge gaps and potential innovations. *Soil Tillage Res.* **2019**, *188*, 41–52. [CrossRef]

4. Rasse, D.P.; Rumpel, C.; Dignac, M.-F. Is soil carbon mostly root carbon? Mechanisms for a specific stabilisation. *Plant Soil* **2005**, *269*, 341–356. [CrossRef]

5. Kätterer, T.; Bolinder, M.A.; Andrén, O.; Kirchmann, H.; Menichetti, L. Roots contribute more to refractory soil organic matter than above-ground crop residues, as revealed by a long-term field experiment. *Agric. Ecosyst. Environ.* **2011**, *141*, 184–192. [CrossRef]

6. Hirte, J.; Leifeld, J.; Abiven, S.; Oberholzer, H.-R.; Mayer, J. Below ground carbon inputs to soil via root biomass and rhizodeposition of field-grown maize and wheat at harvest are independent of net primary productivity. *Agric. Ecosyst. Environ.* **2018**, *265*, 556–566. [CrossRef]

7. Kell, D.B. Breeding crop plants with deep roots: Their role in sustainable carbon, nutrient and water sequestration. *Ann. Bot.* **2011**, *108*, 407–418. [CrossRef]

8. Dolezal, J.; Jandova, V.; Macek, M.; Liancourt, P. Contrasting biomass allocation responses across ontogeny and stress gradients reveal plant adaptations to drought and cold. *Funct. Ecol.* **2020**, *35*, 32–42. [CrossRef]

9. McCarthy, M.C.; Enquist, B. Consistency between an allometric approach and optimal partitioning theory in global patterns of plant biomass allocation. *Funct. Ecol.* **2007**, *21*, 713–720. [CrossRef]

10. Eziz, A.; Yan, Z.; Tian, D.; Han, W.; Tang, Z.; Fang, J. Drought effect on plant biomass allocation: A meta-analysis. *Ecol. Evol.* **2017**, *7*, 11002–11010. [CrossRef]

11. Lin, Y.; Yi, X.; Tang, S.; Chen, W.; Wu, F.; Yang, X.; Jiang, X.; Shi, H.; Ma, J.; Chen, G.; et al. Dissection of Phenotypic and Genetic Variation of Drought-Related Traits in Diverse Chinese Wheat Landraces. *Plant Genome* **2019**, *12*, 190025. [CrossRef] [PubMed]

12. Verma, H.; Borah, J.L.; Sarma, R. Variability Assessment for Root and Drought Tolerance Traits and Genetic Diversity Analysis of Rice Germplasm using SSR Markers. *Sci. Rep.* **2019**, *9*, 16513. [CrossRef] [PubMed]

13. Palta, J.A.; Chen, X.; Milroy, S.P.; Rebetzke, G.J.; Dreccer, M.F.; Watt, M. Large root systems: Are they useful in adapting wheat to dry environments? *Funct. Plant Biol.* **2011**, *38*, 347–354. [CrossRef] [PubMed]

14. Figueroa-Bustos, V.; Palta, J.A.; Chen, Y.; Siddique, K.H. Characterization of Root and Shoot Traits in Wheat Cultivars with Putative Differences in Root System Size. *Agronomy* **2018**, *8*, 109. [CrossRef]

15. Fang, Y.; Xu, B.; Turner, N.C.; Li, F. Does root pruning increase yield and water-use efficiency of winter wheat? *Crop Pasture Sci.* **2010**, *61*, 899–910. [CrossRef]

16. Schneider, H.M.; Lynch, J.P. Should Root Plasticity Be a Crop Breeding Target? *Front. Plant Sci.* **2020**, *11*, 546. [CrossRef]

17. Lynch, J.P. Root phenes that reduce the metabolic costs of soil exploration: Opportunities for 21st century agriculture. *Plant Cell Environ.* **2015**, *38*, 1775–1784. [CrossRef]

18. Hirte, J.; Leifeld, J.; Abiven, S.; Mayer, J. Maize and wheat root biomass, vertical distribution, and size class as affected by fertilization intensity in two long-term field trials. *Field Crop. Res.* **2018**, *216*, 197–208. [CrossRef]

19. Friedli, C.N.; Abiven, S.; Fossati, D.; Hund, A. Modern wheat semi-dwarfs root deep on demand: Response of rooting depth to drought in a set of Swiss era wheats covering 100 years of breeding. *Euphytica* **2019**, *215*, 85. [CrossRef]

20. Manschadi, A.M.; Christopher, J.; Devoil, P.; Hammer, G. The role of root architectural traits in adaptation of wheat to water-limited environments. *Funct. Plant Biol.* **2006**, *33*, 823–837. [CrossRef]

21. Xie, Q.; Fernando, K.M.C.; Mayes, S.; Sparkes, D.L. Identifying seedling root architectural traits associated with yield and yield components in wheat. *Ann. Bot.* **2017**, *119*, 1115–1129. [CrossRef] [PubMed]

22. Bai, C.; Ge, Y.; Ashton, R.W.; Evans, J.; Milne, A.; Hawkesford, M.J.; Whalley, W.R.; Parry, M.A.J.; Melichar, J.; Feuerhelm, D.; et al. The relationships between seedling root screens, root growth in the field and grain yield for wheat. *Plant Soil* **2019**, *440*, 311–326. [CrossRef]
23. Lynch, J.P. Roots of the Second Green Revolution. *Aust. J. Bot.* **2007**, *55*, 493–512. [CrossRef]

24. Dwivedi, S.L.; Britt, A.B.; Tripathi, L.; Sharma, S.; Upadhyaya, H.D.; Ortiz, R. Haploids: Constraints and opportunities in plant breeding. *Biotechnol. Adv.* **2015**, *33*, 812–829. [CrossRef]

25. Yan, G.; Liu, H.; Wang, H.; Lu, Z.; Wang, Y.; Mullan, D.; Hamblin, J.; Liu, C. Accelerated Generation of Selfed Pure Line Plants for Gene Identification and Crop Breeding. *Front. Plant Sci.* **2017**, *8*, 1786. [CrossRef] [PubMed]

26. Reynolds, M.P.; Pellegrineschi, A.; Skovmand, B. Sink-limitation to yield and biomass: A summary of some investigations in spring wheat. *Ann. Appl. Biol.* **2005**, *146*, 39–49. [CrossRef]

27. Shamuyarira, K.W.; Shimelis, H.; Mathew, I.; Zengeni, R.; Chaplot, V. A meta-analysis of combining ability effects in wheat for agronomic traits and drought adaptation: Implications for optimizing biomass allocation. *Crop Sci.* **2022**, *62*, 139–156. [CrossRef]

28. Chen, X.-X.; Zhang, W.; Liang, X.-Y.; Liu, Y.-M.; Xu, S.-J.; Zhao, Q.-Y.; Du, Y.-F.; Zhang, L.; Chen, X.-P.; Zou, C.-Q. Physiological and developmental traits associated with the grain yield of winter wheat as affected by phosphorus fertilizer management. *Sci. Rep.* **2019**, *9*, 16580. [CrossRef]

29. Bastos, L.M.; Carciochi, W.; Lollato, R.P.; Jaenisch, B.R.; Rezende, C.R.; Schwalbert, R.; Prasad, P.V.; Zhang, G.; Fritz, A.K.; Foster, C.; et al. Winter Wheat Yield Response to Plant Density as a Function of Yield Environment and Tillering Potential: A Review and Field Studies. *Front. Plant Sci.* **2020**, *11*, 54. [CrossRef]

30. Figueroa-Bustos, V.; Palta, J.A.; Chen, Y.; Siddique, K.H. Early Season Drought Largely Reduces Grain Yield in Wheat Cultivars with Smaller Root Systems. *Plants* **2019**, *8*, 305. [CrossRef]

31. Středa, T.; Dostál, V.; Horáková, V.; Chloupek, O. Effective use of water by wheat varieties with different root system sizes in rain-fed experiments in Central Europe. *Agric. Water Manag.* **2012**, *104*, 203–209. [CrossRef]

32. Fang, Y.; Du, Y.; Wang, J.; Wu, A.; Qiao, S.; Xu, B.; Zhang, S.; Siddique, K.; Chen, Y. Moderate Drought Stress Affected Root Growth and Grain Yield in Old, Modern and Newly Released Cultivars of Winter Wheat. *Front. Plant Sci.* **2017**, *8*, 672. [CrossRef] [PubMed]

33. Severini, A.D.; Wasson, A.P.; Evans, J.R.; Richards, R.A.; Watt, M. Root phenotypes at maturity in diverse wheat and triticale genotypes grown in three field experiments: Relationships to shoot selection, biomass, grain yield, flowering time, and environment. *Field Crop. Res.* **2020**, *255*, 107807. [CrossRef]

34. Qi, Y.; Wei, W.; Chen, C.; Chen, L. Plant root-shoot biomass allocation over diverse biomes: A global synthesis. *Glob. Ecol. Conserv.* **2019**, *18*, e00606. [CrossRef]

35. Ma, S.-C.; Li, F.-M.; Xu, B.; Huang, Z.-B. Effect of lowering the root/shoot ratio by pruning roots on water use efficiency and grain yield of winter wheat. *Field Crop. Res.* **2010**, *115*, 158–164. [CrossRef]

36. Chang, H.; Chen, D.; Kam, J.; Richardson, T.; Drenth, J.; Guo, X.; McIntyre, C.L.; Chai, S.; Rae, A.L.; Xue, G.-P. Abiotic stress upregulated TaZFP34 represses the expression of type-B response regulator and SHY2 genes and enhances root to shoot ratio in wheat. *Plant Sci.* **2016**, *252*, 88–102. [CrossRef]

37. Bacher, H.; Sharaby, Y.; Walia, H.; Peleg, Z. Modifying root-to-shoot ratio improves root water influxes in wheat under drought stress. *J. Exp. Bot.* **2021**, *73*, 1643–1654. [CrossRef]

38. Lilley, J.; Kirkegaard, J. Benefits of increased soil exploration by wheat roots. *Field Crop. Res.* **2011**, *122*, 118–130. [CrossRef]

39. Feng, F.; Han, Y.; Wang, S.; Yin, S.; Peng, Z.; Zhou, M.; Gao, W.; Wen, X.; Qin, X.; Siddique, K. The Effect of Grain Position on Genetic Improvement of Grain Number and Thousand Grain Weight in Winter Wheat in North China. *Front. Plant Sci.* **2018**, *9*, 129. [CrossRef]

40. Fischer, R. The importance of grain or kernel number in wheat: A reply to Sinclair and Jamieson. *Field Crop. Res.* **2008**, *105*, 15–21. [CrossRef]

41. Shokat, S.; Sehgal, D.; Vikram, P.; Liu, F.; Singh, S. Molecular Markers Associated with Agro-Physiological Traits under Terminal Drought Conditions in Bread Wheat. *Int. J. Mol. Sci.* **2020**, *21*, 3156. [CrossRef] [PubMed]

42. Arjona, J.; Royo, C.; Dreisigacker, S.; Ammar, K.; Villegas, D. Effect of Ppd-A1 and Ppd-B1 Allelic Variants on Grain Number and Thousand Kernel Weight of Durum Wheat and Their Impact on Final Grain Yield. *Front. Plant Sci.* **2018**, *9*, 888. [CrossRef]

43. Stella, T.; Mouratiadou, I.; Gaiser, T.; Berg-Mohnicke, M.; Wallor, E.; Ewert, F.; Nendel, C. Estimating the contribution of crop residues to soil organic carbon conservation. *Environ. Res. Lett.* **2019**, *14*, 094008. [CrossRef]

44. Zadoks, J.C.; Chang, T.T.; Konzak, C.F. A decimal code for the growth stages of cereals. *Weed Res.* **1974**, *14*, 415–421. [CrossRef]

45. Department of Agriculture, Forestry and Fisheries (DAFF). Wheat Production Guideline; Department of Agriculture, Forestry and Fisheries: Pretoria, South Africa, 2016. Available online: https://www.nda.agric.za/Portals/0/Brochures%20and%20Production%20guidelines/Wheat%20-%20Production%20Guideline.pdf (accessed on 19 April 2021).

46. Payne, R.; Murray, D.; Harding, S. *An Introduction to the GenStat Command Language*; VSN International: Hempstead, UK, 2017.

47. Dewey, D.R.; Lu, K.H. A Correlation and Path-Coefficient Analysis of Components of Crested Wheatgrass Seed Production. *Agron. J.* **1959**, *51*, 515–518. [CrossRef]

48. R Development Core Team. R: A Language and Environment for Statistical Computing; R Foundation for Statistical Computing: Vienna, Austria, 2020. Available online: https://www.R-project.org/ (accessed on 10 February 2022).