Genome-wide association study and quantitative trait loci mapping of seed dormancy in common wheat (Triticum aestivum L.)

Linkage and association analysis to identify wheat pre-harvest sprouting resistance genetic regions and develop KASP markers

Pre-harvest sprouting (PHS) of wheat (Triticum aestivum L.) is one of the complex traits that result in rainfall-dependent reductions in grain production and quality worldwide. Breeding new varieties and germplasm with PHS resistance is of great importance to reduce this problem. However, research on markers and genes related to PHS resistance is limited, especially in marker-assisted selection (MAS) wheat breeding. To this end, we studied PHS resistance in recombinant inbred line (RIL) population and in 171 wheat germplasm accessions in different environments and genotyped using the wheat Infinium 50 K/660 K SNP array. Quantitative trait loci (QTL) mapping and genome-wide association studies (GWAS) identified 59 loci controlling PHS. Upon comparison with previously reported QTL affecting PHS, 16 were found to be new QTL, and the remaining 43 loci were co-localized with QTL from previous studies. We also pinpointed 12 candidate genes within these QTL intervals that share functional similarities with genes previously known to influence PHS resistance. In addition, we developed and validated two kompetitive allele-specific PCR (KASP) markers within the chromosome 7B region identified by linkage analysis. These QTL, candidate genes, and the KASP marker identified in this study have the potential to improve PHS resistance of wheat, and they may enhance our understanding of the genetic basis of PHS resistance, thus being useful for MAS breeding.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-024-01526-0.

Functional analysis of a wheat class III peroxidase gene, TaPer12-3A, in seed dormancy and germination

BACKGROUND

Class III peroxidases (PODs) perform crucial functions in various developmental processes and responses to biotic and abiotic stresses. However, their roles in wheat seed dormancy (SD) and germination remain elusive.

RESULTS

Here, we identified a wheat class III POD gene, named TaPer12-3A, based on transcriptome data and expression analysis. TaPer12-3A showed decreasing and increasing expression trends with SD acquisition and release, respectively. It was highly expressed in wheat seeds and localized in the endoplasmic reticulum and cytoplasm. Germination tests were performed using the transgenic Arabidopsis and rice lines as well as wheat mutant mutagenized with ethyl methane sulfonate (EMS) in Jing 411 (J411) background. These results indicated that TaPer12-3A negatively regulated SD and positively mediated germination. Further studies showed that TaPer12-3A maintained H2O2 homeostasis by scavenging excess H2O2 and participated in the biosynthesis and catabolism pathways of gibberellic acid and abscisic acid to regulate SD and germination.

CONCLUSION

These findings not only provide new insights for future functional analysis of TaPer12-3A in regulating wheat SD and germination but also provide a target gene for breeding wheat varieties with high pre-harvest sprouting resistance by gene editing technology.

Genome-Wide Association Study of Preharvest Sprouting in Wheat

Genome-wide association study (GWAS) is widely used to characterize genes or quantitative trait loci (QTLs) associated with preharvest sprouting and seed dormancy. GWAS can identify both previously discovered and novel QTLs across diverse genetic panels. The high-throughput SNP arrays or next-generation sequencing technologies have facilitated the identification of numerous genetic markers, thereby significantly enhancing the resolution of GWAS. Although various methods have been developed, the fundamental principles underlying these techniques remain constant. Here, we provide a basic technological flow to perform seed dormancy assay, followed by GWAS using population structure control, and compared it with previous identified QTLs and genes.

Genetic Improvement of Wheat with Pre-Harvest Sprouting Resistance in China

Wheat pre-harvest sprouting (PHS) refers to the germination of seeds directly on the spike due to rainy weather before harvest, which often results in yield reduction, quality deterioration, and seed value loss. In this study, we reviewed the research progress in the quantitative trait loci (QTL) detection and gene excavation related to PHS resistance in wheat. Simultaneously, the identification and creation of germplasm resources and the breeding of wheat with PHS resistance were expounded in this study. Furthermore, we also discussed the prospect of molecular breeding during genetic improvement of PHS-resistant wheat.

