Shifting the limits in wheat research and breeding using a fully annotated reference genome

Multiple roles of wheat ferritin genes during stress treatment and TaFER5D-1 as a positive regulator in response to drought and salt tolerance

Ferritin not only regulates the plant's iron content but also plays a significant role in the plant's development and resistance to oxidative damage. However, the role of the FER family in wheat has not been systematically elucidated. In this study, 39 FERs identified from wheat and its ancestral species were clustered into two subgroups, and gene members from the same group contain relatively conservative protein models. The structural analyses indicated that the gene members from the same group contained relatively conserved protein models. The cis-acting elements and expression patterns analysis suggested that TaFERs might play an important role combating to abiotic and biotic stresses. In the transcriptional analysis, the TaFER5D-1 gene was found to be significantly up-regulated under drought and salt stresses and was, therefore, selected to further explore the biological functions Moreover, the GFP expression assay revealed the subcellular localization of TaFER5D-1 proteins in the chloroplast, nucleus, membrane and cytoplasm. Over-expression of TaFER5D-1 in transgenic Arabidopsis lines conferred greater tolerance to drought and salt stress. According to the qRT-PCR data, TaFER5D-1 gene over-expression increased the expression of genes related to root development (Atsweet-17 and AtRSL4), iron storage (AtVIT1 and AtYSL1), and stress response (AtGolS1 and AtCOR47). So it is speculated that TaFER5D-1 could improve stress tolerance by promoting root growth, iron storage, and stress-response ability. Thus, the current study provides insight into the role of TaFER genes in wheat.

Identification and Molecular Mapping of the Gene YrFL in Wheat Cultivar Flanders with Durable Resistance to Stripe Rust

Stripe rust, caused by Puccinia striiformis f. sp. tritici (Pst), is one of the most severe diseases of wheat (Triticum aestivum L.) worldwide. Identification and characterization of resistance genes is advantageous to cultivating wheat varieties with durable resistance, which is the most economic and effective strategy to control stripe rust. Flanders, a common wheat cultivar released in France in 1986, confers effective resistance to stripe rust both at the seedling and adult plant stages. To elucidate the genetic basis of resistance in Flanders, F1, F2, and F2:3 generations derived from the cross Mingxian169 × Flanders were evaluated with the most prevalent Chinese Pst race CYR33 at the seedling stage. Inheritance analysis showed that the stripe rust resistance of Flanders was controlled by a single dominant gene, temporarily designated as YrFL. Bulked segregant analysis (BSA) combined with a wheat 660K single-nucleotide polymorphism (SNP) array indicated that polymorphic SNP markers were mainly located in the 0 to 150 Mb on wheat chromosome 5A. One hundred and eleven kompetitive allele-specific PCR (KASP) and 39 simple sequence repeat (SSR) markers on chromosome 5A were used to locate the YrFL. Linkage analysis mapped YrFL with 19 KASP and three SSR markers on wheat chromosome 5AS, and the genetic distances of the closest flanking markers AX108925494 and Xbarc56 to YrFL were 0.6 and 2.0 cM, respectively. Chromosome location, resistance characterization, and molecular marker positions indicated that YrFL is likely a novel stripe rust resistance gene on wheat chromosome 5AS and could be pyramided with other resistance genes to improve resistance in wheat breeding programs.

Genome-wide association mapping in elite winter wheat breeding for yield improvement

Increased grain yield (GY) is the primary breeding target of wheat breeders. We performed the genome-wide association study (GWAS) on 168 elite winter wheat lines from an ongoing breeding program to identify the main determinants of grain yield. Sequencing of Diversity Array Technology fragments (DArTseq) resulted in 19,350 single-nucleotide polymorphism (SNP) and presence-absence variation (PAV) markers. We identified 15 main genomic regions located in ten wheat chromosomes (1B, 2B, 2D, 3A, 3D, 5A, 5B, 6A, 6B, and 7B) that explained from 7.9 to 20.3% of the variation in grain yield and 13.3% of the yield stability. Loci identified in the reduced genepool are important for wheat improvement using marker-assisted selection. We found marker-trait associations between three genes involved in starch biosynthesis and grain yield. Two starch synthase genes (TraesCS2B03G1238800 and TraesCS2D03G1048800) and a sucrose synthase gene (TraesCS3D03G0024300) were found in regions of QGy.rut-2B.2, QGy.rut-2D.1, and QGy.rut-3D, respectively. These loci and other significantly associated SNP markers found in this study can be used for pyramiding favorable alleles in high-yielding varieties or to improve the accuracy of prediction in genomic selection.

Pangenome-based trajectories of intracellular gene transfers in Poaceae unveil high cumulation in Triticeae

Intracellular gene transfers (IGTs) between the nucleus and organelles, including plastids and mitochondria, constantly reshape the nuclear genome during evolution. Despite the substantial contribution of IGTs to genome variation, the dynamic trajectories of IGTs at the pangenomic level remain elusive. Here, we developed an approach, IGTminer, that maps the evolutionary trajectories of IGTs using collinearity and gene reannotation across multiple genome assemblies. We applied IGTminer to create a nuclear organellar gene (NOG) map across 67 genomes covering 15 Poaceae species, including important crops. The resulting NOGs were verified by experiments and sequencing data sets. Our analysis revealed that most NOGs were recently transferred and lineage specific and that Triticeae species tended to have more NOGs than other Poaceae species. Wheat (Triticum aestivum) had a higher retention rate of NOGs than maize (Zea mays) and rice (Oryza sativa), and the retained NOGs were likely involved in photosynthesis and translation pathways. Large numbers of NOG clusters were aggregated in hexaploid wheat during 2 rounds of polyploidization, contributing to the genetic diversity among modern wheat accessions. We implemented an interactive web server to facilitate the exploration of NOGs in Poaceae. In summary, this study provides resources and insights into the roles of IGTs in shaping interspecies and intraspecies genome variation and driving plant genome evolution.

Genetic basis of genome size variation of wheat

Research on various species has revealed a connection between genome size variation and the physiological and ecological characteristics of the species, suggesting that it could be a crucial factor influencing a species' adaptability to different environments. Wheat, being one of the world's three primary grains, holds significance in this regard. Investigating the genome size of wheat and analyzing the genetic factors contributing to its variation could offer valuable insights for enhancing wheat agronomic traits. This project has developed a conservative site ratio calculation approach to determine the size of the wheat genome. Additionally, it employs flow cytometry and k-mer distribution analysis to validate this method. Furthermore, the researchers use re-sequencing data to investigate the impact of environmental selection pressure and transposon dynamics on the variation in the size of the wheat genome. The findings from this study demonstrate a strong relationship between the size of the wheat genome and several environmental factors. These results serve as a valuable reference for understanding the development of variation in the size of the hetero-hexaploid wheat genome. Moreover, they contribute to advancing fundamental research on the genetic mechanisms underlying wheat characteristics. Additionally, the study paves the way for exploring new research directions in wheat breeding, which holds promise for future advancements in this field.

Effects of Allopolyploidization and Homoeologous Chromosomal Segment Exchange on Homoeolog Expression in a Synthetic Allotetraploid Wheat under Variable Environmental Conditions

Allopolyploidy through the combination of divergent genomes into a common nucleus at doubled dosage is known as a potent genetic and evolutionary force. As a macromutation, a striking feature of allopolyploidy in comparison with other mutational processes is that 'genome shock' can be evoked, thereby generating rapid and saltational biological consequences. A major manifestation of genome shock is genome-wide gene expression rewiring, which previously remained to be fully elucidated. Here, using a large set of RNAseq-based transcriptomic data of a synthetic allotetraploid wheat (genome AADD) and its parental species, we performed in-depth analyses of changes in the genome-wide gene expression under diverse environmental conditions at the subgenome (homoeolog) level and investigated the additional effects of homoeologous chromosomal segment exchanges (abbreviated HEs). We show that allopolyploidy caused large-scale changes in gene expression that were variable across the conditions and exacerbated by both stresses and HEs. Moreover, although both subgenomes (A and D) showed clear commonality in the changes, they responded differentially under variable conditions. The subgenome- and condition-dependent differentially expressed genes were enriched for different gene ontology terms implicating different biological functions. Our results provide new insights into the direct impacts of allopolyploidy on condition-dependent changes in subgenome expression and the additional effects of HEs in nascent allopolyploidy.

Breeding from 1891 to 2010 did not increase the content of amylase/trypsin-inhibitors in wheat (Triticum aestivum)

The prevalence of hypersensitivities towards wheat has increased in the last decades. Apart from celiac disease these include allergic and other inflammatory reactions summarized under the term non-celiac wheat sensitivity. One suspected trigger is the family of amylase/trypsin-inhibitors (ATIs), non-gluten proteins that are prominent wheat allergens and that activate the toll-like receptor 4 on intestinal immune cells to promote intestinal and extra-intestinal inflammation. We therefore quantified 13 ATIs in 60 German hexaploid winter wheat cultivars originating from 1891 to 2010 and harvested in three years by targeted liquid chromatography-tandem mass spectrometry combined with stable isotope dilution assay using specific marker peptides as internal standards. The total ATI content and that of the two major ATIs 0.19 and CM3 did not change from old cultivars (first registered from 1891 to 1950) to modern cultivars (1951-2010). There were also no significant changes in ATI distribution.

