Genomewide association study for seeding emergence and tiller number using SNP markers in an elite winter wheat population

Quantitative trait loci detection for three tiller-related traits and the effects on wheat (Triticum aestivum L.) yields

KEY MESSAGE

A total of 38 putative additive QTLs and 55 pairwise putative epistatic QTLs for tiller-related traits were reported, and the candidate genes underlying qMtn-KJ-5D, a novel major and stable QTL for maximum tiller number, were characterized. Tiller-related traits play an important role in determining the yield potential of wheat. Therefore, it is important to elucidate the genetic basis for tiller number when attempting to use genetic improvement as a tool for enhancing wheat yields. In this study, a quantitative trait locus (QTL) analysis of three tiller-related traits was performed on the recombinant inbred lines (RILs) of a mapping population, referred to as KJ-RILs, that was derived from a cross between the Kenong 9204 (KN9204) and Jing 411 (J411) lines. A total of 38 putative additive QTLs and 55 pairwise putative epistatic QTLs for spike number per plant (SNPP), maximum tiller number (MTN), and ear-bearing tiller rate (EBTR) were detected in eight different environments. Among these QTLs with additive effects, three major and stable QTLs were first documented herein. Almost all but two pairwise epistatic QTLs showed minor interaction effects accounting for no more than 3.0% of the phenotypic variance. The genetic effects of two colocated major and stable QTLs, i.e., qSnpp-KJ-5D.1 and qMtn-KJ-5D, for yield-related traits were characterized. The breeding selection effect of the beneficial allele for the two QTLs was characterized, and its genetic effects on yield-related traits were evaluated. The candidate genes underlying qMtn-KJ-5D were predicted based on multi-omics data, and TraesKN5D01HG00080 was identified as a likely candidate gene. Overall, our results will help elucidate the genetic architecture of tiller-related traits and can be used to develop novel wheat varieties with high yields.

Identification of molecular markers and candidate regions associated with grain number per spike in Pubing3228 using SLAF-BSA

Grain number per spike, a pivotal agronomic trait dictating wheat yield, lacks a comprehensive understanding of its underlying mechanism in Pubing3228, despite the identification of certain pertinent genes. Thus, our investigation sought to ascertain molecular markers and candidate regions associated with grain number per spike through a high-density genetic mapping approach that amalgamates site-specific amplified fragment sequencing (SLAF-seq) and bulked segregation analysis (BSA). To facilitate this, we conducted a comparative analysis of two wheat germplasms, Pubing3228 and Jing4839, known to exhibit marked discrepancies in spike shape. By leveraging this methodology, we successfully procured 2,810,474 SLAF tags, subsequently resulting in the identification of 187,489 single nucleotide polymorphisms (SNPs) between the parental strains. We subsequently employed the SNP-index association algorithm alongside the extended distribution (ED) association algorithm to detect regions associated with the trait. The former algorithm identified 24 trait-associated regions, whereas the latter yielded 70. Remarkably, the intersection of these two algorithms led to the identification of 25 trait-associated regions. Amongst these regions, we identified 399 annotated genes, including three genes harboring non-synonymous mutant SNP loci. Notably, the APETALA2 (AP2) transcription factor families, which exhibited a strong correlation with spike type, were also annotated. Given these findings, it is plausible to hypothesize that these genes play a critical role in determining spike shape. In summation, our study contributes significant insights into the genetic foundation of grain number per spike. The molecular markers and candidate regions we have identified can be readily employed for marker-assisted breeding endeavors, ultimately leading to the development of novel wheat cultivars possessing enhanced yield potential. Furthermore, conducting further functional analyses on the identified genes will undoubtedly facilitate a comprehensive elucidation of the underlying mechanisms governing spike development in wheat.

Genome-wide association analysis of type II resistance to Fusarium head blight in common wheat

Background

Fusarium head blight (FHB) is a disease affecting wheat spikes caused by some Fusarium species and leads to cases of severe yield reduction and seed contamination. Identifying resistance genes/QTLs from wheat germplasm may help to improve FHB resistance in wheat production.