GWAS and genomic prediction for pre-harvest sprouting tolerance involving sprouting score and two other related traits in spring wheat

In wheat, a genome-wide association study (GWAS) and genomic prediction (GP) analysis were conducted for pre-harvest sprouting (PHS) tolerance and two of its related traits. For this purpose, an association panel of 190 accessions was phenotyped for PHS (using sprouting score), falling number, and grain color over two years and genotyped with 9904 DArTseq based SNP markers. GWAS for main-effect quantitative trait nucleotides (M-QTNs) using three different models (CMLM, SUPER, and FarmCPU) and epistatic QTNs (E-QTNs) using PLINK were performed. A total of 171 M-QTNs (CMLM, 47; SUPER, 70; FarmCPU, 54) for all three traits, and 15 E-QTNs involved in 20 first-order epistatic interactions were identified. Some of the above QTNs overlapped the previously reported QTLs, MTAs, and cloned genes, allowing delineating 26 PHS-responsive genomic regions that spread over 16 wheat chromosomes. As many as 20 definitive and stable QTNs were considered important for use in marker-assisted recurrent selection (MARS). The gene, TaPHS1, for PHS tolerance (PHST) associated with one of the QTNs was also validated using the KASP assay. Some of the M-QTNs were shown to have a key role in the abscisic acid pathway involved in PHST. Genomic prediction accuracies (based on the cross-validation approach) using three different models ranged from 0.41 to 0.55, which are comparable to the results of previous studies. In summary, the results of the present study improved our understanding of the genetic architecture of PHST and its related traits in wheat and provided novel genomic resources for wheat breeding based on MARS and GP.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-023-01357-5.

Exogenous gibberellic acid shortening after-ripening process and promoting seed germination in a medicinal plant Panax notoginseng

BACKGROUND

Panax notoginseng (Burk) F.H. Chen is an essential plant in the family of Araliaceae. Its seeds are classified as a type of morphophysiological dormancy (MPD), and are characterized by recalcitrance during the after-ripening process. However, it is not clear about the molecular mechanism on the after-ripening in recalcitrant seeds.

RESULTS

In this study, exogenous supply of gibberellic acid (GA₃) with different concentrations shortened after-ripening process and promoted the germination of P. notoginseng seeds. Among the identified plant hormone metabolites, exogenous GA₃ results in an increased level of endogenous hormone GA₃ through permeation. A total of 2971 and 9827 differentially expressed genes (DEGs) were identified in response to 50 mg L^-1 GA₃ (LG) and 500 mg L^-1 GA₃ (HG) treatment, respectively, and the plant hormone signal and related metabolic pathways regulated by GA₃ was significantly enriched. Weighted gene co-expression network analysis (WGCNA) revealed that GA₃ treatment enhances GA biosynthesis and accumulation, while inhibiting the gene expression related to ABA signal transduction. This effect was associated with higher expression of crucial seed embryo development and cell wall loosening genes, Leafy Contyledon1 (LEC1), Late Embryogenesis Abundant (LEA), expansins (EXP) and Pectinesterase (PME).

CONCLUSIONS

Exogenous GA₃ application promotes germination and shorts the after-ripening process of P. notoginseng seeds by increasing GA₃ contents through permeation. Furthermore, the altered ratio of GA and ABA contributes to the development of the embryo, breaks the mechanical constraints of the seed coat and promotes the protrusion of the radicle in recalcitrant P. notoginseng seeds. These findings improve our knowledge of the contribution of GA to regulating the dormancy of MPD seeds during the after-ripening process, and provide new theoretical guidance for the application of recalcitrant seeds in agricultural production and storage.

Identification and validation of coding and non-coding RNAs involved in high-temperature-mediated seed dormancy in common wheat

INTRODUCTION

