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

Deleterious and Adaptive Mutations in Plant Germplasm Conserved Ex Situ

Conserving more than 7 million plant germplasm accessions in 1,750 genebanks worldwide raises the hope of securing the food supply for humanity for future generations. However, there is a genetic cost for such long-term germplasm conservation, which has been largely unaccounted for before. We investigated the extent and variation of deleterious and adaptive mutations in 490 individual plants representing barley, wheat, oat, soybean, maize, rapa, and sunflower collections in a seed genebank using RNA-Seq technology. These collections were found to have a range of deleterious mutations detected from 125 (maize) to 83,695 (oat) with a mean of 13,537 and of the averaged sample-wise mutation burden per deleterious locus from 0.069 to 0.357 with a mean of 0.200. Soybean and sunflower collections showed that accessions acquired earlier had increased mutation burdens. The germplasm with more years of storage in several collections carried more deleterious and fewer adaptive mutations. The samples with more cycles of germplasm regeneration revealed fewer deleterious and more adaptive mutations. These findings are significant for understanding mutational dynamics and genetic cost in conserved germplasm and have implications for long-term germplasm management and conservation.

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.

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

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

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.

High-quality assembly and methylome of a Tibetan wild tree peony genome (Paeonia ludlowii) reveal the evolution of giant genome architecture

Tree peony belongs to one of the Saxifragales families, Paeoniaceae. It is one of the most famous ornamental plants, and is also a promising woody oil plant. Although two Paeoniaceae genomes have been released, their assembly qualities are still to be improved. Additionally, more genomes from wild peonies are needed to accelerate genomic-assisted breeding. Here we assemble a high-quality and chromosome-scale 10.3-Gb genome of a wild Tibetan tree peony, Paeonia ludlowii, which features substantial sequence divergence, including around 75% specific sequences and gene-level differentials compared with other peony genomes. Our phylogenetic analyses suggest that Saxifragales and Vitales are sister taxa and, together with rosids, they are the sister taxon to asterids. The P. ludlowii genome is characterized by frequent chromosome reductions, centromere rearrangements, broadly distributed heterochromatin, and recent continuous bursts of transposable element (TE) movement in peony, although it lacks recent whole-genome duplication. These recent TE bursts appeared during the uplift and glacial period of the Qinghai-Tibet Plateau, perhaps contributing to adaptation to rapid climate changes. Further integrated analyses with methylome data revealed that genome expansion in peony might be dynamically affected by complex interactions among TE proliferation, TE removal, and DNA methylation silencing. Such interactions also impact numerous recently duplicated genes, particularly those related to oil biosynthesis and flower traits. This genome resource will not only provide the genomic basis for tree peony breeding but also shed light on the study of the evolution of huge genome structures as well as their protein-coding genes.

Paternally imprinted LATE-FLOWERING2 transcription factor contributes to paternal-excess interploidy hybridization barriers in wheat

Interploidy hybridization between hexaploid and tetraploid genotypes occurred repeatedly during genomic introgression events throughout wheat evolution, and is commonly employed in wheat breeding programs. Hexaploid wheat usually serves as maternal parent because the reciprocal cross generates progeny with severe defects and poor seed germination, but the underlying mechanism is poorly understood. Here, we performed detailed analysis of phenotypic variation in endosperm between two interploidy reciprocal crosses arising from tetraploid (Triticum durum, AABB) and hexaploid wheat (Triticum aestivum, AABBDD). In the paternal- versus the maternal-excess cross, the timing of endosperm cellularization was delayed and starch granule accumulation in the endosperm was repressed, causing reduced germination percentage. The expression profiles of genes involved in nutrient metabolism differed strongly between these endosperm types. Furthermore, expression patterns of parental alleles were dramatically disturbed in interploidy versus intraploidy crosses, leading to increased number of imprinted genes. The endosperm-specific TaLFL2 showed a paternally imprinted expression pattern in interploidy crosses partially due to allele-specific DNA methylation. Paternal TaLFL2 binds to and represses a nutrient accumulation regulator TaNAC019, leading to reduced storage protein and starch accumulation during endosperm development in paternal-excess cross, as confirmed by interploidy crosses between tetraploid wild-type and clustered regularly interspaced palindromic repeats (CRISPR) - CRISPR-associated protein 9 generated hexaploid mutants. These findings reveal a contribution of genomic imprinting to paternal-excess interploidy hybridization barriers during wheat evolution history and explains why experienced breeders preferentially exploit maternal-excess interploidy crosses in wheat breeding programs.

Whole-genome sequencing of diverse wheat accessions uncovers genetic changes during modern breeding in China and the United States

Breeding has dramatically changed the plant architecture of wheat (Triticum aestivum), resulting in the development of high-yielding varieties adapted to modern farming systems. However, how wheat breeding shaped the genomic architecture of this crop remains poorly understood. Here, we performed a comprehensive comparative analysis of a whole-genome resequencing panel of 355 common wheat accessions (representing diverse landraces and modern cultivars from China and the United States) at the phenotypic and genomic levels. The genetic diversity of modern wheat cultivars was clearly reduced compared to landraces. Consistent with these genetic changes, most phenotypes of cultivars from China and the United States were significantly altered. Of the 21 agronomic traits investigated, 8 showed convergent changes between the 2 countries. Moreover, of the 207 loci associated with these 21 traits, more than half overlapped with genomic regions that showed evidence of selection. The distribution of selected loci between the Chinese and American cultivars suggests that breeding for increased productivity in these 2 regions was accomplished by pyramiding both shared and region-specific variants. This work provides a framework to understand the genetic architecture of the adaptation of wheat to diverse agricultural production environments, as well as guidelines for optimizing breeding strategies to design better wheat varieties.

Physiological, Epigenetic, and Transcriptome Analyses Provide Insights into the Responses of Wheat Seedling Leaves to Different Water Depths under Flooding Conditions

Flooding stress, including waterlogging and submergence, is one of the major abiotic stresses that seriously affects the growth and development of plants. In the present study, physiological, epigenetic, and transcriptomic analyses were performed in wheat seedling leaves under waterlogging (WL), half submergence (HS), and full submergence (FS) treatments. The results demonstrate that FS increased the leaves' hydrogen peroxide (H2O2) and malondialdehyde (MDA) contents and reduced their chlorophyll contents (SPAD), photosynthetic efficiency (Fv/Fm), and shoot dry weight more than HS and WL. In addition, FS increased catalase (CAT) and peroxidase (POD) activities more than HS and WL. However, there were no significant differences in the contents of H2O2, MDA, SPAD, and Fv/Fm, and the activities of superoxide dismutase (SOD) and POD between the HS and WL treatments. The changes in DNA methylation were related to stress types, increasing under the WL and HS treatments and decreasing under the FS treatment. Additionally, a total of 9996, 10,619, and 24,949 genes were differentially expressed under the WL, HS, and FS treatments, respectively, among which the 'photosynthesis', 'phenylpropanoid biosynthesis', and 'plant hormone signal transduction' pathways were extensively enriched under the three flooding treatments. The genes involved in these pathways showed flooding-type-specific expression. Moreover, flooding-type-specific responses were observed in the three conditions, including the enrichment of specific TFs and response pathways. These results will contribute to a better understanding of the molecular mechanisms underlying the responses of wheat seedling leaves to flooding stress and provide valuable genetic and epigenetic information for breeding flood-tolerant varieties of wheat.

Rht12b, a widely used ancient allele of TaGA2oxA13, reduces plant height and enhances yield potential in wheat

KEY MESSAGE

