Next Article in Journal
Development of Enumerative and Binomial Sequential Sampling Plans for Monitoring Lymantria dispar (L.) (Lepidoptera Erebidae) in Mediterranean Oak Forests
Previous Article in Journal
Bread Wheat Productivity in Response to Humic Acid Supply and Supplementary Irrigation Mode in Three Northwestern Coastal Sites of Egypt
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Trait Selection for Yield Improvement in Foxtail Millet (Setaria italica Beauv.) under Climate Change in the North China Plain

1
Dryland Farming Institute, Hebei Academy of Agricultural and Forestry Science, Hengshui 053000, China
2
Key Laboratory of Crop Drouht Tolerance Research of Heibei Province, Hengshui 053000, China
3
College of Agronomy, Yulin Normal University, Yulin 537000, China
4
Key Laboratory for Conservation and Utilization of Subtropical Bio-Resources, Education Department of Guangxi Zhuang Autonomous Region, Yulin Normal University, Yulin 537000, China
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Submission received: 8 May 2022 / Revised: 13 June 2022 / Accepted: 21 June 2022 / Published: 23 June 2022

Abstract

:
Weather factors and drought could impact the yield of foxtail millet, and varieties with traits that could alleviate the negative effects of deteriorating weather factors in the future should be developed. A total of 25 foxtail millet varieties were evaluated in experiments from 2016 to 2020 under well-watered (WW) and water-stressed (WS) treatments. Future climate change might favor an increased temperature that impedes grain yield, so varieties with characters that are less sensitive to temperature change are preferred. Varieties with a high panicle dry-weight per plant, thousand grain weight, leaf area and water productivity in deep soil layer usually gave better grain production under both water treatments. Under the WW treatment, low grain abortion rate, optimal chlorophyll and canopy temperature and more roots in the upper soil layer could favor a high yield and drought resistance. Under the WS treatment, varieties with a high harvest index, low rate of water loss and more roots in the upper soil layer usually produced a high yield.

1. Introduction

Foxtail millet [Setaria italica (L.) P. Beauv.] is one of the earliest cultivated crops in human civilization and plays an important role in food and nutritional security in current society. Foxtail millet ranks as the second millet in the total world production of millets [1,2] in China, 1.74 MT foxtail millet was produced in 2014 and China was the second producer of millets in Asia [3]. Foxtail millet is adaptable to a wide range of soil and climatic conditions, and it is much more tolerant to biotic and abiotic stresses as compared with other major crops, such as maize, wheat and rice (which requires about 30% less water than maize [4]). Foxtail millet has some morphological and anatomical features, including thick cell walls, a dense root structure that grows deeper into the soil profile, and a small leaf area, all of which give it an excellent drought tolerance [5]. The crop is considered as an ideal model crop for plant functional studies and molecular mechanisms of drought tolerance [6,7].
Global climate change has a growing impact on crop production [8,9]. Studies have shown that the average global temperature will increase by roughly 0.2 per decade over the next 30 years [10]. New genomic models that take the multiple traits and multiple environments into consideration [11], along with trait × environment, trait × genotype, and trait × genotype × environment interactions, offer a huge potential for the exploitation of correlations between different variables and for the differentiation between effects [12]. Climate change can potentially influence the cereal yield directly by heat and drought [13]. Under rainfed conditions, drought is expected to continue as one of the main abiotic factors affecting crop yields [14]. Krishna et al. [15] reported drought from earhead emergence to grain filling of finger millet reduced the grain yield by 36.6%. Although foxtail millet is drought resistant compared to other cereal crops, drought is still the most prevalent abiotic constraint that causes yield reduction [16]. Drought-resistant breeding of foxtail millet is an effective measure to achieve higher yield potential under drought conditions. Therefore, the trait-based breeding approach would help to improve yield enhancement. Identifying traits that are related to yield and drought resilience from released varieties would be an appropriate option. Krishna et al. [15] found that the leaf temperature of finger millet was the most important physiological trait under drought stress, and leaf temperature can be used as one of the important selection criteria for drought response; among the agronomic traits, productive tillers and panicle size contributed towards the grain yield under drought conditions.
Canopy temperature has been used to study the genotypic response to drought in wheat. Rashid et al. [17] found a significant correlation between canopy temperature and yield under drought conditions and drought-stress susceptibility index values, which indicated their potential for screening drought-tolerant genotypes. Canopy temperatures also provide an indication for yield performance during drought and could effectively be used as a technique to assess the genotypic response to drought [17]. Reza [18] reported that wheat genotypes with a low canopy temperature can maintain a high transpiration and photosynthetic rate, as well as produce a high yield under water-stressed conditions. There was a significant positive correlation between chlorophyll content and yield in both well-watered and water-stressed environments [18]. The significant correlation between canopy temperature and chlorophyll content with yield under water-stressed conditions indicated the potential for screening wheat genotypes for drought response.
The drought resistance index is closely associated with the panicle harvest index; therefore, it is a trait that is related to a better grain-setting and grain-filling capacity [19]. Roots have a range of functions, including anchorage, mechanical support, nutrient and water uptake, and signaling. Roots are also extremely sensitive to water deficit and high temperatures [20]. Drought has different effects depending on the severity; if a moderate drought occurs, primary root development is increased, while lateral roots are repressed [21].
Foxtail millet is mainly grown in northern China’s arid and semi-arid region, where the development of drought seriously influences foxtail millet production. Although there have been about 870 varieties released since 1950 [22], the development of drought resistant varieties is hampered by the lack of effective selection criteria. In this study, 25 foxtail millet varieties that were collected from 2016 to 2020 were used to investigate the influence of weather and soil water conditions on the grain yield and the beneficial traits of foxtail millet varieties for future foxtail millet cultivation and breeding.

2. Materials and Methods

2.1. Experiment Site

The experiment was conducted at Hengshui Dryland Agricultural Experimental Station (37°13′ N, 114°37′ E; 23 m above sea level), Hebei province, Northern China from 2016 to 2020. The experiment location was in a semi-arid area of China and the soil was classified as silt loam with a pH of 8.0 to 8.1. The soil organic matter was 1.0 to 1.3% and the soil bulk density was 1.2 to 1.5 g cm−3 for the tillage soil layer.

2.2. Experimental Design and Field Management

