Optical maps refine the bread wheat Triticum aestivum cv. Chinese Spring genome assembly

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.

Genetic mapping of loci affecting seedling and adult-plant resistance to powdery mildew derived from two CIMMYT wheat lines

MAIN CONCLUSION

PM3 and PM8 alleles carried by two CIMMYT wheat lines confer powdery mildew resistance in seedlings and/or adult plants. A stage-specific epistatic interaction was observed between PM3 and PM8. Powdery mildew is an important foliar disease of wheat. Major genes for resistance, which have been widely used in wheat breeding programs, are typically effective against only limited numbers of virulence genes of the pathogen. The main aim of this study was to map resistance loci in wheat lines 7HRWSN58 and ZWW09-149 from the International Maize and Wheat Improvement Center (CIMMYT). Doubled haploid populations (Magenta/7HRWSN58 and Emu Rock/ZWW09-149) were developed and grown in controlled environment experiments and inoculated with a composite of Blumeria graminis f.sp. tritici isolates that had been collected at various locations in Western Australia. Plants were assessed for powdery mildew symptoms (percentage leaf area diseased) on seedlings and adult plants. Populations were subjected to genotyping-by-sequencing and assayed for known SNPs in the resistance gene PM3. Linkage maps were constructed, and markers were anchored to the wheat reference genome sequence. In both populations, there were asymptomatic lines that exhibited no symptoms. Among symptomatic lines, disease severity varied widely. In the Magenta/7HRWSN58 population, most of the observed variation was attributed to the PM3 region of chromosome 1A, with the allele from 7HRWSN58 conferring resistance in seedlings and adult plants. In the Emu Rock/ZWW09-149 population, two interacting quantitative trait loci were mapped: one at PM3 and the other on chromosome 1B. The Emu Rock/ZWW09-149 population was confirmed to segregate for a 1BL·1RS translocation that carries the PM8 powdery mildew resistance gene from rye. Consistent with previous reports that PM8-derived resistance can be suppressed by PM3 alleles, the observed interaction between the quantitative trait loci on chromosomes 1A and 1B indicated that the PM3 allele carried by ZWW09-149 suppresses PM8-derived resistance from ZWW09-149, but only at the seedling stage. In adult plants, the PM8 region conferred resistance regardless of the PM3 genotype. The resistance sources and molecular markers that were investigated here could be useful in wheat breeding.

Comprehensive genome-wide analysis of wheat xylanase inhibitor protein (XIP) genes: unveiling their role in Fusarium head blight resistance and plant immune mechanisms

In this comprehensive genome-wide study, we identified and classified 83 Xylanase Inhibitor Protein (XIP) genes in wheat, grouped into five distinct categories, to enhance understanding of wheat's resistance to Fusarium head blight (FHB), a significant fungal threat to global wheat production. Our analysis reveals the unique distribution of XIP genes across wheat chromosomes, particularly at terminal regions, suggesting their role in the evolutionary expansion of the gene family. Several XIP genes lack signal peptides, indicating potential alternative secretion pathways that could be pivotal in plant defense against FHB. The study also uncovers the sequence homology between XIPs and chitinases, hinting at a functional diversification within the XIP gene family. Additionally, the research explores the association of XIP genes with plant immune mechanisms, particularly their linkage with plant hormone signaling pathways like abscisic acid and jasmonic acid. XIP-7A3, in particular, demonstrates a significant increase in expression upon FHB infection, highlighting its potential as a key candidate gene for enhancing wheat's resistance to this disease. This research not only enriches our understanding of the XIP gene family in wheat but also provides a foundation for future investigations into their role in developing FHB-resistant wheat cultivars. The findings offer significant implications for wheat genomics and breeding, contributing to the development of more resilient crops against fungal diseases.

Characterization of a new greenbug resistance gene Gb9 in a synthetic hexaploid wheat

Greenbug [Schizaphis graminum (Rondani)] is a serious insect pest that not only damages cereal crops, but also transmits several destructive viruses. The emergence of new greenbug biotypes in the field makes it urgent to identify novel greenbug resistance genes in wheat. CWI 76364 (PI 703397), a synthetic hexaploid wheat (SHW) line, exhibits greenbug resistance. Evaluation of an F2:3 population from cross OK 14319 × CWI 76364 indicated that a dominant gene, designated Gb9, conditions greenbug resistance in CWI 76364. Selective genotyping of a subset of F2 plants with contrasting phenotypes by genotyping-by-sequencing identified 25 SNPs closely linked to Gb9 on chromosome arm 7DL. Ten of these SNPs were converted to Kompetitive allele-specific polymerase chain reaction (KASP) markers for genotyping the entire F2 population. Genetic analysis delimited Gb9 to a 0.6-Mb interval flanked by KASP markers located at 599,835,668 bp (Stars-KASP872) and 600,471,081 bp (Stars-KASP881) on 7DL. Gb9 was 0.5 cM distal to Stars-KASP872 and 0.5 cM proximal to Stars-KASP881. Allelism tests indicated that Gb9 is a new greenbug resistance gene which confers resistance to greenbug biotypes C, E, H, I, and TX1. TX1 is one of the most widely virulent biotypes and has overcome most known wheat greenbug resistance genes. The introgression of Gb9 into locally adapted wheat cultivars is of economic importance, and the KASP markers developed in this study can be used to tag Gb9 in cultivar development.

Genome-wide association study for seedling heat tolerance under two temperature conditions in bread wheat (Triticum aestivum L.)

BACKGROUND

As the greenhouse effect intensifies, global temperatures are steadily increasing, posing a challenge to bread wheat (Triticum aestivum L.) production. It is imperative to comprehend the mechanism of high temperature tolerance in wheat and implement breeding programs to identify and develop heat-tolerant wheat germplasm and cultivars.

RESULTS

To identify quantitative trait loci (QTL) related to heat stress tolerance (HST) at seedling stage in wheat, a panel of 253 wheat accessions which were re-sequenced used to conduct genome-wide association studies (GWAS) using the factored spectrally transformed linear mixed models (FaST-LMM). For most accessions, the growth of seedlings was found to be inhibited under heat stress. Analysis of the phenotypic data revealed that under heat stress conditions, the main root length, total root length, and shoot length of seedlings decreased by 47.46%, 49.29%, and 15.19%, respectively, compared to those in normal conditions. However, 17 varieties were identified as heat stress tolerant germplasm. Through GWAS analysis, a total of 115 QTLs were detected under both heat stress and normal conditions. Furthermore, 15 stable QTL-clusters associated with heat response were identified. By combining gene expression, haplotype analysis, and gene annotation information within the physical intervals of the 15 QTL-clusters, two novel candidate genes, TraesCS4B03G0152700/TaWRKY74-B and TraesCS4B03G0501400/TaSnRK3.15-B, were responsive to temperature and identified as potential regulators of HST in wheat at the seedling stage.

CONCLUSIONS

This study conducted a detailed genetic analysis and successfully identified two genes potentially associated with HST in wheat at the seedling stage, laying a foundation to further dissect the regulatory mechanism underlying HST in wheat under high temperature conditions. Our finding could serve as genomic landmarks for wheat breeding aimed at improving adaptation to heat stress in the face of climate change.

Future of durum wheat research and breeding: Insights from early career researchers

Durum wheat (Triticum turgidum ssp. durum) is globally cultivated for pasta, couscous, and bulgur production. With the changing climate and growing world population, the need to significantly increase durum production to meet the anticipated demand is paramount. This review summarizes recent advancements in durum research, encompassing the exploitation of existing and novel genetic diversity, exploration of potential new diversity sources, breeding for climate-resilient varieties, enhancements in production and management practices, and the utilization of modern technologies in breeding and cultivar development. In comparison to bread wheat (T. aestivum), the durum wheat community and production area are considerably smaller, often comprising many small-family farmers, notably in African and Asian countries. Public breeding programs such as the International Maize and Wheat Improvement Center (CIMMYT) and the International Center for Agricultural Research in the Dry Areas (ICARDA) play a pivotal role in providing new and adapted cultivars for these small-scale growers. We spotlight the contributions of these and others in this review. Additionally, we offer our recommendations on key areas for the durum research community to explore in addressing the challenges posed by climate change while striving to enhance durum production and sustainability. As part of the Wheat Initiative, the Expert Working Group on Durum Wheat Genomics and Breeding recognizes the significance of collaborative efforts in advancing toward a shared objective. We hope the insights presented in this review stimulate future research and deliberations on the trajectory for durum wheat genomics and breeding.

