Multi-modal deep learning improves grain yield prediction in wheat breeding by fusing genomics and phenomics.

MOTIVATION: Developing new crop varieties with superior performance is highly important to ensure robust and sustainable global food security. The speed of variety development is limited by long field cycles and advanced generation selections in plant breeding programs. While methods to predict yield from genotype or phenotype data have been proposed, improved performance and integrated models are needed. RESULTS: We propose a machine learning model that leverages both genotype and phenotype measurements by fusing genetic variants with multiple data sources collected by unmanned aerial systems. We use a deep multiple instance learning framework with an attention mechanism that sheds light on the importance given to each input during prediction, enhancing interpretability. Our model reaches 0.754 ± 0.024 Pearson correlation coefficient when predicting yield in similar environmental conditions; a 34.8% improvement over the genotype-only linear baseline (0.559 ± 0.050). We further predict yield on new lines in an unseen environment using only genotypes, obtaining a prediction accuracy of 0.386 ± 0.010, a 13.5% improvement over the linear baseline. Our multi-modal deep learning architecture efficiently accounts for plant health and environment, distilling the genetic contribution and providing excellent predictions. Yield prediction algorithms leveraging phenotypic observations during training therefore promise to improve breeding programs, ultimately speeding up delivery of improved varieties. AVAILABILITY AND IMPLEMENTATION: Available at https://github.com/BorgwardtLab/PheGeMIL (code) and https://doi.org/doi:10.5061/dryad.kprr4xh5p (data).

GWAS elucidated grain yield genetics in Indian spring wheat under diverse water conditions.

KEY MESSAGE: Underpinned natural variations and key genes associated with yield under different water regimes, and identified genomic signatures of genetic gain in the Indian wheat breeding program. A novel KASP marker for TKW under water stress was developed and validated. A comprehensive genome-wide association study was conducted on 300 spring wheat genotypes to elucidate the natural variations associated with grain yield and its eleven contributing traits under fully irrigated, restricted water, and simulated no water conditions. Utilizing the 35K Wheat Breeders' Array, we identified 1155 quantitative trait nucleotides (QTNs), with 207 QTNs exhibiting stability across diverse conditions. These QTNs were further delimited into 539 genomic regions using a genome-wide LD value of 3.0 Mbp, revealing pleiotropic control across traits and conditions. Sub-genome A was significantly associated with traits under irrigated conditions, while sub-genome B showed more QTNs under water stressed conditions. Favourable alleles with significantly associated QTNs were delineated, with a notable pyramiding effect for enhancing trait performance. Additionally, allele of only 921 QTNs significantly affected the population mean. Allele profiling highlighted C-306 as a most potential source of drought tolerance. Moreover, 762 genes overlapping significant QTNs were identified, narrowing down to 27 putative candidate genes overlapping 29 novel and functional SNPs expressing (≥ 0.5 tpm) relevance across various growth conditions. A new KASP assay was developed, targeting a gene TraesCS2A03G1123700 regulating thousand kernel weight under severe drought condition. Genomic selection models (GBLUP, BayesB, MxE, and R-Norm) demonstrated an average prediction accuracy of 0.06-0.58 across environments, indicating potential for trait selection. Retrospective analysis of the Indian wheat breeding program supported a genetic gain in GY at the rate of ca. 0.56% per breeding cycle, since 1960, supporting the identification of genomic signatures driving trait selection and genetic gain. These findings offer insight into improving the rate of genetic gain in wheat breeding programs globally.

Genomic selection for spot blotch in bread wheat breeding panels, full-sibs and half-sibs and index-based selection for spot blotch, heading and plant height.

