*Article* **Elevated Atmospheric CO<sup>2</sup> Concentration Influences the Rooting Habits of Winter-Wheat (***Triticum aestivum* **L.) Varieties**

**Balázs Varga \*, Zsuzsanna Farkas, Emese Varga-László \*, Gyula Vida and Ottó Veisz**

Agricultural Institute, Centre for Agricultural Research, Eötvös Loránd Research Network, 2462 Martonvásár, Hungary; farkas.zsuzsanna@atk.hu (Z.F.); vida.gyula@atk.hu (G.V.); veisz.otto@atk.hu (O.V.) **\*** Correspondence: varga.balazs@atk.hu (B.V.); emese.laszlo82@gmail.com (E.V.-L.)

**Abstract:** The intensity and the frequency of extreme drought are increasing worldwide. An elevated atmospheric CO<sup>2</sup> concentration could counterbalance the negative impacts of water shortage; however, wheat genotypes show high variability in terms of CO<sup>2</sup> reactions. The development of the root system is a key parameter of abiotic stress resistance. In our study, biomass and grain production, as well as the root growth of three winter-wheat varieties were examined under optimum watering and simulated drought stress in a combination with ambient and elevated atmospheric CO<sup>2</sup> concentrations. The root growth was monitored by a CI-600 in situ root imager and the photos were analyzed by RootSnap software. As a result of the water shortage, the yield-related parameters decreased, but the most substantial yield reduction was first detected in Mv Karizma. The water shortage influenced the depth of the intensive root development, while under water-limited conditions, the root formation occurred in the deeper soil layers. The most intensive root development was observed until the heading, and the maximum root length was recorded at the beginning of the heading. The period of root development took longer under elevated CO<sup>2</sup> concentration. The elevated CO<sup>2</sup> concentration induced an accelerated root development in almost every soil layer, but generally, the CO<sup>2</sup> fertilization induced in the root length of all genotypes and under each treatment.

**Keywords:** cereals; climate change; water shortage; carbon dioxide; root development

#### **1. Introduction**

One of the most important drivers of climatic changes is the increasing atmospheric CO<sup>2</sup> concentration [1]. The rising temperature is a well-known worldwide phenomenon; however, the changes in precipitation vary greatly by region [2–4]. Still, it can be stated that the frequency and the intensity of drought have been increasing in the most important agricultural areas of the world, and this has resulted in reduced harvested yield and has caused uncertainties in this strategic sector [5]. Many studies have focused on the negative effects of a water-limited environment regarding winter-wheat productivity [2,5,6]. Apart from yield reduction, due to intensive drought, water sources can be utilized less efficiently [7,8]. Several recently published studies confirmed that the elevated CO<sup>2</sup> could more or less counterbalance the negative effects of water shortage through the intensification of photosynthesis in C4 plants [9,10] and the regulation of the stomata closure in C3 plants [11,12]. The majority of these studies focused on the aboveground biomass, especially on the harvested yield; however, it must be highlighted that the root system plays an important role in the water and nutrient uptake of the plants [13–15]. Rooting depth, as well as the structure of the root system and the depth of intensive root development, are the crucial factors that can influence plants' water uptake. Still, a well-developed root structure is not the sole determinant of drought tolerance. Roots are often more varied than shoots and are affected by changes in the climate, soil conditions, plant varieties and soil nutrient and water availability throughout the growing season [16]. The accurate examination of the belowground parts of a plant is more complicated than the analysis of the aboveground

**Citation:** Varga, B.; Farkas, Z.; Varga-László, E.; Vida, G.; Veisz, O. Elevated Atmospheric CO<sup>2</sup> Concentration Influences the Rooting Habits of Winter-Wheat (*Triticum aestivum* L.) Varieties. *Sustainability* **2022**, *14*, 3304. https://doi.org/ 10.3390/su14063304

Academic Editor: Roberto Mancinelli

Received: 24 February 2022 Accepted: 10 March 2022 Published: 11 March 2022

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

**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/).

biomass. Conventionally, destructive and non-repeatable approaches have been applied in order to examine the rooting habits of plants, such as soil coring [17,18], the use of mini rhizotrons [19] or pots of different sizes, but all of these methods have their limitations. First of all, the root development and turnover cannot be measured throughout the vegetation period within the same plant stand. Recently, various methods have been developed to detect the properties of the belowground biomass. They mainly consist of imaging and an image-processing phase. Another option could be to use impedance spectroscopy; however, only a few parameters of the root system can be determined by this approach. This tool is not suitable for measuring the parameters of the individual roots or for inspecting the depth of the intensive root development during the plant-growth period [20,21]. MRI (Magnetic Resonance Imaging) or CT (X-ray Computed Tomography) technologies are highly sophisticated methodologies, but they are exceedingly expensive and could be applied only in pot experiments [22–24]. However, under real field conditions, as well as in a model experiment, root-scanning technology can be efficiently used throughout vegetation. The transparent polycarbonate tubes can be dug or drilled into the soil at different positions and lengths to detect the rooting habits of the plants from germination until harvest. The advantages of this approach are that the measurements can be carried out in the same position and the potential influences of the various environmental factors on the rooting can be excluded [25]. Wheat landraces are better adapted to changing climatic conditions and to stress environments than modern cultivars due to their population's genetic structure, buffering capacity, and morpho-physiological traits, such as rooting habits conferring adaptability to stress environments [26]. However, even among the recently bred cultivars, there is high variability in terms of abiotic stress tolerance [27,28]. Drought resistance is a complex phenomenon, and it develops through the interaction of various plant properties that are determined by several genes, including dwarfing genes, such as *Rht1, Rht2, Rht8*, etc. [29]. Using dwarfing genes to reduce plant height increases the harvest index, improves lodging resistance and increases grain yield. Their application has been one of the major strategies in developing modern bread-wheat cultivars [30]. The presence of these genes in modern wheat varieties is necessary because due to the cultivation technology, shorter plants are more resistant to lodging. The presence of the effective dwarfing genes in the genome not only reduces the plant height but also can negatively affect the intensity of the root development [29].