We identified a new wheat dwarfing allele Rht12b conferring reduced height and higher grain yield, pinpointed its causal variations, developed a breeding-applicable marker, and traced its origin and worldwide distribution. Plant height control is essential to optimize lodging resistance and yield gain in crops. RHT12 is a reduced height (Rht) locus that is identified in a mutationally induced dwarfing mutant and encodes a gibberellin 2-oxidase TaGA2oxA13. However, the artificial dwarfing allele is not used in wheat breeding due to excessive height reduction. Here, we confirmed a stable Rht locus, overlapping with RHT12, in a panel of wheat cultivars and its dwarfing allele reduced plant height by 5.4-8.2 cm, equivalent to Rht12b, a new allele of RHT12. We validated the effect of Rht12b on plant height in a bi-parent mapping population. Importantly, wheat cultivars carrying Rht12b had higher grain yield than those with the contrasting Rht12a allele. Rht12b conferred higher expression level of TaGA2oxA13. Transient activation assays defined SNP-390(C/A) in the promoter of TaGA2oxA13 as the causal variation. An efficient kompetitive allele-specific PCR marker was developed to diagnose Rht12b. Conjoint analysis showed that Rht12b plus the widely used Rht-D1b, Rht8 and Rht24b was the predominant Rht combination and conferred a moderate plant height in tested wheat cultivars. Evolutionary tracking uncovered that RHT12 locus arose from a tandem duplication event with Rht12b firstly appearing in wild emmer. The frequency of Rht12b was approximately 70% (700/1005) in a worldwide wheat panel and comparable to or higher than those of other widely used Rht genes, suggesting it had been subjected to positive selection. These findings not only identify a valuable Rht gene for wheat improvement but also develop a functionally diagnostic tool for marker-assisted breeding.

Combined linkage analysis and association mapping identifies genomic regions associated with yield-related and drought-tolerance traits in wheat (Triticum aestivum L.)

KEY MESSAGE

Combined linkage analysis and association mapping identified genomic regions associated with yield and drought tolerance, providing information to assist breeding for high yield and drought tolerance in wheat. Wheat (Triticum aestivum L.) is one of the most widely grown food crops and provides adequate amounts of protein to support human health. Drought stress is the most important abiotic stress constraining yield during the flowering and grain development periods. Precise targeting of genomic regions underlying yield- and drought tolerance-responsive traits would assist in breeding programs. In this study, two water treatments (well-watered, WW, and rain-fed water stress, WS) were applied, and five yield-related agronomic traits (plant height, PH; spike length, SL; spikelet number per spike, SNPS; kernel number per spike, KNPS; thousand kernel weight, TKW) and drought response values (DRVs) were used to characterize the drought sensitivity of each accession. Association mapping was performed on an association panel of 304 accessions, and linkage analysis was applied to a doubled haploid (DH) population of 152 lines. Eleven co-localized genomic regions associated with yield traits and DRV were identified in both populations. Many previously cloned key genes were located in these regions. In particular, a TKW-associated region on chromosome 2D was identified using both association mapping and linkage analysis and a key candidate gene, TraesCS2D02G142500, was detected based on gene annotation and differences in expression levels. Exonic SNPs were analyzed by sequencing the full length of TraesCS2D02G142500 in the association panel, and a rare haplotype, Hap-2, which reduced TKW to a lesser extent than Hap-1 under drought stress, and the Hap-2 varieties presented drought-insensitive. Altogether, this study provides fundamental insights into molecular targets for high yield and drought tolerance in wheat.

Transposable element-initiated enhancer-like elements generate the subgenome-biased spike specificity of polyploid wheat

Transposable elements (TEs) comprise ~85% of the common wheat genome, which are highly diverse among subgenomes, possibly contribute to polyploid plasticity, but the causality is only assumed. Here, by integrating data from gene expression cap analysis and epigenome profiling via hidden Markov model in common wheat, we detect a large proportion of enhancer-like elements (ELEs) derived from TEs producing nascent noncoding transcripts, namely ELE-RNAs, which are well indicative of the regulatory activity of ELEs. Quantifying ELE-RNA transcriptome across typical developmental stages reveals that TE-initiated ELE-RNAs are mainly from RLG_famc7.3 specifically expanded in subgenome A. Acquisition of spike-specific transcription factor binding likely confers spike-specific expression of RLG_famc7.3-initiated ELE-RNAs. Knockdown of RLG_famc7.3-initiated ELE-RNAs resulted in global downregulation of spike-specific genes and abnormal spike development. These findings link TE expansion to regulatory specificity and polyploid developmental plasticity, highlighting the functional impact of TE-driven regulatory innovation on polyploid evolution.

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

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

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.

CRISPR-mediated acceleration of wheat improvement: advances and perspectives

Common wheat (Triticum aestivum) is one of the most widely cultivated and consumed crops globally. In the face of limited arable land and climate changes, it is a great challenge to maintain current and increase future wheat production. Enhancing agronomic traits in wheat by introducing mutations across all three homoeologous copies of each gene has proven to be a difficult task due to its large genome with high repetition. However, clustered regularly interspaced short palindromic repeat (CRISPR)/CRISPR-associated nuclease (Cas) genome editing technologies offer a powerful means of precisely manipulating the genomes of crop species, thereby opening up new possibilities for biotechnology and breeding. In this review, we first focus on the development and optimization of the current CRISPR-based genome editing tools in wheat, emphasizing recent breakthroughs in precise and multiplex genome editing. We then describe the general procedure of wheat genome editing and highlight different methods to deliver the genome editing reagents into wheat cells. Furthermore, we summarize the recent applications and advancements of CRISPR/Cas technologies for wheat improvement. Lastly, we discuss the remaining challenges specific to wheat genome editing and its future prospects.

Single- and multi-trait genomic prediction and genome-wide association analysis of grain yield and micronutrient-related traits in ICARDA wheat under drought environment

Globally, over 2 billion people suffer from malnutrition due to inadequate intake of micronutrients. Genomic-assisted breeding is identified as a valuable method to facilitate developing new improved plant varieties targeting grain yield and micronutrient-related traits. In this study, a genome-wide association study (GWAS) and single- and multi-trait-based genomic prediction (GP) analysis was conducted using a set of 252 elite wheat genotypes from the International Center for Agricultural Research in Dry Areas (ICARDA). The objective was to identify linked SNP markers, putative candidate genes and to evaluate the genomic estimated breeding values (GEBVs) of grain yield and micronutrient-related traits.. For this purpose, a field trial was conducted at a drought-prone station, Merchouch, Morocco for 2 consecutive years (2018 and 2019) followed by GWAS and genomic prediction analysis with 10,173 quality SNP markers. The studied genotypes exhibited a significant genotypic variation in grain yield and micronutrient-related traits. The GWAS analysis identified highly significantly associated markers and linked putative genes on chromosomes 1B and 2B for zinc (Zn) and iron (Fe) contents, respectively. The genomic predictive ability of selenium (Se) and Fe traits with the multi-trait-based GP GBLUP model was 0.161 and 0.259 improving by 6.62 and 4.44%, respectively, compared to the corresponding single-trait-based models. The identified significantly linked SNP markers, associated putative genes, and developed GP models could potentially facilitate breeding programs targeting to improve the overall genetic gain of wheat breeding for grain yield and biofortification of micronutrients via marker-assisted (MAS) and genomic selection (GS) methods.

Deciphering spike architecture formation towards yield improvement in wheat

Wheat is the most widely grown crop globally, providing 20% of the daily consumed calories and protein content around the world. With the growing global population and frequent occurrence of extreme weather caused by climate change, ensuring adequate wheat production is essential for food security. The architecture of the inflorescence plays a crucial role in determining the grain number and size, which is a key trait for improving yield. Recent advances in wheat genomics and gene cloning techniques have improved our understanding of wheat spike development and its applications in breeding practices. Here, we summarize the genetic regulation network governing wheat spike formation, the strategies used for identifying and studying the key factors affecting spike architecture, and the progress made in breeding applications. Additionally, we highlight future directions that will aid in the regulatory mechanistic study of wheat spike determination and targeted breeding for grain yield improvement.