KEY MESSAGE: Genomic selection is a promising tool to select for spot blotch resistance and index-based selection can simultaneously select for spot blotch resistance, heading and plant height. A major biotic stress challenging bread wheat production in regions characterized by humid and warm weather is spot blotch caused by the fungus Bipolaris sorokiniana. Since genomic selection (GS) is a promising selection tool, we evaluated its potential for spot blotch in seven breeding panels comprising 6736 advanced lines from the International Maize and Wheat Improvement Center. Our results indicated moderately high mean genomic prediction accuracies of 0.53 and 0.40 within and across breeding panels, respectively which were on average 177.6% and 60.4% higher than the mean accuracies from fixed effects models using selected spot blotch loci. Genomic prediction was also evaluated in full-sibs and half-sibs panels and sibs were predicted with the highest mean accuracy (0.63) from a composite training population with random full-sibs and half-sibs. The mean accuracies when full-sibs were predicted from other full-sibs within families and when full-sibs panels were predicted from other half-sibs panels were 0.47 and 0.44, respectively. Comparison of GS with phenotypic selection (PS) of the top 10% of resistant lines suggested that GS could be an ideal tool to discard susceptible lines, as greater than 90% of the susceptible lines discarded by PS were also discarded by GS. We have also reported the evaluation of selection indices to simultaneously select non-late and non-tall genotypes with low spot blotch phenotypic values and genomic-estimated breeding values. Overall, this study demonstrates the potential of integrating GS and index-based selection for improving spot blotch resistance in bread wheat.

Harnessing translational research in wheat for climate resilience.

Despite being the world's most widely grown crop, research investments in wheat (Triticum aestivum and Triticum durum) fall behind those in other staple crops. Current yield gains will not meet 2050 needs, and climate stresses compound this challenge. However, there is good evidence that heat and drought resilience can be boosted through translating promising ideas into novel breeding technologies using powerful new tools in genetics and remote sensing, for example. Such technologies can also be applied to identify climate resilience traits from among the vast and largely untapped reserve of wheat genetic resources in collections worldwide. This review describes multi-pronged research opportunities at the focus of the Heat and Drought Wheat Improvement Consortium (coordinated by CIMMYT), which together create a pipeline to boost heat and drought resilience, specifically: improving crop design targets using big data approaches; developing phenomic tools for field-based screening and research; applying genomic technologies to elucidate the bases of climate resilience traits; and applying these outputs in developing next-generation breeding methods. The global impact of these outputs will be validated through the International Wheat Improvement Network, a global germplasm development and testing system that contributes key productivity traits to approximately half of the global wheat-growing area.

The Aegilops ventricosa 2NvS segment in bread wheat: cytology, genomics and breeding.

KEY MESSAGE: The first cytological characterization of the 2NvS segment in hexaploid wheat; complete de novo assembly and annotation of 2NvS segment; 2NvS frequency is increasing 2NvS and is associated with higher yield. The Aegilops ventricosa 2NvS translocation segment has been utilized in breeding disease-resistant wheat crops since the early 1990s. This segment is known to possess several important resistance genes against multiple wheat diseases including root knot nematode, stripe rust, leaf rust and stem rust. More recently, this segment has been associated with resistance to wheat blast, an emerging and devastating wheat disease in South America and Asia. To date, full characterization of the segment including its size, gene content and its association with grain yield is lacking. Here, we present a complete cytological and physical characterization of this agronomically important translocation in bread wheat. We de novo assembled the 2NvS segment in two wheat varieties, 'Jagger' and 'CDC Stanley,' and delineated the segment to be approximately 33 Mb. A total of 535 high-confidence genes were annotated within the 2NvS region, with > 10% belonging to the nucleotide-binding leucine-rich repeat (NLR) gene families. Identification of groups of NLR genes that are potentially N genome-specific and expressed in specific tissues can fast-track testing of candidate genes playing roles in various disease resistances. We also show the increasing frequency of 2NvS among spring and winter wheat breeding programs over two and a half decades, and the positive impact of 2NvS on wheat grain yield based on historical datasets. The significance of the 2NvS segment in wheat breeding due to resistance to multiple diseases and a positive impact on yield highlights the importance of understanding and characterizing the wheat pan-genome for better insights into molecular breeding for wheat improvement.

Grain yield genetic gains and changes in physiological related traits for CIMMYT's High Rainfall Wheat Screening Nursery tested across international environments.