During the five growing seasons (2016–2020), a total of 25 recently released varieties (Appendix A, Table A1) collected from different breeders were used. Some of the varieties were replaced by the newly released varieties each season. The inclusion of the newly released varieties each season ensured that their characters were considered. The variety Jigu19 was grown in all five seasons to analyze the weather factors on seasonal yield variation.
The experimental design was a randomized block design with three replicates for each variety. The plot size was 2.0 m ×1.5 m and each plot was surrounded by a 20 cm thick concrete wall to avoid water exchange. The depth of each plot was 3.0 m and contained 2 m of soil layer at the top, and 1m of buffer layer at the bottom. There was a movable shed to prevent rainfall. Foxtail millet plants were grown under two irrigation treatments (ITr) during the growing seasons, i.e., water-stressed (WS)—no irrigation; well-watered (WW)—two irrigations. Irrigation was conducted at the jointing stage and the grain-filling stage. For each irrigation, 900 m3 ha−1 (i.e., 90 mm) of water was applied; the amount of irrigation water was controlled by a water meter. The PVC tube was used to irrigate, and this was moved around the plot to ensure that the water content was uniform.
All field managements were the same for all plots, except for the irrigation. The seeds were sown manually. At the seedling, the plants were thinned to a density of 60 plants/m2. The sowing date and harvest date are shown in Table A2. The row spacing was set at 20 cm. Before planting, the soil was irrigated to ensure that the soil moisture at the top 1 m of the soil profile was more than 80% of the field capacity to promote successful germination for all the treatments. Base fertilizers of 225 kg ha−1 for N, 225 kg ha−1 for P2O5 and 225 kg ha−1 for K2O in the form of composite chemicals were incorporated into the soil. No other fertilizer was applied during the growing season.

2.3. Measurements

2.3.1. Crop Growth and Production

The time of sowing, harvesting, and the occurrence of the major growth stages were recorded. During the growing seasons, the plants were managed free of pest and disease and weeds were controlled manually. At the harvest, firstly, 100 plants were randomly collected from each plot to measure the panicle dry-weight per plant (PDW/plant). Next, all the plants in each plot were harvested manually and the grain and straw were dried to a constant weight at 80 °C to obtain the grain yield (GY) and straw yield (SY). Thousand grain weight (TGW) was measured manually.

2.3.2. Canopy Temperature and Root Sampling

At shoot emergence, anthesis and the grain-filling stages, the canopy temperature (CanopyT) was measured for each plot at an angle of 30° by radiation thermometer (HIOKI 3460-50, HIOKI E. E. COPPORATION, Nagano, Japan). The measurements were taken at 12:00–14:00 on clear days. The average of four measurements were taken as one replicate.
The roots of control variety Jigu19 were sampled in 2019 and 2020 seasons at maturity for both the WW and WS treatments. Soil cores were taken with a 9 cm inner-diameter stainless-steel soil drill down to 160 cm with increments of 20 cm. Three cores were taken for each treatment. The soil samples from the cores were collected with a mesh bag and washed to obtain the roots. The roots were scanned by the Expression 12000XL scanner, and the root length was measured by WinRHIZO (WinRHIZO Pro 32-bit 2013e). The root-length density (RLD) at different depths was calculated by dividing the root length by the sampled soil volume (cm/cm3).

2.3.3. Leaf Area, Leaf Chlorophyll Content and Rate of Water Loss

Five plants were selected for the determination of leaf area (LA); all the green leaves of each plant were measured using a measuring tape. The length and width of each leaf were measured and LA was estimated based on the following formula: leaf area = length × width × 0.75, in which the length and width were the maximum length and width [23].
The chlorophyll content (ChlC) was measured by chlorophyll meter (SPAD-502, Minolta Camera Co. Ltd., Osaka, Japan); each flag-leaf SPAD value obtained was the average of 10 readings (5 on each side of leaf at midrib), and 10 plants were measured for each plot from anthesis to the late grain-filling stage.
At 9:00–10:00 local time, three young fully expanded leaves from each plot in the field were sampled randomly and transported to a laboratory within 10 min. the laboratory was at a constant temperature and relative humidity. The leaves were surface dried quickly and weighed; the weight was initial weight (IWt) and the time (t0) was recorded. Then, the leaves were laid flat on the mesh. After 2 h on the mesh, the leaves were weighted; the weight was the final weight (FWt) and the time (t1) was also recorded. Then, the leaves were dried at 80 °C for 48 h to measure the dry weight (DWt).
Rate of water loss (RWL) = (IWt − FWt)/DWt/(t1 − t0)

2.3.4. Weather Factor Monitoring

Daily meteorological data, including daily maximum temperature (Tmax), daily minimum temperature (Tmin), accumulated temperature (AT), relative humidity (RH), and sunshine duration (Shr) were collected at a national weather station that was close to the experimental site during the study period. The diurnal temperature range (DTR) was calculated by subtracting the daily Tmin from the daily Tmax. The reference evapotranspiration (ET0) was calculated with the crop-water program that was developed by FAO using the FAO Penman–Monteith equation [24].
The average Tmin, Tmax, Shr, RH and DTR during the vegetative stage, reproductive stage and the whole growing stage; the total AT and ET0 during the vegetative stage, reproductive stage and the whole growing stage were calculated. The vegetative stage was defined as from sowing to booting, while the reproductive stage was defined as from booting to harvest. The duration for the two periods is shown in Table A2. The date at which more than 50% of the plants reached a particular growth stage (booting) was considered as their booting (growth stage) date.

2.3.5. Calculations

Harvest index is calculated as: Harvest index (HI) = GY/(GY + SY).
Abortive grain rate (AGR) is defined as: 1 − (grain weight per panicle/panicle weight per plant).
Water productivity is defined as grain yield divided by the seasonal evapotranspiration (ET), which was calculated using the following equation: ET = I + P + C + SWD − R − D, where I is irrigation; P is rainfall; C is capillary rise; SWD is soil water depletion (soil water contents at sowing, minus that at harvesting for the 2 m soil profile); R is runoff; and D is drainage from the root-zone profile. Due to the rainfall being prevented by the moveable shelter, no runoff was observed, irrigation did not result in the drainage, and capillary rise was not observed due to the deep groundwater level. P, R, D and C were taken as zero under the conditions of this study. Therefore, ET was calculated as I + SWD.
The drought resistance index (DRI) for each cultivar was calculated using the equation [25]: DRI = (GY2WS/GYWW) × (GYCK-WW/GYCK-WS)2, where GYWS and GYCK-WS are the yield of the newly collected and control variety, respectively, under WS treatment; GYWW and GYCK-WW are the yield of the newly collected and control variety, respectively, under WW treatment.