Deciphering the evolution and complexity of wheat germplasm from a genomic perspective

Bread wheat provides an essential fraction of the daily calorific intake for humanity. Due to its huge and complex genome, progress in studying on the wheat genome is substantially trailed behind those of the other two major crops, rice and maize, for at least a decade. With rapid advances in genome assembling and reduced cost of high-throughput sequencing, emerging de novo genome assemblies of wheat and whole-genome sequencing data are leading to a paradigm shift in wheat research. Here, we review recent progress in dissecting the complex genome and germplasm evolution of wheat since the release of the first high-quality wheat genome. New insights have been gained in the evolution of wheat germplasm during domestication and modern breeding progress, genomic variations at multiple scales contributing to the diversity of wheat germplasm, and complex transcriptional and epigenetic regulations of functional genes in polyploid wheat. Genomics databases and bioinformatics tools meeting the urgent needs of wheat genomics research are also summarized. The ever-increasing omics data, along with advanced tools and well-structured databases, are expected to accelerate deciphering the germplasm and gene resources in wheat for future breeding advances.

Multi-omics atlas of combinatorial abiotic stress responses in wheat

Field-grown crops rarely experience growth conditions in which yield can be maximized. Environmental stresses occur in combination, with advancements in crop tolerance further complicated by its polygenic nature. Strategic targeting of causal genes is required to meet future crop production needs. Here, we employed a systems biology approach in wheat (Triticum aestivum L.) to investigate physio-metabolic adjustments and transcriptome reprogramming involved in acclimations to heat, drought, salinity and all combinations therein. A significant shift in magnitude and complexity of plant response was evident across stress scenarios based on the agronomic losses, increased proline concentrations and 8.7-fold increase in unique differentially expressed transcripts (DETs) observed under the triple stress condition. Transcriptome data from all stress treatments were assembled into an online, open access eFP browser for visualizing gene expression during abiotic stress. Weighted gene co-expression network analysis revealed 152 hub genes of which 32% contained the ethylene-responsive element binding factor-associated amphiphilic repression (EAR) transcriptional repression motif. Cross-referencing against the 31 DETs common to all stress treatments isolated TaWRKY33 as a leading candidate for greater plant tolerance to combinatorial stresses. Integration of our findings with available literature on gene functional characterization allowed us to further suggest flexible gene combinations for future adaptive gene stacking in wheat. Our approach demonstrates the strength of robust multi-omics-based data resources for gene discovery in complex environmental conditions. Accessibility of such datasets will promote cross-validation of candidate genes across studies and aid in accelerating causal gene validation for crop resiliency.

CoreDetector: a flexible and efficient program for core-genome alignment of evolutionary diverse genomes

MOTIVATION

Whole genome alignment of eukaryote species remains an important method for the determination of sequence and structural variations and can also be used to ascertain the representative non-redundant core-genome sequence of a population. Many whole genome alignment tools were first developed for the more mature analysis of prokaryote species with few current tools containing the functionality to process larger genomes of eukaryotes as well as genomes of more divergent species. In addition, the functionality of these tools becomes computationally prohibitive due to the significant compute resources needed to handle larger genomes.

RESULTS

In this research, we present CoreDetector, an easy-to-use general-purpose program that can align the core-genome sequences for a range of genome sizes and divergence levels. To illustrate the flexibility of CoreDetector, we conducted alignments of a large set of closely related fungal pathogen and hexaploid wheat cultivar genomes as well as more divergent fly and rodent species genomes. In all cases, compared to existing multiple genome alignment tools, CoreDetector exhibited improved flexibility, efficiency, and competitive accuracy in tested cases.

AVAILABILITY AND IMPLEMENTATION

CoreDetector was developed in the cross platform, and easily deployable, Java language. A packaged pipeline is readily executable in a bash terminal without any external need for Perl or Python environments. Installation, example data, and usage instructions for CoreDetector are freely available from https://github.com/mfruzan/CoreDetector.

Genome evolution and initial breeding of the Triticeae grass Leymus chinensis dominating the Eurasian Steppe

Leymus chinensis, a dominant perennial grass in the Eurasian Steppe, is well known for its remarkable adaptability and forage quality. Hardly any breeding has been done on the grass, limiting its potential in ecological restoration and forage productivity. To enable genetic improvement of the untapped, important species, we obtained a 7.85-Gb high-quality genome of L. chinensis with a particularly long contig N50 (318.49 Mb). Its allotetraploid genome is estimated to originate 5.29 million years ago (MYA) from a cross between the Ns-subgenome relating to Psathyrostachys and the unknown Xm-subgenome. Multiple bursts of transposons during 0.433-1.842 MYA after genome allopolyploidization, which involved predominantly the Tekay and Angela of LTR retrotransposons, contributed to its genome expansion and complexity. With the genome resource available, we successfully developed a genetic transformation system as well as the gene-editing pipeline in L. chinensis. We knocked out the monocot-specific miR528 using CRISPR/Cas9, resulting in the improvement of yield-related traits with increases in the tiller number and growth rate. Our research provides valuable genomic resources for Triticeae evolutionary studies and presents a conceptual framework illustrating the utilization of genomic information and genome editing to accelerate the improvement of wild L. chinensis with features such as polyploidization and self-incompatibility.

Initiation of B-type starch granules in wheat endosperm requires the plastidial α-glucan phosphorylase PHS1

The plastidial α-glucan phosphorylase (PHS1) can elongate and degrade maltooligosaccharides (MOSs), but its exact physiological role in plants is poorly understood. Here, we discover a specialized role of PHS1 in establishing the unique bimodal characteristic of starch granules in wheat (Triticum spp.) endosperm. Wheat endosperm contains large A-type granules that initiate at early grain development and small B-type granules that initiate in later grain development. We demonstrate that PHS1 interacts with B-GRANULE CONTENT1 (BGC1), a carbohydrate-binding protein essential for normal B-type granule initiation. Mutants of tetraploid durum wheat (Triticum turgidum) deficient in all homoeologs of PHS1 had normal A-type granules but fewer and larger B-type granules. Grain size and starch content were not affected by the mutations. Further, by assessing granule numbers during grain development in the phs1 mutant and using a double mutant defective in both PHS1 and BGC1, we demonstrate that PHS1 is exclusively involved in B-type granule initiation. The total starch content and number of starch granules per chloroplast in leaves were not affected by loss of PHS1, suggesting that its role in granule initiation in wheat is limited to the endosperm. We therefore propose that the initiation of A- and B-type granules occurs via distinct biochemical mechanisms, where PHS1 plays an exclusive role in B-type granule initiation.

Genome-Wide Investigation and Expression Analysis of the Catalase Gene Family in Oat Plants (Avena sativa L.)

Through the degradation of reactive oxygen species (ROS), different antioxidant enzymes, such as catalase (CAT), defend organisms against oxidative stress. These enzymes are crucial to numerous biological functions, like plant development and defense against several biotic and abiotic stresses. However, despite the major economic importance of Avena sativa around the globe, little is known about the CAT gene's structure and organization in this crop. Thus, a genome-wide investigation of the CAT gene family in oat plants has been carried out to characterize the potential roles of those genes under different stressors. Bioinformatic approaches were used in this study to predict the AvCAT gene's structure, secondary and tertiary protein structures, physicochemical properties, phylogenetic tree, and expression profiling under diverse developmental and biological conditions. A local Saudi oat variety (AlShinen) was used in this work. Here, ten AvCAT genes that belong to three groups (Groups I-III) were identified. All identified CATs harbor the two conserved domains (pfam00199 and pfam06628), a heme-binding domain, and a catalase activity motif. Moreover, identified AvCAT proteins were located in different compartments in the cell, such as the peroxisome, mitochondrion, and cytoplasm. By analyzing their promoters, different cis-elements were identified as being related to plant development, maturation, and response to different environmental stresses. Gene expression analysis revealed that three different AvCAT genes belonging to three different subgroups showed noticeable modifications in response to various stresses, such as mannitol, salt, and ABA. As far as we know, this is the first report describing the genome-wide analysis of the oat catalase gene family, and these data will help further study the roles of catalase genes during stress responses, leading to crop improvement.

