Phenotypic and genetic analysis of spike and kernel characteristics in wheat reveals long-term genetic trends of grain yield components

Association analysis of response to take-all disease with agronomic traits and molecular markers and selection ideal genotypes in bread wheat (Triticum aestivum L.) genotypes

Take-all disease, caused by the fungal pathogen Gaeumannomyces tritici, severely impacts the growth and grain yield of wheat. Identifying loci associated with disease resistance can be achieved through molecular methods, along with data on morphological traits and disease severity. This study analyzed 100 bread wheat genotypes using molecular markers (SSR, IRAP technique, and translocation wheat-rye) and agronomical traits to pinpoint loci related to resistance to take-all disease. In this research, we propose a new approach using TOPSIS method for identifying ideal genotypes with resistance to take-all disease and the best in point of other agronomic traits. Genotypes were grouped based on agronomical traits (yield and its components) observed in the field, as well as root weight characteristics, root lignin content, and disease severity. These groupings effectively distinguished between resistant and sensitive genotypes. Stepwise regression techniques unveiled significant loci linked to disease resistance and agronomical traits. The presence of common loci suggests a potential pleiotropic nature of disease resistance. Molecular analysis exposed interactive loci contributing to trait variations and disease resistance, indicating gene-by-gene interactions. Using the IRAP technique, a locus from the LTR retrotransposon marker (LTR14) showed a strong correlation with take-all disease resistance and agronomic traits. This marker can serve as an informative and promising candidate for marker-assisted selection in wheat breeding programs. The TOPSIS method assisted in identifying genotypes showing high yield and resistance to take-all disease.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-025-01554-4.

Enhancing drought resilience in durum wheat: effect of root architecture and genotypic performance in semi-arid rainfed regions

Background

Developing drought-adapted genotypes is a primary goal for achieving resilient agriculture in the Mediterranean region. Durum wheat, a widely grown crop in the drylands of the Mediterranean basin, would significantly benefit from increased drought resistance.

Methods

We investigated a diverse set of 30 durum wheat varieties, including both local landraces and modern cultivars that have proven successful in Algeria. These varieties were evaluated in field trials over two consecutive years with contrasting rainfall patterns (one very dry, the other quite wet). Grain yield (PGY), yield components, and flag leaf characteristics such as area, canopy temperature, or rolling index were evaluated. Data from previous studies of root traits recorded on the same set of genotypes at seedling and adult growth stages were used to search for possible associations with grain yield and other agronomic traits measured in the current work.

Results

Genotypic variation was found for all traits measured under both conditions. Grain yield and aerial biomass were reduced by 76% (from 5.28 to 1.97 Mg ha-1) and 66% (from 15.94 to 3.80 Mg ha-1), respectively in the dry year, whereas the harvest index increased by 32%. The breeding history of the germplasm (cultivar vs. landrace) had a significant effect on the traits studied. Landraces showed higher biomass only under drought (4.27 vs. 3.63 Mg ha-1), whereas modern cultivars out-yielded landraces only under non-drought conditions (5.56 vs. 4.49 Mg ha-1). Promising associations were found between root and agronomic traits, especially with grain yield, indicating that a profuse (large root length) and shallow (wide root angle) root system was related to increased yield of modern cultivars only in the dry year, without penalizing yield in the wet year.

Conclusion

Breeding programs could improve grain yield under Algerian, semi-arid conditions, by making crosses between selected landraces with good growth potential under drought and modern cultivars, with high efficiency of biomass conversion into grain, and searching for lines with acceptable agronomic performance, which combine these desirable traits from landraces and modern cultivars, with the presence of shallow and profuse root systems.

Meta-QTL mapping for wheat thousand kernel weight

Wheat domestication and subsequent genetic improvement have yielded cultivated species with larger seeds compared to wild ancestors. Increasing thousand kernel weight (TKW) remains a crucial goal in many wheat breeding programs. To identify genomic regions influencing TKW across diverse genetic populations, we performed a comprehensive meta-analysis of quantitative trait loci (MQTL), integrating 993 initial QTL from 120 independent mapping studies over recent decades. We refined 242 loci into 66 MQTL, with an average confidence interval (CI) 3.06 times smaller than that of the original QTL. In these 66 MQTL regions, a total of 4,913 candidate genes related to TKW were identified, involved in ubiquitination, phytohormones, G-proteins, photosynthesis, and microRNAs. Expression analysis of the candidate genes showed that 95 were specific to grain and might potentially affect TKW at different seed development stages. These findings enhance our understanding of the genetic factors associated with TKW in wheat, providing reliable MQTL and potential candidate genes for genetic improvement of this trait.

Genetic improvement of important agronomic traits in Chinese wheat breeding over the past 70 years

BACKGROUND

Understanding the genetic improvement patterns of agronomic traits in Chinese wheat (Triticum aestivum L.) breeding is essential for devising future breeding strategies. However, a systematic analysis of the genetic improvement of important traits in Chinese wheat is lacking. This study aimed to provide insights into the improvement progress of yield-related traits in the Chinese wheat breeding process and clarify the selection pressure on important agronomic traits in different agroecological zones. Phenotypic evaluations of 481 wheat accessions including 157 Chinese landraces (CLs) and 324 modern Chinese cultivars (MCCs), were carried out in multiple locations and years.

RESULTS

The population structure analyses showed that all accessions could be basically divided into CLs and MCCs subpopulations. Pearson correlation analysis revealed that the negative correlation between grain number per spike and thousand-grain weight gradually decreased while thousand-grain weight, grain number per spike, and effective tiller number exhibited synergistic improvements during the modern breeding process. Phenotypic differences among MCCs released from the 1950s to the 2000s indicated that grain number per spike and grain weight-related traits increased linearly, whereas plant height and effective tiller number decreased significantly. Furthermore, since the 1950s, the heading date, flowering date, and maturity date have become earlier, while the spike length and spikelet number per spike have not changed significantly with the advancement of breeding years. The annual genetic gain analysis of agronomic traits showed that plant height had the greatest increase (‒0.96%), followed by thousand-grain weight (0.38%), while the lowest for grain number per spike (0.13%). Phenotypic difference analysis of CLs and MCCs with different geographical origins further revealed that heading date, flowering date, plant height, thousand-grain weight, grain width, and grain thickness experienced strong selection with the same trend in seven agroecological zones. Among zones, the northern winter wheat zone experienced the strongest selection pressure, and plant height and thousand-grain weight were strongly selected in all zones.

CONCLUSIONS

This study reveals that CLs and MCCs in China with obvious phenotypic differences, plant height and thousand-grain weight were strongly selected during wheat breeding, and further improvement of wheat in China will inevitably involve a continuous increase in grain number per spike.

Evaluation of the Spike Diversity of Seven Hexaploid Wheat Species and an Artificial Amphidiploid Using a Quadrangle Model Obtained from 2D Images

The spike shape and morphometric characteristics are among the key characteristics of cultivated cereals, being associated with their productivity. These traits are often used for the plant taxonomy and authenticity of hexaploid wheat species. Manual measurement of spike characteristics is tedious and not precise. Recently, the authors of this study developed a method for wheat spike morphometry utilizing 2D image analysis. Here, this method is applied to study variations in spike size and shape for 190 plants of seven hexaploid (2n = 6x = 42) species and one artificial amphidiploid of wheat. Five manually estimated spike traits and 26 traits obtained from digital image analysis were analyzed. Image-based traits describe the characteristics of the base, center and apex of the spike and common parameters (circularity, roundness, perimeter, etc.). Estimates of similar traits by manual measurement and image analysis were shown to be highly correlated, suggesting the practical importance of digital spike phenotyping. The utility of spike traits for classification into types (spelt, normal and compact) and species or amphidiploid is shown. It is also demonstrated that the estimates obtained made it possible to identify the spike characteristics differing significantly between species or between accessions within the same species. The present work suggests the usefulness of wheat spike shape analysis using an approach based on characteristics obtained by digital image analysis.

Identification and characterization of QSFS.sau-MC-5A for sterile florets genetically independent of fertile ones per spike in wheat

KEY MESSAGE

A major and stable QTL for sterile florets per spike and sterile florets per spikelet was identified, it was mapped within a 2.22-Mb interval on chromosome 5AL, and the locus was validated using two segregating populations with different genetic backgrounds. Both the number of fertile florets per spike (FFS) and the number of sterile florets per spike (SFS) significantly influence the final yield of wheat (Triticum aestivum L.), and a trade-off theoretically exists between them. To enhance crop yield, wheat breeders have historically concentrated on easily measurable traits such as FFS, spikelets per spike, and spike length. Other traits of agronomic importance, including SFS and sterile florets per spikelet (SFPs), have been largely overlooked. In the study, reported here, genetic bases of SFS and SFPs were investigated based on the assessment of a population of recombinant inbred lines (RILs) population. The RIL population was developed by crossing a spontaneous mutant with higher SFS (msf) with the cultivar Chuannong 16. A total of 10 quantitative trait loci (QTL) were identified, with QSFS.sau-MC-5A for SFS and QSFPs.sau-MC-5A for SFPs being the major and stable ones, and they were co-located on the long arm of chromosome 5A. The locus was located within a 2.22-Mb interval, and it was further validated in two additional populations based on a tightly linked Kompetitive Allele-Specific PCR (KASP) marker, K_sau_5A_691403852. Expression differences and promoter sequence variations were observed between the parents for both TraesCS5A03G1247300 and TraesCS5A03G1250300. The locus of QSFS.sau-MC-5A/QSFPs.sau-MC-5A showed a significantly positive correlation with spike length, florets in the middle spikelet, and total florets per spike, but it showed no correlation with either kernel number per spike (KNS) or kernel weight per spike. Appropriate nitrogen fertilizer application led to reduced SFS and increased KNS, supporting results from previous reports on the positive effect of nitrogen fertilizer on wheat spike and floret development. Based on these results, we propose a promising approach for breeding wheat cultivars with multiple fertile florets per spike, which could increase the number of kernels per spike and potentially improve yield. Collectively, these findings will facilitate further fine mapping of QSFS.sau-MC-5A/QSFPs.sau-MC-5A and be instrumental in strategies to increase KNS, thereby enhancing wheat yield.