Methods

Our study evaluated 205 elite winter wheat cultivars for FHB resistance. A high-density 90K SNP array was used for genotyping the panel. A genome-wide association study (GWAS) from cultivars from three different environments was performed using a mixed linear model (MLM).

Results

Sixty-six significant marker-trait associations (MTAs) were identified (P < 0.001) on fifteen chromosomes that explained the phenotypic variation ranging from 5.4 to 11.2%. Some important new MTAs in genomic regions involving FHB resistance were found on chromosomes 2A, 3B, 5B, 6A, and 7B. Six MTAs at 92 cM on chromosome 7B were found in cultivars from two different environments. Moreover, there were 11 MTAs consistently associated with diseased spikelet rate and diseased rachis rate as pleiotropic effect loci and D_contig74317_533 on chromosome 5D was novel for FHB resistance. Eight new candidate genes of FHB resistance were predicated in wheat in this study. Three candidate genes, TraesCS5D02G006700, TraesCS6A02G013600, and TraesCS7B02G370700 on chromosome 5DS, 6AS, and 7BL, respectively, were perhaps important in defending against FHB by regulating intramolecular transferase activity, GTP binding, or chitinase activity in wheat, but further validation in needed. In addition, a total of five favorable alleles associated with wheat FHB resistance were discovered. These results provide important genes/loci for enhancing FHB resistance in wheat breeding by marker-assisted selection.

Genome-wide association study to identify genomic loci associated with early vigor in bread wheat under simulated water deficit complemented with quantitative trait loci meta-analysis

A genome-wide association study (GWAS) was used to identify associated loci with early vigor under simulated water deficit and grain yield under field drought in a diverse collection of Iranian bread wheat landraces. In addition, a meta-quantitative trait loci (MQTL) analysis was used to further expand our approach by retrieving already published quantitative trait loci (QTL) from recombinant inbred lines, double haploids, back-crosses, and F2 mapping populations. In the current study, around 16%, 14%, and 16% of SNPs were in significant linkage disequilibrium (LD) in the A, B, and D genomes, respectively, and varied between 5.44% (4A) and 21.85% (6A). Three main subgroups were identified among the landraces with different degrees of admixture, and population structure was further explored through principal component analysis. Our GWAS identified 54 marker-trait associations (MTAs) that were located across the wheat genome but with the highest number found in the B sub-genome. The gene ontology (GO) analysis of MTAs revealed that around 75% were located within or closed to protein-coding genes. In the MQTL analysis, 23 MQTLs, from a total of 215 QTLs, were identified and successfully projected onto the reference map. MQT-YLD4, MQT-YLD9, MQT-YLD13, MQT-YLD17, MQT-YLD18, MQT-YLD19, and MQTL-RL1 contributed to the highest number of projected QTLs and were therefore regarded as the most reliable and stable QTLs under water deficit conditions. These MQTLs greatly facilitate the identification of putative candidate genes underlying at each MQTL interval due to the reduced confidence of intervals associated with MQTLs. These findings provide important information on the genetic basis of early vigor traits and grain yield under water deficit conditions and set the foundation for future investigations into adaptation to water deficit in bread wheat.

Quantitative trait loci identification and breeding value estimation of grain weight-related traits based on a new wheat 50K single nucleotide polymorphism array-derived genetic map