The effects of climate change together with the projected future demand represents a huge challenge for wheat production systems worldwide. Wheat breeding can contribute to global food security through the creation of genotypes exhibiting stress tolerance and higher yield potential. The objectives of our study were to (i) estimate the annual grain yield (GY) genetic gain of High Rainfall Wheat Yield Trials (HRWYT) grown from 2007 (15th HRWYT) to 2016 (24th HRWYT) across international environments, and (ii) determine the changes in physiological traits associated with GY genetic improvement. The GY genetic gains were estimated as genetic progress per se (GYP) and in terms of local checks (GYLC). In total, 239 international locations were classified into two groups: high- and low-rainfall environments based on climate variables and trial management practices. In the high-rainfall environment, the annual genetic gains for GYP and GYLC were 3.8 and 1.17 % (160 and 65.1 kg ha-1 yr-1), respectively. In the low-rainfall environment, the genetic gains were 0.93 and 0.73 % (40 and 33.1 kg ha-1 yr-1), for GYP and GYLC respectively. The GY of the lines included in each nursery showed a significant phenotypic correlation between high- and low-rainfall environments in all the examined years and several of the five best performing lines were common in both environments. The GY progress was mainly associated with increased grain weight (R2 = 0.35 p < 0.001), days to maturity (R2 = 0.20, p < 0.001) and grain filling period (R2 = 0.06, p < 0.05). These results indicate continuous GY genetic progress and yield stability in the HRWYT germplasm developed and distributed by CIMMYT.

High-throughput phenotyping platforms enhance genomic selection for wheat grain yield across populations and cycles in early stage.

Genomic selection (GS) models have been validated for many quantitative traits in wheat (Triticum aestivum L.) breeding. However, those models are mostly constrained within the same growing cycle and the extension of GS to the case of across cycles has been a challenge, mainly due to the low predictive accuracy resulting from two factors: reduced genetic relationships between different families and augmented environmental variances between cycles. Using the data collected from diverse field conditions at the International Wheat and Maize Improvement Center, we evaluated GS for grain yield in three elite yield trials across three wheat growing cycles. The objective of this project was to employ the secondary traits, canopy temperature, and green normalized difference vegetation index, which are closely associated with grain yield from high-throughput phenotyping platforms, to improve prediction accuracy for grain yield. The ability to predict grain yield was evaluated reciprocally across three cycles with or without secondary traits. Our results indicate that prediction accuracy increased by an average of 146% for grain yield across cycles with secondary traits. In addition, our results suggest that secondary traits phenotyped during wheat heading and early grain filling stages were optimal for enhancing the prediction accuracy for grain yield.

Integrating genomic-enabled prediction and high-throughput phenotyping in breeding for climate-resilient bread wheat.

Genomic selection and high-throughput phenotyping (HTP) are promising tools to accelerate breeding gains for high-yielding and climate-resilient wheat varieties. Hence, our objective was to evaluate them for predicting grain yield (GY) in drought-stressed (DS) and late-sown heat-stressed (HS) environments of the International maize and wheat improvement center's elite yield trial nurseries. We observed that the average genomic prediction accuracies using fivefold cross-validations were 0.50 and 0.51 in the DS and HS environments, respectively. However, when a different nursery/year was used to predict another nursery/year, the average genomic prediction accuracies in the DS and HS environments decreased to 0.18 and 0.23, respectively. While genomic predictions clearly outperformed pedigree-based predictions across nurseries, they were similar to pedigree-based predictions within nurseries due to small family sizes. In populations with some full-sibs in the training population, the genomic and pedigree-based prediction accuracies were on average 0.27 and 0.35 higher than the accuracies in populations with only one progeny per cross, indicating the importance of genetic relatedness between the training and validation populations for good predictions. We also evaluated the item-based collaborative filtering approach for multivariate prediction of GY using the green normalized difference vegetation index from HTP. This approach proved to be the best strategy for across-nursery predictions, with average accuracies of 0.56 and 0.62 in the DS and HS environments, respectively. We conclude that GY is a challenging trait for across-year predictions, but GS and HTP can be integrated in increasing the size of populations screened and evaluating unphenotyped large nurseries for stress-resilience within years.

Genome Wide Association Study of Karnal Bunt Resistance in a Wheat Germplasm Collection from Afghanistan.