Seed dormancy (SD) significantly decreases under high temperature (HT) environment during seed maturation, resulting in pre-harvest sprouting (PHS) damage under prolonged rainfall and wet weather during wheat harvest. However, the molecular mechanism underlying HT-mediated SD remains elusiveSeed dormancy (SD) significantly decreases under high temperature (HT) environment during seed maturation, resulting in pre-harvest sprouting (PHS) damage under prolonged rainfall and wet weather during wheat harvest. However, the molecular mechanism underlying HT-mediated SD remains elusive.

METHODS

Here, the wheat landrace 'Waitoubai' with strong SD and PHS resistance was treated with HT from 21 to 35 days post anthesis (DPA). Then, the seeds under HT and normal temperature (NT) environments were collected at 21 DPA, 28 DPA, and 35 DPA and subjected to whole-transcriptome sequencing.

RESULTS

The phenotypic data showed that the seed germination percentage significantly increased, whereas SD decreased after HT treatment compared with NT, consistent with the results of previous studies. In total, 5128 mRNAs, 136 microRNAs (miRNAs), 273 long non-coding RNAs (lncRNAs), and 21 circularRNAs were found to be responsive to HT, and some of them were further verified through qRT-PCR. In particular, the known gibberellin (GA) biosynthesis gene TaGA20ox1 (TraesCS3D02G393900) was proved to be involved in HT-mediated dormancy by using the EMS-mutagenized wheat cultivar Jimai 22. Similarly, a novel gene TaCDPK21 (TraesCS7A02G267000) involved in the calcium signaling pathway was validated to be associated with HT-mediated dormancy by using the EMS mutant. Moreover, TaCDPK21 overexpression in Arabidopsis and functional complementarity tests supported the negative role of TaCDPK21 in SD. We also constructed a co-expression regulatory network based on differentially expressed mRNAs, miRNAs, and lncRNAs and found that a novel miR27319 was located at a key node of this regulatory network. Subsequently, using Arabidopsis and rice lines overexpressing miR27319 precursor or lacking miR27319 expression, we validated the positive role of miR27319 in SD and further preliminarily dissected the molecular mechanism of miR27319 underlying SD regulation through phytohormone abscisic acid and GA biosynthesis, catabolism, and signaling pathways.

DISCUSSION

These findings not only broaden our understanding of the complex regulatory network of HT-mediated dormancy but also provide new gene resources for improving wheat PHS resistance to minimize PHS damage by using the molecular pyramiding approach.

Genome-wide association study reveals a NAC transcription factor TaNAC074 linked to pre-harvest sprouting tolerance in wheat

Twelve QTL associated with pre-harvest sprouting tolerance were identified using association analysis in wheat. Two markers were validated and a candidate gene TaNAC074 for Qgpf.cas-3B.2 was verified using Agrobacterium-mediated transformation. Pre-harvest sprouting (PHS) is a considerable global threat to wheat yield and quality. Due to this threat, breeders must identify quantitative trait loci (QTL) and genes conferring PHS-tolerance (PHST) to reduce the negative effects of PHS caused by low seed dormancy. In this study, we evaluated a panel of 302 diverse wheat genotypes for PHST in four environments and genotyped the panel with a high-density wheat 660 K SNP array. By using a genome-wide association study (GWAS), we identified 12 stable loci significantly associated with PHST (P < 0.0001), explaining 3.34 - 9.88% of the phenotypic variances. Seven of these loci co-located with QTL and genes reported previously. Five loci (Qgpf.cas-3B.2, Qgpf.cas-3B.3, Qgpf.cas-3B.4, Qgpf.cas-7B.2, and Qgpf.cas-7B.3), located in genomic regions with no known PHST QTL or genes, are likely to be new QTL conferring PHST. Additionally, two molecular markers were developed for Qgpf.cas-3A and Qgpf.cas-7B.3, and validated using a different set of 233 wheat accessions. Finally, the PHST-related function of candidate gene TaNAC074 for Qgpf.cas-3B.2 was confirmed by CAPS (cleaved amplified polymorphic sequences) marker association analysis in 233 wheat accessions and by expression and phenotypic analysis of transgenic wheat. Overexpression of TaNAC074 significantly reduced seed dormancy in wheat. This study contributes to broaden the genetic basis and molecular marker-assisted breeding of PHST.