Genetic dissection of value-added quality traits and agronomic parameters through genome-wide association mapping in bread wheat (T. aestivum L.)

Bread wheat (T. aestivum) is one of the world's most widely consumed cereals. Since micronutrient deficiencies are becoming more common among people who primarily depend upon cereal-based diets, a need for better-quality wheat varieties has been felt. An association panel of 154 T. aestivum lines was evaluated for the following quality traits: grain appearance (GA) score, grain hardness (GH), phenol reaction (PR) score, protein percent, sodium dodecyl sulfate (SDS) sedimentation value, and test weight (TWt). In addition, the panel was also phenotyped for grain yield and related traits such as days to heading, days to maturity, plant height, and thousand kernel weight for the year 2017-18 at the Borlaug Institute for South Asia (BISA) Ludhiana and Jabalpur sites. We performed a genome-wide association analysis on this panel using 18,351 genotyping-by-sequencing (GBS) markers to find marker-trait associations for quality and grain yield-related traits. We detected 55 single nucleotide polymorphism (SNP) marker trait associations (MTAs) for quality-related traits on chromosomes 7B (10), 1A (9), 2A (8), 3B (6), 2B (5), 7A (4), and 1B (3), with 3A, 4A, and 6D, having two and the rest, 4B, 5A, 5B, and 1D, having one each. Additionally, 20 SNP MTAs were detected for yield-related traits based on a field experiment conducted in Ludhiana on 7D (4) and 4D (3) chromosomes, while 44 SNP MTAs were reported for Jabalpur on chromosomes 2D (6), 7A (5), 2A (4), and 4A (4). Utilizing these loci in marker-assisted selection will benefit from further validation studies for these loci to improve hexaploid wheat for better yield and grain quality.

Analysis of the Effects of the Vrn-1 and Ppd-1 Alleles on Adaptive and Agronomic Traits in Common Wheat (Triticum aestivum L.)

Wheat heading time is primarily governed by two loci: VRN-1 (response to vernalization) and PPD-1 (response to photoperiod). Five sets of near-isogenic lines (NILs) were studied with the aim of investigating the effect of the aforementioned genes on wheat vegetative period duration and 14 yield-related traits. Every NIL was sown in the hydroponic greenhouse of the Institute of Cytology and Genetics, SB RAS. To assess their allelic composition at the VRN-1 and PPD-1 loci, molecular markers were used. It was shown that HT in plants with the Vrn-A1vrn-B1vrn-D1 genotype was reduced by 29 and 21 days (p < 0.001) in comparison to HT in plants with the vrn-A1Vrn-B1vrn-D1 and the vrn-A1vrn-B1Vrn-D1 genotypes, respectively. In our study, we noticed a decrease in spike length as well as spikelet number per spike parameter for some NIL carriers of the Vrn-A1a allele in comparison to carriers of the Vrn-B1 allele. PCA revealed three first principal components (PC), together explaining more than 70% of the data variance. Among the studied genetic traits, the Vrn-A1a and Ppd-D1a alleles showed significant correlations with PCs. Regarding genetic components, significant correlations were calculated between PC3 and Ppd-B1a (-0.26, p < 0.05) and Vrn-B1 (0.57, p < 0.05) alleles. Thus, the presence of the Vrn-A1a allele affects heading time, while Ppd-D1a is associated with plant height reduction.

Genome-wide association study and genomic selection of spike-related traits in bread wheat

KEY MESSAGE

Identification of 337 stable MTAs for wheat spike-related traits improved model accuracy, and favorable alleles of MTA259 and MTA64 increased grain weight and yield per plant. Wheat (Triticum aestivum L.) is one of the three primary global, staple crops. Improving spike-related traits in wheat is crucial for optimizing spike and plant morphology, ultimately leading to increased grain yield. Here, we performed a genome-wide association study using a dataset of 24,889 high-quality unique single-nucleotide polymorphisms (SNPs) and phenotypic data from 314 wheat accessions across eight diverse environments. In total, 337 stable and significant marker-trait associations (MTAs) related to spike-related traits were identified. MTA259 and MTA64 were consistently detected in seven and six environments, respectively. The presence of favorable alleles associated with MTA259 and MTA64 significantly reduced wheat spike exsertion length and spike length, while enhancing thousand kernel weight and yield per plant. Combined gene expression and network analyses identified TraesCS6D03G0692300 and TraesCS6D03G0692700 as candidate genes for MTA259 and TraesCS2D03G0111700 and TraesCS2D03G0112500 for MTA64. The identified MTAs significantly improved the prediction accuracy of each model compared with using all the SNPs, and the random forest model was optimal for genome selection. Additionally, the eight stable and major MTAs, including MTA259, MTA64, MTA66, MTA94, MTA110, MTA165, MTA180, and MTA164, were converted into cost-effective and efficient detection markers. This study provided valuable genetic resources and reliable molecular markers for wheat breeding programs.

Dissection and validation of a promising QTL controlling spikelet number on 5B in bread wheat

KEY MESSAGE

Five environmentally stable QTLs for spikelet number per spike and days to heading were identified using a high-genetic map containing 95,444 SNPs, among which QSns.ucas-5B was validated using residual heterozygous line at multiple environments. Spikelet number per spike (SNS) and days to heading (DTH) play pivotal roles in the improvement of wheat yield. In this study, a high-density genetic map for a recombinant inbred lines (RILs) population derived from Zhengnong 17 (ZN17) and Yangbaimai (YBM) was constructed using 95,444 single-nucleotide polymorphism (SNP) markers from the Wheat660K SNP array. Our study identified a total of five environmentally stable QTLs for SNS and DTH, one of which was named QSns.ucas-5B, with a physical interval of approximately 545.4-552.1 Mb on the 5BL chromosome arm. Importantly, the elite haplotype within QSns.ucas-5B showed a consistent and positive effect on SNS, grain number and weight per spike, without extending the days to heading. These findings provide a foundation for future efforts to map and clone the gene(s) responsible for QSns.ucas-5B and further indicate the potential application of the developed and validated InDel marker of QSns.ucas-5B for molecular breeding purposes, aimed at improving wheat grain yield.

Assessment of genetic variation among wheat genotypes for drought tolerance utilizing microsatellite markers and morpho-physiological characteristics

Drought is a major abiotic stress that severely limits sustainable wheat (Triticum aestivum L.) productivity via morphological and physio-biochemical alterations of cellular processes. The complex nature and polygenic control of drought tolerance traits make breeding tolerant genotypes quite challenging. However, naturally occurring variabilities among wheat germplasm resources could potentially help combating drought. The present study was conducted to assess the drought tolerance of 18 Bangladeshi hexaploid wheat genotypes, focusing on the identification of potent sources of diversity by combining microsatellite markers, also known as single sequence repeat markers, and morpho-physiological characteristics that might help accelerating wheat crop improvement programs. Initially, the genotypes were evaluated using 25 microsatellite markers followed by an on-field evaluation of 7 morphological traits (plant height, spike number, spike length, grains per spike, 1000-grain weight, grain yield, biological yield) and 6 physiological traits (SPAD value, membrane stability index, leaf relative water content, proline content, canopy temperature depression, and leaf K+ ion content). The field-trial was conducted in a factorial fashion of 18 wheat genotypes and two water regimes (control and drought) following a split-plot randomized complete block design. Regardless of genotype, drought was significantly damaging for all the tested traits; however, substantial variability in drought stress tolerance was evident among the genotypes. Spike length, 1000-grain weight, SPAD value, leaf relative water content, canopy temperature depression, proline content, and potassium (K+) ion content were the most representative of drought-induced growth and yield impairments and also correlated well with the contrasting ability of genotypic tolerance. Microsatellite markers amplified 244 alleles exhibiting 79% genetic diversity. Out of 25 markers, 23 was highly polymorphic showing 77% average polymorphism. Morpho-physiological trait-based hierarchical clustering and microsatellite marker-based neighbor-jointing clustering both revealed three genotypic clusters with 71% co-linearity between them. In both cases, the genotypes Kanchan, BAW-1147, BINA Gom 1, BARI Gom 22, BARI Gom 26, and BARI Gom 33 were found to be comparatively more tolerant than the other tested genotypes, showing potential for cultivation in water-deficit environments. The findings of this study would contribute to the present understanding of drought tolerance in wheat and would provide a basis for future genotype selection for drought-tolerant wheat breeding programs.

