# Cadmium and Zinc Concentration in Grain of Durum Wheat in Relation

Authors: Xiaopeng Gao; Cynthia A. Grant
Journal: Applied and Environmental Soil Science
 (2012)

---

## Abstract

Field experiments were conducted at two locations in Manitoba, Canada, to determine the effect of crop rotation, phosphorus (P) fertilization and tillage on grain yield and grain concentrations of Cd and Zn in durum wheat (Triticum durum L.). Compared to conventional tillage (CT), reduced tillage (RT) management decreased grain Cd and increased grain yield and grain Zn in half of the site-years. The type of preceding crops of spring wheat-flax or canola-flax had little influence. Rate and timing of P application had little effect on grain Cd, but increasing P rate tended to decrease grain Zn. No interactive effect was detected among tested factors. Grain Zn was not related to grain Cd, but positively to other nutrients such as Fe, Mn, P, Ca, K, and Mg. Both grain Zn and Fe correlated positively with grain protein content, suggesting protein may represent a sink for micronutrients. The study suggested that the tillage management may have beneficial effects on both grain yield and quality. Phosphorus fertilizer can remain available for subsequent crops and high annual inputs in the crop sequence may decrease crop grain Zn. Understanding the environment is important in determining the impact of agricultural management on agronomic and nutrient traits.


---

## Body

## 1. Introduction

Cadmium (Cd) accumulation in soils and cereal crops and its transfer to the human diet is a widespread problem around the world. Durum wheat (Triticum  durum L.) is of particular concern because it accumulates more Cd than the other commonly grown cereals with accumulation increasing in the order of rye < barley < oats < bread wheat < durum wheat [1]. Cadmium concentration in durum wheat grain harvested on Canadian prairies have been reported to range from less than 50 to more than 300 μg kg−1 [2], at times exceeding the 200 μg kg−1 limit set by the Codex Alimentarius Commission [3]. In addition, approximately 2.1 × 106 ha durum wheat, which occupies 10% of worldwide durum production area, is grown in the western Prairie region of Canada [4]. Therefore, there is a desire in the Canadian farming industry to control the Cd levels in the durum grain, either by improved agricultural management practices [5] or by breeding low Cd-accumulating cultivars [6, 7].Accumulation of metal elements in crop grains can be regulated by several physiological processes, including uptake from the soil solution, root-to-shoot translocation, and retranslocation into the grain during maturation. Zinc (Zn) and Cd are chemically similar and can compete for common transport mechanisms for uptake and translocation in the crop [5]. In contrast to Cd, Zn concentrations in grains of cereal crops such as wheat are often lower than desired level as sources of minerals for the human diet [8]. High consumption of cereal-based foods with low Zn concentration can result in Zn malnutrition in human. On the Canadian prairies, Zn concentration in wheat grain generally ranges between 20 and 30 mg kg−1 [9]. These concentrations are not adequate for human nutrition in diets with wheat constituting the main source of essential minerals [8]. Therefore, it is desirable to increase the Zn concentration in wheat grain to reduce the risk of Zn deficiency in the human diet.Phosphorus fertilizer is a major input for crop production. The fertilizer efficiency, however, is often quite low during the season of application because of the low mobility and rapid immobilization of phosphate after being applied to the soil. Application of P fertilizer can influence soil Cd availability and Cd accumulation in crops, either by direct addition of Cd as a contaminant in P fertilizers, or by indirect impact on soil properties, plant nutrition, and plant growth [5, 10]. For example, application of high-Cd P fertilizers led to increases in soil and plant Cd in many long-term field studies [11–13]. Phosphorus fertilizer is also frequently reported to reduce Zn concentration in crops and potential mechanisms include interactions between Zn and P in the soil [14], less translocation of Zn from the roots to the shoots [15], metabolic disorder within plant cells [16], and dilution effect [17]. In a short-term growth chamber study [18], the immediate increase of Cd accumulation in durum wheat following phosphate application was mostly due to increased Cd uptake and translocation by decreasing Zn accumulation in durum wheat induced by P fertilization, rather than to a direct addition effect. Similar results were also reported in other studies [19, 20], where using reagent-grade phosphate increased Cd concentration, but decreased Zn concentration of durum wheat. Though these studies showed the long-term or short-term effect of P fertilization on phytoavailability of Cd and Zn in soil-plant system, there is little information available on the effect of the residual P from P fertilizer that was applied to preceding crops.On the Canadian prairies, rotation of wheat with oilseeds or pulse crops is widely adopted and has been reported to increase grain yield and quality of subsequent crops through reduced disease incidence, improved water use efficiency, soil physical and chemical properties, and enhanced soil ecological environments [21–23]. Several studies have addressed the impact of crop rotation on Cd concentration in crops. Oliver et al. [24] reported a greater Cd concentration in spring wheat grown after lupin than after cereal and proposed the effect was due to increased Cd availability by rhizosphere acidification and root release of citric acid in lupin. In a 3-year field trial in Manitoba,Canada, seed Cd concentration in flax was lower when grown after bread wheat (Triticum aestivum L.) than after canola (Brassica napus L.), which was attributed to an increase in mycorrhizal colonization or to the lower Cd concentration in the decomposing wheat residue relative to canola residue [25].Reduced tillage practices are widely adopted by farmers on the Canadian prairies because of the benefits in conserving soil moisture, preventing soil erosion and reducing cost of production [26]. Difference in tillage management may have an impact on pH and nutrient stratification in the soil profile or residue decomposition [27], and consequently affect availability of Cd and Zn. The effect of tillage management on uptake of Cd and Zn in crops has been investigated in several field studies and the results are inconsistent. Tiller et al. [28] reported a 30% greater Cd concentration in spring wheat grain under direct drilling, compared to reduced tillage. Franzluebbers and Hons [29] reported that soil under a no-tillage system contained greater amounts of extractable Zn than under a conventional tillage system. In some other field studies, however, long-term tillage practices did not consistently influence Cd concentration in soils and crops [25, 30, 31].Information on the influence of agricultural management practices on crop yield and on Cd and Zn concentration and accumulation in crops is important in order to select management practices that optimize both crop productivity and quality. Although many studies have evaluated the impact of tillage system, crop sequence and fertilizer practices in the short term, effects of these management practices can persist and have impacts on subsequent crops. For example, P fertilizer not used by a crop in the year of production will remain in the soil and may exert residual effect on the following crop. Effects of tillage may increase over time as changes in soil organic matter and nutrient stratification intensify. However, information on the persistence of management effects on the yield and quality of later crops in a cropping system are limited.Based on field studies conducted at two locations over a four-year period, yield and concentrations of Cd and Zn of flax were shown to be influenced by tillage system, type of preceding crop and P application to the preceding crop [25, 32], but information is lacking about whether such effects can persist to influence subsequent crops. Therefore, to address this question, the field study described by Grant et al. [25, 32] was extended for an additional year to determine if preceding crop, P fertilization and tillage system would have a continuing effect on the yield and quality of the durum wheat grown in the third year of the cropping sequence.