Combining QTL mapping and gene co-expression network analysis for prediction of candidate genes and molecular network related to yield in wheat

BACKGROUND

Wheat (Triticum aestivum L.) is an important cereal crop. Increasing grain yield for wheat is always a priority. Due to the complex genome of hexaploid wheat with 21 chromosomes, it is difficult to identify underlying genes by traditional genetic approach. The combination of genetics and omics analysis has displayed the powerful capability to identify candidate genes for major quantitative trait loci (QTLs), but such studies have rarely been carried out in wheat. In this study, candidate genes related to yield were predicted by a combined use of linkage mapping and weighted gene co-expression network analysis (WGCNA) in a recombinant inbred line population.

RESULTS

QTL mapping was performed for plant height (PH), spike length (SL) and seed traits. A total of 68 QTLs were identified for them, among which, 12 QTLs were stably identified across different environments. Using RNA sequencing, we scanned the 99,168 genes expression patterns of the whole spike for the recombinant inbred line population. By the combined use of QTL mapping and WGCNA, 29, 47, 20, 26, 54, 46 and 22 candidate genes were predicted for PH, SL, kernel length (KL), kernel width, thousand kernel weight, seed dormancy, and seed vigor, respectively. Candidate genes for different traits had distinct preferences. The known PH regulation genes Rht-B and Rht-D, and the known seed dormancy regulation genes TaMFT can be selected as candidate gene. Moreover, further experiment revealed that there was a SL regulatory QTL located in an interval of about 7 Mbp on chromosome 7A, named TaSL1, which also involved in the regulation of KL.

CONCLUSIONS

A combination of QTL mapping and WGCNA was applied to predicted wheat candidate genes for PH, SL and seed traits. This strategy will facilitate the identification of candidate genes for related QTLs in wheat. In addition, the QTL TaSL1 that had multi-effect regulation of KL and SL was identified, which can be used for wheat improvement. These results provided valuable molecular marker and gene information for fine mapping and cloning of the yield-related trait loci in the future.

Genome-Wide Association Study of Kernel Traits Using a 35K SNP Array in Bread Wheat (Triticum aestivum L.)

Kernel size and weight are crucial components of grain yield in wheat. Deciphering their genetic basis is essential for improving yield potential in wheat breeding. In this study, five kernel traits, including kernel length (KL), kernel width (KW), kernel diameter ratio (KDR), kernel perimeter (KP), and thousand-kernel weight (TKW), were evaluated in a panel consisting of 198 wheat accessions under six environments. Wheat accessions were genotyped using the 35K SNP iSelect chip array, resulting in a set of 13,228 polymorphic SNP markers that were used for genome-wide association study (GWAS). A total of 146 significant marker-trait associations (MTAs) were identified for five kernel traits on 21 chromosomes [-log₁₀(P) ≥ 3], which explained 5.91-15.02% of the phenotypic variation. Of these, 12 stable MTAs were identified in multiple environments, and six superior alleles showed positive effects on KL, KP, and KDR. Four potential candidate genes underlying the associated SNP markers were predicted for encoding ML protein, F-box protein, ethylene-responsive transcription factor, and 1,4-α-glucan branching enzyme. These genes were strongly expressed in grain development at different growth stages. The results will provide new insights into the genetic basis of kernel traits in wheat. The associated SNP markers and predicted candidate genes will facilitate marker-assisted selection in wheat breeding.

Comprehensive evaluation of mapping complex traits in wheat using genome-wide association studies