Fighting wheat powdery mildew: from genes to fields

KEY MESSAGE

Host resistance conferred by Pm genes provides an effective strategy to control powdery mildew. The study of Pm genes helps modern breeding develop toward more intelligent and customized. Powdery mildew of wheat is one of the most destructive diseases seriously threatening the crop yield and quality worldwide. The genetic research on powdery mildew (Pm) resistance has entered a new era. Many Pm genes from wheat and its wild and domesticated relatives have been mined and cloned. Meanwhile, modern breeding strategies based on high-throughput sequencing and genome editing are emerging and developing toward more intelligent and customized. This review highlights mining and cloning of Pm genes, molecular mechanism studies on the resistance and avirulence genes, and prospects for genomic-assisted breeding for powdery mildew resistance in wheat.

Identification of KASP markers and candidate genes for drought tolerance in wheat using 90K SNP array genotyping of near-isogenic lines targeting a 4BS quantitative trait locus

KEY MESSAGE

This study identified a novel SNP and developed a highly efficient KASP marker for drought tolerance in wheat by genotyping NILs targeting a major QTL for drought tolerance using an SNP array and validation with commercial varieties. Common wheat (Triticum aestivum L.) is an important winter crop worldwide and a typical allopolyploid with a large and complex genome. With global warming, the environmental volatility and incidence of drought in wheat-producing areas will increase. Molecular markers for drought tolerance are urgently needed to enhance drought tolerance breeding. Here, we genotyped four near-isogenic line (NIL) pairs targeting a major QTL qDSI.4B.1 on wheat chromosome arm 4BS for drought tolerance using the 90K SNP Illumina iSelect array and discovered a single nucleotide polymorphism (SNP) (Excalibur_c100336_106) with consistent genotype-phenotype associations among all four NIL pairs and their parents. Then, we converted the SNP into a Kompetitive Allele-Specific PCR (KASP) marker, with an accuracy of 100% for the four NIL pairs and their parents and as high as 81.8% for the 44 tested wheat lines with known phenotypes collected from Australia and China. Two genes near this SNP were suggested as candidate genes for drought tolerance in wheat after checking the Chinese Spring reference genome annotation version 1.1. One gene, TraesCS4B02G085300, encodes an F-box protein reportedly related to the ABA network, a main pathway for drought tolerance, and another gene, TraesCS4B02G085400, encodes a calcineurin-like metallophos-phoesterase transmembrane protein, which participates in Ca2+-dependent phosphorylation regulatory system. Based on this work and previous research on pre-harvest sprouting, we established a quick and efficient general SQV-based approach for KASP marker development, integrating genotyping by SNP arrays (S) using NILs targeting major QTL for a specific trait (Q) and validating them with commercial varieties (V). The identified SNP and developed KASP marker could be applied to marker-assisted selection in drought breeding, and further study of the candidate genes may improve our understanding of drought tolerance in wheat.

Genetic analysis of stripe rust resistance in the common wheat line Kfa/2*Kachu under a Chinese rust environment

KEY MESSAGE

We mapped a new race-specific seedling stripe rust resistance gene on wheat chromosome 5BL and a new APR locus QYr.hazu-2BS from CIMMYT wheat line Kfa/2*Kachu. Breeding resistant wheat (Triticum aestivum) varieties is the most economical and efficient way to manage wheat stripe rust, but requires the prior identification of new resistance genes and development of associated molecular markers for marker-assisted selection. To map stripe rust resistance loci in wheat, we used a recombinant inbred line population generated by crossing the stripe rust-resistant parent 'Kfa/2*Kachu' and the susceptible parent 'Apav#1'. We employed genotyping-by-sequencing and bulked segregant RNA sequencing to map a new race-specific seedling stripe rust resistance gene, which we designated YrK, to wheat chromosome arm 5BL. TraesCS5B02G330700 encodes a receptor-like kinase and is a high-confidence candidate gene for YrK based on virus-induced gene silencing results and the significant induction of its expression 24 h after inoculation with wheat stripe rust. To assist breeding, we developed functional molecular markers based on the polymorphic single nucleotide polymorphisms in the coding sequence region of YrK. We also mapped four adult plant resistance (APR) loci to wheat chromosome arms 1BL, 2AS, 2BS and 4AL. Among these APR loci, we determined that QYr.hazu-1BL and QYr.hazu-2AS are allelic to the known pleiotropic resistance gene Lr46/Yr29/Pm39 and the race-specific gene Yr17, respectively. However, QYr.hazu-2BS is likely a new APR locus, for which we converted closely linked SNP polymorphisms into breeder-friendly Kompetitive allele-specific PCR (KASP) markers. In the present study, we provided new stripe rust resistance locus/gene and molecular markers for wheat breeder to develop rust-resistant wheat variety.

DMC1 stabilizes crossovers at high and low temperatures during wheat meiosis

Effective chromosome synapsis and crossover formation during meiosis are essential for fertility, especially in grain crops such as wheat. These processes function most efficiently in wheat at temperatures between 17-23 °C, although the genetic mechanisms for such temperature dependence are unknown. In a previously identified mutant of the hexaploid wheat reference variety 'Chinese Spring' lacking the long arm of chromosome 5D, exposure to low temperatures during meiosis resulted in asynapsis and crossover failure. In a second mutant (ttmei1), containing a 4 Mb deletion in chromosome 5DL, exposure to 13 °C led to similarly high levels of asynapsis and univalence. Moreover, exposure to 30 °C led to a significant, but less extreme effect on crossovers. Previously, we proposed that, of 41 genes deleted in this 4 Mb region, the major meiotic gene TaDMC1-D1 was the most likely candidate for preservation of synapsis and crossovers at low (and possibly high) temperatures. In the current study, using RNA-guided Cas9, we developed a new Chinese Spring CRISPR mutant, containing a 39 bp deletion in the 5D copy of DMC1, representing the first reported CRISPR-Cas9 targeted mutagenesis in Chinese Spring, and the first CRISPR mutant for DMC1 in wheat. In controlled environment experiments, wild-type Chinese Spring, CRISPR dmc1-D1 and backcrossed ttmei1 mutants were exposed to either high or low temperatures during the temperature-sensitive period from premeiotic interphase to early meiosis I. After 6-7 days at 13 °C, crossovers decreased by over 95% in the dmc1-D1 mutants, when compared with wild-type plants grown under the same conditions. After 24 hours at 30 °C, dmc1-D1 mutants exhibited a reduced number of crossovers and increased univalence, although these differences were less marked than at 13 °C. Similar results were obtained for ttmei1 mutants, although their scores were more variable, possibly reflecting higher levels of background mutation. These experiments confirm our previous hypothesis that DMC1-D1 is responsible for preservation of normal crossover formation at low and, to a certain extent, high temperatures. Given that reductions in crossovers have significant effects on grain yield, these results have important implications for wheat breeding, particularly in the face of climate change.

Plant genome resequencing and population genomics: Current status and future prospects

Advances in DNA sequencing technology have sparked a genomics revolution, driving breakthroughs in plant genetics and crop breeding. Recently, the focus has shifted from cataloging genetic diversity in plants to exploring their functional significance and delivering beneficial alleles for crop improvement. This transformation has been facilitated by the increasing adoption of whole-genome resequencing. In this review, we summarize the current progress of population-based genome resequencing studies and how these studies affect crop breeding. A total of 187 land plants from 163 countries have been resequenced, comprising 54 413 accessions. As part of resequencing efforts 367 traits have been surveyed and 86 genome-wide association studies have been conducted. Economically important crops, particularly cereals, vegetables, and legumes, have dominated the resequencing efforts, leaving a gap in 49 orders, including Lycopodiales, Liliales, Acorales, Austrobaileyales, and Commelinales. The resequenced germplasm is distributed across diverse geographic locations, providing a global perspective on plant genomics. We highlight genes that have been selected during domestication, or associated with agronomic traits, and form a repository of candidate genes for future research and application. Despite the opportunities for cross-species comparative genomics, many population genomic datasets are not accessible, impeding secondary analyses. We call for a more open and collaborative approach to population genomics that promotes data sharing and encourages contribution-based credit policy. The number of plant genome resequencing studies will continue to rise with the decreasing DNA sequencing costs, coupled with advances in analysis and computational technologies. This expansion, in terms of both scale and quality, holds promise for deeper insights into plant trait genetics and breeding design.

Chromosome-specific painting in Thinopyrum species using bulked oligonucleotides