## 2. Materials and Methods

### 2.1. Experimental Treatment and Cultural Practices

Field plot experiments were conducted at two different locations near Brandon, Manitoba. Both sites were classified as clay loam, Orthic Black Chernozemic soils. Initial characteristics of the soils are provided in Table1 [25]. Sampling and determination methods for soil characteristics were also described by Grant et al. [25]. One site (MZTRF) had been under reduced tillage (RT) management for 6 years before the study was initiated, while the second site (BRC North) had been under conventional tillage (CT) management until the establishment of the study. Neither of the soils was industrially contaminated and as such both are representative of soils commonly used for durum wheat production on the Canadian prairies.Table 1
Soil characteristics of the experimental sites prior to initiation of the study.
Soil propertySoil depth (cm)BRC NorthMZTRF2001*20022003200120022003pH0–157.828.377.838.437.877.5415–30——8.17——7.7130–60——8.43——7.95EC (μS cm−1)0–1529028960326129358515–3021829053722822849530–60299419560226276523N (mg kg−1)0–158.343.782.507.924.197.5015–303.312.562.753.422.033.1930–602.413.561.172.921.811.81P (mg kg−1)0–1512.5311.636.759.3611.349.94K (mg kg−1)0–152672022599.36298290DTPA-Cd (μg kg−1)0–15102.673.372.089.6134.8103.6DTPA-Zn (mg kg−1)0–150.770.600.801.011.081.20*Soil samples were taken in the first phase of the study, prior to seeding the preceding canola and wheat.The study was conducted over a 5 yr period from 1999–2003, with 3 yr cropping cycles being conducted at each of two locations in 1999-2000-2001, 2000-2001-2002, and 2001-2002-2003.In the first phase of the 3 yr cycle, either spring wheat (cv. AC Barrie) or canola (cv. G3295) was seeded under CT or RT, with P fertilizer treatments of 0, 11 and 22 kg P ha−1 as monoammonium phosphate (MAP). In the second phase, each plot where wheat or canola had been grown in the preceding season was divided into two subplots. Flax (cv. AC Emerson) was seeded into each subplot, with P fertilizer treatment of 0 or 11 kg P ha−1. Tillage treatment was kept the same as the preceding season. For the first and second phases, spring wheat, canola, and flax were harvested at maturity with a plot combine and the crop residue was chopped and returned to the plot where it had originated. Details on the cultural practices and sampling are available by Grant et al. [25, 32].In the third phase, each subplot where flax had been grown in the preceding season was seeded with durum wheat (cv. AC Avonlea). Ammonium nitrate was banded in all subplots prior to seeding at a rate of 100 kg N ha−1, adjusted for the amount of N added in the MAP. Each subplot where no P was applied in the preceding season received 11 kg P ha−1 side-banded as MAP, whereas the subplots where 11 kg P ha−1 was applied in the preceding season received no P. This provided a comparison between the effects of P applied in the current season and the residual effects of P applied in the previous year. Again, tillage treatment was continued as in the preceding season. Therefore, as shown in Table 2, the experimental design was comprised of two tillage system, two crop rotation, and six P fertilizer treatments, for a total of 24 combinations. The 3 yr cycle of spring wheat-flax-durum wheat or canola-flax-durum wheat was repeated three times at each site. This paper reports data on durum wheat over three growing seasons (2001–2003) on two sites, for a total of 6 site-years.Table 2
Description of the experimental design in the study.
TreatmentsPhases1st2nd3rdTillageCTCTCTZTZTZTCrop rotationSpring wheatFlaxDurum wheatCanolaFlaxDurum wheatP fertilization (kg MAP ha−1)0025025025025252505002550250The treatments, with four replicates, were arranged as a randomized complete block, with a split-plot layout for a total of 96 subplots at each location. Main blocks were the two tillage treatments and subplots were 12 combinations of crop rotation and P fertilizer treatments. Tillage blocks were initiated in the first year of the study and continued until the completion of the study. The RT consisted of only fertilizer banding and seeding operations and the CT received two tillage passes in the fall and one pass in the spring with a cultivator equipped with tine harrows. The preceding crops were established on a different area within the tillage blocks each year. Registered herbicides were applied as required and according to recommendations of the Manitoba Guild to Crop Protection.


### 2.2. Sampling and Analysis