Genome-wide association studies (GWAS) are effectively applied to detect the marker trait associations (MTAs) using whole genome-wide variants for complex quantitative traits in different crop species. GWAS has been applied in wheat for different quality, biotic and abiotic stresses, and agronomic and yield-related traits. Predictions for marker-trait associations are controlled with the development of better statistical models taking population structure and familial relatedness into account. In this review, we have provided a detailed overview of the importance of association mapping, population design, high-throughput genotyping and phenotyping platforms, advancements in statistical models and multiple threshold comparisons, and recent GWA studies conducted in wheat. The information about MTAs utilized for gene characterization and adopted in breeding programs is also provided. In the literature that we surveyed, as many as 86,122 wheat lines have been studied under various GWA studies reporting 46,940 loci. However, further utilization of these is largely limited. The future breakthroughs in area of genomic selection, multi-omics-based approaches, machine, and deep learning models in wheat breeding after exploring the complex genetic structure with the GWAS are also discussed. This is a most comprehensive study of a large number of reports on wheat GWAS and gives a comparison and timeline of technological developments in this area. This will be useful to new researchers or groups who wish to invest in GWAS.

Phenotypic Evaluation and Genetic Analysis of Seedling Emergence in a Global Collection of Wheat Genotypes (Triticum aestivum L.) Under Limited Water Availability

The challenge in establishing an early-sown wheat crop in southern Australia is the need for consistently high seedling emergence when sowing deep in subsoil moisture (>10 cm) or into dry top-soil (4 cm). However, the latter is strongly reliant on a minimum soil water availability to ensure successful seedling emergence. This study aimed to: (1) evaluate 233 Australian and selected international wheat genotypes for consistently high seedling emergence under limited soil water availability when sown in 4 cm of top-soil in field and glasshouse (GH) studies; (2) ascertain genetic loci associated with phenotypic variation using a genome-wide association study (GWAS); and (3) compare across loci for traits controlling coleoptile characteristics, germination, dormancy, and pre-harvest sprouting. Despite significant (P < 0.001) environment and genotype-by-environment interactions within and between field and GH experiments, eight genotypes that included five cultivars, two landraces, and one inbred line had consistently high seedling emergence (mean value > 85%) across nine environments. Moreover, 21 environment-specific quantitative trait loci (QTL) were detected in GWAS analysis on chromosomes 1B, 1D, 2B, 3A, 3B, 4A, 4B, 5B, 5D, and 7D, indicating complex genetic inheritance controlling seedling emergence. We aligned QTL for known traits and individual genes onto the reference genome of wheat and identified 16 QTL for seedling emergence in linkage disequilibrium with coleoptile length, width, and cross-sectional area, pre-harvest sprouting and dormancy, germination, seed longevity, and anthocyanin development. Therefore, it appears that seedling emergence is controlled by multifaceted networks of interrelated genes and traits regulated by different environmental cues.

Genome-Wide Linkage Mapping for Preharvest Sprouting Resistance in Wheat Using 15K Single-Nucleotide Polymorphism Arrays

Preharvest sprouting (PHS) significantly reduces grain yield and quality. Identification of genetic loci for PHS resistance will facilitate breeding sprouting-resistant wheat cultivars. In this study, we constructed a genetic map comprising 1,702 non-redundant markers in a recombinant inbred line (RIL) population derived from cross Yangxiaomai/Zhongyou9507 using the wheat 15K single-nucleotide polymorphism (SNP) assay. Four quantitative trait loci (QTL) for germination index (GI), a major indicator of PHS, were identified, explaining 4.6-18.5% of the phenotypic variances. Resistance alleles of Qphs.caas-3AL, Qphs.caas-3DL, and Qphs.caas-7BL were from Yangxiaomai, and Zhongyou9507 contributed a resistance allele in Qphs.caas-4AL. No epistatic effects were detected among the QTL, and combined resistance alleles significantly increased PHS resistance. Sequencing and linkage mapping showed that Qphs.caas-3AL and Qphs.caas-3DL corresponded to grain color genes Tamyb10-A and Tamyb10-D, respectively, whereas Qphs.caas-4AL and Qphs.caas-7BL were probably new QTL for PHS. We further developed cost-effective, high-throughput kompetitive allele-specific PCR (KASP) markers tightly linked to Qphs.caas-4AL and Qphs.caas-7BL and validated their association with GI in a test panel of cultivars. The resistance alleles at the Qphs.caas-4AL and Qphs.caas-7BL loci were present in 72.2 and 16.5% cultivars, respectively, suggesting that the former might be subjected to positive selection in wheat breeding. The findings provide not only genetic resources for PHS resistance but also breeding tools for marker-assisted selection.