Bulked segregant RNA-seq reveals complex resistance expression profile to powdery mildew in wild emmer wheat W762

Powdery mildew, caused by Blumeria graminis f. sp. tritici (Bgt), is one of the most destructive fungal diseases threatening global wheat production. Exploring powdery mildew resistance (Pm) gene(s) and dissecting the molecular mechanism of the host resistance are critical to effectively and reasonably control this disease. Durum wheat (Triticum turgidum L. var. durumDesf.) is an important gene donor for wheat improvement against powdery mildew. In this study, a resistant durum wheat accession W762 was used to investigate its potential resistance component(s) and profile its expression pattern in responding to Bgt invasion using bulked segregant RNA-Seq (BSR-Seq) and further qRT-PCR verification. Genetic analysis showed that the powdery mildew resistance in W762 did not meet monogenic inheritance and complex genetic model might exist within the population of W762 × Langdon (susceptible durum wheat). After BSR-Seq, 6,196 consistently different single nucleotide polymorphisms (SNPs) were called between resistant and susceptible parents and bulks, and among them, 763 SNPs were assigned to the chromosome arm 7B. Subsequently, 3,653 differentially expressed genes (DEGs) between resistant and susceptible parents and bulks were annotated and analyzed by Gene Ontology (GO), Cluster of Orthologous Groups (COG), and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment. The potential regulated genes were selected and analyzed their temporal expression patterns following Bgt inoculation. As a result, nine disease-related genes showed distinctive expression profile after Bgt invasion and might serve as potential targets to regulate the resistance against powdery mildew in W762. Our study could lay a foundation for analysis of the molecular mechanism and also provide potential targets for the improvement of durable resistance against powdery mildew.

Deacclimation and reacclimation processes in winter wheat: novel perspectives from time-series transcriptome analysis

Winter wheat achieves freezing tolerance (FT) through cold acclimation (CA) - a process which is induced by low positive temperatures in autumn. The increasing occurrences of temperature fluctuations in winter lead to deacclimation (DEA), causing premature loss of FT, and the cultivars capable of reacclimation (REA) are more likely to survive the subsequent cold spells. The genetic mechanisms of DEA and REA remain poorly understood, necessitating further research to bolster climate resilience in winter wheat. Here, we selected two winter wheat genotypes with contrasting levels of FT and conducted a ten-week-long experiment imitating low-temperature fluctuations after CA under controlled conditions. Crown and leaf tissue samples for RNA-sequencing were collected at CA, DEA, and REA time-points. It is the first transcriptomic study covering both short- and long-term responses to DEA and REA in winter wheat. The study provides novel knowledge regarding CA, DEA, and REA and discusses the gene expression patterns conferring FT under temperature fluctuations. The freezing-tolerant genotype "Lakaja DS" showed elevated photosynthetic activity in leaf tissue and upregulated cryoprotective protein-encoding genes in crowns after CA when compared to the freezing-susceptible "KWS Ferrum". "Lakaja DS" also expressed cold acclimation-associated transcripts at a significantly higher level after 1 week of DEA. Following REA, "Lakaja DS" continued to upregulate dehydrin-related genes in crowns and exhibited significantly higher expression of chitinase transcripts in leaves, when compared to "KWS Ferrum". The findings of this study shed light on the genetic mechanisms governing DEA and REA in winter wheat, thus addressing the gaps in knowledge regarding FT under low-temperature fluctuations. The identified genes should be further examined as potential molecular markers for breeding strategies focused on developing freezing-tolerant winter-type crops. Publicly available datasets generated in this study are valuable resources for further research into DEA and REA, contributing towards the enhancement of winter wheat under global climate change.

Identification of Novel Loci Precisely Modulating Pre-Harvest Sprouting Resistance and Red Color Components of the Seed Coat in T. aestivum L

The association between pre-harvest sprouting (PHS) and seed coat color has long been recognized. Red-grained wheats generally exhibit greater PHS resistance compared to white-grained wheat, although variability in PHS resistance exists within red-grained varieties. Here, we conducted a genome-wide association study on a panel consisting of red-grained wheat varieties, aimed at uncovering genes that modulate PHS resistance and red color components of seed coat using digital image processing. Twelve loci associated with PHS traits were identified, nine of which were described for the first time. Genetic loci marked by SNPs AX-95172164 (chromosome 1B) and AX-158544327 (chromosome 7D) explained approximately 25% of germination index variance, highlighting their value for breeding PHS-resistant varieties. The most promising candidate gene for PHS resistance was TraesCS6B02G147900, encoding a protein involved in aleurone layer morphogenesis. Twenty-six SNPs were significantly associated with grain color, independently of the known Tamyb10 gene. Most of them were related to multiple color characteristics. Prioritization of genes within the revealed loci identified TraesCS1D03G0758600 and TraesCS7B03G1296800, involved in the regulation of pigment biosynthesis and in controlling pigment accumulation. In conclusion, our study identifies new loci associated with grain color and germination index, providing insights into the genetic mechanisms underlying these traits.

Mapping of Leaf Rust Resistance Loci in Two Kenyan Wheats and Development of Linked Markers

Leaf rust caused by the pathogen Puccinia triticina (Pt) is a destructive fungal disease of wheat that occurs in almost all wheat-growing areas across the globe. Genetic resistance has proven to be the best solution to mitigate the disease. Wheat breeders are continuously seeking new diversified and durable sources of resistance to use in developing new varieties. We developed recombinant inbred line (RIL) populations from two leaf rust-resistant genotypes (Kenya Kudu and AUS12568) introduced from Kenya to identify and characterize resistance to Pt and to develop markers linked closely to the resistance that was found. Our studies detected four QTL conferring adult plant resistance (APR) to leaf rust. Two of these loci are associated with known genes, Lr46 and Lr68, residing on chromosomes 1B and 7B, respectively. The remaining two, QLrKK_2B and QLrAus12568_5A, contributed by Kenya Kudu and AUS12568 respectively, are putatively new loci for Pt resistance. Both QLrKK_2B and QLrAus12568_5A were found to interact additively with Lr46 in significantly reducing the disease severity at adult plant growth stages in the field. We further developed a suite of six closely linked markers within the QLrAus12568_5A locus and four within the QLrKK_2B region. Among these, markers sunKASP_522 and sunKASP_524, flanking QLrAus12568_5A, and sunKASP_536, distal to QLrKK_2B, were identified as the most closely linked and reliable for marker-assisted selection. The markers were validated on a selection of 64 Australian wheat varieties and found to be polymorphic and robust, allowing for clear allelic discrimination. The identified new loci and linked molecular markers will enable rapid adoption by breeders in developing wheat varieties carrying diversified and durable resistance to leaf rust.

Chromosome-scale genome assembly of bread wheat's wild relative Triticum timopheevii

Wheat (Triticum aestivum) is one of the most important food crops with an urgent need for increase in its production to feed the growing world. Triticum timopheevii (2n = 4x = 28) is an allotetraploid wheat wild relative species containing the At and G genomes that has been exploited in many pre-breeding programmes for wheat improvement. In this study, we report the generation of a chromosome-scale reference genome assembly of T. timopheevii accession PI 94760 based on PacBio HiFi reads and chromosome conformation capture (Hi-C). The assembly comprised a total size of 9.35 Gb, featuring a contig N50 of 42.4 Mb and included the mitochondrial and plastid genome sequences. Genome annotation predicted 166,325 gene models including 70,365 genes with high confidence. DNA methylation analysis showed that the G genome had on average more methylated bases than the At genome. In summary, the T. timopheevii genome assembly provides a valuable resource for genome-informed discovery of agronomically important genes for food security.

QTL mapping for the flag leaf-related traits using RILs derived from Trititrigia germplasm line SN304 and wheat cultivar Yannong15 in multiple environments

BACKGROUND

Developing and enriching genetic resources plays important role in the crop improvement. The flag leaf affects plant architecture and contributes to the grain yield of wheat (Triticum aestivum L.). The genetic improvement of flag leaf traits faces problems such as a limited genetic basis. Among the various genetic resources of wheat, Thinopyrum intermedium has been utilized as a valuable resource in genetic improvement due to its disease resistance, large spikes, large leaves, and multiple flowers. In this study, a recombinant inbred line (RIL) population was derived from common wheat Yannong15 and wheat-Th. intermedium introgression line SN304 was used to identify the quantitative trait loci (QTL) for flag leaf-related traits.