Predictor bias in genomic and phenomic selection

KEY MESSAGE

NIRS of wheat grains as phenomic predictors for grain yield show inflated prediction ability and are biased toward grain protein content. Estimating the breeding value of individuals using genome-wide marker data (genomic prediction) is currently one of the most important drivers of breeding progress in major crops. Recently, phenomic technologies, including remote sensing and aerial hyperspectral imaging of plant canopies, have made it feasible to predict the breeding value of individuals in the absence of genetic marker data. This is commonly referred to as phenomic prediction. Hyperspectral measurements in the form of near-infrared spectroscopy have been used since the 1980 s to predict compositional parameters of harvest products. Moreover, in recent studies NIRS from grains was used to predict grain yield. The same studies showed that phenomic prediction can outperform genomic prediction for grain yield. The genome is static and not environment dependent, thereby limiting genomic prediction ability. Gene expression is tissue specific and differs under environmental influences, leading to a tissue- and environment-specific phenome, potentially explaining the higher predictive ability of phenomic prediction. Here, we compare genomic prediction and phenomic prediction from hyperspectral measurements of wheat grains for the prediction of a variety of traits including grain yield. We show that phenomic predictions outperform genomic prediction for some traits. However, phenomic predictions are biased toward the information present in the predictor. Future studies on this topic should investigate whether population parameters are retained in phenomic prediction as they are in genomic prediction. Furthermore, we find that unbiased phenomic prediction abilities are considerably lower than previously reported and recommend a method to circumvent this issue.

Genome-Wide Association Study of Asian and European Common Wheat Accessions for Yield-Related Traits and Stripe Rust Resistance

Identifying novel loci of yield-related traits and resistance to stripe rust (caused by Puccinia striiformis f. sp. tritici) in wheat will help in breeding wheat that can meet projected demands in diverse environmental and agricultural practices. We performed a genome-wide association study with 24,767 single nucleotide polymorphisms (SNPs) in 180 wheat accessions that originated in 16 Asian or European countries between latitudes 30°N and 45°N. We detected seven accessions with desirable yield-related traits and 42 accessions that showed stable, high degrees of stripe rust resistance in multienvironment field assessments. A marker-trait association analysis of yield-related traits detected 18 quantitative trait loci (QTLs) in at least two test environments and two QTLs related to stripe rust resistance in at least three test environments. Five of these QTLs were identified as potentially novel QTLs by comparing their physical locations with those of known QTLs in the Chinese Spring (CS) reference genome RefSeq v1.1 published by the International Wheat Genome Sequencing Consortium; two were for spike length, one was for grain number per spike, one was for spike number, and one was for stripe rust resistance at the adult plant stage. We also identified 14 candidate genes associated with the five novel QTLs. These QTLs and candidate genes will provide breeders with new germplasm and can be used to conduct marker-assisted selection in breeding wheat with improved yield and stripe rust resistance.

Exploring the trade-off between individual fitness and community performance of wheat crops using simulated canopy shade

The genetic heritage of wheat (Triticum spp.) crops has been shaped by millions of years of predomestication natural selection, often driven by competition among individuals. However, genetic improvements in yield potential are thought to involve selection towards reduced competitiveness, thus enhancing adaptation to the crop environment. We investigated potential trade-offs between individual plant fitness and community performance using a population of introgression lines carrying chromosome segments of wild emmer (nondomesticated) in the background of an elite durum cultivar. We focused on light as a primary factor affecting plant-plant interactions and assessed morphological and biomass phenotypes of single plants grown in mixtures under sunlight and a simulated canopy shade, and the relevance of these phenotypes for the monoculture community in the field. We found that responses to canopy shade resemble responses to high density and contribute to both the individual and the community. Stepwise regressions suggested that grain number per spike and its persistence under shade are essential attributes of productive communities, advocating their use as a breeding target during early-generation selection. Overall, multiple phenotypes attained under shade could better explain community performance. Our novel, applicable, high-throughput set-up provides new prospects for studying and selecting single-plant phenotypes in a canopy-like environment.

Identification and Validation of a Stable Major-Effect Quantitative Trait Locus for Kernel Number per Spike on Chromosome 2D in Wheat (Triticum aestivum L.)

A recombinant inbred line population including 371 lines was developed by a high kernel number per spike (KNPS) genotype T1208 and a low KNPS genotype Chuannong18 (CN18). A genetic linkage map consisting of 11,583 markers was constructed by the Wheat55K SNP Array. The quantitative trait loci (QTLs) related to KNPS were detected in three years. Eight, twenty-seven, and four QTLs were identified using the ICIM-BIP, ICIM-MET, and ICIM-EPI methods, respectively. One QTL, QKnps.sau-2D.1, which was mapped on chromosome 2D, can explain 18.10% of the phenotypic variation (PVE) on average and be considered a major and stable QTL for KNPS. This QTL was located in a 0.89 Mb interval on chromosome 2D and flanked by the markers AX-109283238 and AX-111606890. Moreover, KASP-AX-111462389, a Kompetitive Allele-Specific PCR (KASP) marker which closely linked to QKnps.sau-2D.1, was designed. The genetic effect of QKnps.sau-2D.1 on KNPS was successfully confirmed in two RIL populations. The results also showed that the significant increase of KNPS and 1000-kernel weight (TKW) was caused by QKnps.sau-2D.1 overcoming the disadvantage due to the decrease of spike number (SN) and finally lead to a significant increase of grain yield. In addition, within the interval in which QKnps.sau-2D.1 is located in Chinese Spring reference genomes, only fifteen genes were found, and two genes that might associate with KNPS were identified. QKnps.sau-2D.1 may provide a new resource for the high-yield breeding of wheat in the future.

Delayed development of basal spikelets in wheat explains their increased floret abortion and rudimentary nature

Large differences exist in the number of grains per spikelet across an individual wheat (Triticum aestivum L.) spike. The central spikelets produce the highest number of grains, while apical and basal spikelets are less productive, and the most basal spikelets are commonly only developed in rudimentary form. Basal spikelets are delayed in initiation, yet they continue to develop and produce florets. The precise timing or the cause of their abortion remains largely unknown. Here, we investigated the underlying causes of basal spikelet abortion using shading applications in the field. We found that basal spikelet abortion is likely to be the consequence of complete floret abortion, as both occur concurrently and have the same response to shading treatments. We detected no differences in assimilate availability across the spike. Instead, we show that the reduced developmental age of basal florets pre-anthesis is strongly associated with their increased abortion. Using the developmental age pre-abortion, we were able to predict final grain set per spikelet across the spike, alongside the characteristic gradient in the number of grains from basal to central spikelets. Future efforts to improve spikelet homogeneity across the spike could thus focus on improving basal spikelet establishment and increasing floret development rates pre-abortion.

A high-resolution genotype-phenotype map identifies the TaSPL17 controlling grain number and size in wheat

BACKGROUND

Large-scale genotype-phenotype association studies of crop germplasm are important for identifying alleles associated with favorable traits. The limited number of single-nucleotide polymorphisms (SNPs) in most wheat genome-wide association studies (GWASs) restricts their power to detect marker-trait associations. Additionally, only a few genes regulating grain number per spikelet have been reported due to sensitivity of this trait to variable environments.

RESULTS

We perform a large-scale GWAS using approximately 40 million filtered SNPs for 27 spike morphology traits. We detect 132,086 significant marker-trait associations and the associated SNP markers are located within 590 associated peaks. We detect additional and stronger peaks by dividing spike morphology into sub-traits relative to GWAS results of spike morphology traits. We propose that the genetic dissection of spike morphology is a powerful strategy to detect signals for grain yield traits in wheat. The GWAS results reveal that TaSPL17 positively controls grain size and number by regulating spikelet and floret meristem development, which in turn leads to enhanced grain yield per plant. The haplotypes at TaSPL17 indicate geographical differentiation, domestication effects, and breeding selection.

CONCLUSION

Our study provides valuable resources for genetic improvement of spike morphology and a fast-forward genetic solution for candidate gene detection and cloning in wheat.

An intron-located single nucleotide variation of TaGS5-3D is related to wheat grain size through accumulating intron retention transcripts

KEY MESSAGE

Thirty-three stable QTL for 13 yield-related traits across ten environments were identified in the PD34/MY47 RIL population, and a candidate gene TaGS5-3D in Qmt.nwafu.3D was preliminarily identified to affect grain-related traits through accumulation of specific transcripts. Dissecting the genetic basis of yield-related traits is pivotal for improvement of wheat yield potential. In this study, a recombinant inbred line (RIL) population genotyped by SNP markers was used to detect quantitative trait loci (QTL) related to yield-related traits in ten environments. A total of 102 QTL were detected, including 33 environmentally stable QTL and 69 putative QTL. Among them, Qmt.nwafu.3D was identified as a pleiotropic QTL in the physical interval of 149.77-154.11 Mb containing a potential candidate gene TaGS5-3D. An SNP (T > C) was detected in its ninth intron, and TaGS5-3D-C was validated as a superior allele associated with larger grains using a CAPS marker. Interestingly, we found that TaGS5-3D-C was closely related to significantly up-regulated expression of intron-retained transcript (TaGS5-3D-PD34.1), while TaGS5-3D-T was related to dominant expression of normal splicing transcript (TaGS5-3D-MY47.1). Our results indicated that alternative splicing associated with the SNP T/C could be involved in the regulation of grain-related traits, laying a foundation for the functional analysis of TaGS5-3D and its greater potential application in high-yield wheat breeding.