Identification and introgression of a novel leaf rust resistance gene from Thinopyrum intermedium chromosome 7Js into wheat

KEY MESSAGE

A novel leaf rust resistance locus located on a terminal segment (0-69.29 Mb) of Thinopyrum intermedium chromosome arm 7JsS has been introduced into wheat genome for disease resistance breeding. Xiaoyan 78829, a wheat-Thinopyrum intermedium partial amphiploid, exhibits excellent resistance to fungal diseases in wheat. To transfer its disease resistance to common wheat (Triticum aestivum), we previously developed a translocation line WTT26 using chromosome engineering. Disease evaluation showed that WTT26 was nearly immune to 14 common races of leaf rust pathogen (Puccinia triticina) and highly resistant to Ug99 race PTKST of stem rust pathogen (P. graminis f. sp. tritici) at the seedling stage. It also displayed high adult plant resistance to powdery mildew (caused by Blumeria graminis f. sp. tritici). Cytogenetic and molecular marker analysis revealed that WTT26 carried a T4BS·7JsS chromosome translocation. Once transferred into the susceptible wheat genetic background, chromosome 7JsS exhibited its resistance to leaf rust, indicating that the resistance locus was located on this alien chromosome. To enhance the usefulness of this locus in wheat breeding, we further developed several new translocation lines with small Th. intermedium segments using irradiation and developed 124 specific markers using specific-locus amplified fragment sequencing, which increased the marker density of chromosome 7JsS. Furthermore, a refined physical map of chromosome 7JsS was constructed with 74 specific markers, and six bins were thus arranged according to the co-occurrence of markers and alien chromosome segments. Combining data from specific marker amplification and resistance evaluation, we mapped a new leaf rust resistance locus in the 0-69.29 Mb region on chromosome 7JsS. The translocation lines carrying the new leaf rust resistance locus and its linked markers will contribute to wheat disease-resistance breeding.

Transcriptomic analysis of wheat reveals possible resistance mechanism mediated by Yr10 to stripe rust

Stripe rust, caused by Puccinia striiformis f. sp. tritici (Pst), is a catastrophic disease that threatens global wheat yield. Yr10 is a race-specific all-stage disease resistance gene in wheat. However, the resistance mechanism of Yr10 is poorly characterized. Therefore, to elucidate the potential molecular mechanism mediated by Yr10, transcriptomic sequencing was performed at 0, 18, and 48 h post-inoculation (hpi) of compatible wheat Avocet S (AvS) and incompatible near-isogenic line (NIL) AvS + Yr10 inoculated with Pst race CYR32. Respectively, 227, 208, and 4050 differentially expressed genes (DEGs) were identified at 0, 18, and 48 hpi between incompatible and compatible interaction. The response of Yr10 to stripe rust involved various processes and activities, as indicated by the results of Gene Ontology (GO) enrichment analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis. Specifically, the response included photosynthesis, defense response to fungus, metabolic processes related to salicylic acid (SA) and jasmonic acid (JA), and activities related to reactive oxygen species (ROS). Ten candidate genes were selected for qRT-PCR verification and the results showed that the transcriptomic data was reliable. Through the functional analysis of candidate genes by the virus-induced gene silencing (VIGS) system, it was found that the gene TaHPPD (4-hydroxyphenylpyruvate dioxygenase) negatively regulated the resistance of wheat to stripe rust by affecting SA signaling, pathogenesis-related (PR) gene expression, and ROS clearance. Our study provides insight into Yr10-mediated resistance in wheat.

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.

Genome-wide linkage mapping of root system architecture-related traits in common wheat (Triticum aestivum L.)

Identifying loci for root system architecture (RSA) traits and developing available markers are crucial for wheat breeding. In this study, RSA-related traits, including total root length (TRL), total root area (TRA), and number of root tips (NRT), were evaluated in the Doumai/Shi4185 recombinant inbred line (RIL) population under hydroponics. In addition, both the RILs and parents were genotyped using the wheat 90K single-nucleotide polymorphism (SNP) array. In total, two quantitative trait loci (QTLs) each for TRL (QTRL.caas-4A.1 and QTRL.caas-4A.2), TRA (QTRA.caas-4A and QTRA.caas-4D), and NRT (QNRT.caas-5B and QNRT.caas-5D) were identified and each explaining 5.94%-9.47%, 6.85%-7.10%, and 5.91%-10.16% phenotypic variances, respectively. Among these, QTRL.caas-4A.1 and QTRA.caas-4A overlapped with previous reports, while QTRL.caas-4A.2, QTRA.caas-4D, QNRT.caas-5B, and QNRT.caas-5D were novel. The favorable alleles of QTRL.caas-4A.1, QTRA.caas-4A, and QTRA.caas-5B were contributed by Doumai, whereas the favorable alleles of QTRL.caas-4A.2, QTRA.caas-4D, and QTRA.caas-5D originated from Shi 4185. Additionally, two competitive allele-specific PCR (KASP) markers, Kasp_4A_RL (QTRA.caas-4A) and Kasp_5D_RT (QNRT.caas-5D), were developed and validated in 165 wheat accessions. This study provides new loci and available KASP markers, accelerating wheat breeding for higher yields.

Developing adapted wheat lines with broad-spectrum resistance to stem rust: Introgression of Sr59 through backcrossing and selections based on genotyping-by-sequencing data

Control of stem rust, caused by Puccinia graminis f.sp. tritici, a highly destructive fungal disease of wheat, faces continuous challenges from emergence of new virulent races across wheat-growing continents. Using combinations of broad-spectrum resistance genes could impart durable stem rust resistance. This study attempted transfer of Sr59 resistance gene from line TA5094 (developed through CSph1bM-induced T2DS·2RL Robertsonian translocation conferring broad-spectrum resistance). Poor agronomic performance of line TA5094 necessitates Sr59 transfer to adapted genetic backgrounds and utility evaluations for wheat improvement. Based on combined stem rust seedling and molecular analyses, 2070 BC1F1 and 1230 BC2F1 plants were derived from backcrossing BAJ#1, KACHU#1, and REEDLING#1 with TA5094. Genotyping-by-sequencing (GBS) results revealed the physical positions of 15,116 SNPs on chromosome 2R. The adapted genotypes used for backcrossing were found not to possess broad-spectrum resistance to selected stem rust races, whereas Sr59-containing line TA5094 showed resistance to all races tested. Stem rust seedling assays combined with kompetitive allele-specific PCR (KASP) marker analysis successfully selected and generated the BC2F2 population, which contained the Sr59 gene, as confirmed by GBS. Early-generation data from backcrossing suggested deviations from the 3:1 segregation, suggesting that multiple genes may contribute to Sr59 resistance reactions. Using GBS marker data (40,584 SNPs in wheat chromosomes) to transfer the recurrent parent background to later-generation populations resulted in average genome recovery of 71.2% in BAJ#1*2/TA5094, 69.8% in KACHU#1*2/TA5094, and 70.5% in REEDLING#1*2/TA5094 populations. GBS data verified stable Sr59 introgression in BC2F2 populations, as evidenced by presence of the Ph1 locus and absence of the 50,936,209 bp deletion in CSph1bM. Combining phenotypic selections, stem rust seedling assays, KASP markers, and GBS data substantially accelerated transfer of broad-spectrum resistance into adapted genotypes. Thus, this study demonstrated that the Sr59 resistance gene can be introduced into elite genetic backgrounds to mitigate stem rust-related yield losses.