2.4. Data Analysis

The average values from three replicates for each variety were used for the statistical analysis; a one-way ANOVA was performed by SPSS 16.0 and a Pearson correlation analysis was used (SPSS Version 16.0) to assess the correlations between traits. A path analysis was performed by SPSS Version 16.0; linear regression was used to calculate the standardized coefficients (direct effect) and correlation coefficients, grain yield was the dependent variable, and agronomic and physiological traits were the independent variables. The indirect effect was the product of the standardized coefficients and correlation coefficients.

3. Results

3.1. Selecting Traits for Improving Yield under Climate Change

3.1.1. The Effects of Weather and Irrigation on Yield

The changes in yield of all foxtail millet varieties for each season under the two water treatments are shown in Figure 1. A large yield difference was observed between the varieties in each season under both water treatments. Under the WS treatment, the smallest difference (0.36 t/ha) and the largest difference (1.43 t/ha) were observed in 2018 and 2020, respectively. Under the WW treatment, the smallest difference (0.38 t/ha) and the largest difference (1.75 t/ha) were observed in 2019 and 2017, respectively. The results indicated that there was a large variation in the yield performance among different varieties. The yield of most varieties was significantly different from the control variety.
The average grain yield across varieties in each season varied from 3.50 (in 2018) to 5.77 (in 2020) t/ha under the WS treatment; from 5.01 (in 2017) to 6.95 (in 2020) t/ha under the WW treatment. The average grain yield of all varieties in each season decreased slightly from 2016 to 2017, then increased from 2018 to 2020, which indicated that the yield of millet had the potential to increase.
Compared with no irrigation, the average yield of all varieties in each season increase was around 25.66%, 13.19%, 43.01%, 40.40% and 20.39% irrigation in 2016, 2017, 2018, 2019 and 2020 seasons, respectively. The yield of foxtail millet was greatly improved with irrigation application compared with that without irrigation.

3.1.2. Weather Factors Related to Yield

The correlation of the yield with different weather factors during the vegetative, reproductive and the whole growth durations for variety Jigu19 that was grown from 2011 to 2020 was analysed (Table 1). Among the different weather factors, only temperature significantly affected GY under the WW treatment. High accumulated temperature and daily Tmax during the whole growth period had significant negative effects on GY, and high accumulated temperature and daily Tmax and daily Tmin during the vegetative growth period also had significant negative effects on GY under the WW treatment. Other weather factors did not significantly affect the GY of Jigu19. The results indicated that temperature was the main weather factor that influenced millet production; high temperature could impede high grain yield.

3.1.3. Climate Change and Trait Selection

The changes in the main weather factors in the past 15 years (from 2006 to 2020) at the experimental site are shown in Figure 2. There were large variations in the main weather factors. DTR, sunshine duration and temperature were in a slight increase trend and ET0 also showed a positive increase trend, while the RH showed a negative decreasing trend. The increased ET0 was related to the increased temperature and decreased relative humidity. The correlation analysis of yield with the main weather factors (Table 1) showed that decreased temperature favored yield formation, especially under well-watered conditions. The increasing temperature trend would negatively influence foxtail millet production in sufficient water supply conditions.
Grain yield in 2018 was the lowest, compared to the grain yield in other years (Figure 1), Table A3 showed that the Tmin in 2018 was the highest among the five seasons, especially during the vegetative stage. The AT and Tmax were not the highest, but they were relatively high. The high temperature could reduce the GY. Figure 2 shows the daily Tmin, daily Tmax and AT increase trends which indicate the possible negative effects of temperature on GY.

3.2. Traits Related to the Yield Improvement of Foxtail Millet under Different Water Supply Conditions

Table 2 lists the correlations of GY and DRI with the physiological and agronomic parameters of those varieties that were tested in 2016–2020. DRI is an important parameter for drought resistant estimation. Table 2 shows that the correlation between GY and DRI was stronger under the WS treatment than under the WW treatment; although the correlations were not statistically significant, DRI improvement could be beneficial to a high grain yield under water stress.
GY significantly positively correlated with PDW/plant, TGW, and LA under both WW and WS treatments, which indicated that high PDW/plant, TGW, and LA played a vital role in improving GY independent of the water conditions. Under the WW treatment, significant positive correlations between GY, ChlC, and CanopyT were found, while significant negative correlations were found between GY and AGR. Thus, high ChlC and CanopyT and low AGR were more important under sufficient water-supply conditions. However, under the WS treatment, there were significant positive correlations between GY and HI, and significant negative correlations between GY and RWL.
In order to obtain a clear picture of the interrelationships between different traits, the direct and indirect effects of different traits were worked out using a path analysis in respect of grain yield (Table 3). Under the WW treatment, the direct effect of TGW, PDW/plant and HI on grain yield were positive, and the direct effect of AGR on grain yield were negative. The indirect effects of TGW on AGR had the greatest positive influence on grain yield. Under the WS treatment, the direct effect of PDW/plant, LA and DRI on grain yield were positive, and the direct effect of ChlC on grain yield were negative. The indirect effects of ChlC on PDW/plant and LA had the greatest positive and negative influence on grain yield, respectively. The indirect effects of DRI on PDW/plant had a great positive influence on grain yield. The indirect effects of LA on PDW/plant and PDW/plant on LA had great negative influences on grain yield.
DRI significantly positively correlated with HI under the WS treatment, which could contribute to a high yield of foxtail millet under water stress conditions. However, DRI significantly negatively affected CanopyT under both WW and WS treatments. Low CanopyT could be beneficial to DRI improvement.
Root distribution in the soil profile is an important factor for efficient soil water utilization, especially under limited irrigations. Most roots occurred in the 0–40 cm soil layer. Under the WS treatment, the roots in the 0–40 cm soil layer accounted for 37.4% and 58.4% in 2019 and 2020, respectively (Figure 3). Under the WW treatment, this accounted for 57.6% and 83.8% in 2019 and 2020, respectively. In the 0–60 cm soil layer, root density was higher under WW than that under the WS treatment, while in the deeper soil layer, root density was higher under WS than under the WW treatment.
The correlations of root density in the upper (0–40 cm), middle (40–100 cm) and deep (100–160 cm) soil layers with water productivity of different varieties in 2019 and 2020 are shown in Table 4. There was a significant negative correlation between root density and water productivity in upper layer soil under both the WW and WS treatments. For root density in the middle and deep soil layers, no significant correlations were found between them and water productivity under both the WW and WS treatments. The results indicated that water productivity of foxtail millet closely related to the root system in the upper soil layer.
For GY, the root density in the upper soil layer significantly correlated with GY under the WS treatment, but it did not significantly correlate with GY under the WW treatment. There were no significant correlations between root density in the middle and deep soil and GY under both the WS and WW treatments. Under the WS treatment, water productivity in the deep soil significantly positively correlated with GY, but water productivity in the upper and middle soil significantly negatively correlated with GY. Under the WW treatment, water productivity in the deep soil significantly positively correlated with GY. High root density in deep soil would contribute to a high grain yield.
In terms of DRI, there were no significant correlations between root density and water productivity in the upper, middle and deep soil and DRI under both the WS and WW treatments. Thus, root density and water productivity did not significantly influence DRI.