Karnal bunt disease of wheat, caused by the fungus Neovossia indica, is one of the most important challenges to the grain industry as it affects the grain quality and also restricts the international movement of infected grain. It is a seed-, soil- and airborne disease with limited effect of chemical control. Currently, this disease is contained through the deployment of host resistance but further improvement is limited as only a few genotypes have been found to carry partial resistance. To identify genomic regions responsible for resistance in a set of 339 wheat accessions, genome-wide association study (GWAS) was undertaken using the DArTSeq® technology, in which 18 genomic regions for Karnal bunt resistance were identified, explaining 5-20% of the phenotypic variation. The identified quantitative trait loci (QTL) on chromosome 2BL showed consistently significant effects across all four experiments, whereas another QTL on 5BL was significant in three experiments. Additional QTLs were mapped on chromosomes 1DL, 2DL, 4AL, 5AS, 6BL, 6BS, 7BS and 7DL that have not been mapped previously, and on chromosomes 4B, 5AL, 5BL and 6BS, which have been reported in previous studies. Germplasm with less than 1% Karnal bunt infection have been identified and can be used for resistance breeding. The SNP markers linked to the genomic regions conferring resistance to Karnal bunt could be used to improve Karnal bunt resistance through marker-assisted selection.

Genome-wide association mapping for resistance to leaf rust, stripe rust and tan spot in wheat reveals potential candidate genes.

KEY MESSAGE: Genome-wide association mapping in conjunction with population sequencing map and Ensembl plants was used to identify markers/candidate genes linked to leaf rust, stripe rust and tan spot resistance in wheat. Leaf rust (LR), stripe rust (YR) and tan spot (TS) are some of the important foliar diseases in wheat (Triticum aestivum L.). To identify candidate resistance genes for these diseases in CIMMYT's (International Maize and Wheat Improvement Center) International bread wheat screening nurseries, we used genome-wide association studies (GWAS) in conjunction with information from the population sequencing map and Ensembl plants. Wheat entries were genotyped using genotyping-by-sequencing and phenotyped in replicated trials. Using a mixed linear model, we observed that seedling resistance to LR was associated with 12 markers on chromosomes 1DS, 2AS, 2BL, 3B, 4AL, 6AS and 6AL, and seedling resistance to TS was associated with 14 markers on chromosomes 1AS, 2AL, 2BL, 3AS, 3AL, 3B, 6AS and 6AL. Seedling and adult plant resistance (APR) to YR were associated with several markers at the distal end of chromosome 2AS. In addition, YR APR was also associated with markers on chromosomes 2DL, 3B and 7DS. The potential candidate genes for these diseases included several resistance genes, receptor-like serine/threonine-protein kinases and defense-related enzymes. However, extensive LD in wheat that decays at about 5 × 107 bps, poses a huge challenge for delineating candidate gene intervals and candidates should be further mapped, functionally characterized and validated. We also explored a segment on chromosome 2AS associated with multiple disease resistance and identified seventeen disease resistance linked genes. We conclude that identifying candidate genes linked to significant markers in GWAS is feasible in wheat, thus creating opportunities for accelerating molecular breeding.

Resistance to Spot Blotch in Two Mapping Populations of Common Wheat Is Controlled by Multiple QTL of Minor Effects.

Spot blotch (SB) is an important fungal disease of wheat in South Asia and South America. Host resistance is regarded as an economical and environmentally friendly approach of controlling SB, and the inheritance of resistance is mostly quantitative. In order to gain a better understanding on the SB resistance mechanism in CIMMYT germplasm, two bi-parental mapping populations were generated, both comprising 232 F2:7 progenies. Elite CIMMYT breeding lines, BARTAI and WUYA, were used as resistant parents, whereas CIANO T79 was used as susceptible parent in both populations. The two populations were evaluated for field SB resistance at CIMMYT's Agua Fria station for three consecutive years, from the 2012⁻2013 to 2014⁻2015 cropping seasons. Phenological traits like plant height (PH) and days to heading (DH) were also determined. Genotyping was performed using the DArTSeq genotyping-by-sequencing (GBS) platform, and a few D-genome specific SNPs and those for phenological traits were integrated for analysis. The most prominent quantitative trait locus (QTL) in both populations was found on chromosome 5AL at the Vrn-A1 locus, explaining phenotypic variations of 7⁻27%. Minor QTL were found on chromosomes 1B, 3A, 3B, 4B, 4D, 5B and 6D in BARTAI and on chromosomes 1B, 2A, 2D and 4B in WUYA, whereas minor QTL contributed by CIANO T79 were identified on chromosome 1B, 1D, 3A, 4B and 7A. In summary, resistance to SB in the two mapping populations was controlled by multiple minor QTL, with strong influence from Vrn-A1.