Dormancy and dormancy release in white-grained wheat (Triticum aestivum L.)

Dormancy in white-grained wheat is conditioned by the cumulative effects of several QTL that delay the onset of the capacity to germinate during ripening and after-ripening. Grain dormancy at harvest-ripeness is a major component of resistance to preharvest sprouting in wheat (Triticum aestivum L.) and an important trait in regions where rain is common during the harvest period. Breeding lines developed in Australia maintained their dormancy phenotype over multiple seasons and during grain ripening, the time between anthesis and the acquisition of the capacity to germinate, dormancy release, increased in line with the strength of dormancy. Genetic dissection of two dormant lines indicated that dormancy was due to the cumulative action of between one and three major genetic loci and several minor loci. This presents a significant challenge for breeders targeting environments with a high risk of sprouting where strong dormancy is desirable. Only around half of the difference in dormancy between the dormant lines and a non-dormant variety could be attributed to the major genetic loci on chromosomes 4A and 3A. A QTL that was mapped on chromosome 5A may be an orthologue of a minor QTL for dormancy in barley. At each locus, the dormancy allele increased the time to dormancy release during ripening. In combination, these alleles had cumulative effects. Embryo sensitivity to abscisic acid was related to the dormancy phenotype of the whole caryopsis, however, changes in concentrations of abscisic acid and gibberellins in embryo sections and de-embryonated grains during ripening and after-ripening could not be linked to the timing of dormancy release.

Exome sequencing of bulked segregants identified a novel TaMKK3-A allele linked to the wheat ERA8 ABA-hypersensitive germination phenotype

Using bulked segregant analysis of exome sequence, we fine-mapped the ABA-hypersensitive mutant ERA8 in a wheat backcross population to the TaMKK3-A locus of chromosome 4A. Preharvest sprouting (PHS) is the germination of mature grain on the mother plant when it rains before harvest. The ENHANCED RESPONSE TO ABA8 (ERA8) mutant increases seed dormancy and, consequently, PHS tolerance in soft white wheat 'Zak.' ERA8 was mapped to chromosome 4A in a Zak/'ZakERA8' backcross population using bulked segregant analysis of exome sequenced DNA (BSA-exome-seq). ERA8 was fine-mapped relative to mutagen-induced SNPs to a 4.6 Mb region containing 70 genes. In the backcross population, the ERA8 ABA-hypersensitive phenotype was strongly linked to a missense mutation in TaMKK3-A-G1093A (LOD 16.5), a gene associated with natural PHS tolerance in barley and wheat. The map position of ERA8 was confirmed in an 'Otis'/ZakERA8 but not in a 'Louise'/ZakERA8 mapping population. This is likely because Otis carries the same natural PHS susceptible MKK3-A-A660^S allele as Zak, whereas Louise carries the PHS-tolerant MKK3-A-C660^R allele. Thus, the variation for grain dormancy and PHS tolerance in the Louise/ZakERA8 population likely resulted from segregation of other loci rather than segregation for PHS tolerance at the MKK3 locus. This inadvertent complementation test suggests that the MKK3-A-G1093A mutation causes the ERA8 phenotype. Moreover, MKK3 was a known ABA signaling gene in the 70-gene 4.6 Mb ERA8 interval. None of these 70 genes showed the differential regulation in wild-type Zak versus ERA8 expected of a promoter mutation. Thus, the working model is that the ERA8 phenotype results from the MKK3-A-G1093A mutation.