4. Discussion

In this study, large seasonal variation in grain yield was observed under both irrigation and non-irrigation treatments (Figure 1), and there was a large variation in the yield of some varieties (Figure 2, Table A4). Selecting a better variety could alleviate some of the negative effects of water deficit and weather factors on foxtail millet production. A better understanding of how climate factors and water stress influence grain yield is important for developing measures to eliminate the negative impacts.
Previous studies have indicated that weather factors significantly affect millet production [26,27,28]. There was a strong negative effect of both rainfall and temperature on millet yield, suggesting that there was variation in the climatic requirement, there is no doubt that farmers increasingly need detailed maps of this kind to plan crop-planting schemes and monitor the rate of yield (smart agriculture) [28]. In this study, the main weather factor that influenced the grain yield of foxtail millet Jigu19 was identified as being temperature only under the WW treatment, and the influence of other weather factors was small. Low temperature was beneficial for grain yield, but future climate change might favor an increased temperature that impedes grain yield in water-sufficient conditions. The weather factors in this study did not significantly impact grain yield in water-deficient conditions. Therefore, weather factors were not the main factors that influence the grain yield of foxtail millet varieties, which is similar to Jigu19 in a dry place without irrigation.
Djanaguiraman et al. [29] found that high temperature stress (≥36/26 °C) imposed at different stages and durations caused a decrease in the number of seeds, individual seed weight and seed yield. Two periods (10–12 days and 2–0 days before anthesis) were identified as the most sensitive to short episodes of stress, causing maximum decreases in the pollen germination percentage and seed numbers of pearl millet. In this study, the grain yield of Jigu19 was more sensitive to the temperature during the vegetative stage than to the temperature during the reproductive stage under the WW treatment. High AT, Tmax and Tmin during the vegetative stage could decrease the grain yield of Jigu19. Therefore, millet varieties with characters that are less sensitive to temperature change during the vegetative stage would be preferred for the regions with sufficient water supply. However, trait × environment, trait × genotype, and trait × genotype × environment interactions are complicated, so further research at different experiment sites is needed.
The correlation analysis between grain yield and drought resistance index with physiological factors and agronomic traits could be used to select new varieties with traits under different water supply conditions, which could reduce the negative effects of water deficiency. Sanjana Reddy et al. [30] pointed out that the interaction between variety and environment was highly significant for the grain yield of pearl millet. In India, from the early 1950s to 2010, the development and diffusion of improved varieties along with suitable agronomic practices facilitated productivity, which increased grain yield from 430 kg ha−1 to 1079 kg ha−1 over a period of four decades [31]. In this study, the panicle dry weight per plant, thousand grain weight, leaf area and grain yield of different varieties of foxtail millet in five seasons had strong positive relationships under both the WW and WS treatments. High panicle dry weight per plant, thousand grain weight and leaf area could improve grain yield. However, abortive grain rate negatively affected grain yield and chlorophyll concentration positively affected grain yield only under the WW treatment; less abortive grains and high chlorophyll concentration could contribute to high grain yield, especially with sufficient water supply. Chlorophyll is a crucial component of crop photosynthesis machinery; high chlorophyll concentration could cause more photosynthate, which causes a low abortive grain rate.. Reza [18] also found a significant positive correlation between chlorophyll concentration and grain yield in both WW and WS conditions in wheat. However, Gu et al. [32] reported that a relatively low chlorophyll content increased photosynthesis and photosynthetic nitrogen-use efficiency. They also reported that it would avoid excessive light absorption, increase PSII efficiency, and result in a higher electron transport rate, and lowering chlorophyll content would not only benefit leaf photosynthesis, but also improve the light distribution in the canopy. Therefore, optimum chlorophyll concentration is better than too much or too little chlorophyll concentration in good water supply conditions.
In this study, CanopyT significantly positively affected grain yield under the WW treatment, while CanopyT significantly negatively affected the drought resistance index under the WW and WS treatment. In wheat, it is considered that wheat genotypes with a low canopy temperature can maintain a high transpiration and photosynthetic rate, as well as produce a high yield under water-stressed conditions [18]. Thus, optimal canopy temperature is preferred for a relatively high yield and drought resistance in sufficient water supply conditions. In addition, under the WS treatment, harvest index significantly positively affected grain yield and drought resistance index, while the rate of water loss significantly negatively affected grain yield. High harvest index and low rate-of-water-loss varieties selection could be considered for high grain yield, and a high harvest index is also beneficial to a high drought resistance index. Therefore, in dry areas, harvest index improvement and rate-of-water-loss reduction are very important for drought-resistant millet breeding or drought-resistant mechanism research.
Root development and distribution in the soil profile are of great importance for the water and nutrients uptake of foxtail millet. In the upper soil layer (0–40 cm), more roots were grown under WW than under the WS treatments, while in the deeper soil layer, more roots were grown under WS than under the WW treatments, because foxtail millet reorients its root growth toward deeper soil layers that retain more water in drought conditions [33]. The enhancement of water productivity in agriculture is very important to achieve efficient crop production [34]. In this study, grain yield significantly positively correlated with water productivity in deep soil layers under both the WS and WW treatments. However, grain yield significantly negatively correlated with water productivity in the upper and middle soil layers under only the WS treatment; grain yield significantly positively correlated with root density in the upper soil layer under only the WS treatment, which indicated that root density and water-use efficiency were more important than weather factors to millet production under water deficient condition, and high root density in the upper soil layer, high water productivity in the deep soil layer and low water productivity in the upper and middle soil layers could improve grain yield. However, in this study, the drought resistance index did not significantly correlate with root density and water productivity. Probably, for millet, root density and water productivity change could not influence drought resistance.