RESULTS

QTL mapping was performed for flag leaf length (FLL), flag leaf width (FLW) and flag leaf area (FLA). A total of 77 QTLs were detected, and among these, 51 QTLs with positive alleles were contributed by SN304. Fourteen major QTLs for flag leaf traits were detected on chromosomes 2B, 3B, 4B, and 2D. Additionally, 28 QTLs and 8 QTLs for flag leaf-related traits were detected in low-phosphorus and drought environments, respectively. Based on major QTLs of positive alleles from SN304, we identified a pair of double-ended anchor primers mapped on chromosome 2B and amplified a specific band of Th. intermedium in SN304. Moreover, there was a major colocated QTL on chromosome 2B, called QFll/Flw/Fla-2B, which was delimited to a physical interval of approximately 2.9 Mb and contained 20 candidate genes. Through gene sequence and expression analysis, four candidate genes associated with flag leaf formation and growth in the QTL interval were identified.

CONCLUSION

These results promote the fine mapping of QFll/Flw/Fla-2B, which have pleiotropic effects, and will facilitate the identification of candidate genes for flag leaf-related traits. Additionally, this work provides a theoretical basis for the application of Th. intermedium in wheat breeding.

Non-B-form DNA is associated with centromere stability in newly-formed polyploid wheat

Non-B-form DNA differs from the classic B-DNA double helix structure and plays a crucial regulatory role in replication and transcription. However, the role of non-B-form DNA in centromeres, especially in polyploid wheat, remains elusive. Here, we systematically analyzed seven non-B-form DNA motif profiles (A-phased DNA repeat, direct repeat, G-quadruplex, inverted repeat, mirror repeat, short tandem repeat, and Z-DNA) in hexaploid wheat. We found that three of these non-B-form DNA motifs were enriched at centromeric regions, especially at the CENH3-binding sites, suggesting that non-B-form DNA may create a favorable loading environment for the CENH3 nucleosome. To investigate the dynamics of centromeric non-B form DNA during the alloploidization process, we analyzed DNA secondary structure using CENH3 ChIP-seq data from newly formed allotetraploid wheat and its two diploid ancestors. We found that newly formed allotetraploid wheat formed more non-B-form DNA in centromeric regions compared with their parents, suggesting that non-B-form DNA is related to the localization of the centromeric regions in newly formed wheat. Furthermore, non-B-form DNA enriched in the centromeric regions was found to preferentially form on young LTR retrotransposons, explaining CENH3's tendency to bind to younger LTR. Collectively, our study describes the landscape of non-B-form DNA in the wheat genome, and sheds light on its potential role in the evolution of polyploid centromeres.

A bird's-eye view: exploration of the flavin-containing monooxygenase superfamily in common wheat

The Flavin Monooxygenase (FMO) gene superfamily in plants is involved in various processes most widely documented for its involvement in auxin biosynthesis, specialized metabolite biosynthesis, and plant microbial defense signaling. The roles of FMOs in defense signaling and disease resistance have recently come into focus as they may present opportunities to increase immune responses in plants including leading to systemic acquired resistance, but are not well characterized. We present a comprehensive catalogue of FMOs found in genomes across vascular plants and explore, in depth, 170 wheat TaFMO genes for sequence architecture, cis-acting regulatory elements, and changes due to Transposable Element insertions. A molecular phylogeny separates TaFMOs into three clades (A, B, and C) for which we further report gene duplication patterns, and differential rates of homoeologue expansion and retention among TaFMO subclades. We discuss Clade B TaFMOs where gene expansion is similarly seen in other cereal genomes. Transcriptome data from various studies point towards involvement of subclade B2 TaFMOs in disease responses against both biotrophic and necrotrophic pathogens, substantiated by promoter element analysis. We hypothesize that certain TaFMOs are responsive to both abiotic and biotic stresses, providing potential targets for enhancing disease resistance, plant yield and other important agronomic traits. Altogether, FMOs in wheat and other crop plants present an untapped resource to be exploited for improving the quality of crops.

CRISPR-based editing of the ω- and γ-gliadin gene clusters reduces wheat immunoreactivity without affecting grain protein quality

Wheat immunotoxicity is associated with abnormal reaction to gluten-derived peptides. Attempts to reduce immunotoxicity using breeding and biotechnology often affect dough quality. Here, the multiplexed CRISPR-Cas9 editing of cultivar Fielder was used to modify gluten-encoding genes, specifically focusing on ω- and γ-gliadin gene copies, which were identified to be abundant in immunoreactive peptides based on the analysis of wheat genomes assembled using the long-read sequencing technologies. The whole-genome sequencing of an edited line showed mutation or deletion of nearly all ω-gliadin and half of the γ-gliadin gene copies and confirmed the lack of editing in the α/β-gliadin genes. The estimated 75% and 64% reduction in ω- and γ-gliadin content, respectively, had no negative impact on the end-use quality characteristics of grain protein and dough. A 47-fold immunoreactivity reduction compared to a non-edited line was demonstrated using antibodies against immunotoxic peptides. Our results indicate that the targeted CRISPR-based modification of the ω- and γ-gliadin gene copies determined to be abundant in immunoreactive peptides by analysing high-quality genome assemblies is an effective mean for reducing immunotoxicity of wheat cultivars while minimizing the impact of editing on protein quality.

Integrated VIS/NIR Spectrum and Genome-Wide Association Study for Genetic Dissection of Cellulose Crystallinity in Wheat Stems

Cellulose crystallinity is a crucial factor influencing stem strength and, consequently, wheat lodging. However, the genetic dissection of cellulose crystallinity is less reported due to the difficulty of its measurement. In this study, VIS/NIR spectra and cellulose crystallinity were measured for a wheat accession panel with diverse genetic backgrounds. We developed a reliable VIS/NIR model for cellulose crystallinity with a high determination coefficient (R2) (0.95) and residual prediction deviation (RPD) (4.04), enabling the rapid screening of wheat samples. A GWAS of the cellulose crystallinity in 326 wheat accessions revealed 14 significant SNPs and 13 QTLs. Two candidate genes, TraesCS4B03G0029800 and TraesCS5B03G1085500, were identified. In summary, this study establishes an efficient method for the measurement of cellulose crystallinity in wheat stems and provides a genetic basis for enhancing lodging resistance in wheat.

Fine mapping of a major QTL, qKl-1BL controlling kernel length in common wheat

KEY MESSAGE

A major stable QTL, qKl-1BL, for kernel length of wheat was narrowed down to a 2.04-Mb interval on chromosome 1BL; the candidate genes were predicated and the genetic effects on yield-related traits were characterized. As a key factor influencing kernel weight, wheat kernel shape is closely related to yield formation, and in turn affects both wheat processing quality and market value. Fine mapping of the major quantitative trait loci (QTL) for kernel shape could provide genetic resources and a theoretical basis for the genetic improvement of wheat yield-related traits. In this study, a major QTL for kernel length (KL) on 1BL, named qKl-1BL, was identified from the recombinant inbred lines (RIL) in multiple environments based on the genetic map and physical map, with 4.76-21.15% of the phenotypic variation explained. To fine map qKl-1BL, the map-based cloning strategy was used. By using developed InDel markers, the near-isogenic line (NIL) pairs and eight key recombinants were identified from a segregating population containing 3621 individuals derived from residual heterozygous lines (RHLs) self-crossing. In combination with phenotype identification, qKl-1BL was finely positioned into a 2.04-Mb interval, KN1B:698.15-700.19 Mb, with eight differentially expressed genes enriched at the key period of kernel elongation. Based on transcriptome analysis and functional annotation information, two candidate genes for qKl-1BL controlling kernel elongation were identified. Additionally, genetic effect analysis showed that the superior allele of qKl-1BL from Jing411 could increase KL, thousand kernel weight (TKW), and yield per plant (YPP) significantly, as well as kernel bulk density and stability time. Taken together, this study identified a QTL interval for controlling kernel length with two possible candidate genes, which provides an important basis for qKl-1BL cloning, functional analysis, and application in molecular breeding programs.

Genome-wide association mapping for pre-harvest sprouting in European winter wheat detects novel resistance QTL, pleiotropic effects, and structural variation in multiple genomes