Genomic and pedigree-based prediction for leaf, stem, and stripe rust resistance in wheat.

KEY MESSAGE: Genomic prediction for seedling and adult plant resistance to wheat rusts was compared to prediction using few markers as fixed effects in a least-squares approach and pedigree-based prediction. The unceasing plant-pathogen arms race and ephemeral nature of some rust resistance genes have been challenging for wheat (Triticum aestivum L.) breeding programs and farmers. Hence, it is important to devise strategies for effective evaluation and exploitation of quantitative rust resistance. One promising approach that could accelerate gain from selection for rust resistance is 'genomic selection' which utilizes dense genome-wide markers to estimate the breeding values (BVs) for quantitative traits. Our objective was to compare three genomic prediction models including genomic best linear unbiased prediction (GBLUP), GBLUP A that was GBLUP with selected loci as fixed effects and reproducing kernel Hilbert spaces-markers (RKHS-M) with least-squares (LS) approach, RKHS-pedigree (RKHS-P), and RKHS markers and pedigree (RKHS-MP) to determine the BVs for seedling and/or adult plant resistance (APR) to leaf rust (LR), stem rust (SR), and stripe rust (YR). The 333 lines in the 45th IBWSN and the 313 lines in the 46th IBWSN were genotyped using genotyping-by-sequencing and phenotyped in replicated trials. The mean prediction accuracies ranged from 0.31-0.74 for LR seedling, 0.12-0.56 for LR APR, 0.31-0.65 for SR APR, 0.70-0.78 for YR seedling, and 0.34-0.71 for YR APR. For most datasets, the RKHS-MP model gave the highest accuracies, while LS gave the lowest. GBLUP, GBLUP A, RKHS-M, and RKHS-P models gave similar accuracies. Using genome-wide marker-based models resulted in an average of 42% increase in accuracy over LS. We conclude that GS is a promising approach for improvement of quantitative rust resistance and can be implemented in the breeding pipeline.

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.

Sparse kernel models provide optimization of training set design for genomic prediction in multiyear wheat breeding data.

The success of genomic selection (GS) in breeding schemes relies on its ability to provide accurate predictions of unobserved lines at early stages. Multigeneration data provides opportunities to increase the training data size and thus, the likelihood of extracting useful information from ancestors to improve prediction accuracy. The genomic best linear unbiased predictions (GBLUPs) are performed by borrowing information through kinship relationships between individuals. Multigeneration data usually becomes heterogeneous with complex family relationship patterns that are increasingly entangled with each generation. Under these conditions, historical data may not be optimal for model training as the accuracy could be compromised. The sparse selection index (SSI) is a method for training set (TRN) optimization, in which training individuals provide predictions to some but not all predicted subjects. We added an additional trimming process to the original SSI (trimmed SSI) to remove less important training individuals for prediction. Using a large multigeneration (8 yr) wheat (Triticum aestivum L.) grain yield dataset (n = 68,836), we found increases in accuracy as more years are included in the TRN, with improvements of ∼0.05 in the GBLUP accuracy when using 5 yr of historical data relative to when using only 1 yr. The SSI method showed a small gain over the GBLUP accuracy but with an important reduction on the TRN size. These reduced TRNs were formed with a similar number of subjects from each training generation. Our results suggest that the SSI provides a more stable ranking of genotypes than the GBLUP as the TRN becomes larger.

Identification of Genomic Regions and Sources for Wheat Blast Resistance through GWAS in Indian Wheat Genotypes.

