Resequencing of 145 Landmark Cultivars Reveals Asymmetric Sub-genome Selection and Strong Founder Genotype Effects on Wheat Breeding in China

A system genetics analysis uncovers the regulatory variants controlling drought response in wheat

Plants activate a variable response to drought stress by modulating transcription of key genes. However, our knowledge of genetic variations governing gene expression in response to drought stress remains limited in natural germplasm. Here, we performed a comprehensive analysis of the transcriptional variability of 200 wheat accessions in response to drought stress by using a systems genetics approach integrating pan-transcriptome, co-expression networks, transcriptome-wide association study (TWAS), and expression quantitative trait loci (eQTLs) mapping. We identified 1621 genes and eight co-expression modules significantly correlated with wheat drought tolerance. We also defined 620 664 and 654 798 independent eQTLs associated with the expression of 17 429 and 18 080 eGenes under normal and drought stress conditions. Focusing on dynamic regulatory variants, we further identified 572 eQTL hotspots and constructed transcription factors governed drought-responsive network by the XGBoost model. Subsequently, by combining with genome-wide association study (GWAS), we uncovered a 369-bp insertion variant in the TaKCS3 promoter containing multiple cis-regulatory elements recognized by eQTL hotspot-associated transcription factors that enhance its transcription. Further functional analysis indicated that elevating TaKCS3 expression affects cuticular wax composition to reduce water loss during drought stress, and thereby increase drought tolerance. This study sheds light on the genome-wide genetic variants that influence dynamic transcriptional changes during drought stress and provides a valuable resource for the mining of drought-tolerant genes in the future.

Genomic structural variation in an alpha/beta hydrolase triggers hybrid necrosis in wheat

Hybrid necrosis, a century-old mystery in wheat, is caused by complementary genes Ne1 and Ne2. Ne2, encoding a nucleotide-binding leucine-rich repeat (NLR) immune receptor, has been cloned, yet Ne1 remains elusive. Here, we report that Ne1, which encodes an alpha/beta hydrolase (ABH) protein generated by structural variation, triggers hybrid necrosis with Ne2 by activating autoimmune responses. We further verify that not only allelic variation but also copy number variation (CNV) of Ne1 are pivotal for hybrid necrosis diversity in wheat. Ne1 likely originates from wild emmer wheat, potentially through duplication and ectopic recombination events. Unlike Ne2, which is frequently selected for rust resistance in wheat breeding, the lower prevalence of Ne1 in modern wheat cultivars is attributed to its association with hybrid necrosis. Altogether, these findings illuminate the co-evolution of the NLR/ABH gene pair in plant development and innate immunity, offering potential benefits for wheat breeding.

TaFT-D1 positively regulates grain weight by acting as a coactivator of TaFDL2 in wheat

FLOWERING LOCUS T (FT), a multifunctional regulator in crops, modulates multiple key agronomic traits such as flowering time or heading date and plant height; however, its role in grain development regulation is unclear. Herein, through genome-wide association studies (GWAS), we identified TaFT-D1, which encodes a phosphatidylethanolamine-binding protein (PEBP), as a candidate gene for grain weight in wheat. A one-bp insertion/deletion (InDel) (G/-) in the third exon of TaFT-D1, resulting in different protein lengths, was significantly associated with grain weight. TaFT-D1 knockout via the CRISPR-Cas9 system reduced grain size and weight, and TaFT-D1 increased grain size by promoting cell proliferation and starch synthesis. Transcriptome analysis revealed a significant decrease in the expression of cell cycle- and starch synthesis-related genes, including TaNAC019-3A, TaSWEET15-like-7B, TaCYCD4;1 and TaCYCD3;2, in the taft-d1 knockout line. TaFT-D1 interacted with the bZIP transcription factor TaFDL2, and the tafdl2 mutant presented relatively small grains, suggesting that TaFDL2 is a positive regulator of grain size. Moreover, TaFDL2 bound to the promoters of downstream cell cycle- and starch synthesis-related genes, activating their expression, whereas TaFT-D1 increased this activation via TaFDL2. Interaction assays demonstrated that TaFT-D1, Ta14-3-3A and TaFDL2 formed a regulatory complex. Furthermore, the TaFT-D1(G) allele was significantly correlated with greater thousand-grain weight and earlier heading. This favourable allele has undergone strong positive selection during wheat breeding in China. Our findings provide novel insights into how TaFT-D1 regulates grain weight and highlight its potential application for yield improvement in wheat.

Genomic analysis of Zhou8425B, a key founder parent, reveals its genetic contributions to elite agronomic traits in wheat breeding

High-quality genome information is essential for efficiently deciphering and improving crop traits. Here, we report a highly contiguous and accurate hexaploid genome assembly for the key wheat breeding parent Zhou8425B, an elite 1BL/1RS translocation line with durable adult plant resistance (APR) against yellow rust (YR) disease. By integrating HiFi and Hi-C sequencing reads, we have generated a 14.75-Gb genome assembly for Zhou8425B with a contig N50 of 70.94 and a scaffold N50 of 735.11 Mb. Comparisons with previously sequenced common wheat cultivars shed light on structural changes in the 1RS chromosome arm, which has been extensively used in wheat improvement. Interestingly, Zhou8425B 1RS carries more genes encoding AP2/ERF-ERF or B3 transcription factors than its counterparts in four previously sequenced wheat and rye genotypes. The Zhou8425B genome assembly aided in the fine mapping of a new APR locus (YrZH3BS) that confers resistance to YR disease and promotes grain yield under field conditions. Notably, pyramiding YrZH3BS with two previously characterized APR loci (YrZH22 and YrZH84) can further reduce YR severity and enhance grain yield, with the triple combination (YrZH3B + YrZH22 + YrZH84) having the greatest effect. Finally, the founder genotype effects of Zhou8425B were explored using publicly available genome resequencing data, which reveals the presence of important Zhou8425B genomic blocks in its derivative cultivars. Our data demonstrate the value of the Zhou8425B genome assembly for further study of the structural and functional characteristics of 1RS, the genetic basis of durable YR resistance, and founder genotype effects in wheat breeding. Our resources will facilitate the development of elite wheat cultivars through genomics-assisted breeding.

Wheat2035: Integrating pan-omics and advanced biotechnology for future wheat design

Wheat (Triticum aestivum) production is vital for global food security, providing energy and protein to millions of people worldwide. Recent advancements in wheat research have led to significant increases in production, fueled by technological and scientific innovation. Here, we summarize the major advancements in wheat research, particularly the integration of biotechnologies and a deeper understanding of wheat biology. The shift from multi-omics to pan-omics approaches in wheat research has greatly enhanced our understanding of the complex genome, genomic variations, and regulatory networks to decode complex traits. We also outline key scientific questions, potential research directions, and technological strategies for improving wheat over the next decade. Since global wheat production is expected to increase by 60% in 2050, continued innovation and collaboration are crucial. Integrating biotechnologies and a deeper understanding of wheat biology will be essential for addressing future challenges in wheat production, ensuring sustainable practices and improved productivity.

The location and genome origin of alien chromatin in wheat founder parent Xiaoyan 6

KEY MESSAGE

The location of alien chromatin in Xiaoyan 6 was identified using mc-GISH analysis, genetic mapping and whole genome re-sequencing, and its possible origin was discussed. As a founder parent, Xiaoyan 6 has played an important role in distant hybridization breeding in China. Although it came from the cross between common wheat and Thinopyrum ponticum (Podp.) Barkworth and D.R. Dewey, the location of its alien chromatin has not been determined using traditional genomic in situ hybridization (GISH). In the present study, chromosome variation in Xiaoyan 6 was discovered by multicolor GISH analysis. Four alien-specific markers were developed by specific-locus amplified fragment sequencing technique. Their amplified sequences were analyzed by basic local alignment search tool with the reference genome sequences of common wheat Chinese Spring (CS) and Th. elongatum, and the whole genome re-sequencing reads of Th. ponticum and CS. Furthermore, the four markers were mapped on three different chromosomes in two RIL populations. By dissecting the mapped reads depth of the whole genome re-sequencing of Xiaoyan 6, we found that the depth of nine chromosome regions was obviously lower than the average. Among these, three regions on 1A, 3A and 7B were demonstrated as the alien introgressions in Xiaoyan 6 by multiple methods. Finally, the genetic transmission of the alien chromatin was analyzed in a set of wheat-Th. ponticum introgression lines. Some stable QTLs for morphological and physiological traits have been mapped near the alien chromatin.

Genetic dissection of flag leaf morphology traits and fine mapping of a novel QTL (Qflw.sxau-6BL) in bread wheat (Triticum aestivum L.)