Pre-harvest sprouting (PHS), germination of seeds before harvest, is a major problem in global wheat (Triticum aestivum L.) production, and leads to reduced bread-making quality in affected grain. Breeding for PHS resistance can prevent losses under adverse conditions. Selecting resistant lines in years lacking pre-harvest rain, requires challenging of plants in the field or in the laboratory or using genetic markers. Despite the availability of a wheat reference and pan-genome, linking markers, genes, allelic, and structural variation, a complete understanding of the mechanisms underlying various sources of PHS resistance is still lacking. Therefore, we challenged a population of European wheat varieties and breeding lines with PHS conditions and phenotyped them for PHS traits, grain quality, phenological and agronomic traits to conduct genome-wide association mapping. Furthermore, we compared these marker-trait associations to previously reported PHS loci and evaluated their usefulness for breeding. We found markers associated with PHS on all chromosomes, with strong evidence for novel quantitative trait locus/loci (QTL) on chromosome 1A and 5B. The QTL on chromosome 1A lacks pleiotropic effect, for the QTL on 5B we detected pleiotropic effects on phenology and grain quality. Multiple peaks on chromosome 4A co-located with the major resistance locus Phs-A1, for which two causal genes, TaPM19 and TaMKK3, have been proposed. Mapping markers and genes to the pan-genome and chromosomal alignments provide evidence for structural variation around this major PHS-resistance locus. Although PHS is controlled by many loci distributed across the wheat genome, Phs-A1 on chromosome 4A seems to be the most effective and widely deployed source of resistance, in European wheat varieties.

Genome-wide association studies reveal stable loci for wheat grain size under different sowing dates

Background

Wheat (Tritium aestivum L.) production is critical for global food security. In recent years, due to climate change and the prolonged growing period of rice varieties, the delayed sowing of wheat has resulted in a loss of grain yield in the area of the middle and lower reaches of the Yangtze River. It is of great significance to screen for natural germplasm resources of wheat that are resistant to late sowing and to explore genetic loci that stably control grain size and yield.

Methods

A collection of 327 wheat accessions from diverse sources were subjected to genome-wide association studies using genotyping-by-sequencing. Field trials were conducted under normal, delayed, and seriously delayed sowing conditions for grain length, width, and thousand-grain weight at two sites. Additionally, the additive main effects and multiplicative interaction (AMMI) model was applied to evaluate the stability of thousand-grain weight of 327 accessions across multiple sowing dates.

Results

Four wheat germplasm resources have been screened, demonstrating higher stability of thousand-grain weight. A total of 43, 35, and 39 significant MTAs were determined across all chromosomes except for 4D under the three sowing dates, respectively. A total of 10.31% of MTAs that stably affect wheat grain size could be repeatedly identified in at least two sowing dates, with PVE ranging from 0.03% to 38.06%. Among these, six were for GL, three for GW, and one for TGW. There were three novel and stable loci (4A_598189950, 4B_307707920, 2D_622241054) located in conserved regions of the genome, which provide excellent genetic resources for pyramid breeding strategies of superior loci. Our findings offer a theoretical basis for cultivar improvement and marker-assisted selection in wheat breeding practices.

Genome-wide identification of microsatellites for mapping, genetic diversity and cross-transferability in wheat (Triticum spp)

Wheat (Triticum aestivum L.) is a crucial global staple crop, and is consistently being improved to enhance yield, disease resistance, and quality traits. However, the development of molecular markers is a challenging task due to its hexaploid genome. Molecular marker system such as simple sequence repeat (SSR) and single nucleotide polymorphism (SNP) are helpful for breeding, but SNP has limitations due to its development cost and its conversion to breeder markers. The study proposed an in-silico approach, by utilizing the low-cost transcriptome sequencing of two parental lines, 'TAC 75' and 'WH 1105', to identify polymorphic SSRs for mapping in a recombinant inbred line (RIL) population. This study introduces a new approach to bridge wheat genetics intricacies and next-generation sequencing potential. It presents a comprehensive genome-wide SSR distribution using IWGSC CS RefSeq v2.1 genome assembly and to identify 189 polymorphic loci through in-silico strategy. Of these, 54.76% showed polymorphism between parents, surpassing the traditional low polymorphic success rate. A RIL population screening validated these markers, demonstrating the fitness of identified markers through chi-square tests. The designed SSRs were also validated for genetic diversity analysis in a subset of 37 Indian wheat genotypes and cross-transferability in the wild/relative wheat species. In diversity analysis, a subset of 38 markers revealed 95 alleles (2.5 allele/locus), indicating substantial genetic variation. Population structure analysis unveiled three distinct groups, supported by phylogenetic and PCoA analyses. Further the polymorphic SSRs were also analyzed for SSR-gene association using gene ontology analysis. By utilizing the developing seed transcriptome data within parental lines, the study has enhanced the polymorphic SSR identification precision and facilitated in the RIL population. The undertaken study pioneers the use of transcriptome sequencing and genetic mapping to overcome challenges posed by the intricate wheat genome. This approach offers a cost-effective, less labour-intensive alternative to conventional methods, providing a platform for advancing wheat breeding research.

Different wheat loci are associated to heritable free asparagine content in grain grown under different water and nitrogen availability

KEY MESSAGE

Different wheat QTLs were associated to the free asparagine content of grain grown in four different conditions. Environmental effects are a key factor when selecting for low acrylamide-forming potential. The amount of free asparagine in grain of a wheat genotype determines its potential to form harmful acrylamide in derivative food products. Here, we explored the variation in the free asparagine, aspartate, glutamine and glutamate contents of 485 accessions reflecting wheat worldwide diversity to define the genetic architecture governing the accumulation of these amino acids in grain. Accessions were grown under high and low nitrogen availability and in water-deficient and well-watered conditions, and plant and grain phenotypes were measured. Free amino acid contents of grain varied from 0.01 to 1.02 mg g-1 among genotypes in a highly heritable way that did not correlate strongly with grain yield, protein content, specific weight, thousand-kernel weight or heading date. Mean free asparagine content was 4% higher under high nitrogen and 3% higher in water-deficient conditions. After genotyping the accessions, single-locus and multi-locus genome-wide association study models were used to identify several QTLs for free asparagine content located on nine chromosomes. Each QTL was associated with a single amino acid and growing environment, and none of the QTLs colocalised with genes known to be involved in the corresponding amino acid metabolism. This suggests that free asparagine content is controlled by several loci with minor effects interacting with the environment. We conclude that breeding for reduced asparagine content is feasible, but should be firmly based on multi-environment field trials.

KEY MESSAGE

Different wheat QTLs were associated to the free asparagine content of grain grown in four different conditions. Environmental effects are a key factor when selecting for low acrylamide-forming potential.

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.

Characterization of a partially saturated and glycosylated apocarotenoid from wheat that is depleted upon leaf rust infection

Recent semi-targeted metabolomics studies have highlighted a number of metabolites in wheat that associate with leaf rust resistance genes and/or rust infection. Here, we report the structural characterization of a novel glycosylated and partially saturated apocarotenoid, reminiscent of a reduced form of mycorradicin, (6E,8E,10E)-4,9-dimethyl-12-oxo-12-((3,4,5-trihydroxy-6-(2-hydroxyethoxy)tetrahydro-2H-pyran-2-yl)methoxy)-3-((3,4,5-trihydroxy-6-(hydroxymethyl)tetrahydro-2H-pyran-2-yl)oxy)dodeca-6,8,10-trienoic acid, isolated from Triticum aestivum L. (Poaceae) variety 'Thatcher' (Tc) flag leaves. While its accumulation was not associated with any of Lr34, Lr67 or Lr22a resistance genes, infection of Tc with leaf rust was found to deplete it, consistent with the idea of this metabolite being a glycosylated-storage form of an apocarotenoid of possible relevance to plant defense. A comparative analysis of wheat transcriptomic changes shows modulation of terpenoid, carotenoid, UDP-glycosyltransferase and glycosylase -related gene expression profiles, consistent with anticipated biosynthesis and degradation mechanisms. However, details of the exact nature of the relevant pathways remain to be validated in the future. Together these findings highlight another example of the breadth of unique metabolites underlying plant host-fungal pathogen interactions.

Bread wheat satellitome: a complex scenario in a huge genome