Wheat blast (WB) is a devastating fungal disease that has recently spread to Bangladesh and poses a threat to the wheat production in India, which is the second-largest wheat producing country in the world. In this study, 350 Indian wheat genotypes were evaluated for WB resistance in 12 field experiments in three different locations, namely Jashore in Bangladesh and Quirusillas and Okinawa in Bolivia. Single nucleotide polymorphisms (SNPs) across the genome were obtained using DArTseq® technology, and 7554 filtered SNP markers were selected for a genome-wide association study (GWAS). All the three GWAS approaches used identified the 2NS translocation as the only major source of resistance, explaining up to 32% of the phenotypic variation. Additional marker-trait associations were located on chromosomes 2B, 3B, 4D, 5A and 7A, and the combined effect of three SNPs (2B_180938790, 7A_752501634 and 5A_618682953) showed better resistance, indicating their additive effects on WB resistance. Among the 298 bread wheat genotypes, 89 (29.9%) carried the 2NS translocation, the majority of which (60 genotypes) were CIMMYT introductions, and 29 were from India. The 2NS carriers with a grand mean WB index of 6.6 showed higher blast resistance compared to the non-2NS genotypes with a mean index of 46.5. Of the 52 durum wheats, only one genotype, HI 8819, had the 2NS translocation and was the most resistant, with a grand mean WB index of 0.93. Our study suggests that the 2NS translocation is the only major resistance source in the Indian wheat panel analysed and emphasizes the urgent need to identify novel non-2NS resistance sources and genomic regions.

Dissecting the Genetic Architecture of Phenology Affecting Adaptation of Spring Bread Wheat Genotypes to the Major Wheat-Producing Zones in India.

Spring bread wheat adaptation to diverse environments is supported by various traits such as phenology and plant architecture. A large-scale genome-wide association study (GWAS) was designed to investigate and dissect the genetic architecture of phenology affecting adaptation. It used 48 datasets from 4,680 spring wheat lines. For 8 years (2014-2021), these lines were evaluated for days to heading (DH) and maturity (DM) at three sites: Jabalpur, Ludhiana, and Samastipur (Pusa), which represent the three major Indian wheat-producing zones: the Central Zone (CZ), North-Western Plain Zone (NWPZ), and North-Eastern Plain Zone (NEPZ), respectively. Ludhiana had the highest mean DH of 103.8 days and DM of 148.6 days, whereas Jabalpur had the lowest mean DH of 77.7 days and DM of 121.6 days. We identified 119 markers significantly associated with DH and DM on chromosomes 5B (76), 2B (18), 7D (10), 4D (8), 5A (1), 6B (4), 7B (1), and 3D (1). Our results clearly indicated the importance of the photoperiod-associated gene (Ppd-B1) for adaptation to the NWPZ and the Vrn-B1 gene for adaptation to the NEPZ and CZ. A maximum variation of 21.1 and 14% was explained by markers 2B_56134146 and 5B_574145576 linked to the Ppd-B1 and Vrn-B1 genes, respectively, indicating their significant role in regulating DH and DM. The results provide important insights into the genomic regions associated with the two phenological traits that influence adaptation to the major wheat-producing zones in India.

Genome-Wide Association Mapping Indicates Quantitative Genetic Control of Spot Blotch Resistance in Bread Wheat and the Favorable Effects of Some Spot Blotch Loci on Grain Yield.

Spot blotch caused by the fungus Bipolaris sorokiniana poses a serious threat to bread wheat production in warm and humid wheat-growing regions of the world. Hence, the major objective of this study was to identify consistent genotyping-by-sequencing (GBS) markers associated with spot blotch resistance using genome-wide association mapping on a large set of 6,736 advanced bread wheat breeding lines from the International Maize and Wheat Improvement Center. These lines were phenotyped as seven panels at Agua Fria, Mexico between the 2013-2014 and 2019-2020 crop cycles. We identified 214 significant spot blotch associated GBS markers in all the panels, among which only 96 were significant in more than one panel, indicating a strong environmental effect on the trait and highlights the need for multiple phenotypic evaluations to identify lines with stable spot blotch resistance. The 96 consistent GBS markers were on chromosomes 1A, 1B, 1D, 2A, 3B, 4A, 5B, 5D, 6B, 7A, 7B, and 7D, including markers possibly linked to the Lr46, Sb1, Sb2 and Sb3 genes. We also report the association of the 2NS translocation from Aegilops ventricosa with spot blotch resistance in some environments. Moreover, the spot blotch favorable alleles at the 2NS translocation and two markers on chromosome 3BS (3B_2280114 and 3B_5601689) were associated with increased grain yield evaluated at several environments in Mexico and India, implying that selection for favorable alleles at these loci could enable simultaneous improvement for high grain yield and spot blotch resistance. Furthermore, a significant relationship between the percentage of favorable alleles in the lines and their spot blotch response was observed, which taken together with the multiple minor effect loci identified to be associated with spot blotch in this study, indicate quantitative genetic control of resistance. Overall, the results presented here have extended our knowledge on the genetic basis of spot blotch resistance in bread wheat and further efforts to improve genetic resistance to the disease are needed for reducing current and future losses under climate change.