Identification and validation of new quantitative trait loci for spike-related traits in two RIL populations

Wheat (Triticum aestivum L.) is one of the most important cereal crops for ensuring food security worldwide. Identification of major quantitative trait loci (QTL) for spike-related traits is important for improvement of yield potential in wheat breeding. In this study, by using the wheat 55K single nucleotide polymorphism (SNP) array and diversity array technology (DArT), two recombinant inbred line populations derived from crosses avocet/chilero and avocet/huites were used to map QTL for kernel number per spike (KNS), total spikelet number per spike (TSS), fertile spikelet number per spike (FSS), and spike compactness (SC). Forty-two QTLs were identified on chromosomes 2A (4), 2B (3), 3A (2), 3B (7), 5A (11), 6A (4), 6B, and 7A (10), explaining 3.13-21.80% of the phenotypic variances. Twelve QTLs were detected in multi-environments on chromosomes 2A, 3B (2), 5A (4), 6A (3), 6B, and 7A, while four QTL clusters were detected on chromosomes 3A, 3B, 5A, and 7A. Two stable and new QTL clusters, QKns/Tss/Fss/SC.haust-5A and QKns/Tss/Fss.haust-7A, were detected in the physical intervals of 547.49-590.46 Mb and 511.54-516.15 Mb, accounting for 7.53-14.78% and 7.01-20.66% of the phenotypic variances, respectively. High-confidence annotated genes for QKns/Tss/Fss/SC.haust-5A and QKns/Tss/Fss.haust-7A were more highly expressed in spike development. The results provide new QTL and molecular markers for marker-assisted breeding in wheat.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-023-01401-4.

Performance and Stability Analysis of Selected Durum Wheat Genotypes Differing in Their Kernel Characteristics

Breeding of high-yielding and stable durum wheat varieties with improved kernel characteristics is needed for dry regions around the globe. The aim of this study was to investigate the performance and stability of eight durum wheat genotypes varying in their kernel characteristics across 15 contrasting environments. The tested material included three recombinant inbred lines (NUR-072, NUR-106 and NUR-238) derived from a cross between Norsi, a Jordanian landrace with special kernel characteristics and UC1113 Yr36+Gpc-B1, an elite line from USA. Field trials were carried out for three constitutive growing seasons under rainfed conditions, except for three environments where supplementary irrigation was provided. After the harvest, grain yield (GY), total yield (TW), and harvest index (HI) were recorded. Additionally, several kernel-related traits, including thousand kernel weight (TKW), kernel area (KA), kernel width (KW), kernel length (KL), kernel circularity (KC), and kernel length-width ratio (KL:KW) were evaluated. Analysis of variance for all tested traits revealed high significant variations (p ≤ 0.01) between the genotype (except for TW) and the genotype × environment (G × E) interaction. Genotype effect contributed to substantial percentage of variation (>75%) for KA, KL, KC and KL:KW, whereas KW showed a lower percentage similar to GY. Regarding the G × E effect, explained variation was highest for the TW (67.79%), and lowest for KL (6.47%). For GY, Norsi produced significantly the lowest mean value (249.99 g.m^-2) while, Bolenga produced the highest mean value (377.85 g.m^-2) although no significant differences were observed with the remaining genotypes. On the other hand, Norsi, NUR-072 and NUR-106 showed best performance for TKW and kernel-related traits with NUR-106 producing the highest mean value for KL (9.07 mm). The GGE biplot and AMMI analysis of GY identified Bolenga, Um Qais and NUR-106 as good performers across several environments, while Norsi exhibited the poorest performance. For TKW, Norsi was the best performer across different environments followed by NUR-106, which showed excellent performance under irrigated and saline conditions. For stability analysis, NUR-106 emerged as the most stable genotype in this study for GY and several kernel-related traits, particularly for KL and KC. In conclusion, the results of this study offer valuable insights for durum wheat breeders seeking to develop high-yielding and stable varieties with special kernel characteristics suitable for cultivation in dry areas.

Genomics for Yield and Yield Components in Durum Wheat

In recent years, many efforts have been conducted to dissect the genetic basis of yield and yield components in durum wheat thanks to linkage mapping and genome-wide association studies. In this review, starting from the analysis of the genetic bases that regulate the expression of yield for developing new durum wheat varieties, we have highlighted how, currently, the reductionist approach, i.e., dissecting the yield into its individual components, does not seem capable of ensuring significant yield increases due to diminishing resources, land loss, and ongoing climate change. However, despite the identification of genes and/or chromosomal regions, controlling the grain yield in durum wheat is still a challenge, mainly due to the polyploidy level of this species. In the review, we underline that the next-generation sequencing (NGS) technologies coupled with improved wheat genome assembly and high-throughput genotyping platforms, as well as genome editing technology, will revolutionize plant breeding by providing a great opportunity to capture genetic variation that can be used in breeding programs. To date, genomic selection provides a valuable tool for modeling optimal allelic combinations across the whole genome that maximize the phenotypic potential of an individual under a given environment.

Identification and validation of a major quantitative trait locus for spike length and compactness in the wheat (Triticum aestivum L.) line Chuanyu12D7

Spike length (SL) and spike compactness (SC) are crucial traits related to wheat (Triticum aestivum L.) yield potential. In this study, a backcrossed inbred lines (BILs) population segregating for SL/SC was developed by using a commercial variety chuanyu25 as recurrent parent and a backbone parent Chuanyu12D7. Bulked segregant analysis (BSA) combined with the Wheat 660K SNP array was performed to conduct quantitative trait locus (QTL) mapping. A major and stable SL/SC QTL (designated as QSl/Sc.cib-2D.1) was identified on chromosome 2DS, explaining 45.63-59.72% of the phenotypic variation. QSl/Sc.cib-2D.1 was mapped to a 102.29-Kb interval by flanking SNPs AX-110276364 and AX-111593853 using a BC₄F_2:3 population. Since QSl/Sc.cib-2D.1 is linked to the Rht8 gene, their additive effects on plant type and spike type were analysed. Remarkably, the superior allele of QSl/Sc.cib-2D.1 combined with Rht8 can increase SL and TGW, and decrese SC without any apparent trade-offs in other yield-related traits. In addition, the closely linked kompetitive allele-specific PCR (KASP) markers of this locus were developed for marker-assisted selection (MAS) breeding. Four genes within the physical interval were considered as potential candidates based on expression patterns as well as orthologous gene functions. These results laid the foundation for map-based cloning of the gene(s) underlying QSl/Sc.cib-2D.1 and its potential application in wheat ideotype breeding.

Chlorophyll Meter: A Precision Agricultural Decision-Making Tool for Nutrient Supply in Durum Wheat (Triticum turgidum L.) Cultivation under Drought Conditions

How crop biodiversity adapts to drought conditions and enhances grain yield became the most important issue facing agronomists and plant breeders at the turn of the century. Variations in genetic response, inadequacy of nutrients in the soil, and insufficient access to nutrients are factors that aggravate drought stressors. The development of screening tools for identifying drought tolerance is important in the deployment of durum wheat varieties suited to drought-prone environments. An experiment was conducted to evaluate durum wheat varieties under a range of nutrient supplies in naturally imposed drought conditions. The treatments consisting of two nitrogen regimes (i.e., control and 60 kg ha−1), four durum wheat varieties, and three types of nutrients (control, sulfur, and zinc) that were arranged in a split-split plot design with three replications. Both foliar-based sulfur and zinc fertilization were employed at the flag leaf stage, at a rate of 4 and 3-L ha−1, respectively. The results showed a significant (p < 0.05) genetic variation in chlorophyll concentration, grain protein content, tillering potential, and leaf area index. Varieties that contained better leaf chlorophyll content had improved grain yield by about 8.33% under 60 kg/ha nitrogen. A combined application of nitrogen and zinc at flag leaf stage significantly improved grain yield of Duragold by about 21.3%. Leaf chlorophyll content was found to be a more important trait than spikes per m2 to discriminate durum wheat varieties. Foliar application of sulfur increased the grain yield of drought-stressed plants by about 12.23%. Grain yield and protein content were strongly correlated with late-season SPAD readings. Significant (p < 0.05) correlation coefficients were obtained between normalized difference vegetation index, leaf area index, grain yield, and protein content with late-season chlorophyll content, revealing the importance of chlorophyll content in studying and identifying drought-tolerant varieties.

Amand P, Bai G, Baker J, Baker S, Liu S. QTL mapping of yield components and kernel traits in wheat cultivars TAM 112 and Duster