5. Conclusions

A large difference in yield among some varieties was observed, indicating the possibility of selecting better varieties to mitigate the negative effects of fluctuations in weather factors. For foxtail millet that was grown without water stress, temperature was the main weather factor influencing the grain yield of Jigu19; varieties with characters that are less sensitive to temperature change would be preferred.
The high panicle dry weight per plant, thousand grain weight and leaf area would improve the grain yield of foxtail millet under both irrigation treatments. Low abortion grain rate and optimal chlorophyll and canopy temperature could favor high grain yield and drought resistance under good water supply conditions. A high harvest index and low rate of water loss could favor high grain yield under limited water supply conditions.
The varieties with more roots in the upper soil layer could be beneficial to high grain yield, especially under water stress conditions. High water productivity in the deep soil layer was beneficial to high grain yield under good water supply and water stress conditions.

Author Contributions

Conceptualization, W.Z., B.W. and G.L.; formal analysis, C.B.; investigation, B.L. and Z.C.; resources, W.Z. and B.W.; writing—original draft preparation, C.B.; writing—review and editing, C.B. and Y.G.; supervision, W.Z. and B.W.; project administration, W.Z. and B.W.; funding acquisition, W.Z., B.W., C.B. and Y.G. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Key Research and Development Program, grant number 2018-YFD1000704-6; China Agricultural Research System, grant number CARS-06-14.5-A34; the HAAFS Agriculture Science and Technology Innovation Project, grant number 2022KJCXZX-HZS-3; Research Foundation for Talented Scholar, grant number G2019ZK40; G2019ZK45; and Guangxi Science and Technology Foundation for Talented Scholar, grant number GuikeAD20238077, GuikeAD20238066.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Table A1. The information of 25 varieties.
Table A1. The information of 25 varieties.
VarietyPedigreeRegions of AdaptionYear of ReleaseYear GrownYield
(t/ha)
Thousand Grain Weight (g)
Canggu9Heng968 × K492Hebei201820185.302.77
Heng1475(Henggu12 × Yugu18) × kn2009-2Hebei201920175.222.60
Heng2011-2Jigu17 × Henggu15Hebei, Jilin201820175.082049
Henggu13Yugu15 × Jigu31Hebei, Shandong, Henan, Liaoning, Jilin20182016, 20196.502.91
Henggu15Not clearHenan, Liaoning, Jilin201820167.532.88
Henggu16Shagu × K721-1Henan, Liaoning, Jilin201820175.402.68
Henggu17Shagu × K720-3Henan, Liaoning, Jilin201820175.622.69
Henggu18Not clearHenan, Liaoning, Jilin 20173.882.57
Henggu2018-1Henggu11 × Ji0506Hebei201820185.512.32
Henggu2018-2Henggu11 × (Hengyan5 × SR3522)Hebei201920185.562.32
Henggu2018-3(Hengyan5 × SR3522) × M1508Hebei201920184.492.68
Hengsi1Not clearHebei201820167.622.86
Jigu22Yugu9 × Jigu25Hebei, Shandong, Henannot clear20207.082.81
Jigu168Yugu18 × 1310-2Hebei, Liaoning, Neimongol, Shaanxi202020206.372.99
Jigu19Ai88 × QingfengguShanxi, Shandong, Henan, Shaanxi20042016–20205.482.49
Jigu34Not clearShanxi, Shandong, Henan, Shaanxi201320168.592.94
Jigu37Not clearShanxi, Shandong, Henan, Shaanxi 20168.913.13
Jigu39An09-8525 × [An4585 × (Jigu24 × 2010-M1445)]Hebei, Henan, Shandong, Shanxi, Xinjiang, Beijing, Liaoning, Jilin20182019, 20207.113.10
Jigu42An4585 × (Shi98622 × 1310-2)Hebei, Henan, Shandong, Xinjiang, Liaoning, Jilin, Neimongol, Shanxi, Shaanxi, Heilongjiang201820184.002.48
Lvmi16Not clearNot clearNot clear20167.473.27
Zhonggu10Zhonggu2 × Chuang877Beijing, Tianjin, Hebei, Shanxi, Neimongol202020206.923.05
Zhonggu2Yugu1 × Ai88Beijing, Tianjin, Hebei, Shanxi, Neimongol20152017–20206.002.72
Zhonggu7Not clearBeijing, Hebei, Henan, Shandong202020174.262.31
Zhonggu8Jigu20 × Q31Beijing, Hebei, Henan, Shandong 202020175.642.51
Zhonggu18Not clearBeijing, Hebei, Henan, Shandong 20182018–20206.072.44
Table A2. The recorded dates of sowing, booting and harvest.
Table A2. The recorded dates of sowing, booting and harvest.
YearSowing DateBooting DateHarvest Date
201621 June4 August27 September
201719 June5 August25 September
201818 June3 August23 September
201916 June4 August25 September
202023 June6 August29 September
Table A3. The major weather factors during the five seasons at different growing stages of foxtail millet.
Table A3. The major weather factors during the five seasons at different growing stages of foxtail millet.
Weather FactorsSeasonsWhole
Stage
Vegetative
Stage
Reproductive
Stage
Accumulated temperature (°C)20162508.801314.701194.10
20172546.001350.001196.00
20182603.401396.101207.30
20192616.601408.401208.20
20202541.001312.901228.10
Minimum temperature (°C)201621.0722.3119.83
201721.4622.7720.20
201822.3023.9220.69
201922.0623.7020.41
202021.6022.1621.03
Maximum temperature (°C)201630.9432.0329.86
201731.4032.9029.94
201831.7033.8129.66
201932.0434.5029.58
202031.0731.8030.34
Diurnal temperature range (°C)20169.879.7210.02
20179.9410.139.74
20189.419.858.97
20199.9910.809.17
20209.479.649.31
Reference evaporation (mm)2016481.38263.91217.47
2017514.67292.71221.95
2018511.89296.15215.74
2019493.78288.57205.21
2020487.74281.94205.80
Sunshine duration (hr/d)20166.776.277.27
20177.477.087.87
20186.677.336.02
20196.436.925.93
20207.468.506.43
Relatively humidity (%)201678.1575.6780.63
201774.6870.1279.24
201872.8870.0675.69
201964.8260.1469.49
202069.8266.2273.41
Table A4. The maximum, minimum, mean and coefficient of variation of grain yield of different varieties in five seasons.
Table A4. The maximum, minimum, mean and coefficient of variation of grain yield of different varieties in five seasons.
TreatmentYearMaximumMinimumMeanCoefficient of Variation (%)
WW20166.685.606.028.07
20175.643.895.0113.06
20185.564.005.0112.49
20196.606.226.412.89
20207.366.376.955.23
WS20165.293.894.7912.97
20175.234.044.4311.78
20183.673.303.504.26
20194.864.114.577.75
20206.194.765.779.50