Mining novel and less utilized thousand grain weight (TGW) related genes are useful for improving wheat yield. In this study, a recombinant inbred line population from a cross between Zhongkemai 138 (ZKM138, high TGW) and Chuanmai 44 (CM44, low TGW) was used to construct a new Wheat 50K SNP array-derived genetic map that spanned 1,936.59 cM and contained 4, 139 markers. Based on this map, ninety-one quantitative trait loci (QTL) were detected for eight grain-related traits in six environments. Among 58 QTLs, whose superior alleles were contributed by ZKM138, QTgw.cib-6A was a noticeable major stable QTL and was also highlighted by bulked segregant analysis with RNA sequencing (BSR-Seq). It had a pyramiding effect on TGW enhancement but no significant trade-off effect on grain number per spike or tiller number, with two other QTLs (QTgw.cib-2A.2 and QTgw.cib-6D), possibly explaining the excellent grain performance of ZKM138. After comparison with known loci, QTgw.cib-6A was deduced to be a novel locus that differed from nearby TaGW2 and TaBT1. Seven simple sequence repeat (SSR) and thirty-nine kompetitive allele-specific PCR (KASP) markers were finally developed to narrow the candidate interval of QTgw.cib-6A to 4.1 Mb. Only six genes in this interval were regarded as the most likely candidate genes. QTgw.cib-6A was further validated in different genetic backgrounds and presented 88.6% transmissibility of the ZKM138-genotype and a 16.4% increase of TGW in ZKM138 derivatives. And the geographic pattern of this locus revealed that its superior allele is present in only 6.47% of 433 Chinese modern wheat varieties, indicating its potential contribution to further high-yield breeding.

Four-Year and Five-Developing-Stage Dynamic QTL Mapping for Tiller Number in the Hybrid Population of Agropyron Gaertn

Tiller number (TN) is an important agronomic trait affecting gramineous crop yield. To understand the static and dynamic information of quantitative trait locus (QTLs) controlling TN of Agropyron Gaertn., both the unconditional and conditional quantitative trait loci (QTL) mapping of TN were conducted using a cross-pollinated (CP) hybrid population with a total of 113 plant lines from the cross between Agropyron cristatum (L.) Gaertn. Z1842 and Allium mongolicum Keng Z2098, based on the phenotypic data of TN at five developmental stages [i.e., recovering stage (RS), jointing stage (JS), heading stage (HS), flowering stage (FS), and maturity stage (MS)] in 4 years (i.e., 2017, 2018, 2020, and 2021) and the genetic map constructed of 1,023 single-nucleotide polymorphism (SNP) markers. Thirty-seven QTLs controlling TN were detected using two analysis methods in 4 years, which were distributed in six linkage groups. Each QTL explained 2.96-31.11% of the phenotypic variation, with a logarithum of odds (LOD) value of 2.51-13.95. Nine of these loci detected both unconditional and conditional QTLs. Twelve unconditional major QTLs and sixteen conditional major QTLs were detected. Three relatively major stable conditional QTLs, namely, cQTN1-3, cQTN1-5, and cQTN4-1, were expressed in 2020 and 2021. Meantime, two pairs of major QTLs cQTN1-5 and qTN1-4 and also cQTN2-4 and qTN2-3 were located at the same interval but in different years. Except for qTN2-2 and qTN3-5/cQTN3-5, other thirty-four QTLs were first detected in this study. This study provides a better interpretation of genetic factors that selectively control tiller at different developmental stages and a reference for molecular marker-assisted selection in the related plant improvement.

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.

Genome wide association study of plant height and tiller number in hulless barley

Hulless barley (Hordeum vulgare L. var. nudum), also called naked barley, is a unique variety of cultivated barley. The genome-wide specific length amplified fragment sequencing (SLAF-seq) method is a rapid deep sequencing technology that is used for the selection and identification of genetic loci or markers. In this study, we collected 300 hulless barley accessions and used the SLAF-seq method to identify candidate genes involved in plant height (PH) and tiller number (TN). We obtained a total of 1407 M paired-end reads, and 228,227 SLAF tags were developed. After filtering using an integrity threshold of >0.8 and a minor allele frequency of >0.05, 14,504,892 single-nucleotide polymorphisms (SNP) loci were screened out. The remaining SNPs were used for the construction of a neighbour-joining phylogenetic tree, and the three subcluster members showed no obvious differentiation among regional varieties. We used a genome wide association study approach to identify 1006 and 113 SNPs associated with TN and PH, respectively. Based on best linear unbiased predictors (BLUP), 41 and 29 SNPs associated with TN and PH, respectively. Thus, several of genes, including Hd3a and CKX5, may be useful candidates for the future genetic breeding of hulless barley. Taken together, our results provide insight into the molecular mechanisms controlling barley architecture, which is important for breeding and yield.