In the Southern Great Plains, wheat cultivars have been selected for a combination of outstanding yield and drought tolerance as a long-term breeding goal. To understand the underlying genetic mechanisms, this study aimed to dissect the quantitative trait loci (QTL) associated with yield components and kernel traits in two wheat cultivars `TAM 112' and `Duster' under both irrigated and dryland environments. A set of 182 recombined inbred lines (RIL) derived from the cross of TAM 112/Duster were planted in 13 diverse environments for evaluation of 18 yield and kernel related traits. High-density genetic linkage map was constructed using 5,081 single nucleotide polymorphisms (SNPs) from genotyping-by-sequencing (GBS). QTL mapping analysis detected 134 QTL regions on all 21 wheat chromosomes, including 30 pleiotropic QTL regions and 21 consistent QTL regions, with 10 QTL regions in common. Three major pleiotropic QTL on the short arms of chromosomes 2B (57.5 - 61.6 Mbps), 2D (37.1 - 38.7 Mbps), and 7D (66.0 - 69.2 Mbps) colocalized with genes Ppd-B1, Ppd-D1, and FT-D1, respectively. And four consistent QTL associated with kernel length (KLEN), thousand kernel weight (TKW), plot grain yield (YLD), and kernel spike^-1 (KPS) (Qklen.tamu.1A.325, Qtkw.tamu.2B.137, Qyld.tamu.2D.3, and Qkps.tamu.6A.113) explained more than 5% of the phenotypic variation. QTL Qklen.tamu.1A.325 is a novel QTL with consistent effects under all tested environments. Marker haplotype analysis indicated the QTL combinations significantly increased yield and kernel traits. QTL and the linked markers identified in this study will facilitate future marker-assisted selection (MAS) for pyramiding the favorable alleles and QTL map-based cloning.

Soil Fertility Clock-Crop Rotation as a Paradigm in Nitrogen Fertilizer Productivity Control

The Soil Fertility Clock (SFC) concept is based on the assumption that the critical content (range) of essential nutrients in the soil is adapted to the requirements of the most sensitive plant in the cropping sequence (CS). This provides a key way to effectively control the productivity of fertilizer nitrogen (Nf). The production goals of a farm are set for the maximum crop yield, which is defined by the environmental conditions of the production process. This target can be achieved, provided that the efficiency of Nf approaches 1.0. Nitrogen (in fact, nitrate) is the determining yield-forming factor, but only when it is balanced with the supply of other nutrients (nitrogen-supporting nutrients; N-SNs). The condition for achieving this level of Nf efficiency is the effectiveness of other production factors, including N-SNs, which should be set at ≤1.0. A key source of N-SNs for a plant is the soil zone occupied by the roots. N-SNs should be applied in order to restore their content in the topsoil to the level required by the most sensitive crop in a given CS. Other plants in the CS provide the timeframe for active controlling the distance of the N-SNs from their critical range.

Identification of stable QTL controlling multiple yield components in a durum × cultivated emmer wheat population under field and greenhouse conditions

Crop yield gains are needed to keep pace with a growing global population and decreasing resources to produce food. Cultivated emmer wheat is a progenitor of durum wheat and a useful source of genetic variation for trait improvement in durum. Here, we evaluated a recombinant inbred line population derived from a cross between the North Dakota durum wheat variety Divide and the cultivated emmer wheat accession PI 272527 consisting of 219 lines. The population was evaluated in 3 field environments and 2 greenhouse experiments to identify quantitative trait locus associated with 11 yield-related traits that were expressed in a consistent manner over multiple environments. We identified 27 quantitative trait locus expressed in at least 2 field environments, 17 of which were also expressed under greenhouse conditions. Seven quantitative trait locus regions on chromosomes 1B, 2A, 2B, 3A, 3B, 6A, and 7B had pleiotropic effects on multiple yield-related traits. The previously cloned genes Q and FT-B1, which are known to be associated with development and morphology, were found to consistently be associated with multiple traits across environments. PI 272527 contributed beneficial alleles for quantitative trait locus associated with multiple traits, especially for seed morphology quantitative trait locus on chromosomes 1B, 2B, and 6A. Three recombinant inbred lines with increased grain size and weight compared to Divide were identified and demonstrated the potential for improvement of durum wheat through deployment of beneficial alleles from the cultivated emmer parent. The findings from this study provide knowledge regarding stable and robust quantitative trait locus that breeders can use for improving yield in durum wheat.

Transcriptional signatures of wheat inflorescence development

In order to maintain global food security, it will be necessary to increase yields of the cereal crops that provide most of the calories and protein for the world's population, which includes common wheat (Triticum aestivum L.). An important wheat yield component is the number of grain-holding spikelets which form on the spike during inflorescence development. Characterizing the gene regulatory networks controlling the timing and rate of inflorescence development will facilitate the selection of natural and induced gene variants that contribute to increased spikelet number and yield. In the current study, co-expression and gene regulatory networks were assembled from a temporal wheat spike transcriptome dataset, revealing the dynamic expression profiles associated with the progression from vegetative meristem to terminal spikelet formation. Consensus co-expression networks revealed enrichment of several transcription factor families at specific developmental stages including the sequential activation of different classes of MIKC-MADS box genes. This gene regulatory network highlighted interactions among a small number of regulatory hub genes active during terminal spikelet formation. Finally, the CLAVATA and WUSCHEL gene families were investigated, revealing potential roles for TtCLE13, TtWOX2, and TtWOX7 in wheat meristem development. The hypotheses generated from these datasets and networks further our understanding of wheat inflorescence development.

A comprehensive meta-QTL analysis for yield-related traits of durum wheat (Triticum turgidum L. var. durum) grown under different water regimes

Abiotic stress strongly affects yield-related traits in durum wheat, in particular drought is one of the main environmental factors that have effect on grain yield and plant architecture. In order to obtain new genotypes well adapted to stress conditions, the highest number of desirable traits needs to be combined in the same genotype. In this context, hundreds of quantitative trait loci (QTL) have been identified for yield-related traits in different genetic backgrounds and environments. Meta-QTL (MQTL) analysis is a useful approach to combine data sets and for creating consensus positions for the QTL detected in independent studies for the reliability of their location and effects. MQTL analysis is a useful method to dissect the genetic architecture of complex traits, which provide an extensive allelic coverage, a higher mapping resolution and allow the identification of putative molecular markers useful for marker-assisted selection (MAS). In the present study, a complete and comprehensive MQTL analysis was carried out to identify genomic regions associated with grain-yield related traits in durum wheat under different water regimes. A total of 724 QTL on all 14 chromosomes (genomes A and B) were collected for the 19 yield-related traits selected, of which 468 were reported under rainfed conditions, and 256 under irrigated conditions. Out of the 590 QTL projected on the consensus map, 421 were grouped into 76 MQTL associated with yield components under both irrigated and rainfed conditions, 12 genomic regions containing stable MQTL on all chromosomes except 1A, 4A, 5A, and 6B. Candidate genes associated to MQTL were identified and an in-silico expression analysis was carried out for 15 genes selected among those that were differentially expressed under drought. These results can be used to increase durum wheat grain yields under different water regimes and to obtain new genotypes adapted to climate change.

Whole-genome analysis of hard winter wheat germplasm identifies genomic regions associated with spike and kernel traits

Genetic dissection of yield component traits including spike and kernel characteristics is essential for the continuous improvement in wheat yield. Genome-wide association studies (GWAS) have been frequently used to identify genetic determinants for spike and kernel-related traits in wheat, though none have been employed in hard winter wheat (HWW) which represents a major class in US wheat acreage. Further, most of these studies relied on assembled diversity panels instead of adapted breeding lines, limiting the transferability of results to practical wheat breeding. Here we assembled a population of advanced/elite breeding lines and well-adapted cultivars and evaluated over four environments for phenotypic analysis of spike and kernel traits. GWAS identified 17 significant multi-environment marker-trait associations (MTAs) for various traits, representing 12 putative quantitative trait loci (QTLs), with five QTLs affecting multiple traits. Four of these QTLs mapped on three chromosomes 1A, 5B, and 7A for spike length, number of spikelets per spike (NSPS), and kernel length are likely novel. Further, a highly significant QTL was detected on chromosome 7AS that has not been previously associated with NSPS and putative candidate genes were identified in this region. The allelic frequencies of important quantitative trait nucleotides (QTNs) were deduced in a larger set of 1,124 accessions which revealed the importance of identified MTAs in the US HWW breeding programs. The results from this study could be directly used by the breeders to select the lines with favorable alleles for making crosses, and reported markers will facilitate marker-assisted selection of stable QTLs for yield components in wheat breeding.

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.

Overexpression of the WAPO-A1 gene increases the number of spikelets per spike in bread wheat

Two homoeologous QTLs for number of spikelets per spike (SPS) were mapped on chromosomes 7AL and 7BL using two wheat MAGIC populations. Sets of lines contrasting for the QTL on 7AL were developed which allowed for the validation and fine mapping of the 7AL QTL and for the identification of a previously described candidate gene, WHEAT ORTHOLOG OF APO1 (WAPO1). Using transgenic overexpression in both a low and a high SPS line, we provide a functional validation for the role of this gene in determining SPS also in hexaploid wheat. We show that the expression levels of this gene positively correlate with SPS in multiple MAGIC founder lines under field conditions as well as in transgenic lines grown in the greenhouse. This work highlights the potential use of WAPO1 in hexaploid wheat for further yield increases. The impact of WAPO1 and SPS on yield depends on other genetic and environmental factors, hence, will require a finely balanced expression level to avoid the development of detrimental pleiotropic phenotypes.

