Genome-Wide Mapping of Spike-Related and Agronomic Traits in a Common Wheat Population Derived from a Supernumerary Spikelet Parent and an Elite Parent

Identification and validation of two major QTL for grain number per spike on chromosomes 2B and 2D in bread wheat (Triticum aestivum L.)

KEY MESSAGE

Major QTL for grain number per spike were identified on chromosomes 2B and 2D. Haplotypes and candidate genes of QGns.cib-2B.1 were analyzed. Grain number per spike (GNS) is one of the main components of wheat yield. Genetic dissection of their regulatory factors is essential to improve the yield potential. In present study, a recombinant inbred line population comprising 180 lines developed from the cross between a high GNS line W7268 and a cultivar Chuanyu12 was employed to identify quantitative trait loci (QTL) associated with GNS across six environments. Two major QTL, QGns.cib-2B.1 and QGns.cib-2D.1, were detected in at least four environments with the phenotypic variations of 12.99-27.07% and 8.50-13.79%, respectively. And significant interactions were observed between the two major QTL. In addition, QGns.cib-2B.1 is a QTL cluster for GNS, grain number per spikelet and fertile tiller number, and they were validated in different genetic backgrounds using Kompetitive Allele Specific PCR (KASP) markers. QGns.cib-2B.1 showed pleotropic effects on other yield-related traits including plant height, spike length, and spikelet number per spike, but did not significantly affect thousand grain weight which suggested that it might be potentially applicable in breeding program. Comparison analysis suggested that QGns.cib-2B.1 might be a novel QTL. Furthermore, haplotype analysis of QGns.cib-2B.1 indicated that it is a hot spot of artificial selection during wheat improvement. Based on the expression patterns, gene annotation, orthologs analysis and sequence variations, the candidate genes of QGns.cib-2B.1 were predicted. Collectively, the major QTL and KASP markers reported here provided a wealth of information for the genetic basis of GNS and grain yield improvement.

Genomic insights into glume pubescence in durum wheat: GWAS and haplotype analysis implicates TdELD1-1A as a candidate gene

Glume pubescence is an important morphological trait for the characterization of wheat cultivars. It shows tolerance to biotic and abiotic stresses to some extent. Hg1 (formerly named Hg) locus on chromosome 1AS controls glume pubescence in wheat. Its genetic analysis, fine-mapping and candidate gene analysis have been widely studied recently, however, the cloning of Hg1 has not yet been reported. Here, we conducted a GWAS between a dense panel of 171,103 SNPs and glume pubescence (Gp) in a durum wheat population of 145 lines, and further analyzed the candidate genes of Hg1 combined with the gene expression, functional annotation, and haplotype analysis. As a results, TRITD0Uv1G104670 (TdELD1-1A), encoding glycosyltransferase-like ELD1/KOBITO 1, was detected as the most promising candidate gene of Hg1 for glume pubescence in durum wheat. Our findings not only contribute to a deeper understanding of its cloning and functional validation but also underscore the significance of accurate genome sequences and annotations. Additionally, our study highlights the relevance of unanchored sequences in chrUn and the application of bioinformatics analysis for gene discovery in durum wheat.

Image-based classification of wheat spikes by glume pubescence using convolutional neural networks

Introduction

Pubescence is an important phenotypic trait observed in both vegetative and generative plant organs. Pubescent plants demonstrate increased resistance to various environmental stresses such as drought, low temperatures, and pests. It serves as a significant morphological marker and aids in selecting stress-resistant cultivars, particularly in wheat. In wheat, pubescence is visible on leaves, leaf sheath, glumes and nodes. Regarding glumes, the presence of pubescence plays a pivotal role in its classification. It supplements other spike characteristics, aiding in distinguishing between different varieties within the wheat species. The determination of pubescence typically involves visual analysis by an expert. However, methods without the use of binocular loupe tend to be subjective, while employing additional equipment is labor-intensive. This paper proposes an integrated approach to determine glume pubescence presence in spike images captured under laboratory conditions using a digital camera and convolutional neural networks.