In bread wheat (Triticum aestivum L.), chromosome associations during meiosis are extremely regulated and initiate at the telomeres and subtelomeres, which are enriched in satellite DNA (satDNA). We present the study and characterization of the bread wheat satellitome to shed light on the molecular organization of wheat subtelomeres. Our results revealed that the 2.53% of bread wheat genome is composed by satDNA and subtelomeres are particularly enriched in such DNA sequences. Thirty-four satellite DNA (21 for the first time in this work) have been identified, analyzed and cytogenetically validated. Many of the satDNAs were specifically found at particular subtelomeric chromosome regions revealing the asymmetry in subtelomere organisation among the wheat subgenomes, which might play a role in proper homologous recognition and pairing during meiosis. An integrated physical map of the wheat satellitome was also constructed. To the best of our knowledge, our results show that the combination of both cytogenetics and genome research allowed the first comprehensive analysis of the wheat satellitome, shedding light on the complex wheat genome organization, especially on the polymorphic nature of subtelomeres and their putative implication in chromosome recognition and pairing during meiosis.

Transcriptomic insights into shared responses to Fusarium crown rot infection and drought stresses in bread wheat (Triticum aestivum L.)

KEY MESSAGE

Shared changes in transcriptomes caused by Fusarium crown rot infection and drought stress were investigated based on a single pair of near-isogenic lines developed for a major locus conferring tolerance to both stresses. Fusarium crown rot (FCR) is a devastating disease in many areas of cereal production worldwide. It is well-known that drought stress enhances FCR severity but possible molecular relationship between these two stresses remains unclear. To investigate their relationships, we generated several pairs of near isogenic lines (NILs) targeting a locus conferring FCR resistance on chromosome 2D in bread wheat. One pair of these NILs showing significant differences between the two isolines for both FCR resistance and drought tolerance was used to investigate transcriptomic changes in responsive to these two stresses. Our results showed that the two isolines likely deployed different strategies in dealing with the stresses, and significant differences in expressed gene networks exist between the two time points of drought stresses evaluated in this study. Nevertheless, results from analysing Gene Ontology terms and transcription factors revealed that similar regulatory frameworks were activated in coping with these two stresses. Based on the position of the targeted locus, changes in expression following FCR infection and drought stresses, and the presence of non-synonymous variants between the two isolines, several candidate genes conferring resistance or tolerance to these two types of stresses were identified. The NILs generated, the large number of DEGs with single-nucleotide polymorphisms detected between the two isolines, and the candidate genes identified would be invaluable in fine mapping and cloning the gene(s) underlying the targeted locus.

A co-located QTL for seven spike architecture-related traits shows promising breeding use potential in common wheat (Triticum aestivum L.)

KEY MESSAGE

A co-located novel QTL for TFS, FPs, FMs, FFS, FFPs, KWS, and KWPs with potential of improving wheat yield was identified and validated. Spike-related traits, including fertile florets per spike (FFS), kernel weight per spike (KWS), total florets per spike (TFS), florets per spikelet (FPs), florets in the middle spikelet (FMs), fertile florets per spikelet (FFPs), and kernel weight per spikelet (KWPs), are key traits in improving wheat yield. In the present study, quantitative trait loci (QTL) for these traits evaluated under various environments were detected in a recombinant inbred line population (msf/Chuannong 16) mainly genotyped using the 16 K SNP array. Ultimately, we identified 60 QTL, but only QFFS.sau-MC-1A for FFS was a major and stably expressed QTL. It was located on chromosome arm 1AS, where loci for TFS, FPs, FMs, FFS, FFPs, KWS, and KWPs were also simultaneously co-mapped. The effect of QFFS.sau-MC-1A was further validated in three independent segregating populations using a Kompetitive Allele-Specific PCR marker. For the co-located QTL, QFFS.sau-MC-1A, the presence of a positive allele from msf was associate with increases for all traits: + 12.29% TFS, + 10.15% FPs, + 13.97% FMs, + 17.12% FFS, + 14.75% FFPs, + 22.17% KWS, and + 19.42% KWPs. Furthermore, pleiotropy analysis showed that the positive allele at QFFS.sau-MC-1A simultaneously increased the spike length, spikelet number per spike, and thousand-kernel weight. QFFS.sau-MC-1A represents a novel QTL for marker-assisted selection with the potential for improving wheat yield. Four genes, TraesCS1A03G0012700, TraesCS1A03G0015700, TraesCS1A03G0016000, and TraesCS1A03G0016300, which may affect spike development, were predicted in the physical interval harboring QFFS.sau-MC-1A. Our results will help in further fine mapping QFFS.sau-MC-1A and be useful for improving wheat yield.

Genetic analysis of late-maturity α-amylase in twelve wheat populations

MAIN CONCLUSION

Genetic loci, particularly those with an effect in the independent panel, could be utilised to further reduce LMA expression when used with favourable combinations of genes known to affect LMA. Late maturity α-amylase (LMA) is a grain quality defect involving elevated α-amylase within the aleurone of wheat (Triticum aestivum L.) grains. The genes known to affect expression are the reduced height genes Rht-B1 (chromosome 4B) and Rht-D1 (chromosome 4D), and an ent-copalyl diphosphate synthase gene (LMA-1) on chromosome 7B. Other minor effect loci have been reported, but these are poorly characterised and further genetic understanding is needed. In this study, twelve F4-derived populations were created through single seed descent, genotyped and evaluated for LMA. LMA-1 haplotype C and the Rht-D1b allele substantially reduced LMA expression. The alternative dwarfing genes Rht13 and Rht18 had no significant effect on LMA expression. Additional quantitative trait loci (QTL) were mapped at 16 positions in the wheat genome. Effects on LMA expression were detected for four of these QTL in a large independent panel of Australian wheat lines. The QTL detected in mapping populations and confirmed in the large independent panel provide further opportunity for selection against LMA, especially if combined with Rht-D1b and/or favourable haplotypes of LMA-1.

Enrichment and Diversification of the Wheat Genome via Alien Introgression

Wheat, including durum and common wheat, respectively, is an allopolyploid with two or three homoeologous subgenomes originating from diploid wild ancestral species. The wheat genome's polyploid origin consisting of just three diploid ancestors has constrained its genetic variation, which has bottlenecked improvement. However, wheat has a large number of relatives, including cultivated crop species (e.g., barley and rye), wild grass species, and ancestral species. Moreover, each ancestor and relative has many other related subspecies that have evolved to inhabit specific geographic areas. Cumulatively, they represent an invaluable source of genetic diversity and variation available to enrich and diversify the wheat genome. The ancestral species share one or more homologous genomes with wheat, which can be utilized in breeding efforts through typical meiotic homologous recombination. Additionally, genome introgressions of distant relatives can be moved into wheat using chromosome engineering-based approaches that feature induced meiotic homoeologous recombination. Recent advances in genomics have dramatically improved the efficacy and throughput of chromosome engineering for alien introgressions, which has served to boost the genetic potential of the wheat genome in breeding efforts. Here, we report research strategies and progress made using alien introgressions toward the enrichment and diversification of the wheat genome in the genomics era.

Heat stress during male meiosis impairs cytoskeletal organization, spindle assembly and tapetum degeneration in wheat

The significance of heat stress in agriculture is ever-increasing with the progress of global climate changes. Due to a negative effect on the yield of staple crops, including wheat, the impairment of plant reproductive development triggered by high ambient temperature became a restraint in food production. Although the heat sensitivity of male meiosis and the following gamete development in wheat has long been recognized, a detailed structural characterization combined with a comprehensive gene expression analysis has not been done about this phenomenon. We demonstrate here that heat stress severely alters the cytoskeletal configuration, triggers the failure of meiotic division in wheat. Moreover, it changes the expression of genes related to gamete development in male meiocytes and the tapetum layer in a genotype-dependent manner. 'Ellvis', a heat-tolerant winter wheat cultivar, showed high spikelet fertility rate and only scarce structural aberrations upon exposure to high temperature. In addition, heat shock genes and genes involved in scavenging reactive oxygen species were significantly upregulated in 'Ellvis', and the expression of meiosis-specific and major developmental genes showed high stability in this cultivar. In the heat-sensitive 'Mv 17-09', however, genes participating in cytoskeletal fiber nucleation, the spindle assembly checkpoint genes, and tapetum-specific developmental regulators were downregulated. These alterations may be related to the decreased cytoskeleton content, frequent micronuclei formation, and the erroneous persistence of the tapetum layer observed in the sensitive genotype. Our results suggest that understanding the heat-sensitive regulation of these gene functions would be an essential contribution to the development of new, heat-tolerant cultivars.

Dissection of a rapidly evolving wheat resistance gene cluster by long-read genome sequencing accelerated the cloning of Pm69