Prior to seeding of the durum wheat, 0–15 cm depth soil core samples were collected in each plot. Soil samples were air-dried and sieved. Available P was extracted from the soil with 0.5 mol L−1 NaHCO3 [33]. Cadmium and Zn were extracted using DTPA [34]. Soil pH and EC were determined by a glass electrode using a 1 : 2 soil : water ratio. Concentrations of P and Zn in the extract were measured with an ARL 3410 ICP unit, and Cd on a Varian 300/400 atomic absorption spectrophotometer at a wavelength of 228.8 nm using a graphite furnace with deuterium correction (detection limit 0.01 μg Cd L−1). Reliability of the analysis was assessed by including certified soil reference materials and repeated samples in each set of soil extracts. Measured concentrations matched the stated ranges in the standards.At crop maturity, grain yield of durum wheat was measured by harvesting the center four rows of the plot using a plot combine. Grain was dried at 30°C to a constant moisture level and the weight recorded for each plot. Grain samples were ground and digested with a boiling acid mixture (HNO3 + HClO4) for Cd and Zn [35]. The concentrations of Zn in the digest were determined on an ARL 3520 inductively coupled plasma (ICP) and Cd concentrations were determined on a Varian 300/400 atomic absorption spectrophotometer (AAS) at a wavelength of 228.8 nm using a graphite furnace with deuterium correction (detection limit 0.01 μg Cd L−1). The reliability of the analysis was assessed by including certified plant reference materials and duplicate samples in each set of digests.


### 2.3. Statistical Analysis

Analysis of variance (Proc Mixed of SAS) was conducted on data separated by site-year. All interactions among fixed effect variables (tillage, crop rotation and P fertilization) were tested in the models and results showed that all interactions for grain yield and grain concentrations of elements were not significant. Therefore, only main effects of treatments were presented in this study. Means of the treatments were compared using protected least significance difference (LSD) at the 5% level of probability. Relationships between parameters were determined using linear regression analysis. Principal component analysis (PCA) was used to describe relationships among variables and was performed using the correlation matrix method. Loading plots were generated using principal components 1 and 2 as axes, and variables were plotted along the axes. PCA was conducted with Minitab Statistical Software (Minitab 15) and other analyses were performed with SAS program, release 9.1 (SAS Institute Inc., Cary, USA).


## 2.1. Experimental Treatment and Cultural Practices

Field plot experiments were conducted at two different locations near Brandon, Manitoba. Both sites were classified as clay loam, Orthic Black Chernozemic soils. Initial characteristics of the soils are provided in Table1 [25]. Sampling and determination methods for soil characteristics were also described by Grant et al. [25]. One site (MZTRF) had been under reduced tillage (RT) management for 6 years before the study was initiated, while the second site (BRC North) had been under conventional tillage (CT) management until the establishment of the study. Neither of the soils was industrially contaminated and as such both are representative of soils commonly used for durum wheat production on the Canadian prairies.Table 1
Soil characteristics of the experimental sites prior to initiation of the study.
Soil propertySoil depth (cm)BRC NorthMZTRF2001*20022003200120022003pH0–157.828.377.838.437.877.5415–30——8.17——7.7130–60——8.43——7.95EC (μS cm−1)0–1529028960326129358515–3021829053722822849530–60299419560226276523N (mg kg−1)0–158.343.782.507.924.197.5015–303.312.562.753.422.033.1930–602.413.561.172.921.811.81P (mg kg−1)0–1512.5311.636.759.3611.349.94K (mg kg−1)0–152672022599.36298290DTPA-Cd (μg kg−1)0–15102.673.372.089.6134.8103.6DTPA-Zn (mg kg−1)0–150.770.600.801.011.081.20*Soil samples were taken in the first phase of the study, prior to seeding the preceding canola and wheat.The study was conducted over a 5 yr period from 1999–2003, with 3 yr cropping cycles being conducted at each of two locations in 1999-2000-2001, 2000-2001-2002, and 2001-2002-2003.In the first phase of the 3 yr cycle, either spring wheat (cv. AC Barrie) or canola (cv. G3295) was seeded under CT or RT, with P fertilizer treatments of 0, 11 and 22 kg P ha−1 as monoammonium phosphate (MAP). In the second phase, each plot where wheat or canola had been grown in the preceding season was divided into two subplots. Flax (cv. AC Emerson) was seeded into each subplot, with P fertilizer treatment of 0 or 11 kg P ha−1. Tillage treatment was kept the same as the preceding season. For the first and second phases, spring wheat, canola, and flax were harvested at maturity with a plot combine and the crop residue was chopped and returned to the plot where it had originated. Details on the cultural practices and sampling are available by Grant et al. [25, 32].In the third phase, each subplot where flax had been grown in the preceding season was seeded with durum wheat (cv. AC Avonlea). Ammonium nitrate was banded in all subplots prior to seeding at a rate of 100 kg N ha−1, adjusted for the amount of N added in the MAP. Each subplot where no P was applied in the preceding season received 11 kg P ha−1 side-banded as MAP, whereas the subplots where 11 kg P ha−1 was applied in the preceding season received no P. This provided a comparison between the effects of P applied in the current season and the residual effects of P applied in the previous year. Again, tillage treatment was continued as in the preceding season. Therefore, as shown in Table 2, the experimental design was comprised of two tillage system, two crop rotation, and six P fertilizer treatments, for a total of 24 combinations. The 3 yr cycle of spring wheat-flax-durum wheat or canola-flax-durum wheat was repeated three times at each site. This paper reports data on durum wheat over three growing seasons (2001–2003) on two sites, for a total of 6 site-years.Table 2
Description of the experimental design in the study.
TreatmentsPhases1st2nd3rdTillageCTCTCTZTZTZTCrop rotationSpring wheatFlaxDurum wheatCanolaFlaxDurum wheatP fertilization (kg MAP ha−1)0025025025025252505002550250The treatments, with four replicates, were arranged as a randomized complete block, with a split-plot layout for a total of 96 subplots at each location. Main blocks were the two tillage treatments and subplots were 12 combinations of crop rotation and P fertilizer treatments. Tillage blocks were initiated in the first year of the study and continued until the completion of the study. The RT consisted of only fertilizer banding and seeding operations and the CT received two tillage passes in the fall and one pass in the spring with a cultivator equipped with tine harrows. The preceding crops were established on a different area within the tillage blocks each year. Registered herbicides were applied as required and according to recommendations of the Manitoba Guild to Crop Protection.


## 2.2. Sampling and Analysis