The objectives of the study were (1) to determine the dynamics of root development during the vegetation period of winter-wheat varieties, (2) to examine how plants react through the modification of the root development to the drought stress in the different soil layers, and (3) to quantify the genotypic responses to elevated atmospheric CO<sup>2</sup> concentrations through the root development under optimum and limited water availability. The results of the experiments could contribute to the efficiency of plant-breeding activities for improved drought tolerance. Moreover, from a practical point of view, the experiments emphasize that the rooting habits of the different varieties should also be considered by the farmers.

#### **2. Materials and Methods**

#### *2.1. Experimental Layout*

Three registered Hungarian winter-wheat (*Triticum aestivum* L.) varieties, Mv Pálma, Mv Karéj and Mv Karizma carrying *Rht8, Rht1* and *Rht2* dwarfing genes (30), respectively, were selected for the model experiment, which was carried out in two similar climatecontrolled greenhouse chambers of the Agricultural Institute, Centre for Agricultural Research, Hungary. Plants were vernalized at 4 ◦C for 6 weeks and the germinated seeds were planted into plastic containers (120 cm× 90 cm × 100 cm) filled with ca. 1000 liters of a 3:1:1 (*v*/*v*) mixture of soil, sand and humus. Water-soluble fertilizer (14% N, 7% P2O5, 21% K2O, 1% Mg, 1% B, Cu, Mn, Fe, Zn; Volldünger Classic; Kwizda Agro Ltd., Vienna, Austria) was added bi-weekly to both water treatments based on the manufacturer's recommendations. Plant density was 450 plants/m<sup>2</sup> , which is similar to the commonly

applied sowing rate in local agricultural practice (Figures 1b and 2). Each container consisted of 8 rows and each row of 28 plants. ppm and the gas concentration was enhanced to 750 ppm in the other chamber by using a network of perforated pipes placed at a height of 0.5 m above the plants. The uniform distribution was achieved through ventilation. atmospheric CO2 concentration in the control (ambient) chamber was maintained at ~400 ppm and the gas concentration was enhanced to 750 ppm in the other chamber by using a network of perforated pipes placed at a height of 0.5 m above the plants. The uniform

a 3:1:1 (*v/v*) mixture of soil, sand and humus. Water‐soluble fertilizer (14% N, 7% P2O5, 21% K2O, 1% Mg, 1% B, Cu, Mn, Fe, Zn; Volldünger Classic; Kwizda Agro Ltd., Vienna, Austria) was added bi‐weekly to both water treatments based on the manufacturer's rec‐ ommendations. Plant density was 450 plants/m2, which is similar to the commonly ap‐ plied sowing rate in local agricultural practice (Figures 1b and 2). Each container consisted

a 3:1:1 (*v/v*) mixture of soil, sand and humus. Water‐soluble fertilizer (14% N, 7% P2O5, 21% K2O, 1% Mg, 1% B, Cu, Mn, Fe, Zn; Volldünger Classic; Kwizda Agro Ltd., Vienna, Austria) was added bi‐weekly to both water treatments based on the manufacturer's rec‐ ommendations. Plant density was 450 plants/m2, which is similar to the commonly ap‐ plied sowing rate in local agricultural practice (Figures 1b and 2). Each container consisted

The air temperature and the additional light intensity of the greenhouse chambers were automatically regulated. The air temperature was increased from the initial 10–12 °C to 24–26 °C over 16 weeks, while air humidity was maintained between 60% and 80% and was regulated by ventilating the greenhouse chambers' air [31]. Whenever it was neces‐ sary, the natural‐light intensity was enhanced by artificial illumination to 500 μmol m–2s–1 at the beginning of the vegetation period, which was gradually increased to 700 μmol m– 2s–1. Thiovit Jet fungicide (Syngenta AG, Basel, Switzerland) (active ingredient: sulfur) and Karate 2.5 WG insecticide (Syngenta AG, Basel, Switzerland) (active ingredient: lambda‐ cyhalothrin) were applied two times (BBCH 23 and 37) [32] following the distributor's recommendations by the dosage against powdery mildew and aphids, respectively. The atmospheric CO2 concentration in the control (ambient) chamber was maintained at ~400

The air temperature and the additional light intensity of the greenhouse chambers were automatically regulated. The air temperature was increased from the initial 10–12 °C to 24–26 °C over 16 weeks, while air humidity was maintained between 60% and 80% and was regulated by ventilating the greenhouse chambers' air [31]. Whenever it was neces‐ sary, the natural‐light intensity was enhanced by artificial illumination to 500 μmol m–2s–1 at the beginning of the vegetation period, which was gradually increased to 700 μmol m– 2s–1. Thiovit Jet fungicide (Syngenta AG, Basel, Switzerland) (active ingredient: sulfur) and Karate 2.5 WG insecticide (Syngenta AG, Basel, Switzerland) (active ingredient: lambda‐ cyhalothrin) were applied two times (BBCH 23 and 37) [32] following the distributor's recommendations by the dosage against powdery mildew and aphids, respectively. The

*Sustainability* **2022**, *14*, x FOR PEER REVIEW 3 of 14

*Sustainability* **2022**, *14*, x FOR PEER REVIEW 3 of 14

of 8 rows and each row of 28 plants.

of 8 rows and each row of 28 plants.

distribution was achieved through ventilation.

**Figure 1.** (**a**) The position of the polycarbonate tubes in the container; (**b**) The plant stand after the planting. **Figure 1.** (**a**) The position of the polycarbonate tubes in the container; (**b**) The plant stand after the planting. **Figure 1.** (**a**) The position of the polycarbonate tubes in the container; (**b**) The plant stand after the planting.

the polycarbonate tubes for scanning the root system. **Figure 2.** The experimentaldesign inthe ambient chamber. The tubes with the black end caps are the polycarbonate tubes for scanning the root system. **Figure 2.** The experimental design in the ambient chamber. The tubes with the black end caps are the polycarbonate tubes for scanning the root system.

The air temperature and the additional light intensity of the greenhouse chambers were automatically regulated. The air temperature was increased from the initial 10–12 ◦C to 24–26 ◦C over 16 weeks, while air humidity was maintained between 60% and 80% and was regulated by ventilating the greenhouse chambers' air [31]. Whenever it was necessary, the natural-light intensity was enhanced by artificial illumination to 500 µmol m−<sup>2</sup> s <sup>−</sup><sup>1</sup> at the beginning of the vegetation period, which was gradually increased to 700 µmol m−<sup>2</sup> s −1 . Thiovit Jet fungicide (Syngenta AG, Basel, Switzerland) (active ingredient: sulfur) and Karate 2.5 WG insecticide (Syngenta AG, Basel, Switzerland) (active ingredient: lambdacyhalothrin) were applied two times (BBCH 23 and 37) [32] following the distributor's recommendations by the dosage against powdery mildew and aphids, respectively. The atmospheric CO<sup>2</sup> concentration in the control (ambient) chamber was maintained at ~400 ppm and the gas concentration was enhanced to 750 ppm in the other chamber by using a network of perforated pipes placed at a height of 0.5 m above the plants. The uniform distribution was achieved through ventilation.

The containers were separated into two parts of equal size (optimum watered and drought-stressed) by using water-insulating PVC foil (thickness 1 mm).

The water-holding capacity of the soil was determined by using the gravimetric method before starting the experiment and the control treatments were watered until optimum (60%) soil-water content [SWC]. The water content of the soil was monitored by 5 TA sensors (Decagon Devices Ltd., Pullman, WA, USA) at 3 depths (30, 60 and 90 cm). Water-stressed plants did not receive additional watering after the planting until the volumetric soil-water content dropped below 8–10% (average of the three depths). The plants in the drought-stress treatment were irrigated first at the BBCH 51 stage, 68 days after the planting. Afterwards, halved water doses compared to the control were applied to the stress-treated stands. TA sensors (Decagon Devices Ltd., Pullman, WA, USA) at 3 depths (30, 60 and 90 cm). Water‐stressed plants did not receive additional watering after the planting until the vol‐ umetric soil‐water content dropped below 8–10% (average of the three depths). The plants in the drought‐stress treatment were irrigated first at the BBCH 51 stage, 68 days after the planting. Afterwards, halved water doses compared to the control were applied to the stress‐treated stands.

The containers were separated into two parts of equal size (optimum watered and

The water‐holding capacity of the soil was determined by using the gravimetric method before starting the experiment and the control treatments were watered until op‐ timum (60%) soil‐water content [SWC]. The water content of the soil was monitored by 5

#### *2.2. Measurements*

Without overlapping, three transparent polycarbonate tubes were set up in the containers in a horizontal position, at 30, 60 and 90 cm soil depths (Figure 1a). The root development/turnover was monitored every two weeks in the same positions of the tubes by a CI-600 in situ root imager (CID-Bioscience Ltd., Camas, WA, USA). The plant phenophases were ranked according to the BBCH scale (Table 1) [32]. *2.2. Measurements* Without overlapping, three transparent polycarbonate tubes were set up in the con‐ tainers in a horizontal position, at 30, 60 and 90 cm soil depths (Figure 1a). The root de‐ velopment/turnover was monitored every two weeks in the same positions of the tubes by a CI‐600 in situ root imager (CID‐Bioscience Ltd., Camas, WA, USA). The plant phe‐ nophases were ranked according to the BBCH scale (Table 1) [32].


**Table 1.** List of phenophases when the root length measurements were carried. **Table 1.** List of phenophases when the root length measurements were carried

BBCH 83 Ripening, early dough

*Sustainability* **2022**, *14*, x FOR PEER REVIEW 4 of 14

drought‐stressed) by using water‐insulating PVC foil (thickness 1 mm).

RootSnap software (CID-Bioscience Ltd., Camas, WA, USA) was applied for image processing and to determine the root length (Figure 3). Simultaneously, soil temperature and soil-water content were monitored continuously in the three layers by 5 TE sensors and EM50 data loggers (Decagon Devices Ltd., Pullman, WA, USA). RootSnap software (CID‐Bioscience Ltd., Camas, WA, USA) was applied for image processing and to determine the root length (Figure 3). Simultaneously, soil temperature and soil‐water content were monitored continuously in the three layers by 5 TE sensors and EM50 data loggers (Decagon Devices Ltd., Pullman, WA, USA).

**Figure 3.** An image taken by the CI-600 root scanner under analysis by the RootSnap software.

At the end of the vegetation period, all plants were manually harvested row by row (8 replications). The weight of the total aboveground biomass, grain weight and

thousand-kernel weight was measured by a digital balance (ME 1002E, Mettler-Toledo Ltd., Worthington, OH, USA), and the harvest index (Equation (1)) was calculated.

$$\text{Harvest index} = \frac{\text{grain yield (g)}}{\text{total aboveground biomass} \, (\text{g})} \ast 100 \tag{1}$$

Relative changes in the root length in response to elevated carbon dioxide levels were calculated as:

$$\text{CO}\_2\text{ responses} = \frac{\text{Ex}}{\text{A}}\tag{2}$$

where A is the root length at the 400 ppm CO<sup>2</sup> level, Ex is the root length at the 700 ppm CO<sup>2</sup> level.

#### *2.3. Analysis*

The experimental design (2 × 2 × 3) consisted of two CO<sup>2</sup> levels, two watering levels and three genotypes. The rooting habits of a closed plant stand were analyzed. The scanning was always carried out in three replications in the same position. The effects of the tested factors on yield parameters were determined by multi-way ANOVA and means were compared by Tukey's HSD test (*p* ≤ 0.05). The significant differences in the length parameters between the water treatments and CO<sup>2</sup> levels were evaluated by Student's t-test and the significant differences between the root lengths measured in the plant growth stages were analyzed by one-way ANOVA followed by Tukey's HSD test.

#### **3. Results**

*3.1. Effects of Water Shortage and Elevated Atmospheric CO<sup>2</sup> Concentration on Yield Parameters*

The biomass (BM), grain yield (GY), thousand-kernel weight (TKW) and harvest index (HI) of the individual treatments are summarized in Table 2.


**Table 2.** Responses of winter-wheat varieties to elevated CO<sup>2</sup> and water shortage.

BM: aboveground biomass; GY: grain yield; TKW thousand-kernel weight; HI: harvest index; NC: ambient CO2; EC: 750 ppm CO2; C: controlled watering; D: drought stress, n.s. the effects of the factor are not significant (the presented data are means of eight replications).

The ANOVA shows that the effects of the watering levels, genotypes and the CO<sup>2</sup> concentrations significantly influenced the grain yield, but their interactions were not statistically significant. The effects of the watering and the genotypes were significant on the aboveground biomass, but the ANOVA shows that the interactions of the factors remained insignificant (Tables S1 and S2). Different CO<sup>2</sup> levels and irrigation regimes were assumed to have contradictory effects on the yield-related parameters. Water shortage resulted in decreased BM and GY both under ambient and elevated CO<sup>2</sup> concentrations. The ratio of the yield reduction did not differ between the two CO<sup>2</sup> concentrations. Compared to the well-watered plant stands, the biomass of Mv Pálma and Mv Karéj decreased by 8.02% and 10.03% under ambient CO<sup>2</sup> concentration and 11.2% and 10.45% under elevated CO2, respectively, in the drought-stress treatment (Table 2).

Generally, the simulated drought stress did not significantly affect the TKW and the HI, but the effects of the genotypes and watering levels were statistically significant in terms of both parameters (Tables S3 and S4). In the case of Mv Pálma, water shortage led to a 12.75% and 13.75% increase in TKW compared with the non-stressed controls under 400 ppm and 750 ppm CO<sup>2</sup> concentration, respectively. The highest yielding capacity was observed for Mv Karizma, and this genotype showed the lowest yield reduction under the simulated drought stress. In Mv Karizma, the elevated CO<sup>2</sup> concentration played a role in counterbalancing the negative effects of the water shortage. Taking the average of the three examined varieties, the CO<sup>2</sup> enrichment increased the biomass (by 4.57% and 4.97% under normal watering and drought stress, respectively) as well as the grain yield (by 8.1% and 8.96% under normal watering and drought stress, respectively) (Table 2). Improved harvest indices (4.3% and 4.0% under normal watering and drought-stress conditions, respectively) could be determined as a result of the observed tendencies in biomass and grain weight, which would be highly favorable for nutrient- and water-utilization efficiency. In Mv Karéj, no differences were observed between the two water treatments in terms of CO<sup>2</sup> reactions. The interaction of water shortage and CO<sup>2</sup> fertilization showed opposite tendencies in the other two varieties. The increase in the yield parameters was more intense under the control water supply in Mv Pálma, while the positive effect of the elevated CO<sup>2</sup> compared with the control was more significant under drought-stress conditions in Mv Karizma (Table 2).

### *3.2. Dynamics of the Root Development of Winter-Wheat Varieties under Optimum Watering and Drought-Stressed Conditions Grown at Ambient and Elevated CO<sup>2</sup> Concentrations*

Under ambient CO<sup>2</sup> concentration, the root length of Mv Pálma continuously developed in the well-watered treatment after planting until it reached its maximum level at the BBCH 51 stage in the two upper soil layers (30 and 60 cm), while a faster root development was detected in the drought-stressed treatment when the root length did not significantly increase after the BBCH 29 and BBCH 37 stages at 30 and 60 cm, respectively. (Table 3). Significantly higher root length was observed at 30 cm between the BBCH 29 and 51 stages under drought-stressed conditions, but the root length did not significantly differ between the water treatment at 60 cm except at the BBCH 21 and 83 stages. At 90 cm, the root length of Mv Pálma reached its maximum (does not significantly increase further) at BBCH 69 under optimum irrigation and at BBCH 51 under drought-stressed conditions, and significantly higher root length was measured under drought stress between BBCH 37 and BBCH 83 than in the control treatment (Table 2). Under elevated CO2, the root length of Mv Pálma was consequently higher in the well-irrigated treatment at 30 cm between BBCH 21 and BBCH 77 than under drought stress (Table 3). The root length did not significantly increase at 30 or 90 cm after the BBCH 37 stage, but at 60 cm, the root length did not significantly increase after the BBCH 29 stage under drought-stressed conditions.

Under ambient CO<sup>2</sup> concentration, the root length of Mv Karéj reached its maximum in the upper soil layer (30 cm) at the BBCH 29 stage in both water treatments, and the root length was significantly higher under optimum watering at BBCH 51 and BBCH 77 than in the drought-stressed treatment (Table 4). In the middle layer (60 cm), the root length did not significantly increase after BBCH 37 and BBCH 29 in the control and drought-stressed treatments, respectively, and a significant difference between the water treatments can be observed only at the BBCH 29 phenophase. In Mv Karéj, the maximum root length was observed at the BBCH 51 stage at 90 cm in both water treatments, but significantly higher values were measured under drought-stressed conditions from BBCH 37 compared to the control treatment.

**Table 3.** Root length of Mv Pálma in different stages of vegetation under well-watered and droughtstressed conditions in different soil layers by ambient and elevated atmospheric CO<sup>2</sup> concentrations.


C30, D30, C60, D60, C90 and D90 indicate the soil layer at 30 cm under controlled watering, the soil layer at 30 cm under drought stress, the soil layer at 60 cm under controlled watering, the soil layer at 60 cm under drought stress, the soil layer at 90 cm under controlled watering and the soil layer at 90 cm under drought stress, respectively. Lowercase letters indicate significant differences between drought and control treatments within the same soil layer and CO<sup>2</sup> treatment (Student's *t*-test) (*p* < 0.05) and the numbers in superscripts indicate significant differences between phenophases (HSD test (*p* ≤ 0.05) (n = 3), n.r., no roots were observed.

**Table 4.** Root length of Mv Karéj at different stages of vegetation under well-watered and droughtstressed conditions in different soil layers by ambient and elevated atmospheric CO<sup>2</sup> concentrations.


C30, D30, C60, D60, C90 and D90 indicate the soil layer at 30 cm under controlled watering, the soil layer at 30 cm under drought stress, the soil layer at 60 cm under controlled watering, the soil layer at 60 cm under drought stress, the soil layer at 90 cm under controlled watering and the soil layer at 90 cm under drought stress, respectively. Lowercase letters indicate significant differences between drought and control treatments within the same soil layer and CO<sup>2</sup> treatment (Student's *t*-test) (*p* < 0.05) and the numbers in superscripts indicate significant differences between phenophases (HSD test (*p* ≤ 0.05) (n = 3), n.r., no roots were observed.

Under 750 ppm CO2, the root length of Mv Karéj reached its maximum at the BBCH 37 stage at each soil layer under optimum and drought-stressed conditions. No significant differences were observed between the water treatments at the 30 cm soil layer after the BBCH 21 stage, but the CO<sup>2</sup> enrichment resulted in a significant decrease in root length in each phenophase at 60 and 90 cm (Table 4).

The root development of Mv Karizma showed the highest variability among the studied varieties under ambient CO<sup>2</sup> concentration in different soil layers and irrigation regimes. Intensive root development was observed from planting until the BBCH 29 stage at 30 cm, but the root length did not increase further during vegetation. Oppositely, a significant reduction in root length was determined from the BBCH 69 stage (Table 5). The root length at 30 cm was significantly higher under optimum watering than under drought stress in each phenophase. At 60 cm, the root-length development was faster under drought-stress conditions than under optimum irrigation, and the measured data did not significantly increase after BBCH 37 and BBCH 51 under drought stress and optimum

watering, respectively. The water shortage induced a significantly more developed root system in Mv Karizma for each phenophase at 60 cm (Table 5). The root system developed intensively until the BBCH 51 stage at 90 cm under both watering regimes, then a root turnover can be observed in the BBCH 69 and 77 stages, but another intensive root formation was detected at the BBCH 83 stage. Significantly higher root-length values were measured in Mv Karizma under drought-stressed conditions at 90 cm than in the control treatment from BBCH 29 until BBCH 83, except BBCH 37 at 400 ppm CO<sup>2</sup> concentration. Under elevated CO<sup>2</sup> levels (750 ppm), a delay can be observed at 30 cm in the time that the root length takes to reach its maximum compared to the ambient treatment. The maximum values were observed at the BBCH 37 stage in both water treatments. The water supply had no significant effects on the root length of Mv Karizma during vegetation at 30 cm under 750 ppm CO<sup>2</sup> between the BBCH29 and BBCH77 stages. At 60 cm, the highest root length was observed at the end of vegetation under optimum watering, while the highest data were measured at the BBCH 37 and BBCH 51 stages under limited water supply. The water shortage induced a significant increase in root length in each phenophase of Mv Karizma under CO<sup>2</sup> enrichment. The trends in root development at 90 cm were similar to that at 60 cm, and the water shortage resulted in an increase in root length between the BBCH 37 and BBCH 77 stages (Table 5).

**Table 5.** Root length of Mv Karizma at different stages of vegetation under well-watered and droughtstressed conditions in different soil layers by ambient and elevated atmospheric CO<sup>2</sup> concentrations.


C30, D30, C60, D60, C90 and D90 indicate the soil layer at 30 cm under controlled watering, the soil layer at 30 cm under drought stress, the soil layer at 60 cm under controlled watering, the soil layer at 60 cm under drought stress, the soil layer at 90 cm under controlled watering and the soil layer at 90 cm under drought stress, respectively. Lowercase letters indicate significant differences between drought and control treatments within the same soil layer and CO<sup>2</sup> treatment (Student's *t*-test) (*p* < 0.05) and the numbers in superscripts indicate significant differences between phenophases (HSD test (*p* ≤ 0.05) (n = 3), n.r., no roots were observed.

#### *3.3. CO<sup>2</sup> Reactions of Winter-Wheat Genotypes during Vegetation at Different Soil Layers under Well-Watered and Drought-Stressed Conditions*

The responses of Mv Pálma to the elevated CO<sup>2</sup> showed a variability during the vegetation period and significant alterations were determined between the well-watered and drought-stressed treatments (Figure 4). The CO<sup>2</sup> response was positive and significant from the sowing until the BBCH 21 stage at 30 and 60 cm, and this tendency can be observed for 90 cm at the BBCH 29 stage. Generally, the CO<sup>2</sup> enrichment induced faster root development in Mv Pálma.

Positive and statistically significant CO<sup>2</sup> reactions were observed between the BBCH 37 and BBCH 83 growth stages under optimum watering at 90 cm, and between BBCH 51 and BBCH 69 at 30 cm (Figure 4). The CO<sup>2</sup> fertilization resulted in a significant decrease in the root length of Mv Pálma under drought-stressed conditions at the BBCH 37 stage, and at the later stages of vegetation in each soil layer except in the BBCH 51 stage at 60 cm (Figure 4).

Overall, the CO<sup>2</sup> responses of Mv Karéj were negative in each phenophase, and the reactions were not influenced by the watering (Figure 5). The negative CO<sup>2</sup> responses were more intensive at the end of vegetation, which indicates that the CO<sup>2</sup> fertilization influenced the water balance of plants and induced faster root turnover. The unfavorable impacts of the CO<sup>2</sup> enrichment were more intensive under drought-stressed conditions on the root length, but this process could be a component of the survival strategy of this genotype (Figure 5). drought‐stressed treatments (Figure 4). The CO2 response was positive and significant from the sowing until the BBCH 21 stage at 30 and 60 cm, and this tendency can be ob‐ served for 90 cm at the BBCH 29 stage. Generally, the CO2 enrichment induced faster root development in Mv Pálma.

layer at 90 cm under drought stress, respectively. Lowercase letters indicate significant differences between drought and control treatments within the same soil layer and CO2 treatment (Student's *t*‐test) (*p* < 0.05) and the numbers in superscripts indicate significant differences between phe‐

*3.3. CO2 Reactions of Winter‐Wheat Genotypes during Vegetation at Different Soil Layers under*

The responses of Mv Pálma to the elevated CO2 showed a variability during the veg‐ etation period and significant alterations were determined between the well‐watered and

*Sustainability* **2022**, *14*, x FOR PEER REVIEW 9 of 14

nophases (HSD test (*p* ≤ 0.05) (n = 3), n.r., no roots were observed.

*Well‐Watered and Drought‐Stressed Conditions*

**Figure 4.** Relative changes in root length of Mv Pálma to elevated CO2 concentration during the vegetation period at 30, 60 and 90 cm soil layer under optimum watering (control) and drought‐ stressed conditions. Full bars represent significant differences compared to the control (400 ppm) *p* ≤ 0.05 level (n = 3). **Figure 4.** Relative changes in root length of Mv Pálma to elevated CO<sup>2</sup> concentration during the vegetation period at 30, 60 and 90 cm soil layer under optimum watering (control) and droughtstressed conditions. Full bars represent significant differences compared to the control (400 ppm) *p* ≤ 0.05 level (n = 3). dentially negative between the BBCH29 and BBCH 83 stages under drought‐stressed con‐ ditions and positive at 60 cm under controlled irrigation. Negative CO2 reactions were determined at 90 cm at both watering levels, but this tendency was significant at first when plants reached the BBCH 51 stage (Figure 6).

**Figure 5.** Relative changes in root length of Mv Karéj to elevated CO2 concentration during the veg‐ etation period at 30, 60 and 90 cm soil layer under optimum watering (control) and drought‐stressed conditions. Full bars represent significant differences compared to the control (400 ppm) *p ≤* 0.05 level (n = 3). **Figure 5.** Relative changes in root length of Mv Karéj to elevated CO<sup>2</sup> concentration during the vegetation period at 30, 60 and 90 cm soil layer under optimum watering (control) and droughtstressed conditions. Full bars represent significant differences compared to the control (400 ppm) *p* ≤ 0.05 level (n = 3).

The CO<sup>2</sup> reactions of Mv Karizma showed variability during vegetation, and this parameter was influenced significantly by the watering (Figure 6). Roots can be observed only at 30 cm at the BBCH 17 stage, and significant positive CO<sup>2</sup> responses were observed in terms of the root development under drought-stressed conditions, while the opposite tendency was observed under optimum irrigation. The negative impacts of the CO<sup>2</sup> enrichment were detected at 30 cm under optimum irrigation in the BBCH 21 stage, but significant positive responses to both watering levels were observed at 60 cm. The reactions of Mv Karizma to CO<sup>2</sup> enrichment was consequently negative at 30 cm, and this trend was not influenced by the intensity of watering, while the CO<sup>2</sup> reactions were tendentially negative between the BBCH29 and BBCH 83 stages under drought-stressed conditions and positive at 60 cm under controlled irrigation. Negative CO<sup>2</sup> reactions were determined at 90 cm at both watering levels, but this tendency was significant at first when plants reached the BBCH 51 stage (Figure 6). **Figure 5.** Relative changes in root length of Mv Karéj to elevated CO2 concentration during the veg‐ etation period at 30, 60 and 90 cm soil layer under optimum watering (control) and drought‐stressed conditions. Full bars represent significant differences compared to the control (400 ppm) *p ≤* 0.05 level (n = 3).

*Sustainability* **2022**, *14*, x FOR PEER REVIEW 10 of 14

when plants reached the BBCH 51 stage (Figure 6).

tendency was observed under optimum irrigation. The negative impacts of the CO2 en‐ richment were detected at 30 cm under optimum irrigation in the BBCH 21 stage, but significant positive responses to both watering levels were observed at 60 cm . The reac‐ tions of Mv Karizma to CO2 enrichment was consequently negative at 30 cm, and this trend was not influenced by the intensity of watering, while the CO2 reactions were ten‐ dentially negative between the BBCH29 and BBCH 83 stages under drought‐stressed con‐ ditions and positive at 60 cm under controlled irrigation. Negative CO2 reactions were determined at 90 cm at both watering levels, but this tendency was significant at first

**Figure 6.** Relative changes in root length of Mv Karizma to elevated CO<sup>2</sup> concentration during the vegetation period at 30, 60 and 90 cm soil layer under optimum watering (control) and droughtstressed conditions. Full bars represent significant differences compared to the control (400 ppm) *p* ≤ 0.05 level (n = 3).

#### **4. Discussion**

CO<sup>2</sup> molecules in the atmosphere are essential substrates of photosynthesis; therefore, in general, increasing the concentration of CO<sup>2</sup> leads to improved assimilation and crop productivity. Previously, the effects of elevated CO<sup>2</sup> have been widely investigated in various field crops, such as wheat [33,34], maize [35,36], rice [37], sorghum [35], etc. Primarily, these studies determined that CO<sup>2</sup> fertilization improved biomass and grain production, which is also confirmed by our experiment. It has already been correspondingly established that there are differences in the CO<sup>2</sup> reactions of cereal varieties, especially under stress conditions [37–39]. In our experiment, Mv Karizma showed the most intense responses to the increased CO<sup>2</sup> level. Previously conducted studies also highlighted that the increasing CO<sup>2</sup> concentration could partly counterbalance the negative impact of abiotic stresses, such as drought or heat [7,40]. In our experiment, the CO<sup>2</sup> response was variety dependent and the most favorable reactions were observed for Mv Karizma. The yield reduction under drought stress dropped by approximately one half under elevated CO<sup>2</sup> compared to ambient conditions. The physiological background of this result might be that the high CO<sup>2</sup>

concentration modifies the intensity of the photosynthesis [39,41], decreases the stomatal conductance and increases the evaporation [42]. Based on this, it can be concluded that although the influence of elevated CO<sup>2</sup> on the aboveground parts of the plants is well-known, the development of the root system should also offset the increased water loss and nutrient uptake. There are differences in the rooting habits of wheat varieties, even under ambient conditions [43,44]; furthermore, the cultivars also give diverse responses to CO<sup>2</sup> [45,46]. Considering the above, it would be highly important to reveal how these effects interact under various irrigation regimes.

The novel aspect of our study is that the effects of the elevated CO<sup>2</sup> concentration were combined with the drought stress, and the differences in the CO<sup>2</sup> responses of the tested varieties were also determined.

Asseng et al. [47] described in a minirhizotron study that the fastest root-growth development was observed between the visible terminal spikelet stage (BBCH 51) and anthesis (BBCH 65), and the maximum root length was measured beyond 90 days after sowing. Based on our results, the most intensive root growth occurred in the vegetative phase, while the date of the highest root length was influenced by the atmospheric CO<sup>2</sup> as well. Under ambient conditions, the maximum root length could be observed at the BBCH 37 stage, but the elevated CO<sup>2</sup> resulted in a constant increase in the root length until the BBCH 51 phase. Uddin et al. [48] described that elevated CO<sup>2</sup> resulted in the intensive rootlength formation of wheat plants, especially in the upper soil layer. In our study, however, significant differences were observed between the examined varieties in terms of CO<sup>2</sup> reactions. Our study partly confirmed this trend: although accelerated root development was detected in each soil layer induced by the CO<sup>2</sup> fertilization, overall, the root length was reduced under elevated CO<sup>2</sup> concentration. Mitchell et al. [49] proposed that the potential reason for this phenomenon could be that the surplus assimilates availability for extra root growth under elevated CO2. Manschadi et al. [50] found in a wheat experiment, in root-observation chambers, that genotypes that were developed for dry conditions had longer root systems in deep soil layers. In our study, the intensive root formation in the deeper soil layers (60 and 90 cm) under drought-stressed conditions indicate the good adaptive capacity of the examined genotypes, especially that of Mv Pálma and Mv Karizma, to a water-limited environment.

In the root characteristics, no differences were detected in near-isogenic lines in the case of the reduced height genes (*Rht1* and *Rht2*) in Western Australia [51], whereas in Argentina, an increase in the total root length and root weight was associated with reduced height [52]. In our experiment, Mv Karizma had the greatest root length, especially under dry conditions, in the deepest soil layer. This might be due to the presence of the lessefficient dwarfing gene (*Rht2*) in its genome. Mv Pálma had the less-developed root system, even under elevated CO<sup>2</sup> concentration, which could be in accordance with the presence of the *Rht8* gene.

#### **5. Conclusions**

As a consequence of the water shortage, the depth of the intensive root development remained in the deeper soil layers and under elevated CO<sup>2</sup> concentration, the distribution of the root system was more homogeneous in the whole soil profile than under ambient conditions. The elevated CO<sup>2</sup> concentration induced an accelerated root formation, but considering the whole vegetation period, the CO<sup>2</sup> fertilization had a reducing effect on the root length.

The most intensive root development was detected in the vegetative stage of the plants. The maximum root length was observed at the beginning of the heading, and this period took longer under elevated CO<sup>2</sup> concentration. After the heading, the development of the new roots became slower and intensive root turnover was observed. Increased root formation was determined at the maturity stage in almost every treatment, which could be in association with the re-growth of the stubble under field conditions.

The water shortage generally stimulated root growth, but the depth of the intensive root development proved to be variety dependent and this phenomenon was more intensive under ambient CO<sup>2</sup> concentration.

**Supplementary Materials:** The following supporting information can be downloaded at: https: //www.mdpi.com/article/10.3390/su14063304/s1, Table S1: The effects of the tested factors and their interactions on the aboveground biomass, Table S2: The effects of the tested factors and their interactions on the grain yield, Table S3: The effects of the tested factors and their interactions on the thousand-kernel weight, Table S4: The effects of the tested factors and their interactions on the harvest index.

**Author Contributions:** Conceptualization, B.V.; methodology, B.V.; validation, G.V., O.V. data curation, B.V., E.V.-L.; writing—original draft preparation, B.V., Z.F.; writing—review and editing G.V., O.V.; visualization, E.V.-L.; funding acquisition, O.V. All authors have read and agreed to the published version of the manuscript.

**Funding:** This research was funded by the Hungarian Government and the European Union, with the co-funding of the European Regional Development Fund within the framework of the Széchenyi 2020 Program's grant number GINOP-2.3.2-15-2016-00029.

**Institutional Review Board Statement:** Not applicable.

**Informed Consent Statement:** Not applicable.

**Data Availability Statement:** The data presented in this study are available on request from the corresponding author.

**Conflicts of Interest:** The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

#### **References**