Genetic analysis of morphological traits in spring wheat from the Northeast of China by a genome-wide association study

Identification of the gene for agronomic traits is important for the wheat marker-assisted selection (MAS) breeding. To identify the new and stable loci for agronomic traits, including flag leaf length (FLL), flag leaf width (FLW), uppermost internode length (UIL), and plant morphology (PM, including prostrate, semi-prostrate, and erect). A total of 251 spring wheat accessions collected from the Northeast of China were used to conduct genome-wide association study (GWAS) by 55K SNP arrays. A total of 30 loci for morphological traits were detected, and each explained 4.8–17.9% of the phenotypic variations. Of these, 13 loci have been reported by previous studies, and the other 17 are novel. We have identified seven genes involved in the signal transduction, cell-cycle progression, and plant development pathway as candidate genes. This study provides new insights into the genetic basis of morphological traits. The associated SNPs and accessions with more of favorable alleles identified in this study could be used to promote the wheat breeding progresses.

Genetic architecture and QTL selection response for Kernza perennial grain domestication traits

Analysis of multi-year breeding program data revealed that the genetic architecture of an intermediate wheatgrass population was highly polygenic for both domestication and agronomic traits, supporting the use of genomic selection for new crop domestication. Perennial grains have the potential to provide food for humans and decrease the negative impacts of annual agriculture. Intermediate wheatgrass (IWG, Thinopyrum intermedium, Kernza®) is a promising perennial grain candidate that The Land Institute has been breeding since 2003. We evaluated four consecutive breeding cycles of IWG from 2016 to 2020 with each cycle containing approximately 1100 unique genets. Using genotyping-by-sequencing markers, quantitative trait loci (QTL) were mapped for 34 different traits using genome-wide association analysis. Combining data across cycles and years, we found 93 marker-trait associations for 16 different traits, with each association explaining 0.8-5.2% of the observed phenotypic variance. Across the four cycles, only three QTL showed an F_ST differentiation > 0.15 with two corresponding to a decrease in floret shattering. Additionally, one marker associated with brittle rachis was 216 bp from an ortholog of the btr2 gene. Power analysis and quantitative genetic theory were used to estimate the effective number of QTL, which ranged from a minimum of 33 up to 558 QTL for individual traits. This study suggests that key agronomic and domestication traits are under polygenic control and that molecular methods like genomic selection are needed to accelerate domestication and improvement of this new crop.

Genetic Incorporation of the Favorable Alleles for Three Genes Associated With Spikelet Development in Wheat

The number of spikelets per spike is an important trait that directly affects grain yield in wheat. Three quantitative trait loci (QTLs) associated with spikelet nodes per spike (SNS) were mapped in a population of recombinant inbred lines generated from a cross between two advanced breeding lines of winter wheat based on the phenotypic variation evaluated over six locations/years. Two of the three QTLs are QSns.sxau-2A at the WHEATFRIZZY PANICLE (WFZP) loci and QSns.sxau-7A at the WHEAT ORTHOLOG OF APO1 (WAPO1) loci. The WFZP-A1b allele with a 14-bp deletion at QSns.sxau-2A was associated with increased spikelets per spike. WAPO-A1e, as a novel allele at WAPO1, were regulated at the transcript level that was associated with the SNS trait. The third SNS QTL, QSns.sxau-7D on chromosome 7D, was not associated with homoeologous WAPO-D1 or any other genes known to regulate SNS. The favorable alleles for each of WZFP-A1, WAPO-A1, and QSns.sxau-7D are identified and incorporated to increase up to 3.4 spikelets per spike in the RIL lines. Molecular markers for the alleles were developed. This study has advanced our understanding of the genetic basis of natural variation in spikelet development in wheat.

Identification and validation of QTLs for kernel number per spike and spike length in two founder genotypes of wheat

BACKGROUND

Kernel number per spike (KNS) and spike length (SL) are important spike-related traits in wheat variety improvement. Discovering genetic loci controlling these traits is necessary to elucidate the genetic basis of wheat yield traits and is very important for marker-assisted selection breeding.

RESULTS

In the present study, we used a recombinant inbred line population with 248 lines derived from the two founder genotypes of wheat, Bima4 and BainongAK58, to construct a high-density genetic map using wheat 55 K genotyping assay. The final genetic linkage map consists of 2356 bin markers (14,812 SNPs) representing all 21 wheat chromosomes, and the entire map spanned 4141.24 cM. A total of 7 and 18 QTLs were identified for KNS and SL, respectively, and they were distributed on 11 chromosomes. The allele effects of the flanking markers for 12 stable QTLs, including four QTLs for KNS and eight QTLs for SL, were estimated based on phenotyping data collected from 15 environments in a diverse wheat panel including 384 elite cultivars and breeding lines. The positive alleles at seven loci, namely, QKns.his-7D2-1, QKns.his-7D2-2, QSl.his-4A-1, QSl.his-5D1, QSl.his-4D2-2, QSl.his-5B and QSl.his-5A-2, significantly increased KNS or SL in the diverse panel, suggesting they are more universal in their effects and are valuable for gene pyramiding in breeding programs. The transmission of Bima4 allele indicated that the favorite alleles at five loci (QKns.his-7D2-1, QSl.his-5A-2, QSl.his-2D1-1, QSl.his-3A-2 and QSl.his-3B) showed a relatively high frequency or an upward trend following the continuity of generations, suggesting that they underwent rigorous selection during breeding. At two loci (QKns.his-7D2-1 and QSl.his-5A-2) that the positive effects of the Bima4 alleles have been validated in the diverse panel, two and one kompetitive allele-specific PCR (KASP) markers were further developed, respectively, and they are valuable for marker-assisted selection breeding.

CONCLUSION

Important chromosome regions controlling KNS and SL were identified in the founder parents. Our results are useful for knowing the molecular mechanisms of founder parents and future molecular breeding in wheat.

Genomic variants affecting homoeologous gene expression dosage contribute to agronomic trait variation in allopolyploid wheat

Allopolyploidy greatly expands the range of possible regulatory interactions among functionally redundant homoeologous genes. However, connection between the emerging regulatory complexity and expression and phenotypic diversity in polyploid crops remains elusive. Here, we use diverse wheat accessions to map expression quantitative trait loci (eQTL) and evaluate their effects on the population-scale variation in homoeolog expression dosage. The relative contribution of cis- and trans-eQTL to homoeolog expression variation is strongly affected by both selection and demographic events. Though trans-acting effects play major role in expression regulation, the expression dosage of homoeologs is largely influenced by cis-acting variants, which appear to be subjected to selection. The frequency and expression of homoeologous gene alleles showing strong expression dosage bias are predictive of variation in yield-related traits, and have likely been impacted by breeding for increased productivity. Our study highlights the importance of genomic variants affecting homoeolog expression dosage in shaping agronomic phenotypes and points at their potential utility for improving yield in polyploid crops.

GWAS identifies genetic loci underlying nitrogen responsiveness in the climate resilient C4 model Setaria italica (L.)

INTRODUCTION

N responsiveness is the capacity to perceive and induce morpho-physiological adaptation to external and internal Nitrogen (N). Crop productivity is propelled by N fertilizer and requires the breeding/selection of cultivars with intrinsically high N responsiveness. This trait has many advantages in being more meaningful in commercial/environmental context, facilitating in-season N management and not being inversely correlated with N availability over processes regulating NUE. Current lack of its understanding at the physio-genetic basis is an impediment to select for cultivars with a predictably high N response.

OBJECTIVES

To dissect physio-genetic basis of N responsiveness in 142 diverse population of foxtail millet, Setaria italica (L.) by employing contrasting N fertilizer nutrition regimes.

METHODS

We phenotyped S. italica accessions for major yield related traits under low (N10, N25) and optimal (N100) growth conditions and genotyped them to subsequently perform a genome-wide association study to identify genetic loci associated with nitrogen responsiveness trait. Groups of accessions showing contrasting trait performance and allelic forms of specific linked genetic loci (showing haplotypes) were further accessed for N dependent transcript abundances of their proximal genes.

RESULTS

Our study show that N dependent yield rise in S. italica is driven by grain number whose responsiveness to N availability is genetically underlined. We identify 22 unique SNP loci strongly associated with this trait out of which six exhibit haplotypes and consistent allelic variation between lines with contrasting N dependent grain number response and panicle architectures. Furthermore, differential transcript abundances of specific genes proximally linked to these SNPs in same lines is indicative of their N dependence in a genotype specific manner.

CONCLUSION

The study demonstrates the value/ potential of N responsiveness as a selection trait and identifies key genetic components underlying the trait in S. italica. This has major implications for improving crop N sustainability and food security.

WAPO-A1 is the causal gene of the 7AL QTL for spikelet number per spike in wheat