Prior to seeding of the durum wheat, 0–15 cm depth soil core samples were collected in each plot. Soil samples were air-dried and sieved. Available P was extracted from the soil with 0.5 mol L−1 NaHCO3 [33]. Cadmium and Zn were extracted using DTPA [34]. Soil pH and EC were determined by a glass electrode using a 1 : 2 soil : water ratio. Concentrations of P and Zn in the extract were measured with an ARL 3410 ICP unit, and Cd on a Varian 300/400 atomic absorption spectrophotometer at a wavelength of 228.8 nm using a graphite furnace with deuterium correction (detection limit 0.01 μg Cd L−1). Reliability of the analysis was assessed by including certified soil reference materials and repeated samples in each set of soil extracts. Measured concentrations matched the stated ranges in the standards.At crop maturity, grain yield of durum wheat was measured by harvesting the center four rows of the plot using a plot combine. Grain was dried at 30°C to a constant moisture level and the weight recorded for each plot. Grain samples were ground and digested with a boiling acid mixture (HNO3 + HClO4) for Cd and Zn [35]. The concentrations of Zn in the digest were determined on an ARL 3520 inductively coupled plasma (ICP) and Cd concentrations were determined on a Varian 300/400 atomic absorption spectrophotometer (AAS) at a wavelength of 228.8 nm using a graphite furnace with deuterium correction (detection limit 0.01 μg Cd L−1). The reliability of the analysis was assessed by including certified plant reference materials and duplicate samples in each set of digests.


## 2.3. Statistical Analysis

Analysis of variance (Proc Mixed of SAS) was conducted on data separated by site-year. All interactions among fixed effect variables (tillage, crop rotation and P fertilization) were tested in the models and results showed that all interactions for grain yield and grain concentrations of elements were not significant. Therefore, only main effects of treatments were presented in this study. Means of the treatments were compared using protected least significance difference (LSD) at the 5% level of probability. Relationships between parameters were determined using linear regression analysis. Principal component analysis (PCA) was used to describe relationships among variables and was performed using the correlation matrix method. Loading plots were generated using principal components 1 and 2 as axes, and variables were plotted along the axes. PCA was conducted with Minitab Statistical Software (Minitab 15) and other analyses were performed with SAS program, release 9.1 (SAS Institute Inc., Cary, USA).


## 3. Results

### 3.1. Grain Yield

Grain yield was greater with RT than CT in 3 site-years and less in 1 site-year (Table3). Crop rotation and P fertilization, however, had little effect on grain yield. Grain yield was greater when canola-flax rather than wheat-flax were the preceding crops in only one of 6 site-years, being BRC-North site in 2002. Differences in grain yield were only observed among P treatments at MZTRF in 2002. There was no significant interaction among the treatments in their effect on grain yield.Table 3
Grain yield (kg ha−1) of durum wheat as affected by tillage system, crop rotation, and P fertilization on at two sites over 3 years.
BRC NorthMZTRF200120022003200120022003TillageConventional179030862696328232212912Reduced188828353042323535662870LSD7113726979139101Crop rotationwheat-flax-durum182328512888323433932933Canola-flax-durum185430702851328333942849LSD7013726979139101P fertilization0-0-251894299329273291330029000-25-017332801267732943131289525-0-2518392919286032173552286925-25-018323043290031423284283750-0-2518712931293932653546294550-25-0186530772911334335462901LSD123237466137241174Analysis of varianceSourceDFPr ≥ FTillage10.00790.00050.0127n.s.<0.0001n.s.Crop rotation1n.s.0.0022n.s.n.s.n.s.n.s.P fertilization5n.s.n.s.n.s.n.s.0.0015n.s.


### 3.2. Grain Cd Concentration and Accumulation

Grain Cd concentration ranged from 66.6 to 150.4μg kg−1 dry weight (56.6 to 127.5 μg kg−1 fresh weight based on approximately 15% water content). Grain Cd was less with RT than CT in 3 site-years and greater in 1 site-year (Table 4). Accumulation of Cd in the grain was also less with RT in 2 of 6 site-years and greater in 1 site-year. Effects of crop rotation and P fertilizer on grain Cd concentration and accumulation were insignificant in all but one site-year, being MZTRF-2001, where grain Cd concentration and accumulation was greater when following canola-flax than wheat-flax. Accumulation of Cd was increased by P fertilizer only at MZTRF in 2002, and was not affected in other site-years. There was no significant interaction among the treatments in their effect on Cd concentration and accumulation in the grain.Table 4
Grain Cd concentration (μg kg−1) and accumulation (mg ha−1) of durum wheat as affected by tillage system, crop rotation, and P fertilization on at two sites over 3 years.
Grain Cd concentration (μg kg−1)Grain Cd accumulation (mg ha−1)BRC NorthMZTRFBRC NorthMZTRF200120022003200120022003200120022003200120022003TillageConventional120.1122.8144.282.586.183.5216380401269279247Reduced120.470.1119.881.861.399.8228205380267221289LSD9.113.915.94.88.711.3204880183137Crop rotationwheat-flax-durum117.392.0131.178.670.590.8215270394255239270Canola-flax-durum123.3100.9132.985.676.992.5229316387282260266LSD9.113.915.94.88.711.3204880183137P fertilization0-0-25117.189.9132.283.168.488.62232724002742262580-25-0110.397.5116.280.372.987.919128731826522625725-0-25116.994.8133.982.381.894.621628339326529127725-25-0120.798.9126.078.568.383.122330736924822623950-0-25131.199.7150.487.184.3111.124430445928529532950-25-0125.597.8133.281.566.684.6234302404272237247LSD15.724.127.58.315.119.53584139325464Analysis of varianceSourceDFPr ≥ FTillage1n.s.<0.0010.003n.s.<0.0010.005n.s.<0.001n.s.n.s.0.0040.025Crop rotation1n.s.n.s.n.s.0.005n.s.n.s.n.s.n.s.n.s.0.005n.s.n.s.P fertilization5n.s.n.s.n.s.n.s.n.s.n.s.n.s.n.s.n.s.n.s.0.014n.s.


### 3.3. Grain Zn Concentration and Accumulation