Chromosome-specific painting probes were developed to identify the individual chromosomes from 1 to 7E in Thinopyrum species and detect alien genetic material of the E genome in a wheat background. The E genome of Thinopyrum is closely related to the ABD genome of wheat (Triticum aestivum L.) and harbors genes conferring beneficial traits to wheat, including high yield, disease resistance, and unique end-use quality. Species of Thinopyrum vary from diploid (2n = 2x = 14) to decaploid (2n = 10x = 70), and chromosome structural variation and differentiation have arisen during polyploidization. To investigate the variation and evolution of the E genome, we developed a complete set of E genome-specific painting probes for identification of the individual chromosomes 1E to 7E based on the genome sequences of Th. elongatum (Host) D. R. Dewey and wheat. By using these new probes in oligonucleotide-based chromosome painting, we showed that Th. bessarabicum (PI 531711, EbEb) has a close genetic relationship with diploid Th. elongatum (EeEe), with five chromosomes (1E, 2E, 3E, 6E, and 7E) maintaining complete synteny in the two species except for a reciprocal translocation between 4 and 5Eb. All 14 pairs of chromosomes of tetraploid Th. elongatum have maintained complete synteny with those of diploid Th. elongatum (Thy14), but the two sets of E genomes have diverged. This study also demonstrated that the E genome-specific painting probes are useful for rapid and effective detection of the alien genetic material of E genome in wheat-Thinopyrum derived lines.

Centromere Plasticity With Evolutionary Conservation and Divergence Uncovered by Wheat 10+ Genomes

Centromeres (CEN) are the chromosomal regions that play a crucial role in maintaining genomic stability. The underlying highly repetitive DNA sequences can evolve quickly in most eukaryotes, and promote karyotype evolution. Despite their variability, it is not fully understood how these widely variable sequences ensure the homeostasis of centromere function. In this study, we investigated the genetics and epigenetics of CEN in a population of wheat lines from global breeding programs. We captured a high degree of sequences, positioning, and epigenetic variations in the large and complex wheat CEN. We found that most CENH3-associated repeats are Cereba element of retrotransposons and exhibit phylogenetic homogenization across different wheat lines, but the less-associated repeat sequences diverge on their own way in each wheat line, implying specific mechanisms for selecting certain repeat types as functional core CEN. Furthermore, we observed that CENH3 nucleosome structures display looser wrapping of DNA termini on complex centromeric repeats, including the repositioned CEN. We also found that strict CENH3 nucleosome positioning and intrinsic DNA features play a role in determining centromere identity among different lines. Specific non-B form DNAs were substantially associated with CENH3 nucleosomes for the repositioned centromeres. These findings suggest that multiple mechanisms were involved in the adaptation of CENH3 nucleosomes that can stabilize CEN. Ultimately, we proposed a remarkable epigenetic plasticity of centromere chromatin within the diverse genomic context, and the high robustness is crucial for maintaining centromere function and genome stability in wheat 10+ lines as a result of past breeding selections.

QTL Pyramiding Provides Marginal Improvement in 2NvS-Based Wheat Blast Resistance

Wheat blast, caused by the fungus Magnaporthe oryzae Triticum pathotype (MoT), is a devastating disease affecting South America, Bangladesh, and Zambia. Resistance to wheat blast has strongly relied on the 2NvS translocation; however, newer MoT isolates have increased aggressiveness, threatening the 2NvS translocation's effectiveness and durability. To identify genomic regions associated with wheat blast resistance, we performed a quantitative trait loci (QTL) mapping study using 187 double-haploid (DH) lines from a cross between the Brazilian wheat cultivars 'TBIO Alvorada' and 'TBIO Sossego', which are moderately resistant and susceptible to blast, respectively. The DH population was evaluated in a greenhouse in Brazil and Bolivia, and field conditions in Bolivia. Contrasting models best explained the relationship between traits evaluated according to differences in disease levels and the presence of the 2NvS. A large effect-locus, derived from 'TBIO Sossego', was identified on chromosome 2AS, which was confirmed to be 2NvS translocation and explained 33.5 to 82.4% of the phenotypic variance. Additional significant loci were identified on 5AL, 1DS, 4DS, 5DL, and 6DL chromosome arms with phenotypic variance <6%, but they were not consistent across trait-environment combinations. QTL pyramiding analyses showed that some specific loci had an additive effect when combined with the 2NvS, suggesting that stacking multiple loci may be an effective strategy to help manage wheat blast. The markers associated with the 2NvS can be used as dominant diagnostic markers for this alien translocation. Additional characterization of these loci using a broader set of MoT isolates is critical to validate their effectiveness against current MoT populations.

Mapping of Major Resistance Gene Qse.xn-2BL for Sharp Eyespot in the Wheat-Psathyrostachys huashanica Introgression Line H83

Sharp eyespot, a soil-borne disease of wheat (Triticum aestivum L.), is one of the most devastating diseases and severely affects grain production. The most efficient and economical method of controlling the disease is the utilization of genetic resistance. In this study, the wheat-Psathyrostachys huashanica introgression line H83 processed the enhanced resistance to Rhizoctonia cerealis isolate R0301 than its wheat parent 7182. A resistance locus in the 600 to 800 Mb interval of chromosome 2BL was screened using 244 segregation population F2 plants of H83×Huixianhong with bulked segregant analysis and wheat axiom 660K genotyping array. Furthermore, by using 12 kompetitive allele-specific PCR markers, a major resistance gene, designated as Qse.xn-2BL, was identified in a secondary segregating population with 138 F3:4 lines and initially mapped to a 765.6 to 775.5 Mb interval on chromosome 2BL. Molecular cytology analysis revealed that H83 probably has an alien introgression at the distal of chromosome 2BL, where it overlapped with the mapping target gene. Above all, H83 showed great potential to improve wheat resistance to sharp eyespot and can be expected to improve resistance in wheat breeding.

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

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

Supplementary Information

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

Einkorn genomics sheds light on history of the oldest domesticated wheat

Einkorn (Triticum monococcum) was the first domesticated wheat species, and was central to the birth of agriculture and the Neolithic Revolution in the Fertile Crescent around 10,000 years ago1,2. Here we generate and analyse 5.2-Gb genome assemblies for wild and domesticated einkorn, including completely assembled centromeres. Einkorn centromeres are highly dynamic, showing evidence of ancient and recent centromere shifts caused by structural rearrangements. Whole-genome sequencing analysis of a diversity panel uncovered the population structure and evolutionary history of einkorn, revealing complex patterns of hybridizations and introgressions after the dispersal of domesticated einkorn from the Fertile Crescent. We also show that around 1% of the modern bread wheat (Triticum aestivum) A subgenome originates from einkorn. These resources and findings highlight the history of einkorn evolution and provide a basis to accelerate the genomics-assisted improvement of einkorn and bread wheat.

Mapping QTL conferring flag leaf senescence in durum wheat cultivars

Flag leaf senescence is a critical factor affecting the yield and quality of wheat. The aim of this study was to identify QTLs associated with flag leaf senescence in an F10 recombinant inbred line population derived from durum wheats UC1113 and Kofa. Bulked segregant analysis using the wheat 660K SNP array identified 3225 SNPs between extreme-phenotype bulks, and the differential SNPs were mainly clustered on chromosomes 1A, 1B, 3B, 5A, 5B, and 7A. BSR-Seq indicated that the significant SNPs were mainly located in two intervals of 354.0-389.0 Mb and 8.0-15.0 Mb on 1B and 3B, respectively. Based on the distribution of significant SNPs on chromosomes 1B and 3B, a total of 109 insertion/deletion (InDel) markers were developed, and 8 of them were finally used to map QTL in UC1113/Kofa population for flag leaf senescence. Inclusive composite interval mapping identified two major QTL in marker intervals Mar2005-Mar2116 and Mar207-Mar289, explaining 14.2-15.4% and 31.4-68.6% of the phenotypic variances across environments, respectively. Using BSR-Seq, gene expression and sequence analysis, the TraesCS1B02G211600 and TraesCS3B02G023000 were identified as candidate senescence-associated genes. This study has potential to be used in cloning key genes for flag leaf senescence and provides available molecular markers for genotyping and marker-assisted selection breeding.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-023-01410-3.

Identification and characterization of RuvBL DNA helicase genes for tolerance against abiotic stresses in bread wheat (Triticum aestivum L.) and related species

Recombination UVB (sensitivity) like (RuvBL) helicase genes represent a conserved family of genes, which are known to be involved in providing tolerance against abiotic stresses like heat and drought. We identified nine wheat RuvBL genes, one each on nine different chromosomes, belonging to homoeologous groups 2, 3, and 4. The lengths of genes ranged from 1647 to 2197 bp and exhibited synteny with corresponding genes in related species including Ae. tauschii, Z. mays, O. sativa, H. vulgare, and B. distachyon. The gene sequences were associated with regulatory cis-elements and transposable elements. Two genes, namely TaRuvBL1a-4A and TaRuvBL1a-4B, also carried targets for a widely known miRNA, tae-miR164. Gene ontology revealed that these genes were closely associated with ATP-dependent formation of histone acetyltransferase complex. Analysis of the structure and function of RuvBL proteins revealed that the proteins were localized mainly in the cytoplasm. A representative gene, namely TaRuvBL1a-4A, was also shown to be involved in protein-protein interactions with ten other proteins. On the basis of phylogeny, RuvBL proteins were placed in two sub-divisions, namely RuvBL1 and RuvBL2, which were further classified into clusters and sub-clusters. In silico studies suggested that these genes were differentially expressed under heat/drought. The qRT-PCR analysis confirmed that expression of TaRuvBL genes differed among wheat cultivars, which differed in the level of thermotolerance. The present study advances our understanding of the biological role of wheat RuvBL genes and should help in planning future studies on RuvBL genes in wheat including use of RuvBL genes in breeding thermotolerant wheat cultivars.

