Identification of Novel Genomic Regions for Biofortification Traits Using an SNP Marker-Enriched Linkage Map in Wheat (Triticum aestivum L.)

Prospects for mineral biofortification of wheat: classical breeding and agronomy

Low intake of micro- and macroelements and vitamins in food negatively affects the health of more than two billion people around the world provoking chronic diseases. For the majority of the world's population, these are soft and durum wheats that provide beneficial nutrients, however their modern high-yielding varieties have a significantly depleted grain mineral composition that have reduced mineral intake through food. Biofortification is a new research trend, whose main goal is to improve the nutritional qualities of agricultural crops using a set of classical (hybridization and selection) methods as well and the modern ones employing gene/QTL mapping, bioinformatic analysis, transgenesis, mutagenesis and genome editing. Using the classical breeding methods, biofortified varieties have been bred as a part of various international programs funded by HarvestPlus, CIMMYT, ICARDA. Despite the promise of transgenesis and genome editing, these labor-intensive methods require significant investments, so these technologies, when applied to wheat, are still at the development stage and cannot be applied routinely. In recent years, the interest in wheat biofortification has increased due to the advances in mapping genes and QTLs for agronomically important traits. The new markers obtained from wheat genome sequencing and application of bioinformatic methods (GWAS, meta-QTL analysis) has expanded our knowledge on the traits that determine the grain mineral concentration and has identified the key gene candidates. This review describes the current research on genetic biofortification of wheat in the world and in Russia and provides information on the use of cultivated and wild-relative germplasms to expand the genetic diversity of modern wheat varieties.

Improving wheat grain composition for human health by constructing a QTL atlas for essential minerals

Wheat is an important source of minerals for human nutrition and increasing grain mineral content can contribute to reducing mineral deficiencies. Here, we identify QTLs for mineral micronutrients in grain of wheat by determining the contents of six minerals in a total of eleven sample sets of three biparental populations from crosses between A.E. Watkins landraces and cv. Paragon. Twenty-three of the QTLs are mapped in two or more sample sets, with LOD scores above five in at least one set with the increasing alleles for sixteen of the QTLs being present in the landraces and seven in Paragon. Of these QTLs, the number for each mineral varies between three and five and they are located on 14 of the 21 chromosomes, with clusters on chromosomes 5A (four), 6A (three), and 7A (three). The gene content within 5 megabases of DNA on either side of the marker for the QTL with the highest LOD score is determined and the gene responsible for the strongest QTL (chromosome 5A for Ca) identified as an ATPase transporter gene (TraesCS5A02G543300) using mutagenesis. The identification of these QTLs, together with associated SNP markers and candidate genes, will facilitate the improvement of grain nutritional quality.

Genome-wide association study of carotenoids in maize kernel

In this study, the contents of four carotenoids in 244 maize inbred lines were detected and about three million single nucleotide polymorphisms (SNPs) for genome-wide association study to preliminarily analyze the genetic mechanism of maize kernel carotenoids. We identified 826 quantitative trait loci (QTLs) were significantly associated with carotenoids contents, and two key candidate genes Zm00001d029526 (CYP18) and Zm00001d023336 (wrky91) were obtained. In addition, we found a germplasm IL78 with higher carotenoids. The results of this study can provide a theoretical basis for screening genes that guide kernel carotenoids selection breeding.

Genome-wide association study and expression of candidate genes for Fe and Zn concentration in sorghum grains

Sorghum germplasm showed grain Fe and Zn genetic variability, but a few varieties were biofortified with these minerals. This work contributes to narrowing this gap. Fe and Zn concentrations along with 55,068 high-quality GBS SNP data from 140 sorghum accessions were used in this study. Both micronutrients exhibited good variability with respective ranges of 22.09-52.55 ppm and 17.92-43.16 ppm. Significant marker-trait associations were identified on chromosomes 1, 3, and 5. Two major effect SNPs (S01_72265728 and S05_58213541) explained 35% and 32% of Fe and Zn phenotypic variance, respectively. The SNP S01_72265728 was identified in the cytochrome P450 gene and showed a positive effect on Fe accumulation in the kernel, while S05_58213541 was intergenic near Sobic.005G134800 (zinc-binding ribosomal protein) and showed negative effect on Zn. Tissue-specific in silico expression analysis resulted in higher levels of Sobic.003G350800 gene product in several tissues such as leaf, root, flower, panicle, and stem. Sobic.005G188300 and Sobic.001G463800 were expressed moderately at grain maturity and anthesis in leaf, root, panicle, and seed tissues. The candidate genes expressed in leaves, stems, and grains will be targeted to improve grain and stover quality. The haplotypes identified will be useful in forward genetics breeding.