Grain Zn concentration ranged from 23.0 to 42.3 mg kg−1 dry weight. Grain Zn concentration was higher with RT than CT in 3 of 6 site-years (Table 5). Compared to CT, RT increased grain Zn accumulation in 2001 and 2003, but decreased it in 2002 on BRC-North site. In contrast, grain Zn accumulation was not affected by tillage management on MZTRF site. Grain Zn concentration and accumulation in durum wheat was higher when following wheat-flax than canola-flax only at MZTRF site in 2003, and was not affected in other site-years. Application of P fertilizer to either the previous crop or the current durum wheat crop did not affect grain Zn concentration. With increased application rate, however, grain Zn concentration showed a significantly decreasing trend in 2 of 6 site-years. Accumulation of Zn was not affected by P management in any site-years. Similar to Cd, grain Zn concentration and accumulation was not affected by any interactions among the treatments.Table 5
Grain Zn concentration (mg kg−1) and accumulation (g ha−1) of durum wheat as affected by tillage system, crop rotation, and P fertilization on at two sites over 3 years.
Grain Zn concentration (mg kg−1)Grain Zn accumulation (g ha−1)BRC-NorthMZTRFBRC-NorthMZTRF200120022003200120022003200120022003200120022003TillageConventional38.026.623.735.637.025.368.282.465.9116.9119.374.6Reduced40.126.026.337.435.624.975.574.282.9120.9126.672.1LSD1.91.32.51.82.31.54.76.513.95.98.96.4Crop rotationwheat-flax-durum39.526.225.736.936.226.172.074.977.1119.6123.177.3Canola-flax-durum38.726.524.336.136.424.171.781.771.7118.3122.869.3LSD1.91.32.51.82.31.54.76.513.95.98.96.4P fertilization0-0-2537.927.424.938.035.725.671.982.376.1125.1118.174.80-25-042.327.426.537.438.727.473.277.673.4123.0121.080.625-0-2539.825.623.737.134.923.872.874.771.1119.4124.169.325-25-039.426.926.735.737.026.172.282.979.3112.0122.574.750-0-2537.925.023.435.734.323.071.473.571.8116.5121.168.450-25-037.325.624.835.237.024.669.578.874.8117.6130.772.1LSD3.32.24.33.14.12.68.211.324.110.315.411.2Analysis of varianceSourceDFPr ≥ FTillage10.036n.s.0.0340.043n.s.n.s.0.0030.0140.018n.s.n.s.n.s.Crop rotation1n.s.n.s.n.s.n.s.n.s.0.013n.s.0.042n.s.n.s.n.s.0.015P fertilization50.034n.s.n.s.n.s.n.s.0.015n.s.n.s.n.s.n.s.n.s.n.s.


### 3.4. Correlation Analysis

Correlation analysis (Table6) showed that grain Cd was positively correlated with grain Fe, P, K, Mg, and negatively correlated with grain Mn, Cu, and Ca. A slightly negative but significant relationship was observed between yield and grain Cd (r=-0.135, P<0.001). Both grain Fe and Zn were highly significantly positively correlated (P<0.001) with other elements including Mn, P, Ca, K, and Mg and with protein content. Grain Fe, but not Zn, was negatively correlated with grain yield and positively correlated with grain Cd. Grain Fe and Zn were also positively correlated. Grain P was positively correlated with other measured element concentrations and negatively correlated with grain yield. The full dataset used for correlation matrix analysis is presented in Supplementary Table 1 (available online at doi:10.1155/2012/817107).Table 6
Correlation coefficients (r) of grain concentrations of Cd, Zn, Fe and P to other elements concentrations, grain yield, grain protein concentration, and soil characteristics (The full dataset is presented in Supplementary Table 1).
CdZnFePGrainZn−0.022—0.507***0.688***Fe0.294***0.507***—0.598***Mn−0.150***0.239***0.441***0.504***Cu−0.134***0.096*0.0670.101*P0.223***0.688***0.598***—Ca−0.365***0.408***−0.0170.168***K0.081*0.755***0.579***0.833***Mg0.191***0.491***0.390***0.691***Protein0.0360.425***0.529***0.636***Yield−0.135***−0.014−0.302***−0.292***SoilpH−0.080−0.635***−0.193**−0.255***EC−0.094−0.047−0.041−0.095Olsen-P0.309***0.324***0.537***0.553***DTPA-Cd0.339***0.213***0.235***0.012DTPA-Zn0.0440.523***0.193***0.201***, **, *** indicate significant at 0.05, 0.01, and 0.001, respectively.Grain concentrations of Cd, Zn and P across site-years correlated significantly and positively with their extractable concentrations in soils (Table6). Durum grain Cd was significantly related to soil extractable Cd concentration (r=0.339, P<0.001), as well as to soil extractable P concentration (r=0.309, P<0.001). Grain concentrations of Zn, Fe, and P, but not Cd, were negatively related with soil pH.PCA of the dataset across 6 site-years extracted two major principal components, explaining 78% of the total variation in data (Figure1). PC1 explained 54.4% of the variation, and was loaded positively with grain Zn, Mn, P, K, Mg, Fe, and soil P, and negatively with grain yield and soil Cd. PC2 explained 23.6% of the variation, and was loaded positively with grain Ca and negatively with grain Cd and Cu. The PCA showed strong associations between grain protein and grain concentrations of minerals such as Zn, Mn, P, K, Mg, and Fe. In the biplot, the long and outside-scattered distribution of the six site-year vectors indicated that the environments were likely to have a great influence on the measured parameters. Grain Cd concentration was greatest in BRC-2003 and least in MZTRF-2002. In contrast, grain Zn concentration on either site was generally greater in 2001 than in other years. The least grain yield was recorded in BRC-2001 even though the soil had a relatively high P availability.Figure 1
Principal component analysis (based on correlation matrix) of grain yield, grain concentrations of Cd, Zn, Cu, Fe, Mg, K, P, Mn, Ca, grain protein concentration, and selected soil properties (soil P and soil Cd) across 6 site-years. Biplot vectors are trait factor loadings for principal component 1 (PC1) and PC2. Data at each site-year are averaged across P fertilizer treatments.