Assessment of Indian wheat germplasm for Septoria nodorum blotch and tan spot reveals new QTLs conferring resistance along with recessive alleles of Tsn1 and Snn3

The leaf blight diseases, Septoria nodorum blotch (SNB), and tan spot (TS) are emerging due to changing climatic conditions in the northern parts of India. We screened 296 bread wheat cultivars released in India over the past 20 years for seedling resistance against SNB (three experiments) and TS (two experiments). According to a genome-wide association study, six QTLs on chromosome arms 1BL, 2AS, 5BL, and 6BL were particularly significant for SNB across all three years, of which Q.CIM.snb.1BL, Q.CIM.snb.2AS1, Q.CIM.snb.2AS.2, and Q.CIM.snb.6BL appeared novel. In contrast, those on 5BS and 5BL may correspond to Snn3 and Tsn1, respectively. The allelic combination of tsn1/snn3 conferred resistance to SNB, whereas that of Tsn1/Snn3 conferred high susceptibility. As for TS, Tsn1 was the only stably significant locus identified in this panel. Several varieties like PBW 771, DBW 277, and HD 3319, were identified as highly resistant to both diseases that can be used in future wheat improvement programs as resistant donors.

Identification of genetic loci for early maturity in spring bread wheat using the association analysis and gene dissection

Background

Early maturity in spring bread wheat is highly desirable in the regions where it enables the plants to evade high temperatures and plant pathogens at the end of the growing season.

Methods

To reveal the genetic loci responsible for the maturity time association analysis was carried out based on phenotyping for an 11-year period and high-throughput SNP genotyping of a panel of the varieties contrasting for this trait. The expression of candidate genes was verified using qPCR. The association between the SNP markers and the trait was validated using the biparental F2:3 population.

Results

Our data showed that under long-day conditions, the period from seedling to maturity is mostly influenced by the time from heading to maturity, rather than the heading time. The QTLs associated with the trait were located on 2A, 3B, 4A, 5B, 7A and 7B chromosomes with the 7BL locus being the most significant and promising for its SNPs accelerated the maturity time by about 9 days. Gene dissection in this locus detected a number of candidates, the best being TraesCS7B02G391800 (bZIP9) and TraesCS7B02G412200 (photosystem II reaction center). The two genes are predominantly expressed in the flag leaf while flowering. The effect of the SNPs was verified in F2:3 population and confirmed the association of the 4A, 5B and 7BL loci with the maturity time.

Field phenotyping for African crops: overview and perspectives

Improvements in crop productivity are required to meet the dietary demands of the rapidly-increasing African population. The development of key staple crop cultivars that are high-yielding and resilient to biotic and abiotic stresses is essential. To contribute to this objective, high-throughput plant phenotyping approaches are important enablers for the African plant science community to measure complex quantitative phenotypes and to establish the genetic basis of agriculturally relevant traits. These advances will facilitate the screening of germplasm for optimum performance and adaptation to low-input agriculture and resource-constrained environments. Increasing the capacity to investigate plant function and structure through non-invasive technologies is an effective strategy to aid plant breeding and additionally may contribute to precision agriculture. However, despite the significant global advances in basic knowledge and sensor technology for plant phenotyping, Africa still lags behind in the development and implementation of these systems due to several practical, financial, geographical and political barriers. Currently, field phenotyping is mostly carried out by manual methods that are prone to error, costly, labor-intensive and may come with adverse economic implications. Therefore, improvements in advanced field phenotyping capabilities and appropriate implementation are key factors for success in modern breeding and agricultural monitoring. In this review, we provide an overview of the current state of field phenotyping and the challenges limiting its implementation in some African countries. We suggest that the lack of appropriate field phenotyping infrastructures is impeding the development of improved crop cultivars and will have a detrimental impact on the agricultural sector and on food security. We highlight the prospects for integrating emerging and advanced low-cost phenotyping technologies into breeding protocols and characterizing crop responses to environmental challenges in field experimentation. Finally, we explore strategies for overcoming the barriers and maximizing the full potential of emerging field phenotyping technologies in African agriculture. This review paper will open new windows and provide new perspectives for breeders and the entire plant science community in Africa.

Molecular characterization of stable QTL and putative candidate genes for grain zinc and iron concentrations in two related wheat populations

KEY MESSAGE

Major QTL for grain zinc and iron concentrations were identified on the long arm of chromosomes 2D and 6D. Gene-based KASP markers were developed for putative candidate genes TaIPK1-2D and TaNAS10-6D. Micronutrient malnutrition is one of the most common public health problems in the world. Biofortification, the most attractive and sustainable solution to surmount malnutrition requires the development of micronutrient enriched new crop cultivars. In this study, two recombinant inbred line (RIL) populations, ZM175/XY60 and ZM175/LX987, were used to identify QTL for grain zinc concentration (GZnC), grain iron concentration (GFeC) and thousand grain weight (TGW). Eight QTL for GZnC, six QTL for GFeC and five QTL for TGW were detected. Three QTL on chromosomes 2DL and 4BS and chromosome 6A showed pleiotropic effects on all three traits. The 4BS and 6A QTL also increased plant height and might be Rht-B1a and Rht25a, respectively. The 2DL locus within a suppressed recombination region was identified in both RIL populations and the favorable allele simultaneously increasing GZnC, GFeC and TGW was contributed by XY60 and LX987. A QTL on chromosome 6DL associated only with GZnC was detected in ZM175/XY60 and was validated in JD8/AK58 RILs using kompetitive allele-specific PCR (KASP) marker K_AX-110119937. Both the 2DL and 6DL QTL were new loci for GZnC. Based on gene annotations, sequence variations and expression profiles, the phytic acid biosynthesis gene TaIPK1-2D and nicotianamine synthase gene TaNAS10-6D were predicted as candidate genes. Their gene-based KASP markers were developed and validated in a cultivar panel of 343 wheat accessions. This study investigated the genetic basis of GZnC and GFeC and provided valuable candidate genes and markers for breeding Zn- and Fe-enriched wheat.

Efficient proteome-wide identification of transcription factors targeting Glu-1: A case study for functional validation of TaB3-2A1 in wheat

High-molecular-weight glutenin subunits (HMW-GS), a major component of seed storage proteins (SSP) in wheat, largely determine processing quality. HMW-GS encoded by GLU-1 loci are mainly controlled at the transcriptional level by interactions between cis-elements and transcription factors (TFs). We previously identified a conserved cis-regulatory module CCRM1-1 as the most essential cis-element for Glu-1 endosperm-specific high expression. However, the TFs targeting CCRM1-1 remained unknown. Here, we built the first DNA pull-down plus liquid chromatography-mass spectrometry platform in wheat and identified 31 TFs interacting with CCRM1-1. TaB3-2A1 as proof of concept was confirmed to bind to CCRM1-1 by yeast one hybrid and electrophoretic mobility shift assays. Transactivation experiments demonstrated that TaB3-2A1 repressed CCRM1-1-driven transcription activity. TaB3-2A1 overexpression significantly reduced HMW-GS and other SSP, but enhanced starch content. Transcriptome analyses confirmed that enhanced expression of TaB3-2A1 down-regulated SSP genes and up-regulated starch synthesis-related genes, such as TaAGPL3, TaAGPS2, TaGBSSI, TaSUS1 and TaSUS5, suggesting that it is an integrator modulating the balance of carbon and nitrogen metabolism. TaB3-2A1 also had significant effects on agronomic traits, including heading date, plant height and grain weight. We identified two major haplotypes of TaB3-2A1 and found that TaB3-2A1-Hap1 conferred lower seed protein content, but higher starch content, plant height and grain weight than TaB3-2A1-Hap2 and was subjected to positive selection in a panel of elite wheat cultivars. These findings provide a high-efficiency tool to detect TFs binding to targeted promoters, considerable gene resources for dissecting regulatory mechanisms underlying Glu-1 expression, and a useful gene for wheat improvement.

