Genome-wide association mapping and genomic prediction for pre‑harvest sprouting resistance, low α-amylase and seed color in Iranian bread wheat

Identification of Novel Loci Precisely Modulating Pre-Harvest Sprouting Resistance and Red Color Components of the Seed Coat in T. aestivum L

The association between pre-harvest sprouting (PHS) and seed coat color has long been recognized. Red-grained wheats generally exhibit greater PHS resistance compared to white-grained wheat, although variability in PHS resistance exists within red-grained varieties. Here, we conducted a genome-wide association study on a panel consisting of red-grained wheat varieties, aimed at uncovering genes that modulate PHS resistance and red color components of seed coat using digital image processing. Twelve loci associated with PHS traits were identified, nine of which were described for the first time. Genetic loci marked by SNPs AX-95172164 (chromosome 1B) and AX-158544327 (chromosome 7D) explained approximately 25% of germination index variance, highlighting their value for breeding PHS-resistant varieties. The most promising candidate gene for PHS resistance was TraesCS6B02G147900, encoding a protein involved in aleurone layer morphogenesis. Twenty-six SNPs were significantly associated with grain color, independently of the known Tamyb10 gene. Most of them were related to multiple color characteristics. Prioritization of genes within the revealed loci identified TraesCS1D03G0758600 and TraesCS7B03G1296800, involved in the regulation of pigment biosynthesis and in controlling pigment accumulation. In conclusion, our study identifies new loci associated with grain color and germination index, providing insights into the genetic mechanisms underlying these traits.

Genetic analyses and prediction for lodging‑related traits in a diverse Iranian hexaploid wheat collection

Lodging is one of the most important limiting environmental factors for achieving the maximum yield and quality of grains in cereals, including wheat. However, little is known about the genetic foundation underlying lodging resistance (LR) in wheat. In this study, 208 landraces and 90 cultivars were phenotyped in two cropping seasons (2018-2019 and 2019-2020) for 19 LR-related traits. A genome-wide association study (GWAS) and genomics prediction were carried out to dissect the genomic regions of LR. The number of significant marker pairs (MPs) was highest for genome B in both landraces (427,017) and cultivars (37,359). The strongest linkage disequilibrium (LD) between marker pairs was found on chromosome 4A (0.318). For stem lodging-related traits, 465, 497, and 478 marker-trait associations (MTAs) and 45 candidate genes were identified in year 1, year 2, and pooled. Gene ontology exhibited genomic region on Chr. 2B, 6B, and 7B control lodging. Most of these genes have key roles in defense response, calcium ion transmembrane transport, carbohydrate metabolic process, nitrogen compound metabolic process, and some genes harbor unknown functions that, all together may respond to lodging as a complex network. The module associated with starch and sucrose biosynthesis was highlighted. Regarding genomic prediction, the GBLUP model performed better than BRR and RRBLUP. This suggests that GBLUP would be a good tool for wheat genome selection. As a result of these findings, it has been possible to identify pivotal QTLs and genes that could be used to improve stem lodging resistance in Triticum aestivum L.

Identification and estimation of lodging in bread wheat genotypes using machine learning predictive algorithms

BACKGROUND

Lodging or stem bending decreases wheat yield quality and quantity. Thus, the traits reflected in early lodging wheat are helpful for early monitoring to some extent. In order to identify the superior genotypes and compare multiple linear regression (MLR) with support vector regression (SVR), artificial neural network (ANN), and random forest regression (RF) for predicting lodging in Iranian wheat accessions, a total of 228 wheat accessions were cultivated under field conditions in an alpha-lattice experiment, randomized incomplete block design, with two replications in two cropping seasons (2018-2019 and 2019-2020). To measure traits, a total of 20 plants were isolated from each plot and were measured using image processing.

RESULTS

The lodging score index (LS) had the highest positive correlation with plant height (r = 0.78**), Number of nodes (r = 0.71**), and internode length 1 (r = 0.70**). Genotypes were classified into four groups based on heat map output. The most lodging-resistant genotypes showed a lodging index of zero or close to zero. The findings revealed that the RF algorithm provided a more accurate estimate (R2 = 0.887 and RMSE = 0.091 for training data and R2 = 0.768 and RMSE = 0.124 for testing data) of wheat lodging than the ANN and SVR algorithms, and its robustness was as good as ANN but better than SVR.

CONCLUSION

Overall, it seems that the RF model can provide a helpful predictive and exploratory tool to estimate wheat lodging in the field. This work can contribute to the adoption of managerial approaches for precise and non-destructive monitoring of lodging.

Analysis of genetic diversity and genome-wide association study for drought tolerance related traits in Iranian bread wheat

BACKGROUND