Gene cloning in repeat-rich polyploid genomes remains challenging. Here, we describe a strategy for overcoming major bottlenecks in cloning of the powdery mildew resistance gene (R-gene) Pm69 derived from tetraploid wild emmer wheat. A conventional positional cloning approach was not effective owing to suppressed recombination. Chromosome sorting was compromised by insufficient purity. A Pm69 physical map, constructed by assembling Oxford Nanopore Technology (ONT) long-read genome sequences, revealed a rapidly evolving nucleotide-binding leucine-rich repeat (NLR) R-gene cluster with structural variations. A single candidate NLR was identified by anchoring RNA sequencing reads from susceptible mutants to ONT contigs and was validated by virus-induced gene silencing. Pm69 is likely a newly evolved NLR and was discovered in only one location across the wild emmer wheat distribution range in Israel. Pm69 was successfully introgressed into cultivated wheat, and a diagnostic molecular marker was used to accelerate its deployment and pyramiding with other R-genes.

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.

Ameliorating the effects of multiple stresses on agronomic traits in crops: modern biotechnological and omics approaches

While global climate change poses a significant environmental threat to agriculture, the increasing population is another big challenge to food security. To address this, developing crop varieties with increased productivity and tolerance to biotic and abiotic stresses is crucial. Breeders must identify traits to ensure higher and consistent yields under inconsistent environmental challenges, possess resilience against emerging biotic and abiotic stresses and satisfy customer demands for safer and more nutritious meals. With the advent of omics-based technologies, molecular tools are now integrated with breeding to understand the molecular genetics of genotype-based traits and develop better climate-smart crops. The rapid development of omics technologies offers an opportunity to generate novel datasets for crop species. Identifying genes and pathways responsible for significant agronomic traits has been made possible by integrating omics data with genetic and phenotypic information. This paper discusses the importance and use of omics-based strategies, including genomics, transcriptomics, proteomics and phenomics, for agricultural and horticultural crop improvement, which aligns with developing better adaptability in these crop species to the changing climate conditions.

Genome-Wide Identification of the CBF Gene Family and ICE Transcription Factors in Walnuts and Expression Profiles under Cold Conditions

Cold stress impacts woody tree growth and perennial production, especially when the temperature rapidly changes in late spring. To address this issue, we conducted the genome-wide identification of two important transcription factors (TFs), CBF (C-repeat binding factors) and ICE (inducers of CBF expression), in three walnut (Juglans) genomes. Although the CBF and ICE gene families have been identified in many crops, very little systematic analysis of these genes has been carried out in J. regia and J. sigillata. In this study, we identified a total of 16 CBF and 12 ICE genes in three Juglans genomes using bioinformatics analysis. Both CBF and ICE had conserved domains, motifs, and gene structures, which suggests that these two TFs were evolutionarily conserved. Most ICE genes are located at both ends of the chromosomes. The promoter cis-regulatory elements of CBF and ICE genes are largely involved in light and phytohormone responses. Based on 36 RNA sequencing of leaves from four walnut cultivars ('Zijing', 'Lvling', 'Hongren', and 'Liao1') under three temperature conditions (8 °C, 22 °C, and 5 °C) conditions in late spring, we found that the ICE genes were expressed more highly than CBFs. Both CBF and ICE proteins interacted with cold-related proteins, and many putative miRNAs had interactions with these two TFs. These results determined that CBF1 and ICE1 play important roles in the tolerance of walnut leaves to rapid temperature changes. Our results provide a useful resource on the function of the CBF and ICE genes related to cold tolerance in walnuts.

A new wild emmer wheat panel allows to map new loci associated with resistance to stem rust at seedling stage

Wheat stem rust, caused by Puccinia graminis f. sp. tritici (Pgt), is a major wheat disease worldwide. A collection of 283 wild emmer wheat [Triticum turgidum L. subsp. dicoccoides (Körn. ex Asch. & Graebn.) Thell] accessions, representative of the entire Fertile Crescent region where wild emmer naturally occurs, was assembled, genotyped, and characterized for population structure, genetic diversity, and rate of linkage disequilibrium (LD) decay. Then, the collection was employed for mapping Pgt resistance genes, as a proof of concept of the effectiveness of genome-wide association studies in wild emmer. The collection was evaluated in controlled conditions for reaction to six common Pgt pathotypes (TPMKC, TTTTF, JRCQC, TRTTF, TTKSK/Ug99, and TKTTF). Most resistant accessions originated from the Southern Levant wild emmer lineage, with some showing a resistance reaction toward three to six tested races. Association analysis was conducted considering a 12K polymorphic single-nucleotide polymorphisms dataset, kinship relatedness between accessions, and population structure. Eleven significant marker-trait associations (MTA) were identified across the genome, which explained from 17% to up to 49% of phenotypic variance with an average 1.5 additive effect (based on the 1-9 scoring scale). The identified loci were either effective against single or multiple races. Some MTAs colocalized with known Pgt resistance genes, while others represent novel resistance loci useful for durum and bread wheat prebreeding. Candidate genes with an annotated function related to plant response to pathogens were identified at the regions linked to the resistance and defined according to the estimated small LD (about 126 kb), as typical of wild species.

Integrative physiological, critical plant endogenous hormones, and transcriptomic analyses reveal the difenoconazole stress response mechanism in wheat (Triticum aestivum L.)

Difenoconazole (DFN) is widely utilized as a fungicide in wheat production. However, its accumulation in plant tissues has a profound impact on the physiological functions of wheat plants, thus severely threatening wheat growth and even jeopardizing human health. This study aims to comprehensively analyze the dynamic dissipation patterns of DFN, along with an investigation into the physiological, hormonal, and transcriptomic responses of wheat seedlings exposed to DFN. The results demonstrated that exposure of wheat roots to DFN (10 mg/kg in soil) led to a significant accumulation of DFN in wheat plants, with the DFN content in roots being notably higher than that in leaves. Accumulating DFN triggered an increase in reactive oxygen species content, malonaldehyde content, and antioxidant enzyme activities, while concurrently inhibiting photosynthesis. Transcriptome analysis further revealed that the number of differentially expressed genes was greater in roots compared with leaves under DFN stress. Key genes in roots and leaves that exhibited a positive response to DFN-induced stress were identified through weighted gene co-expression network analysis. Metabolic pathway analysis indicated that these key genes mainly encode proteins involved in glutathione metabolism, plant hormone signaling, amino acid metabolism, and detoxification/defense pathways. Further results indicated that abscisic acid and salicylic acid play vital roles in the detoxification of leaf and root DFN, respectively. In brief, the abovementioned findings contribute to a deeper understanding of the detrimental effects of DFN on wheat seedlings, while shedding light on the molecular mechanisms underlying the responses of wheat root and leaves to DFN exposure.

Three-dimensional chromatin architecture in plants - General features and novelties

Research on the three-dimensional (3D) structure of the genome and its distribution within the nuclear space has made a big leap in the last two decades. Work in the animal field has led to significant advances in our general understanding on eukaryotic genome organization. This did not only bring along insights into how the 3D genome interacts with the epigenetic landscape and the transcriptional machinery but also how 3D genome architecture is relevant for fundamental developmental processes, such as cell differentiation. In parallel, the 3D organization of plant genomes have been extensively studied, which resulted in both congruent and novel findings, contributing to a more complete view on how eukaryotic genomes are organized in multiple dimensions. Plant genomes are remarkably diverse in size, composition, and ploidy. Furthermore, as intrinsically sessile organisms without the possibility to relocate to more favorable environments, plants have evolved an elaborate epigenetic repertoire to rapidly respond to environmental challenges. The diversity in genome organization and the complex epigenetic programs make plants ideal study subjects to acquire a better understanding on universal features and inherent constraints of genome organization. Furthermore, considering a wide range of species allows us to study the evolutionary crosstalk between the various levels of genome architecture. In this article, we aim at summarizing important findings on 3D genome architecture obtained in various plant species. These findings cover many aspects of 3D genome organization on a wide range of levels, from gene loops to topologically associated domains and to global 3D chromosome configurations. We present an overview on plant 3D genome organizational features that resemble those in animals and highlight facets that have only been observed in plants to date.

Skim exome capture genotyping in wheat