Triticeae crop genome biology: an endless frontier

Triticeae, the wheatgrass tribe, includes several major cereal crops and their wild relatives. Major crops within the Triticeae are wheat, barley, rye, and oat, which are important for human consumption, animal feed, and rangeland protection. Species within this tribe are known for their large genomes and complex genetic histories. Powered by recent advances in sequencing technology, researchers worldwide have made progress in elucidating the genomes of Triticeae crops. In addition to assemblies of high-quality reference genomes, pan-genome studies have just started to capture the genomic diversities of these species, shedding light on our understanding of the genetic basis of domestication and environmental adaptation of Triticeae crops. In this review, we focus on recent signs of progress in genome sequencing, pan-genome analyses, and resequencing analysis of Triticeae crops. We also propose future research avenues in Triticeae crop genomes, including identifying genome structure variations, the association of genomic regions with desired traits, mining functions of the non-coding area, introgression of high-quality genes from wild Triticeae resources, genome editing, and integration of genomic resources.

Comparative genomic analysis of 5Mg chromosome of Aegilops geniculata and 5Uu chromosome of Aegilops umbellulata reveal genic diversity in the tertiary gene pool

Wheat is one of the most important cereal crops for the global food security. Due to its narrow genetic base, modern bread wheat cultivars face challenges from increasing abiotic and biotic stresses. Since genetic improvement is the most sustainable approach, finding novel genes and alleles is critical for enhancing the genetic diversity of wheat. The tertiary gene pool of wheat is considered a gold mine for genetic diversity as novel genes and alleles can be identified and transferred to wheat cultivars. Aegilops geniculata and Ae. umbellulata are the key members of the tertiary gene pool of wheat and harbor important genes against abiotic and biotic stresses. Homoeologous-group five chromosomes (5U^u and 5M^g) have been extensively studied from Ae. geniculata and Ae. umbellulata as they harbor several important genes including Lr57, Lr76, Yr40, Yr70, Sr53 and chromosomal pairing loci. In the present study, using chromosome DNA sequencing and RNAseq datasets, we performed comparative analysis to study homoeologous gene evolution in 5M^g, 5U^u, and group 5 wheat chromosomes. Our findings highlight the diversity of transcription factors and resistance genes, resulting from the differential expansion of the gene families. Both the chromosomes were found to be enriched with the "response to stimulus" category of genes providing resistance against biotic and abiotic stress. Phylogenetic study positioned the M genome closer to the D genome, with higher proximity to the A genome than the B genome. Over 4000 genes were impacted by SNPs on 5D, with 4-5% of those genes displaying non-disruptive variations that affect gene function.

Fine mapping and characterization of a major QTL for plant height on chromosome 5A in wheat

KEY MESSAGE

We precisely mapped QPH.caas-5AL for plant height in wheat, predicted candidate genes and confirmed genetic effects in a panel of wheat cultivars. Plant height is an important agronomic trait, and appropriately reduced height can improve yield potential and stability in wheat, usually combined with sufficient water and fertilizer. We previously detected a stable major-effect quantitative trait locus QPH.caas-5AL for plant height on chromosome 5A in a recombinant inbred line population of the cross 'Doumai × Shi 4185' using the wheat 90 K SNP assay. Here , QPH.caas-5AL was confirmed using new phenotypic data in additional environment and new-developed markers. We identified nine heterozygous recombinant plants for fine mapping of QPH.caas-5AL and developed 14 breeder-friendly kompetitive allele-specific PCR markers in the region of QPH.caas-5AL based on the genome re-sequencing data of parents. Phenotyping and genotyping analyses of secondary populations derived from the self-pollinated heterozygous recombinant plants delimited QPH.caas-5AL into an approximate 3.0 Mb physical region (521.0-524.0 Mb) according to the Chinese Spring reference genome. This region contains 45 annotated genes, and six of them were predicted as the candidates of QPH.caas-5AL based on genome and transcriptome sequencing analyses. We further validated that QPH.caas-5AL has significant effects on plant height but not yield component traits in a diverse panel of wheat cultivars; its dwarfing allele is frequently used in modern wheat cultivars. These findings lay a solid foundation for the map-based cloning of QPH.caas-5AL and also provide a breeding-applicable tool for its marker-assisted selection. Keymessage We precisely mapped QPH.caas-5AL for plant height in wheat, predicted candidate genes and confirmed genetic effects in a panel of wheat cultivars.

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

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

Genome wide association and haplotype analyses for the crease depth trait in bread wheat (Triticum aestivum L.)

Wheat grain has a complex structure that includes a crease on one side, and tissues within the crease region play an important role in nutrient transportation during wheat grain development. However, the genetic architecture of the crease region is still unclear. In this study, 413 global wheat accessions were resequenced and a method was developed for evaluating the phenotypic data of crease depth (CD). The CD values exhibited continuous and considerable large variation in the population, and the broad-sense heritability was 84.09%. CD was found to be positively correlated with grain-related traits and negatively with quality-related traits. Analysis of differentiation of traits between landraces and cultivars revealed that grain-related traits and CD were simultaneously improved during breeding improvement. Moreover, 2,150.8-Mb genetic segments were identified to fall within the selective sweeps between the landraces and cultivars; they contained some known functional genes for quality- and grain-related traits. Genome-wide association study (GWAS) was performed using around 10 million SNPs generated by genome resequencing and 551 significant SNPs and 18 QTLs were detected significantly associated with CD. Combined with cluster analysis of gene expression, haplotype analysis, and annotated information of candidate genes, two promising genes TraesCS3D02G197700 and TraesCS5A02G292900 were identified to potentially regulate CD. To the best of our knowledge, this is the first study to provide the genetic basis of CD, and the genetic loci identified in this study may ultimately assist in wheat breeding programs.

Genetic dissection of grain morphology and yield components in a wheat line with defective grain filling

KEY MESSAGE

We identified stable QTL for grain morphology and yield component traits in a wheat defective grain filling line and validated genetic effects in a panel of cultivars using breeding-relevant markers. Grain filling capacity is essential for grain yield and appearance quality in cereal crops. Identification of genetic loci for grain filling is important for wheat improvement. However, there are few genetic studies on grain filling in wheat. Here, a defective grain filling (DGF) line wdgf1 characterized by shrunken grains was identified in a population derived from multi-round crosses involving nine parents and a recombinant inbreed line (RIL) population was generated from the cross between wdgf1 and a sister line with normal grains. We constructed a genetic map of the RIL population using the wheat 15K single nucleotide polymorphism chip and detected 25 stable quantitative trait loci (QTL) for grain morphology and yield components, including three for DGF, eleven for grain size, six for thousand grain weight, three for grain number per spike and two for spike number per m². Among them, QDGF.caas-7A is co-located with QTGW.caas-7A and can explain 39.4-64.6% of the phenotypic variances, indicating that this QTL is a major locus controlling DGF. Sequencing and linkage mapping showed that TaSus2-2B and Rht-B1 were candidate genes for QTGW.caas-2B and the QTL cluster (QTGW.caas-4B, QGNS.caas-4B, and QSN.caas-4B), respectively. We developed kompetitive allele-specific PCR markers tightly linked to the stable QTL without corresponding to known yield-related genes, and validated their genetic effects in a diverse panel of wheat cultivars. These findings not only lay a solid foundation for genetic dissection underlying grain filling and yield formation, but also provide useful tools for marker-assisted breeding.

Genetic basis of resistance against powdery mildew in the wheat cultivar "Tabasco"

European winter wheat cultivar "Tabasco" was reported to have resistance to powdery mildew disease caused by Blumeria graminis f. sp. tritici (Bgt) in China. In previous studies, Tabasco was reported to have the resistance gene designated as Pm48 on the short arm of chromosome 5D when a mapping population was phenotyped with pathogen isolate Bgt19 collected in China and was genotyped with simple sequence repeat (SSR) markers. In this study, single-nucleotide polymorphism (SNP) chips were used to rapidly determine the resistance gene by mapping a new F2 population that was developed from Tabasco and a susceptible cultivar "Ningmaizi119" and inoculated with pathogen isolate NCF-D-1-1 that was collected in the USA. The segregation of resistance in the population was found to link with Pm2 which was identified in Tabasco. Therefore, it was concluded that the previously reported Pm48 on chromosome arm 5DS in Tabasco should be the Pm2 gene on the same chromosome. The Pm2 was also found in European cultivars "Mattis" and "Claire" but not in any of the accessions from diploid wheat Aegilops tauschii or modern cultivars such as "Gallagher," "Smith's Gold," and "OK Corral" being used in the Great Plains in the USA. A KASP marker was developed to track the resistance allele Pm2 in wheat breeding.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-023-01402-3.