Multivariate Genome-Wide Association Study of Concentrations of Seven Elements in Seeds Reveals Four New Loci in Russian Wheat Lines

Wheat is a cereal grain that plays an important role in the world's food industry. The identification of the loci that change the concentration of elements in wheat seeds is an important challenge nowadays especially for genomic selection and breeding of novel varieties. In this study, we performed a multivariate genome-wide association study (GWAS) of the seven traits-concentrations of Zn, Mg, Mn, Ca, Cu, Fe, and K in grain-of the Russian collection of common wheat Triticum aestivum (N = 149 measured in two years in two different fields). We replicated one known locus associated with the concentration of Zn (IAAV1375). We identified four novel loci-BS00022069_51 (associated with concentrations of Ca and K), RFL_Contig6053_3082 (associated with concentrations of Fe and Mn), Kukri_rep_c70864_329 (associated with concentrations of all elements), and IAAV8416 (associated with concentrations of Fe and Mn)-three of them were located near the genes TraesCS6A02G375400, TraesCS7A02G094800, and TraesCS5B02G325400. Our result adds novel information on the loci involved in wheat grain element contents and may be further used in genomic selection.

Consensus genomic regions for grain quality traits in wheat revealed by Meta-QTL analysis and in silico transcriptome integration

Grain quality traits are the key factors that determine the economic value of wheat and are largely influenced by genetics and the environment. In this study, using a meta-analysis of quantitative trait loci (QTLs) and a comprehensive in silico transcriptome assessment, we identified key genomic regions and putative candidate genes for the grain quality traits protein content, gluten content, and test weight. A total of 508 original QTLs were collected from 41 articles on QTL mapping for the three quality traits in wheat published from 2003 to 2021. When these original QTLs were projected onto a high-density consensus map consisting of 14,548 markers, 313 QTLs resulted in the identification of 64 MQTLs distributed across 17 of the 21 chromosomes. Most of the meta-QTLs (MQTLs) were distributed on sub-genomes A and B. Compared with the original QTLs, the confidence interval (CI) of the MQTLs was smaller, with an average CI of 4.47 cM, while the projected QTLs CI was 11.13 cM (2.49-fold lower). The corresponding physical length of the MQTL ranged from 0.45 to 239.01 Mb. Thirty-one of these 64 MQTLs were validated in at least one genome-wide association study. In addition, five of the 64 MQTLs were selected and designated as core MQTLs. The 211 quality-related genes from rice were used to identify wheat homologs in MQTLs. In combination with transcriptional and omics analyses, 135 putative candidate genes were identified from 64 MQTL regions. The findings should contribute to a better understanding of the molecular genetic mechanisms underlying grain quality and the improvement of these traits in wheat breeding.

Nutritional Genomic Approach for Improving Grain Protein Content in Wheat

Grain protein content (GPC) is a key aspect of grain quality, a major determinant of the flour functional properties and grain nutritional value of bread wheat. Exploiting diverse germplasms to identify genes for improving crop performance and grain nutritional quality is needed to enhance food security. Here, we evaluated GPC in a panel of 255 Triticum aestivum L. accessions from 27 countries. GPC determined in seeds from three consecutive crop seasons varied from 8.6 to 16.4% (11.3% on average). Significant natural phenotypic variation in GPC among genotypes and seasons was detected. The population was evaluated for the presence of the trait-linked single nucleotide polymorphism (SNP) markers via a genome-wide association study (GWAS). GWAS analysis conducted with calculated best linear unbiased estimates (BLUEs) of phenotypic data and 90 K SNP array using the fixed and random model circulating probability unification (FarmCPU) model identified seven significant genomic regions harboring GPC-associated markers on chromosomes 1D, 3A, 3B, 3D, 4B and 5A, of which those on 3A and 3B shared associated SNPs with at least one crop season. The verified SNP-GPC associations provide new promising genomic signals on 3A (SNPs: Excalibur_c13709_2568 and wsnp_Ku_c7811_13387117) and 3B (SNP: BS00062734_51) underlying protein improvement in wheat. Based on the linkage disequilibrium for significant SNPs, the most relevant candidate genes within a 4 Mbp-window included genes encoding a subtilisin-like serine protease; amino acid transporters; transcription factors; proteins with post-translational regulatory functions; metabolic proteins involved in the starch, cellulose and fatty acid biosynthesis; protective and structural proteins, and proteins associated with metal ions transport or homeostasis. The availability of molecular markers within or adjacent to the sequences of the detected candidate genes might assist a breeding strategy based on functional markers to improve genetic gains for GPC and nutritional quality in wheat.