## 3.1. Grain Yield

Grain yield was greater with RT than CT in 3 site-years and less in 1 site-year (Table3). Crop rotation and P fertilization, however, had little effect on grain yield. Grain yield was greater when canola-flax rather than wheat-flax were the preceding crops in only one of 6 site-years, being BRC-North site in 2002. Differences in grain yield were only observed among P treatments at MZTRF in 2002. There was no significant interaction among the treatments in their effect on grain yield.Table 3
Grain yield (kg ha−1) of durum wheat as affected by tillage system, crop rotation, and P fertilization on at two sites over 3 years.
BRC NorthMZTRF200120022003200120022003TillageConventional179030862696328232212912Reduced188828353042323535662870LSD7113726979139101Crop rotationwheat-flax-durum182328512888323433932933Canola-flax-durum185430702851328333942849LSD7013726979139101P fertilization0-0-251894299329273291330029000-25-017332801267732943131289525-0-2518392919286032173552286925-25-018323043290031423284283750-0-2518712931293932653546294550-25-0186530772911334335462901LSD123237466137241174Analysis of varianceSourceDFPr ≥ FTillage10.00790.00050.0127n.s.<0.0001n.s.Crop rotation1n.s.0.0022n.s.n.s.n.s.n.s.P fertilization5n.s.n.s.n.s.n.s.0.0015n.s.


## 3.2. Grain Cd Concentration and Accumulation

Grain Cd concentration ranged from 66.6 to 150.4μg kg−1 dry weight (56.6 to 127.5 μg kg−1 fresh weight based on approximately 15% water content). Grain Cd was less with RT than CT in 3 site-years and greater in 1 site-year (Table 4). Accumulation of Cd in the grain was also less with RT in 2 of 6 site-years and greater in 1 site-year. Effects of crop rotation and P fertilizer on grain Cd concentration and accumulation were insignificant in all but one site-year, being MZTRF-2001, where grain Cd concentration and accumulation was greater when following canola-flax than wheat-flax. Accumulation of Cd was increased by P fertilizer only at MZTRF in 2002, and was not affected in other site-years. There was no significant interaction among the treatments in their effect on Cd concentration and accumulation in the grain.Table 4
Grain Cd concentration (μg kg−1) and accumulation (mg ha−1) of durum wheat as affected by tillage system, crop rotation, and P fertilization on at two sites over 3 years.
Grain Cd concentration (μg kg−1)Grain Cd accumulation (mg ha−1)BRC NorthMZTRFBRC NorthMZTRF200120022003200120022003200120022003200120022003TillageConventional120.1122.8144.282.586.183.5216380401269279247Reduced120.470.1119.881.861.399.8228205380267221289LSD9.113.915.94.88.711.3204880183137Crop rotationwheat-flax-durum117.392.0131.178.670.590.8215270394255239270Canola-flax-durum123.3100.9132.985.676.992.5229316387282260266LSD9.113.915.94.88.711.3204880183137P fertilization0-0-25117.189.9132.283.168.488.62232724002742262580-25-0110.397.5116.280.372.987.919128731826522625725-0-25116.994.8133.982.381.894.621628339326529127725-25-0120.798.9126.078.568.383.122330736924822623950-0-25131.199.7150.487.184.3111.124430445928529532950-25-0125.597.8133.281.566.684.6234302404272237247LSD15.724.127.58.315.119.53584139325464Analysis of varianceSourceDFPr ≥ FTillage1n.s.<0.0010.003n.s.<0.0010.005n.s.<0.001n.s.n.s.0.0040.025Crop rotation1n.s.n.s.n.s.0.005n.s.n.s.n.s.n.s.n.s.0.005n.s.n.s.P fertilization5n.s.n.s.n.s.n.s.n.s.n.s.n.s.n.s.n.s.n.s.0.014n.s.


## 3.3. Grain Zn Concentration and Accumulation

Grain Zn concentration ranged from 23.0 to 42.3 mg kg−1 dry weight. Grain Zn concentration was higher with RT than CT in 3 of 6 site-years (Table 5). Compared to CT, RT increased grain Zn accumulation in 2001 and 2003, but decreased it in 2002 on BRC-North site. In contrast, grain Zn accumulation was not affected by tillage management on MZTRF site. Grain Zn concentration and accumulation in durum wheat was higher when following wheat-flax than canola-flax only at MZTRF site in 2003, and was not affected in other site-years. Application of P fertilizer to either the previous crop or the current durum wheat crop did not affect grain Zn concentration. With increased application rate, however, grain Zn concentration showed a significantly decreasing trend in 2 of 6 site-years. Accumulation of Zn was not affected by P management in any site-years. Similar to Cd, grain Zn concentration and accumulation was not affected by any interactions among the treatments.Table 5
Grain Zn concentration (mg kg−1) and accumulation (g ha−1) of durum wheat as affected by tillage system, crop rotation, and P fertilization on at two sites over 3 years.
Grain Zn concentration (mg kg−1)Grain Zn accumulation (g ha−1)BRC-NorthMZTRFBRC-NorthMZTRF200120022003200120022003200120022003200120022003TillageConventional38.026.623.735.637.025.368.282.465.9116.9119.374.6Reduced40.126.026.337.435.624.975.574.282.9120.9126.672.1LSD1.91.32.51.82.31.54.76.513.95.98.96.4Crop rotationwheat-flax-durum39.526.225.736.936.226.172.074.977.1119.6123.177.3Canola-flax-durum38.726.524.336.136.424.171.781.771.7118.3122.869.3LSD1.91.32.51.82.31.54.76.513.95.98.96.4P fertilization0-0-2537.927.424.938.035.725.671.982.376.1125.1118.174.80-25-042.327.426.537.438.727.473.277.673.4123.0121.080.625-0-2539.825.623.737.134.923.872.874.771.1119.4124.169.325-25-039.426.926.735.737.026.172.282.979.3112.0122.574.750-0-2537.925.023.435.734.323.071.473.571.8116.5121.168.450-25-037.325.624.835.237.024.669.578.874.8117.6130.772.1LSD3.32.24.33.14.12.68.211.324.110.315.411.2Analysis of varianceSourceDFPr ≥ FTillage10.036n.s.0.0340.043n.s.n.s.0.0030.0140.018n.s.n.s.n.s.Crop rotation1n.s.n.s.n.s.n.s.n.s.0.013n.s.0.042n.s.n.s.n.s.0.015P fertilization50.034n.s.n.s.n.s.n.s.0.015n.s.n.s.n.s.n.s.n.s.n.s.