Drought is most likely the most significant abiotic stress affecting wheat yield. The discovery of drought-tolerant genotypes is a promising strategy for dealing with the world's rapidly diminishing water resources and growing population. A genome-wide association study (GWAS) was conducted on 298 Iranian bread wheat landraces and cultivars to investigate the genetic basis of yield, yield components, and drought tolerance indices in two cropping seasons (2018-2019 and 2019-2020) under rainfed and well-watered environments.

RESULTS

A heatmap display of hierarchical clustering divided cultivars and landraces into four categories, with high-yielding and drought-tolerant genotypes clustering in the same group. The results of the principal component analysis (PCA) demonstrated that selecting genotypes based on the mean productivity (MP), geometric mean productivity (GMP), harmonic mean (HM), and stress tolerance index (STI) can help achieve high-yield genotypes in the environment. Genome B had the highest number of significant marker pairs in linkage disequilibrium (LD) for both landraces (427,017) and cultivars (370,359). Similar to cultivars, marker pairs on chromosome 4A represented the strongest LD (r2 = 0.32). However, the genomes D, A, and B have the highest LD, respectively. The single-locus mixed linear model (MLM) and multi-locus random-SNP-effect mixed linear model (mrMLM) identified 1711 and 1254 significant marker-trait association (MTAs) (-log10 P > 3) for all traits, respectively. A total of 874 common quantitative trait nucleotides (QTNs) were simultaneously discovered by both MLM and mrMLM methods. Gene ontology revealed that 11, 18, 6, and 11 MTAs were found in protein-coding regions (PCRs) for spike weight (SW), thousand kernel weight (TKW), grain number per spike (GN), and grain yield (GY), respectively.

CONCLUSION

The results identified rich regions of quantitative trait loci (QTL) on Ch. 4A and 5A suggest that these chromosomes are important for drought tolerance and could be used in wheat breeding programs. Furthermore, the findings indicated that landraces studied in Iranian bread wheat germplasm possess valuable alleles, that are responsive to water-limited conditions. This GWAS experiment is one of the few types of research conducted on drought tolerance that can be exploited in the genome-mediated development of novel varieties of wheat.

Uncovering the genetic basis for quality traits in the Mediterranean old wheat germplasm and phenotypic and genomic prediction assessment by cross-validation test

The release of new wheat varieties is based on two main characteristics, grain yield and quality, to meet the consumer's demand. Identifying the genetic architecture for yield and key quality traits has wide attention for genetic improvement to meet the global requirement. In this sense, the use of landraces represents an impressive source of natural allelic variation. In this study, a genome-wide association analysis (GWAS) with PCA and kinship matrix was performed to detect QTLs in bread wheat for fifteen quality and agronomic traits using 170 diverse landraces from 24 Mediterranean countries in two years of field trials. A total of 53 QTL hotspots containing 165 significant marker-trait associations (MTAs) were located across the genome for quality and agronomical traits except for chromosome 2D. The major specific QTL hotspots for quality traits were QTL_3B.3 (13 MTAs with a mean PVE of 8.2%) and QTL_4A.3 (15 MTAs, mean PVE of 11.0%), and for yield-related traits were QTL_2B.1 (8 MTAs, mean PVE of 7.4%) and QTL_4B.2 (5 MTAs, mean PVE of 10.0%). A search for candidate genes (CG) identified 807 gene models within the QTL hotspots. Ten of these CGs were expressed specifically in grain supporting the role of identified QTLs in Landraces, associated to bread wheat quality traits and grain formation. A cross-validation approach within the collection was performed to calculate the accuracies of genomic prediction for quality and agronomical traits, ranging from -0.03 to 0.64 for quality and 0.46 to 0.65 for agronomic traits. In addition, five prediction equations using the phenotypic data were developed to predict bread loaf volume in landraces. The prediction ability varied from 0.67 to 0.82 depending on the complexity of the traits considered to predict loaf volume.

Genome-wide association mapping and genomic prediction of agronomical traits and breeding values in Iranian wheat under rain-fed and well-watered conditions

BACKGROUND

The markers detected by genome-wide association study (GWAS) make it possible to dissect genetic structure and diversity at many loci. This can enable a wheat breeder to reveal and used genomic loci controlling drought tolerance. This study was focused on determining the population structure of Iranian 208 wheat landraces and 90 cultivars via genotyping-by-sequencing (GBS) and also on detecting marker-trait associations (MTAs) by GWAS and genomic prediction (GS) of wheat agronomic traits for drought-tolerance breeding. GWASs were conducted using both the original phenotypes (pGWAS) and estimated breeding values (eGWAS). The bayesian ridge regression (BRR), genomic best linear unbiased prediction (gBLUP), and ridge regression-best linear unbiased prediction (rrBLUP) approaches were used to estimate breeding values and estimate prediction accuracies in genomic selection.