Resistance to Powdery Mildew Is Conferred by Different Genetic Loci at the Adult-Plant and Seedling Stages in Winter Wheat Line Tianmin 668

Winter wheat line Tianmin 668 was crossed with susceptible cultivar Jingshuang 16 to develop 216 recombinant inbred lines (RILs) for dissecting its adult-plant resistance (APR) and all-stage resistance (ASR) against powdery mildew. The RIL population was genotyped on a 16K genotyping by target sequencing single-nucleotide polymorphism array and phenotyped in six field trials and in the greenhouse. Three loci-QPmtj.caas-2BL, QPmtj.caas-2AS, and QPmtj.caas-5AL-conferring APR to powdery mildew were detected on chromosomes 2BL, 2AS, and 5AL, respectively, of Tianmin 668. The effect of resistance to powdery mildew for QPmtj.caas-2BL was greater than that of the other two loci. A Kompetitive allele-specific PCR marker specific for QPmtj.caas-2BL was developed and verified on 402 wheat cultivars or breeding lines. Results of virulence and avirulence patterns to 17 Blumeria graminis f. sp. tritici isolates, bulked segregant analysis-RNA-sequencing, and a genetic linkage mapping identified a resistance allele at locus Pm4 in Tianmin 668 based on the seedling phenotypes of the RIL population. The PCR-based DNA sequence alignment and cosegregation of the functional marker with the phenotypes of the RIL population demonstrated that Pm4d was responsible for the ASR to isolate Bgt1 in Tianmin 668. The dissection of genetic loci for APR and ASR may facilitate the application of Tianmin 668 in developing powdery mildew-resistant wheat cultivars.

A major yellow rust resistance QTL on chromosome 6A shows increased frequency in recent Norwegian spring wheat cultivars and breeding lines

KEY MESSAGE

A major yellow rust resistance QTL, QYr.nmbu.6A, contributed consistent adult plant resistance in field trials across Europe, China, Kenya and Mexico. Puccinia striiformis f. sp. tritici, causing wheat yellow rust (YR), is one of the most devastating biotrophic pathogens affecting global wheat yields. Owing to the recent epidemic of the PstS10 race group in Europe, yellow rust has become a reoccurring disease in Norway since 2014. As all stage resistances (ASR) (or seedling resistances) are usually easily overcome by pathogen evolution, deployment of durable adult plant resistance (APR) is crucial for yellow rust resistance breeding. In this study, we assessed a Nordic spring wheat association mapping panel (n = 301) for yellow rust field resistance in seventeen field trials from 2015 to 2021, including nine locations in six countries across four different continents. Nine consistent QTL were identified across continents by genome-wide association studies (GWAS). One robust QTL on the long arm of chromosome 6A, QYr.nmbu.6A, was consistently detected in nine out of the seventeen trials. Haplotype analysis of QYr.nmbu.6A confirmed significant QTL effects in all tested environments and the effect was also validated using an independent panel of new Norwegian breeding lines. Increased frequency of the resistant haplotype was found in new varieties and breeding lines in comparison to older varieties and landraces, implying that the resistance might have been selected for due to the recent changes in the yellow rust pathogen population in Europe.

Association mapping reveals novel genes and genomic regions controlling grain size architecture in mini core accessions of Indian National Genebank wheat germplasm collection

Wheat (Triticum aestivum L.) is a staple food crop for the global human population, and thus wheat breeders are consistently working to enhance its yield worldwide. In this study, we utilized a sub-set of Indian wheat mini core germplasm to underpin the genetic architecture for seed shape-associated traits. The wheat mini core subset (125 accessions) was genotyped using 35K SNP array and evaluated for grain shape traits such as grain length (GL), grain width (GW), grain length, width ratio (GLWR), and thousand grain weight (TGW) across the seven different environments (E₁, E₂, E₃, E₄, E₅, E₅, E₆, and E₇). Marker-trait associations were determined using a multi-locus random-SNP-effect Mixed Linear Model (mrMLM) program. A total of 160 non-redundant quantitative trait nucleotides (QTNs) were identified for four grain shape traits using two or more GWAS models. Among these 160 QTNs, 27, 36, 38, and 35 QTNs were associated for GL, GW, GLWR, and TGW respectively while 24 QTNs were associated with more than one trait. Of these 160 QTNs, 73 were detected in two or more environments and were considered reliable QTLs for the respective traits. A total of 135 associated QTNs were annotated and located within the genes, including ABC transporter, Cytochrome450, Thioredoxin_M-type, and hypothetical proteins. Furthermore, the expression pattern of annotated QTNs demonstrated that only 122 were differentially expressed, suggesting these could potentially be related to seed development. The genomic regions/candidate genes for grain size traits identified in the present study represent valuable genomic resources that can potentially be utilized in the markers-assisted breeding programs to develop high-yielding varieties.

Physical mapping of the wheat genes in low-recombination regions: radiation hybrid mapping of the C-locus

KEY MESSAGE

This work reports the physical mapping of an important gene affecting spike compactness located in a low-recombination region of hexaploid wheat. This work paves the way for the eventual isolation and characterization of the factor involved but also opens up possibilities to use this approach to precisely map other wheat genes located on proximal parts of wheat chromosomes that show highly reduced recombination. Mapping wheat genes, in the centromeric and pericentromeric regions (~ 2/3rd of a given chromosome), poses a formidable challenge due to highly suppressed recombination. Using an example of compact spike locus (C-locus), this study provides an approach to precisely map wheat genes in the pericentromeric and centromeric regions that house ~ 30% of wheat genes. In club-wheat, spike compactness is controlled by the dominant C-locus, but previous efforts have failed to localize it, on a particular arm of chromosome 2D. We integrated radiation hybrid (RH) and high-resolution genetic mapping to locate C-locus on the short arm of chromosome 2D. Flanking markers of the C-locus span a physical distance of 11.0 Mb (231.0-242 Mb interval) and contain only 11 high-confidence annotated genes. This work demonstrates the value of this integrated strategy in mapping dominant genes in the low-recombination regions of the wheat genome. A comparison of the mapping resolutions of the RH and genetic maps using common anchored markers indicated that the RH map provides ~ 9 times better resolution that the genetic map even with much smaller population size. This study provides a broadly applicable approach to fine map wheat genes in regions of suppressed recombination.

Symbiotic Variations among Wheat Genotypes and Detection of Quantitative Trait Loci for Molecular Interaction with Auxin-Producing Azospirillum PGPR

Crop varieties differ in their ability to interact with Plant Growth-Promoting Rhizobacteria (PGPR), but the genetic basis for these differences is unknown. This issue was addressed with the PGPR Azospirillum baldaniorum Sp245, using 187 wheat accessions. We screened the accessions based on the seedling colonization by the PGPR and the expression of the phenylpyruvate decarboxylase gene ppdC (for synthesis of the auxin indole-3-acetic acid), using gusA fusions. Then, the effects of the PGPR on the selected accessions stimulating Sp245 (or not) were compared in soil under stress. Finally, a genome-wide association approach was implemented to identify the quantitative trait loci (QTL) associated with PGPR interaction. Overall, the ancient genotypes were more effective than the modern genotypes for Azospirillum root colonization and ppdC expression. In non-sterile soil, A. baldaniorum Sp245 improved wheat performance for three of the four PGPR-stimulating genotypes and none of the four non-PGPR-stimulating genotypes. The genome-wide association did not identify any region for root colonization but revealed 22 regions spread on 11 wheat chromosomes for ppdC expression and/or ppdC induction rate. This is the first QTL study focusing on molecular interaction with PGPR bacteria. The molecular markers identified provide the possibility to improve the capacity of modern wheat genotypes to interact with Sp245, as well as, potentially, other Azospirillum strains.

Mapping and validation of a novel major QTL for resistance to stripe rust in four wheat populations derived from landrace Qishanmai