Improving our understanding of the genes regulating grain yield can contribute to the development of more productive wheat varieties. Previously, a highly significant QTL affecting spikelet number per spike (SNS), grain number per spike (GNS) and grain yield was detected on chromosome arm 7AL in multiple genome-wide association studies. Using a high-resolution genetic map, we established that the A-genome homeolog of WHEAT ORTHOLOG OF APO1 (WAPO-A1) was a leading candidate gene for this QTL. Using mutants and transgenic plants, we demonstrate in this study that WAPO-A1 is the causal gene underpinning this QTL. Loss-of-function mutants wapo-A1 and wapo-B1 showed reduced SNS in tetraploid wheat, and the effect was exacerbated in wapo1 combining both mutations. By contrast, spikes of transgenic wheat plants carrying extra copies of WAPO-A1 driven by its native promoter had higher SNS, a more compact spike apical region and a smaller terminal spikelet than the wild type. Taken together, these results indicate that WAPO1 affects SNS by regulating the timing of terminal spikelet formation. Both transgenic and wapo1 mutant plants showed a wide range of floral abnormalities, indicating additional roles of WAPO1 on wheat floral development. Previously, we found three widespread haplotypes in the QTL region (H1, H2 and H3), each associated with particular WAPO-A1 alleles. Results from this and our previous study show that the WAPO-A1 allele in the H1 haplotype (115-bp deletion in the promoter) is expressed at significantly lower levels in the developing spikes than the alleles in the H2 and H3 haplotypes, resulting in reduced SNS. Field experiments also showed that the H2 haplotype is associated with the strongest effects in increasing SNS and GNS (H2>H3>H1). The H2 haplotype is already present in most modern common wheat varieties but is rare in durum wheat, where it might be particularly useful to improve grain yield.

Deciphering the environmental impact on spike architectural traits for grain yield consolidation in bread wheat (T. aestivum L.)

The realization of grain yield in wheat is decided by source-sink balance under prevailing environmental conditions. Management conditions like changing the sowing time influence the source-sink capacity through modification in agronomic traits. Therefore, this experiment was conducted to decipher the influence of spike architectural traits (SATs) on grain yield and to open avenues for further grain yield enhancement. Comparatively early sowing over timely sowing gives the advantage of realizing higher grain yield with a positive relationship with SATs namely spike length, spikelets per spike, individual spike weight, individual grain weight, number of grains per spikelet, grain length, and grain width of upper and lower spike portion. Confirmatory factorial analysis revealed that spike length, spikelets per spike, individual spike weight, grains per spikelet were having a significant effect in deciding grain yield in early sown. The presence of a significant effect of genotype by environment interaction over grain yield and SATs allows the exploitation of available genotypic and environmental variability for further yield enhancement. GGE analysis on transformed and standardized grain yield-trait (GY-trait) combinations was used in the selection of genotypes having high GY-trait combinations for both sowing times. In early sowing, WG 11 was the best for high GY with high individual spike weight; grain length and grain width at lower and upper parts of the spike; and shorter days to 50% flowering. Genotypes exclusively having the high GY-trait combination along with low values of remaining GY-trait combinations were also selected with genotype focused GGE approach.

Ectopic expression of TaBG1 increases seed size and alters nutritional characteristics of the grain in wheat but does not lead to increased yields

BACKGROUND

Grain size is thought to be a major component of yield in many plant species. Here we set out to understand if knowledge from other cereals such as rice could translate to increased yield gains in wheat and lead to increased nitrogen use efficiency. Previous findings that the overexpression of OsBG1 in rice increased yields while increasing seed size suggest translating gains from rice to other cereals may help to increase yields.

RESULTS

The orthologous genes of OsBG1 were identified in wheat. One homoeologous wheat gene was cloned and overexpressed in wheat to understand its role in controlling seed size. Potential alteration in the nutritional profile of the grains were also analyzed in wheat overexpressing TaBG1. It was found that increased TaBG1-A expression could indeed lead to larger seed size but was linked to a reduction in seed number per plant leading to no significant overall increase in yield. Other important components of yield such as biomass or tillering did not change significantly with increased TaBG1-A expression. The nutritional profile of the grain was altered, with a significant decrease in the Zn levels in the grain associated with increased seed size, but Fe and Mn concentrations were unchanged. Protein content of the wheat grain also fell under moderate N fertilization levels but not under deficient or adequate levels of N.

CONCLUSIONS

TaBG1 does control seed size in wheat but increasing the seed size per se does not increase yield and may come at the cost of lower concentrations of essential elements as well as potentially lower protein content. Nevertheless, TaBG1 could be a useful target for further breeding efforts in combination with other genes for increased biomass.

Genomic selection of parents and crosses beyond the native gene pool of a breeding program

Genomic selection has become a valuable tool for selecting cultivar candidates in many plant breeding programs. Genomic selection of elite parents and crossing combinations with germplasm developed outside a breeding program has, however, hardly been explored until now. The aim of this study was to assess the potential of this method for commonly ranking and selecting elite germplasm developed within and beyond a given breeding program. A winter wheat (Triticum aestivum L.) population consisting of 611 in-house and 87 externally developed lines was used to compare training population compositions and statistical models for genomically predicting baking quality in this framework. Augmenting training populations with lines from other breeding programs had a larger influence on the prediction ability than adding in-house generated lines when aiming to commonly rank both germplasm sets. Exploiting preexisting information of secondary correlated traits resulted likewise in more accurate predictions both in empirical analyses and simulations. Genotyping germplasm developed beyond a given breeding program is moreover a convenient way to clarify its relationships with a breeder's own germplasm because pedigree information is oftentimes not available for this purpose. Genomic predictions can thus support a more informed diversity management, especially when integrating simply to phenotype correlated traits to partly circumvent resource reallocations for a costly phenotyping of germplasm from other programs.

Training set design in genomic prediction with multiple biparental families

Genomic selection is a powerful tool to reduce the cycle length and enhance the genetic gain of complex traits in plant breeding. However, questions remain about the optimum design and composition of the training set. In this study, we used 944 soybean [Glycine max (L.) Merr.] recombinant inbred lines from eight families derived through a partial-diallel mating design among five parental lines. The cross-validated prediction accuracies for the six traits seed yield, 1,000-seed weight, protein yield, plant height, protein content, and oil content were high, ranging from 0.79 to 0.87. We investigated among-family predictions, making use of the special mating design with different degrees of relatedness among families. Generally, the prediction accuracy decreased from full-sibs to half-sib families to unrelated families. However, half-sib and unrelated families also showed substantial variation in their prediction accuracy for a given family, which appeared to be caused at least in part by the shared segregation of quantitative trait loci in both the training and prediction sets. Combining several half-sib families in composite training sets generally led to an increase in the prediction accuracy compared with the best family alone. The prediction accuracy increased with the size of the training set, but for comparable prediction accuracy, substantially more half-sibs were required than full-sibs. Collectively, our results highlight the potential of genomic selection for soybean breeding and, in a broader context, illustrate the importance of the targeted design of the training set.

Increased Predictive Accuracy of Multi-Environment Genomic Prediction Model for Yield and Related Traits in Spring Wheat (Triticum aestivum L.)

Genomic selection (GS) has the potential to improve the selection gain for complex traits in crop breeding programs from resource-poor countries. The GS model performance in multi-environment (ME) trials was assessed for 141 advanced breeding lines under four field environments via cross-predictions. We compared prediction accuracy (PA) of two GS models with or without accounting for the environmental variation on four quantitative traits of significant importance, i.e., grain yield (GRYLD), thousand-grain weight, days to heading, and days to maturity, under North and Central Indian conditions. For each trait, we generated PA using the following two different ME cross-validation (CV) schemes representing actual breeding scenarios: (1) predicting untested lines in tested environments through the ME model (ME_CV1) and (2) predicting tested lines in untested environments through the ME model (ME_CV2). The ME predictions were compared with the baseline single-environment (SE) GS model (SE_CV1) representing a breeding scenario, where relationships and interactions are not leveraged across environments. Our results suggested that the ME models provide a clear advantage over SE models in terms of robust trait predictions. Both ME models provided 2-3 times higher prediction accuracies for all four traits across the four tested environments, highlighting the importance of accounting environmental variance in GS models. While the improvement in PA from SE to ME models was significant, the CV1 and CV2 schemes did not show any clear differences within ME, indicating the ME model was able to predict the untested environments and lines equally well. Overall, our results provide an important insight into the impact of environmental variation on GS in smaller breeding programs where these programs can potentially increase the rate of genetic gain by leveraging the ME wheat breeding trials.

Grain yield and correlated traits of bread wheat lines: Implications for yield improvement

Global wheat yields are suffering due to differences in regional climatic conditions and soil fertility. Plant breeders are continuously working to improve the yield per unit area of wheat crop through selecting superior lines as parents. The screening and field evaluation of available lines allow the selection of superior ones and subsequently improved varieties. Therefore, heritable distinctions among 33 bread wheat lines for yield and related attributes were assessed under field conditions. The experiment included thirty lines and three check varieties. Data relating to different plant characteristics was collected at maturity. Significant differences were recorded for yield and related traits of tested wheat lines and check varieties. Wheat lines V6, V12 and V20 proved better with reduced number of days to reach anthesis and other desirable traits compared to check varieties. Days to start heading had strong correlation with spike length and number of spikelets spike-1. Flag leaf area had positive relationship with peduncle length and yield related traits. The 1000-garin weight and grain yield were also correlated with each other. It is concluded that V6, V10 and V20 proved better for all studied traits than the rest of the lines. Therefore, these lines could be used in wheat breeding program as parents to improve yield.