References

  1. Bala-Ravi, S. Neglected millets that save the poor from starvation. LEISA India 2004, 6, 1–8. [Google Scholar]
  2. Amadou, I.; Gounga, M.E.; Shi, Y.H.; Le, G.W. Fermentation and heat-moisture treatment induced changes on the physicochemical properties of foxtail millet (Setaria italica) flour. Food Bioprod. Process. 2014, 92, 38–45. [Google Scholar] [CrossRef]
  3. FAOSTAT (Food and Agriculture Organization Statistics). Available online: http://faostat.fao.org/site/567/DesktopDefault.aspx?PageID=567 (accessed on 20 July 2016).
  4. Koch, D.W. Foxtail Millet: Management for Supplemental and Emergency Forage; B-1122.3; Cooperative Extension Service; University of Wyoming: Laramie, WY, USA, 2002. [Google Scholar]
  5. Liu, T.; Ye, N.; Song, T.; Cao, Y.; Gao, B.; Zhang, D.; Zhu, F.; Chen, M.; Zhang, Y.; Xu, W.; et al. Rhizosheath formation and involvement in foxtail millet (Setaria italica) root growth under drought stress. J. Integr. Plant Biol. 2019, 61, 449–462. [Google Scholar] [CrossRef] [PubMed]
  6. Pan, J.; Li, Z.; Wang, Q.; Pan, J.; Garrell, A.K.; Liu, M.; Guan, Y.; Zhou, W.; Liu, W. Comparative proteomic investigation of drought responses in foxtail millet. BMC Plant Biol. 2018, 18, 315. [Google Scholar] [CrossRef] [PubMed]
  7. Yu, T.; Zhao, W.; Fu, J.; Liu, Y.; Chen, M.; Zhou, Y.; Ma, Y.; Xu, Z.; Xi, Y. Genome-wide analysis of cdpk family in foxtail millet and determination of SiCDPK24 functions in drought stress. Front. Plant Sci. 2018, 9, 651. [Google Scholar] [CrossRef] [PubMed]
  8. Tripathi, A.; Tripathi, D.K.; Chauhan, D.; Kumar, N.; Singh, G. Paradigms of climate change impacts on some major food sources of the world: A review on current knowledge and future prospects. Agric. Ecosyst. Environ. 2016, 216, 356–373. [Google Scholar] [CrossRef]
  9. Howden, S.M.; Soussana, J.F.; Tubiello, F.N.; Chhetri, N.; Dunlop, M.; Meinke, H. Adapting agriculture to climate change. Proc. Natl. Acad. Sci. USA 2007, 104, 19691–19696. [Google Scholar] [CrossRef] [Green Version]
  10. Solomon, S.; Qin, D.; Manning, M.; Marquis, M.; Averyt, K.; Tignor, M.M.B.; Miller, H.L., Jr. Climate Change 2007: The Physical Science Basis. Working Group I Contribution to the Fourth Assessment Report of the IPCC; Cambridge University Press: Cambridge, UK, 2007; Volume 4. [Google Scholar]
  11. Resende, R.T.; Piepho, H.P.; Siva-Junior, B.O.; Silva, F.F.; Resende, M.D.V.; Grattapaglia, D. Enviromics in breeding: Applications and perspectives on envirotypic-assisted selection. Theor. Appl. Genet. 2020, 134, 95–112. [Google Scholar] [CrossRef]
  12. Crossa, J.; Fritsche-Neto, R.; Montesinos-Lopez, O.A.; Costa-Neto, G.; Dreisigacker, S.; Montesinos-Lopez, A.; Bentley, R.A. The modern plant breeding triangle: Optimizing the use of genomics, phenomics, and enviromics data. Front. Plant Sci. 2021, 12, 651480. [Google Scholar] [CrossRef]
  13. Porter, J.R.; Xie, L.; Challinor, A.J.; Cochrane, K.; Howden, S.M.; Iqbal, M.M.; Lobell, D.B.; Travasso, M.I. Food security and food production systems. In Climate Change 2014: Impacts, Adaptation, and Vulnerability. Part A: Global and Sectoral Aspects. Contribution of Working Group II to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change; Cambridge University Press: Cambridge, UK; New York, NY, USA, 2014; pp. 485–533. [Google Scholar] [CrossRef] [Green Version]
  14. Jalihal, C.; Srinivasan, J.; Chakraborty, A. Modulation of Indian monsoon by water vapor and cloud feedback over the past 22,000 years. Nat. Commun. 2019, 10, 5701. [Google Scholar] [CrossRef] [Green Version]
  15. Krishna, S.S.; Reddy, Y.A.N.; Kumar, R.L.R. Assessment of traits for grain yield under drought in finger millet. Plant Physiol. Rep. 2021, 26, 84–94. [Google Scholar] [CrossRef]
  16. Robin, S.; Pathan, M.S.; Courtois, B.; Lafitte, R.; Carandang, S.; Lanceras, S.; Amante, M.; Nguyen, T.H.; Li, Z. Mapping osmotic adjustment in an advanced back-cross inbred population of rice. Theor. Appl. Genet. 2003, 107, 1288–1296. [Google Scholar] [CrossRef] [PubMed]
  17. Rashid, A.; Stark, J.C.; Tanveer, A.; Mustafa, T. Use of canopy temperature measurements as a screening tool for drought tolerance in spring wheat. J. Agron. Crop Sci. 1999, 182, 231–237. [Google Scholar] [CrossRef]
  18. Talebi, R. Evaluation of chlorophyll content and canopy temperature as indicators for drought tolerance in durum wheat (triticum durum desf.). Aust. J. Basic Appl. Sci. 2011, 5, 1457–1462. [Google Scholar] [CrossRef]
  19. Vadez, V.; Hash, T.; Bidinger, R.F.; Kholova, J., II. 1.5 Phenotyping pearl millet for adaptation to drought. Front. Physiol. 2012, 3, 386. [Google Scholar] [CrossRef] [Green Version]
  20. Porter, J.R.; Gawith, M. Temperatures and the growth and development of wheat: A review. Eur. J. Agron. 1999, 10, 23–36. [Google Scholar] [CrossRef]
  21. Smith, S.; De Smet, I. Root system architecture: Insights from Arabidopsis and cereal crops. In: Philosophical Transactions of the Royal Society of London. Biol. Sci. 2012, 367, 1441–1452. [Google Scholar] [CrossRef] [Green Version]
  22. Liu, Z.; Zhang, T.; Li, C.; Bai, G. Genetic diversity and classification of cytoplasm of Chinese elite foxtail millet [Setaria italica (L.) P. Beauv.] germplasm. Crop Sci. 2014, 54, 659–666. [Google Scholar] [CrossRef] [Green Version]
  23. Zhang, D.; Panhwar, R.B.; Liu, J.; Gong, X.; Liang, J.; Minxuan, L.I.U.; Lu, P.; Gao, X.; Feng, B. Morphological diversity and correlation analysis of phenotypes and quality traits of proso millet (Panicum miliaceum L.) core collections. J. Integr. Agric. 2019, 18, 958–969. [Google Scholar] [CrossRef]
  24. Allen, R.G.; Pereira, L.S.; Raes, D.; Smith, M. Crop Evapotranspiration: Guidelines for computing crop water requirements. In FAO Irrigation and Drainage 56; FAO: Rome, Italy, 1998; Volume 300. [Google Scholar]
  25. Zhang, W.; Diao, X.; Jia, G.; Liu, B.; Wang, B.; Zhi, H.; Pang, J. National Agriculture Standards of P.R. of China: Technical Cose of Practice for Identification of Drought Resistance Varieties in Foxtail Millet; China Agriculture Press: Beijing, China, 2021; NY/T 3851-2021. [Google Scholar]
  26. Qu, Y.; Zhou, Y.; Chao, G.; Zhang, P.; Song, H.; Feng, B. Water use efficiency of foxtail millet (Panicum italicum L.) under climate change conditions in Northwest regions of China. Agrociencia 2016, 50, 665–676. [Google Scholar]
  27. Guo, D.; Chen, C.; Guo, P.; Yuan, X.; Xing, X.; Ma, X. College of wat evaluation of AquaCrop model for foxtail millet (setaria italica) growth and water use with plastic film mulching and no mulching under different weather conditions. Water 2018, 10, 836. [Google Scholar] [CrossRef] [Green Version]
  28. Ojo, O.I.; Olaniyan, A.O.; Gbadegesin, A.S.; Ilunga, M.F. Assessment of Climatic Variability Effect on Millet Production and Yield. In Handbook of Climate Change Resilience, 1st ed.; Filho, L.W., Ed.; Springer: Cham, Switzerland, 2020; pp. 1269–1289. [Google Scholar] [CrossRef]
  29. Djanaguiraman, M.; Perumal, R.; Ciampitti, I.A.; Gupta, S.K.; Prasad, P.V.V. Quantifying pearl millet response to high temperature stress: Thresholds, sensitive stages, genetic variability and relative sensitivity of pollen and pistil. Plant Cell Environ. 2018, 41, 993–1007. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  30. Sanjana Reddy, P.; Satyavathi, C.T.; Khandelwal, V.; Patil, H.T.; Gupta, P.C.; Sharma, L.D.; Mungra, K.D.; Singh, S.P.; Narasimhulu, R.; Bhadarge, H.H.; et al. Performance and stability of pearl millet varieties for grain yield and micronutrients in arid and semi-arid regions of India. Front Plant Sci. 2021, 12, 670201. [Google Scholar] [CrossRef] [PubMed]
  31. Kumara, C.D.; Bantilan, C.; Rajalaxmi, A.; Rai, K.N.; Yadav, O.P.; Gupta, S.K.; Singh, P.N. Development and Diffusion of Pearl Millet Improved Cultivars in India: Impact on Growth and Yield Stability; ICRISAT: Patancheru, India, 2014. [Google Scholar] [CrossRef]
  32. Gu, J.; Zhou, Z.; Li, Z.; Chen, Y.; Wang, Z.; Zhang, H. Rice (Oryza sativa L.) with reduced chlorophyll content exhibit higher photosynthetic rate and efficiency, improved canopy light distribution, and greater yields than normally pigmented plants. Field Crops Res. 2017, 200, 58–70. [Google Scholar] [CrossRef]
  33. Faye, A.; Sine, B.; Chopart, J.-L.; Grondin, A.; Lucas, M.; Diedhiou, A.G.; Pascal, G.; Laurent, C.; Doohong, M.; Alain, A.; et al. Development of a model estimating root length density from root impacts on a soil profile in pearl millet (Pennisetum glaucum (L.) R. Br). Application to measure root system response to water stress in field conditions. PLoS ONE 2019, 14, e0214182. [Google Scholar] [CrossRef] [Green Version]
  34. Ali, M.H.; Talukder, M.S.U. Increasing water productivity in crop production—A synthesis. Agric. Water Manag. 2008, 95, 1201–1213. [Google Scholar] [CrossRef]
