Variation in grain Zn concentration, and the grain ionome, in field-grown Indian wheat

Wheat is an important dietary source of zinc (Zn) and other mineral elements in many countries. Dietary Zn deficiency is widespread, especially in developing countries, and breeding (genetic biofortification) through the HarvestPlus programme has recently started to deliver new wheat varieties to help alleviate this problem in South Asia. To better understand the potential of wheat to alleviate dietary Zn deficiency, this study aimed to characterise the baseline effects of genotype (G), site (E), and genotype by site interactions (GxE) on grain Zn concentration under a wide range of soil conditions in India. Field experiments were conducted on a diverse panel of 36 Indian-adapted wheat genotypes, grown on a range of soil types (pH range 4.5–9.5), in 2013–14 (five sites) and 2014–15 (six sites). Grain samples were analysed using inductively coupled plasma-mass spectrometry (ICP-MS). The mean grain Zn concentration of the genotypes ranged from 24.9–34.8 mg kg-1, averaged across site and year. Genotype and site effects were associated with 10% and 6% of the overall variation in grain Zn concentration, respectively. Whilst G x E interaction effects were evident across the panel, some genotypes had consistent rankings between sites and years. Grain Zn concentration correlated positively with grain concentrations of iron (Fe), sulphur (S), and eight other elements, but did not correlate negatively with grain yield, i.e. no yield dilution was observed. Despite a relatively small contribution of genotype to the overall variation in grain Zn concentration, due to experiments being conducted across many contrasting sites and two years, our data are consistent with reports that biofortifying wheat through breeding is likely to be effective at scale given that some genotypes performed consistently across diverse soil types. Notably, all soils in this study were probably Zn deficient and interactions between wheat genotypes and soil Zn availability/management (e.g. the use of Zn-containing fertilisers) need to be better-understood to improve Zn supply in food systems.