Wheat yield has been constrained by stripe rust disease globally. A wheat landrace (Qishanmai, QSM) consistently showed lower stripe rust severities in multiple year studies than susceptible check varieties including Suwon11 (SW) at the adult plant stage. To detect QTL for reducing the severity in QSM, 1218 recombinant inbred lines (RILs) were developed from SW × QSM. QTL detection was conducted firstly using 112 RILs selected for similarity in pheno-morphological characters. The 112 RILs were assessed for stripe rust severity at the 2nd leaf, 6th leaf and flag leaf stages under field and greenhouse conditions, and genotyping was done primarily with a single nucleotide polymorphism (SNP) array. On the basis of these phenotypic and genotypic data, a major QTL (QYr.cau-1DL) was detected on chromosome 1D at the 6th leaf and flag leaf stages. Further mapping was conducted by genotyping 1218 RILs using new simple sequence repeat (SSR) markers, which were developed by referring to the sequences of the wheat line Chinese Spring (IWGSC RefSeq v1.0). QYr.cau-1DL was mapped within a 0.5 cM (5.2 Mb) interval delimited by the SSR markers 1D-320.58 and 1D-325.79. These markers were applied to select for QYr.cau-1DL by screening F₂ or BC₄F₂ plants of the wheat crosses RL6058 × QSM, Lantian10 × QSM and Yannong21 × QSM. F_2:3 or BC₄F_2:3 families derived from the selected plants were assessed for stripe rust resistance in the fields of two locations and in a greenhouse. Wheat plants carrying the resistant marker haplotype in homozygous state for QYr.cau-1DL showed lower stripe rust severities (by 44% to 48%) than plants lacking this QTL. The trial of RL6058 (a carrier of Yr18) × QSM also indicated that QYr.cau-1DL had larger effect than Yr18 on reducing severity; they acted synergistically, yielding an elevated level of stripe rust resistance.

Optimizing genomic selection of agricultural traits using K-wheat core collection

The agricultural traits that constitute basic plant breeding information are usually quantitative or complex in nature. This quantitative and complex combination of traits complicates the process of selection in breeding. This study examined the potential of genome-wide association studies (GWAS) and genomewide selection (GS) for breeding ten agricultural traits by using genome-wide SNPs. As a first step, a trait-associated candidate marker was identified by GWAS using a genetically diverse 567 Korean (K)-wheat core collection. The accessions were genotyped using an Axiom^® 35K wheat DNA chip, and ten agricultural traits were determined (awn color, awn length, culm color, culm length, ear color, ear length, days to heading, days to maturity, leaf length, and leaf width). It is essential to sustain global wheat production by utilizing accessions in wheat breeding. Among the traits associated with awn color and ear color that showed a high positive correlation, a SNP located on chr1B was significantly associated with both traits. Next, GS evaluated the prediction accuracy using six predictive models (G-BLUP, LASSO, BayseA, reproducing kernel Hilbert space, support vector machine (SVM), and random forest) and various training populations (TPs). With the exception of the SVM, all statistical models demonstrated a prediction accuracy of 0.4 or better. For the optimization of the TP, the number of TPs was randomly selected (10%, 30%, 50% and 70%) or divided into three subgroups (CC-sub 1, CC-sub 2 and CC-sub 3) based on the subpopulation structure. Based on subgroup-based TPs, better prediction accuracy was found for awn color, culm color, culm length, ear color, ear length, and leaf width. A variety of Korean wheat cultivars were used for validation to evaluate the prediction ability of populations. Seven out of ten cultivars showed phenotype-consistent results based on genomics-evaluated breeding values (GEBVs) calculated by the reproducing kernel Hilbert space (RKHS) predictive model. Our research provides a basis for improving complex traits in wheat breeding programs through genomics assisted breeding. The results of our research can be used as a basis for improving wheat breeding programs by using genomics-assisted breeding.

Genome-wide characterization of the VQ genes in Triticeae and their functionalization driven by polyploidization and gene duplication events in wheat

Valine-glutamine motif-containing (VQ) proteins are transcriptional cofactors widely involved in plant growth, development, and response to various stresses. Although the VQ family has been genome-wide identified in some species, but the knowledge regarding duplication-driven functionalization of VQ genes among evolutionarily related species is still lacking. Here, 952 VQ genes have been identified from 16 species, emphasizing seven Triticeae species including the bread wheat. Comprehensive phylogenetic and syntenic analyses allow us to establish the orthologous relationship of VQ genes from rice (Oryza sativa) to bread wheat (Triticum aestivum). The evolutionary analysis revealed that whole-genome duplication (WGD) drives the expansion of OsVQs, while TaVQs expansion is associated with a recent burst of gene duplication (RBGD). We also analyzed the motif composition and molecular properties of TaVQ proteins, enriched biological functions, and expression patterns of TaVQs. We demonstrate that WGD-derived TaVQs have become divergent in both protein motif composition and expression pattern, while RBGD-derived TaVQs tend to adopt specific expression patterns, suggesting their functionalization in certain biological processes or in response to specific stresses. Furthermore, some RBGD-derived TaVQs are found to be associated with salt tolerance. Several of the identified salt-related TaVQ proteins were located in the cytoplasm and nucleus and their salt-responsive expression patterns were validated by qPCR analysis. Yeast-based functional experiments confirmed that TaVQ27 may be a new regulator to salt response and regulation. Overall, this study lays the foundation for further functional validation of VQ family members within the Triticeae species.

Population Structure of Modern Winter Wheat Accessions from Central Asia

Despite the importance of winter wheat in Central Asian countries, there are limited reports describing their diversity within this region. In this study, the population structures of 115 modern winter wheat cultivars from four Central Asian countries were compared to germplasms from six other geographic origins using 10,746 polymorphic single-nucleotide polymorphism (SNP) markers. After applying the STRUCTURE package, we found that in terms of the most optimal K steps, samples from Kazakhstan and Kyrgyzstan were grouped together with samples from Russia, while samples from Tajikistan and Uzbekistan were grouped with samples from Afghanistan. The mean value of Nei's genetic diversity index for the germplasm from four groups from Central Asia was 0.261, which is comparable to that of the six other groups studied: Europe, Australia, the USA, Afghanistan, Turkey, and Russia. The Principal Coordinate Analysis (PCoA) showed that samples from Kyrgyzstan, Tajikistan, and Uzbekistan were close to samples from Turkey, while Kazakh accessions were located near samples from Russia. The evaluation of 10,746 SNPs in Central Asian wheat suggested that 1006 markers had opposing allele frequencies. Further assessment of the physical positions of these 1006 SNPs in the Wheat Ensembl database indicated that most of these markers are constituents of genes associated with plant stress tolerance and adaptability. Therefore, the SNP markers identified can be effectively used in regional winter wheat breeding projects for facilitating plant adaptation and stress resistance.

Uncovering the transcriptional regulatory network involved in boosting wheat regeneration and transformation

Genetic transformation is important for gene functional study and crop improvement. However, it is less effective in wheat. Here we employed a multi-omic analysis strategy to uncover the transcriptional regulatory network (TRN) responsible for wheat regeneration. RNA-seq, ATAC-seq and CUT&Tag techniques were utilized to profile the transcriptional and chromatin dynamics during early regeneration from the scutellum of immature embryos in the wheat variety Fielder. Our results demonstrate that the sequential expression of genes mediating cell fate transition during regeneration is induced by auxin, in coordination with changes in chromatin accessibility, H3K27me3 and H3K4me3 status. The built-up TRN driving wheat regeneration was found to be dominated by 446 key transcription factors (TFs). Further comparisons between wheat and Arabidopsis revealed distinct patterns of DNA binding with one finger (DOF) TFs in the two species. Experimental validations highlighted TaDOF5.6 (TraesCS6A02G274000) and TaDOF3.4 (TraesCS2B02G592600) as potential enhancers of transformation efficiency in different wheat varieties.

Wheat derived glucuronokinase as a potential target for regulating ascorbic acid and phytic acid content with increased root length under drought and ABA stresses in Arabidopsis thaliana

Glucuronokinase (GlcAK) converts glucuronic acid into glucuronic acid-1-phosphate, which is then converted into UDP-glucuronic acid (UDP-GlcA) via myo-inositol oxygenase (MIOX) pathway. UDP-GlcA acts as a precursor in the synthesis of nucleotide-sugar moieties forming cell wall biomass. GlcAK being present at the bifurcation point between UDP-GlcA and ascorbic acid (AsA) biosyntheses, makes it necessary to study its role in plants. In this study, the three homoeologs of GlcAK gene from hexaploid wheat were overexpressed in Arabidopsis thaliana. The GlcAK overexpressing transgenic lines showed decreased contents of AsA and phytic acid (PA) as compared to control plants. Root length and seed germination analyses under abiotic stress (drought and abscisic acid) conditions revealed enhanced root length in transgenic lines as compared to control plants. These results indicate that the MIOX pathway might be contributing towards AsA biosynthesis as evident by the decreased AsA content in the GlcAK overexpressing transgenic Arabidopsis thaliana plants. Findings of the present study will enhance the understanding of the involvement of GlcAK gene in MIOX pathway and subsequent physiological effects in plants.

TaTPP-7A positively feedback regulates grain filling and wheat grain yield through T6P-SnRK1 signalling pathway and sugar-ABA interaction