Methods

Initially, image segmentation is conducted to extract the contour of the spike body, followed by cropping of the spike images to an equal size. These images are then classified based on glume pubescence (pubescent/glabrous) using various convolutional neural network architectures (Resnet-18, EfficientNet-B0, and EfficientNet-B1). The networks were trained and tested on a dataset comprising 9,719 spike images.

Results

For segmentation, the U-Net model with EfficientNet-B1 encoder was chosen, achieving the segmentation accuracy IoU = 0.947 for the spike body and 0.777 for awns. The classification model for glume pubescence with the highest performance utilized the EfficientNet-B1 architecture. On the test sample, the model exhibited prediction accuracy parameters of F1 = 0.85 and AUC = 0.96, while on the holdout sample it showed F1 = 0.84 and AUC = 0.89. Additionally, the study investigated the relationship between image scale, artificial distortions, and model prediction performance, revealing that higher magnification and smaller distortions yielded a more accurate prediction of glume pubescence.

Unraveling the genetic basis of grain number-related traits in a wheat-Agropyron cristatum introgressed line through high-resolution linkage mapping

BACKGROUND

Grain number per spike (GNS) is a pivotal determinant of grain yield in wheat. Pubing 3228 (PB3228), a wheat-Agropyron cristatum germplasm, exhibits a notably higher GNS.

RESULTS

In this study, we developed a recombinant inbred line (RIL) population derived from PB3228/Gao8901 (PG-RIL) and constructed a high-density genetic map comprising 101,136 loci, spanning 4357.3 cM using the Wheat 660 K SNP array. The genetic map demonstrated high collinearity with the wheat assembly IWGSC RefSeq v1.0. Traits related to grain number and spikelet number per spike were evaluated in seven environments for quantitative trait locus (QTL) analysis. Five environmentally stable QTLs were detected in at least three environments. Among these, two major QTLs, QGns-4A.2 and QGns-1A.1, associated with GNS, exhibited positive alleles contributed by PB3228. Further, the conditional QTL analysis revealed a predominant contribution of PB3228 to the GNS QTLs, with both grain number per spikelet (GNSL) and spikelet number per spike (SNS) contributing to the overall GNS trait. Four kompetitive allele-specific PCR (KASP) markers that linked to QGns-4A.2 and QGns-1A.1 were developed and found to be effective in verifying the QTL effect within a diversity panel. Compared to previous studies, QGns-4A.2 exhibited stability across different trials, while QGns-1A.1 represents a novel QTL. The results from unconditional and conditional QTL analyses are valuable for dissecting the genetic contribution of the component traits to GNS at the individual QTL level and for understanding the genetic basis of the superior grain number character in PB3228. The KASP markers can be utilized in marker-assisted selection for enhancing GNS.

CONCLUSIONS

Five environmentally stable QTLs related to grain number and spikelet number per spike were identified. PB3228 contributed to the majority of the QTLs associated with GNS.

QTL mapping of agronomic traits in wheat using the UK Avalon × Cadenza reference mapping population grown in Kazakhstan

Background

The success of wheat production is largely dependent on local breeding projects that focus on the development of high-yielding cultivars with the use of novel molecular tools. One strategy for improving wheat productivity involves the deployment of diverse germplasms with a high potential yield. An important factor for achieving success involves the dissection of quantitative trait loci (QTLs) for complex agronomic traits, such as grain yield components, in targeted environments for wheat growth.

Methods

In this study, we tested the United Kingdom (UK) spring set of the doubled haploid (DH) reference population derived from the cross between two British cultivars, Avalon (winter wheat) and Cadenza (spring wheat), in the Northern, Central, and Southern regions (Karabalyk, Karaganda, Kyzylorda) of Kazakhstan over three years (2013–2015). The DH population has previously been genotyped by UK scientists using 3647 polymorphic DNA markers. The list of tested traits includes the heading time, seed maturation time, plant height, spike length, productive tillering, number of kernels per spike, number of kernels per meter, thousand kernel weight, and yield per square meter. Windows QTL Cartographer was applied for QTL mapping using the composite interval mapping method.