KEY MESSAGE

Total 60-QRC for FLM traits were detected by meta-genomics analysis, nine major and stable QTL identified by DH population and validated, and a novel QTL  Qflw.sxau-6BL was fine mapped. The flag leaf is an "ideotypic" morphological trait providing photosynthetic assimilates in wheat. Although flag leaf morphology (FLM) traits had been extensively investigated through genetic mapping, there is a desire for FLM-related loci to be validated in multi-environments and fine mapping. In order to identify the stable genomic regions for FLM traits, we conducted a meta-genomic analysis based on reports from 2008 to 2024. Experimentally, a doubled haploid (DH) population was used to assess the genetic regions associated with FLM traits in nine environments. The meta-genomic analysis extracted 60 QTL-rich clusters (QRC), 45 of which were verified in marker-trait association (MTA) study. Nine major and stable QTL were found being associated with FLM traits across three-to-seven environments including BLUP, with phenotypic variance explained (PVE) ranging from 5.05 to 34.95%. The KASP markers of the nine QTL were validated (P < 0.005) in more than three environments using a panel of diverse wheat collections from Shanxi Province in China. Two co-located major and stable QTL viz. Qflw.sxau-6B.5 and Qfla.sxau-6B.4 were found novel and contributed to increase FLW by 12.09-19.21% and FLA by 5.45-13.28%. They also demonstrated high recombination rates in LD analysis based on the resequencing of 145 wheat landmark cultivars. The fine mapping of Qflw.sxau-6BL narrowed it down to a 1.27 Mb region as a result of the combined genotypic and phenotypic analysis for secondary mapping population. Comparing to NIL-ND3338, the NIL-LF5064 showed higher FLW by 20.45-27.37%, thousand-grain weight by 1.88-2.57% and grain length by 0.47-2.30% across all environments. The expression analysis of 11 tissues revealed seven highly expressed genes within the fine map region. This study provides a genetic basis for the FLM traits for further map-based cloning of FLW genes in wheat.

Refining polyploid breeding in sweet potato through allele dosage enhancement

Allele dosage plays a key role in the phenotypic variation of polyploids. Here we present a genome-wide variation map of hexaploid sweet potato that captures allele dosage information, constructed from deep sequencing of 294 hexaploid accessions. Genome-wide association studies identified quantitative trait loci with dosage effects on 23 agronomic traits. Our analyses reveal that sweet potato breeding has progressively increased the dosage of favourable alleles to enhance trait performance. Notably, the Mesoamerican gene pool has evolved towards higher dosages of favourable alleles at multiple loci, which have been increasingly introgressed into modern Chinese cultivars. We substantiated the breeding-driven dosage accumulation through transgenic validation of IbEXPA4, an expansin gene influencing tuberous root weight. In addition, we explored causative sequence variations that alter the expression of the Orange gene, which regulates flesh colour. Our findings illuminate the breeding history of sweet potato and establish a foundation for leveraging allele dosages in polyploid breeding practices.

Pan-genome bridges wheat structural variations with habitat and breeding

Wheat is the second largest food crop with a very good breeding system and pedigree record in China. Investigating the genomic footprints of wheat cultivars will unveil potential avenues for future breeding efforts1,2. Here we report chromosome-level genome assemblies of 17 wheat cultivars that chronicle the breeding history of China. Comparative genomic analysis uncovered a wealth of structural rearrangements, identifying 249,976 structural variations with 49.03% (122,567) longer than 5 kb. Cultivars developed in 1980s displayed significant accumulations of structural variations, a pattern linked to the extensive incorporation of European and American varieties into breeding programmes of that era. We further proved that structural variations in the centromere-proximal regions are associated with a reduction of crossover events. We showed that common wheat evolved from spring to winter types via mutations and duplications of the VRN-A1 gene as an adaptation strategy to a changing environment. We confirmed shifts in wheat cultivars linked to dietary preferences, migration and cultural integration in Northwest China. We identified large presence or absence variations of pSc200 tandem repeats on the 1RS terminal, suggesting its own rapid evolution in the wheat genome. The high-quality genome assemblies of 17 representatives developed and their good complementarity to the 10+ pan-genomes offer a robust platform for future genomics-assisted breeding in wheat.

Cytogenetic Identification and Molecular Marker Analysis of Two Wheat-Thinopyrum ponticum Translocations with Stripe Rust Resistance

Stripe rust, induced by Puccinia striiformis f. sp. tritici (Pst), is one of the most destructive fungal diseases of wheat worldwide. Thinopyrum ponticum, a significant wild relative for wheat improvement, exhibits innate immunity to this disease. To transfer the stripe rust resistance gene from Th. ponticum to wheat, two translocation lines, SN21171 and SN52684, were produced through distant hybridization techniques. Disease evaluation results showed that these two lines were immune to Pst species CYR32 at the adult plant stage. Molecular cytogenetic analyses and specific intron-targeting markers amplification results revealed that SN21171 and SN52684 harbor several T3Eb-3DS·3DL and T1Eb-1BS·1BL translocation chromosomes. Furthermore, the comparison of the chromosome karyotype from two translocation lines and their recurrent parent YN15, revealed that structural variation occurred in chromosomes 2A, 5A, 2B, 4B, 5B, and 6B in SN21171 and chromosomes 5A, 3B, 4B, 5B, 6B, and 7B in SN52684. Agronomic trait assessments uncovered advantageous properties in both lines, with SN21171 matching the recurrent parent and SN52684 exhibiting elevated higher grain number per main spike and increased thousand grain weight. These two translocation lines and specific markers may apply to wheat stripe rust-resistance breeding.

Single- and multi-locus genome-wide association study reveals genomic regions of thirteen yield-related traits in common wheat

Genetic dissection of yield-related traits can be used to improve wheat yield through molecular design breeding. In this study, we genotyped 245 wheat varieties and measured 13 yield-related plant height-, grain-, and spike-related traits, in seven environments, and identified 778 loci for these traits by genome-wide association study (GWAS) using single- and multi-locus models. Among these loci, nine were major, of which seven were novel, including Qph/lph.ahau-7A for plant height (PH) and leaf pillow height (LPH), Qngps/sps.ahau-1A for number of grains per spike (NGPS) and spikelet number per spike (SPS), Qsd.ahau-2B.1 and Qsd.ahau-5A.2 for spikelet density (SD), Qlph.ahau-7B.2 for LPH, Qgl.ahau-7B.3 for grain length (GL), and Qsl.ahau-3A.3 for spike length (SL). Through marker development, re-GWAS, gene annotation and cloning, and sequence variation, haplotype, and expression analyses, we confirmed two novel major loci and identified potential candidate genes, TraesCS7A02G118000 (named TaF-box-7A) and TraesCS1A02G190200 (named TaBSK2-1A) underlying Qph/lph.ahau-7A for PH-related traits and Qngps/sps.ahau-1A for spike-related traits. We also reported two favorable haplotypes, including TaF-box-Hap1 associated with low PH and LPH and TaBSK2-Hap3 associated with high NGPS and SPS. In summary, these findings can be applied to improve wheat yield and enrich our understanding of the complex genetic mechanisms of yield-related traits.

Exploring the genetic diversity and population structure of an ancient hexaploid wheat species Triticum sphaerococcum using SNP markers

BACKGROUND

Understanding genetic diversity and population structure is crucial for strategizing and enhancing breeding efficiency. Wheat, a globally cultivated crop, is a significant source of daily calories for humans. To overcome challenges such as extreme climatic fluctuations, stagnant yields, and diminishing genetic variation, it is essential to develop diverse germplasms with new alleles. Triticum sphaerococcum, an underutilized ancient hexaploid wheat species, shows promise for contributing beneficial alleles. However, the genetic diversity of its germplasms remains unstudied. This is the first report where we have examined the genetic diversity and population structure of 116 T. sphaerococcum accessions using a 35 K SNP Array. The objective of this study is to apply these findings to improve wheat breeding programs.

RESULTS

Analysis of the population's genetic structure identified four potential subpopulations, which was supported by principal coordinate analysis. Allele neutrality tests showed an abundance of intermediate genotypes, suggesting that many beneficial alleles are maintained through balancing selection. Among the three subgenomes, subgenome B exhibited the highest genetic diversity. AMOVA (Analysis of Molecular Variance) revealed significant variation both among (35%) and within (65%) the four subpopulations. The high genetic differentiation between subpopulations was corroborated by a moderate level of haploid migrant numbers (Nm = 1.286), indicating sufficient gene flow. SP4 emerged as the most diverse subpopulation, showing the highest values for allelic pattern indices due to its larger size and higher percentage of polymorphic loci. The D subgenome displayed a faster linkage disequilibrium (LD) decay rate compared to the A and B subgenomes. Haplotype block analysis identified 260 haplotype blocks of varying sizes distributed across the genome.