Genome-Wide Association Study for Grain Protein, Thousand Kernel Weight, and Normalized Difference Vegetation Index in Bread Wheat (Triticum aestivum L.)

Genomic regions governing grain protein content (GPC), 1000 kernel weight (TKW), and normalized difference vegetation index (NDVI) were studied in a set of 280 bread wheat genotypes. The genome-wide association (GWAS) panel was genotyped using a 35K Axiom array and phenotyped in three environments. A total of 26 marker-trait associations (MTAs) were detected on 18 chromosomes covering the A, B, and D subgenomes of bread wheat. The GPC showed the maximum MTAs (16), followed by NDVI (6), and TKW (4). A maximum of 10 MTAs was located on the B subgenome, whereas, 8 MTAs each were mapped on the A and D subgenomes. In silico analysis suggest that the SNPs were located on important putative candidate genes such as NAC domain superfamily, zinc finger RING-H2-type, aspartic peptidase domain, folylpolyglutamate synthase, serine/threonine-protein kinase LRK10, pentatricopeptide repeat, protein kinase-like domain superfamily, cytochrome P450, and expansin. These candidate genes were found to have different roles including regulation of stress tolerance, nutrient remobilization, protein accumulation, nitrogen utilization, photosynthesis, grain filling, mitochondrial function, and kernel development. The effects of newly identified MTAs will be validated in different genetic backgrounds for further utilization in marker-aided breeding.

Mapping of the QTLs governing grain micronutrients and thousand kernel weight in wheat (Triticum aestivum L.) using high density SNP markers

Biofortification is gaining importance globally to improve human nutrition through enhancing the micronutrient content, such as vitamin A, iron, and zinc, in staple food crops. The present study aims to identify the chromosomal regions governing the grain iron concentration (GFeC), grain zinc concentration (GZnC), and thousand kernel weight (TKW) using recombinant inbred lines (RILs) in wheat, developed from a cross between HD3086 and HI1500. The experiment was conducted in four different production conditions at Delhi viz., control, drought, heat, and combined heat and drought stress and at Indore under drought stress. Grain iron and zinc content increased under heat and combined stress conditions, while thousand kernel weight decreased. Medium to high heritability with a moderate correlation between grain iron and zinc was observed. Out of 4,106 polymorphic markers between the parents, 3,407 SNP markers were used for linkage map construction which spanned over a length of 14791.18 cm. QTL analysis identified a total of 32 chromosomal regions governing the traits under study, which includes 9, 11, and 12 QTLs for GFeC, GZnC, and TKW, respectively. A QTL hotspot was identified on chromosome 4B which is associated with grain iron, grain zinc, and thousand kernel weight explaining the phenotypic variance of 29.28, 10.98, and 17.53%, respectively. Similarly, common loci were identified on chromosomes 4B and 4D for grain iron, zinc, and thousand kernel weight. In silico analysis of these chromosomal regions identified putative candidate genes that code for proteins such as Inositol 1,3,4-trisphosphate 5/6-kinase, P-loop containing nucleoside triphosphate hydrolase, Pleckstrin homology (PH) domains, Serine-threonine/tyrosine-protein kinase and F-box-like domain superfamily proteins which play role in many important biochemical or physiological process. The identified markers linked to QTLs can be used in MAS once successfully validated.

Molecular Mapping of Biofortification Traits in Bread Wheat (Triticum aestivum L.) Using a High-Density SNP Based Linkage Map