Grain size and filling are two key determinants of grain thousand-kernel weight (TKW) and crop yield, therefore they have undergone strong selection since cereal was domesticated. Genetic dissection of the two traits will improve yield potential in crops. A quantitative trait locus significantly associated with wheat grain TKW was detected on chromosome 7AS flanked by a simple sequence repeat marker of Wmc17 in Chinese wheat 262 mini-core collection by genome-wide association study. Combined with the bulked segregant RNA-sequencing (BSR-seq) analysis of an F₂ genetic segregation population with extremely different TKW traits, a candidate trehalose-6-phosphate phosphatase gene located at 135.0 Mb (CS V1.0), designated as TaTPP-7A, was identified. This gene was specifically expressed in developing grains and strongly influenced grain filling and size. Overexpression (OE) of TaTPP-7A in wheat enhanced grain TKW and wheat yield greatly. Detailed analysis revealed that OE of TaTPP-7A significantly increased the expression levels of starch synthesis- and senescence-related genes involved in abscisic acid (ABA) and ethylene pathways. Moreover, most of the sucrose metabolism and starch regulation-related genes were potentially regulated by SnRK1. In addition, TaTPP-7A is a crucial domestication- and breeding-targeted gene and it feedback regulates sucrose lysis, flux, and utilization in the grain endosperm mainly through the T6P-SnRK1 pathway and sugar-ABA interaction. Thus, we confirmed the T6P signalling pathway as the central regulatory system for sucrose allocation and source-sink interactions in wheat grains and propose that the trehalose pathway components have great potential to increase yields in cereal crops.

Suppressing ASPARTIC PROTEASE 1 prolongs photosynthesis and increases wheat grain weight

The elongation of photosynthesis, or functional staygreen, represents a feasible strategy to propel metabolite flux towards cereal kernels. However, achieving this goal remains a challenge in food crops. Here we report the cloning of wheat CO2 assimilation and kernel enhanced 2 (cake2), the mechanism underlying the photosynthesis advantages and natural alleles amenable to breeding elite varieties. A premature stop mutation in the A-genome copy of the ASPARTIC PROTEASE 1 (APP-A1) gene increased the photosynthesis rate and yield. APP1 bound and degraded PsbO, the protective extrinsic member of photosystem II critical for increasing photosynthesis and yield. Furthermore, a natural polymorphism of the APP-A1 gene in common wheat reduced APP-A1's activity and promoted photosynthesis and grain size and weight. This work demonstrates that the modification of APP1 increases photosynthesis, grain size and yield potentials. The genetic resources could propel photosynthesis and high-yield potentials in elite varieties of tetraploid and hexaploid wheat.

ZIP4 is required for normal progression of synapsis and for over 95% of crossovers in wheat meiosis

Tetraploid (AABB) and hexaploid (AABBDD) wheat have multiple sets of similar chromosomes, with successful meiosis and preservation of fertility relying on synapsis and crossover (CO) formation only taking place between homologous chromosomes. In hexaploid wheat, the major meiotic gene TaZIP4-B2 (Ph1) on chromosome 5B, promotes CO formation between homologous chromosomes, whilst suppressing COs between homeologous (related) chromosomes. In other species, ZIP4 mutations eliminate approximately 85% of COs, consistent with loss of the class I CO pathway. Tetraploid wheat has three ZIP4 copies: TtZIP4-A1 on chromosome 3A, TtZIP4-B1 on 3B and TtZIP4-B2 on 5B. Here, we have developed single, double and triple zip4 TILLING mutants and a CRISPR Ttzip4-B2 mutant, to determine the effect of ZIP4 genes on synapsis and CO formation in the tetraploid wheat cultivar 'Kronos'. We show that disruption of two ZIP4 gene copies in Ttzip4-A1B1 double mutants, results in a 76-78% reduction in COs when compared to wild-type plants. Moreover, when all three copies are disrupted in Ttzip4-A1B1B2 triple mutants, COs are reduced by over 95%, suggesting that the TtZIP4-B2 copy may also affect class II COs. If this is the case, the class I and class II CO pathways may be interlinked in wheat. When ZIP4 duplicated and diverged from chromosome 3B on wheat polyploidization, the new 5B copy, TaZIP4-B2, could have acquired an additional function to stabilize both CO pathways. In tetraploid plants deficient in all three ZIP4 copies, synapsis is delayed and does not complete, consistent with our previous studies in hexaploid wheat, when a similar delay in synapsis was observed in a 59.3 Mb deletion mutant, ph1b, encompassing the TaZIP4-B2 gene on chromosome 5B. These findings confirm the requirement of ZIP4-B2 for efficient synapsis, and suggest that TtZIP4 genes have a stronger effect on synapsis than previously described in Arabidopsis and rice. Thus, ZIP4-B2 in wheat accounts for the two major phenotypes reported for Ph1, promotion of homologous synapsis and suppression of homeologous COs.

Wnt binding to Coatomer proteins directs secretion on exosomes independently of palmitoylation

Wnt proteins are secreted hydrophobic glycoproteins that act over long distances through poorly understood mechanisms. We discovered that Wnt7a is secreted on extracellular vesicles (EVs) following muscle injury. Structural analysis identified the motif responsible for Wnt7a secretion on EVs that we term the Exosome Binding Peptide (EBP). Addition of the EBP to an unrelated protein directed secretion on EVs. Disruption of palmitoylation, knockdown of WLS, or deletion of the N-terminal signal peptide did not affect Wnt7a secretion on purified EVs. Bio-ID analysis identified Coatomer proteins as candidates responsible for loading Wnt7a onto EVs. The crystal structure of EBP bound to the COPB2 coatomer subunit, the binding thermodynamics, and mutagenesis experiments, together demonstrate that a dilysine motif in the EBP mediates binding to COPB2. Other Wnts contain functionally analogous structural motifs. Mutation of the EBP results in a significant impairment in the ability of Wnt7a to stimulate regeneration, indicating that secretion of Wnt7a on exosomes is critical for normal regeneration in vivo . Our studies have defined the structural mechanism that mediates binding of Wnt7a to exosomes and elucidated the singularity of long-range Wnt signalling.

The translational landscape of bread wheat during grain development

The dynamics of gene expression in crop grains has typically been investigated at the transcriptional level. However, this approach neglects translational regulation, a widespread mechanism that rapidly modulates gene expression to increase the plasticity of organisms. Here, we performed ribosome profiling and polysome profiling to obtain a comprehensive translatome data set of developing bread wheat (Triticum aestivum) grains. We further investigated the genome-wide translational dynamics during grain development, revealing that the translation of many functional genes is modulated in a stage-specific manner. The unbalanced translation between subgenomes is pervasive, which increases the expression flexibility of allohexaploid wheat. In addition, we uncovered widespread previously unannotated translation events, including upstream open reading frames (uORFs), downstream open reading frames (dORFs), and open reading frames (ORFs) in long noncoding RNAs, and characterized the temporal expression dynamics of small ORFs. We demonstrated that uORFs act as cis-regulatory elements that can repress or even enhance the translation of mRNAs. Gene translation may be combinatorially modulated by uORFs, dORFs, and microRNAs. In summary, our study presents a translatomic resource that provides a comprehensive and detailed overview of the translational regulation in developing bread wheat grains. This resource will facilitate future crop improvements for optimal yield and quality.

Genome-wide association study reveals the genetic variation and candidate gene for grain calcium content in bread wheat

KEY MESSAGE

This study provides important information on the genetic basis of GCaC in wheat, thus contributing to breeding efforts to improve the nutrient quality of wheat. Calcium (Ca) plays important roles in the human body. Wheat grain provides the main diet for billions of people worldwide but is low in Ca content. Here, grain Ca content (GCaC) of 471 wheat accessions was determined in four field environments. A genome-wide association study (GWAS) was performed to reveal the genetic basis of GCaC using the phenotypic data form four environments and a wheat 660 K single nucleotide polymorphism (SNP) array. Twelve quantitative trait locus (QTLs) for GCaC were identified on chromosomes 1A, 1D, 2A, 3B, 6A, 6D, 7A, and 7D, which was significant in at least two environments. Haplotype analysis revealed that the phenotypic difference between the haplotypes of TraesCS6D01G399100 was significant (P ≤ 0.05) across four environments, suggesting it as an important candidate gene for GCaC. This research enhances our understanding of the genetic architecture of GCaC for further improving the nutrient quality of wheat.

The wheat stem rust resistance gene Sr43 encodes an unusual protein kinase

To safeguard bread wheat against pests and diseases, breeders have introduced over 200 resistance genes into its genome, thus nearly doubling the number of designated resistance genes in the wheat gene pool¹. Isolating these genes facilitates their fast-tracking in breeding programs and incorporation into polygene stacks for more durable resistance. We cloned the stem rust resistance gene Sr43, which was crossed into bread wheat from the wild grass Thinopyrum elongatum^2,3. Sr43 encodes an active protein kinase fused to two domains of unknown function. The gene, which is unique to the Triticeae, appears to have arisen through a gene fusion event 6.7 to 11.6 million years ago. Transgenic expression of Sr43 in wheat conferred high levels of resistance to a wide range of isolates of the pathogen causing stem rust, highlighting the potential value of Sr43 in resistance breeding and engineering.

ASYNAPSIS 1 ensures crossover fidelity in polyploid wheat by promoting homologous recombination and suppressing non-homologous recombination

