- Open Access
Genome-wide association mapping and genomic prediction for pre‑harvest sprouting resistance, low α-amylase and seed color in Iranian bread wheat
BMC Plant Biology volume 22, Article number: 300 (2022)
Pre-harvest sprouting (PHS) refers to a phenomenon, in which the physiologically mature seeds are germinated on the spike before or during the harvesting practice owing to high humidity or prolonged period of rainfall. Pre-harvest sprouting (PHS) remarkably decreases seed quality and yield in wheat; hence it is imperative to uncover genomic regions responsible for PHS tolerance to be used in wheat breeding. A genome-wide association study (GWAS) was carried out using 298 bread wheat landraces and varieties from Iran to dissect the genomic regions of PHS tolerance in a well-irrigated environment. Three different approaches (RRBLUP, GBLUP and BRR) were followed to estimate prediction accuracies in wheat genomic selection.
Genomes B, A, and D harbored the largest number of significant marker pairs (MPs) in both landraces (427,017, 328,006, 92,702 MPs) and varieties (370,359, 266,708, 63,924 MPs), respectively. However, the LD levels were found the opposite, i.e., genomes D, A, and B have the highest LD, respectively. Association mapping by using GLM and MLM models resulted in 572 and 598 marker-trait associations (MTAs) for imputed SNPs (− log10 P > 3), respectively. Gene ontology exhibited that the pleitropic MPs located on 1A control seed color, α-Amy activity, and PHS. RRBLUP model indicated genetic effects better than GBLUP and BRR, offering a favorable tool for wheat genomic selection.
Gene ontology exhibited that the pleitropic MPs located on 1A can control seed color, α-Amy activity, and PHS. The verified markers in the current work can provide an opportunity to clone the underlying QTLs/genes, fine mapping, and genome-assisted selection.Our observations uncovered key MTAs related to seed color, α-Amy activity, and PHS that can be exploited in the genome-mediated development of novel varieties in wheat.
Wheat (Triticum aestivum L.) has gradually become the global pioneer in supplying human nutrition and calories [1, 2]. The seeds of this crop are prone to sprouting at maturity when reiterated rainfall happens in the time of harvest in the field, leading to a remarkable decrease in flour quality and grain yield . As a result, pre-harvest sprouting (PHS) is known as a detrimental restricting factor in wheat productivity . Given this challenge, genetic improvements in PHS tolerance have become a serious focus of wheat breeders.
PHS tolerance depends on several factors, including i) environmental factors, such as relative humidity and temperature ; ii) biophysiological traits, such as germination-inhibitory compounds in the glumes, α-amylase (α-Amy) activity, grain structure and color, phytohormones, and seed dormancy ; iii) morphological traits, such as awn and spike structure . Of these factors, grain color is genetically related to PHS tolerance, the red-grained genotypes are more tolerant to PHS than white ones . Genes coding MYB transcriptional factors responsible for the flavonoid biosynthesis, i.e., Tamyb10-1, have been reported as candidates that determine grain color . Myb10 confers PHS resistance in wheat, which activates 9-cis-epoxycarotenoid dioxygenase (NCED) by biding the secondary wall MYB-responsive element (SMRE) to promote ABA biosynthesis in early wheat seed development stages [9,10,11]. Moreover, experimental evidence highlight seed dormancy is a key genetic component that determines PHS tolerance in wheat genotypes .
To date, numerous quantitative trait loci (QTLs) associated with PHS tolerance in wheat have been recognized in previous studies . These works have either assayed PHS tolerance indirectly by germination testing of harvest-ripe grains in a controlled environment [13,14,15] and/or directly by evaluating spikes in the field or in misting chambers [15, 16]. Most identified genomic segments are mapped on chromosome 4A [17,18,19], followed by 3A, 3B, and 3D . The PHS tolerance genes located on the chromosomes 3D, 3B, and 3A are known to be pleiotropic or closely linked with red coat controlled by allele R . Several resistant genes such as MKK3 , Vp1[22, 23], PM19 [24, 25], MFT , PHS1 , PHS-3D , ABI5 , FUS3  and DOG1 , were characterized from wheat for grain dormancy. Recently, Torada et al.  cloned MKK3 as the causal gene for grain dormancy. Further development of functional markers related to PHS tolerance is critical in wheat.
Genome-wide association study (GWAS) is an alternative tool to determine QTLs in natural populations . The establishment of genotyping technologies, from SSRs to SNPs, could facilitate association studies for accurate and efficient exploring of potential loci involved in complex traits, including PHS resistance in wheat [7, 13, 31] and grain-associated traits [32, 33]. However, the molecular mechanisms of PHS resistance remain unclear. Genomic selection (GS) along with GWAS can dramatically accelerate genetic gain in breeding [34, 35]. Several methods, including SNP-BLUP, have been suggested for genomic prediction .
In this study, a total of 298 Iranian wheat genotypes were evaluated for genotyping-by-sequencing (GBS)-based GWAS to achieve two objectives: i) uncovering genetic loci associated with PHS resistance; (2) identifying the best model for estimating prediction accuracies in genomic selection.
Phenotypic data summary
The results of descriptive statistics of traits related to pre‑harvest sprouting are shown in Table 1. Germination percentage occurred among Iranian wheat cultivars and landraces were ranged from zero to %100. The averages of germination percentage in landraces and cultivars were 71.31% and are 79.67%, respectively, which shows that native populations harbor more value of this trait. Sprouting index, sprouting score, and sprouting spike also confirm the lower pre‑harvest sprouting rate of native populations than cultivated varieties. The α-Amy enzyme activities in native populations and cultivars were 9.38 and 10.76, respectively, which indicates less activity of the enzyme in landraces than that of varieties. Color indices including L, a, and b do not differ significantly between cultivars and landraces.
From Fig. 1, wheat germination percentage (GP) indicated significant, negative correlations with most seed traits with coefficients ranging from 0.99 to 0.31 (P < 0.01). GP had the highest correlation with sprouting score (SS) (r = 0.99), followed by sprouting index (SI) (r = 0.98), α-Amy (A.amy) (r = 0.89), percentage of sprouted spike (SSp) (r = 0.68), color index L (r = 0.38), brightness index (WI) (r = 0.38), color index b, and Chroma (r = 0.25 and 0.24).
Assessment of SNPs
After eight Ion Proton runs, a total of 566,439,207 reads were identified with 458,363,607 (about 81%) high-quality barcoded reads. A total of 133,039 unique SNPs were called after filtering out duplicated reads. After imputation and discarding the SNPs with > 20% missing values, > 10% heterozygosity, and < 5% miner allele frequency, 43,525 SNPs were identified across all 21 wheat chromosomes. Out of them, 15,951, 21,864, and 5,710 SNPs were mapped to A, B, and D genomes, respectively, which included 36.7%, 50.2%, and 13.1% of total SNPs (Fig. 2). The highest and lowest numbers of SNPs were located on 3A (4034 SNPs) and 4D (270 SNPs), respectively.
Population structure and kinship matrix
In order to determine the appropriate number of subpopulations, the number of clusters was plotted (K) against ΔK. The largest ΔK value was observed at K = 3 suggesting the presence of three subpopulations (Fig. 3a). Using the structure software, the population of 298 accessions was structured into three subpopulations, Sub1, Sub2, and Sub_3 (Fig. 2). Sub_1 contains 113 accessions with 107 landraces and 6 varieties, Sub_2 contains 111 accessions with 97 landraces and 14 varieties; Sub_3 contains 74 studies with 70 varieties and 4 landraces (Fig. 3b). Molecular markers-based PCA showed that the first and second components justified 16.9% and 6.3% of total genetic variance occurred between wheat accessions. Thus, our study can distinguish favorably cultivars and native populations (Fig. 4). As expected, a population structure was identified in the Iranian wheat landraces, with the first five eigenvalues accounting for 30.5% of genetic diversity. From the clustering results, the native populations were divided into two subgroups. Clustering based on the nearest neighbor also indicated that cultivars and landraces were appropriately separated by using the imputed markers (Fig. 5).
Linkage disequilibrium (LD)
The levels of LD in genomes A, B, and D were 2279, 1707, and 5135, respectively. This reflects that genomes D, A, and B have the highest LD, respectively (Fig. 6). An analysis on landraces identified a total of 1,867,575 marker pairs with r2 = 0.182, of which 847,725 (45.39%) harbored significant linkages at P < 0.001.
Similar to cultivars, marker pairs on chromosome 4A showed the strongest LD (r2 = 0.369). Moreover, most of the significant marker pairs were found at distance of < 10 cM. Genomes D and B possessed the lowest and highest number of marker pairs (92,702 and 427,017), respectively. A total of 1,858,425 marker pairs with r2 = 0.211 were identified in cultivars, of which 700,991 (37.72%) harbored significant linkages at P < 0.001. Based on the observations, most of the significant marker pairs were found at distance of < 10 cM. Genomes D and B possessed the lowest and highest number of significant marker pairs (63,924 and 370,359), respectively (Table 2; Fig. 6).
MTAs for morphometric seed traits
In total, 566 and 598 significant marker pairs (MTAs) were identified by using GLM and MLM approaches, respectively, for PHS-related traits (–log10 P > 3). Of the total number of MTAs in the GLM method, 204, 271, and 97 MTAs were assigned to genomes A, B, and D, respectively. Of 598 MTAs in the MLM method, 220, 273, and 105 MTAs belonged to genomes A, B, and D, respectively. Genome B with 47.9% (GLM) and 45.7% (MLM) harbored the highest significant marker pairs and genome D with 16.1% (GLM) and 17.6% (MLM) possessed the lowest marker pairs, respectively. The number of significant markers for GP, SS, SI, SSp, A.amy, L, a, b, Hue, Chroma, and WI traits using the GLM method were 60, 65, 72, 120, 40, 30, 50, 35, 39, 35, and 20, as well as using the MLM method were 65, 66, 64, 170, 34, 30, 41, 35, 36, 37, and 20, respectively. The highest and lowest numbers of significant marker pairs using GLM and MLM methods were related to SSp (120 and 170 marker pairs) and WI (20 and 20 marker pairs), respectively. The most significant markers for PHS were on genome B, which has a greater effect on seed dormancy when compared to other genomes. However, the seed brightness (L and WI)-associated markers were located on genome A (Fig. 7). Manhattan diagrams for common areas associated with each seed trait are shown in Fig. 8.
The markers with the highest significance (P < 0.0001) and pleiotropy were studied in more detail. A total of 41 markers with high significance and pleitropic were identified, most of which were on 1A, 1B, 2A, 3B, 6D, and 7A. The marker pairs located on 1A were found to be able to control seed color, α-Amy activity, and germination percentage. Some of the significant MTAs were responsible for important molecular and biological processes, including protein kinase, G protein-coupled receptor signaling, signal transduction, intracellular transport, oxidoreductase activity, Fe ion binding, oxidation–reduction process, monooxygenase activity, protein binding, regulation of transcription, and double-stranded DNA binding (Table 3).
Based on the rice reference genome, the following pathways were discovered: hormone signal transduction (Fig. 9), metabolic pathways (Supplementary Fig. 1), MAPK signaling pathway (Supplementary Fig. 2), purine metabolism (Supplementary Fig. 3), spliceosome (Supplementary Fig. 4), and glycolysis/gluconeogenesis (Supplementary Fig. 5) ([37,38,39], www.kegg.jp/kegg/kegg1.html).
BRR, RR-BLUP, and GBLUP models using imputed SNPs exhibited the highest prediction accuracy for phenotypes 6, 3, and 2. The highest prediction accuracy by the GBLUP was achieved for SSp, Hue, and WI; by the RR-BLUP method for SS, SI, A.amy, a, L, and b; as well as by the BRR for GP and L traits (Fig. 10). BRR, RR-BLUP, and GBLUP models using significant SNPs indicated the highest prediction accuracy for phenotypes 2, 7, and 2. The highest prediction accuracy by the GBLUP was achieved for L and WI; by the RR-BLUP method for GP, SS, SI, SSp, Hue, a, and b; as well as by the BRR for A.amy trait. Overall, the RR-BLUP showed higher prediction accuracy and the BRR had a slight difference in accuracy with the RR-BLUP.
PHS tolerance in wheat is a complicated quantitative trait influenced by genetic background and environment . Thus, reliable phenotyping and genotyping for monitoring PHS tolerance can enhance the accuracy of QTL mapping. In this study, a total of 298 Iranian wheat accessions including 208 landraces and 90 cultivars were assembled as a natural population for mapping QTLs related to α-Amy enzyme activity, seed color, PHS using GWAS.
A high level of variation was uncovered in the studied traits for Iranian wheat accessions, suggesting the potential of the GWAS technique for exploring QTLs, as reported by Rahimi et al. . The α-Amy enzyme activity was lower in native populations than that of cultivars. Moreover, the seeds of landraces were exposed to longer dormancy when compared to cultivars. From correlation analysis, the below facts were concluded based on the associations among α-Amy activity, grain color, and pr-harvest sprouting, i) the lower the α-Amy activity, the higher the resistance of accessions to PHS; ii) the darker the seed, the more dormant it is; and iii) the more dormant the seed, the more PHS resistant it is. Similar results were reported by Zhou et al. , Zhou et al. , and Albrecht et al. .
The possibility for false associations can be raised in mapping studies if population structure is not suitably accounted for . Two kinds of kinships lead to a high rate of false positives in GWAS: cryptic relatedness and ancestry difference. Cryptic relatedness appears when some plant accessions are closely related; however, these shared ancestries are undisclosed to breeders . Large populations inevitably consist of accessions having common ancestry from various populations. Ancestry difference also refers to various ancestries among accessions in research . To evaluate the population structure in Iranian wheat accessions, PCA analysis and clustering were performed. Of results, the panel of accessions was stratified into three groups. The selection effects in breeding programs are considered as the reasons for such a genetic separation . Rahimi et al.  observed the same grouping on these Iranian wheat accessions. Cultivars made up one group, while landraces made up the other two groups, regardless of their geographic origins. This mixture can be originated from grain exchanges between farmers in different local markets throughout the country .
Of the results, the detected SNPs could cover the wheat genome well. The SNPs were higher in genome B and lower in genome D. Therefore, it seems there is a direct correlation between chromosome size and SNP density , because of the smaller size of chromosomes B compared to A ones. The higher frequency of SNPs in genome B resulted from the evolutionary processes. This inference was also stated by Alipour et al.  and Mourad et al. .
Genomes D, A, and B have the highest LD, respectively. The strongest LD was recorded between marker pairs on chromosome 4A. The fact that cultivars exhibited higher LD in contrast to landraces, particularly in genome D, is presumably a consequence of selection throughout the time of breeding efforts for PHS-related traits . The differences in LD occurred between genomes and accessions, in addition to the evolutionary processes, indicate the impact of breeding schedules. Similarly, Liu et al.  observed that the distance of LD decay in the native populations is less than cultivated varieties in wheat Pakistan/China collections.
Of the results, 1A, 2A, 4A, 1B, 2B, 6B, 4B, 3B, 5B, 7B, 6D, 5D, 4D, and 2D harbor genomic regions controlling PHS-related traits. Genome B possessed the highest number of MTAs, suggesting the potential of this genome in wheat adaptability to PHS. The most significant markers for PHS were on genome B, which has a greater effect on seed dormancy when compared to other genomes. However, the seed brightness-associated markers were located on genome A. These observations are in agreement with previous studies. For instance, Zhu et al.  mapped three key loci for PHS tolerance on chromosomes 6BL, 3BS, and 1AL, as well as validated one dCAPS and two CAPS markers for implementation in wheat genomics-based selection.
Genomic regions controlling PHS were detected in most wheat chromosomes in this study. To date, seven PHS QTLs were detected, including Qphs.ahau-7A.2, Qphs.ahau- 6A, Qphs.ahau-5D, Qphs.ahau-5B.4, Qphs.ahau-3B, Qphs. ahau-2B.3 and Qphs.ahau-2A.1 . Our observations showed that the darker the seed, the more dormant it will be and thus the more resistant it will be to PHS. Of justifying the cause, some associations were observed between grain color and PHS tolerance. Zhu et al. , for instance, discovered the positive correlations between PHS tolerance and seed color and suggested that this association occurs because the red-colored populations harbor more tolerant Qphs.ahau-1A and Qphs.ahau-3B alleles. Therefore, the authors stated that wheat seed color may be modulated collectively via Tamyb10-1 and other QTLs. In this work, MTAs related to grain color were found on 7B, 2A, etc. In this regard, the Psy1 gene coding phytoene synthase 1, responsible for yellow pigment, is co-segregated with seed brightness on 7B . A major QTL controlling both a* (redness) and L* (brightness) was also reported on 2A . Therefore, it seems that QTLs located on 7B and 2A are involved in wheat seed brightness, and thereby PHS tolerance. In this work, MTAs related to seed dormancy were found in some chromosomes, such as 4A. Similarly, Torada et al.  mapped TaMKK3-A as a candidate gene for the wheat seed dormancy, namely Phs1, on chromosomes 4A. They suggested that a single amino acid substitution in the kinase domain of this protein is related to the length of seed dormancy. From our findings, α-Amy-related genomic regions were found on 6B, 6D, 7B, etc. This is in line with previous studies. Lazarus et al.  demonstrated that α-Amy-related genomic regions are multigene families located on the chromosomes 7A, 7B, 7D (α-Amy2) and 6A, 6B, 6D (α-Amy1).
The flanking sequences of imputed SNPs were searched and aligned versus the RefSeq v1.0 (, https://urgi.versailles.inra.fr/blast_iwgsc/). Interestingly, output indicated that most marker pairs are in the protein-coding regions, which control the transcription process. DNA-binding, transcription factor activity, and transmembrane transport are other examples that are likely responsible for PHS tolerance. These findings are similar to the earlier researches . Based on the rice reference genome, the following pathways were discovered: metabolic pathways, hormone signal transduction, MAPK signaling pathway, purine metabolism, spliceosome, and glycolysis/gluconeogenesis. Liu et al.  observed that the slowed glycolysis leads to down-regulate glycerate-3-phosphate and inhibits seed germination (i.e., PHS). Torada et al.  uncovered a MKK3 by a map-based approach as a candidate gene for the locus Phs1 on 4A in wheat. Liu et al.  revealed that water status changes transcript levels of key genes involved in auxin, JA, and ethylene biosynthesis and their metabolic pathways, suggesting roles in regulating seed dormancy and germination. Nonogaki et al.  showed that seed germination and dormancy, the two main factors around PHS, are controled by endogenous hormone balance, especially between GA and ABA, reflecting their vital roles in PHS. Wang et al.  indicated that MAPK signaling and hormone signal transduction are associated with PHS. Zhang et al.  also highlighted that transcripts of spliceosome-related genes are abundant in the early stage of seed germination, suggesting the role of spliceosome in PHS process.
The highest prediction accuracy by GBLUP was achieved for SSp, Hue, and WI; by RR-BLUP method for SS, SI, A.amy, a, L, and b; as well as by BRR for GP and L traits. Shabannejad et al.  revealed BRR and RR-BLUP are superior to other GP models, which are utilized in well-irrigated and rain-fed environments, respectively. Overall, obtaining the highest GP accuracy is depend on the genomic selection method, level of LD, genetic architecture, and genetic variation . In this study, RRBLUP model indicated genetic effects better than GBLUP and BRR, offering a favorable tool for wheat genomic selection. It was reported that high genetic variation would be achieved by the GBLUP if markers were closely associated with the trait of interest and/or plant populations were advanced. RR-BLUP can work well for genetic architecture consisting of numerous loci with small impacts. BRR is similar to RR-BLUP however its shrinkage depends on the size of the studied population .
In the current study, GWAS for PHS in Iranian bread wheat accessions revealed the lowest LD decay distance and the highest number of marker pairs on genome B due to evolutionary events. The loci controlling the traits of interest were mapped on 1A, 2A, 4A, 1B, 2B, 7B, 3B, 5B, 6B, 4B, 6D, 2D, 5D, and 4D. Gene ontology exhibited that the pleitropic MPs located on 1A can control seed color, α-Amy activity, and PHS. The verified markers in the current work can provide an opportunity to clone the underlying QTLs/genes, fine mapping, and genome-assisted selection.
Material and methods
Plant material and field trial
To monitor PHS resistance, 208 native landraces and 90 cultivars were cultured in an alpha-lattice with two repeats in three crop seasons (2017–18, 2018–19, and 2019–20) under well-irrigated conditions (Table 4). The sizes of plots were adjusted as 1 m2. After physiological maturing, a total of 50 spikes were chosen from each plot and stored at -20 °C. After about a month, the spikes were taken out of the refrigerator and kept at 25 °C for 48 h. From each repeat, 10 healthy spikes were selected and soaked in distilled water for 3 h. Spikes immersed under 100% humidity were placed inside the controlled chambers, in which the steam and mist systems are utilized to spray and to maintain the moisture of the spikes, with a 16 h light/ 8 h dark photoperiod at 25 °C . The authors declare that all study complies with relevant institutional, national, and international guidelines and legislation for plant ethics in the methods section. Samples are provided from the Gene Bank of Agronomy and Plant Breeding Group and these samples are available at USDA and CIMMYT with USDA PI number and CIMMYT number (Supplementary Table 1), respectively. The authors declare that all that permissions or licenses were obtained to collect the wheat plant.
Pre harvest sprouting and α-Amy activity
After 7 days, PHS resistance was measured (Supplementary Table 2) based on the sprouting score as follows: the wheat spikes were given a score of one to nine, including one (germinated), two (less than 5%), three (5 to 15%), four (16 to 25%), five (26 to 45%), six (46 to 65%), seven (66 to 85%), eight (86 to 95%), and nine (more than 95%). The sprouting Index (SI, Eq. 1) was given a score of zero to five for each spike, in which zero was considered as the non-germinated spikes and five as 100% germinated spikes. The germination percentage (GP) and sprouted spikes (SS) were estimated from Eq. (2) and (3), respectively [59, 60].
Where n represents the number of clusters,
Where ni and N are the numbers of germinated and total seeds, respectively,
Where mi is the number of sprouted spikes and M is the total number of spikes.
To estimate α-Amy activity, the spikes of all accessions were taken out of the refrigerator, threshing was conducted by hand to avoid damaging the seed coat or embryos. Therefore, seeds were imbibed in a petri dish for a duration of 24 h at 25 °C and then prepared for enzyme extraction . 0.5 ml of the seed extract (60 mM phosphate buffer (pH 8.6) and 0.5 ml of starch solution were incubated at 37 °C for 30 min. The reaction was ceased by adding 1 ml of hydrochloric acid (0.1 N), and then 1 ml of the iodine reagent was added to the solution. The color absorption was recorded using a plate reader at 620 nm .
Evaluation of seed color with digital images
The digital images of wheat grains in the current work were provided by a camera (Canon SX540 HS) equipped with 800 dpi resolution. The captured images were analyzed and processed via Python 3.7 . For calibration, the regression between L, a, and b indices calculated with the Japanese CR_400 colorimeter and a photo box of 17 standard colors printed on 8 cm squares were used. Chroma saturation or index was calculated by Eq. (4), Hue Angle by Eq. (5), and Whiteness Index by Eq. (6).
Where, L, a, and b are color indices.
GBS and imputation
The GBS libraries were established and sequenced for the Iranian wheat genotypes following the procedure as explained by Alipour et al. . SNPs were discovered via internal alignments after trimming reads to 64 bp and categorizing them into tags. SNP calling was carried out using the UNEAK GBS pipeline, where SNPs with low allele frequency < 1% and reads with a low-quality score < 15 were removed to keep away from false-positive outputs. The imputation was accomplished according to available allele frequencies in BEAGLE version 3.3.2 . The distance of LD decay was determined by the ggplot2 package in RStudio . The W7984 reference genome was used because it fulfills the highest accuracy of imputation among various wheat reference genomes .
Population structure and kinship matrix
Population structure in the Iranian wheat accessions was assayed via STRUCTURE version 2.3.4. An admixture model was exploited along with a simulation phase consisting of 10,000 steps for K = 1–10. In this study, ΔK was exerted to estimate the most likely number of subpopulations . To measure LD among markers, the expected and observed allele frequencies were introduced into TASSEL. To determine the relationships among the Iranian wheat accessions, a neighbor-joining tree was constructed according to a pairwise distance matrix by TASSEL version 5 .
Genome-wide association study
The general linear model (GLM) and mixed linear model (MLM) approaches were accomplished to obtain the marker effect estimations. The GLM was performed with population structure (Q matrix) integrated as covariate to correct for the effects of population substructure. The MLM was employed with accounting for both population structure and family structure matrix (Kinship) to control both Type I and Type II errors. The association mapping was carried out using GLM and MLM functions of TASSEL5 [65, 68]. To correct for multiple testing, a false discovery rate (FDR) method described was used to declare significant marker-trait associations with relevant grain phenotype descriptor. A Manhattan plot was obtained using the CMplot package to explore associations between genotypes and phenotypes.
Annotation of genes
Sequences harboring associated SNP markers were exploited for the gene annotation by aligning to the IWGSC-RefSeq V1.0 (IWGSC) using Gramene (http://www.gramene.org/), an integrated database for comparative genomics in plant species. The overlapping genes with the highest blast score were picked out for further analysis. The ensemble-gramene database was utilized to extract the molecular functions and biological processes of genes in the gene ontology. Moreover, the significant SNPs were utilized in the enrichment analysis of gene ontology via KOBAS version 2.0 for testing in the KEGG (https://www.genome.jp/kegg/).
Genomic prediction strategies
GP was calculated by various approaches, including BRR , GBLUP , and RR-BLUP  based on whole 43,525 marker set and GWAS on the training set. All of the analyses were performed by iPat Tool . The GP accuracy was determined as Pearson’s correlations (r) between GEBVs and BLUPs over the validation and training sets .
The descriptive statistics and correlation analysis was conducted by R 4.1 using the ggplot2, dplyr, ggpubr and psych packages to reveal the distribution of wheat traits. To classify wheat genotypes, heatmap analysis was carried out in RStudio.
Availability of data and materials
The datasets generated and analyzed during the current study are available in the Figshare repository [https://doi.org/10.6084/m9.figshare.18774476.v1].
Rabieyan E, Alipour H. NGS-based multiplex assay of trait-linked molecular markers revealed the genetic diversity of Iranian bread wheat landraces and cultivars. Crop Pasture Sci. 2021;72(3):173–82. https://doi.org/10.1071/CP20362.
Rabieyan E, Bihamta MR, Moghaddam ME, Mohammadi V, Alipour H. Imaging-based screening of wheat seed characteristics towards distinguishing drought-responsive Iranian landraces and cultivars. Crop Pasture Sci. 2022;73(4):337–55. https://doi.org/10.1071/CP21500.
Zhu Y, Wang S, Wei W, Xie H, Liu K, Zhang C, Ma C. Genome-wide association study of pre-harvest sprouting tolerance using a 90K SNP array in common wheat (Triticum aestivum L.). Theor Appl Genet. 2019;132(11):2947–63. https://doi.org/10.1007/s00122-019-03398-x.
Zhou Y, Tang H, Cheng MP, Dankwa KO, Chen ZX, Li ZY, et al. Genome-wide association study for pre-harvest sprouting resistance in a large germplasm collection of Chinese wheat landraces. Front Plant Sci. 2017;8:401. https://doi.org/10.3389/fpls.2017.00401.
Mares DJ, Mrva K. Wheat grain pre-harvest sprouting and late maturity alpha-amylase. Planta. 2014;240(6):1167–78. https://doi.org/10.1007/s00425-014-2172-5.
Cao L, Hayashi K, Tokui M, Mori M, Miura H, Onishi K. Detection of QTLs for traits associated with pre-harvest sprouting resistance in bread wheat (Triticum aestivum L.). Breed Sci. 2016;66(2):260–70. https://doi.org/10.1270/jsbbs.66.260.
Lin M, Zhang DD, Liu SB, Zhang GR, Yu JM, Fritz AK, Bai GH. Genome-wide association analysis on pre-harvest sprouting resistance and grain color in U.S. winter wheat. BMC Genomics. 2016;17(1):794–810. https://doi.org/10.1186/s12864-016-3148-6.
Himi E, Maekawa M, Miura H, Noda K. Development of PCR markers for Tamyb10 related to R-1. red grain color gene in wheat. Theor Appl Genet. 2011;122(8):1561–76. https://doi.org/10.1007/s00122-011-1555-2.
Lang J, Fu Y, Zhou Y, Cheng M, Deng M, Li M, Wang J. Myb10-D confers PHS-3D resistance to pre-harvest sprouting by regulating NCED in ABA biosynthesis pathway of wheat. New Phytol. 2021;230(5):1940–52.
Yang J, Tan C, Lang J, Tang H, Hao M, Tan Z, Wang J. Identification of qPHS. sicau-1B and qPHS. sicau-3D from synthetic wheat for pre-harvest sprouting resistance wheat improvement. Mol Breeding. 2019;39(9):1–12.
Himi E, Noda K. Red grain colour gene (R) of wheat is a Myb-type transcription factor. Euphytica. 2005;143(3):239–42.
Kulwal P, Ishikawa G, Benscher D, Feng Z, Yu LX, Jadhav A, Sorrells ME. Association mapping for pre-harvest sprouting resistance in white winter wheat. Theor Appl Genet. 2012;125(4):793–805. https://doi.org/10.1007/s00122-012-1872-0.
Lin M, Cai SH, Wang S, Liu SB, Zhang GR, Bai GH. Genotyping- by-sequencing (GBS) identified SNP tightly linked to QTL for pre-harvest sprouting resistance. Theor Appl Genet. 2015;128(7):1385–95. https://doi.org/10.1007/s00122-015-2513-1.
Zhang YJ, Miao XL, Xia XC, He ZH. Cloning of seed dormancy genes (TaSdr) associated with tolerance to pre-harvest sprouting in common wheat and development of a functional marker. Theor Appl Genet. 2014;127(4):855–66. https://doi.org/10.1007/s00122-014-2262-6.
Somyong S, Ishikawa G, Munkvold JD, Tanaka J, Benscher D, Cho YG, Sorrells ME. Fine mapping of a preharvest sprouting QTL interval on chromosome 2B in white wheat. Theor Appl Genet. 2014;127(8):1843–55. https://doi.org/10.1007/s00122-014-2345-4.
Albrecht T, Oberforster M, Kempf H, Ramgraber L, Schacht J, Kazman E, et al. Genome wide association mapping of pre-harvest sprouting resistance in a diversity panel of European winter wheat. J Appl Genet. 2015;56(3):277–85. https://doi.org/10.1007/s13353-015-0286-5.
Cabral AL, Jordan MC, McCartney CA, You FM, Humphreys DG, MacLachlan R, Pozniak CJ. Identification of candidate genes. regions and markers for pre-harvest sprouting resistance in wheat (Triticum aestivum L.). BMC Plant Biol. 2014;14(1):1–12. https://doi.org/10.1186/s12870-014-0340-1.
Liu S, Bai G, Cai S, Chen C. Dissection of genetic components of preharvest sprouting resistance in white wheat. Mol Breed. 2011;27(4):511–23. https://doi.org/10.1007/s11032-010-9448-7.
Ogbonnaya FC, Imtiaz M, Ye G, Hearnden PR, Hernandez E, Eastwood RF. Genetic and QTL analyses of seed dormancy and preharvest sprouting resistance in the wheat germplasm CN10955. Theor Appl Genet. 2008;116(7):891–902. https://doi.org/10.1007/s00122-008-0712-8.
Liu S, Bai G. Dissection and fine mapping of a major QTL for preharvest sprouting resistance in white wheat Rio Blanco. Theor Appl Genet. 2010;121(8):1395–404. https://doi.org/10.1007/s00122-010-1396-4.
Torada A, Koike M, Ogawa T, Takenouchi Y, Tadamura K, Wu J, et al. A causal gene for seed dormancy on wheat chromosome 4A encodes a MAP kinase kinase. Curr Biol. 2016;26(6):782–7. https://doi.org/10.1016/j.cub.2016.01.063.
Chang C, Feng JM, Si HQ, Yin B, Zhang HP, Ma CX. Validating a novel allele of viviparous-1 (Vp-1Bf) associated with high seed dormancy of Chinese wheat landrace. Wanxianbaimaizi Mol Breeding. 2010;25(3):517–25.
Yang Y, Zhao XL, Xia LQ, Chen XM, Xia XC, Yu Z, Röder M. Development and validation of a Viviparous-1 STS marker for pre-harvest sprouting tolerance in Chinese wheats. Theor Appl Genet. 2007;115(7):971–80.
Barrero JM, Cavanagh C, Verbyla KL, Tibbits JF, Verbyla AP, Huang BE, Gubler F. Transcriptomic analysis of wheat near-isogenic lines identifies PM19-A1 and A2 as candidates for a major dormancy QTL. Genome Biol. 2015;16(1):1–18.
Torada A, Koike M, Ogawa T, Takenouchi Y, Tadamura K, Wu J, Ogihara Y. A causal gene for seed dormancy on wheat chromosome 4A encodes a MAP kinase kinase. Curr Biol. 2016;26(6):782–7.
Lei L, Zhu X, Wang S, Zhu M, Carver BF, Yan L. TaMFT-A1 is associated with seed germination sensitive to temperature in winter wheat. PLoS ONE. 2013;8(9):e73330.
Liu S, Sehgal SK, Lin M, Li J, Trick HN, Gill BS, Bai G. Independent mis-splicing mutations in Ta PHS 1 causing loss of preharvest sprouting (PHS) resistance during wheat domestication. New Phytol. 2015;208(3):928–35.
Yang Y, Ma YZ, Xu ZS, Chen XM, He ZH, Yu Z, Xia LQ. Isolation and characterization of Viviparous-1 genes in wheat cultivars with distinct ABA sensitivity and pre-harvest sprouting tolerance. J Exp Bot. 2007;58(11):2863–71.
Chiu RS, Nahal H, Provart NJ, Gazzarrini S. The role of the Arabidopsis FUSCA3transcription factor during inhibition of seed germination at high temperature. BMC Plant Biol. 2012;12(1):1–16.
Ashikawa I, Mori M, Nakamura S, Abe F. A transgenic approach to controlling wheat seed dormancy level by using Triticeae DOG1-like genes. Transgenic Res. 2014;23(4):621–9.
Rehman-Arif MA, Neumann K, Nagel M, Kobiljski B, Lohwasser U, Brner A. An association mapping analysis of dormancy and pre-harvest sprouting in wheat. Euphytica. 2012;88(3):409–17. https://doi.org/10.1007/s10681-012-0705-1.
Wang SX, Zhu YL, Zhang DX, Shao H, Liu P, Hu JB, et al. Genome wide association study for grain yield and related traits in elite wheat varieties and advanced lines using SNP markers. PLoS ONE. 2017;12(11):e0188662. https://doi.org/10.1371/journal.pone.0188662.
Chen GF, Zhang H, Deng ZY, Wu RG, Li DM, Wang MY, Tian JC. Genome-wide association study for kernel weight-related traits using SNPs in a Chinese winter wheat population. Euphytica. 2016;212(2):173–85. https://doi.org/10.1007/s10681-016-1750-y.
Lozada DN, Mason RE, Sarinelli JM, Brown-Guedira G. Accuracy of genomic selection for grain yield and agronomic traits in soft red winter wheat. BMC Genet. 2019;20(1):1–12. https://doi.org/10.1186/s12863-019-0785-1.
Juliana P, Poland J, Huerta-Espino J, Shrestha S, Crossa J, Crespo-Herrera L, Singh RP. Improving grain yield. stress resilience and quality of bread wheat using large-scale genomics. Nat Genet. 2019;51(10):1530–9. https://doi.org/10.1038/s41588-019-0496-6.
Tsai HY, Janss LL, Andersen JR, Orabi J, Jensen JD, Jahoor A, Jensen J. Genomic prediction and GWAS of yield quality and disease-related traits in spring barley and winter wheat. Sci Rep. 2020;10(1):1–15. https://doi.org/10.1038/s41598-020-60203-2.
Kanehisa M, Goto S. KEGG: Kyoto Encyclopedia of Genes and Genomes. Nucleic Acids Res. 2000;28:27–30. https://doi.org/10.1093/nar/28.1.27.
Kanehisa M. Toward understanding the origin and evolution of cellular organisms. Protein Sci. 2019;28:1947–51. https://doi.org/10.1002/pro.3715.
Kanehisa M, Furumichi M, Sato Y, Ishiguro-Watanabe M, Tanabe M. KEGG: integrating viruses and cellular organisms. Nucleic Acids Res. 2021;49:D545–51. https://doi.org/10.1093/nar/gkaa970.
Rahimi Y, Bihamta MR, Taleei A, Alipour H, Ingvarsson PK. Genome-wide association study of agronomic traits in bread wheat reveals novel putative alleles for future breeding programs. BMC Plant Biol. 2019;19(1):1–19. https://doi.org/10.1186/s12870-019-2165-4.
Wang X, Liu H, Siddique KH, Yan G. Transcriptomic profiling of wheat near-isogenic lines reveals candidate genes on chromosome 3A for pre-harvest sprouting resistance. BMC Plant Biol. 2021;21(1):1–14. https://doi.org/10.1186/s12870-021-02824-x.
Sul JH, Martin LS, Eskin E. Population structure in genetic studies: Confounding factors and mixed models. PLoS Genet. 2018;14(12):e1007309. https://doi.org/10.1371/journal.pgen.1007309.
Alipour H, Abdi H, Rahimi Y, Bihamta MR. Dissection of the genetic basis of genotype-by-environment interactions for grain yield and main agronomic traits in Iranian bread wheat landraces and cultivars. Sci Rep. 2021;11(1):1–17. https://doi.org/10.1038/s41598-021-96576-1.
Alemu A, Feyissa T, Tuberosa R, Maccaferri M, Sciara G, Letta T, Abeyo B. Genome-wide association mapping for grain shape and color traits in Ethiopian durum wheat (Triticum turgidum ssp. durum). The Crop Journal. 2020;8(5):757–68. https://doi.org/10.1016/j.cj.2020.01.001.
Sabzehzari M, Zeinali M, Naghavi MR. Alternative sources and metabolic engineering of Taxol: Advances and future perspectives. Biotechnol Adv. 2020;43: 107569. https://doi.org/10.1016/j.biotechadv.2020.107569.
Alipour H, Bihamta MR, Mohammadi V, Peyghambari SA, Bai G, Zhang G. Genotyping-by-sequencing (GBS) revealed molecular genetic diversity of Iranian wheat landraces and cultivars. Front Plant Sci. 2017;8:1293. https://doi.org/10.3389/fpls.2017.01293.
Mourad AMI, Belamkar V, Baenziger PS. Molecular genetic analysis of spring wheat core collection using genetic diversity population structure and linkage disequilibrium. BMC Genomics. 2020;21(1):1–12. https://doi.org/10.1186/s12864-020-06835-0.
Liu J, Rasheed A, He Z, Imtiaz M, Arif A, Mahmood T, Xia X. Genome-wide variation patterns between landraces and cultivars uncover divergent selection during modern wheat breeding. Theor Appl Genet. 2019;132(9):2509–23. https://doi.org/10.1007/s00122-019-03367-4.
Pozniak CJ, Knox RE, Clarke FR, Clarke JM. Identification of QTL and association of a phytoene synthase gene with endosperm colour in durum wheat. Theor Appl Genet. 2007;114(3):525–37. https://doi.org/10.1007/s00122-006-0453-5.
Lazarus CM, Baulcombe DC, Martienssen RA. α-amylase genes of wheat are two multigene families which are differentially expressed. Plant Mol Biol. 1985;5(1):13–24. https://doi.org/10.1007/BF00017869.
Appels R, Eversole K, Stein N, Feuillet C, Keller B, Rogers J, et al. Shifting the limits in wheat research and breeding using a fully annotated reference genome. Science. 2018;361:(6403). https://doi.org/10.1126/science.aar7191.
Liu C, Ding F, Hao F, Yu M, Lei H, Wu X, Tang H. Reprogramming of seed metabolism facilitates pre-harvest sprouting resistance of wheat. Sci Rep. 2016;6(1):1–11. https://doi.org/10.1038/srep20593.
Liu A, Gao F, Kanno Y, Jordan MC, Kamiya Y, Seo M, Ayele BT. Regulation of Wheat Seed Dormancy by After-Ripening Is Mediated by Specific Transcriptional Switches That Induce Changes in Seed Hormone Metabolism and Signaling. PLoS ONE. 2013;8(2):e56570. https://doi.org/10.1371/journal.pone.0056570.
Nonogaki H, Barrero JM, Li C. Editorial Seed dormancy germination and pre-harvest sprouting. Front Plant Sci. 2018;9:1783. https://doi.org/10.3389/fpls.2018.01783.
Zhang Q, Zhang X, Wang S, Tan C, Zhou G, Li C. Involvement of alternative splicing in barley seed germination. PLoS ONE. 2016;11(3):e0152824. https://doi.org/10.1371/journal.pone.0152824.
Shabannejad M, Bihamta MR, Majidi-Hervan E, Ebrahimi A. A simple. cost-effective high-throughput image analysis pipeline improves genomic prediction accuracy for days to maturity in wheat. Plant Methods. 2020;16(1):1–11. https://doi.org/10.1186/s13007-020-00686-2.
Sabzehzari M, Zeinali M, Naghavi MR. CRISPR-based metabolic editing: Next-generation metabolic engineering in plants. Gene. 2020;759:144993. https://doi.org/10.1016/j.gene.2020.144993.
Shabannejad M, Bihamta MR, Majidi-Hervan E, Alipour H, Ebrahimi A. A classic approach for determining genomic prediction accuracy under terminal drought stress and well-watered conditions in wheat landraces and cultivars. PLoS One. 2021;16(3):e0247824. https://doi.org/10.1371/journal.pone.0247824.
Clarke JM, DePauw RM, Mcleod JG, McCraig TN. Variation for preharvest sprouting resistance in durum wheat. Crop Sci. 1994;34(6):1632–5. https://doi.org/10.2135/cropsci1994.0011183X003400060038x.
Hucl P. Divergent selection for sprouting resistance in spring wheat. Plant Breeding. 1995;114(3):199–204. https://doi.org/10.1111/j.1439-0523.1995.tb00793.x.
Mori M, Uchino N, Chono M, Kato K, Miura H. Mapping QTLs for grain dormancy on wheat chromosome 3A and the group 4 chromosomes and their combined effect. Theor Appl Genet. 2005;110(7):1315–23. https://doi.org/10.1007/s00122-005-1972-1.
Zhang C, Si Y, Lamkey J, Boydston RA, Garland-Campbell KA, Sankaran S. High-throughput phenotyping of seed/seedling evaluation using digital image analysis. Agronomy. 2018;8(5):63. https://doi.org/10.3390/agronomy8050063.
Browning BL, Browning SR. A unified approach to genotype imputation and haplotype-phase inference for large data sets of trios and unrelated individuals. Am J Hum Genet. 2009;84(2):210–23. https://doi.org/10.1016/j.ajhg.2009.01.005.
Team R. RStudio: integrated development for R. RStudio. Inc. Boston. MA URL http://www.rstudio.com. 2015;42:14.
Alipour H, Bai G, Zhang G, Bihamta MR, Mohammadi V, Peyghambari SA. Imputation accuracy of wheat genotyping-by-sequencing (GBS) data using barley and wheat genome references. PLoS ONE. 2019;14(1):e0208614. https://doi.org/10.1371/journal.pone.0208614.
Pritchard JK, Stephens M, Donnelly P. Inference of population structure using multilocus genotype data. Genetics. 2000;155(2):945–59. https://doi.org/10.1093/genetics/155.2.945.
Bradbury PJ, Zhang Z, Kroon DE, Casstevens TM, Ramdoss Y, Buckler ES. TASSEL: software for association mapping of complex traits in diverse samples. Bioinformatics. 2007;23(19):2633–5. https://doi.org/10.1093/bioinformatics/btm308.
Lipka AE, Tian F, Wang Q, Peiffer J, Li M, Bradbury PJ, et al. GAPIT: genome association and prediction integrated tool. Bioinformatics. 2012;28(18):2397–9. https://doi.org/10.1093/bioinformatics/bts444.
Pérez P, de Los CG. Genome-wide regression and prediction with the BGLR statistical package. Genetics. 2014;198(2):483–95. https://doi.org/10.1534/genetics.114.164442.
Clark SA, van der Werf J. Genomic best linear unbiased prediction (gBLUP) for the estimation of genomic breeding values. Methods Mol Biol. 2013;1019:321–30. https://doi.org/10.1007/978-1-62703-447-0_13.
Endelman JB. Ridge regression and other kernels for genomic selection with R package rrBLUP. The Plant Genome. 2011;4(3):250–5. https://doi.org/10.3835/plantgenome2011.08.0024.
Chen CJ, Zhang Z. iPat: intelligent prediction and association tool for genomic research. Bioinformatics. 2018;34(11):1925–7. https://doi.org/10.1093/bioinformatics/bty015.
Resende MF, Munoz P, Resende MD, Garrick DJ, Fernando RL, Davis JM, Kirst M. Accuracy of genomic selection methods in a standard data set of loblolly pine (Pinus taeda L). Genetics. 2012;190(4):1503–10. https://doi.org/10.1534/genetics.111.137026.
This research did not receive any specific funding.
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Overview on the landraces and cultivars of Iranian wheat studied. Supplementary Table 2. Phenotypic (germination, grain color, alpha-amylase) data measured in landraces and cultivars of Iranian wheat studied. Supplementary Fig 1. The KEGG pathway of metabolic pathways. Supplementary Fig 2. The KEGG pathway of MAPK signaling. Supplementary Fig 3. The KEGG pathway of purine metabolism. Supplementary Fig 4. The KEGG pathway of spliceosome. Supplementary Fig 5. The KEGG pathway of glycolysis/gluconeogenesis.
About this article
Cite this article
Rabieyan, E., Bihamta, M.R., Moghaddam, M.E. et al. Genome-wide association mapping and genomic prediction for pre‑harvest sprouting resistance, low α-amylase and seed color in Iranian bread wheat. BMC Plant Biol 22, 300 (2022). https://doi.org/10.1186/s12870-022-03628-3
- Pre-harvest sprouting