1. Introduction
Wheat (
Triticum aestivum L.), as one of the world’s important food crops, accounts for about one third of human daily food demand [
1,
2,
3]. Facing the challenge of continuous population growth, undoubtedly, much more wheat production is required to satisfy the increasing demand [
4,
5]. To date, the Guinness World Record for the highest wheat grain yield is 17.95 t/ha [
6]. However, large yield gaps, i.e., the difference between the practical farmers’ yields achieved in large areas and the potential highest yields that can be achieved by using the best soil-crop management practices with the best adapted cultivar, exist in in many regions of the world [
7]. For example, the average wheat grain yield of China in 2023 is only 5.83 t/ha [
8]. Zinc (Zn) is an essential micronutrient for wheat growth and human health [
9]. However, wheat grains generally contain low Zn, which can’t meet the daily human nutrition requirement [
10]. Chen et al. [
11] investigated wheat grain Zn concentrations in seven major wheat production provinces of China and found an average of only 23.3 mg/kg, which is far lower than the recommended biofortification target value of 40-50 mg/kg [
12]. Therefore, while wheat grain yield needs to increase, it must simultaneously increase its Zn concentration to achieve both food and nutrition security.
Unfortunately, negative yield-nutrition trade-offs in wheat grains were frequently reported. High-yielding wheat cultivars generally had low gluten/protein concentrations [
13], and wheat lines with low phytate, an antinutritional compound that reduces grain Zn bioavailability, were always along with low grain yields [
14]. With the variety release years increasing from 1838 to 2012 in the UK, concentrations of Zn, iron (Fe) and protein in wheat grains were all unintentionally decreased, indicating modern wheat varieties with higher grain yields than old ones have lower nutritional quality in minerals [
15,
16]. Thus, the “Green Revolution” was actually constrained by the yield-nutrition “dilution effect” [
17,
18,
19], and achieving high yield and high nutritional quality simultaneously was an unavoidable dilemma and challenge for wheat breeders [
20,
21]. Actually, the “dilution effect” was reported to be mainly due to wheat breeding but seldom due to agronomic management practices, which were rarely considered in the studies above-mentioned [
22]. Whether agronomic practices, e.g., fertilization, can effectively overcome the yield-nutrition dilution effect or not is an emerging important question worthy to be answered [
23,
24].
Increasing Zn supply via soil fertilization and/or foliar spraying could effectively correct or prevent the symptom occurrence of Zn deficiency, ensure sufficient Zn uptake by wheat, and improve grain yields and Zn concentrations, particularly in calcareous soils with high pH, and low soil organic matter and moisture [
10]. Notably, many previous studies have shown that foliar Zn spraying is much more effective and much more economically efficient (low dose) than soil Zn application in increasing Zn concentration and bioavailability in the wheat grain/flour [
23]. In the HarvestZinc international study involving 7 countries (China, India, Kazakhstan, Mexico, Pakistan, Turkey, and Zambia), 23 experimental sites, 10 different wheat varieties, and 3 years, the grain Zn concentration was increased by an average of 83.5% achieved by foliar Zn spraying, and by only 12.0% by soil Zn application [
25]. Numerous studies have shown that increasing nitrogen (N) supply moderately (not excessively) increases wheat grain yields as well as improves grain Zn and N concentrations, i.e., the “N-Zn synergism”, but the input of phosphorus (P) fertilizer generally reduces crop grain Zn concentrations, being termed the “P-Zn” antagonism [
12,
23,
26].
Therefore, the optimized N and P fertilizers in combination with appropriate soil and foliar Zn applications may act as a double-win strategy to simultaneously achieve high wheat grain yields and high grain Zn concentrations/bioavailability [
23]. However, most previous researches have focused on the dose effects of a single element (N, P, or Zn), but less on their chemical forms and the integration or interactions among these elements, and most studies on foliar Zn spraying were conducted using small watering cans or knapsack sprayers with high labor costs and low efficiency, resulting in a clear separation between the experimental work and the farmers’ practices [
12]. There is a lack of research on compound, slow/controlled release, and organic/microbial fertilizers, and water and fertilizer integration technology adopted by farmers, and a lack of research on the development of new types of high-efficient foliar Zn fertilizers and foliar spraying techniques using modern agricultural drones.
In addition to Zn, other micronutrient elements, such as Fe, manganese (Mn) and copper (Cu), and macronutrients including N, P, potassium (K), calcium (Ca) and magnesium (Mg), are also essential nutrient elements determining wheat grain yields and nutritional quality and human dietary health [
22,
27,
28,
29,
30]. Most studies on wheat grain Zn biofortification focus on only a few of these elements (no more than 3 in most situations), but less on the whole suite of these mineral nutrients in wheat grains. Consequently, there is a lack of systematic understanding on the effects of different crop management practices on changes of these micro- and macro- nutrients in wheat grains, their cross-talks among each other and relationships with grain yields.
In this study, we changed the component ratio of N, P2O5, K2O and micronutrient in compound fertilizers through enlarging N or micronutrient and shrinking P in the soil fertilization experiment, sprayed foliar solutions with different chemical forms of Zn or without Zn using an agricultural drone, and investigated their effects on (1) wheat grain yields, yield components and other agronomic traits including the plant height (PH), spike length (SL), spike number (SN) per 666.7 m2, kernel number per spike (KNPS), thousand kernel weight (TKW), grain yield (Y), and harvest index (HI); (2) grain micronutrient accumulation including Zn, Fe, Mn and Cu; (3) changes of grain macronutrients (N, P, K, Ca, Mg); and (4) relationships among the agronomic and nutritional traits across different soil fertilization or foliar spraying treatments. The bioavailability of Zn and Fe, estimated by molar ratios of phytic acid (PA)/Zn, PA/Fe, PA × Ca/Zn and PA × Ca/Fe in wheat grains in the soil fertilization experiment was also investigated. These findings would provide a better understanding of the yield-nutrition relationship among agronomic practices, wheat yields and grain nutritional quality to achieve food security in quantity and quality, and biofortification of wheat grains with micronutrient (especially for Zn) to alleviate malnutrition.
2. Materials and Methods
2.1. Study Site
Field experiments were conducted at two sites/years during the winter wheat growing season from October to June in a winter wheat-summer maize rotation system. The soil fertilization treatments were conducted at Jiyang Experimental Station, Shandong Academy of Agricultual Sciences, Jinan, China, during 2020-2021. The foliar spraying treatments were conducted at Liuyuan planting base, Maifeng Wheat Planting Professional Cooperative in Dongming County, Heze, China, during 2022-2023. The two sites are located in west of Shandong Province of China, the area has a typical continental and warm-temperate monsoon climate, with cold and dry spring and winter, and a hot and rainy summer. The annual mean temperature is 12.0-14.0 °C and the annual frost-free period is 195-220 days. The annual precipitation is 500-700 mm, with 70% rainfall occurring during June-September. Detailed site information including geographic coordinates and soil basal properties in 0-20 cm soil layers of the two experimental sites before wheat sowing are presented in
Table 1.
2.2. Experimental Design and Crop Management
For the soil fertilization experiment, the single-factor randomized block design was applied with 5 treatments and 4 replicates. Corresponding to 5 treatments, five different compound fertilizers, (NPK 15-15-15 as a control, 15-15-15+Micronutrient, 17-17-17, 26-10-15 and 30-10-11, Stanley Agricultural Group Co., Ltd., Linshu, China), were evenly distributed and incorporated into the upper 20 cm of the soil prior to wheat planting, respectively, with a quantity of 750 kg·ha-1. The planting area of each treatment was 400 m2 (10 m × 40 m).
For the foliar spraying experiment, the single-factor randomized block design was applied with 3 treatments (
Table 2) and 3 replicates. These three treatments included: (1) spraying of deionized water as a control (CK); (2) spraying of a mixed solution with deionized water and ZnO (2.0%,
w/
v); and (3) spraying of a mixed solution with deionized water and ZnSO
4·7H
2O (2.0%,
w/
v), labeled Zn. All solutions contained 0.01% (
v/
v) TWEEN 20 as a surfactant and each time a dosage of 10 L/ha was sprayed to the area of 660 m
2 (44 m × 15 m) for each treatment using an agricultural drone (DJI AGRAS T40, DJI Agriculture, Shenzhen, China). All foliar spraying treatments were conducted three times. The spraying occurred for the first time 5 days after wheat flowering on May 8, and was repeated on May 17 and May 23, respectively, at 6-9 day intervals (
Table 2).
In addition, for each treatment of the two experiments above-mentioned, a 112.5 kg of N ha-1 (supplied as urea) was top-dressed with irrigation/rainfall during the regreening-jointing stage. Winter wheat (Triticum aestivum L.) variety “Jimai 22” was sown around 20 October in autumn and harvested around 10 June in the following year. All plots were adequately irrigated and no obvious biotic (weeds, pests, disease, etc.) and abiotic (drought, cold damage, etc.) stress was observed during the wheat growing season. To control aphids, omethoate (2-dimethoxyphosphinoylthio-N-methylacetamide) was sprayed at the booting stage, and no fungicide was applied.
2.3. Plant Sampling and Nutrient Analysis
At maturity, a 1-m2 area of spikes in each plot of the two experiments was manually harvested to determine the wheat grain yield and yield components (SN: spike number, KNPS: kernel number per spike, and TKW: thousand kernel weight). In addition, 10 random plants of each plot were used to determine the average above-ground plant height (PH) and spike length (SL).
After determination of the grain yield and yield components, grain samples were rapidly washed with deionized water, oven-dried at 65 °C for 72 h until constant weight, and then ground with a stainless steel grinder. Sub-samples after ground were digested with HNO
3-H
2O
2 in a microwave accelerated reaction system (CEM Corp., Matthews, North Carolina, USA). The concentrations of Zn, Fe, Mn, Cu, P, K, Ca and Mg in the digested solutions were determined by inductively coupled plasma optical emission spectroscopy (ICP-OES, Avio
TM 200, PerkinElmer, Waltham, Massachusetts, USA). A certified reference grain material (IPE556, Wageningen University) and two blanks were included in each digestion batch to ensure analytical quality. The grain N concentration was determined by the H
2SO
4-H
2O
2 digestion-micro-Kjeldahl method. Phytate-P concentration was determined by the method of Haug and Lantzsch [
31], and converted to PA by dividing by 0.282. The molar ratios of PA/Zn, PA × Ca/Zn, PA/Fe and PA × Ca/Fe were calculated to predict the bioavailability of Zn and Fe in wheat grains.
2.4. Statistical Analysis
Data from the single-factor randomized block design experiments were subjected to one-way analysis of variance (ANOVA) using The SAS System for Windows V8 (SAS Institute Inc., Cary, North Carolina, USA) and the mean value was compared according to Fisher’s protected least significant difference (LSD, p ≤ 0.05). Pearson’s correlation analysis and principal component analysis (PCA) were performed by OriginPro 2021 (OriginLab Corp., Northampton, Massachusetts, USA).
Figure 1.
Principle component analysis (PCA) of the effects of different soil fertilization treatments (a) and foliar spraying treatments (b) on various investigated parameters of wheat plants. 15-15-15, 17-17-17, 26-10-15, and 30-10-11 are ratios of N-P2O5-K2O in compound fertilizers (a). In panel (b), CK: spraying of deionized water; ZnO: spraying of a mixed solution with deionized water and ZnO; Zn: spraying of a mixed solution with deionized water and ZnSO4·7H2O. The abbreviations of various parameters investigated are as follows: yield (Y), plant height (PH), spike length (SL), spike number per 666.7 m2 (SN), kernel number per spike (KNPS), thousand kernel weight (TKW), HI (harvest index), concentrations of Zn, Fe, Mn, Cu, N, P, K, Ca, Mg and phytate-P, ratios of phytate-P/P, and molar ratios of phytic acid (PA)/Zn, PA × Ca/Zn, PA/Fe and PA × Ca/Fe in wheat grains.
Figure 1.
Principle component analysis (PCA) of the effects of different soil fertilization treatments (a) and foliar spraying treatments (b) on various investigated parameters of wheat plants. 15-15-15, 17-17-17, 26-10-15, and 30-10-11 are ratios of N-P2O5-K2O in compound fertilizers (a). In panel (b), CK: spraying of deionized water; ZnO: spraying of a mixed solution with deionized water and ZnO; Zn: spraying of a mixed solution with deionized water and ZnSO4·7H2O. The abbreviations of various parameters investigated are as follows: yield (Y), plant height (PH), spike length (SL), spike number per 666.7 m2 (SN), kernel number per spike (KNPS), thousand kernel weight (TKW), HI (harvest index), concentrations of Zn, Fe, Mn, Cu, N, P, K, Ca, Mg and phytate-P, ratios of phytate-P/P, and molar ratios of phytic acid (PA)/Zn, PA × Ca/Zn, PA/Fe and PA × Ca/Fe in wheat grains.
Figure 2.
Corrplot representing correlations among measured grain yields, yield components, other agronomic traits, and grain nutritional parameters of wheat crop across different soil fertilization treatments (a) and across different foliar spraying treatments (b). Positive correlations are displayed in blue and negative correlations in red. The color legend on the right-hand side of corrplot shows correlation coefficients and the corresponding colors. The intensity of the color is proportional to the correlation coefficient, and the ellipse size indicates the range of scattered experimental data points. “*”, “**” and “***” indicate significant correlations at p ≤ 0.05, 0.01 and 0.001, respectively. The abbreviations are as follows: yield (Y), plant height (PH), spike length (SL), spike number per 666.7 m2 (SN), kernel number per spike (KNPS), thousand kernel weight (TKW), harvest index (HI), concentrations of Zn, Fe, Mn, Cu, N, P, K, Ca, Mg and phytate-P, ratios of phytate-P/P, and molar ratios of phytic acid (PA)/Zn, PA × Ca/Zn, PA/Fe and PA × Ca/Fe in wheat grains.
Figure 2.
Corrplot representing correlations among measured grain yields, yield components, other agronomic traits, and grain nutritional parameters of wheat crop across different soil fertilization treatments (a) and across different foliar spraying treatments (b). Positive correlations are displayed in blue and negative correlations in red. The color legend on the right-hand side of corrplot shows correlation coefficients and the corresponding colors. The intensity of the color is proportional to the correlation coefficient, and the ellipse size indicates the range of scattered experimental data points. “*”, “**” and “***” indicate significant correlations at p ≤ 0.05, 0.01 and 0.001, respectively. The abbreviations are as follows: yield (Y), plant height (PH), spike length (SL), spike number per 666.7 m2 (SN), kernel number per spike (KNPS), thousand kernel weight (TKW), harvest index (HI), concentrations of Zn, Fe, Mn, Cu, N, P, K, Ca, Mg and phytate-P, ratios of phytate-P/P, and molar ratios of phytic acid (PA)/Zn, PA × Ca/Zn, PA/Fe and PA × Ca/Fe in wheat grains.
Figure 3.
A schematic diagram showing integrative strategies for simultaneously achieving yield increase and wheat grain Zn biofortification.
Figure 3.
A schematic diagram showing integrative strategies for simultaneously achieving yield increase and wheat grain Zn biofortification.
Table 1.
Detailed site information including geographic coordinates and soil basal properties in 0-20 cm soil layers before wheat sowing.
Table 1.
Detailed site information including geographic coordinates and soil basal properties in 0-20 cm soil layers before wheat sowing.
Experimental year |
Site |
Geographic coordinates |
Soil type |
pH (2.5:1 Water: soil ratio) |
Organic matter (g/kg) |
Total nitrogen (g/kg) |
Olsen-P (mg/kg) |
Exchangeable K (mg/kg) |
DTPA-extractable Zn (mg/kg) |
2020-2021 |
Jiyang |
116°58'E, 36°58'N |
Calcareous alluvial soil |
8.1 |
12.7 |
0.94 |
23.1 |
103.5 |
1.7 |
2022-2023 |
Liuyuan |
115°07'E, 35°02'N |
Calcareous alluvial soil |
8.3 |
13.5 |
0.82 |
15.3 |
98.7 |
1.2 |
Table 2.
Treatments of drone-based foliar spraying.
Table 2.
Treatments of drone-based foliar spraying.
Treatment |
Number of sprays/ spraying date |
Deionized water (L/ha) |
ZnSO4·7H2O (kg/ha) |
ZnO (kg/ha) |
CK 1
|
First/2023-05-08 |
150 L/ha |
0 |
0 |
ZnO |
150 L/ha |
0 |
3 kg/ha |
Zn |
150 L/ha |
3 kg/ha |
0 |
CK |
Second/2023-05-17 |
150 L/ha |
0 |
0 |
ZnO |
150 L/ha |
0 |
3 kg/ha |
Zn |
150 L/ha |
3 kg/ha |
0 |
CK |
Third/2023-05-23 |
150 L/ha |
0 |
0 |
ZnO |
150 L/ha |
0 |
3 kg/ha |
Zn |
150 L/ha |
3 kg/ha |
0 |
Table 3.
Grain yields, yield components and other agronomic traits of wheat as affected by different soil fertilization treatments and different drone-based foliar spraying treatments.
Table 3.
Grain yields, yield components and other agronomic traits of wheat as affected by different soil fertilization treatments and different drone-based foliar spraying treatments.
Experiment
|
Treatment
|
Plant height (cm) |
Spike length (cm) |
Spike number (10000/666.7 m2) |
KNPS
|
TKW (g) |
Grain yield (t/ha) |
HI (%) |
Soil Fertilization
|
N-P2O5-K2O (15-15-15) |
72.4a 1
|
6.5c |
34.4a |
34.5ab |
48.2a |
7.1b |
63.2a |
15-15-15+Micronutrient |
74.5a |
6.7bc |
36.5a |
32.3b |
48.2a |
7.5a |
61.1ab |
17-17-17 |
72.3a |
6.9abc |
39.2a |
36.5a |
47.3a |
7.6a |
62.0ab |
26-10-15 |
75.5a |
7.1ab |
35.7a |
33.4ab |
48.2a |
7.5a |
59.8b |
30-10-11 |
75.6a |
7.3a |
36.8a |
36.2ab |
46.7a |
7.4a |
61.5ab |
Foliar Spraying
|
CK |
77.7a |
8.3a |
39.8a |
36.8a |
43.9b |
7.5b |
59.4a |
ZnO |
84.5a |
8.4a |
41.4a |
37.5a |
45.7a |
8.6a |
61.5a |
Zn |
78.7a |
8.1a |
41.8a |
39.9a |
46.2a |
8.8a |
61.8a |
Table 4.
Grain nutrient concentration and acquisition of wheat as affected by different soil fertilization treatments.
Table 4.
Grain nutrient concentration and acquisition of wheat as affected by different soil fertilization treatments.
Parameter |
Treatment |
Zn |
Fe |
Mn |
Cu |
|
N |
P |
K |
Ca |
Mg |
Concentration |
|
mg/kg |
|
g/kg |
N-P2O5-K2O (15-15-15) |
19.9b 1
|
26.4a |
26.5a |
2.4ab |
|
19.2a |
3.1a |
4.5a |
0.33a |
1.39ab |
15-15-15+ Micronutrient |
23.5ab |
26.0ab |
25.2ab |
2.6ab |
|
18.6a |
3.0a |
4.7a |
0.34a |
1.41ab |
17-17-17 |
19.4b |
25.0ab |
26.1ab |
2.3b |
|
17.0a |
3.0a |
4.6a |
0.34a |
1.37b |
26-10-15 |
22.5ab |
22.2c |
23.9b |
2.4ab |
|
18.6a |
3.0a |
4.6a |
0.32a |
1.38b |
30-10-11 |
27.0a |
24.2b |
25.7ab |
2.8a |
|
19.9a |
3.2a |
4.8a |
0.33a |
1.45a |
Acquisition |
|
g/ha |
|
kg/ha |
N-P2O5-K2O (15-15-15) |
141.2b |
187.7ab |
188.1a |
17.4b |
|
136.9a |
12.6b |
32.2a |
2.3a |
9.9a |
15-15-15+ Micronutrient |
176.2ab |
195.0a |
189.0a |
19.3ab |
|
139.7a |
14.6ab |
35.0a |
2.5a |
10.6a |
17-17-17 |
144.5b |
189.1ab |
196.6a |
17.4b |
|
128.1a |
14.1ab |
34.5a |
2.6a |
10.3a |
26-10-15 |
168.9ab |
167.6b |
180.1a |
17.8ab |
|
140.4a |
15.2a |
34.3a |
2.4a |
10.4a |
30-10-11 |
197.6a |
178.3ab |
188.5a |
20.8a |
|
145.9a |
14.7ab |
35.2a |
2.5a |
10.7a |
Table 5.
Grain nutrient concentration and acquisition of wheat as affected by different drone-based foliar spraying treatments.
Table 5.
Grain nutrient concentration and acquisition of wheat as affected by different drone-based foliar spraying treatments.
Parameter |
Treatment |
Zn |
Fe |
Mn |
Cu |
|
P |
K |
Ca |
Mg |
Concentration |
|
mg/kg |
|
g/kg |
|
CK |
33.5c 1
|
25.9b |
13.3c |
4.5c |
|
3.5a |
4.3a |
0.392ab |
1.5a |
|
ZnO |
41.9b |
28.3a |
15.9a |
4.8b |
|
3.6a |
4.0a |
0.389b |
1.5a |
|
Zn |
43.6a |
29.0a |
14.4b |
5.4a |
|
3.7a |
4.2a |
0.430a |
1.5a |
Acquisition |
|
g/ha |
|
kg/ha |
|
CK |
252.5b |
195.0c |
100.0c |
33.7c |
|
26.5b |
32.1c |
3.0c |
11.0b |
|
ZnO |
360.8a |
243.2b |
136.7a |
41.3b |
|
31.0ab |
34.7ab |
3.3b |
12.8a |
|
Zn |
383.5a |
254.9a |
126.8b |
47.3a |
|
32.5a |
36.9a |
3.8a |
13.0a |
Table 6.
Concentrations of phytic acid (PA) and phytate-P, ratios of phytate-P/P, and molar ratios of PA/Zn, PA × Ca/Zn, PA/Fe and PA × Ca/Fe in wheat grains as affected by different soil fertilization treatments.
Table 6.
Concentrations of phytic acid (PA) and phytate-P, ratios of phytate-P/P, and molar ratios of PA/Zn, PA × Ca/Zn, PA/Fe and PA × Ca/Fe in wheat grains as affected by different soil fertilization treatments.
Treatment |
PA (g/kg) |
Phytate-P (g/kg) |
Phytate-P/P |
PA/Zn |
PA × Ca/Zn |
PA/Fe |
PA × Ca/Fe |
N-P2O5-K2O (15-15-15) |
8.1a 1
|
2.3a |
0.74a |
40.2ab |
327.1ab |
26.0c |
212.5b |
15-15-15+Micronutrient |
8.1a |
2.3a |
0.75a |
34.3b |
290.3b |
26.6bc |
225.9ab |
17-17-17 |
8.9a |
2.5a |
0.82a |
46.5a |
394.7a |
30.1abc |
257.5ab |
26-10-15 |
8.8a |
2.5a |
0.82a |
38.4ab |
307.6b |
33.4a |
268.6a |
30-10-11 |
8.8a |
2.5a |
0.79a |
33.5b |
278.5b |
31.0ab |
259.2ab |