A set of 188 recombinant inbred lines (RILs) derived from a cross between a high-yielding Indian bread wheat cultivar HD2932 and a synthetic hexaploid wheat (SHW) Synthetic 46 derived from tetraploid Triticum turgidum (AA, BB 2n = 28) and diploid Triticum tauschii (DD, 2n = 14) was used to identify novel genomic regions associated in the expression of grain iron concentration (GFeC), grain zinc concentration (GZnC), grain protein content (GPC) and thousand kernel weight (TKW). The RIL population was genotyped using SNPs from 35K Axiom® Wheat Breeder's Array and 34 SSRs and phenotyped in two environments. A total of nine QTLs including five for GPC (QGpc.iari_1B, QGpc.iari_4A, QGpc.iari_4B, QGpc.iari_5D, and QGpc.iari_6B), two for GFeC (QGfec.iari_5B and QGfec.iari_6B), and one each for GZnC (QGznc.iari_7A) and TKW (QTkw.iari_4B) were identified. A total of two stable and co-localized QTLs (QGpc.iari_4B and QTkw.iari_4B) were identified on the 4B chromosome between the flanking region of Xgwm149-AX-94559916. In silico analysis revealed that the key putative candidate genes such as P-loop containing nucleoside triphosphatehydrolase, Nodulin-like protein, NAC domain, Purine permease, Zinc-binding ribosomal protein, Cytochrome P450, Protein phosphatase 2A, Zinc finger CCCH-type, and Kinesin motor domain were located within the identified QTL regions and these putative genes are involved in the regulation of iron homeostasis, zinc transportation, Fe, Zn, and protein remobilization to the developing grain, regulation of grain size and shape, and increased nitrogen use efficiency. The identified novel QTLs, particularly stable and co-localized QTLs are useful for subsequent use in marker-assisted selection (MAS).

Mapping QTL for Phenological and Grain-Related Traits in a Mapping Population Derived from High-Zinc-Biofortified Wheat

Genomic regions governing days to heading (DH), days to maturity (DM), plant height (PH), thousand-kernel weight (TKW), and test weight (TW) were investigated in a set of 190 RILs derived from a cross between a widely cultivated wheat-variety, Kachu (DPW-621-50), and a high-zinc variety, Zinc-Shakti. The RIL population was genotyped using 909 DArTseq markers and phenotyped in three environments. The constructed genetic map had a total genetic length of 4665 cM, with an average marker density of 5.13 cM. A total of thirty-seven novel quantitative trait loci (QTL), including twelve for PH, six for DH, five for DM, eight for TKW and six for TW were identified. A set of 20 stable QTLs associated with the expression of DH, DM, PH, TKW, and TW were identified in two or more environments. Three novel pleiotropic genomic-regions harboring co-localized QTLs governing two or more traits were also identified. In silico analysis revealed that the DArTseq markers were located on important putative candidate genes such as MLO-like protein, Phytochrome, Zinc finger and RING-type, Cytochrome P450 and pentatricopeptide repeat, involved in the regulation of pollen maturity, the photoperiodic modulation of flowering-time, abiotic-stress tolerance, grain-filling duration, thousand-kernel weight, seed morphology, and plant growth and development. The identified novel QTLs, particularly stable and co-localized QTLs, will be validated to estimate their effects in different genetic backgrounds for subsequent use in marker-assisted selection (MAS).

Ensuring Nutritional Security in India through Wheat Biofortification: A Review

Undernourishment of nutrients, also known as hidden hunger, affects over 2 billion populace globally. Even though stunting among children below five years of age has decreased in India in the last ten years, India is home to roughly thirty percent of the world's population of stunted pre-schoolers. A significant improvement has been witnessed in the targeted development and deployment of biofortified crops; approximately 20 million farm households from developing counties benefit from cultivating and consuming biofortified crops. There is ample scope for including biofortified varieties in the seed chain, ensuring nutritional security. Wheat is a dietary staple in India, typically consumed as wholemeal flour in the form of flatbreads such as chapatti and roti. Wheat contributes to nearly one fifth of global energy requirements and can also provide better amounts of iron (Fe) and zinc (Zn). As a result, biofortified wheat can serve as a medium for delivery of essential micronutrients such as Fe and Zn to end users. This review discusses wheat biofortification components such as Fe and Zn dynamics, its uptake and movement in plants, the genetics of their buildup, and the inclusion of biofortified wheat varieties in the seed multiplication chain concerning India.