QTL Mapping Reveals Both All-Stage and Adult-Plant Resistance to Powdery Mildew in Chinese Elite Wheat Cultivars

Powdery mildew caused by Blumeria graminis f. sp. tritici is a threat to wheat production in China. Mapping quantitative trait loci (QTL) for resistance to powdery mildew and developing breeder-friendly markers are important initial steps in breeding resistant cultivars. An all-stage resistance gene and several QTL were identified using a population of 254 recombinant inbred lines developed from a Jingdong 8/Aikang 58 cross. The population was evaluated for powdery mildew resistance across six field environments over three consecutive growing seasons utilizing two different mixtures of B. graminis f. sp. tritici isolates, named #Bgt-HB and #Bgt-BJ. Using genotypic data obtained from the Wheat TraitBreed 50K single-nucleotide polymorphism array, seven stable QTL were identified on chromosome arms 1DL, 2AL, 2DS, 4DL, 5AL, 6BL.1, and 6BL.2. The QTL on 2AL conferred all-stage resistance to B. graminis f. sp. tritici race E20 in greenhouse tests and explained up to 52% of the phenotypic variance in field trials but was resistant only against #Bgt-HB. The gene involved in this QTL was predicted to be Pm4a based on genome location and gene sequence. QPmja.caas-1DL, QPmja.caas-4DL, and QPmja.caas-6BL.1 were identified as potentially new QTL for powdery mildew resistance. QPmja.caas-2DS and QPmja.caas-6BL.1 were effective against both B. graminis f. sp. tritici mixtures, indicating their probable broad-spectrum resistance. A Kompetitive allele-specific PCR marker closely linked to QPmja.caas-2DS was developed and validated in a panel of 286 wheat cultivars. Because both Jingdong 8 and Aikang 58 have been leading cultivars and breeding parents, the QTL and marker reported represent valuable resources for wheat researchers and breeders.

Grain yield improvement by genome editing of TaARF12 that decoupled peduncle and rachis development trajectories via differential regulation of gibberellin signalling in wheat

Plant breeding is constrained by trade-offs among different agronomic traits by the pleiotropic nature of many genes. Genes that contribute to two or more favourable traits with no penalty on yield are rarely reported, especially in wheat. Here, we describe the editing of a wheat auxin response factor TaARF12 by using CRISPR/Cas9 that rendered shorter plant height with larger spikes. Changes in plant architecture enhanced grain number per spike up to 14.7% with significantly higher thousand-grain weight and up to 11.1% of yield increase under field trials. Weighted Gene Co-Expression Network Analysis (WGCNA) of spatial-temporal transcriptome profiles revealed two hub genes: RhtL1, a DELLA domain-free Rht-1 paralog, which was up-regulated in peduncle, and TaNGR5, an organ size regulator that was up-regulated in rachis, in taarf12 plants. The up-regulation of RhtL1 in peduncle suggested the repression of GA signalling, whereas up-regulation of TaNGR5 in spike may promote GA response, a working model supported by differential expression patterns of GA biogenesis genes in the two tissues. Thus, TaARF12 complemented plant height reduction with larger spikes that gave higher grain yield. Manipulation of TaARF12 may represent a new strategy in trait pyramiding for yield improvement in wheat.

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.

Increasing amyloplast size in wheat endosperm through mutation of PARC6 affects starch granule morphology

The determination of starch granule morphology in plants is poorly understood. The amyloplasts of wheat endosperm contain large discoid A-type granules and small spherical B-type granules. To study the influence of amyloplast structure on these distinct morphological types, we isolated a mutant in durum wheat (Triticum turgidum) defective in the plastid division protein PARC6, which had giant plastids in both leaves and endosperm. Endosperm amyloplasts of the mutant contained more A- and B-type granules than those of the wild-type. The mutant had increased A- and B-type granule size in mature grains, and its A-type granules had a highly aberrant, lobed surface. This morphological defect was already evident at early stages of grain development and occurred without alterations in polymer structure and composition. Plant growth and grain size, number and starch content were not affected in the mutants despite the large plastid size. Interestingly, mutation of the PARC6 paralog, ARC6, did not increase plastid or starch granule size. We suggest TtPARC6 can complement disrupted TtARC6 function by interacting with PDV2, the outer plastid envelope protein that typically interacts with ARC6 to promote plastid division. We therefore reveal an important role of amyloplast structure in starch granule morphogenesis in wheat.

The dynamics of touch-responsive gene expression in cereals

Wind, rain, herbivores, obstacles, neighbouring plants, etc. provide important mechanical cues to steer plant growth and survival. Mechanostimulation to stimulate yield and stress resistance of crops is of significant research interest, yet a molecular understanding of transcriptional responses to touch is largely absent in cereals. To address this, we performed whole-genome transcriptomics following mechanostimulation of wheat, barley, and the recent genome-sequenced oat. The largest transcriptome changes occurred ±25 min after touching, with most of the genes being upregulated. While most genes returned to basal expression level by 1-2 h in oat, many genes retained high expression even 4 h post-treatment in barley and wheat. Functional categories such as transcription factors, kinases, phytohormones, and Ca2+ regulation were affected. In addition, cell wall-related genes involved in (hemi)cellulose, lignin, suberin, and callose biosynthesis were touch-responsive, providing molecular insight into mechanically induced changes in cell wall composition. Furthermore, several cereal-specific transcriptomic footprints were identified that were not observed in Arabidopsis. In oat and barley, we found evidence for systemic spreading of touch-induced signalling. Finally, we provide evidence that both the jasmonic acid-dependent and the jasmonic acid-independent pathways underlie touch-signalling in cereals, providing a detailed framework and marker genes for further study of (a)biotic stress responses in cereals.

Identification of a Novel Pm65 Allele Conferring a Wide Spectrum of Resistance to Powdery Mildew in Wheat Accession PI 351817

Powdery mildew is caused by the highly adaptive biotrophic fungus Blumeria graminis f. sp. tritici infecting wheat worldwide. Novel powdery mildew resistance genes are urgently needed that can be used rapidly in wheat cultivar development with minimal disruption of trait advances elsewhere. PI 351817 is a German cultivar exhibiting a wide spectrum of resistance to B. graminis f. sp. tritici isolates collected from different wheat-growing regions of the United States. Evaluation of an F2 population and 237 F2:3 lines derived from OK1059060-2C14 × PI 351817 for responses to B. graminis f. sp. tritici isolate OKS(14)-B-3-1 identified a single dominant gene, designated Pm351817, for powdery mildew resistance in PI 351817. Using bulked segregant analysis (BSA) and simple sequence repeat (SSR) markers, Pm351817 was mapped in the terminal region of the long arm of chromosome 2A. Deep sequencing of the genotyping-by-sequencing libraries of the two parental lines identified a set of single-nucleotide polymorphism (SNP) markers in the 2AL candidate gene region. Those SNP markers was subsequently converted to Kompetitive allele-specific PCR (KASP) markers for genotyping the mapping population. Linkage analysis delimited Pm351817 to a 634-kb interval between Stars-KASP656 (771,207,512 bp) and Stars-KASP662 (771,841,609 bp) on 2AL, based on the Chinese Spring reference sequence IWGSC RefSeq v 2.1. Tests of allelism indicated that Pm351817 is located at the Pm65 locus. Pm351817 shows resistance to all B. graminis f. sp. tritici isolates used in this study and can be used to enhance powdery mildew resistance in the United States. KASP markers flanking Pm351817 can be used to select Pm351817 in wheat breeding programs after further tests for polymorphism.

The meso-octoploid Heliophila variabilis genome sheds a new light on the impact of polyploidization and diploidization on the diversity of the Cape flora

