Jia Liu1,2†
Lin Huang2†
Tingxuan Li3
Yaxi Liu1,2,4
Zehong Yan1,2,4
Guan Tang2
Youliang Zheng1,2,4
Dengcai Liu1,2,4
Bihua Wu1,2,4*- 1State Key Laboratory of Crop Gene Exploration and Utilization in Southwest China, Sichuan Agricultural University, Wenjiang, China
- 2Triticeae Research Institute, Sichuan Agricultural University, Wenjiang, China
- 3College of Resources, Sichuan Agricultural University, Wenjiang, China
- 4Key Laboratory of Crop Genetic Resources and Improvement, Ministry of Education, Sichuan Agricultural University, Ya'an, China
Wheat is one of the important staple crops as the resources of both food and micronutrient for most people of the world. However, the levels of micronutrients (especially Fe and Zn) in common wheat are inherently low. Biofortification is an effective way to increase the micronutrient concentration of wheat. Wild emmer wheat (Triticum turgidum ssp. dicoccoides, AABB, 2n = 4x = 28) is an important germplasm resource for wheat micronutrients improvement. In the present study, a genome-wide association study (GWAS) was performed to characterize grain iron, zinc, and manganese concentration (GFeC, GZnC, and GMnC) in 161 advanced lines derived from wild emmer. Using both the general linear model and mixed linear model, we identified 14 high-confidence significant marker-trait associations (MTAs) that were associated with GFeC, GZnC, and GMnC of which nine MTAs were novel. Six MTAs distributed on chromosomes 3B, 4A, 4B, 5A, and 7B were significantly associated with GFeC. Three MTAs on 1A and 2A were significantly associated with GZnC and five MTAs on 1B were significantly associated with GMnC. These MTAs show no negative effects on thousand kernel weight (TKW), implying the potential value for simultaneous improvement of micronutrient concentrations and TKW in breeding. Meanwhile, the GFeC, GZnC and GMnC are positively correlated, suggesting that these traits could be simultaneously improved. Genotypes containing high-confidence MTAs and 61 top genotypes with a higher concentration of grain micronutrients were recommended for wheat biofortification breeding. A total of 38 candidate genes related to micronutrient concentrations were identified. These candidates can be classified into four main groups: enzymes, transporter proteins, MYB transcription factor, and plant defense responses proteins. The MTAs and associated candidate genes provide essential information for wheat biofortification breeding through marker-assisted selection (MAS).
Introduction
Mineral elements are essential micronutrients for animals and human beings. It is estimated that approximately half of the world's population are at risk of micronutrient deficiency (Gregory et al., 2017; Tabbita et al., 2017), resulting in serious health problems (Cakmak and Kutman, 2018). For example, the deficiency in iron (Fe) and zinc (Zn) may cause overall poor health, anemia, increase morbidity and mortality rates, and low worker productivity (Peleg et al., 2008; Gibson, 2012; Borrill et al., 2014), whereas the manganese (Mn) shortfall may affect reproductive ability, mental development, and bone development (Li and Tan, 2012). In addition, the previous study showed that crops grown in high CO2 environments not only have low nutritional value, but also lack Fe, Zn, and other important micronutrients (Myers et al., 2014). Developing micronutrient-enriched agricultural crops (biofortification), agronomically and/or genetically, is considered the most cost-effective and sustainable approach to alleviate malnutrition and related health problems (Peleg et al., 2008; Bouis and Saltzman, 2017).
Wheat (Triticum aestivum, 2n = 6x = 42, AABBDD) is a major staple food crop worldwide. Therefore, the composition and nutritional quality of wheat grains have a significant impact on human health and well-being. Modern wheat cultivars are, however, inherently low in micronutrients and show a narrow genetic variation for micronutrients to be exploited in breeding programs (Velu et al., 2017). One effective strategy to improve the amount of mineral elements in wheat grains is to exploit the “left behind” genetic variation in wild relatives for grain micronutrients (Peleg et al., 2008, 2009).
Wild emmer wheat (Triticum turgidum ssp. dicoccoides, 2n = 4x = AABB) is the tetraploid progenitor of cultivated wheat, harboring a rich allelic repertoire for improvement of various economically important traits in wheat (Gomez-Becerra et al., 2010), including micronutrient concentrations (Yan et al., 2018). Substantially high concentrations of up to 190 mg kg−1 for Zn and 109 mg kg−1 for Fe in grains of wild emmer have been reported (Cakmak et al., 2004, 2010).
The accumulation of minerals in seeds depends on a plethora of processes that are controlled by several genes (Peleg et al., 2009). Till now, dissection of the genetic basis of grain micronutrients content in wild emmer wheat was performed using tetraploid wheat recombinant inbred lines (RILs), derived from a cross between durum wheat (T. turgidum ssp. durum, cv. Langdon) and wild emmer wheat (accession G18-16) (Peleg et al., 2009; Fatiukha et al., 2020). The introgression and identification of wild emmer quantitative trait locus (QTL) linked to micronutrient accumulation in a hexaploid wheat background have been less reported. The wild emmer gene Gpc-B1 that belongs to the NAC-domain transcription factor affects grain protein content (GPC), grain Zn, Fe, and Mn concentrations (GZnC, GFeC, and GMnC) in wheat (Uauy et al., 2006; Distelfeld et al., 2007). However, this gene was associated with reductions in grain weight and yield in some environments and wheat cultivars (Uauy et al., 2006; Brevis and Dubcovsky, 2010; Tabbita et al., 2013). Thus, the exploration of genomic regions associated with high GZnC, GFeC, and GMnC that have less negative effect on grain yield is desirable.
Genome-wide association analysis (GWAS) has been widely used for deciphering the genetic basis of multiple traits in crops (Su et al., 2016; Wang et al., 2016; Ates et al., 2018; Liu et al., 2019). It was available to study QTLs related to agronomically important traits in large sets of germplasm resources such as landraces (Liu et al., 2017), elite cultivars (Sukumaran et al., 2015), and advanced breeding lines (Wang et al., 2017) as well as backbone parents and their derived lines (Yu and Tian, 2012; Yu et al., 2014; Xiao et al., 2016; Liu et al., 2019). In wheat, GWAS has been extensively applied to reveal genomic regions controlling traits such as GPC (Liu et al., 2019), grain ionome (Fatiukha et al., 2020), and yield-related traits (Tadesse et al., 2015).
In our previous study, the agronomically stable advanced lines were obtained from a cross between common wheat cultivar Chuannong16 (CN16 hereafter) and wild emmer accession D1 (Liu et al., 2019). We found that most of the tested advanced lines without Gpc-B1 showed thousand kernel weight (TKW), GPC and grain mineral concentrations (e.g. GZnC and GFeC) simultaneous improvement (Wang, 2015). In this study, GWAS was used to analyze the genetic basis of GFeC, GZnC, and GMnC in a multi-parent population which consisted of wild emmer as backbone parent and its derived advanced lines. The objectives of the current study were to identify genomic regions associated with high GFeC, GZnC, and GMnC in a hexaploid wheat population, and to scan promising candidate genes using the whole genome assembly of wild emmer wheat (Avni et al., 2017) and Chinese Spring (Appels et al., 2018). The identified markers and genomic regions will provide essential information for cloning genes related to high GFeC, GZnC, and GMnC and be useful in biofortification breeding programs.
Materials and Methods
Plant Materials and Experimental Design
The same set of 161 advanced lines (Liu et al., 2019) derived from a cross between CN16 and D1 were used in the present study. Wheat plants were arranged in the field using a randomized complete block design with three replicates over two growing seasons (2015 and 2016) at the Chongzhou (2015CZ and 2016CZ) and Wenjiang (2015WJ and 2016WJ) experimental stations of Sichuan Agricultural University, China (Liu et al., 2019). Individual plants were spaced 10 cm apart within a 2 m row with 30 cm between rows. Each replicate contained 20 plants in a 2 m row. Cultivation and management followed local field production conditions. The soil types at fields of Wenjiang and Chongzhou experimental stations are paddy soil and yellow soil, respectively. Mature seeds were harvested from the middle six wheat plants of each row and measured for GFeC, GZnC, GMnC, and TKW.
Phenotypic Measurements
Wheat grains were dried to a constant weight and then ground to a fine powder using Chopin CD1 AUTO (Renault, Boulogne-Billancourt, France). A 0.3 g powder samples were digested with a diacid mixture [HNO3 (4): HClO4 (1) (v/v)]. Atomic absorption Spectrometer (PinAAcle 900T, USA) was used to measure the GFeC, GZnC, and GMnC. The weight of 300 randomly sampled kernels (GB T 5519-2008, 2008) was recorded with an electronic balance to represent the TKW.
Genotyping and Population Structure and Linkage Disequilibrium (LD) Analysis
The wheat genomic DNA was isolated from seedlings using a CTAB method (Murray and Thompson, 1980) with minor modification and then, genotyped using the DArT markers (https://www.diversityarrays.com/). The obtained data were filtered according to call rate (minimum threshold value of 85%) and reproducibility (minimum threshold value of 95%) (Liu et al., 2019). The recalled marker data with missing data >10% and minor allele frequency (MAF) <5% were further filtered (Liu et al., 2019). The retained 13,116 markers were used for further study. The linkage disequilibrium (LD) and population structure analyses were conducted as described by Liu et al. (2019).
GWAS for GZnC, GFeC, and GMnC
The best linear unbiased prediction (BLUP) across four tested environments was performed using the META-R (Alvarado et al., 2015). The association analysis of markers and GZnC, GFeC, and GMnC was conducted using the general linear model (GLM) and mixed linear model (MLM) in TASSEL. The population structure (Q matrix) was used as a covariate to adjust population stratification. Kinship matrix (K) was calculated using the Scaled IBS method (Yu et al., 2006). A threshold p-value of –log10 (p) ≥ 3 (Alomari et al., 2018) was used to define significant marker-trait associations (MTAs). MTAs were displayed with a Manhattan plot using Haploview4.2 software (Barrett et al., 2004). Important p-value distributions were shown with a quantile-quantile plot. The prediction of candidate genes linked to MTAs was performed as described by Liu et al. (2019).
Statistical Analysis
Analysis of variance (ANOVA) and estimates of heritability were performed using the META-R (Alvarado et al., 2015). T-test and the phenotypic Pearson correlation coefficients were calculated using SPSS version 22.0 (SPSS Inc., Chicago, IL, USA).
Results
Phenotypic Variation for GFeC, GZnC, and GMnC
The variation for GFeC, GZnC, and GMnC in the advanced lines and their core parents were showed in Table 1. A wide range of phenotypic variation was found among the wheat advanced lines. These traits were significantly different among genotypes and showed medium-to-high heritability (0.67–0.98) (Table 1). Most genotypes displayed relatively stable GFeC, GZnC, and GMnC across four environments (2015WJ, 2015CZ, 2016WJ, and 2016CZ) (Figure 1).
Table 1. The variation for GFeC, GZnC, GMnC, and TKW in advanced lines under four environments.
Figure 1. Heat map of GFeC, GZnC, and GMnC of the advanced lines grown at 2015WJ, 2015CZ, 2016WJ, and 2016CZ.
The wild emmer D1 had significantly (P < 0.05) higher GFeC (mean range 92.69–115.47 mg/kg) compared to that of CN16 (mean range 38.15–49.20 mg/kg) across all test environments (Table 1). The GFeC of the advanced lines were ranged from 42.44 to 168.00 mg/kg (mean range 83.71–109.68 mg/kg) across four environments. The lowest GFeC (42.44 mg/kg) was recorded at 2016CZ, whereas the highest GFeC (168.00 mg/kg) was recorded at 2015 WJ.
The D1 showed significantly (P < 0.05) higher GZnC (mean range 80.91–97.87 mg/kg) compared to that of CN16 (mean range 36.67–41.64 mg/kg) across four environments (Table 1). The GZnC of the advanced lines ranged from 55.99 to 68.75 mg/kg across environments. The lowest (30.17 mg/kg) and highest GZnC (120.52 mg/kg) were recorded at 2015CZ.
There is no significant difference in GMnC between D1 and CN16, except in two environments (2015CZ and 2015WJ). The GMnC of the advanced lines were ranged from 16.28 to 53.85 mg/kg (mean range 30.36–37.14 mg/kg) across four environments. The lowest GMnC (16.28 mg/kg) was recorded in 2016 at Wenjiang, whereas the highest GMnC (53.85 mg/kg) was recorded in 2015 at Wenjiang (Table 1).
The mean TKW of the advanced lines under 2015WJ, 2015CZ, 2016WJ, and 2016CZ environments was 46.08, 52.90, 49.88, and 49.06 g, respectively, which were all higher than those of female parent CN16 (Table 1).
Phenotypic Frequency Distributions and Pearson's Correlation Analysis
Frequency distributions of advanced lines for GFeC, GZnC, and GMnC with an indication of the bi-parents controls were presented in Figure 2. The advanced lines exhibited normal distributions for all micronutrient variables under four environments. The highest frequency distribution of GFeC showed a large difference in four environments, implying this variable was subjected to environmental factors. In contrast, the highest frequency distributions of GZnC, GMnC, and TKW were clustered in a narrow range of 50–70 mg/kg, 30–40 mg/kg, and 45–55 g, respectively (Figure 2). Most wheat lines had both higher grain micronutrient concentrations and TKW than those of CN16 in four tested environments (Supplementary Table 1).
Figure 2. Frequency distribution of GFeC, GZnC, GMnC, and TKW under four environments in the advanced lines. Solid arrows represent the values of D1. Dashed arrows represent the values of CN16.
Pearson's correlation analysis revealed highly significant positive correlations among the GFeC, GZnC, and GMnC base on BLUP. The TKW had a significant positive correlation (P < 0.05) with GMnC and no negative correlation with GFeC and GZnC (Table 2).
Table 2. Pearson's correlation coefficients of GFeC, GZnC, GMnC, and TKW in advanced lines based on BLUP data.
Population Structure and Linkage Disequilibrium (LD) Analysis
Based on population structure and LD analysis, the 161 RILs could be divided into three subgroups and the LD decay distances for the A, B, and D subgenomes and the entire genome were about 9, 12, 13, and 12 cM, respectively, which were rapidly decreased with increasing pairwise distance as described previously (Liu et al., 2019).
GWAS for GFeC, GZnC, and GMnC in Advanced Lines
Based on GLM, 556 MTAs relative to grain micronutrient concentrations were distributed on all wheat chromosomes with phenotypic variance explained (PVE) of 5.45–18.70% (Table 3). In the MLM, 22 MTAs were detected with PVEs of 7.18–11.84% (Table 3). The Manhattan and quantile-quantile charts based on the GWAS were presented in Figure 3. The GLM and MLM models as well as quantile-quantile charts were considered together to remove the potential false-positive associations. Fourteen high-confidence MTAs were detected by both GLM and MLM (Table 4). Of these, six MTAs relative to GFeC were found on chromosomes 3B, 4A, 4B, 5A, and 7B (Table 4). Three MTAs for GZnC were clustered on chromosomes 1A and 2A (Table 4). Five MTAs for GMnC were on chromosome 1B (Table 4). According to the results of GWAS, the top-ranking 61 (37.9%) genotypes with a high number of MTAs that associated with GFeC, GZnC, and GMnC (Figure 4, Supplementary Table 2) were identified. A wide range of phenotypic variation for GFeC (range 62.73–125.58 mg/kg), GZnC (range 47.00–75.09 mg/kg), GMnC (range 21.71–49.50 mg/kg), and TKW (range 42.64–53.87 g) was observed among these genotypes (Supplementary Table 2). Moreover, t-test showed that the average values of GFeC (100.62 mg/kg), GZnC (60.74 mg/kg), GMnC (34.56 mg/k), and TKW (49.54 g) for the 61 genotypes were significantly higher than those of the parent CN16 (p < 0.001) (Supplementary Table 2).
Table 3. Significant MTAs identified by GWAS using GLM and MLM model.
Figure 3. Manhattan plots of GWAS scan for GFeC, GZnC, and GMnC in four environments. Significant MTAs identified by MLM (A) and GLM (B). Red lines: the –log10(p) threshold of 3.00. The quantile-quantile charts of MLM (C) and GLM (D). Black line: the expected values.
Table 4. Significant high-confidence MTAs for GFeC, GZnC, and GMnC identified by MLM and GLM.
Figure 4. Potential genotypes with high micronutrient concentrations selected from the RILs based on high-confidence MTAs associated with GFeC, GZnC, and GMnC.
Candidate Genes That May Be Associated With GFeC, GZnC, and GMnC
The high-confidence MTAs were selected to predicate candidate genes using the annotated Chinese Spring genome sequence (RefSeq v.1.0) and wild emmer genome sequence (WEWSeq v.1.0). A total of 38 putative candidate genes for micronutrient concentrations were predicted (Supplementary Table 3). These genes can be roughly classified into four groups based on the types of proteins they encoded. The first group of candidates encoded enzymes related to iron-dependent dioxygenase, metal-dependent hydrolase, phosphatase, and methyltransferases. The second group consisted of transporter proteins, such as sugar transporter, ABC transporter, and heavy metal transport proteins. The third group included MYB transcription factor, zinc finger family, and transmembrane proteins. The last group consisted of proteins related to plant defense responses, such as receptor-like kinase, leucine-rich repeat-containing protein, and mitogen-activated kinase (MAPK) (Supplementary Table 3).
Discussion
The poor bioavailability and low concentration of grain micronutrients (especially Zn and Fe) in cultivated wheat are common problems (Xu et al., 2012). Enhancement of the wheat grain nutritional value through genetic biofortification is available (Borrill et al., 2014). The wild emmer accessions are a valuable resource for micronutrient improvement in wheat (Chatzav et al., 2010). In the present study, the same set of 161 advanced lines derived from wild emmer (Liu et al., 2019) were phenotyped in four environments and genotyped using the DArT markers to understand the genetic basis of GFeC, GZnC, and GMnC accumulation in wheat. GWAS identified 14 high-confidence MTAs for GFeC, GZnC, and GMnC by both GLM and MLM (Table 4), and 38 putative candidate genes were predicate using the wheat sequences annotation (RefSeq v.1.0; WEWSeq v.1.0).
We have found that the GFeC, GZnC and GMnC of advanced lines derived from wild emmer were significantly higher than those of CN16, with the exception of GMnC in 2016WJ and 2016CZ. Medium-to-high heritability (0.67–0.98) were estimated for GFeC, GZnC, GMnC, and TKW, which indicated these traits were largely governed by genetic factors. Similar heritability for these traits has been reported in previous studies (Velu et al., 2017; Bhatta et al., 2018). These results demonstrated the potential value of wild emmer for grain micronutrient improvement in wheat. On the other hand, we have observed significant positive correlations among GFeC, GZnC, and GMnC and no significant negative correlations between grain micronutrient concentrations and TKW. Our results were consistent with a previous study that the wild emmer GZnC and GFeC were positively correlated (Peleg et al., 2009). These results indicate that the alleles for GZnC and GFeC are either co-segregated or have a pleiotropic effect.
Six MTAs on chromosomes 3B, 4A, 4B, 5A (2), and 7B are related to GFeC. Two MTAs (2258533 and 3033960) at position 0–0.38 cM on chromosome 5A are located in the same GFeC-QTL region as reported by Peleg et al. (2009). The MTA 2255722 (46.59 cM) was located in a GFeC-QTL region between markers wmc468-barc170 at position 32.6–50.8 cM on chromosome 4A (Pu et al., 2014). The remaining 3 MTAs on 3B, 4B, and 7B in our studies might be novel.
Among the three MTAs related to high GZnC, two MTAs (1077698 and 1234362) on chromosome 2A at position 139.90–145.15 cM were located in the same interval of wild emmer GZnC-QTL (112.4 ± 35.0 cM) (Peleg et al., 2009). The position of MTA 1136167 (359.46 cM) in the current study is different to the GZnC-QTL (39 cM) on chromosome 1A (Roshanzamir et al., 2013).
The five MATs related to high GMnC located on chromosome 1B. Several GMnC-QTLs have been reported on different wheat chromosomes such as 1A, 2A, 2B, 2D, 3A, 3B, 4B, 4D, 5A, 5D, 6B, and 7B (Shi et al., 2013; Pu et al., 2014; Bhatta et al., 2018). To our best knowledge, there have no reports of GMnC QTLs on chromosome 1B. Noteworthy, we have found that some wheat genotypes carrying MTAs for high GFeC, GZnC, and GMnC, but without reduction in TKW (Figure 4, Supplementary Table 2). The GFeC, GZnC, GMnC, and TKW of these wheat genotypes were significantly higher than those of the parent CN16 (p < 0.001) (Supplementary Table 2). These results indicate the presence of novel wild emmer alleles that would be used for grain minerals biofortification.
In the present study, the high-confidence MTAs were used to predicate putative candidate genes. We found that MYB transcription factor, dioxygenase, and hydrolase genes may play important roles in conferring high GFeC. For example, the dioxygenase genes in higher plants were reported to be involved in the absorption and transport of Fe (Kobayashi and Nishizawa, 2012). The MYB transcription factor in Malus xiaojinensis was responded to Fe deficiency stress (Han et al., 2020). In addition, we have identified hydrolase family protein/HAD-superfamily protein, thioredoxin-like protein AAED1, kelch repeat-containing protein, receptor-like protein kinase, and HD domain-containing metal-dependent phosphohydrolase family protein on chromosome 4A and 3B. The previous study showed that the kelch repeat-containing protein was associated with Fe acquisition and metal homeostasis (Kawahara et al., 2017).
Our results indicated that zinc ion binding (FAR1), heavy metal transport, zinc finger (C3HC4-type RING finger) family, and mitogen-activated protein kinase (MAPK) may be promising candidates for GZnC (Supplementary Table 3). Previous studies found that the metal transporter (ZIP1) could transport Zn in wild emmer wheat (Durmaz et al., 2011). A recent study reported that a MAPK-related gene may be involved in Zn accumulation in chickpea seeds (Upadhyaya et al., 2016; Alomari et al., 2018).
We have identified S-adenosyl-L-methionine-dependent methyltransferases, ABC transporter, phosphatase, and argininosuccinate synthase as the candidates for GMnC. For example, the ABC transporter was proved to be involved in the uptake and transport of manganese ions (Li and Tan, 2012). Arginase is a type of manganese enzyme, and therefore argininosuccinate synthase can lead to Mn accumulation (Li and Tan, 2012).
Previous studies have demonstrated that the grain yield and micronutrients (Fe, Zn, and Mn) are negatively correlated in wheat, making the simultaneous improvement of these two traits challenging (Shi et al., 2013; Bhatta et al., 2018). However, Peleg et al. (2008) reported neither negative nor positive associations between grain micronutrient concentrations and yield in wild emmer wheat. In our previous studies, we have introgressed the high GFeC, GZnC, and TKW traits of wild emmer into common wheat through wide hybridization (Wu et al., 2008; Wang, 2015). In this study, we found no significant negative correlation between grain micronutrient concentrations (GFeC, GZnC, and GMnC) and TKW in the advanced lines. Our findings demonstrate the improvement of grain micronutrient concentrations (GFeC, GZnC, and GMnC) without sacrificing grain yield could be possible by the exploitation of wild emmer MTAs.
Conclusion
Our results show the importance of wild emmer as valuable germplasm for the improvement of grain micronutrient concentrations in wheat. The advanced lines derived from wild emmer showed high grain micronutrients and a weak correlation between grain micronutrients and TKW, indicating simultaneous improvement of grain micronutrient concentrations (especially GFeC and GZnC) and TKW could be possible. The top-ranking 61 wheat genotypes carrying MTAs for high GFeC, GZnC, and GMnC without reduction in TKW could be valuable donor parents for wheat biofortification breeding.
A GWAS identified 14 high-confidence significant MTAs that were associated with GFeC, GZnC, and GMnC, of which nine MTAs were novel. In total, 38 putative candidate genes for these MTAs were predicated. The MTAs and associated candidate genes provide information for fine mapping and cloning of genes related to GFeC, GZnC, and GMnC and for wheat biofortification breeding programs.
Data Availability Statement
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author.
Author Contributions
BW, JL, and LH designed the experiment. JL, TL, and GT performed experiments. JL analyzed data and drafted the manuscript. YL and ZY provided the genotyping data. YZ, DL, and BW provided their valuable suggestions. LH and BW revised the manuscript. All authors read and approved the final manuscript.
Funding
This work was financially supported by grants from the National Natural Science Foundation of China (Grant No. 31571668) and the National Key Research and Development Program of China (Grant No. 2017YFD0100900).
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
The reviewer VG declared a past collaboration with two of the authors JL and BW to the handling editor.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2021.651283/full#supplementary-material
References
Alomari, D. Z., Eggert, K., Von Wirén, N., Alqudah, A. M., Polley, A., Plieske, J., et al. (2018). Identifying candidate genes for enhancing grain Zn concentration in wheat. Front. Plant Sci. 9:1313. doi: 10.3389/fpls.2018.01313
Alvarado, G., López, M., Vargas, M., Pacheco, Á., Rodríguez, F., Burgueño, J., et al. (2015). META-R (Multi Environment Trail Analysis with R for Windows) Version 5.0. Texcoco: International Maize and Wheat Improvement Center.
Appels, R., Eversole, K., Feuillet, C., Keller, B., Rogers, J., Stein, N., et al. (2018). Shifting the limits in wheat research and breeding using a fully annotated reference genome. Science 361:eaar7191. doi: 10.1126/science.aar7191
Ates, D., Asciogul, T. K., Nemli, S., Erdogmus, S., Esiyok, D., and Tanyolac, M. B. (2018). Association mapping of days to flowering in common bean (Phaseolus vulgaris L.) revealed by DArT markers. Mol. Breed. 38:113. doi: 10.1007/s11032-018-0868-0
Avni, R., Nave, M., Barad, O., Baruch, K., Twardziok, S. O., Gundlach, H., et al. (2017). Wild emmer genome architecture and diversity elucidate wheat evolution and domestication. Science 357, 93–97. doi: 10.1126/science.aan0032
Barrett, J. C., Fry, B., Maller, J., and Daly, M. J. (2004). Haploview: analysis and visualization of LD and haplotype maps. Bioinformatics 21, 263–265. doi: 10.1093/bioinformatics/bth457
Bhatta, M., Baenziger, P., Waters, B., Poudel, R., Belamkar, V., Poland, J., et al. (2018). Genome-wide association study reveals novel genomic regions associated with 10 grain minerals in synthetic hexaploid wheat. Int. J. Mol. Sci. 19:E3237. doi: 10.3390/ijms19103237
Borrill, P., Connorton, J., Balk, J., Miller, T., Sanders, D., and Uauy, C. (2014). Biofortification of wheat grain with iron and zinc: integrating novel genomic resources and knowledge from model crops. Front. Plant Sci. 5:53. doi: 10.3389/fpls.2014.00053
Bouis, H. E., and Saltzman, A. (2017). Improving nutrition through biofortification: a review of evidence from HarvestPlus, 2003 through 2016. Glob. Food Secur. 12, 49–58. doi: 10.1016/j.gfs.2017.01.009
Brevis, J. C., and Dubcovsky, J. (2010). Effects of the chromosome region including the Gpc-B1 locus on wheat grain and protein yield. Crop Sci. 50, 93–104. doi: 10.2135/cropsci2009.02.0057
Cakmak, I., and Kutman, U. Á. (2018). Agronomic biofortification of cereals with zinc: a review. Eur. J. Soil Sci. 69, 172–180. doi: 10.1111/ejss.12437
Cakmak, I., Pfeiffer, W. H., and McClafferty, B. (2010). Biofortification of durum wheat with zinc and iron. Cereal Chem. 87, 10–20. doi: 10.1094/CCHEM-87-1-0010
Cakmak, I., Torun, A., Millet, E., Feldman, M., Fahima, T., Korol, A., et al. (2004). Triticum dicoccoides: an important genetic resource for increasing zinc and iron concentration in modern cultivated wheat. Soil Sci. Plant Nutr. 50, 1047–1054. doi: 10.1080/00380768.2004.10408573
Chatzav, M., Peleg, Z., Ozturk, L., Yazici, A., Fahima, T., Cakmak, I., et al. (2010). Genetic diversity for grain nutrients in wild emmer wheat: potential for wheat improvement. Ann. Bot. 105, 1211–1220. doi: 10.1093/aob/mcq024
Distelfeld, A., Cakmak, I., Peleg, Z., Ozturk, L., Yazici, A. M., Budak, H., et al. (2007). Multiple QTL effects of wheat Gpc-B1 locus on grain protein and micronutrient concentrations. Physiol. Plantarum 129, 635–643. doi: 10.1111/j.1399-3054.2006.00841.x
Durmaz, E., Coruh, C., Dinler, G., Grusak, M. A., Peleg, Z., Saranga, Y., et al. (2011). Expression and cellular localization of ZIP1 transporter under zinc deficiency in wild emmer wheat. Plant Mol. Biol. Rep. 29, 582–596. doi: 10.1007/s11105-010-0264-3
Fatiukha, A., Klymiuk, V., Peleg, Z., Saranga, Y., Cakmak, I., Krugman, T., et al. (2020). Variation in phosphorus and sulfur content shapes the genetic architecture and phenotypic associations within the wheat grain ionome. Plant J. 101, 555–572. doi: 10.1111/tpj.14554
GB T 5519-2008 (2008). Cereals and Pulses-Determination of the Mas of 1000 Grains. Beijing: Standards Press of China.
Gibson, R. S. (2012). Zinc deficiency and human health: etiology, health consequences, and future solutions. Plant Soil 361, 291–299. doi: 10.1007/s11104-012-1209-4
Gomez-Becerra, H. F., Yazici, A., Ozturk, L., Budak, H., Peleg, Z., Morgounov, A., et al. (2010). Genetic variation and environmental stability of grain mineral nutrient concentrations in Triticum dicoccoides under five environments. Euphytica 171, 39–52. doi: 10.1007/s10681-009-9987-3
Gregory, P. J., Wahbi, A., Adu-Gyamfi, J., Heiling, M., Gruber, R., Joy, E. J., et al. (2017). Approaches to reduce zinc and iron deficits in food systems. Glob. Food Secur. 15, 1–10. doi: 10.1016/j.gfs.2017.03.003
Han, D., Zhou, Z., Du, M., Li, T., Wu, X., Yu, J., et al. (2020). Overexpression of a Malus xiaojinensis WRKY transcription factor gene (MxWRKY55) increased iron and high salinity stress tolerance in Arabidopsis thaliana. In Vitro Cell. Dev-Pl. 56, 600–609. doi: 10.1007/s11627-020-10129-1
Kawahara, Y., Hashimoto, T., Nakayama, H., and Kitamura, Y. (2017). Galactose oxidase/kelch repeat-containing protein is involved in the iron deficiency stress response in the roots of Hyoscyamus albus. Plant Root 11, 58–63. doi: 10.3117/plantroot.11.58
Kobayashi, T., and Nishizawa, N. K. (2012). Iron uptake, translocation, and regulation in higher plants. Annu. Rev. Plant Biol. 63, 131–152. doi: 10.1146/annurev-arplant-042811-105522
Li, W., and Tan, X. (2012). Manganese trafficking, metabolism and homeostasis. Chin. Bull. Life Sci. 24, 867–880. doi: 10.1111/jnc.13170
Liu, J., Huang, L., Wang, C., Liu, Y., Yan, Z., Wang, Z., et al. (2019). Genome-wide association study reveals novel genomic regions associated with high grain protein content in wheat lines derived from wild emmer wheat. Front. Plant Sci. 10, 464–464. doi: 10.3389/fpls.2019.00464
Liu, Y., Lin, Y., Gao, S., Li, Z., Ma, J., Deng, M., et al. (2017). A genome-wide association study of 23 agronomic traits in Chinese wheat landraces. Plant J. 91, 861–873. doi: 10.1111/tpj.13614
Murray, M. G., and Thompson, W. F. (1980). Rapid isolation of high molecular weight plant DNA. Nucleic Acids Res. 8, 4321–4326. doi: 10.1093/nar/8.19.4321
Myers, S. S., Zanobetti, A., Kloog, I., Huybers, P., Leakey, A. D., Bloom, A. J., et al. (2014). Increasing CO2 threatens human nutrition. Nature 510, 139–142. doi: 10.1038/nature13179
Peleg, Z., Cakmak, I., Ozturk, L., Yazici, A., Jun, Y., Budak, H., et al. (2009). Quantitative trait loci conferring grain mineral nutrient concentrations in durum wheat × wild emmer wheat RIL population. Theor. Appl. Genet. 119, 353–369. doi: 10.1007/s00122-009-1044-z
Peleg, Z., Saranga, Y., Yazici, A., Fahima, T., Ozturk, L., and Cakmak, I. (2008). Grain zinc, iron and protein concentrations and zinc-efficiency in wild emmer wheat under contrasting irrigation regimes. Plant Soil 306, 57–67. doi: 10.1007/s11104-007-9417-z
Pu, Z. E., Ma, Y. U., He, Q. Y., Chen, G. Y., Wang, J. R., Liu, Y. X., et al. (2014). Quantitative trait loci associated with micronutrient concentrations in two recombinant inbred wheat lines. J. Integr. Agr. 13, 2322–2329. doi: 10.1016/S2095-3119(13)60640-1
Roshanzamir, H., Kordenaeej, A., and Bostani, A. (2013). Mapping QTLs related to Zn and Fe concentrations in bread wheat (Triticum aestivum) grain using microsatellite markers. Iranian J. Genet. Plant. Breed. 2, 10–17.
Shi, R. L., Tong, Y. P., Jing, R. L., Zhang, F. S., and Zou, C. Q. (2013). Characterization of quantitative trait loci for grain minerals in hexaploid wheat (Triticum aestivum L.). J. Integr. Agr. 12, 1512–1521. doi: 10.1016/S2095-3119(13)60559-6
Su, J., Fan, S., Li, L., Wei, H., Wang, C., Wang, H., et al. (2016). Detection of favorable QTL alleles and candidate genes for lint percentage by GWAS in Chinese upland cotton. Front. Plant Sci. 7:1576. doi: 10.3389/fpls.2016.01576
Sukumaran, S., Dreisigacker, S., Lopes, M., Chavez, P., and Reynolds, M. P. (2015). Genome-wide association study for grain yield and related traits in an elite spring wheat population grown in temperate irrigated environments. Theor. Appl. Genet. 128, 353–363. doi: 10.1007/s00122-014-2435-3
Tabbita, F., Lewis, S., Vouilloz, J. P., Ortega, M. A., Kade, M., Abbate, P. E., et al. (2013). Effects of the Gpc-B1 locus on high grain protein content introgressed into a rgentinean wheat germplasm. Plant Breed. 132, 48–52. doi: 10.1111/pbr.12011
Tabbita, F., Pearce, S., and Barneix, A. J. (2017). Breeding for increased grain protein and micronutrient content in wheat: ten years of the GPC-B1 gene. J. Cereal Sci. 73, 183–191. doi: 10.1016/j.jcs.2017.01.003
Tadesse, W., Ogbonnaya, F. C., Jighly, A., Sanchez-Garcia, M., Sohail, Q., Rajaram, S., et al. (2015). Genome-wide association mapping of yield and grain quality traits in winter wheat genotypes. PLoS ONE 10:e0141339. doi: 10.1371/journal.pone.0141339
Uauy, C., Distelfeld, A., Fahima, T., Blechl, A., and Dubcovsky, J. (2006). A NAC gene regulating senescence improves grain protein, zinc, and iron content in wheat. Science 314, 1298–1301. doi: 10.1126/science.1133649
Upadhyaya, H. D., Bajaj, D., Das, S., Kumar, V., Gowda, C. L., Sharma, S., et al. (2016). Genetic dissection of seed-iron and zinc concentrations in chickpea. Sci. Rep. 6:24050. doi: 10.1038/srep24050
Velu, G., Tutus, Y., Gomez-Becerra, H. F., Hao, Y., Demir, L., Kara, R., et al. (2017). QTL mapping for grain zinc and iron concentrations and zinc efficiency in a tetraploid and hexaploid wheat mapping populations. Plant Soil 411, 81–99. doi: 10.1007/s11104-016-3025-8
Wang, D. (2015). The Expression Level of NAM-B1 Gene and its Effect on Grain Nutritional Quality in Derivatives between Triticum turgidum ssp. dicoccoides and Common Wheat Cultivar Chuannong 16 (in Chinese). Master dissertation, Sichuan Agricultural University, Ya'an, China.
Wang, N., Chen, B., Xu, K., Gao, G., Li, F., Qiao, J., et al. (2016). Association mapping of flowering time QTLs and insight into their contributions to rapeseed growth habits. Front. Plant Sci. 7:338. doi: 10.3389/fpls.2016.00338
Wang, S. X., Zhu, Y. L., Zhang, D. X., Shao, H., Liu, P., Hu, J. B., et al. (2017). Genome-wide association study for grain yield and related traits in elite wheat varieties and advanced lines using SNP markers. PLoS ONE 12:e0188662. doi: 10.1371/journal.pone.0188662
Wu, B., Lu, J., Li, P., Hu, X., and Zheng, Y. (2008). Genetic diversity and potential utilization of grain protein content and sink traits in Triticum dicoccoides. J. Sichuan Agric. Univ. 26, 221–225. doi: 10.16036/j.issn.1000-2650.2008.03.004
Xiao, Y., Tong, H., Yang, X., Xu, S., Pan, Q., Qiao, F., et al. (2016). Genome-wide dissection of the maize ear genetic architecture using multiple populations. New Phytol. 210, 1095–1106. doi: 10.1111/nph.13814
Xu, Y., An, D., Liu, D., Zhang, A., Xu, H., and Li, B. (2012). Molecular mapping of QTLs for grain zinc, iron and protein concentration of wheat across two environments. Field Crops Res. 138, 57–62. doi: 10.1016/j.fcr.2012.09.017
Yan, J., Xue, W. T., Yang, R. Z., Qin, H. B., Zhao, G., Tzion, F., et al. (2018). Quantitative trait loci conferring grain selenium nutrient in durum wheat × wild emmer wheat RIL population. Czech J. Genet. Plant Breed. 54, 52–58. doi: 10.17221/112/2016-CJGPB
Yu, H., and Tian, J. (2012). Association between starch pasting properties and DArT markers in common wheat. Acta Agron. Sin. 38, 1997–2006. doi: 10.3724/SP.J.1006.2012.01997
Yu, H., Xiao, J., and Tian, J. (2014). Genome-wide association analysis of flour color (Whiteness) using DArT markers in common wheat. Acta Agron. Sin. 40, 2198–2202. doi: 10.3724/SP.J.1006.2014.02198
Keywords: wild emmer, common wheat, grain micronutrient concentrations, wide hybridization, GWAS, biofortification
Citation: Liu J, Huang L, Li T, Liu Y, Yan Z, Tang G, Zheng Y, Liu D and Wu B (2021) Genome-Wide Association Study for Grain Micronutrient Concentrations in Wheat Advanced Lines Derived From Wild Emmer. Front. Plant Sci. 12:651283. doi: 10.3389/fpls.2021.651283
Received: 09 January 2021; Accepted: 08 April 2021;
Published: 14 May 2021.
Edited by:
Soren K. Rasmussen, University of Copenhagen, DenmarkReviewed by:
Hikmet Budak, Montana Bioagriculture, Inc., United StatesVelu Govindan, International Maize and Wheat Improvement Center, Mexico
Antonio Blanco, University of Bari Aldo Moro, Italy
Copyright © 2021 Liu, Huang, Li, Liu, Yan, Tang, Zheng, Liu and Wu. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Bihua Wu, wubihua2017@126.com
†These authors have contributed equally to this work