Results

In total, 83 out of 232 QTLs were identified as stable QTLs from at least two environments. A literature survey suggests that 40 QTLs had previously been reported elsewhere, indicating that this study identified 43 QTLs that are presumably novel marker-trait associations (MTA) for these environments. Hence, the phenotyping of the DH population in new environments led to the discovery of novel MTAs. The identified SNP markers associated with agronomic traits in the DH population could be successfully used in local Kazakh breeding projects for the improvement of wheat productivity.

Characterization of the Durum Wheat-Aegilops tauschii 4D(4B) Disomic Substitution Line YL-443 With Superior Characteristics of High Yielding and Stripe Rust Resistance

Durum wheat is one of the important food and cash crops. The main goals in current breeding programs are improving its low yield potential, kernel characteristics, and lack of resistance or tolerance to some biotic and abiotic stresses. In this study, a nascent synthesized hexaploid wheat Lanmai/AT23 is used as the female parent in crosses with its AB genome donor Lanmai. A tetraploid line YL-443 with supernumerary spikelets and high resistance to stripe rust was selected out from the pentaploid F₇ progeny. Somatic analysis using multicolor fluorescence in situ hybridization (mc-FISH) revealed that this line is a disomic substitution line with the 4B chromosome pair of Lanmai replaced by the 4D chromosome pair of Aegilops tauschii AT23. Comparing with Lanmai, YL-443 shows an increase in the number of spikelets and florets per spike by 36.3 and 75.9%, respectively. The stripe rust resistance gene Yr28 carried on the 4D chromosome was fully expressed in the tetraploid background. The present 4D(4B) disomic substitution line YL-443 was distinguished from the previously reported 4D(4B) lines with the 4D chromosomes from Chinese Spring (CS). Our study demonstrated that YL-443 can be used as elite germplasm for durum wheat breeding targeting high yield potential and stripe rust resistance. The Yr28-specific PCR marker and the 4D chromosome-specific KASP markers together with its unique features of pubescent leaf sheath and auricles can be utilized for assisting selection in breeding.

Genome wide genetic dissection of wheat quality and yield related traits and their relationship with grain shape and size traits in an elite × non-adapted bread wheat cross

The genetic gain in yield and quality are two major targets of wheat breeding programs around the world. In this study, a high density genetic map consisting of 10,172 SNP markers identified a total of 43 genomic regions associated with three quality traits, three yield traits and two agronomic traits in hard red spring wheat (HRSW). When compared with six grain shape and size traits, the quality traits showed mostly independent genetic control (~18% common loci), while the yield traits showed moderate association (~53% common loci). Association of genomic regions for grain area (GA) and thousand-grain weight (TGW), with yield suggests that targeting an increase in GA may help enhancing wheat yield through an increase in TGW. Flour extraction (FE), although has a weak positive phenotypic association with grain shape and size, they do not share any common genetic loci. A major contributor to plant height was the Rht8 locus and the reduced height allele was associated with significant increase in grains per spike (GPS) and FE, and decrease in number of spikes per square meter and test weight. Stable loci were identified for almost all the traits. However, we could not find any QTL in the region of major known genes like GPC-B1, Ha, Rht-1, and Ppd-1. Epistasis also played an important role in the genetics of majority of the traits. In addition to enhancing our knowledge about the association of wheat quality and yield with grain shape and size, this study provides novel loci, genetic information and pre-breeding material (combining positive alleles from both parents) to enhance the cultivated gene pool in wheat germplasm. These resources are valuable in facilitating molecular breeding for wheat quality and yield improvement.

Genetic insights into morphometric inflorescence traits of wheat

KEY MESSAGE