## 3.4. Correlation Analysis

Correlation analysis (Table6) showed that grain Cd was positively correlated with grain Fe, P, K, Mg, and negatively correlated with grain Mn, Cu, and Ca. A slightly negative but significant relationship was observed between yield and grain Cd (r=-0.135, P<0.001). Both grain Fe and Zn were highly significantly positively correlated (P<0.001) with other elements including Mn, P, Ca, K, and Mg and with protein content. Grain Fe, but not Zn, was negatively correlated with grain yield and positively correlated with grain Cd. Grain Fe and Zn were also positively correlated. Grain P was positively correlated with other measured element concentrations and negatively correlated with grain yield. The full dataset used for correlation matrix analysis is presented in Supplementary Table 1 (available online at doi:10.1155/2012/817107).Table 6
Correlation coefficients (r) of grain concentrations of Cd, Zn, Fe and P to other elements concentrations, grain yield, grain protein concentration, and soil characteristics (The full dataset is presented in Supplementary Table 1).
CdZnFePGrainZn−0.022—0.507***0.688***Fe0.294***0.507***—0.598***Mn−0.150***0.239***0.441***0.504***Cu−0.134***0.096*0.0670.101*P0.223***0.688***0.598***—Ca−0.365***0.408***−0.0170.168***K0.081*0.755***0.579***0.833***Mg0.191***0.491***0.390***0.691***Protein0.0360.425***0.529***0.636***Yield−0.135***−0.014−0.302***−0.292***SoilpH−0.080−0.635***−0.193**−0.255***EC−0.094−0.047−0.041−0.095Olsen-P0.309***0.324***0.537***0.553***DTPA-Cd0.339***0.213***0.235***0.012DTPA-Zn0.0440.523***0.193***0.201***, **, *** indicate significant at 0.05, 0.01, and 0.001, respectively.Grain concentrations of Cd, Zn and P across site-years correlated significantly and positively with their extractable concentrations in soils (Table6). Durum grain Cd was significantly related to soil extractable Cd concentration (r=0.339, P<0.001), as well as to soil extractable P concentration (r=0.309, P<0.001). Grain concentrations of Zn, Fe, and P, but not Cd, were negatively related with soil pH.PCA of the dataset across 6 site-years extracted two major principal components, explaining 78% of the total variation in data (Figure1). PC1 explained 54.4% of the variation, and was loaded positively with grain Zn, Mn, P, K, Mg, Fe, and soil P, and negatively with grain yield and soil Cd. PC2 explained 23.6% of the variation, and was loaded positively with grain Ca and negatively with grain Cd and Cu. The PCA showed strong associations between grain protein and grain concentrations of minerals such as Zn, Mn, P, K, Mg, and Fe. In the biplot, the long and outside-scattered distribution of the six site-year vectors indicated that the environments were likely to have a great influence on the measured parameters. Grain Cd concentration was greatest in BRC-2003 and least in MZTRF-2002. In contrast, grain Zn concentration on either site was generally greater in 2001 than in other years. The least grain yield was recorded in BRC-2001 even though the soil had a relatively high P availability.Figure 1
Principal component analysis (based on correlation matrix) of grain yield, grain concentrations of Cd, Zn, Cu, Fe, Mg, K, P, Mn, Ca, grain protein concentration, and selected soil properties (soil P and soil Cd) across 6 site-years. Biplot vectors are trait factor loadings for principal component 1 (PC1) and PC2. Data at each site-year are averaged across P fertilizer treatments.


## 4. Discussion