Genetic dissection of grain iron and zinc, and thousand kernel weight in wheat (Triticum aestivum L.) using genome-wide association study

Genetic biofortification is recognized as a cost-effective and sustainable strategy to reduce micronutrient malnutrition. Genomic regions governing grain iron concentration (GFeC), grain zinc concentration (GZnC), and thousand kernel weight (TKW) were investigated in a set of 280 diverse bread wheat genotypes. The genome-wide association (GWAS) panel was genotyped using 35 K Axiom Array and phenotyped in five environments. The GWAS analysis showed a total of 17 Bonferroni-corrected marker-trait associations (MTAs) in nine chromosomes representing all the three wheat subgenomes. The TKW showed the highest MTAs (7), followed by GZnC (5) and GFeC (5). Furthermore, 14 MTAs were identified with more than 10% phenotypic variation. One stable MTA i.e. AX-95025823 was identified for TKW in both E4 and E5 environments along with pooled data, which is located at 68.9 Mb on 6A chromosome. In silico analysis revealed that the SNPs were located on important putative candidate genes such as Multi antimicrobial extrusion protein, F-box domain, Late embryogenesis abundant protein, LEA-18, Leucine-rich repeat domain superfamily, and C3H4 type zinc finger protein, involved in iron translocation, iron and zinc homeostasis, and grain size modifications. The identified novel MTAs will be validated to estimate their effects in different genetic backgrounds for subsequent use in marker-assisted selection. The identified SNPs will be valuable in the rapid development of biofortified wheat varieties to ameliorate the malnutrition problems.

Wheat Biofortification: Utilizing Natural Genetic Diversity, Genome-Wide Association Mapping, Genomic Selection, and Genome Editing Technologies

Alleviating micronutrients associated problems in children below five years and women of childbearing age, remains a significant challenge, especially in resource-poor nations. One of the most important staple food crops, wheat attracts the highest global research priority for micronutrient (Fe, Zn, Se, and Ca) biofortification. Wild relatives and cultivated species of wheat possess significant natural genetic variability for these micronutrients, which has successfully been utilized for breeding micronutrient dense wheat varieties. This has enabled the release of 40 biofortified wheat cultivars for commercial cultivation in different countries, including India, Bangladesh, Pakistan, Bolivia, Mexico and Nepal. In this review, we have systematically analyzed the current understanding of availability and utilization of natural genetic variations for grain micronutrients among cultivated and wild relatives, QTLs/genes and different genomic regions regulating the accumulation of micronutrients, and the status of micronutrient biofortified wheat varieties released for commercial cultivation across the globe. In addition, we have also discussed the potential implications of emerging technologies such as genome editing to improve the micronutrient content and their bioavailability in wheat.

Genome-Wide Association Mapping Identifies Key Genomic Regions for Grain Zinc and Iron Biofortification in Bread Wheat

Accelerating breeding efforts for developing biofortified bread wheat varieties necessitates understanding the genetic control of grain zinc concentration (GZnC) and grain iron concentration (GFeC). Hence, the major objective of this study was to perform genome-wide association mapping to identify consistently significant genotyping-by-sequencing markers associated with GZnC and GFeC using a large panel of 5,585 breeding lines from the International Maize and Wheat Improvement Center. These lines were grown between 2018 and 2021 in an optimally irrigated environment at Obregon, Mexico, while some of them were also grown in a water-limiting drought-stressed environment and a space-limiting small plot environment and evaluated for GZnC and GFeC. The lines showed a large and continuous variation for GZnC ranging from 27 to 74.5 ppm and GFeC ranging from 27 to 53.4 ppm. We performed 742,113 marker-traits association tests in 73 datasets and identified 141 markers consistently associated with GZnC and GFeC in three or more datasets, which were located on all wheat chromosomes except 3A and 7D. Among them, 29 markers were associated with both GZnC and GFeC, indicating a shared genetic basis for these micronutrients and the possibility of simultaneously improving both. In addition, several significant GZnC and GFeC associated markers were common across the irrigated, water-limiting drought-stressed, and space-limiting small plots environments, thereby indicating the feasibility of indirect selection for these micronutrients in either of these environments. Moreover, the many significant markers identified had minor effects on GZnC and GFeC, suggesting a quantitative genetic control of these traits. Our findings provide important insights into the complex genetic basis of GZnC and GFeC in bread wheat while implying limited prospects for marker-assisted selection and the need for using genomic selection.