Modifying morphometric inflorescence traits is important for increasing grain yield in wheat. Mapping revealed nine QTL, including new QTL and a new allele for the q locus, controlling wheat spike morphometric traits. To identify loci controlling spike morphometric traits, namely spike length (SL), internode length (IL), node number per spike (NPS), and node density (ND), we studied 146 Recombinant Inbred Lines of tetraploid wheat (Triticum turgidum L.) derived from standard spike and spike-branching mutant parents. Phenotypic analyses of spike morphometric traits showed low genetic coefficients of variation, resulting in high heritabilities. The phenotypic correlation between NPS with growing degree days (GDD) suggested the importance of GDD in the determination of node number in wheat. The major effect QTL for GDD or heading date was mapped to chromosome 7BS carrying the flowering time gene, Vrn3-B1. Mapping also identified nine QTL controlling spike morphometric traits. Most of these loci controlled more than a single trait, suggesting a close genetic interrelationship among spike morphometric traits. For example, this study identified a new QTL, QND.ipk-4AL, controlling ND (up to 17.6% of the phenotypic variance), IL (up to 11% of the phenotypic variance), and SL (up to 20.8% of the phenotypic variance). Similarly, the major effect QTL for IL was mapped to the q locus. Sequencing of the Q/q gene further revealed a new q allele, qdel-5A, in spike-branching accessions possessing a six base pair deletion close to the miR172 target site. The identification of qdel-5A suggested that the spike-branching tetraploid wheats are double mutants for the spikelet meristem (SM) identity gene, i.e., branched headt (TtBHt), and the q gene, which is believed to be involved in the SM indeterminacy complex in wheat.

Uncovering Genomic Regions Associated With 36 Agro-Morphological Traits in Indian Spring Wheat Using GWAS

Wheat genetic improvement by integration of advanced genomic technologies is one way of improving productivity. To facilitate the breeding of economically important traits in wheat, SNP loci and underlying candidate genes associated with the 36 agro-morphological traits were studied in a diverse panel of 404 genotypes. By using Breeders' 35K Axiom array in a comprehensive genome-wide association study covering 4364.79 cM of the wheat genome and applying a compressed mixed linear model, a total of 146 SNPs (-log10 P ≥ 4) were found associated with 23 traits out of 36 traits studied explaining 3.7-47.0% of phenotypic variance. To reveal this a subset of 260 genotypes was characterized phenotypically for six quantitative traits [days to heading (DTH), days to maturity (DTM), plant height (PH), spike length (SL), awn length (Awn_L), and leaf length (Leaf_L)] under five environments. Gene annotations mined ∼38 putative candidate genes which were confirmed using tissue and stage specific gene expression data from RNA Seq. We observed strong co-localized loci for four traits (glume pubescence, SL, PH, and awn color) on chromosome 1B (24.64 cM) annotated five putative candidate genes. This study led to the discovery of hitherto unreported loci for some less explored traits (such as leaf sheath wax, awn attitude, and glume pubescence) besides the refined chromosomal regions of known loci associated with the traits. This study provides valuable information of the genetic loci and their potential genes underlying the traits such as awn characters which are being considered as important contributors toward yield enhancement.

Genome-Wide Association Study Reveals Novel Genomic Regions for Grain Yield and Yield-Related Traits in Drought-Stressed Synthetic Hexaploid Wheat

Synthetic hexaploid wheat (SHW; 2n = 6x = 42, AABBDD, Triticum aestivum L.) is produced from an interspecific cross between durum wheat (2n = 4x = 28, AABB, T. turgidum L.) and goat grass (2n = 2x = 14, DD, Aegilops tauschii Coss.) and is reported to have significant novel alleles-controlling biotic and abiotic stresses resistance. A genome-wide association study (GWAS) was conducted to unravel these loci [marker⁻trait associations (MTAs)] using 35,648 genotyping-by-sequencing-derived single nucleotide polymorphisms in 123 SHWs. We identified 90 novel MTAs (45, 11, and 34 on the A, B, and D genomes, respectively) and haplotype blocks associated with grain yield and yield-related traits including root traits under drought stress. The phenotypic variance explained by the MTAs ranged from 1.1% to 32.3%. Most of the MTAs (120 out of 194) identified were found in genes, and of these 45 MTAs were in genes annotated as having a potential role in drought stress. This result provides further evidence for the reliability of MTAs identified. The large number of MTAs (53) identified especially on the D-genome demonstrate the potential of SHWs for elucidating the genetic architecture of complex traits and provide an opportunity for further improvement of wheat under rapidly changing climatic conditions.