Genomic Selection: A Tool for Accelerating the Efficiency of Molecular Breeding for Development of Climate-Resilient Crops.

Since the inception of the theory and conceptual framework of genomic selection (GS), extensive research has been done on evaluating its efficiency for utilization in crop improvement. Though, the marker-assisted selection has proven its potential for improvement of qualitative traits controlled by one to few genes with large effects. Its role in improving quantitative traits controlled by several genes with small effects is limited. In this regard, GS that utilizes genomic-estimated breeding values of individuals obtained from genome-wide markers to choose candidates for the next breeding cycle is a powerful approach to improve quantitative traits. In the last two decades, GS has been widely adopted in animal breeding programs globally because of its potential to improve selection accuracy, minimize phenotyping, reduce cycle time, and increase genetic gains. In addition, given the promising initial evaluation outcomes of GS for the improvement of yield, biotic and abiotic stress tolerance, and quality in cereal crops like wheat, maize, and rice, prospects of integrating it in breeding crops are also being explored. Improved statistical models that leverage the genomic information to increase the prediction accuracies are critical for the effectiveness of GS-enabled breeding programs. Study on genetic architecture under drought and heat stress helps in developing production markers that can significantly accelerate the development of stress-resilient crop varieties through GS. This review focuses on the transition from traditional selection methods to GS, underlying statistical methods and tools used for this purpose, current status of GS studies in crop plants, and perspectives for its successful implementation in the development of climate-resilient crops.

Genome-Wide Association Study of Phytic Acid in Wheat Grain Unravels Markers for Improving Biofortification.

Biofortification of cereal grains offers a lasting solution to combat micronutrient deficiency in developing countries where it poses developmental risks to children. Breeding efforts thus far have been directed toward increasing the grain concentrations of iron (Fe) and zinc (Zn) ions. Phytic acid (PA) chelates these metal ions, reducing their bioavailability in the digestive tract. We present a high-throughput assay for quantification of PA and its application in screening a breeding population. After extraction in 96-well megatiter plates, PA content was determined from the phosphate released after treatment with a commercially available phytase enzyme. In a set of 330 breeding lines of wheat grown in the field over 3 years as part of a HarvestPlus breeding program for high grain Fe and Zn, our assay unraveled variation for PA that ranged from 0.90 to 1.72% with a mean of 1.24%. PA content was not associated with grain yield. High yielding lines were further screened for low molar PA/Fe and PA/Zn ratios for increased metal ion bioavailability, demonstrating the utility of our assay. Genome-wide association study revealed 21 genetic associations, six of which were consistent across years. Five of these associations mapped to chromosomes 1A, 2A, 2D, 5A, and 7D. Additivity over four of these haplotypes accounted for an ∼10% reduction in PA. Our study demonstrates it is possible to scale up assays to directly select for low grain PA in forward breeding programs.

Bayesian multitrait kernel methods improve multienvironment genome-based prediction.

When multitrait data are available, the preferred models are those that are able to account for correlations between phenotypic traits because when the degree of correlation is moderate or large, this increases the genomic prediction accuracy. For this reason, in this article, we explore Bayesian multitrait kernel methods for genomic prediction and we illustrate the power of these models with three-real datasets. The kernels under study were the linear, Gaussian, polynomial, and sigmoid kernels; they were compared with the conventional Ridge regression and GBLUP multitrait models. The results show that, in general, the Gaussian kernel method outperformed conventional Bayesian Ridge and GBLUP multitrait linear models by 2.2-17.45% (datasets 1-3) in terms of prediction performance based on the mean square error of prediction. This improvement in terms of prediction performance of the Bayesian multitrait kernel method can be attributed to the fact that the proposed model is able to capture nonlinear patterns more efficiently than linear multitrait models. However, not all kernels perform well in the datasets used for evaluation, which is why more than one kernel should be evaluated to be able to choose the best kernel.