Physical mapping of QTL associated with agronomic and end-use quality traits in spring wheat under conventional and organic management systems

Using phenotypic data of four biparental spring wheat populations evaluated at multiple environments under two management systems, we discovered 152 QTL and 22 QTL hotspots, of which two QTL accounted for up to 37% and 58% of the phenotypic variance, consistently detected in all environments, and fell within genomic regions harboring known genes. Identification of the physical positions of quantitative trait loci (QTL) would be highly useful for developing functional markers and comparing QTL results across multiple independent studies. The objectives of the present study were to map and characterize QTL associated with nine agronomic and end-use quality traits (tillering ability, plant height, lodging, grain yield, grain protein content, thousand kernel weight, test weight, sedimentation volume, and falling number) in hard red spring wheat recombinant inbred lines (RILs) using the International Wheat Genome Sequencing Consortium (IWGSC) RefSeq v2.0 physical map. We evaluated a total of 698 RILs from four populations derived from crosses involving seven parents at 3-8 conventionally (high N) and organically (low N) managed field environments. Using the phenotypic data combined across all environments per management, and the physical map between 1058 and 6526 markers per population, we identified 152 QTL associated with the nine traits, of which 29 had moderate and 2 with major effects. Forty-nine of the 152 QTL mapped across 22 QTL hotspot regions with each region coincident to 2-6 traits. Some of the QTL hotspots were physically located close to known genes. QSv.dms-1A and QPht.dms-4B.1 individually explained up to 37% and 58% of the variation in sedimentation volume and plant height, respectively, and had very large LOD scores that varied from 19.0 to 35.7 and from 16.7 to 55.9, respectively. We consistently detected both QTL in the combined and all individual environments, laying solid ground for further characterization and possibly for cloning.

Identification of novel genomic regions associated with nine mineral elements in Chinese winter wheat grain

BACKGROUND

Mineral elements are important for maintaining good human health besides heavy metals. Mining genes that control mineral elements are paramount for improving their accumulation in the wheat grain. Although previous studies have reported some loci for beneficial trace elements, they have mainly focused on Zn and Fe content. However, little information is available regarding the genetic loci differences in dissecting synchronous accumulation of multiple mineral elements in wheat grains, including beneficial and heavy elements. Therefore, a genome-wide association study (GWAS) was conducted on 205 wheat accessions with 24,355 single nucleotide polymorphisms (SNPs) to identify important loci and candidate genes for controlling Ca, Fe, Zn, Se, Cu, Mn, Cd, As, and Pb accumulation in wheat grains.

RESULTS

A total of 101 marker-trait associations (MTAs) (P < 10^-5) loci affecting the content of nine mineral elements was identified on chromosomes 1B, 1D, 2A, 2B, 3A, 3B, 3D, 4A, 4B, 5A, 5B, 5D, 6B, 7A, 7B, and 7D. Among these, 17 major MTAs loci for the nine mineral elements were located, and four MTAs loci (P < 10^-5) were found on chromosomes 1B, 6B, 7B, and 7D. Eight multi-effect MTAs loci were detected that are responsible for the control of more than one trait, mainly distributed on chromosomes 3B, 7B, and 5A. Furthermore, sixteen candidate genes controlling Ca, Fe, Zn, Se, Cd, and Pb were predicted, whose functions were primarily related to ion binding, including metals, Fe, Ca, Cu, Mg, and Zn, ATP binding, ATPase activity, DNA binding, RNA binding, and protein kinase activity.

CONCLUSIONS