Environmental and Genetic Variation for Yield-Related Traits of Durum Wheat as Affected by Development

Phenology has a profound effect on adaptation and productivity of crops. The impact of phenology on tillering and fertility traits of durum wheat (Triticum turgidum L. subsp. durum Desf.) was evaluated with the aim of specifying which group of flowering genes (Vrn, Ppd, or eps) was involved in their control. A recombinant inbred line population was grown under four contrasting conditions of vernalization and daylength. Phenotyping was carried out according to robust phenological models dissecting both phenological and yield related traits. Whole-genome profiling was performed using the DArT-Seq technology. The genetic variability for tillering was mainly related to the genetic variability for vernalization sensitivity, as shown by the many quantitative trait loci (QTLs) identified in non-vernalized plants associated to both tillering and phenological traits. No effects of photoperiod sensitivity on spikelet number were detected in short-day-grown plants, apparently because of limited genetic variability in photoperiod sensitivity of the population. Earliness per se was involved in control of spikelet number via final leaf number, with these traits genetically correlated and sharing some QTLs. Chaff weight and number of kernels per g chaff were negatively associated and related to anthesis date under most conditions. QTL mapping uncovered novel loci involved in phenological control of tillering and fertility traits, and confirmed the presence of several well-established loci. Phenotyping of both phenology and kernel number according to a robust physiological model amplified the possibility of identifying genetic factors underlying their variations. Also, isolating known flowering gene cues by manipulation of environmental conditions provided the opportunity for each group of genes to be expressed without confounding effects of the others. This information helps to predict the consequences of either genetic manipulation of flowering genes and changes in environmental conditions on the potential yield of durum wheat.

Genome-wide association study and genetic diversity analysis on nitrogen use efficiency in a Central European winter wheat (Triticum aestivum L.) collection

To satisfy future demands, the increase of wheat (Triticum aestivum L.) yield is inevitable. Simultaneously, maintaining high crop productivity and efficient use of nutrients, especially nitrogen use efficiency (NUE), are essential for sustainable agriculture. NUE and its components are inherently complex and highly influenced by environmental factors, nitrogen management practices and genotypic variation. Therefore, a better understanding of their genetic basis and regulation is fundamental. To investigate NUE-related traits and their genetic and environmental regulation, field trials were evaluated in a Central European wheat collection of 93 cultivars at two nitrogen input levels across three seasons. This elite germplasm collection was genotyped on DArTseq® genotypic platform to identify loci affecting N-related complex agronomic traits. To conduct robust genome-wide association mapping, the genetic diversity, population structure and linkage disequilibrium were examined. Population structure was investigated by various methods and two subpopulations were identified. Their separation is based on the breeding history of the cultivars, while analysis of linkage disequilibrium suggested that selective pressures had acted on genomic regions bearing loci with remarkable agronomic importance. Besides NUE, genetic basis for variation in agronomic traits indirectly affecting NUE and its components, moreover genetic loci underlying response to nitrogen fertilisation were also determined. Altogether, 183 marker-trait associations (MTA) were identified spreading over almost the entire genome. We found that most of the MTAs were environmental-dependent. The present study identified several associated markers in those genomic regions where previous reports had found genes or quantitative trait loci influencing the same traits, while most of the MTAs revealed new genomic regions. Our data provides an overview of the allele composition of bread wheat varieties anchored to DArTseq® markers, which will facilitate the understanding of the genetic basis of NUE and agronomically important traits.

Mapping QTLs controlling kernel dimensions in a wheat inter-varietal RIL mapping population