Introduction
Zinc is an essential element for plants and people [1,2,3].Approximately 20% percent of the world's population is at risk of having insufficient Zn in their diets, based on food supply and composition, with much greater proportions at risk in sub-Saharan Africa and South Asia [4,5,6].Dietary Zn deficiency is linked to diarrhoea, respiratory infections and stunting (low height to weight ratios) in children [7,8].Wheat (Triticum aestivum L.) is the third most widely consumed cereal crop by weight, after rice and maize, providing ~20% of total global energy and protein supply [9].This figure rises to >50% of energy in some countries, for example, in India (www.dwr.in).Wheat is an important dietary source of Zn [10] and supplies ~50% Zn of daily needs in India [5] (http://harvestplus.org/where-we-work/india).Over the past 20 years, agronomic and breeding-based approaches have been developed to increase the Zn concentration of wheat grain [11,12].The most effective agronomic approaches to increasing wheat grain Zn concentrations involve foliar applications of ZnSO 4 solution in the presence of adequate nitrogen [13].This approach is likely to be cost-effective in terms of yield and health benefits but may require changes to current farming practices [14][15][16][17][18]. Breeding for increased grain Zn concentration is likely to be a more cost-effective approach in the longer term provided there is adequate plant-available Zn in the soil [16,[19][20][21].
The 'HarvestPlus' programme is developing and deploying new wheat varieties with higher grain Zn concentrations in South Asia [22].The HarvestPlus target is to enhance the grain Zn concentration in locally-adapted wheat varieties by 8-12 mg kg -1 , from a baseline of ~25 mg kg -1 , without compromising grain yield or other grain quality attributes [23,24].Through collaboration with Banaras Hindu University (BHU), Uttar Pradesh (UP), India, and several private Indian seed companies including Sai Seeds, Nirmal Seeds, and Astha Beej, three varieties with increased grain Zn concentration have been multiplied and released in Eastern UP (North Eastern Plain Zone, NEPZ): Akshai (BHU-3), Abhay (Zinc Shakthi), and BHU-6 [25].A high zinc wheat variety 'Zincol-2016' has been recommended for release and multiplication in Pakistan by the Pakistan Agriculture Research Council (PARC; http://www.harvestplus.org/node/1647).The wheat genotypes used in the present study have been selected to represent a wide diversity of crop types which span all of the agro-climatic wheat growing zones of India for good yield, yield components and quality traits.Most of these genotypes are used as reference lines within advanced variety trials for screening new wheat lines, including under saline/ acid soils [26,27].
It is likely that optimal grain Zn concentrations will ultimately be achieved using a combination of crop improvement and agronomy [28,15,16,29,30].However, to date, little work has reported baseline relative contributions of genotype (G), environment (E) and G Ã E interactions on wheat grain Zn concentration under widely differing field conditions.Therefore, the aim of this study was to quantify baseline variation in grain Zn concentrations of fieldgrown wheat, using a diverse panel of Indian-adapted wheat genotypes grown on a wide range of soil types but without addition of Zn fertilisers.Simultaneously, the grain concentration of a wider range of elements was determined (NB, the mineral composition of biological tissues is typically referred to as the 'ionome' [31]).

Wheat genotypes, field sites and management
A panel of 36 elite wheat genotypes comprising Triticum aestivum L. (n = 34) and Triticum durum (n = 2) was selected for study under field conditions.These genotypes represent a diverse genetic background with adaptations to a range of climatic and soil environments.All genotypes and experimental conditions are described in Khokhar et al. [26], including data for yield and yield components.In brief, the panel was grown at six field sites in major wheat growing areas of India, during rabi (winter) seasons 2013-14 and 2014-15.The six sites were, (1)

Determining plant-available soil Zn
Top-and sub-soil were sampled with a hand-auger during both growing seasons at all six sites.Five samples were taken at each experimental site at each depth from the vertices of a 'W-pattern'.A single composite sample of 1 kg was retained for each soil depth.To determine plantavailable soil Zn concentrations, air-dried soil samples were sieved to 2 mm and a 1.5 g subsample extracted with 10 mL of 0.01 M Ca(NO 3 ) 2 by shaking for 1 h on a rotary shaker.Soil suspensions were then centrifuged at 3300 rpm for 20 minutes and the supernatant liquid was filtered through 0.22 μm polyether sulfone (PES) syringe filters and acidified to a final concentration of 2% HNO 3 (Primar Plus™, Fisher).Multi-elemental analyses was conducted on (technical) triplicate samples using inductively coupled plasma-mass spectrometry (ICP-MS), as described for the grain samples below.A total of 30 and 36 soil samples were analysed in 2013-14 and 2014-15, respectively; at Kumarganj-reclaimed and Kumarganj-sodic sites, soil samples were taken only from one depth, 0-30 cm, in 2013-14.

Grain sampling and digestion
Grains of each plot from 5 sites (all sites except Pundibari) in 2013-14, and from plots on all 6 sites in 2014-15 were harvested at physiological maturity and sun-dried.A representative subsample of grain from each plot was carefully hand-cleaned by discarding broken grains and foreign material.In 2013-14, ~10 grains of each genotype were weighed and digested using a microwave system comprising a Multiwave 3000 platform with a 48-vessel MF50 rotor (Anton Paar Gmbh, Graz, Austria); digestion vessels were perfluoroalkoxy (PFA) tubes in polyethylethylketone (PEEK) pressure jackets (Anton Paar GmbH).Grains were digested in 2 mL 70% Trace Analysis Grade (TAG) HNO 3 , 1 mL Milli-Q water (18.2MΩ cm; Fisher Scientific UK Ltd, Loughborough, UK), and 1 mL H 2 O 2 .Prior to digestion, the samples were soaked in this solution for 16 h at room temperature.In 2014-15, grain samples were crushed inside a paper bag, and subsamples (~0.200 g DW) were then used without the soaking step.The microwave settings were: power = 1400 W, temperature = 140˚C, pressure = 2 MPa, time = 45 minutes.Two operational blanks were included in each digestion run.Duplicate samples of a certified reference material (CRM: Wheat flour SRM 1567b, NIST, Gaithersburg, MD, USA) were included approximately every fourth digestion run.A laboratory reference material (LRM), Paragon was used for each digestion run.Following digestion, each tube was made up to a final volume of 15 mL by adding 11 mL Milli-Q water, then transferred to a 25 mL universal tube (Sarstedt Ltd., Nu ¨mbrecht, Germany) and stored at room temperature.Samples were further diluted 1:5 with Milli-Q water into 13 ml tubes (Sarstedt Ltd.) prior to analysis.

Zinc and multi-element analysis
The grain and plant-available soil concentrations of Zn and 30 other elements were determined by ICP-MS (Thermo Fisher Scientific iCAPQ, Thermo Fisher Scientific, Bremen, Germany).The other elements were Ag, Al, As, B, Ba, Be, Ca, Cd, Cr, Co, Cs, Cu, Fe, K, Li, Mg, Mn, Mo, Na, Ni, P, Pb, Rb, S, Se, Sr, Ti, Tl, U, and V. Three operational ICP-MS modes were used: (i) a helium collision-cell (He-cell) with kinetic energy discrimination to remove polyatomic interferences, (ii) standard mode (STD) in which the collision cell was evacuated, and (iii) a hydrogen collision-reaction cell (H 2 -cell).Samples were introduced from an auto sampler incorporating an ASXpress™ rapid uptake module (Cetac ASX-520, Teledyne Technologies Inc., Omaha, NE, USA) through a PEEK nebuliser (Burgener Mira Mist, Mississauga, Burgener Research Inc., Canada).Internal standards were introduced to the sample stream on a separate line via the ASXpress unit and included Sc (20 μg L -1 ), Rh (10 μg L -1 ), Ge (10 μg L -1 ) and Ir (5 μg L -1 ) in 2% TAG HNO3 (Fisher Scientific UK Ltd).External multi-element calibration standards (Claritas-PPT grade CLMS-2; SPEX Certiprep Inc., Metuchen, NJ, USA) included Ag, Al, As, B, Ba, Be, Cd, Ca, Co, Cr, Cs, Cu, Fe, K, Li, Mg, Mn, Mo, Na, Ni, P, Pb, Rb, S, Se, Sr, Ti, Tl (semi-quant), U, V and Zn, in the range 0-100 μg L -1 (0, 20, 40, 100 μg L -1 ).A bespoke external multi-element calibration solution (PlasmaCAL, SCP Science, Courtaboeuf, France) was used to create Ca, K, Mg and Na standards in the range 0-30 mg L -1 .Boron, P and S calibration utilised in-house standard solutions (KH 2 PO 4 , K 2 SO 4 and H 3 BO 3 ).In-sample switching was used to measure B and P in STD mode, Se in H 2 -cell mode and all other elements in He-cell mode.Sample processing was undertaken using Qtegra™ software (Thermo Fisher Scientific) with external cross-calibration between pulse-counting and analogue detector modes when required.In total, 288 grain samples were analysed in 7 runs in 2013-14, and 432 grain samples in 12 runs in 2014-15, excluding blank, CRM and LRM samples.The Zn-specific recovery from CRMs was 96% in 2013-14 and 94% in 2014-15 compared with certified CRM values.

Data analysis of grain Zn concentration and the grain ionome
For each data-point, an element-specific operational blank concentration (mean of each ICP-MS run) was subtracted.Data were then multiplied by initial sample volume, divided by the initial dry mass of material, and converted to mg element kg -1 dry grain material.Elementspecific limits of detection (LODs) were reported as 3 times the standard deviation (SD) of the ten operational blank concentrations, assuming a notional starting dry weight of 0.4 g in 2013-14 and 0.2 g in 2014-15 (S1 Table ).Out of 31 mineral elements, 15 (Ag, Al, B, Be, Cd, Co, Cr, Li, Na, Ni, Pb, Tl, Ti, U, V) were removed from further analysis because their mean grain concentration across all plots was less than the LOD (S2 Table ).For the remaining elements, when an individual sample had a grain element concentration less than the LOD, actual values were replaced with 50% LOD.Grain element concentrations >5 standard deviation (SDs) greater than the global arithmetic mean for each element were also removed from the analysis as a precaution against using contaminated samples (15 and 31 data points were removed in 2013-14 and 2014-15, respectively).All data are provided in S3 Table.

Statistical analyses
Variance components associated with grain concentration of Zn and 15 other grain mineral elements (As, Ba, Ca, Cs, Cu, Fe, K, Mg, Mn, Mo, P, Rb, S, Se, and Sr) were calculated.Analyses of Variance (ANOVA) and Least Significant Difference (LSD) tests were used to test for differences in mean grain elemental concentration between sites and genotypes.Broad-sense heritability (H 2 ) was calculated from ANOVA for all the grain mineral elements across all sites as H 2 = V g /V p where V g is the genotypic variance and V p is the phenotypic variance.A genotype plus genotype by environment (site) interaction (GGE) analysis was performed to construct a biplot for grain Zn concentration to study the genotype and genotype by environment interaction effect, to explore the adaptation of genotypes to the specific sites.Differences between sites and genotypes for grain mineral composition traits were considered significant at P<0.01.Pearson correlation coefficients were calculated for the 16 grain mineral elements, along with grain yield, for each site in each year.All analyses were conducted using GenStat 17 th Edition (VSN International Ltd, Hemel Hempstead, UK).

Results
This study reports the contribution of genotype and site factors to variation in grain Zn concentration, and the grain ionome, of 36 wheat genotypes in India.Primary data are provided in S4 Table.

Variation in wheat grain Zn concentration and the grain ionome
The mean grain Zn concentration across all plots and years (n = 719) was 29.3 mg kg -1 and varied from 3.4 to 60.9 mg kg -1 , across all plots (Table 1).Among the macronutrients, the mean grain concentration of Ca, K, Mg, P and S were 369, 4190, 1160, 3610 and 1490 mg kg -1 , respectively and varied 20-fold for Ca, 14-fold for K, 15-fold for Mg, 189-fold for P and 17-fold for S.Among the micronutrients, the mean grain concentration of Cu, Fe, Mn and Mo across all plots were 4.3, 37.8, 37.5 and 0.74 mg kg -1 , respectively (n = 719).Grain concentration varied 12-fold for Cu, 13-fold for Fe, 11-fold for Mn and 41-fold for Mo.Among trace elements, the mean grain As and Se concentrations were 0.02 and 0.16 mg kg -1 , respectively, and varied from 0.003 to 0.089 for As and from 0.02 to 0.8 mg kg -1 for Se (Table 1).Cadmium (Cd) and lead (Pb) are trace element of interest in the context of Zn biofortification [32] because they are potentially toxic to humans [33][34][35][36].In this study, the mean grain Cd and Pb concentrations across all plots were less than the LODs of 0.0341 and 0.0660 mg kg -1 for Cd and Pb, respectively, and were therefore not included in further data analyses.

Variance components associated with grain Zn concentration and the grain ionome
Genotype (G) and site terms (E) were associated with 10% and 6% of the variation in grain Zn concentration, respectively (Table 2).The G x E interaction term was associated with 20% of the variation in grain Zn concentration.The residual term (R) was associated with 64% of the variation in grain Zn concentration, and represents variation due to year, plot-to-plot variation between replicates, and technical/measurement errors.Site (E) was typically associated with a greater proportion of the variation in the grain concentration of other minerals compared to Zn (Table 2).Among the macronutrients, 28, 18, 36, 8 and 12% of the variation in grain Ca, K, Mg, P and S concentration were associated with E; 18, 18, 11, 12 and 8% of the variation were associated with G, respectively.Among the micronutrients, E was associated with 14, 19, 55 and 79% of the variation in grain Cu, Fe, Mn and Mo concentration, respectively; 14, 9, 11 and 4% of the variation were associated with G, respectively.Among non-essential trace elements, 70 and 67% of the variation in grain As and Se concentration, respectively, were associated with E, with just 2 and 1% of the variation associated with G, respectively.
Most of the variation in yield and yield components was associated with site; 62, 55 and 75% of the variation in GYD, PHT and DTM, respectively (Table 2).Among the macronutrients and micronutrients, the grain concentrations of Ca, Mg, S, Fe, Cu and Mn differed between genotypes in both 2013-14 and 2014-15 (P<0.001).Among the non-essential trace elements, the grain As and Se concentration did not differ significantly between genotypes in either year.

Differences in grain Zn concentration and the grain ionome between genotypes
The mean grain Ca concentration of 36 genotypes, averaged across sites, varied from 285.2 to 437 mg kg -1 in 2013-14 and from 317 to 453 mg kg -1 in 2014-15.The KRL 3-4 and NW

Variation in grain Zn concentration due to variation in plant-available soil Zn
In the present study, plant-available soil Zn concentrations were reported in top and sub-surface soils at six sites in 2013-14 and 2014-15 based on the use of a non-conventional extract of Ca(NO 3 ) 2 [37,38] (Table 3).Plant available soil Zn concentration varied from 0.04-1.10mg kg -1 in 0-15 cm and from 0.03-2.14mg kg -1 in 15-30 cm soil depths in 2013-14 while in 2014-15, it varied from 0.03-0.22mg kg -1 in 0-15 cm and from 0.003-0.06mg kg -1 in 15-30 cm soil depths, at six sites.All sites showed very low concentrations of immediately plant available Zn in the soil in both years except at Pundibari in 2013-14.
At Karnal, soil available Zn concentration was greater in 2014-15 than in 2013-14 at both soil depths and also showed greater grain Zn concentration in 2014-15 than in 2013-14.Contrastingly, at Malda, the soil available Zn concentration was greater in 2013-14 than 2014-15 in both the soil depths, but, grain Zn concentration was greater in 2014-15 than in 2013-14 (Table 3).

Discussion
Wheat is an important dietary source of Zn for the Indian population and there is scope to increase grain Zn concentrations in wheat through breeding to help alleviate dietary Zn deficiency.In the present study, the grain Zn concentration of 36 wheat genotypes varied from 26-32 mg kg -1 .Our results are consistent with Zhao et al. [39] who reported grain Zn concentration from 14 to 35 mg kg -1 in 150 field grown bread wheat lines.A similar range of grain Zn concentration were reported among 66 advanced wheat genotypes (20-39 mg kg -1 ) selected from Central Asia Breeding Programs by Morgounov et al. [40], and grown in Kazakhstan, Kyrgyzstan and Tajikistan.Among 40 bread wheat lines from CIMMYT breeding programmes grown at nine locations in South Asia and Mexico, grain Zn concentration from 29-40 mg kg -1 [23].A summary of grain Zn concentrations was recently published at global scale in field grown wheat and reported grain Zn concentration from 20-31 mg kg -1 , showed a difference of 10 to 20 mg kg -1 from the biofortification target of 40 mg kg -1 for human diet [41].
Because experiments were conducted at a wide range of soil conditions, and over two seasons, genotype (G) accounted for a relatively small proportion of the total variation in grain Zn concentration.This is not surprising given that grain Zn concentration is a complex trait, under the control of many genes that can be modified under different environments.Joshi et al. [42], Velu et al. [23,43] and Swamy et al. 2016 [44] also reported a significant effect of different environments on variation in grain Zn concentration.However, some genotypes still had a consistent ranking across sites and years for grain Zn concentration.For example, Kharchia 65 and BH 1146 genotypes showed comparatively greater grain Zn concentration across six sites.Genotype, Kharchia 65 is a salt tolerant genotype [45] and BH 1146 is an aluminium tolerant genotype [46].Therefore, Kharchia 65 and BH 1146 might be potential new sources of background variation for crossing with new Zn-biofortified wheat varieties, including for saline/sodic and acidic soil areas.In addition to Kharchia 65 and BH 1146, several genotypes achieved more than 31 mg kg -1 of grain Zn concentration at different locations, which is equivalent to 50% of the HarvestPlus target of 12 mg kg -1 above a notional baseline grain Zn concentration of 25 mg kg -1 .Given that these lines already have implicitly good grain yield and other quality attributes for India, such lines could also be useful for crossing with high-Zn lines derived from the HarvestPlus, to support site-specific approaches as adopted previously for yield and yield component traits [26].
The high-Zn wheat varieties were developed by HarvestPlus using the high Zn containing accessions of synthetic wheats, spelt wheats and, further crossing with well adapted high-yielding hexaploid wheat varieties.Synthetic wheat lines were developed by crossing of Triticum durum and high-Zn containing wild accessions of tetraploid Triticum dicoccon with Aegilops squarrosa [22].The T. durum based synthetic hexaploid wheat was developed by crossing of T. durum wheat cultivar with Aegilops squarrosa, a D genome donor of wheat, and further selections have developed a high-Zn variety, the Zinc Shakti (CROC1_/AE.SQUARROSA (210)// INQALAB91 Ã 2/KUKUNA/3/PBW343 Ã 2/KUKUNA) for North-Eastern Plain Zones of India and showed 14 mg kg -1 greater grain Zn concentration than a notional baseline grain Zn concentration of 25 mg kg -1 .The T. dicoccon based synthetic hexaploid wheat was developed by crossing of a high Zn containing accession of wild T. dicoccon with Aegilops squarrosa and further selections have developed the WB02 and HPBW-01 (T.DICOCCON, CI9309/AE.SQUARROSA (409)//MUTUS/3/2 Ã MUTUS) wheat varieties for North-Western Plain Zones of India and showed 7 mg kg -1 greater grain Zn concentration than the notional baseline grain Zn concentration of 25 mg kg -1 .Similarly, Zincol-2016 has been released for the Pakistan and developed by crossing of a well-adapted NARC2011 wheat variety with a high Zn containing accession of T. spelta and showed 9 mg kg -1 greater grain Zn concentration than the notional baseline grain Zn concentration of 25 mg kg -1 [22].
There was no evidence of a yield dilution effect on grain Zn concentration in this study.Our result is therefore consistent with findings of Welch and Graham [47] and Joshi et al. [42], who reported no trade-off between grain Zn and grain yield concentration in wheat.This observation supports the argument that locally adapted, high yielding, genotypes have potential value for developing Zn-biofortified varieties for a wide range of soil conditions than reported previously.In contrast, other studies reported negative correlation of grain Zn with grain yield concentration in wheat [48,49].Oury et al. [49], reported lower grain Zn concentrations in high yielding genotypes which may have been due to yield dilution.It is important to continue to monitor grain Zn concentration in the wider context of yield and yield component traits in breeding programmes, especially given that grain yields in the region are relatively low.
Variation in plant-available Zn in soils is likely to represent a large proportion of environmental sources of variation in grain Zn concentration [11,50,51,52].Most of the soils in the present study were Zn deficient, albeit based on the use of a Ca(NO 3 ) 2 soil extraction to assess Zn-availability rather than the standard diethylene triamine pentaacetic acid (DTPA) method.Karami et al. [50], Joshi et al. [42] and Velu et al. [43] reported lower grain Zn concentrations in wheat genotypes when grown under Zn-deficient soil conditions.Soil Zn concentrations can be increased through the use of fertilisers [11,15,16].Zinc can be applied with granular fertilisers to soil or as a foliar spray.Application of Zn fertiliser in granular form has produced increases in both grain yield and grain Zn concentration.However, application of Zn fertiliser as a foliar spray can reduce the amount Zn fertiliser needed in granular form and cereals are more responsive in terms of increase in grain Zn concentration, albeit at higher application costs.It is possible to reduce the cost of application by combining foliar application of Zn with pesticide applications [6,17,53].Based on a wide survey of farmers, application of basal Zn fertilisers increased grain yield in wheat in Indo-Gangetic area of Sindh and Punjab Provinces of Pakistan [16].A novel aspect of this study is that the same panel was grown in an extremely diverse range of hostile soils that will have complex trace element geochemistry.For example, sodic soils at the Kumarganj site had a pH value of 9.5, which is much higher than in the arable soils routinely studied in crop trials.At such sites, it is likely that reclamation using addition of gypsum (calcium sulphate, CaSO 4 .2H 2 O), and farm yard manure (FYM) coupled with measures to offset soil structural damage and impeded drainage would be needed to optimise grain yields and quality.
Increases in grain Zn concentration may also help to increase the grain Ca, Cu, Fe, K, Mg, Mn, Mo, P, and S concentrations given the observed positive correlation of grain Zn with these elements.The ionome data of high-Zn wheat varieties released by HarvestPlus in India and Pakistan is not available because the standard HarvestPlus method to measure the whole grain Zn concentration is X-ray fluorescence (XRF) [54].However, our results are consistent with those of Pandey et al. [55], who reported a positive relationship between grain Zn concentrations and Ca, Cu, Fe, K, Mg, P and S concentrations in 150 bread wheat lines collected from India and Turkey.Many researchers also reported relationships between grain Zn and Fe [56,23,57], P [58,59] and S [60].Ozkan et al. [61] reported a major QTL on chromosome 5 which increased the grain Fe, Zn, Cu and Mn concentration in Triticum monococcum wheat genotypes.Velu et al. [57] reported co-localisation of grain Zn and Fe QTLs on 2B chromosomes which may help to simultaneous improvement of both the micronutrients.Similarly, McDonald and Mousavvi [60] showed that slightly increase in sulphur concentration increased the grain Zn concentration.It may be due to increase in S containing amino acid methionine, which further increased the level of phytosiderophores and nicotinamide, involved in uptake and translocation process of zinc.However, P is an important component of phytic acid (PA) which affects the bioavailability of Zn [59].Therefore, the role of PA should be considered concomitantly in breeding programs intended to increase the grain Zn concentration in wheat genotypes.
Taken together, the results from this study show considerable potential for genetic improvement for Indian wheat genotypes.Some elite varieties (e.g.Kharchia 65) could be useful sources of variation that can be used at a wide range of sites to cross with new Zn-biofortified lines.This study also reveals the wealth of information that can be gleaned from studying the wider grain ionome across sites with a wide range of conditions, so that both breeding and spatial approaches (e.g.site-specific crop selection, fertiliser-use, etc.) to improving diet quality can deployed most effectively.excluded from analysis (values below than LOD) and values higher than mean+(5 Ã SD) (highlighted in yellow) are removed from data before analysis.(PDF) S4 Table.Data of grain mineral composistion traits (mg kg -1 ) of 36 wheat genotypes grown at five sites in 2013-14 and at six sites in 2014-15 in India.Data below limit of detection (LOD) have been replaced by half LOD values and values above than mean+(5 Ã SD) have been removed.(PDF) the Indian Institute of Wheat & Barley Research (IIWBR), Karnal (29.70˚N; 76.99˚E), Haryana (HR); (2) IIWBR, Hisar (29.18˚N; 75.70˚E) (HR); (3) Narendra Deva University of Agriculture and Technology (NDUA&T), Kumarganj (26.43˚N; 82.17˚E), Faizabad, Uttar Pradesh (UP), reclaimed site; (4) NDUA&T, Kumarganj (26.43˚N; 82.17˚E), Faizabad (UP), sodic site; (5) Uttar Banga Krishi Viswavidyalaya (UBKV), Regional Research Sub-Station (RRSS), Mathurapur (25.57˚N; 87.10˚E), Malda (WB) and (6) UBKV, Pundibari (26.32˚N; 89.45˚E), Cooch Behar, West Bengal (WB).Seeds were manually placed in four rows per plot at 25 cm spacing between rows.The plot length was 2.5 m.The plots were arranged in a simple lattice design (6x6) with two replicates according to standard IIWBR practices.The field experiments were conducted at research farms of IIWBR, Karnal; NDUA&T, Kumarganj and UBKV, Pundibari, with their research director's permission and study received the permission from the Indian Council of Agricultural Research (ICAR) and Department of Biotechnology (DBT) under a collaborative project between India and University of Nottingham, UK; grant number BT/IN/UK/12/IS/2012.

Fig 2 .
Fig 2. The correlation coefficients in grain ionome of a panel of 36 genotypes between 2013-14 and 2014-15.Data represent the means of two replicate plots per genotype at Karnal, Hisar and Malda, and one replicate per genotype at Kumarganj-reclaimed and Kumarganj-sodic sites in 2013-14, and the means of two replicate plots of each genotype at all six sites in 2014-15.Colour represents strength of correlation from strongly negative (dark blue) to strongly positive (dark red).https://doi.org/10.1371/journal.pone.0192026.g002

Fig 3 .
Fig 3. Grain element concentrations of 36 genotypes of wheat at six sites over two years.Data represent the means of two replicate plots per genotype at Karnal, Hisar, Malda and one replicate plot per genotype at Kumarganj-reclaimed and Kumarganj-sodic sites in 2013-14 and two replicate plots per genotype at all six sites in 2014-15.Boxes represent the two mid-quartiles with the median drawn; whiskers are the 95% confidence limits; circles are the outliers.https://doi.org/10.1371/journal.pone.0192026.g003

Fig 5 .
Fig 5.The relationship of grain Zn concentration (mg kg -1 ) with GYD (t ha -1 ) and other grain mineral elements in a panel of 36 wheat genotypes grown at five sites in 2013-14 and six sites in 2014-15.Data are means of two replicate per plot at Karnal, Hisar and Malda and one replicate at Kumarganj-reclaimed and Kumarganj-sodic sites in 2013-14 and two replicate per plot at six sites in 2004-15 (n = 719).Colour represents strength of correlation from strongly negative (dark blue) to strongly positive (dark red).https://doi.org/10.1371/journal.pone.0192026.g005

Table 2 . The contribution of G, E, G Ã E and residual factors to variation percentage (%) in grain Zn concentration and the grain ionome of a panel of 36 wheat geno- types grown at five sites in 2013-14 and at six sites in 2014-15 and in grain yield and yield components of a panel of 36 genotypes grown at six sites over two years.
https://doi.org/10.1371/journal.pone.0192026.t002

Table 3 . Plant available soil Zn (Ca(NO 3 ) 2 -extractable) and grain Zn concentration at six sites in 2013-14 and 2014-15.
Soil Zn data are means of three replicate soil samples of each depth at each site.Grain Zn data are mean of two replicate at Karnal, Hisar, Malda and one replicate at Kumarganj-reclaimed and Kumarganj-sodic sites in 2013-14 and two replicates at all six sites in 2014-15.