CONCLUSIONS

This research demonstrates that Indian dwarf wheat accessions, sourced from three distinct gene banks and local collections, possess considerable genetic diversity. These germplasm collections offer valuable opportunities to investigate their unexplored genetic potential. They can be utilized in wheat improvement initiatives to tackle both present and future breeding challenges. Furthermore, these accessions can introduce new alleles to broaden the genetic base of modern wheat varieties, enhancing their overall diversity.

QTL analysis in the Mianmai 902×Taichang 29 RIL population reveals the genetic basis for the high-yield of wheat cultivars Mianmai 902 in terms of spike and plant architecture

Spike-related traits and plant height (PH) are greatly associated with wheat yield. Identification of stable quantitative trait loci (QTL) for these traits is crucial for understanding the genetic basis for yield and their further application in breeding. In this study, QTL analysis for spikelet number per spike (SNS), spike length (SL), spike compactness (SC) and PH was performed using a recombinant inbred line (RIL) population derived from a cross between wheat cultivars Mianmai902 (MM902) and Taichang29 (TC29). Thirteen stably expressed QTLs were identified, and the most favorable alleles were contributed by MM902. The mjaor QTL, QSNS-MT.cib-2D, QSL-MT.cib-2D, QSC-MT.cib-2D, QSC-MT.cib-6 A, QPH-MT.cib-4B and QPH-MT.cib-4D, were repeatedly detected in multiple environments and explained 5.77-47.11% of the phenotypic variations. By employing the Kompetitive Allele Specific PCR markers, the most major QTLs were successfully validated in multiple populations derived from different genetic backgrounds. Moreover, the individual and pyramiding effects of major QTLs on SNS, SL, SC, PH, grain number per spike (GNS) and thousand-grain weight (TGW) were investigated and their potential utilization value in breeding was showcased. Ppd-D1 was predicted as the candidate for QSNS/SL/SC/PH-MT.cib-2D, and MM902 carried the Ppd-D1d allele. Ppd-D1d is rare in Chinese winter wheat cultivars and may be an elite alternative allele of the Ppd-D1a allele. In summary, these QTLs revealed the genetic basis of the spike-related traits and PH of MM902, which partially contributed to the high yield of MM902 and have application potential in wheat breeding by optimizing spike morphology and PH to improve yield.

Genomic insights into the modifications of spike morphology traits during wheat breeding

Over the past century, environmental changes have significantly impacted wheat spike morphology, crucial for adaptation and grain yield. However, the changes in wheat spike modifications during this period remain largely unknown. This study examines 16 spike morphology traits in 830 accessions released from 1900 to 2020. It finds that spike weight, grain number per spike (GN), and thousand kernel weight have significantly increased, while spike length has no significant change. The increase in fertile spikelets is due to fewer degenerated spikelets, resulting in a higher GN. Genome-wide association studies identified 49,994 significant SNPs, grouped into 293 genomic regions. The accumulation of favorable alleles in these genomic regions indicates the genetic basis for modification in spike morphology traits. Genetic network analysis of these genomic regions reveals the genetic basis for phenotypic correlations among spike morphology traits. The haplotypes of the identified genomic regions display obvious geographical differentiation in global accessions and environmental adaptation over the past 120 years. In summary, we reveal the genetic basis of adaptive evolution and the interactions of spike morphology, offering valuable resources for the genetic improvement of spike morphology to enhance environmental adaptation.

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.

Unveiling RNA structure-mediated regulations of RNA stability in wheat

Despite the critical role of mRNA stability in post-transcriptional gene regulation, research on this topic in wheat, a vital agricultural crop, remains unclear. Our study investigated the mRNA decay landscape of durum wheat (Triticum turgidum L. ssp. durum, BBAA), revealing subgenomic asymmetry in mRNA stability and its impact on steady-state mRNA abundance. Our findings indicate that the 3' UTR structure and homoeolog preference for RNA structural motifs can influence mRNA stability, leading to subgenomic RNA decay imbalance. Furthermore, single-nucleotide variations (SNVs) selected for RNA structural motifs during domestication can cause variations in subgenomic mRNA stability and subsequent changes in steady-state expression levels. Our research on the transcriptome stability of polyploid wheat highlights the regulatory role of non-coding region structures in mRNA stability, and how domestication shaped RNA structure, altering subgenomic mRNA stability. These results illustrate the importance of RNA structure-mediated post-transcriptional gene regulation in wheat and pave the way for its potential use in crop improvement.

On the evolution and genetic diversity of the bread wheat D genome

Bread wheat (Triticum aestivum) became a globally dominant crop after incorporating the D genome from the donor species Aegilops tauschii, but the evolutionary history that shaped the D genome during this process remains to be clarified. Here, we propose a renewed evolutionary model linking Ae. tauschii and the hexaploid wheat D genome by constructing an ancestral haplotype map covering 762 Ae. tauschii and hexaploid wheat accessions. We dissected the evolutionary trajectories of Ae. tauschii lineages and reported a few independent intermediate accessions, demonstrating that low-frequency inter-sublineage gene flow had enriched the diversity of Ae. tauschii. We discovered that the D genome of hexaploid wheat was inherited from a unified ancestral template, but with a mosaic composition that was highly mixed and derived mainly from three Ae. tauschii L2 sublineages located in the Caspian coastal region. This result suggests that early agricultural activities facilitated innovations in D-genome composition and finalized the success of hexaploidization. We found that the majority (51.4%) of genetic diversity was attributed to novel mutations absent in Ae. tauschii, and we identified large Ae. tauschii introgressions from various lineages, which expanded the diversity of the wheat D genome and introduced beneficial alleles. This work sheds light on the process of wheat hexaploidization and highlights the evolutionary significance of the multi-layered genetic diversity of the bread wheat D genome.

Genome wide association study and transcriptome analysis identify candidate genes regulating wheat coleoptile length

Coleoptile length, in wheat, is a significant agronomic trait impacting yield by facilitating the successful establishment of seedlings. In arid regions, varieties possessing longer coleoptile can evade harsh conditions by deep sowing, paving the way for improved yield. However, the study of genes involved in coleoptile development is insufficient. In this study, a high-density 660 K SNP array was used for genome-wide association study (GWAS) on coleoptile length in 150 wheat varieties. The findings revealed the detection of 353 significantly associated SNPs across all environments. The integration of linkage disequilibrium analysis and haplotype analysis mined 23 core QTLs capable responsible for the stable regulating coleoptile length in wheat. In wheat varieties characterized by extended coleoptile length, 6,600, 11,524, and 6,059 genes were found to be differentially expressed at three distinct developmental stages within the coleoptile, respectively. Through GWAS, gene expression levels, and functional annotation, we concluded the identification of two candidate genes (TraesCS2B02G423500, TraesCS2B02G449200) regulating wheat coleoptile length. By employing WGCNA and protein interactions prediction, discovered that the 19 genes were found to interact with candidate genes and participate in plant hormone metabolism and signaling, cell elongation or proliferation, which collectively contributing to coleoptile elongation. Additionally, two KASP markers were developed which can be used in breeding. These results offer a basis for understanding the genetic regulatory network responsible for wheat coleoptile length formation. The QTLs and candidate genes identified in this study can be further utilized for genetic improvement of wheat coleoptile length.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-024-01520-6.

Development and Characterization of Two Wheat-Rye Introgression Lines with Resistance to Stripe Rust and Powdery Mildew

Rye (Secale cereale L.) genes, which contribute to the tertiary gene pool of wheat, include multiple disease resistance genes useful for the genetic improvement of wheat. Introgression lines are the most valuable materials for wheat breeding because of their small alien segments and limited or lack of linkage drag. In the present study, wheat-rye derivative lines SN21627-2 and SN21627-6 were produced via distant hybridization. A genomic in situ hybridization analysis revealed that SN21627-2 and SN21627-6 lack alien segments, while a multi-color fluorescence in situ hybridization analysis detected structural changes in both introgression lines. At the seedling and adult plant stages, SN21627-2 and SN21627-6 were highly resistant to stripe rust and powdery mildew. Primers for 86 PCR-based landmark unique gene markers and 345 rye-specific SLAF markers were used to amplify SN21627-2 and SN21627-6 genomic DNA. Eight markers specific to rye chromosome 2R were detected in both introgression lines, implying these lines carry chromosome 2R segments with genes conferring stripe rust and powdery mildew resistance. Therefore, SN21627-2 and SN21627-6 are resistant to more than one major wheat disease, making them promising bridging parents for breeding disease-resistant wheat lines.