Next-generation sequencing (NGS) technology advancements continue to reduce the cost of high-throughput genome-wide genotyping for breeding and genetics research. Skim sequencing, which surveys the entire genome at low coverage, has become feasible for quantitative trait locus (QTL) mapping and genomic selection in various crops. However, the genome complexity of allopolyploid crops such as wheat (Triticum aestivum L.) still poses a significant challenge for genome-wide genotyping. Targeted sequencing of the protein-coding regions (i.e., exome) reduces sequencing costs compared to whole genome re-sequencing and can be used for marker discovery and genotyping. We developed a method called skim exome capture (SEC) that combines the strengths of these existing technologies and produces targeted genotyping data while decreasing the cost on a per-sample basis compared to traditional exome capture. Specifically, we fragmented genomic DNA using a tagmentation approach, then enriched those fragments for the low-copy genic portion of the genome using commercial wheat exome baits and multiplexed the sequencing at different levels to achieve desired coverage. We demonstrated that for a library of 48 samples, ∼7-8× target coverage was sufficient for high-quality variant detection. For higher multiplexing levels of 528 and 1056 samples per library, we achieved an average coverage of 0.76× and 0.32×, respectively. Combining these lower coverage SEC sequencing data with genotype imputation using a customized wheat practical haplotype graph database that we developed, we identified hundreds of thousands of high-quality genic variants across the genome. The SEC method can be used for high-resolution QTL mapping, genome-wide association studies, genomic selection, and other downstream applications.

Multi-trait genomic selection improves the prediction accuracy of end-use quality traits in hard winter wheat

Improvement of end-use quality remains one of the most important goals in hard winter wheat (HWW) breeding. Nevertheless, the evaluation of end-use quality traits is confined to later development generations owing to resource-intensive phenotyping. Genomic selection (GS) has shown promise in facilitating selection for end-use quality; however, lower prediction accuracy (PA) for complex traits remains a challenge in GS implementation. Multi-trait genomic prediction (MTGP) models can improve PA for complex traits by incorporating information on correlated secondary traits, but these models remain to be optimized in HWW. A set of advanced breeding lines from 2015 to 2021 were genotyped with 8725 single-nucleotide polymorphisms and was used to evaluate MTGP to predict various end-use quality traits that are otherwise difficult to phenotype in earlier generations. The MTGP model outperformed the ST model with up to a twofold increase in PA. For instance, PA was improved from 0.38 to 0.75 for bake absorption and from 0.32 to 0.52 for loaf volume. Further, we compared MTGP models by including different combinations of easy-to-score traits as covariates to predict end-use quality traits. Incorporation of simple traits, such as flour protein (FLRPRO) and sedimentation weight value (FLRSDS), substantially improved the PA of MT models. Thus, the rapid low-cost measurement of traits like FLRPRO and FLRSDS can facilitate the use of GP to predict mixograph and baking traits in earlier generations and provide breeders an opportunity for selection on end-use quality traits by culling inferior lines to increase selection accuracy and genetic gains.

TaSINA2B, interacting with TaSINA1D, positively regulates drought tolerance and root growth in wheat (Triticum aestivum L.)

Wheat (Triticum aestivum L.) is an important food crop mainly grown in arid and semiarid regions worldwide, whose productivity is severely limited by drought stress. Although various E3 ubiquitin (Ub) ligases regulate drought stress, only a few SINA-type E3 Ub ligases are known to participate in such responses. Herein, we identified and cloned 15 TaSINAs from wheat. The transcription level of TaSINA2B was highly induced by drought, osmotic and abscisic acid treatments. Two-type promoters of TaSINA2B were found in 192 wheat accessions; furthermore wheat accessions with promoter TaSINA2BII showed a considerably higher level of drought tolerance and gene expression levels than those characterizing accessions with promoter TaSINA2BI that was mainly caused by a 64 bp insertion in its promoter. Enhanced drought tolerance of TaSINA2B-overexpressing (OE) transgenic wheat lines was found to be associated with root growth promotion. Further, an interaction between TaSINA2B and TaSINA1D was detected through yeast two-hybrid and bimolecular fluorescence complementation assays. And TaSINA1D-OE transgenic wheat lines showed similar drought tolerance and root growth phenotype to those observed when TaSINA2B was overexpressed. Therefore, the variation of TaSINA2B promoter contributed to the drought stress regulation, while TaSINA2B, interacting with TaSINA1D, positively regulated drought tolerance by promoting root growth.

Wheat stripe rust resistance locus YR63 is a hot spot for evolution of defence genes - a pangenome discovery

BACKGROUND

Stripe rust, caused by Puccinia striiformis f. sp. tritici (Pst), poses a threat to global wheat production. Deployment of widely effective resistance genes underpins management of this ongoing threat. This study focused on the mapping of stripe rust resistance gene YR63 from a Portuguese hexaploid wheat landrace AUS27955 of the Watkins Collection.

RESULTS

YR63 exhibits resistance to a broad spectrum of Pst races from Australia, Africa, Asia, Europe, Middle East and South America. It was mapped to the short arm of chromosome 7B, between two single nucleotide polymorphic (SNP) markers sunCS_YR63 and sunCS_67, positioned at 0.8 and 3.7 Mb, respectively, in the Chinese Spring genome assembly v2.1. We characterised YR63 locus using an integrated approach engaging targeted genotyping-by-sequencing (tGBS), mutagenesis, resistance gene enrichment and sequencing (MutRenSeq), RNA sequencing (RNASeq) and comparative genomic analysis with tetraploid (Zavitan and Svevo) and hexaploid (Chinese Spring) wheat genome references and 10+ hexaploid wheat genomes. YR63 is positioned at a hot spot enriched with multiple nucleotide-binding and leucine rich repeat (NLR) and kinase domain encoding genes, known widely for defence against pests and diseases in plants and animals. Detection of YR63 within these gene clusters is not possible through short-read sequencing due to high homology between members. However, using the sequence of a NLR member we were successful in detecting a closely linked SNP marker for YR63 and validated on a panel of Australian bread wheat, durum and triticale cultivars.

CONCLUSIONS

This study highlights YR63 as a valuable source for resistance against Pst in Australia and elsewhere. The closely linked SNP marker will facilitate rapid introgression of YR63 into elite cultivars through marker-assisted selection. The bottleneck of this study reinforces the necessity for a long-read sequencing such as PacBio or Oxford Nanopore based techniques for accurate detection of the underlying resistance gene when it is part of a large gene cluster.

Subgenome phasing for complex allopolyploidy: case-based benchmarking and recommendations

Accurate subgenome phasing is crucial for understanding the origin, evolution and adaptive potential of polyploid genomes. SubPhaser and WGDI software are two common methodologies for subgenome phasing in allopolyploids, particularly in scenarios lacking known diploid progenitors. Triggered by a recent debate over the subgenomic origins of the cultivated octoploid strawberry, we examined four well-documented complex allopolyploidy cases as benchmarks, to evaluate and compare the accuracy of the two software. Our analysis demonstrates that the subgenomic structure phased by both software is in line with prior research, effectively tracing complex allopolyploid evolutionary trajectories despite the limitations of each software. Furthermore, using these validated methodologies, we revisited the controversial issue regarding the progenitors of the octoploid strawberry. The results of both methodologies reaffirm Fragaria vesca and Fragaria iinumae as progenitors of the octoploid strawberry. Finally, we propose recommendations for enhancing the accuracy of subgenome phasing in future studies, recognizing the potential of integrated tools for advanced complex allopolyploidy research and offering a new roadmap for robust subgenome-based phylogenetic analysis.

An online database for einkorn wheat to aid in gene discovery and functional genomics studies

Diploid A-genome wheat (einkorn wheat) presents a nutrition-rich option as an ancient grain crop and a resource for the improvement of bread wheat against abiotic and biotic stresses. Realizing the importance of this wheat species, reference-level assemblies of two einkorn wheat accessions were generated (wild and domesticated). This work reports an einkorn genome database that provides an interface to the cereals research community to perform comparative genomics, applied genetics and breeding research. It features queries for annotated genes, the use of a recent genome browser release, and the ability to search for sequence alignments using a modern BLAST interface. Other features include a comparison of reference einkorn assemblies with other wheat cultivars through genomic synteny visualization and an alignment visualization tool for BLAST results. Altogether, this resource will help wheat research and breeding. Database URL  https://wheat.pw.usda.gov/GG3/pangenome.

Genome wide association in Spanish bread wheat landraces identifies six key genomic regions that constitute potential targets for improving grain yield related traits

KEY MESSAGE