Although the South African Cape flora is one of the most remarkable biodiversity hotspots, its high diversity has not been associated with polyploidy. Here, we report the chromosome-scale genome assembly of an ephemeral cruciferous species Heliophila variabilis (~334 Mb, n = 11) adapted to South African semiarid biomes. Two pairs of differently fractionated subgenomes suggest an allo-octoploid origin of the genome at least 12 million years ago. The ancestral octoploid Heliophila genome (2n = 8x = ~60) has probably originated through hybridization between two allotetraploids (2n = 4x = ~30) formed by distant, intertribal, hybridization. Rediploidization of the ancestral genome was marked by extensive reorganization of parental subgenomes, genome downsizing, and speciation events in the genus Heliophila. We found evidence for loss-of-function changes in genes associated with leaf development and early flowering, and over-retention and sub/neofunctionalization of genes involved in pathogen response and chemical defense. The genomic resources of H. variabilis will help elucidate the role of polyploidization and genome diploidization in plant adaptation to hot arid environments and origin of the Cape flora. The sequenced H. variabilis represents the first chromosome-scale genome assembly of a meso-octoploid representative of the mustard family.

Multi-scale analysis provides insights into the roles of ureide permeases in wheat nitrogen use efficiency

The ureides allantoin and allantoate serve as nitrogen (N) transport compounds in plants, and more recently, allantoin has been shown to play a role in signaling. In planta, tissue ureide levels are controlled by the activity of enzymes of the purine degradation pathway and by ureide transporters called ureide permeases (UPS). Little is known about the physiological function of UPS proteins in crop plants, and especially in monocotyledon species. Here, we identified 13 TaUPS genes in the wheat (Triticum aestivum L.) genome. Phylogenetic and genome location analyses revealed a close relationship of wheat UPSs to orthologues in other grasses and a division into TaUPS1, TaUPS2.1, and TaUPS2.2 groups, each consisting of three homeologs, with a total of four tandem duplications. Expression, localization, and biochemical analyses resolved spatio-temporal expression patterns of TaUPS genes, transporter localization at the plasma membrane, and a role for TaUPS2.1 proteins in cellular import of ureides and phloem and seed loading. In addition, positive correlations between TaUPS1 and TaUPS2.1 transcripts and ureide levels were found. Together the data support that TaUPSs function in regulating ureide pools at source and sink, along with source-to-sink transport. Moreover, comparative studies between wheat cultivars grown at low and high N strengthened a role for TaUPS1 and TaUPS2.1 transporters in efficient N use and in controlling primary metabolism. Co-expression, protein-protein interaction, and haplotype analyses further support TaUPS involvement in N partitioning, N use efficiency, and domestication. Overall, this work provides a new understanding on UPS transporters in grasses as well as insights for breeding resilient wheat varieties with improved N use efficiency.

Cloning of the wheat leaf rust resistance gene Lr47 introgressed from Aegilops speltoides

Leaf rust, caused by Puccinia triticina Eriksson (Pt), is one of the most severe foliar diseases of wheat. Breeding for leaf rust resistance is a practical and sustainable method to control this devastating disease. Here, we report the identification of Lr47, a broadly effective leaf rust resistance gene introgressed into wheat from Aegilops speltoides. Lr47 encodes a coiled-coil nucleotide-binding leucine-rich repeat protein that is both necessary and sufficient to confer Pt resistance, as demonstrated by loss-of-function mutations and transgenic complementation. Lr47 introgression lines with no or reduced linkage drag are generated using the Pairing homoeologous1 mutation, and a diagnostic molecular marker for Lr47 is developed. The coiled-coil domain of the Lr47 protein is unable to induce cell death, nor does it have self-protein interaction. The cloning of Lr47 expands the number of leaf rust resistance genes that can be incorporated into multigene transgenic cassettes to control this devastating disease.

Thermosensitive SUMOylation of TaHsfA1 defines a dynamic ON/OFF molecular switch for the heat stress response in wheat

Dissecting genetic components in crop plants associated with heat stress (HS) sensing and adaptation will facilitate the design of modern crop varieties with improved thermotolerance. However, the molecular mechanisms underlying the ON/OFF switch controlling HS responses (HSRs) in wheat (Triticum aestivum) remain largely unknown. In this study, we focused on the molecular action of TaHsfA1, a class A heat shock transcription factor, in sensing dynamically changing HS signals and regulating HSRs. We show that the TaHsfA1 protein is modified by small ubiquitin-related modifier (SUMO) and that this modification is essential for the full transcriptional activation activity of TaHsfA1 in triggering downstream gene expression. During sustained heat exposure, the SUMOylation of TaHsfA1 is suppressed, which partially reduces TaHsfA1 protein activity, thereby reducing the intensity of downstream HSRs. In addition, we demonstrate that TaHsfA1 interacts with the histone acetyltransferase TaHAG1 in a thermosensitive manner. Together, our findings emphasize the importance of TaHsfA1 in thermotolerance in wheat. In addition, they define a highly dynamic SUMOylation-dependent "ON/OFF" molecular switch that senses temperature signals and contributes to thermotolerance in crops.

Characterization and fine mapping analysis of a major stable QTL qKnps-4A for kernel number per spike in wheat

KEY MESSAGE

A major stable QTL for kernel number per spike was narrowed down to a 2.19-Mb region containing two potential candidate genes, and its effects on yield-related traits were characterized. Kernel number per spike (KNPS) in wheat is a key yield component. Dissection and characterization of major stable quantitative trait loci (QTLs) for KNPS would be of considerable value for the genetic improvement of yield potential using molecular breeding technology. We had previously reported a major stable QTL controlling KNPS, qKnps-4A. In the current study, primary fine-mapping analysis, based on the primary mapping population, located qKnps-4A to an interval of approximately 6.8-Mb from 649.0 to 655.8 Mb on chromosome 4A refering to 'Kenong 9204' genome. Further fine-mapping analysis based on a secondary mapping population narrowed qKnps-4A to an approximately 2.19-Mb interval from 653.72 to 655.91 Mb. Transcriptome sequencing, gene function annotation analysis and homologous gene related reports showed that TraesKN4A01HG38570 and TraesKN4A01HG38590 were most likely to be candidate genes of qKnps-4A. Phenotypic analysis based on paired near-isogenic lines in the target region showed that qKnps-4A increased KNPS mainly by increasing the number of central florets per spike. We also evaluated the effects of qKnps-4A on other yield-related traits. Moreover, we dissected the QTL cluster of qKnps-4A and qTkw-4A and proved that the phenotypic effects were probably due to close linkage of two or more genes rather than pleiotropic effects of a single gene. This study provides molecular marker resource for wheat molecular breeding designed to improve yield potential, and lay the foundation for gene functional analysis of qKnps-4A.

Exploring GWAS and genomic prediction to improve Septoria tritici blotch resistance in wheat

Septoria tritici blotch (STB) is a destructive foliar diseases threatening wheat grain yield. Wheat breeding for STB disease resistance has been identified as the most sustainable and environment-friendly approach. In this work, a panel of 316 winter wheat breeding lines from a commercial breeding program were evaluated for STB resistance at the seedling stage under controlled conditions followed by genome-wide association study (GWAS) and genomic prediction (GP). The study revealed a significant genotypic variation for STB seedling resistance, while disease severity scores exhibited a normal frequency distribution. Moreover, we calculated a broad-sense heritability of 0.62 for the trait. Nine single- and multi-locus GWAS models identified 24 marker-trait associations grouped into 20 quantitative trait loci (QTLs) for STB seedling-stage resistance. The seven QTLs located on chromosomes 1B, 2A, 2B, 5B (two), 7A, and 7D are reported for the first time and could potentially be novel. The GP cross-validation analysis in the RR-BLUP model estimated the genomic-estimated breeding values (GEBVs) of STB resistance with a prediction accuracy of 0.49. Meanwhile, the GWAS assisted wRR-BLUP model improved the accuracy to 0.58. The identified QTLs can be used for marker-assisted backcrossing against STB in winter wheat. Moreover, the higher prediction accuracy recorded from the GWAS-assisted GP analysis implies its power to successfully select superior candidate lines based on their GEBVs for STB resistance.