Our study indicated the existence of gene interactions among mineral elements based on multi-effect MTAs loci and candidate genes. Meanwhile this study provided new insights into the genetic control of mineral element concentrations, and the important loci and genes identified may contribute to the rapid development of beneficial mineral elements and a reduced content of harmful heavy metals in wheat grain.

Genome wide association study of the whiteness and colour related traits of flour and dough sheets in common wheat

Flour whiteness and colour are important factors that influence the quality of wheat flour and end-use products. In this study, a genome wide association study focusing on flour and dough sheet colour using a high density genetic map constructed with 90K single nucleotide polymorphism arrays in a panel of 205 elite winter wheat accessions was conducted in two different locations in 2 years. Eighty-six significant marker-trait associations (MTAs) were detected for flour whiteness and the brightness index (L* value), the redness index (a* value), and the yellowness index (b* value) of flour and dough sheets (P < 10^-4) on homologous group 1, 2, 5 and 7, and chromosomes 3A, 3B, 4A, 6A and 6B. Four, three, eleven, eleven MTAs for the flour whiteness, L* value, a* value, b* value, and one MTA for the dough sheet L* value were identified in more than one environment. Based on MATs, some important new candidate genes were identified. Of these, two candidate genes, TraesCS5D01G004300 and Gsp-1D, for BS00000020_51 were found in wheat, relating to grain hardness. Other candidate genes were associated with proteins, the fatty acid biosynthetic process, the ketone body biosynthetic process, etc.

QTL Analysis of Five Morpho-Physiological Traits in Bread Wheat Using Two Mapping Populations Derived from Common Parents

Traits such as plant height (PH), juvenile growth habit (GH), heading date (HD), and tiller number are important for both increasing yield potential and improving crop adaptation to climate change. In the present study, these traits were investigated by using the same bi-parental population at early (F₂ and F₂-derived F₃ families) and late (F₆ and F₇, recombinant inbred lines, RILs) generations to detect quantitative trait loci (QTLs) and search for candidate genes. A total of 176 and 178 lines were genotyped by the wheat Illumina 25K Infinium SNP array. The two genetic maps spanned 2486.97 cM and 3732.84 cM in length, for the F₂ and RILs, respectively. QTLs explaining the highest phenotypic variation were found on chromosomes 2B, 2D, 5A, and 7D for HD and GH, whereas those for PH were found on chromosomes 4B and 4D. Several QTL detected in the early generations (i.e., PH and tiller number) were not detected in the late generations as they were due to dominance effects. Some of the identified QTLs co-mapped to well-known adaptive genes (i.e., Ppd-1, Vrn-1, and Rht-1). Other putative candidate genes were identified for each trait, of which PINE1 and PIF4 may be considered new for GH and TTN in wheat. The use of a large F₂ mapping population combined with NGS-based genotyping techniques could improve map resolution and allow closer QTL tagging.

Genetic analyses of nitrogen assimilation enzymes in Brassica juncea (L.) Czern & Coss

Nitrogen (N) is a critical input for plant growth and development. A better understanding of N uptake and utilization is important to develop plant breeding strategies for improving nitrogen use efficiency (NUE). With that objective in mind, we assayed a SNP-genotyped association panel comprising 92 inbred lines for the activities of nitrate reductase (NR), nitrite reductase (NIR), glutamine synthetase (GS) and glutamate synthase (GOGAT). All these enzymes are associated with N assimilation. The experiments were carried out at two levels of N application: no added N (N₀) and agrnomically recommened dose (100 kg/ha) of N application (N₁₀₀). Genome wide association studies (GWAS) helped to identify several marker-trait associations (MTAs), involving chromosomes A01, A06, A08, B02, B04, B05 and B08. These explained high phenotypic variation (up to 32%). Annotation of the genomic region(s) in or around significant SNPs allowed prediction of genes encoding high affinity nitrate transporters, glutamine synthetase 1.3, myb-like transcription factor family protein, bidirectional amino acid transporter 1, auxin signaling F-box 3 and oxidoreductases. This is the first attempt to use GWAS for identification of enzyme QTLs to explain variation for nitrogen assimilation enzymes in Brassica juncea.