Association mapping conducted in 189 Spanish bread wheat landraces revealed six key genomic regions that constitute stable QTLs for yield and include 15 candidate genes. Genetically diverse landraces provide an ideal population to conduct association analysis. In this study, association mapping was conducted in a collection of 189 Spanish bread wheat landraces whose genomic diversity had been previously assessed. These genomic data were combined with characterization for yield-related traits, including grain size and shape, and phenological traits screened across five seasons. The association analysis revealed a total of 881 significant marker trait associations, involving 434 markers across the genome, that could be grouped in 366 QTLs based on linkage disequilibrium. After accounting for days to heading, we defined 33 high density QTL genomic regions associated to at least four traits. Considering the importance of detecting stable QTLs, 6 regions associated to several grain traits and thousand kernel weight in at least three environments were selected as the most promising ones to harbour targets for breeding. To dissect the genetic cause of the observed associations, we studied the function and in silico expression of the 413 genes located inside these six regions. This identified 15 candidate genes that provide a starting point for future analysis aimed at the identification and validation of wheat yield related genes.

Genetic Mapping of a Novel Gene PmAege7M from Aegilops geniculata Conferring Resistance to Wheat Powdery Mildew

Powdery mildew, caused by Blumeria graminis f. sp. tritici (Bgt), is one of the most damaging foliage diseases of wheat across the world. Aegilops geniculata Roth is a valuable gene resource for enhancing wheat resistance to powdery mildew. This study identified Ae. geniculata accession PI 487224 as immune and PI 487228 as susceptible to powdery mildew. Genetic analysis of the F1, F2, and F2:3 progeny derived from PI 487224 × PI 487228 showed that powdery mildew resistance in PI 487224 was controlled by two independent dominant genes located on two different nonhomologous chromosomes. By combing bulked segregant RNA-Seq, genetic linkage analysis of a single resistance gene segregation population, and marker analysis of a set of 14 wheat-Ae. geniculata chromosome addition lines, one of the resistance genes, temperately designated PmAege7M, was mapped to a 4.9-cM interval flanked by markers STS7-55926 and SNP7-45792/STS7-65911 on the long arm of chromosome 7 Mg of PI 487224, spanning 604.73 to 622.82 Mb on the 7D long arm based on the Ae. tauschii reference genome (Aet_v4.0). The map and closely linked markers of PmAege7M from Ae. geniculata in this study will facilitate the transfer of PmAege7M into common wheat and fine mapping of the gene.

MOMS: A pipeline for scaffolding using multi-optical maps

Here, we report a new multi-optical maps scaffolder (MOMS) aiming at utilizing complementary information among optical maps labelled by distinct enzymes. This pipeline was designed for data structure organization, scaffolding by path traversal, gap-filling and molecule reuse of optical maps. Our testing showed that this pipeline has uncapped enzyme tolerance in scaffolding. This means that there are no inbuilt limits as to the number of maps generated by different enzymes that can be utilized by MOMS. For the genome assembly of the human GM12878 cell line, MOMS significantly improved the contiguity and completeness with an up to 144-fold increase of scaffold N50 compared with initial assemblies. Benchmarking on the genomes of human and O. sativa showed that MOMS is more effective and robust compared with other optical-map-based scaffolders. We believe this pipeline will contribute to high-fidelity chromosome assembly and chromosome-level evolutionary analysis.

Discovery of isoflavone phytoalexins in wheat reveals an alternative route to isoflavonoid biosynthesis

Isoflavones are a group of phenolic compounds mostly restricted to plants of the legume family, where they mediate important interactions with plant-associated microbes, including in defense from pathogens and in nodulation. Their well-studied health promoting attributes have made them a prime target for metabolic engineering, both for bioproduction of isoflavones as high-value molecules, and in biofortification of food crops. A key gene in their biosynthesis, isoflavone synthase, was identified in legumes over two decades ago, but little is known about formation of isoflavones outside of this family. Here we identify a specialized wheat-specific isoflavone synthase, TaCYP71F53, which catalyzes a different reaction from the leguminous isoflavone synthases, thus revealing an alternative path to isoflavonoid biosynthesis and providing a non-transgenic route for engineering isoflavone production in wheat. TaCYP71F53 forms part of a biosynthetic gene cluster that produces a naringenin-derived O-methylated isoflavone, 5-hydroxy-2',4',7-trimethoxyisoflavone, triticein. Pathogen-induced production and in vitro antimicrobial activity of triticein suggest a defense-related role for this molecule in wheat. Genomic and metabolic analyses of wheat ancestral grasses further show that the triticein gene cluster was introduced into domesticated emmer wheat through natural hybridization ~9000 years ago, and encodes a pathogen-responsive metabolic pathway that is conserved in modern bread wheat varieties.

Chromosome-scale assembly of the wild wheat relative Aegilops umbellulata

Wild wheat relatives have been explored in plant breeding to increase the genetic diversity of bread wheat, one of the most important food crops. Aegilops umbellulata is a diploid U genome-containing grass species that serves as a genetic reservoir for wheat improvement. In this study, we report the construction of a chromosome-scale reference assembly of Ae. umbellulata accession TA1851 based on corrected PacBio HiFi reads and chromosome conformation capture. The total assembly size was 4.25 Gb with a contig N50 of 17.7 Mb. In total, 36,268 gene models were predicted. We benchmarked the performance of hifiasm and LJA, two of the most widely used assemblers using standard and corrected HiFi reads, revealing a positive effect of corrected input reads. Comparative genome analysis confirmed substantial chromosome rearrangements in Ae. umbellulata compared to bread wheat. In summary, the Ae. umbellulata assembly provides a resource for comparative genomics in Triticeae and for the discovery of agriculturally important genes.

Genomic characterization and gene bank curation of Aegilops: the wild relatives of wheat

Genetic diversity found in crop wild relatives is critical to preserve and utilize for crop improvement to achieve sustainable food production amid climate change and increased demand. We genetically characterized a large collection of 1,041 Aegilops accessions distributed among 23 different species using more than 45K single nucleotide polymorphisms identified by genotyping-by-sequencing. The Wheat Genetics Resource Center (WGRC) Aegilops germplasm collection was curated through the identification of misclassified and redundant accessions. There were 49 misclassified and 28 sets of redundant accessions within the four diploid species. The curated germplasm sets now have improved utility for genetic studies and wheat improvement. We constructed a phylogenetic tree and principal component analysis cluster for all Aegilops species together, giving one of the most comprehensive views of Aegilops. The Sitopsis section and the U genome Aegilops clade were further scrutinized with in-depth population analysis. The genetic relatedness among the pair of Aegilops species provided strong evidence for the species evolution, speciation, and diversification. We inferred genome symbols for two species Ae. neglecta and Ae. columnaris based on the sequence read mapping and the presence of segregating loci on the pertinent genomes as well as genetic clustering. The high genetic diversity observed among Aegilops species indicated that the genus could play an even greater role in providing the critical need for untapped genetic diversity for future wheat breeding and improvement. To fully characterize these Aegilops species, there is an urgent need to generate reference assemblies for these wild wheats, especially for the polyploid Aegilops.

A large-scale whole-exome sequencing mutant resource for functional genomics in wheat

Hexaploid wheat (Triticum aestivum), a major staple crop, has a remarkably large genome of ~14.4 Gb (containing 106 913 high-confidence [HC] and 159 840 low-confidence [LC] genes in the Chinese Spring v2.1 reference genome), which poses a major challenge for functional genomics studies. To overcome this hurdle, we performed whole-exome sequencing to generate a nearly saturated wheat mutant database containing 18 025 209 mutations induced by ethyl methanesulfonate (EMS), carbon (C)-ion beams, or γ-ray mutagenesis. This database contains an average of 47.1 mutations per kb in each gene-coding sequence: the potential functional mutations were predicted to cover 96.7% of HC genes and 70.5% of LC genes. Comparative analysis of mutations induced by EMS, γ-rays, or C-ion beam irradiation revealed that γ-ray and C-ion beam mutagenesis induced a more diverse array of variations than EMS, including large-fragment deletions, small insertions/deletions, and various non-synonymous single nucleotide polymorphisms. As a test case, we combined mutation analysis with phenotypic screening and rapidly mapped the candidate gene responsible for the phenotype of a yellow-green leaf mutant to a 2.8-Mb chromosomal region. Furthermore, a proof-of-concept reverse genetics study revealed that mutations in gibberellic acid biosynthesis and signalling genes could be associated with negative impacts on plant height. Finally, we built a publically available database of these mutations with the corresponding germplasm (seed stock) repository to facilitate advanced functional genomics studies in wheat for the broad plant research community.