Wheat powdery mildew resistance: from gene identification to immunity deployment

Powdery mildew is one of the most devastating diseases on wheat and is caused by the obligate biotrophic phytopathogen Blumeria graminis f. sp. tritici (Bgt). Due to the complexity of the large genome of wheat and its close relatives, the identification of powdery mildew resistance genes had been hampered for a long time until recent progress in large-scale sequencing, genomics, and rapid gene isolation techniques. Here, we describe and summarize the current advances in wheat powdery mildew resistance, emphasizing the most recent discoveries about the identification of genes conferring powdery mildew resistance and the similarity, diversity and molecular function of those genes. Multilayered resistance to powdery mildew in wheat could be used for counteracting Bgt, including durable, broad spectrum but partial resistance, as well as race-specific and mostly complete resistance mediated by nucleotide-binding and leucine rich repeat domain (NLR) proteins. In addition to the above mentioned layers, manipulation of susceptibility (S) and negative regulator genes may represent another layer that can be used for durable and broad-spectrum resistance in wheat. We propose that it is promising to develop effective and durable strategies to combat powdery mildew in wheat by simultaneous deployment of multilayered immunity.

Multi-target genome editing reduces polyphenol oxidase activity in wheat (Triticum aestivum L.) grains

Introduction

Polyphenol oxidases (PPO) are dual activity metalloenzymes that catalyse the production of quinones. In plants, PPO activity may contribute to biotic stress resistance and secondary metabolism but is undesirable for food producers because it causes the discolouration and changes in flavour profiles of products during post-harvest processing. In wheat (Triticum aestivum L.), PPO released from the aleurone layer of the grain during milling results in the discolouration of flour, dough, and end-use products, reducing their value. Loss-of-function mutations in the PPO1 and PPO2 paralogous genes on homoeologous group 2 chromosomes confer reduced PPO activity in the wheat grain. However, limited natural variation and the proximity of these genes complicates the selection of extremely low-PPO wheat varieties by recombination. The goal of the current study was to edit all copies of PPO1 and PPO2 to drive extreme reductions in PPO grain activity in elite wheat varieties.

Results

A CRISPR/Cas9 construct with one single guide RNA (sgRNA) targeting a conserved copper binding domain was used to edit all seven PPO1 and PPO2 genes in the spring wheat cultivar 'Fielder'. Five of the seven edited T1 lines exhibited significant reductions in PPO activity, and T2 lines had PPO activity up to 86.7% lower than wild-type. The same construct was transformed into the elite winter wheat cultivars 'Guardian' and 'Steamboat', which have five PPO1 and PPO2 genes. In these varieties PPO activity was reduced by >90% in both T1 and T2 lines. In all three varieties, dough samples from edited lines exhibited reduced browning.

Discussion

This study demonstrates that multi-target editing at late stages of variety development could complement selection for beneficial alleles in crop breeding programs by inducing novel variation in loci inaccessible to recombination.

Lifting of the 1,000 wheat exome project SNPs from Triticum aestivum cv. Chinese Spring assembly RefSeq v1.0 to RefSeq v2.1

OBJECTIVE

The 1,000 wheat exome project captured the single nucleotide variants in the coding regions of a diverse set of 890 wheat accessions to analyse the contribution of introgression to adaptation of wheat. However, this highly useful single nucleotide polymorphism (SNP) dataset is based on RefSeq v1.0 of the International Wheat Genome Sequencing Consortium (IWGSC) assembly of the bread wheat genome of Chinese Spring. This reference sequence has recently been updated using optical maps and long-read sequencing to produce the improved RefSeq v2.1. Our objective was to develop a reliable high-density SNP dataset positioned onto RefSeq v2.1 because it is the current standard reference sequence used by wheat researchers.

RESULTS

The 3,039,822 SNPs originally positioned on RefSeq v1.0 were projected to v2.1 using Liftoff with four different flanking regions, and 2,946,536 SNPs were consistently lifted to the same location irrespective of the flanking region lengths. Of these, 2,799,166 were located on the '+' ve strand. The distribution of the SNPs across the 21 chromosomes on RefSeq v2.1 was similar to that of RefSeq v1.0. Among the SNPs that were based on unanchored scaffolds in RefSeq v1.0, 11,938 were projected to one of the 21 pseudomolecules in the upgraded assembly. This SNP dataset constitutes a much-needed standardized resource for the wheat research community.

LysM Proteins TaCEBiP and TaLYK5 are Involved in Immune Responses Mediated by Chitin Coreceptor TaCERK1 in Wheat

Plant lysin motif (LysM) ectodomain receptors interact with pathogen-associated molecular patterns (PAMPs) and have critical functions in plant-microbe interactions. In this study, 65 LysM family genes were identified using the recent version of the reference sequence of bread wheat (Triticum aestivum), in which 23, 16, 20, and 6 members belonged to LysM-containing receptor-like kinases (LYKs), LysM-containing receptor-like proteins (LYPs), extracellular LysM proteins (LysMes), and intracellular nonsecretory LysM proteins (LysMns), respectively. The study found that TaCEBiP, TaLYK5, and TaCERK1 were highly responsive to PAMP elicitors and phytopathogens, with TaCEBiP and TaLYK5 binding directly to chitin. TaCERK1 acted as a coreceptor with TaCEBiP and TaLYK5 at the plasma membrane. Overexpression of TaCEBiP, TaLYK5, and TaCERK1 in Nicotiana benthamiana leaves exhibited enhanced resistance to Sclerotinia sclerotiorum. Subsequently, knocking down TaCEBiP, TaLYK5, and TaCERK1 genes with barley stripe mosaic virus-VIGS compromised the wheat defense response to an avirulent strain of Puccinia striiformis. The study concluded that wheat has two synergistic chitin perception systems for detecting pathogen elicitors, with the activated CERK1 intracellular kinase domain leading to signaling transduction. This research provides valuable insights into the functional roles and regulatory mechanisms of wheat LysM members under biotic stress.

Wheat (Triticum aestivum) chromosome 6D harbours the broad spectrum common bunt resistance gene Bt11

KEY MESSAGE

A major QTL on chromosome 6DL corresponding to bunt resistance gene Bt11 was identified in four mapping populations generated through crosses with Bt11-carriers PI 166910 and M822123. Common bunt in wheat has witnessed a renaissance with the rise of organic agriculture that began in the 1980s. The abandonment of systemic fungicides in organic farming, together with a lack of resistant cultivars, has led to wide-spread problems due to common bunt infections. Knowledge about genetic sources for resistance is still scarce and only few of the known bunt resistance factors are currently used in breeding. We therefore aimed to map the resistance factor harboured by the Turkish landrace PI 166910, which is the resistance donor for the Bt11 bunt differential line. Four mapping populations (MPs) with 96-132 recombinant inbred lines (RILs) were phenotyped for common bunt resistance over 2, 3 or 4 years with one or two local bunt populations and genotyped with the 25K SNP array. A major bunt resistance locus on the distal end of chromosome 6D designated QBt.ifa-6DL was identified in all MPs and experiments. Additional QTL contributing to resistance were detected on chromosomes 4B, 1A, 1B, 2A and 7B. QBt.ifa-6DL mapped to a region overlapping with the Bt9-locus identified in previous studies, but results indicate that QBt.ifa-6DL is different from Bt9 and convincing evidence from haplotype comparisons suggests that it represents the Bt11 resistance allele. Markers for the distal region of chromosome 6D between 492.6 and 495.2 Mbp can be used to select for QBt.ifa-6DL. This resistance factor confers high and stable resistance against common bunt and should be integrated into organic and low-input wheat breeding programs.