Grain yield was relatively low at the BRC-North site in 2001 (Table3, Figure 1), due to adverse weather conditions, but was average for the area in the other site-years. The range of Cd concentrations in durum grain in this study was also within the expected range, being comparable to that reported in our previous studies on the Canadian prairies [31, 36], and well below the proposed Codex maximum limit in wheat grain [3]. This suggests the durum wheat production in the area is generally safe for human consumption. However, Zn concentrations in the grain were generally less than the biofortification target of 40–60 mg kg−1, which is required to counteract human Zn deficiency in many parts of the world [37]. Therefore, it is desirable to increase Zn concentration in durum grain, either by genetic improvement or agricultural practices.Of the factors evaluated, tillage management showed the greatest influence on grain yield, as well as on grain Cd and Zn concentrations and accumulations in durum wheat. Compared to CT, RT significantly increased grain yield and Zn, and decreased grain Cd in half of the site-years (Tables3, 4, 5). The decrease in Cd, however, was neither necessarily related to the increase in Zn, nor to the increase in grain yield, as suggested by the nonrelationship or weak relationship among these three variables (Table 6). It should be noted that the observed effects of tillage on grain Cd and Zn were not expected because our previous studies suggested that tillage system had very limited influence on either Cd [29] or Zn [38] concentrations in durum grain. Also, in the second phase of the current study where the flax was the tested crop, tillage management also did not exert any consistent effect on flaxseed concentrations of Cd and Zn [25]. It is unclear why there was an effect of tillage in this growing season but not in the previous seasons. Some studies suggested that tillage management can have an impact on metal availability in soils by its influence on soil properties such as pH, organic carbon content and cation exchange capacity [29, 39]. This assumption, however, is not supported in this study because soil concentrations of extractable Cd and Zn were not consistently affected by tillage management (data not present). Instead, the difference in mycorrhizal association could play a role. Compared to RT, annual soil disturbances produced by CT may reduce mycorrhizal colonization in crops [40] and thus influence root uptake of Cd [41] or Zn [42]. In addition, tillage effects may vary from year to year because of differences in weather conditions. Under prairie conditions, RT tends to have greater benefits under warmer, drier conditions where moisture conservation become important. Similarly, seasonal differences in weather may interact with tillage effects on both root distribution in different soil layers and rate of decomposition of crop residues, to influence metal uptake. While the impact of tillage is not consistent across different studies on the Canadian prairies [25, 31, 32, 36], the RT practices may be beneficial to farmers in many regions not only because of the benefits in conserving soil moisture, preventing soil erosion and reducing the cost of production [26], but also for potential improvements in grain yield and grain quality, as has been shown in the current study.There was no consistent difference in grain yield or grain concentration of Cd or Zn between durum wheat following spring wheat-flax as compared to canola-flax. In contrast, in the second phase of this study, flaxseed following canola had significantly higher Cd concentration and lower Zn concentration than following spring wheat [25]. These results suggest that the impact of spring wheat or canola as a preceding crop was significant only on the following crop, and did not persist to affect the third crop in the sequence. Potential mechanisms explaining preceding crop effect on flax were discussed in Grant et al. [25], including difference in Cd and Zn content in the crop residue which was returned to the soil after harvest, difference in mycorrhizal association between mycorrhizal spring wheat and nonmycorrhizal canola, as well as difference in root-induced rhizosphere chemical changes between canola and wheat.In spite of the low initial soil P concentration at the two locations, there was no statistically significant main effect on grain yield from the different amounts or timings of P fertilization through the three-year cropping sequence (Table3). However, in 4 of 6 site-years, there was a statistical or numerical advantage from 160 to 250 kg ha−1 of grain yield with application of the P to the durum wheat rather than in the preceding crop at the lowest rate of application (Table 3), a difference that was significant at the P<0.02–0.28 rate using contrast analysis (data not present). This differential was much lower or nonexistent when higher amounts of P were applied through the cropping sequence. It appears that if the level of input through the cropping sequence is low, there may be a small yield advantage to providing some P as a starter near the seed row in the current crop. The benefit of the use of starter P diminishes as the P input through the rotation decreases.In other field studies, that is, in Sweden [43] and in the Canadian prairies [13], Cd concentration in grain and seed of several crops consistently increased with increasing P application, which was attributed to either the direct addition of the Cd contained in the fertilizer or the indirect effect on soil chemical properties. The MAP used in this study, however, contained only 3.7 mg Cd kg−1 and application of 75 kg fertilizer over 3 growing seasons would add only 278 mg of Cd to the soil per hectare, making it unlikely to have a strong effect on Cd concentration in crop grain. Similarly, in studies in the USA, long-term application of fertilizers containing less than 5 mg Cd kg−1 did not increase Cd concentration in soil and crops [44]. In 2 of 6 site-years, P fertilizer treatment exerted a significant effect on grain Zn concentration, with the trend being that grain Zn decreased with increased application rate over the three growing seasons (Table 5). Similarly, in the second phase of this study, seed Zn concentration of flax decreased with application of P in half of the site-years [25]. The P-induced inhibition of Zn accumulation is frequently reported and may be related to the P-Zn interaction in soil [14] and depression in mycorrhizal association caused by increased P supply [45].The impacts of the agricultural management on grain yield and grain quality were highly unstable. Their performance varied across locations and years. For example, the reduced tillage significantly improved grain yield and grain quality as compared to the conventional tillage in half of the site-years. The effect, however, was absent or even opposite in other site-years. Also, increasing P fertilizer rate decreased grain Zn in 2 site-years, but not in the others. Therefore, the environmental conditions such as location or meteorological factors are important determinants for agronomic and nutrient traits. This assumption is also confirmed by the wide range of the site-year criteria in the PCA biplot (Figure1). Still, only one wheat genotype was tested in the current study, neglecting the potential genotype by environment interactions, as been shown in other studies concerning grain Cd [46] and grain Zn [47]. Hence, greater efforts taking care of genotype by environment interactions are needed to produce improved grain quality for human health.Results of the linear regression and PCA showed grain Zn concentration was correlated positively to other nutrients such as Fe, Mn, P, Ca, K, and Mg, but was not correlated to grain Cd (Table6, Figure 1). The relationships favor the possibility of producing durum grain with moderately high nutrients, while maintaining low concentration of Cd. Also, both grain Zn and Fe concentrations correlated positively with grain protein content. This finding is consistent with other field trials [48] and suggested a possible link between grain protein and the levels of the two micronutrients.


## 5. Conclusions

In summary, of the factors evaluated, tillage had the most consistent effect on grain concentrations of Cd and Zn. Compared to CT, RT increased grain yield, grain Zn, and decreased grain Cd in half of the site-years and should therefore be recommended in tested area. The effect of growing wheat rather than canola as a preceding crop on crop yield or trace element concentration did not persist to affect durum wheat grown as the third crop in the sequence. There was evidence of a small impact of starter P on grain yield, but only if the levels of P applied through the cropping sequence were low, indicating that P fertilizer applied in preceding years can remain plant-available, reducing the requirements for P input in following crops. Rate and timing of P application had little effect on grain concentration of Cd, but increasing P rate tended to decrease grain Zn concentration. The linear and multivariate regressions revealed that grain Zn was not related to grain Cd, but was positively correlated with other nutrients such as Fe, Mn, P, Ca, K, and Mg, suggesting the possibility of producing durum grain with moderately high nutrients, while maintaining low concentration of Cd. Grain protein may represent a sink for micronutrients such as Zn and Fe. No interactive effect was detected among tested factors. Results of the study suggest that tillage management can have persistent effects on both grain yield and quality, that impacts of preceding crop do not persist past the first subsequent season and that P fertilizer can remain available for subsequent crops, reducing the response to annual inputs if P levels are maintained at adequate levels throughout the cropping sequence.



---

*Source: 817107-2012-03-19.xml*