Genome-wide association study identifies loci and candidate genes for grain micronutrients and quality traits in wheat (Triticum aestivum L.)

Malnutrition due to micronutrients and protein deficiency is recognized among the major global health issues. Genetic biofortification of wheat is a cost-effective and sustainable strategy to mitigate the global micronutrient and protein malnutrition. Genomic regions governing grain zinc concentration (GZnC), grain iron concentration (GFeC), grain protein content (GPC), test weight (TW), and thousand kernel weight (TKW) were investigated in a set of 184 diverse bread wheat genotypes through genome-wide association study (GWAS). The GWAS panel was genotyped using Breeders' 35 K Axiom Array and phenotyped in three different environments during 2019-2020. A total of 55 marker-trait associations (MTAs) were identified representing all three sub-genomes of wheat. The highest number of MTAs were identified for GPC (23), followed by TKW (15), TW (11), GFeC (4), and GZnC (2). Further, a stable SNP was identified for TKW, and also pleiotropic regions were identified for GPC and TKW. In silico analysis revealed important putative candidate genes underlying the identified genomic regions such as F-box-like domain superfamily, Zinc finger CCCH-type proteins, Serine-threonine/tyrosine-protein kinase, Histone deacetylase domain superfamily, and SANT/Myb domain superfamily proteins, etc. The identified novel MTAs will be validated to estimate their effects in different genetic backgrounds for subsequent use in marker-assisted selection.

Applications of Genomic Tools in Plant Breeding: Crop Biofortification

Crop breeding has mainly been focused on increasing productivity, either directly or by decreasing the losses caused by biotic and abiotic stresses (that is, incorporating resistance to diseases and enhancing tolerance to adverse conditions, respectively). Quite the opposite, little attention has been paid to improve the nutritional value of crops. It has not been until recently that crop biofortification has become an objective within breeding programs, through either conventional methods or genetic engineering. There are many steps along this long path, from the initial evaluation of germplasm for the content of nutrients and health-promoting compounds to the development of biofortified varieties, with the available and future genomic tools assisting scientists and breeders in reaching their objectives as well as speeding up the process. This review offers a compendium of the genomic technologies used to explore and create biodiversity, to associate the traits of interest to the genome, and to transfer the genomic regions responsible for the desirable characteristics into potential new varieties. Finally, a glimpse of future perspectives and challenges in this emerging area is offered by taking the present scenario and the slow progress of the regulatory framework as the starting point.

Harnessing the Wild Relatives and Landraces for Fe and Zn Biofortification in Wheat through Genetic Interventions—A Review

Micronutrient deficiencies, particularly iron (Fe) and zinc (Zn), in human diets are affecting over three billion people globally, especially in developing nations where diet is cereal-based. Wheat is one of several important cereal crops that provide food calories to nearly one-third of the population of the world. However, the bioavailability of Zn and Fe in wheat is inherently low, especially under Zn deficient soils. Although various fortification approaches are available, biofortification, i.e., development of mineral-enriched cultivars, is an efficient and sustainable approach to alleviate malnutrition. There is enormous variability in Fe and Zn in wheat germplasm, especially in wild relatives, but this is not utilized to the full extent. Grain Fe and Zn are quantitatively inherited, but high-heritability and genetic correlation at multiple locations indicate the high stability of Fe and Zn in wheat. In the last decade, pre-breeding activities have explored the potential of wild relatives to develop Fe and Zn rich wheat varieties. Furthermore, recent advances in molecular biology have improved the understanding of the uptake, storage, and bioavailability of Fe and Zn. Various transportation proteins encoding genes like YSL 2, IRT 1, OsNAS 3, VIT 1, and VIT 2 have been identified for Fe and Zn uptake, transfer, and accumulation at different developing stages. Hence, the availability of major genomic regions for Fe and Zn content and genome editing technologies are likely to result in high-yielding Fe and Zn biofortified wheat varieties. This review covers the importance of wheat wild relatives for Fe and Zn biofortification, progress in genomics-assisted breeding, and transgenic breeding for improving Fe and Zn content in wheat.