Seven kernel dimension QTLs were identified in wheat, and kernel thickness was found to be the most important dimension for grain weight improvement. Kernel morphology and weight of wheat (Triticum aestivum L.) affect both yield and quality; however, the genetic basis of these traits and their interactions has not been fully understood. In this study, to investigate the genetic factors affecting kernel morphology and the association of kernel morphology traits with kernel weight, kernel length (KL), width (KW) and thickness (KT) were evaluated, together with hundred-grain weight (HGW), in a recombinant inbred line population derived from Nanda2419 × Wangshuibai, with data from five trials (two different locations over 3 years). The results showed that HGW was more closely correlated with KT and KW than with KL. A whole genome scan revealed four QTLs for KL, one for KW and two for KT, distributed on five different chromosomes. Of them, QKl.nau-2D for KL, and QKt.nau-4B and QKt.nau-5A for KT were newly identified major QTLs for the respective traits, explaining up to 32.6 and 41.5% of the phenotypic variations, respectively. Increase of KW and KT and reduction of KL/KT and KW/KT ratios always resulted in significant higher grain weight. Lines combining the Nanda 2419 alleles of the 4B and 5A intervals had wider, thicker, rounder kernels and a 14% higher grain weight in the genotype-based analysis. A strong, negative linear relationship of the KW/KT ratio with grain weight was observed. It thus appears that kernel thickness is the most important kernel dimension factor in wheat improvement for higher yield. Mapping and marker identification of the kernel dimension-related QTLs definitely help realize the breeding goals.

Radiation hybrid maps of the D-genome of Aegilops tauschii and their application in sequence assembly of large and complex plant genomes

BACKGROUND

The large and complex genome of bread wheat (Triticum aestivum L., ~17 Gb) requires high resolution genome maps with saturated marker scaffolds to anchor and orient BAC contigs/ sequence scaffolds for whole genome assembly. Radiation hybrid (RH) mapping has proven to be an excellent tool for the development of such maps for it offers much higher and more uniform marker resolution across the length of the chromosome compared to genetic mapping and does not require marker polymorphism per se, as it is based on presence (retention) vs. absence (deletion) marker assay.

METHODS

In this study, a 178 line RH panel was genotyped with SSRs and DArT markers to develop the first high resolution RH maps of the entire D-genome of Ae. tauschii accession AL8/78. To confirm map order accuracy, the AL8/78-RH maps were compared with:1) a DArT consensus genetic map constructed using more than 100 bi-parental populations, 2) a RH map of the D-genome of reference hexaploid wheat 'Chinese Spring', and 3) two SNP-based genetic maps, one with anchored D-genome BAC contigs and another with anchored D-genome sequence scaffolds. Using marker sequences, the RH maps were also anchored with a BAC contig based physical map and draft sequence of the D-genome of Ae. tauschii.

RESULTS

A total of 609 markers were mapped to 503 unique positions on the seven D-genome chromosomes, with a total map length of 14,706.7 cR. The average distance between any two marker loci was 29.2 cR which corresponds to 2.1 cM or 9.8 Mb. The average mapping resolution across the D-genome was estimated to be 0.34 Mb (Mb/cR) or 0.07 cM (cM/cR). The RH maps showed almost perfect agreement with several published maps with regard to chromosome assignments of markers. The mean rank correlations between the position of markers on AL8/78 maps and the four published maps, ranged from 0.75 to 0.92, suggesting a good agreement in marker order. With 609 mapped markers, a total of 2481 deletions for the whole D-genome were detected with an average deletion size of 42.0 Mb. A total of 520 markers were anchored to 216 Ae. tauschii sequence scaffolds, 116 of which were not anchored earlier to the D-genome.

CONCLUSION

This study reports the development of first high resolution RH maps for the D-genome of Ae. tauschii accession AL8/78, which were then used for the anchoring of unassigned sequence scaffolds. This study demonstrates how RH mapping, which offered high and uniform resolution across the length of the chromosome, can facilitate the complete sequence assembly of the large and complex plant genomes.