RESULTS

Population structure analysis using 2,174,975 SNPs revealed four genetically distinct sub-populations from wheat accessions. D-Genome harbored the lowest number of significant marker pairs and the highest linkage disequilibrium (LD), reflecting different evolutionary histories of wheat genomes. From pGWAS, BRR, gBLUP, and rrBLUP, 284, 363, 359 and 295 significant MTAs were found under normal and 195, 365, 362 and 302 under stress conditions, respectively. The gBLUP with the most similarity (80.98 and 71.28% in well-watered and rain-fed environments, correspondingly) with the pGWAS method in the terms of discovered significant SNPs, suggesting the potential of gBLUP in uncovering SNPs. Results from gene ontology revealed that 29 and 30 SNPs in the imputed dataset were located in protein-coding regions for well-watered and rain-fed conditions, respectively. gBLUP model revealed genetic effects better than other models, suggesting a suitable tool for genome selection in wheat.

CONCLUSION

We illustrate that Iranian landraces of bread wheat contain novel alleles that are adaptive to drought stress environments. gBLUP model can be helpful for fine mapping and cloning of the relevant QTLs and genes, and for carrying out trait introgression and marker-assisted selection in both normal and drought environments in wheat collections.

Genome-Wide Association Study (GWAS) and genome prediction of seedling salt tolerance in bread wheat (Triticum aestivum L.)

BACKGROUND

Salinity tolerance in wheat is imperative for improving crop genetic capacity in response to the expanding phenomenon of soil salinization. However, little is known about the genetic foundation underlying salinity tolerance at the seedling growth stage of wheat. Herein, a GWAS analysis was carried out by the random-SNP-effect mixed linear model (mrMLM) multi-locus model to uncover candidate genes responsible for salt tolerance at the seedling stage in 298 Iranian bread wheat accessions, including 208 landraces and 90 cultivars.

RESULTS

A total of 29 functional marker-trait associations (MTAs) were detected under salinity, 100 mM NaCl (sodium chloride). Of these, seven single nucleotide polymorphisms (SNPs) including rs54146, rs257, rs37983, rs18682, rs55629, rs15183, and rs63185 with R² ≥ 10% were found to be linked with relative water content, root fresh weight, root dry weight, root volume, shoot high, proline, and shoot potassium (K⁺), respectively. Further, a total of 27 candidate genes were functionally annotated to be involved in response to the saline environment. Most of these genes have key roles in photosynthesis, response to abscisic acid, cell redox homeostasis, sucrose and carbohydrate metabolism, ubiquitination, transmembrane transport, chromatin silencing, and some genes harbored unknown functions that all together may respond to salinity as a complex network. For genomic prediction (GP), the genomic best linear unbiased prediction (GBLUP) model reflected genetic effects better than both bayesian ridge regression (BRR) and ridge regression-best linear unbiased prediction (RRBLUP), suggesting GBLUP as a favorable tool for wheat genomic selection.

CONCLUSION

The SNPs and candidate genes identified in the current work can be used potentially for developing salt-tolerant varieties at the seedling growth stage by marker-assisted selection.

Genome-wide association mapping for wheat morphometric seed traits in Iranian landraces and cultivars under rain-fed and well-watered conditions

Seed traits in bread wheat are valuable to breeders and farmers, thus it is important exploring putative QTLs responsible for key traits to be used in breeding programs. GWAS was carried out using 298 bread wheat landraces and cultivars from Iran to uncover the genetic basis of seed characteristics in both rain-fed and well-watered environments. The analyses of linkage disequilibrium (LD) between marker pairs showed that the largest number of significant LDs in landraces (427,017) and cultivars (370,359) was recorded in genome B, and the strongest LD was identified on chromosome 4A (0.318). LD decay was higher in the B and A genomes, compared to the D genome. Mapping by using mrMLM (LOD > 3) and MLM (0.05/m, Bonferroni) led to 246 and 67 marker-trait associations (MTAs) under rain-fed, as well as 257 and 74 MTAs under well-watered conditions, respectively. The study found that 3VmrMLM correctly detected all types of loci and estimated their effects in an unbiased manner, with high power and accuracy and a low false positive rate, which led to the identification of 140 MTAs (LOD > 3) in all environments. Gene ontology revealed that 10 and 10 MTAs were found in protein-coding regions for rain-fed and well-watered conditions, respectively. The findings suggest that landraces studied in Iranian bread wheat germplasm possess valuable alleles, which are responsive to water-limited conditions. MTAs uncovered in this study can be exploited in the genome-mediated development of novel wheat cultivars.