Genetic dissection of grain architecture-related traits in a winter wheat population

BACKGROUND

The future productivity of wheat (T. aestivum L.) as the most grown crop worldwide is of utmost importance for global food security. Thousand kernel weight (TKW) in wheat is closely associated with grain architecture-related traits, e.g. kernel length (KL), kernel width (KW), kernel area (KA), kernel diameter ratio (KDR), and factor form density (FFD). Discovering the genetic architecture of natural variation in these traits, identifying QTL and candidate genes are the main aims of this study. Therefore, grain architecture-related traits in 261 worldwide winter accessions over three field-year experiments were evaluated.

RESULTS

Genome-wide association analysis using 90K SNP array in FarmCPU model revealed several interesting genomic regions including 17 significant SNPs passing false discovery rate threshold and strongly associated with the studied traits. Four of associated SNPs were physically located inside candidate genes within LD interval e.g. BobWhite_c5872_589 (602,710,399 bp) found to be inside TraesCS6A01G383800 (602,699,767-602,711,726 bp). Further analysis reveals the four novel candidate genes potentially involved in more than one grain architecture-related traits with a pleiotropic effects e.g. TraesCS6A01G383800 gene on 6A encoding oxidoreductase activity was associated with TKW and KA. The allelic variation at the associated SNPs showed significant differences betweeen the accessions carying the wild and mutated alleles e.g. accessions carying C allele of BobWhite_c5872_589, TraesCS6A01G383800 had significantly higher TKW than the accessions carying T allele. Interestingly, these genes were highly expressed in the grain-tissues, demonstrating their pivotal role in controlling the grain architecture.

CONCLUSIONS

These results are valuable for identifying regions associated with kernel weight and dimensions and potentially help breeders in improving kernel weight and architecture-related traits in order to increase wheat yield potential and end-use quality.

Crop adaptation to climate change as a consequence of long-term breeding

Major global crops in high-yielding, temperate cropping regions are facing increasing threats from the impact of climate change, particularly from drought and heat at critical developmental timepoints during the crop lifecycle. Research to address this concern is frequently focused on attempts to identify exotic genetic diversity showing pronounced stress tolerance or avoidance, to elucidate and introgress the responsible genetic factors or to discover underlying genes as a basis for targeted genetic modification. Although such approaches are occasionally successful in imparting a positive effect on performance in specific stress environments, for example through modulation of root depth, major-gene modifications of plant architecture or function tend to be highly context-dependent. In contrast, long-term genetic gain through conventional breeding has incrementally increased yields of modern crops through accumulation of beneficial, small-effect variants which also confer yield stability via stress adaptation. Here we reflect on retrospective breeding progress in major crops and the impact of long-term, conventional breeding on climate adaptation and yield stability under abiotic stress constraints. Looking forward, we outline how new approaches might complement conventional breeding to maintain and accelerate breeding progress, despite the challenges of climate change, as a prerequisite to sustainable future crop productivity.

Identification of eight QTL controlling multiple yield components in a German multi-parental wheat population, including Rht24, WAPO-A1, WAPO-B1 and genetic loci on chromosomes 5A and 6A

KEY MESSAGE

Quantitative trait locus (QTL) mapping of 15 yield component traits in a German multi-founder population identified eight QTL each controlling ≥2 phenotypes, including the genetic loci Rht24, WAPO-A1 and WAPO-B1. Grain yield in wheat (Triticum aestivum L.) is a polygenic trait representing the culmination of many developmental processes and their interactions with the environment. Toward maintaining genetic gains in yield potential, 'reductionist approaches' are commonly undertaken by which the genetic control of yield components, that collectively determine yield, are established. Here we use an eight-founder German multi-parental wheat population to investigate the genetic control and phenotypic trade-offs between 15 yield components. Increased grains per ear was significantly positively correlated with the number of fertile spikelets per ear and negatively correlated with the number of infertile spikelets. However, as increased grain number and fertile spikelet number per ear were significantly negatively correlated with thousand grain weight, sink strength limitations were evident. Genetic mapping identified 34 replicated quantitative trait loci (QTL) at two or more test environments, of which 24 resolved into eight loci each controlling two or more traits-termed here 'multi-trait QTL' (MT-QTL). These included MT-QTL associated with previously cloned genes controlling semi-dwarf plant stature, and with the genetic locus Reduced height 24 (Rht24) that further modulates plant height. Additionally, MT-QTL controlling spikelet number traits were located to chromosome 7A encompassing the gene WHEAT ORTHOLOG OF APO1 (WAPO-A1), and to its homoeologous location on chromosome 7B containing WAPO-B1. The genetic loci identified in this study, particularly those that potentially control multiple yield components, provide future opportunities for the targeted investigation of their underlying genes, gene networks and phenotypic trade-offs, in order to underpin further genetic gains in yield.

The da1 mutation in wheat increases grain size under ambient and elevated CO2 but not grain yield due to trade-off between grain size and grain number

Grain size is potentially yield determining in wheat, controlled by the ubiquitin pathway and negatively regulated by ubiquitin receptor DA1. We analyzed whether increased thousand grain weight in wheat da1 mutant is translated into higher grain yield and whether additional carbon provided by elevated (e)CO2 can be better used by the da1, displaying higher grain sink strength and size. Yield-related, biomass, grain quality traits, and grain dimensions were analyzed by two-factorial mixed-model analysis, regarding genotype and eCO2. da1 increased grain size but reduced spikes and grains per plant, grains per spike, and spikelets per spike, independent of eCO2 treatment, leaving total grain yield unchanged. eCO2 increased yield and grain number additively and independently of da1 but did not overcome the trade-off between grain size and number observed for da1. eCO2 but not da1 impaired grain quality, strongly decreasing concentrations of several macroelement and microelement. In conclusion, intrinsic stimulation of grain sink strength and grain size, achieved by da1, is not benefitting total yield unless trade-offs between grain size and numbers can be overcome. The results reveal interactions of yield components in da1-wheat under ambient and eCO2, thereby uncovering limitations enhancing wheat yield potential.

Wheat Developmental Traits as Affected by the Interaction between Eps-7D and Temperature under Contrasting Photoperiods with Insensitive Ppd-D1 Background

Earliness per se (Eps) genes are important to fine tune adaptation, and studying their probable pleiotropic effect on wheat yield traits is worthwhile. In addition, it has been shown that some Eps genes interact with temperature and therefore determining the likely Eps × temperature interaction is needed for each newly identified Eps gene. We studied two NILs differing in the newly identified Eps-7D (carrying insensitive Ppd-D1 in the background) under three temperature regimes (9, 15 and 18 °C) and two photoperiods (12 and 24 h). Eps-7D affected time to anthesis as expected and the Eps-7D-late allele extended both the period before and after terminal spikelet. The interaction effect of Eps-7D × temperature was significant but not cross-over: the magnitude and level of significance of the difference between NILs with the late or early allele was affected by the growing temperature (i.e., difference was least at 18 °C and largest at 9 °C), and the differences caused due to temperature sensitivity were influenced by photoperiod. The rate of leaf initiation was faster in NIL with Eps-7D-early than with the late allele which compensated for the shorter duration of leaf initiation resulting in similar final leaf number between two NILs. Eps-7D-late consistently increased spike fertility through improving floret primordia survival as a consequence of extending the late reproductive phase.

Determination of wheat spike and spikelet architecture and grain traits using X-ray Computed Tomography imaging

BACKGROUND

Wheat spike architecture is a key determinant of multiple grain yield components and detailed examination of spike morphometric traits is beneficial to explain wheat grain yield and the effects of differing agronomy and genetics. However, quantification of spike morphometric traits has been very limited because it relies on time-consuming manual measurements.

RESULTS

In this study, using X-ray Computed Tomography imaging, we proposed a method to efficiently detect the 3D architecture of wheat spikes and component spikelets by clustering grains based on their Euclidean distance and relative positions. Morphometric characteristics of wheat spikelets and grains, e.g., number, size and spatial distribution along the spike can be determined. Two commercial wheat cultivars, one old, Maris Widgeon, and one modern, Siskin, were studied as examples. The average grain volume of Maris Widgeon and Siskin did not differ, but Siskin had more grains per spike and therefore greater total grain volume per spike. The spike length and spikelet number were not statistically different between the two cultivars. However, Siskin had a higher spikelet density (number of spikelets per unit spike length), with more grains and greater grain volume per spikelet than Maris Widgeon. Spatial distribution analysis revealed the number of grains, the average grain volume and the total grain volume of individual spikelets varied along the spike. Siskin had more grains and greater grain volumes per spikelet from spikelet 6, but not spikelet 1-5, compared with Maris Widgeon. The distribution of average grain volume along the spike was similar for the two wheat cultivars.

CONCLUSION

The proposed method can efficiently extract spike, spikelet and grain morphometric traits of different wheat cultivars, which can contribute to a more detailed understanding of the sink of wheat grain yield.