Pleiotropic phenotypic effects of the TaCYP78A family on multiple yield-related traits in wheat

Increasing crop yield depends on selecting and utilizing pleiotropic genes/alleles to improve multiple yield-related traits (YRTs) during crop breeding. However, synergistic improvement of YRTs is challenging due to the trade-offs between YRTs in breeding practices. Here, the favourable haplotypes of the TaCYP78A family are identified by analysing allelic variations in 1571 wheat accessions worldwide, demonstrating the selection and utilization of pleiotropic genes to improve yield and related traits during wheat breeding. The TaCYP78A family members, including TaCYP78A3, TaCYP78A5, TaCYP78A16, and TaCYP78A17, are organ size regulators expressed in multiple organs, and their allelic variations associated with various YRTs. However, due to the trade-offs between YRTs, knockdown or overexpression of TaCYP78A family members does not directly increase yield. Favourable haplotypes of the TaCYP78A family, namely A3/5/16/17Ap-Hap II, optimize the expression levels of TaCYP78A3/5/16/17-A across different wheat organs to overcome trade-offs and improve multiple YRTs. Different favourable haplotypes have both complementary and specific functions in improving YRTs, and their aggregation in cultivars under strong artificial selection greatly increase yield, even under various planting environments and densities. These findings provide new support and valuable genetic resources for molecular breeding of wheat and other crops in the era of Breeding 4.0.

Tracing post-domestication historical events and screening pre-breeding germplasm from large gene pools in wheat in the absence of phenotype data

Wheat, particularly common wheat (Triticum aestivum L.), is a major crop accounting for 25% of the world cereal production and thriving in diverse ecogeographic regions. Its adaptation to diverse environments arises from its three distinct genomes adapted to different environments and post-domestication anthropogenic interventions. In search of key genomic regions revealing historic events and breeding significance to common wheat, we performed genome scan and genome-environment association (GEA) analyses using high-marker density genotype datasets. Whole-genome scans revealed highly differentiated regions on chromosomes 2A, 3B, and 4A. In-depth analyses corroborated our previous prediction of the 4A differentiated region signifying the separation between Spelt/Macha and other wheat types. Individual chromosome scans captured key introgressions, including one from T. timopheevii and one from Thinopyrum ponticum on 2B and 3D, respectively, as well as known genes such as Vrn-A1 on 5A. GEA highlighted loci linked to latitude-induced environmental variations, influencing traits such as photoperiodism and responses to abiotic stress. Variation at the Vrn-A1 locus on 5A assigned accessions to two haplotypes (6% and 94%). Further analysis on Vrn-A1 coding gene revealed four subgroups of the major haplotype, while the minor haplotype remained undifferentiated. Analyses at differentiated loci mostly dichotomized the population, illustrating the possibility of isolating pre-breeding materials with desirable traits from large gene pools in the absence of phenotype data. Given the current availability of broad genetic data, the genome-scan-GEA hybrid can be an efficient and cost-effective approach for pinpointing environmentally resilient pre-breeding germplasm from vast gene pools, including gene banks regardless of trait characterization.

Genotype selection identified elite lines through quantitative trait loci mapping of agronomically important traits in wheat

Wheat is one of the most important staple foods in the world. Genetic characterization of wheat agronomically important traits is crucial for yield improvement through molecular breeding. In this study, a recombinant inbred line (RIL) population was developed by crossing a local adapted high yield variety Jimai 22 (JM22) with an external variety Cunmai no.1 (CM1). A high-density genetic map containing 7,359 single nucleotide polymorphism (SNP) markers was constructed. Quantitative trait loci (QTL) mapping identified 61 QTL for eight yield-related traits under six environments (years). Among them, 17 QTL affecting spike number per plant, grain number per spike and thousand grain weight showed high predictability for theoretical yield per plant (TYP), of which, 12 QTL alleles positively contributed to TYP. Nine promising candidate genes for seven of the 12 QTL were identified including three known wheat genes and six rice orthologs. Four elite lines with TYP increased by 5.6%-15.2% were identified through genotype selection which carried 7-9 favorable alleles from JM22 and 2-3 favorable alleles from CM1 of the 12 QTL. Moreover, the linked SNPs of the 12 QTL were converted to high-throughput kompetitive allele-specific PCR (KASP) markers and validated in the population. The mapped QTL, identified promising candidate genes, developed elite lines and KASP markers are highly valuable in future genotype selection to improve wheat yield.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-024-01496-3.

Genomic-wide analysis reveals seven in absentia genes regulating grain development in wheat (Triticum aestivum L.)

Seven in absentia proteins, which contain a conserved SINA domain, are involved in regulating various aspects of wheat (Triticum aestivum L.) growth and development, especially in response to environmental stresses. However, it is unclear whether TaSINA family members are involved in regulating grain development until now. In this study, the expression pattern, genomic polymorphism, and relationship with grain-related traits were analyzed for all TaSINA members. Most of the TaSINA genes identified showed higher expression levels in young wheat spikes or grains than other organs. The genomic polymorphism analysis revealed that at least 62 TaSINA genes had different haplotypes, where the haplotypes of five genes were significantly correlated with grain-related traits. Kompetitive allele-specific PCR markers were developed to confirm the single nucleotide polymorphisms in TaSINA101 and TaSINA109 among the five selected genes in a set of 292 wheat accessions. The TaSINA101-Hap II and TaSINA109-Hap II haplotypes had higher grain weight and width compared to TaSINA101-Hap I and TaSINA109-Hap I in at least three environments, respectively. The qRT-PCR assays revealed that TaSINA101 was highly expressed in the palea shell, seed coat, and embryo in young wheat grains. The TaSINA101 protein was unevenly distributed in the nucleus when transiently expressed in the protoplast of wheat. Three homozygous TaSINA101 transgenic lines in rice (Oryza sativa L.) showed higher grain weight and size compared to the wild type. These findings provide valuable insight into the biological function and elite haplotype of TaSINA family genes in wheat grain development at a genomic-wide level.

Change of Flavonoid Content in Wheatgrass in a Historic Collection of Wheat Cultivars

Wheatgrass is recognized for its nutritional and medicinal properties, partly attributed to its flavonoid content. The objective of this study was to assess the flavonoid content and antioxidant properties of wheatgrass obtained from a wide range of 145 wheat cultivars, which included Chinese landraces (CL), modern Chinese cultivars (MCC), and introduced modern cultivars (IMC). The flavonoids were extracted using a solution of 80% methanol, and their content was evaluated using ultra-high-performance liquid chromatography-mass spectrometry (UHPLC-MS). The results revealed the assessed cultivars showed significant variation in their total flavonoid content (TFC), with MCCs generally having higher amounts compared to CLs. PCA analysis demonstrated clear variations in flavonoid profiles between different cultivar groups, emphasizing the evolutionary inconsistencies in wheat breeding. The antioxidant assays, ABTS, DPPH, and FRAP, exhibited robust abilities for eliminating radicals, which were found to be directly associated with the amounts of flavonoids. In addition, this study investigated the correlation between the content of flavonoids and the ability to resist powdery mildew in a collection of mutated wheat plants. Mutants exhibiting heightened flavonoid accumulation demonstrated a decreased severity of powdery mildew, suggesting that flavonoids play a protective role against fungal infections. The results highlight the potential of wheatgrass as a valuable source of flavonoids that have antioxidant and protective effects. This potential is influenced by the genetic diversity and breeding history of wheatgrass. Gaining insight into these connections can guide future wheat breeding endeavors aimed at improving nutritional value and in strengthening disease resistance. The current finding provides critical information for developing wheatgrass with high flavonoid content and antioxidant activity.

Beyond pathways: Accelerated flavonoids candidate identification and novel exploration of enzymatic properties using combined mapping populations of wheat

Although forward-genetics-metabolomics methods such as mGWAS and mQTL have proven effective in providing myriad loci affecting metabolite contents, they are somehow constrained by their respective constitutional flaws such as the hidden population structure for GWAS and insufficient recombinant rate for QTL. Here, the combination of mGWAS and mQTL was performed, conveying an improved statistical power to investigate the flavonoid pathways in common wheat. A total of 941 and 289 loci were, respectively, generated from mGWAS and mQTL, within which 13 of them were co-mapped using both approaches. Subsequently, the mGWAS or mQTL outputs alone and their combination were, respectively, utilized to delineate the metabolic routes. Using this approach, we identified two MYB transcription factor encoding genes and five structural genes, and the flavonoid pathway in wheat was accordingly updated. Moreover, we have discovered some rare-activity-exhibiting flavonoid glycosyl- and methyl-transferases, which may possess unique biological significance, and harnessing these novel catalytic capabilities provides potentially new breeding directions. Collectively, we propose our survey illustrates that the forward-genetics-metabolomics approaches including multiple populations with high density markers could be more frequently applied for delineating metabolic pathways in common wheat, which will ultimately contribute to metabolomics-assisted wheat crop improvement.