During meiosis, the chromosome axes and synaptonemal complex mediate chromosome pairing and homologous recombination to maintain genomic stability and accurate chromosome segregation. In plants, ASYNAPSIS 1 (ASY1) is a key component of the chromosome axis that promotes inter-homolog recombination, synapsis and crossover formation. Here, the function of ASY1 has been cytologically characterized in a series of hypomorphic wheat mutants. In tetraploid wheat, asy1 hypomorphic mutants experience a reduction in chiasmata (crossovers) in a dosage-specific manner, resulting in failure to maintain crossover (CO) assurance. In mutants with only one functional copy of ASY1, distal chiasmata are maintained at the expense of proximal and interstitial chiasmata, indicating that ASY1 is required to promote chiasma formation away from the chromosome ends. Meiotic prophase I progression is delayed in asy1 hypomorphic mutants and is arrested in asy1 null mutants. In both tetraploid and hexaploid wheat, single asy1 mutants exhibit a high degree of ectopic recombination between multiple chromosomes at metaphase I. To explore the nature of the ectopic recombination, Triticum turgidum asy1b-2 was crossed with wheat-wild relative Aegilops variabilis. Homoeologous chiasmata increased 3.75-fold in Ttasy1b-2/Ae. variabilis compared to wild type/Ae. variabilis, indicating that ASY1 suppresses chiasma formation between divergent, but related chromosomes. These data suggest that ASY1 promotes recombination along the chromosome arms of homologous chromosomes whilst suppressing recombination between non-homologous chromosomes. Therefore, asy1 mutants could be utilized to increase recombination between wheat wild relatives and elite varieties for expediting introgression of important agronomic traits.

An unusual tandem kinase fusion protein confers leaf rust resistance in wheat

The introgression of chromosome segments from wild relatives is an established strategy to enrich crop germplasm with disease-resistance genes¹. Here we use mutagenesis and transcriptome sequencing to clone the leaf rust resistance gene Lr9, which was introduced into bread wheat from the wild grass species Aegilops umbellulata². We established that Lr9 encodes an unusual tandem kinase fusion protein. Long-read sequencing of a wheat Lr9 introgression line and the putative Ae. umbellulata Lr9 donor enabled us to assemble the ~28.4-Mb Lr9 translocation and to identify the translocation breakpoint. We likewise cloned Lr58, which was reportedly introgressed from Aegilops triuncialis³, but has an identical coding sequence compared to Lr9. Cytogenetic and haplotype analyses corroborate that the two genes originate from the same translocation event. Our work sheds light on the emerging role of kinase fusion proteins in wheat disease resistance, expanding the repertoire of disease-resistance genes for breeding.

FHB resistance conferred by Fhb1 is under inhibitory regulation of two genetic loci in wheat (Triticum aestivum L.)

KEY MESSAGE

Two loci inhibiting Fhb1 resistance to Fusarium head blight were identified through genome-wide association mapping and validated in biparental populations. Fhb1 confers Fusarium head blight (FHB) resistance by limiting fungal spread within spikes in wheat (type II resistance). However, not all lines with Fhb1 display the expected resistance. To identify genetic factors regulating Fhb1 effect, a genome-wide association study for type II resistance was first performed with 72 Fhb1-carrying lines using the Illumina 90 K iSelect SNP chip. Of 84 significant marker-trait associations detected, more than half were repeatedly detected in at least two environments, with the SNPs distributed in one region on chromosome 5B and one on chromosome 6A. This result was validated in a collection of 111 lines with Fhb1 and 301 lines without Fhb1. We found that these two loci caused significant resistance variations solely among lines with Fhb1 by compromising the resistance. In1, the inhibitory gene on chromosome 5B, was in close linkage with Xwgrb3860 in a recombinant inbred line population derived from Nanda2419 × Wangshuibai and a double haploid (DH) population derived from R-43 (Fhb1 near isogenic line) × Biansui7 (with Fhb1 and In1); and In2, the inhibitory gene on chromosome 6A, was mapped to the Xwgrb4113-Xwgrb4034 interval using a DH population derived from R-43 × PH8901 (with Fhb1 and In2). In1 and In2 are present in all wheat-growing areas worldwide. Their frequencies in China's modern cultivars are high but have significantly decreased in comparison with landraces. These findings are of great significance for FHB resistance breeding using Fhb1.

Mapping of quantitative trait loci for leaf rust resistance in the wheat population 'Xinmai 26/Zhoumai 22'

Leaf rust, caused by the fungal pathogen Puccinia triticina (Pt), is one of the major and dangerous diseases of wheat, and has caused serious yield loss of wheat worldwide. Here, we investigated adult-plant resistance (APR) to leaf rust in a recombinant inbred line (RIL) population derived from 'Xinmai 26' and 'Zhoumai 22' over 3 years. Linkage mapping for APR to leaf rust revealed four quantitative trait loci (QTL) in this RIL population. Two QTL, QLr.hnau-2BS and QLr.hnau-3BS were contributed by 'Zhoumai22', whereas QLr.hnau-2DS and QLr.hnau-5AL were contributed by 'Xinmai 26'. The QLr.hnau-2BS covering a race-specific resistance gene Lr13 showed the most stable APR to leaf rust. Overexpression of Lr13 significantly increased APR to leaf rust. Interestingly, we found that a CNL(coiled coil-nucleotide-binding site-leucine-rich repeat)-like gene, TaCN, in QLr.hnau-2BS completely co-segregated with leaf rust resistance. The resistant haplotype TaCN-R possessed half the sequence of the coiled-coil domain of TaCN protein. Lr13 strongly interacted with TaCN-R, but did not interact with the full-length TaCN (TaCN-S). In addition, TaCN-R was significantly induced after Pt inoculation and changed the sub-cellular localization of Lr13 after interaction. Therefore, we hypothesized that TaCN-R mediated leaf rust resistance possibly by interacting with Lr13. This study provides important QTL for APR to leaf rust, and new insights into understanding how a CNL gene modulates disease resistance in common wheat.

Proteomic Pathways across Ejection Fraction Spectrum in Heart Failure: an EXSCEL Substudy

Background

Ejection fraction (EF) is a key component of heart failure (HF) classification, including the increasingly codified HF with mildly reduced EF (HFmrEF) category. However, the biologic basis of HFmrEF as an entity distinct from HF with preserved EF (HFpEF) and reduced EF (HFrEF) has not been well characterized.

Methods

The EXSCEL trial randomized participants with type 2 diabetes (T2DM) to once-weekly exenatide (EQW) vs. placebo. For this study, profiling of ∼5000 proteins using the SomaLogic SomaScan platform was performed in baseline and 12-month serum samples from N=1199 participants with prevalent HF at baseline. Principal component analysis (PCA) and ANOVA (FDR p<0.1) were used to determine differences in proteins between three EF groups, as previously curated in EXSCEL (EF>55% [HFpEF], EF 40-55% [HFmrEF], EF<40% [HFrEF]). Cox proportional hazards was used to assess association between baseline levels of significant proteins, and changes in protein level between baseline and 12-month, with time-to-HF hospitalization. Mixed models were used to assess whether significant proteins changed differentially with exenatide vs. placebo therapy.

Results

Of N=1199 EXSCEL participants with prevalent HF, 284 (24%), 704 (59%) and 211 (18%) had HFpEF, HFmrEF and HFrEF, respectively. Eight PCA protein factors and 221 individual proteins within these factors differed significantly across the three EF groups. Levels of the majority of proteins (83%) demonstrated concordance between HFmrEF and HFpEF, but higher levels in HFrEF, predominated by the domain of extracellular matrix regulation, e.g. COL28A1 and tenascin C [TNC]; p<0.0001. Concordance between HFmrEF and HFrEF was observed in a minority of proteins (1%) including MMP-9 (p<0.0001). Biologic pathways of epithelial mesenchymal transition, ECM receptor interaction, complement and coagulation cascades, and cytokine receptor interaction demonstrated enrichment among proteins with the dominant pattern, i.e. HFmrEF-HFpEF concordance. Baseline levels of 208 (94%) of the 221 proteins were associated with time-to-incident HF hospitalization including domains of extracellular matrix (COL28A1, TNC), angiogenesis (ANG2, VEGFa, VEGFd), myocyte stretch (NT-proBNP), and renal function (cystatin-C). Change in levels of 10 of the 221 proteins from baseline to 12 months (including increase in TNC) predicted incident HF hospitalization (p<0.05). Levels of 30 of the 221 significant proteins (including TNC, NT-proBNP, ANG2) were reduced differentially by EQW compared with placebo (interaction p<0.0001).

Conclusions

In this HF substudy of a large clinical trial of people with T2DM, we found that serum levels of most proteins across multiple biologic domains were similar between HFmrEF and HFpEF. HFmrEF may be more biologically similar to HFpEF than HFrEF, and specific related biomarkers may offer unique data on prognosis and pharmacotherapy modification with variability by EF.