Genome-wide association studies of seven agronomic traits under two sowing conditions in bread wheat

BACKGROUND

Wheat is a cool seasoned crop requiring low temperature during grain filling duration and therefore increased temperature causes significant yield reduction. A set of 125 spring wheat genotypes from International Maize and Wheat Improvement Centre (CIMMYT-Mexico) was evaluated for phenological and yield related traits at three locations in Pakistan under normal sowing time and late sowing time for expose to prolonged high temperature. With the help of genome-wide association study using genotyping-by-sequencing, marker trait associations (MTAs) were observed separately for the traits under normal and late sown conditions.

RESULTS

Significant reduction ranging from 9 to 74% was observed in all traits under high temperature. Especially 30, 25, 41 and 66% reduction was observed for days to heading (DH), plant height (PH), spikes per plant (SPP) and yield respectively. We identified 55,954 single nucleotide polymorphisms (SNPs) using genotyping by sequencing of these 125 hexaploid spring wheat genotypes and conducted genome-wide association studies (GWAS) for days to heading (DH), grain filled duration (GFD), plant height (PH), spikes per plant (SPP), grain number per spike (GNS), thousand kernel weight (TKW) and grain yield per plot (GY). Genomic regions identified through GWAS explained up to 13% of the phenotypic variance, on average. A total of 139 marker-trait associations (MTAs) across three wheat genomes (56 on genome A, 55 on B and 28 on D) were identified for all the seven traits studied. For days to heading, 20; grain filled duration, 21; plant height, 23; spikes per plant, 13; grain numbers per spike, 8; thousand kernel weight, 21 and for grain yield, 33 MTAs were detected under normal and late sown conditions.

CONCLUSIONS

This study identifies the essential resource of genetics research and underpins the chromosomal regions of seven agronomic traits under normal and high temperature. Significant relationship was observed between the number of favored alleles and trait observations. Fourteen protein coding genes with their respective annotations have been searched with the sequence of seven MTAs which were identified in this study. These findings will be helpful in the development of a breeder friendly platform for the selection of high yielding wheat lines at high temperature areas.

Utilization of a Wheat55K SNP Array for Mapping of Major QTL for Temporal Expression of the Tiller Number

Maximum tiller number and productive tiller number are important traits for wheat grain yield, but research involving the temporal expression of tiller number at different quantitative trait loci (QTL) levels is limited. In the present study, a population set of 371 recombined inbred lines derived from a cross between Chuan-Nong18 and T1208 was used to construct a high-density genetic map using a Wheat55K SNP Array and to perform dynamic QTL analysis of the tiller number at four growth stages. A high-density genetic map containing 11,583 SNP markers and 59 SSR markers that spanned 4,513.95 cM and was distributed across 21 wheat chromosomes was constructed. A total of 28 single environmental QTL were identified in the recombined inbred lines population, and among these, seven QTL were stable and used for multi-environmental and dynamic analysis. These QTL were mapped to chromosomes 2D, 4A, 4D, 5A, 5D, and 7D, respectively. Each QTL explained 1.63-21.22% of the observed phenotypic variation, with an additive effect from -20.51 to 11.59. Dynamic analysis showed that cqTN-2D.2 can be detected at four growth stages of tillering, explaining 4.92-17.16% of the observed phenotypic variations and spanning 13.71 Mb (AX-109283238-AX-110544009: 82189047-95895626) according to the physical location of the flanking markers. The effects of the stable QTL were validated in the recombined inbred lines population, and the beneficial alleles could be utilized in future marker-assisted selection. Several candidate genes for MTN and PTN were predicted. The results provide a better understanding of the QTL selectively expressing the control of tiller number and will facilitate future map-based cloning. 9.17% SNP markers showed best hits to the Chinese Spring contigs. It was indicated that Wheat55K Array was efficient and valid to construct a high-density wheat genetic map.