Tagging large CNV blocks in wheat boosts digitalization of germplasm resources by ultra-low-coverage sequencing

BACKGROUND

The massive structural variations and frequent introgression highly contribute to the genetic diversity of wheat, while the huge and complex genome of polyploid wheat hinders efficient genotyping of abundant varieties towards accurate identification, management, and exploitation of germplasm resources.

RESULTS

We develop a novel workflow that identifies 1240 high-quality large copy number variation blocks (CNVb) in wheat at the pan-genome level, demonstrating that CNVb can serve as an ideal DNA fingerprinting marker for discriminating massive varieties, with the accuracy validated by PCR assay. We then construct a digitalized genotyping CNVb map across 1599 global wheat accessions. Key CNVb markers are linked with trait-associated introgressions, such as the 1RS·1BL translocation and 2NvS translocation, and the beneficial alleles, such as the end-use quality allele Glu-D1d (Dx5 + Dy10) and the semi-dwarf r-e-z allele. Furthermore, we demonstrate that these tagged CNVb markers promote a stable and cost-effective strategy for evaluating wheat germplasm resources with ultra-low-coverage sequencing data, competing with SNP array for applications such as evaluating new varieties, efficient management of collections in gene banks, and describing wheat germplasm resources in a digitalized manner. We also develop a user-friendly interactive platform, WheatCNVb ( http://wheat.cau.edu.cn/WheatCNVb/ ), for exploring the CNVb profiles over ever-increasing wheat accessions, and also propose a QR-code-like representation of individual digital CNVb fingerprint. This platform also allows uploading new CNVb profiles for comparison with stored varieties.

CONCLUSIONS

The CNVb-based approach provides a low-cost and high-throughput genotyping strategy for enabling digitalized wheat germplasm management and modern breeding with precise and practical decision-making.

A k-mer-based pangenome approach for cataloging seed-storage-protein genes in wheat to facilitate genotype-to-phenotype prediction and improvement of end-use quality

Wheat is a staple food for more than 35% of the world's population, with wheat flour used to make hundreds of baked goods. Superior end-use quality is a major breeding target; however, improving it is especially time-consuming and expensive. Furthermore, genes encoding seed-storage proteins (SSPs) form multi-gene families and are repetitive, with gaps commonplace in several genome assemblies. To overcome these barriers and efficiently identify superior wheat SSP alleles, we developed "PanSK" (Pan-SSP k-mer) for genotype-to-phenotype prediction based on an SSP-based pangenome resource. PanSK uses 29-mer sequences that represent each SSP gene at the pangenomic level to reveal untapped diversity across landraces and modern cultivars. Genome-wide association studies with k-mers identified 23 SSP genes associated with end-use quality that represent novel targets for improvement. We evaluated the effect of rye secalin genes on end-use quality and found that removal of ω-secalins from 1BL/1RS wheat translocation lines is associated with enhanced end-use quality. Finally, using machine-learning-based prediction inspired by PanSK, we predicted the quality phenotypes with high accuracy from genotypes alone. This study provides an effective approach for genome design based on SSP genes, enabling the breeding of wheat varieties with superior processing capabilities and improved end-use quality.

Combined genome and transcriptome analysis of elite fiber quality in Gossypium barbadense

Gossypium barbadense, which is one of several species of cotton, is well known for its superior fiber quality. However, the genetic basis of its high-quality fiber remains largely unexplored. Here, we resequenced 269 G. barbadense accessions. Phylogenetic structure analysis showed that the set of accessions was clustered into 3 groups: G1 and G2 mainly included modern cultivars from Xinjiang, China, and G3 was related to widely introduced accessions in different regions worldwide. A genome-wide association study of 5 fiber quality traits across multiple field environments identified a total of 512 qtls (main-effect QTLs) and 94 qtlEs (QTL-by-environment interactions) related to fiber quality, of which 292 qtls and 57 qtlEs colocated with previous studies. We extracted the genes located in these loci and performed expression comparison, local association analysis, and introgression segment identification. The results showed that high expression of hormone-related genes during fiber development, introgressions from Gossypium hirsutum, and the recombination of domesticated elite allelic variation were 3 major contributors to improve the fiber quality of G. barbadense. In total, 839 candidate genes with encoding region variations associated with elite fiber quality were mined. We confirmed that haplotype GB_D03G0092H traced to G. hirsutum introgression, with a 1-bp deletion leading to a frameshift mutation compared with GB_D03G0092B, significantly improved fiber quality. GB_D03G0092H is localized in the plasma membrane, while GB_D03G0092B is in both the nucleus and plasma membrane. Overexpression of GB_D03G0092H in Arabidopsis (Arabidopsis thaliana) significantly improved the elongation of longitudinal cells. Our study systematically reveals the genetic basis of the superior fiber quality of G. barbadense and provides elite segments and gene resources for breeding high-quality cotton cultivars.

Molecular characterization of QTL for grain zinc and iron concentrations in wheat landrace Chinese Spring

KEY MESSAGE

Three stable QTL for grain zinc concentration were identified in wheat landrace Chinese Spring. Favorable alleles were more frequent in landraces than in modern wheat cultivars. Wheat is a major source of dietary energy for the growing world population. Developing cultivars with enriched zinc and iron can potentially alleviate human micronutrient deficiency. In this study, a recombinant inbred line (RIL) population with 245 lines derived from cross Zhou 8425B/Chinese Spring was used to detect quantitative trait loci (QTL) for grain zinc concentration (GZnC) and grain iron concentration (GFeC) across four environments. Three stable QTL for GZnC with all favorable alleles from Chinese Spring were identified on chromosomes 3BL, 5AL, and 5BL. These QTL explaining maxima of 8.7%, 5.8%, and 7.1% of phenotypic variances were validated in 125 resequenced wheat accessions encompassing both landraces and modern cultivars using six kompetitive allele specific PCR (KASP) assays. The frequencies of favorable alleles for QGZnCzc.caas-3BL, QGZnCzc.caas-5AL and QGZnCzc.caas-5BL were higher in landraces (90.4%, 68.0%, and 100.0%, respectively) compared to modern cultivars (45.9%, 35.4%, and 40.9%), suggesting they were not selected in breeding programs. Candidate gene association studies on GZnC in the cultivar panel further delimited the QTL into 8.5 Mb, 4.1 Mb, and 47.8 Mb regions containing 46, 4, and 199 candidate genes, respectively. The 5BL QTL located in a region where recombination was suppressed. Two stable and three less stable QTL for GFeC with favorable alleles also from Chinese Spring were identified on chromosomes 4BS (Rht-B1a), 4DS (Rht-D1a), 1DS, 3AS, and 6DS. This study sheds light on the genetic basis of GZnC and GFeC in Chinese Spring and provides useful molecular markers for wheat biofortification.

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.

Degenerate oligonucleotide primer MIG-seq: an effective PCR-based method for high-throughput genotyping

Next-generation sequencing (NGS) library construction often involves using restriction enzymes to decrease genome complexity, enabling versatile polymorphism detection in plants. However, plant leaves frequently contain impurities, such as polyphenols, necessitating DNA purification before enzymatic reactions. To overcome this problem, we developed a PCR-based method for expeditious NGS library preparation, offering flexibility in number of detected polymorphisms. By substituting a segment of the simple sequence repeat sequence in the MIG-seq primer set (MIG-seq being a PCR method enabling library construction with low-quality DNA) with degenerate oligonucleotides, we introduced variability in detectable polymorphisms across various crops. This innovation, named degenerate oligonucleotide primer MIG-seq (dpMIG-seq), enabled a streamlined protocol for constructing dpMIG-seq libraries from unpurified DNA, which was implemented stably in several crop species, including fruit trees. Furthermore, dpMIG-seq facilitated efficient lineage selection in wheat and enabled linkage map construction and quantitative trait loci analysis in tomato, rice, and soybean without necessitating DNA concentration adjustments. These findings underscore the potential of the dpMIG-seq protocol for advancing genetic analyses across diverse plant species.

Genome-wide identification, gene expression and haplotype analysis of the rhomboid-like gene family in wheat (Triticum aestivum L.)

The rhomboid-like (RBL) gene encodes serine protease, which plays an important role in the response to cell development and diverse stresses. However, genome-wide identification, expression profiles, and haplotype analysis of the RBL family genes have not been performed in wheat (Triticum aestivum L.). This study investigated the phylogeny and diversity of the RBL family genes in the wheat genome through various approaches, including gene structure analysis, evolutionary relationship analysis, promoter cis-acting element analysis, expression pattern analysis, and haplotype analysis. The 41 TaRBL genes were identified and divided into five subfamilies in the wheat genome. RBL family genes were expanded through segmented duplication and purification selection. The cis-element analysis revealed their involvement in various stress responses and plant development. The results of RNA-seq and quantitative real-time-PCR showed that TaRBL genes displayed higher expression levels in developing spike/grain and were differentially regulated under polyethylene glycol, NaCl, and abscisic acid treatments, indicating their roles in grain development and abiotic stress response. A kompetitive allele-specific PCR molecular marker was developed to confirm the single nucleotide polymorphism of TaRBL14a gene in 263 wheat accessions. We found that the elite haplotype TaRBL14a-Hap2 showed a significantly higher 1000-grain weight than TaRBL14a-Hap11 in at least three environments, and the TaRBL14a-Hap2 was positively selected in wheat breeding. The findings will provide a good insight into the evolutionary and functional characteristics of the TaRBL genes family in wheat and lay the foundation for future exploration of the regulatory mechanisms of TaRBL genes in plant growth and development, as well as their response to abiotic stresses.

Variation of TaMyb10 and their function on grain color and pre-harvest sprouting resistance of wheat

Pre-harvest sprouting (PHS) is a significant threat to global food security due to its association with losses in both yield and quality. Among the genes involved in PHS resistance in wheat, PHS-3D (TaMyb10-D) plays a crucial role. Here, we characterized the sequence variations of TaMyb10 genes in 416 bread wheat and 302 Aegilops tauschii accessions. Within TaMyb10-A sequences, we identified a deletion ranging from 214 to 305 bp in the signal and amino acid coding region, present in 61.3% of the accessions. Similarly, 79.3% of the TaMyb10-B sequences within the third exon region exhibited a 19 bp deletion. Additionally, 40.8% of the accessions lacked the 2.4 Mb fragment (in/del mutations) on Chr3D, where TaMyb10-D/PHS-3D was located. Interestingly, the geographical distribution of accessions showed little correlation with the divergence of TaMyb10. TaMyb10-A-IIIDele, TaMyb10-B-IVDele, and TaMyb10-D-VDele genotypes were prevalent in wheat populations across continents. Despite their structural variations, the five distinct protein types exhibited comparable ability to bind the promoters of downstream genes in the flavonoid and ABA pathways, such as CHS, DFR, and NCED. Furthermore, the combination of TaMyb10 homologs was significantly associated with grain color and germination percentages. Accessions exclusively harboring TaMyb10-D displayed red seed color and reduced germination percentages, indicating the predominant role of TaMyb10-D compared to TaMyb10-A and TaMyb10-B. This comprehensive investigation enhances our understanding of the structural variations and functional divergence of TaMyb10, providing valuable insights and resources for improving PHS resistance in wheat.

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.

Dissecting the molecular basis of spike traits by integrating gene regulatory networks and genetic variation in wheat

Spike architecture influences both grain weight and grain number per spike, which are the two major components of grain yield in bread wheat (Triticum aestivum L.). However, the complex wheat genome and the influence of various environmental factors pose challenges in mapping the causal genes that affect spike traits. Here, we systematically identified genes involved in spike trait formation by integrating information on genomic variation and gene regulatory networks controlling young spike development in wheat. We identified 170 loci that are responsible for variations in spike length, spikelet number per spike, and grain number per spike through genome-wide association study and meta-QTL analyses. We constructed gene regulatory networks for young inflorescences at the double ridge stage and the floret primordium stage, in which the spikelet meristem and the floret meristem are predominant, respectively, by integrating transcriptome, histone modification, chromatin accessibility, eQTL, and protein-protein interactome data. From these networks, we identified 169 hub genes located in 76 of the 170 QTL regions whose polymorphisms are significantly associated with variation in spike traits. The functions of TaZF-B1, VRT-B2, and TaSPL15-A/D in establishment of wheat spike architecture were verified. This study provides valuable molecular resources for understanding spike traits and demonstrates that combining genetic analysis and developmental regulatory networks is a robust approach for dissection of complex traits.

Using UAV-Based Temporal Spectral Indices to Dissect Changes in the Stay-Green Trait in Wheat

Stay-green (SG) in wheat is a beneficial trait that increases yield and stress tolerance. However, conventional phenotyping techniques limited the understanding of its genetic basis. Spectral indices (SIs) as non-destructive tools to evaluate crop temporal senescence provide an alternative strategy. Here, we applied SIs to monitor the senescence dynamics of 565 diverse wheat accessions from anthesis to maturation stages over 2 field seasons. Four SIs (normalized difference vegetation index, green normalized difference vegetation index, normalized difference red edge index, and optimized soil-adjusted vegetation index) were normalized to develop relative stay-green scores (RSGS) as the SG indicators. An RSGS-based genome-wide association study identified 47 high-confidence quantitative trait loci (QTL) harboring 3,079 single-nucleotide polymorphisms associated with SG and 1,085 corresponding candidate genes. Among them, 15 QTL overlapped or were adjacent to known SG-related QTL/genes, while the remaining QTL were novel. Notably, a set of favorable haplotypes of SG-related candidate genes such as TraesCS2A03G1081100, TracesCS6B03G0356400, and TracesCS2B03G1299500 are increasing following the Green Revolution, further validating the feasibility of the pipeline. This study provided a valuable reference for further quantitative SG and genetic research in diverse wheat panels.

Comprehensive Comparative Analysis of the JAZ Gene Family in Common Wheat (Triticum aestivum) and Its D-Subgenome Donor Aegilops tauschii

JASMONATE-ZIM DOMAIN (JAZ) repressor proteins work as co-receptors in the jasmonic acid (JA) signalling pathway and are essential for plant development and environmental adaptation. Despite wheat being one of the main staple food crops, until recently, comprehensive analysis of its JAZ gene family has been limited due to the lack of complete and high-quality reference genomes. Here, using the latest reference genome, we identified 17 JAZ genes in the wheat D-genome donor Aegilops tauschii. Then, 54 TaJAZs were identified in common wheat. A systematic examination of the gene structures, conserved protein domains, and phylogenetic relationships of this gene family was performed. Five new JAZ genes were identified as being derived from tandem duplication after wheat divergence from other species. We integrated RNA-seq data and yield QTL information and found that tandemly duplicated TaJAZ genes were prone to association with spike-related traits. Moreover, 12 TaJAZ genes were located within breeding selection sweeps, including 9 tandemly duplicated ones. Haplotype variation analysis of selected JAZ genes showed significant association of TaJAZ7A and TaJAZ13A with thousand-grain weight. Our work provides a clearer picture of wheat JAZ gene evolution and puts forward the possibility of using these genes for wheat yield improvement.

Genome‑wide identification and expression analysis of the UBC gene family in wheat (Triticum aestivum L.)

BACKGROUND

Ubiquitination is an important regulatory step of selective protein degradation in the plant UPS (ubiquitin-proteasome system), which is involved in various biological processes in eukaryotes. Ubiquitin-conjugating enzymes play an intermediate role in the process of protein ubiquitination reactions and thus play an essential role in regulating plant growth and response to adverse environmental conditions. However, a genome-wide analysis of the UBC gene family in wheat (Triticum aestivum L.) has not yet been performed.

RESULTS

In this study, the number, physiochemical properties, gene structure, collinearity, and phylogenetic relationships of TaUBC family members in wheat were analyzed using bioinformatics methods. The expression pattern of TaUBC genes in different tissues/organs and developmental periods, as well as the transcript levels under abiotic stress treatment, were analyzed using RNA-Seq data and qRT-PCR. Meanwhile, favorable haplotypes of TaUBC25 were investigated based on wheat resequencing data of 681 wheat cultivars from the Wheat Union Database. The analyses identified a total of 93 TaUBC family members containing a UBC domain in wheat genome. These genes were unevenly distributed across 21 chromosomes, and numerous duplication events were observed between gene members. Based on phylogenetic analysis, the TaUBC family was divided into 13 E2 groups and a separate UEV group. We investigated the expression of TaUBC family genes under different tissue/organ and stress conditions by quantitative real-time PCR (qRT-PCR) analysis. The results showed that some TaUBC genes were specifically expressed in certain tissues/organs and that most TaUBC genes responded to NaCl, PEG6000, and ABA treatment with different levels of expression. In addition, we performed association analysis for the two haplotypes based on key agronomic traits such as thousand-kernel weight (TKW), kernel length (KL), kernel weight (KW), and kernel thickness (KT), examining 122 wheat accessions at three environmental sites. The results showed that TaUBC25-Hap II had significantly higher TKW, KL, KW, and KT than TaUBC25-Hap I. The distribution analysis of haplotypes showed that TaUBC25-Hap II was preferred in the natural population of wheat.

CONCLUSION

Our results identified 93 members of the TaUBC family in wheat, and several genes involved in grain development and abiotic stress response. Based on the SNPs detected in the TaUBC sequence, two haplotypes, TaUBC25-Hap I and TaUBC25-Hap II, were identified among wheat cultivars, and their potential value for wheat breeding was validated by association analysis. The above results provide a theoretical basis for elucidating the evolutionary relationships of the TaUBC gene family and lay the foundation for studying the functions of family members in the future.

Natural variation of STKc_GSK3 kinase TaSG-D1 contributes to heat stress tolerance in Indian dwarf wheat

Heat stress threatens global wheat (Triticum aestivum) production, causing dramatic yield losses worldwide. Identifying heat tolerance genes and comprehending molecular mechanisms are essential. Here, we identify a heat tolerance gene, TaSG-D1E286K, in Indian dwarf wheat (Triticum sphaerococcum), which encodes an STKc_GSK3 kinase. TaSG-D1E286K improves heat tolerance compared to TaSG-D1 by enhancing phosphorylation and stability of downstream target TaPIF4 under heat stress condition. Additionally, we reveal evolutionary footprints of TaPIF4 during wheat selective breeding in China, that is, InDels predominantly occur in the TaPIF4 promoter of Chinese modern wheat cultivars and result in decreased expression level of TaPIF4 in response to heat stress. These sequence variations with negative effect on heat tolerance are mainly introduced from European germplasm. Our study provides insight into heat stress response mechanisms and proposes a potential strategy to improve wheat heat tolerance in future.

Systematic identification of wheat spike developmental regulators by integrated multi-omics, transcriptional network, GWAS, and genetic analyses

The spike architecture of wheat plays a crucial role in determining grain number, making it a key trait for optimization in wheat breeding programs. In this study, we used a multi-omic approach to analyze the transcriptome and epigenome profiles of the young spike at eight developmental stages, revealing coordinated changes in chromatin accessibility and H3K27me3 abundance during the flowering transition. We constructed a core transcriptional regulatory network (TRN) that drives wheat spike formation and experimentally validated a multi-layer regulatory module involving TaSPL15, TaAGLG1, and TaFUL2. By integrating the TRN with genome-wide association studies, we identified 227 transcription factors, including 42 with known functions and 185 with unknown functions. Further investigation of 61 novel transcription factors using multiple homozygous mutant lines revealed 36 transcription factors that regulate spike architecture or flowering time, such as TaMYC2-A1, TaMYB30-A1, and TaWRKY37-A1. Of particular interest, TaMYB30-A1, downstream of and repressed by WFZP, was found to regulate fertile spikelet number. Notably, the excellent haplotype of TaMYB30-A1, which contains a C allele at the WFZP binding site, was enriched during wheat breeding improvement in China, leading to improved agronomic traits. Finally, we constructed a free and open access Wheat Spike Multi-Omic Database (http://39.98.48.156:8800/#/). Our study identifies novel and high-confidence regulators and offers an effective strategy for dissecting the genetic basis of wheat spike development, with practical value for wheat breeding.

The TaSOC1-TaVRN1 module integrates photoperiod and vernalization signals to regulate wheat flowering

Wheat needs different durations of vernalization, which accelerates flowering by exposure to cold temperature, to ensure reproductive development at the optimum time, as that is critical for adaptability and high yield. TaVRN1 is the central flowering regulator in the vernalization pathway and encodes a MADS-box transcription factor (TF) that usually works by forming hetero- or homo-dimers. We previously identified that TaVRN1 bound to an MADS-box TF TaSOC1 whose orthologues are flowering activators in other plants. The specific function of TaSOC1 and the biological implication of its interaction with TaVRN1 remained unknown. Here, we demonstrated that TaSOC1 was a flowering repressor in the vernalization and photoperiod pathways by overexpression and knockout assays. We confirmed the physical interaction between TaSOC1 and TaVRN1 in wheat protoplasts and in planta, and further validated their genetic interplay. A Flowering Promoting Factor 1-like gene TaFPF1-2B was identified as a common downstream target of TaSOC1 and TaVRN1 through transcriptome and chromatin immunoprecipitation analyses. TaSOC1 competed with TaVRT2, another MADS-box flowering regulator, to bind to TaVRN1; their coding genes synergistically control TaFPF1-2B expression and flowering initiation in response to photoperiod and low temperature. We identified major haplotypes of TaSOC1 and found that TaSOC1-Hap1 conferred earlier flowering than TaSOC1-Hap2 and had been subjected to positive selection in wheat breeding. We also revealed that wheat SOC1 family members were important domestication loci and expanded by tandem and segmental duplication events. These findings offer new insights into the regulatory mechanism underlying flowering control along with useful genetic resources for wheat improvement.

Innovative computational tools provide new insights into the polyploid wheat genome

Bread wheat (Triticum aestivum) is an important crop and serves as a significant source of protein and calories for humans, worldwide. Nevertheless, its large and allopolyploid genome poses constraints on genetic improvement. The complex reticulate evolutionary history and the intricacy of genomic resources make the deciphering of the functional genome considerably more challenging. Recently, we have developed a comprehensive list of versatile computational tools with the integration of statistical models for dissecting the polyploid wheat genome. Here, we summarize the methodological innovations and applications of these tools and databases. A series of step-by-step examples illustrates how these tools can be utilized for dissecting wheat germplasm resources and unveiling functional genes associated with important agronomic traits. Furthermore, we outline future perspectives on new advanced tools and databases, taking into consideration the unique features of bread wheat, to accelerate genomic-assisted wheat breeding.

A novel variation of TaGW2-6B increases grain weight without penalty in grain protein content in wheat (Triticum aestivum L.)

Yield and quality are two crucial breeding objects of wheat therein grain weight and grain protein content (GPC) are two key relevant factors correspondingly. Investigations of their genetic mechanisms represent special significance for breeding. In this study, 199 F2 plants and corresponding F2:3 families derived from Nongda3753 (ND3753) and its EMS-generated mutant 564 (M564) were used to investigate the genetic basis of larger grain and higher GPC of M564. QTL analysis identified a total of 33 environmentally stable QTLs related to thousand grain weight (TGW), grain area (GA), grain circle (GC), grain length (GL), grain width (GW), and GPC on chromosomes 1B, 2A, 2B, 4D, 6B, and 7D, respectively, among which QGw.cau-6B.1, QTgw.cau-6B.1, QGa.cau-6B.1, and QGc.cau-6B.1 shared overlap confidence interval on chromosome 6B. This interval contained the TaGW2 gene playing the same role as the QTLs, so TaGW2-6B was cloned and sequenced. Sequence alignment revealed two G/A SNPs between two parents, among which the SNP in the seventh exon led to a premature termination in M564. A KASP marker was developed based on the SNP, and single-marker analysis on biparental populations showed that the mutant allele could significantly increase GW and TGW, but had no effect on GPC. Distribution detection of the mutant allele through KASP marker genotyping and sequence alignment against databases ascertained that no materials harbored this allele within natural populations. This allele was subsequently introduced into three different varieties through molecular marker-assisted backcrossing, and it was revealed that the allele had a significant effect on simultaneously increasing GW, TGW, and even GPC in all of three backgrounds. Summing up the above, it could be concluded that a novel elite allele of TaGW2-6B was artificially created and might play an important role in wheat breeding for high yield and quality.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-024-01455-y.

Genetic dissection of major QTL for grain number per spike on chromosomes 5A and 6A in bread wheat (Triticum aestivum L.)

Grain number per spike (GNS) is a crucial component of grain yield and plays a significant role in improving wheat yield. To identify quantitative trait loci (QTL) associated with GNS, a recombinant inbred line (RIL) population derived from the cross of Zhongkemai 13F10 and Chuanmai 42 was employed to conduct QTL mapping across eight environments. Based on the bulked segregant exome sequencing (BSE-Seq), genomic regions associated with GNS were detected on chromosomes 5A and 6A. According to the constructed genetic maps, two major QTL QGns.cib-5A (LOD = 4.35-8.16, PVE = 8.46-14.43%) and QGns.cib-6A (LOD = 3.82-30.80, PVE = 5.44-12.38%) were detected in five and four environments, respectively. QGns.cib-6A is a QTL cluster for other seven yield-related traits. QGns.cib-5A and QGns.cib-6A were further validated using linked Kompetitive Allele Specific PCR (KASP) markers in different genetic backgrounds. QGns.cib-5A exhibited pleiotropic effects on productive tiller number (PTN), spike length (SL), fertile spikelet number per spike (FSN), and ratio of grain length to grain width (GL/GW) but did not significantly affect thousand grain weight (TGW). Haplotype analysis revealed that QGns.cib-5A and QGns.cib-6A were the targets of artificial selection during wheat improvement. Candidate genes for QGns.cib-5A and QGns.cib-6A were predicted by analyzing gene annotation, spatiotemporal expression patterns, and orthologous and sequence differences. These findings will be valuable for fine mapping and map-based cloning of genes underlying QGns.cib-5A and QGns.cib-6A.

Genome-wide identification and analysis of the GGCT gene family in wheat

BACKGROUND

γ-glutamylcyclotransferase (GGCT), an enzyme to maintain glutathione homeostasis, plays a vital role in the response to plant growth and development as well as the adaptation to various stresses. Although the GGCT gene family analysis has been conducted in Arabidopsis and rice, the family genes have not yet been well identified and analyzed at the genome-wide level in wheat (Triticum aestivum L.).

RESULTS

In the present study, 20 TaGGCT genes were identified in the wheat genome and widely distributed on chromosomes 2A, 2B, 2D, 3A, 4A, 5A, 5B, 5D, 6A, 6B, 6D, 7A, 7B, and 7D. Phylogenetic and structural analyses showed that these TaGGCT genes could be classified into three subfamilies: ChaC, GGGACT, and GGCT-PS. They exhibited similar motif compositions and distribution patterns in the same subgroup. Gene duplication analysis suggested that the expansion of TaGGCT family genes was facilitated by segmental duplications and tandem repeats in the wheat evolutionary events. Identification of diverse cis-acting response elements in TaGGCT promoters indicated their potential fundamental roles in response to plant development and abiotic stresses. The analysis of transcriptome data combined with RT-qPCR results revealed that the TaGGCTs genes exhibited ubiquitous expression across plant organs, with highly expressed in roots, stems, and developing grains. Most TaGGCT genes were up-regulated after 6 h under 20% PEG6000 and ABA treatments. Association analysis revealed that two haplotypes of TaGGCT20 gene displayed significantly different Thousand-kernel weight (TKW), Kernel length (KL), and Kernel width (KW) in wheat. The geographical and annual distribution of the two haplotypes of TaGGCT20 gene further revealed that the frequency of the favorable haplotype TaGGCT20-Hap-I was positively selected in the historical breeding process of wheat.

CONCLUSION

This study investigated the genome-wide identification, structure, evolution, and expression analysis of TaGGCT genes in wheat. The motifs of TaGGCTs were highly conserved throughout the evolutionary history of wheat. Most TaGGCT genes were highly expressed in roots, stems, and developing grains, and involved in the response to drought stresses. Two haplotypes were developed in the TaGGCT20 gene, where TaGGCT20-Hap-I, as a favorable haplotype, was significantly associated with higher TKW, KL, and KW in wheat, suggesting that the haplotype is used as a function marker for the selection in grain yield in wheat breeding.

Genome-Wide Association Study of Preharvest Sprouting in Wheat

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

In silico analysis of the Seven IN Absentia (SINA) genes in bread wheat sheds light on their structure in plants

Seven IN Absentia (SINA) is a small family of genes coding for ubiquitin-ligases that play major roles in regulating various plant growth and developmental processes, as well as in plant response to diverse biotic and abiotic stresses. Here, we studied the SINA genes family in bread wheat Triticum aestivum which is a culture of major importance for food security worldwide. One hundred and forty-one SINA family genes have been identified in bread wheat and showed that their number is very high compared to other plant species such as A. thaliana or rice. The expansion of this family seems to have been more important in monocots than in eudicots. In bread wheat, the chromosome 3 distal region is the site of a massive amplification of the SINA family, since we found that 83 of the 141 SINA genes are located on this chromosome in the Chinese Spring variety. This amplification probably occurred as a result of local duplications, followed by sequences divergence. The study was then extended to 4856 SINA proteins from 97 plant species. Phylogenetic and structural analyses identified a group of putative ancestral SINA proteins in plants containing a 58 aminoacid specific signature. Based on sequence homology and the research of that "Ancestral SINA motif" of 58 amino acids, a methodological process has been proposed and lead to the identification of functional SINA genes in a large family such as the Triticae that might be used for other species. Finally, tis paper gives a comprehensive overview of wheat gene family organization and functionalization taken the SINA genes as an example.

Epigenetic modifications regulate cultivar-specific root development and metabolic adaptation to nitrogen availability in wheat

The breeding of crops with improved nitrogen use efficiency (NUE) is crucial for sustainable agriculture, but the involvement of epigenetic modifications remains unexplored. Here, we analyze the chromatin landscapes of two wheat cultivars (KN9204 and J411) that differ in NUE under varied nitrogen conditions. The expression of nitrogen metabolism genes is closely linked to variation in histone modification instead of differences in DNA sequence. Epigenetic modifications exhibit clear cultivar-specificity, which likely contributes to distinct agronomic traits. Additionally, low nitrogen (LN) induces H3K27ac and H3K27me3 to significantly enhance root growth in KN9204, while remarkably inducing NRT2 in J411. Evidence from histone deacetylase inhibitor treatment and transgenic plants with loss function of H3K27me3 methyltransferase shows that changes in epigenetic modifications could alter the strategy preference for root development or nitrogen uptake in response to LN. Here, we show the importance of epigenetic regulation in mediating cultivar-specific adaptation to LN in wheat.

Genome resources for the elite bread wheat cultivar Aikang 58 and mining of elite homeologous haplotypes for accelerating wheat improvement

Despite recent progress in crop genomics studies, the genomic changes brought about by modern breeding selection are still poorly understood, thus hampering genomics-assisted breeding, especially in polyploid crops with compound genomes such as common wheat (Triticum aestivum). In this work, we constructed genome resources for the modern elite common wheat variety Aikang 58 (AK58). Comparative genomics between AK58 and the landrace cultivar Chinese Spring (CS) shed light on genomic changes that occurred through recent varietal improvement. We also explored subgenome diploidization and divergence in common wheat and developed a homoeologous locus-based genome-wide association study (HGWAS) approach, which was more effective than single homoeolog-based GWAS in unraveling agronomic trait-associated loci. A total of 123 major HGWAS loci were detected using a genetic population derived from AK58 and CS. Elite homoeologous haplotypes (HHs), formed by combinations of subgenomic homoeologs of the associated loci, were found in both parents and progeny, and many could substantially improve wheat yield and related traits. We built a website where users can download genome assembly sequence and annotation data for AK58, perform blast analysis, and run JBrowse. Our work enriches genome resources for wheat, provides new insights into genomic changes during modern wheat improvement, and suggests that efficient mining of elite HHs can make a substantial contribution to genomics-assisted breeding in common wheat and other polyploid crops.

Integration of genome-wide association study, linkage analysis, and population transcriptome analysis to reveal the TaFMO1-5B modulating seminal root growth in bread wheat

Bread wheat, one of the keystone crops for global food security, is challenged by climate change and resource shortage. The root system plays a vital role in water and nutrient absorption, making it essential for meeting the growing global demand. Here, using an association-mapping population composed of 406 accessions, we identified QTrl.Rs-5B modulating seminal root development with a genome-wide association study and validated its genetic effects with two F5 segregation populations. Transcriptome-wide association study prioritized TaFMO1-5B, a gene encoding the flavin-containing monooxygenases, as the causal gene for QTrl.Rs-5B, whose expression levels correlate negatively with the phenotyping variations among our population. The lines silenced for TaFMO1-5B consistently showed significantly larger seminal roots in different genetic backgrounds. Additionally, the agriculture traits measured in multiple environments showed that QTrl.Rs-5B also affects yield component traits and plant architecture-related traits, and its favorable haplotype modulates these traits toward that of modern cultivars, suggesting the application potential of QTrl.Rs-5B for wheat breeding. Consistently, the frequency of the favorable haplotype of QTrl.Rs-5B increased with habitat expansion and breeding improvement of bread wheat. In conclusion, our findings identified and demonstrated the effects of QTrl.Rs-5B on seminal root development and illustrated that it is a valuable genetic locus for wheat root improvement.