Figure 1. The changes in the highest and the lowest yield among the foxtail millet varieties each season and the average yield of the varieties under well-watered and water-stressed treatments for the five growing seasons (the Difference was defined as the difference between the varieties with the highest and lowest yield).
Figure 1. The changes in the highest and the lowest yield among the foxtail millet varieties each season and the average yield of the varieties under well-watered and water-stressed treatments for the five growing seasons (the Difference was defined as the difference between the varieties with the highest and lowest yield).
Agronomy 12 01500 g001
Figure 2. Changes in weather factors during the 5 growing seasons of foxtail millet and the recent 15 years from 2006 to 2020 (Sunshine: sunshine duration; Tmax: maximum temperature; Tmin: minimum temperature; accumulated temperature; DTR: average daily diurnal temperature range; (A). RH: relative humidity; ET0: reference evaporation (B)).
Figure 2. Changes in weather factors during the 5 growing seasons of foxtail millet and the recent 15 years from 2006 to 2020 (Sunshine: sunshine duration; Tmax: maximum temperature; Tmin: minimum temperature; accumulated temperature; DTR: average daily diurnal temperature range; (A). RH: relative humidity; ET0: reference evaporation (B)).
Agronomy 12 01500 g002
Figure 3. The distribution of the total root length (TRL) along the soil profile in 2019 for four varieties and 2020 for six varieties under well-watered (WW) and water-stressed (WS) treatments. (Bars represent the range of the RLD among the cultivars sampled).
Figure 3. The distribution of the total root length (TRL) along the soil profile in 2019 for four varieties and 2020 for six varieties under well-watered (WW) and water-stressed (WS) treatments. (Bars represent the range of the RLD among the cultivars sampled).
Agronomy 12 01500 g003
Table 1. Correlation analysis of grain yield of foxtail millet under water-stressed and well-watered conditions with weather factors during vegetative, reproductive and whole growing stages for variety Jigu19.
Table 1. Correlation analysis of grain yield of foxtail millet under water-stressed and well-watered conditions with weather factors during vegetative, reproductive and whole growing stages for variety Jigu19.
Weather FactorsGY (Well-Watered)GY (Water-Stressed)
WhS VS RS WhS VS RS
AT−0.697 *−0.843 **−0.212−0.471−0.109−0.499
Tmin−0.609−0.809 **−0.273−0.468−0.279−0.457
Tmax−0.637 *−0.774 **−0.152−0.366−0.029−0.492
DTR−0.225−0.4340.1030.0080.231−0.229
ET0−0.434−0.517−0.258−0.294−0.084−0.411
Shr−0.157−0.219−0.075−0.1830.035−0.264
RH0.6080.5320.513−0.008−0.2790.277
Note: AT: accumulated temperature; DTR: diurnal temperature range; Tmin: minimum temperature; Tmax: maximum temperature; ET0: reference evapotranspiration; Shr: sunshine duration; RH: relative humidity; WhS: the whole stage; VS: vegetative stage; RS: reproductive stage. * and **: significant at p < 0.05 and p < 0.01, respectively.
Table 2. Correlation analysis of agronomic and physiological traits with the grain yield and drought resistance for foxtail millet varieties grown under two irrigation treatments for five seasons.
Table 2. Correlation analysis of agronomic and physiological traits with the grain yield and drought resistance for foxtail millet varieties grown under two irrigation treatments for five seasons.
ITr GYPDW/PlantTGWAGRHILAChlCCanopyTRWL
WWGY 0.405 *0.745 **−0.729 **0.0540.563 **0.374 *0.365 *−0.311
DRI−0.1510.008−0.1940.3090.320−0.3160.198−0.366 **0.095
WSGY 0.495 **0.629 **0.1710.430 *0.530 **−0.006−0.029−0.505 **
DRI0.3540.347−0.0380.0740.552 **−0.2330.085−0.365 *−0.059
Note: ITr: irrigation treatment; WW: well-watered; WS: water-stressed; GY: grain yield; DRI: drought resistance index; PDW/plant: panicle dry weight/plant; TGW: thousand grain weight; HI: harvest index; AGR: abortive grain rate; LA: leaf area; ChlC: chlorophyll content; CanopyT: average canopy temperature; RWL: rate of water loss. * and **: significant at p < 0.05 and p < 0.01, respectively.
Table 3. Direct and indirect effects of factors influencing grain yield in foxtail millet under two irrigation treatments for five seasons.
Table 3. Direct and indirect effects of factors influencing grain yield in foxtail millet under two irrigation treatments for five seasons.
ITrTraitDirect Effect Indirect Effect through
TGWAGRPDW/PlantHI
WWTGW0.274-0.4030.073−0.021
AGR−0.632−0.175-0.081−0.005
PDW/plant0.5350.037−0.096-0.058
HI0.221−0.0060.013−0.174-
WS PDW/PlantLADRIChlC
PDW/plant0.725-−0.2130.096−0.113
LA0.696−0.222-−0.0640.120
DRI0.2760.252−0.162-−0.012
ChlC−0.2280.360−0.3680.108-
Note: ITr: irrigation treatment; WW: well-watered; WS: water-stressed; GY: grain yield; DRI: drought resistance index; PDW/plant: panicle dry weight/plant; TGW: thousand grain weight; HI: harvest index; AGR: abortive grain rate; LA: leaf area; ChlC: chlorophyll content.
Table 4. Correlation analysis of root density with water productivity and grain yield for foxtail millet varieties grown under two treatments in 2019 and 2020.
Table 4. Correlation analysis of root density with water productivity and grain yield for foxtail millet varieties grown under two treatments in 2019 and 2020.
Treatment Water ProductivityGrain
Yield
Drought
Resistance
Index
Upper
Soil
Middle
Soil
Deep
Soil
Water-
stressed
Root densityupper soil−0.804 ** 0.633 *0.218
middle soil 0.350 −0.454−0.348
deep soil 0.4300.2880.128
Grain yield −0.826 **−0.799 **0.745 *
Drought resistance index−0.440−0.4990.398
Well-
watered
Root densityupper soil−0.722 * 0.5190.271
middle soil 0.096 −0.512−0.211
deep soil −0.003−0.1600.553
Grain yield −0.592−0.0060.789 **
Drought resistance index0.234−0.0970.380
* and **: significant at p < 0.05 and p < 0.01, respectively.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Zhang, W.; Wang, B.; Liu, B.; Chen, Z.; Lu, G.; Ge, Y.; Bai, C. Trait Selection for Yield Improvement in Foxtail Millet (Setaria italica Beauv.) under Climate Change in the North China Plain. Agronomy 2022, 12, 1500. https://0-doi-org.brum.beds.ac.uk/10.3390/agronomy12071500

AMA Style

Zhang W, Wang B, Liu B, Chen Z, Lu G, Ge Y, Bai C. Trait Selection for Yield Improvement in Foxtail Millet (Setaria italica Beauv.) under Climate Change in the North China Plain. Agronomy. 2022; 12(7):1500. https://0-doi-org.brum.beds.ac.uk/10.3390/agronomy12071500

Chicago/Turabian Style

Zhang, Wenying, Bianyin Wang, Binhui Liu, Zhaoyang Chen, Guanli Lu, Yaoxiang Ge, and Caihong Bai. 2022. "Trait Selection for Yield Improvement in Foxtail Millet (Setaria italica Beauv.) under Climate Change in the North China Plain" Agronomy 12, no. 7: 1500. https://0-doi-org.brum.beds.ac.uk/10.3390/agronomy12071500

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop