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

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.

Plant and pathogen genomics: essential approaches for stem rust resistance gene stacks in wheat

The deployment of disease resistance genes is currently the most economical and environmentally sustainable method of crop protection. However, disease resistance genes can rapidly break down because of constant pathogen evolution, particularly when they are deployed singularly. Polygenic resistance is, therefore, considered the most durable, but combining and maintaining these genes by breeding is a laborious process as effective genes are usually unlinked. The deployment of polygenic resistance with single-locus inheritance is a promising innovation that overcomes these difficulties while enhancing resistance durability. Because of major advances in genomic technologies, increasing numbers of plant resistance genes have been cloned, enabling the development of resistance transgene stacks (RTGSs) that encode multiple genes all located at a single genetic locus. Gene stacks encoding five stem rust resistance genes have now been developed in transgenic wheat and offer both breeding simplicity and potential resistance durability. The development of similar genomic resources in phytopathogens has advanced effector gene isolation and, in some instances, enabled functional validation of individual resistance genes in RTGS. Here, the wheat stem rust pathosystem is used as an illustrative example of how host and pathogen genomic advances have been instrumental in the development of RTGS, which is a strategy applicable to many other agricultural crop species.

Post-hybridization introgression and natural selection promoted genomic divergence of Aegilops speltoides and the four S*-genome diploid species

Understanding how different driving forces have promoted biological divergence and speciation is one of the central issues in evolutionary biology. The Triticum/Aegilops species complex contains 13 diploid species belonging to the A-, B- and D-lineages and offers an ideal system to address the evolutionary dynamics of lineage fusion and splitting. Here, we sequenced the whole genomes of one S-genome species (Aegilops speltoides) of the B-lineage and four S*-genome diploid species (Aegilops bicornis, Aegilops longissima, Aegilops sharonensis and Aegilops searsii) of the D-lineage at the population level. We performed detailed comparisons of the five species and with the other four representative A-, B- and D-lineage species. Our estimates identified frequent genetic introgressions from A- and B-lineages to the D-lineage species. A remarkable observation is the contrasting distributions of putative introgressed loci by the A- and B-lineages along all the seven chromosomes to the extant D-lineage species. These genetic introgressions resulted in high levels of genetic divergence at centromeric regions between Ae. speltoides (B-lineage) and the other four S*-genome diploid species (D-lineage), while natural selection is a potential contributor to divergence among the four S*-genome species at telomeric regions. Our study provides a genome-wide view on how genetic introgression and natural selection acted together yet chromosome-regionally divided to promote genomic divergence among the five S- and S*-genome diploid species, which provides new and nuanced insights into the evolutionary history of the Triticum/Aegilops species complex.

A low-coverage 3' RNA-seq to detect homeolog expression in polyploid wheat

Although allopolyploid species are common among natural and crop species, it is not easy to distinguish duplicated genes, known as homeologs, during their genomic analysis. Yet, cost-efficient RNA sequencing (RNA-seq) is to be developed for large-scale transcriptomic studies such as time-series analysis and genome-wide association studies in allopolyploids. In this study, we employed a 3' RNA-seq utilizing 3' untranslated regions (UTRs) containing frequent mutations among homeologous genes, compared to coding sequence. Among the 3' RNA-seq protocols, we examined a low-cost method Lasy-Seq using an allohexaploid bread wheat, Triticum aestivum. HISAT2 showed the best performance for 3' RNA-seq with the least mapping errors and quick computational time. The number of detected homeologs was further improved by extending 1 kb of the 3' UTR annotation. Differentially expressed genes in response to mild cold treatment detected by the 3' RNA-seq were verified with high-coverage conventional RNA-seq, although the latter detected more differentially expressed genes. Finally, downsampling showed that even a 2 million sequencing depth can still detect more than half of expressed homeologs identifiable by the conventional 32 million reads. These data demonstrate that this low-cost 3' RNA-seq facilitates large-scale transcriptomic studies of allohexaploid wheat and indicate the potential application to other allopolyploid species.

All families of transposable elements were active in the recent wheat genome evolution and polyploidy had no impact on their activity

Bread wheat (Triticum aestivum L.) is a major crop and its genome is one of the largest ever assembled at reference-quality level. It is 15 Gb, hexaploid, with 85% of transposable elements (TEs). Wheat genetic diversity was mainly focused on genes and little is known about the extent of genomic variability affecting TEs, transposition rate, and the impact of polyploidy. Multiple chromosome-scale assemblies are now available for bread wheat and for its tetraploid and diploid wild relatives. In this study, we computed base pair-resolved, gene-anchored, whole genome alignments of A, B, and D lineages at different ploidy levels in order to estimate the variability that affects the TE space. We used assembled genomes of 13 T. aestivum cultivars (6x = AABBDD) and a single genome for Triticum durum (4x = AABB), Triticum dicoccoides (4x = AABB), Triticum urartu (2x = AA), and Aegilops tauschii (2x = DD). We show that 5%-34% of the TE fraction is variable, depending on the species divergence. Between 400 and 13,000 novel TE insertions per subgenome were detected. We found lineage-specific insertions for nearly all TE families in di-, tetra-, and hexaploids. No burst of transposition was observed and polyploidization did not trigger any boost of transposition. This study challenges the prevailing idea of wheat TE dynamics and is more in agreement with an equilibrium model of evolution.

Comparative transcriptome analysis reveals the key genes and pathways involved in drought stress response of two wheat (Triticum aestivum L) varieties

Drought stress results in significant yield losses in wheat production. Although studies have reported a number of wheat drought tolerance genes, a deeper understanding of the tolerance mechanisms is required for improving wheat tolerance against drought stress. In this study, we found that "Deguo 2" exhibited higher tolerance to drought than "Truman". Transcriptomics analysis enabled identification of 6084 and 7146 differentially expressed genes (DEGs), mainly mapping flavonoid biosynthesis, plant hormone, phenolamides and antioxidant pathways and revealed altered expression levels of about 700 genes. Exogenous melatonin application enhanced wheat tolerance against drought stress. Co-expression analysis showed that bHLH and bZIP transcription factors may be involved in the regulation of various pathway genes. Take together, these results provide new insights for us on exploring the crosstalk between phytohormones and secondary metabolites, and will deepen the understanding of the complex tolerance mechanisms against drought stress in wheat.

TaTIP41 and TaTAP46 positively regulate drought tolerance in wheat by inhibiting PP2A activity

Drought is a major environmental stress limiting global wheat (Triticum aestivum) production. Exploring drought tolerance genes is important for improving drought adaptation in this crop. Here, we cloned and characterized TaTIP41, a novel drought tolerance gene in wheat. TaTIP41 is a putative conserved component of target of rapamycin (TOR) signaling, and the TaTIP41 homoeologs were expressed in response to drought stress and abscisic acid (ABA). The overexpression of TaTIP41 enhanced drought tolerance and the ABA response, including ABA-induced stomatal closure, while its downregulation using RNA interference (RNAi) had the opposite effect. Furthermore, TaTIP41 physically interacted with TaTAP46, another conserved component of TOR signaling. Like TaTIP41, TaTAP46 positively regulated drought tolerance. Furthermore, TaTIP41 and TaTAP46 interacted with type-2A protein phosphatase (PP2A) catalytic subunits, such as TaPP2A-2, and inhibited their enzymatic activities. Silencing TaPP2A-2 improved drought tolerance in wheat. Together, our findings provide new insights into the roles of TaTIP41 and TaTAP46 in the drought tolerance and ABA response in wheat, and their potential application in improving wheat environmental adaptability.

Genome-wide identification of TaCIPK gene family members in wheat and their roles in host response to Blumeria graminis f. sp. tritici infection

Powdery mildew, caused by Blumeria graminis f. sp. tritici (Bgt), is a destructive disease affecting wheat crops worldwide. Functional genes can be activated in response to Bgt inoculations. Calcineurin B-like protein (CBL) together with CBL-interacting protein kinase (CIPK) forms the CBL-CIPK protein complex that participates in Ca2+ sensor kinase-related signaling pathways responding to abiotic and biotic stresses. In this study, we performed a genome-wide screening and identified 27 CIPK subfamilies (123 CIPK transcripts, TaCIPKs) including 55 new and 47 updated TaCIPKs in wheat. Phylogenetic analysis revealed that 123 TaCIPKs could be divided into four groups. Segmental duplications and tandem repeats promoted the expansion of the TaCIPK family. Gene function was further evidenced by differences in gene structure, cis-elements, and protein domains. TaCIPK15-4A was cloned in this study. TaCIPK15-4A contained 17 serine, seven tyrosine, and 15 threonine phosphorylation sites and localized in the plasma membrane and cytoplasm. TaCIPK15-4A expression was induced after Bgt inoculation. Virus-induced gene silencing and overexpression experiments indicated that TaCIPK15-4A could play a positive role in wheat disease resistance to Bgt. Overall, these results provide insights into the role of the TaCIPK gene family in wheat resistance and could be beneficial for further research to prevent Bgt infection.

Genome-wide circadian gating of a cold temperature response in bread wheat

Circadian rhythms coordinate the responses of organisms with their daily fluctuating environments, by establishing a temporal program of gene expression. This schedules aspects of metabolism, physiology, development and behaviour according to the time of day. Circadian regulation in plants is extremely pervasive, and is important because it underpins both productivity and seasonal reproduction. Circadian regulation extends to the control of environmental responses through a regulatory process known as circadian gating. Circadian gating is the process whereby the circadian clock regulates the response to an environmental cue, such that the magnitude of response to an identical cue varies according to the time of day of the cue. Here, we show that there is genome-wide circadian gating of responses to cold temperatures in plants. By using bread wheat as an experimental model, we establish that circadian gating is crucial to the programs of gene expression that underlie the environmental responses of a crop of major socioeconomic importance. Furthermore, we identify that circadian gating of cold temperature responses are distributed unevenly across the three wheat subgenomes, which might reflect the geographical origins of the ancestors of modern wheat.

Chromosome translocation affects multiple phenotypes, causes genome-wide dysregulation of gene expression, and remodels metabolome in hexaploid wheat

Chromosomal rearrangements (CRs) may occur in newly formed polyploids due to compromised meiotic fidelity. Moreover, CRs can be more readily tolerated in polyploids allowing their longer-term retention and hence potential spreading/fixation within a lineage. The direct functional consequences of CRs in plant polyploids remain unexplored. Here, we identified a heterozygous individual from a synthetic allohexaploid wheat in which the terminal parts of the long-arms of chromosomes 2D (approximately 193 Mb) and 4A (approximately 167 Mb) were reciprocally translocated. Five homogeneous translocation lines including both unbalanced and balanced types were developed by selfing fertilization of the founder mutant (RT [2DL; 4AL]-ter/1, reciprocal translocation). We investigated impacts of these translocations on phenotype, genome-wide gene expression and metabolome. We find that, compared with sibling wild-type, CRs in the form of both unbalanced and balanced translocations induced substantial changes of gene expression primarily via trans-regulation in the nascent allopolyploid wheat. The CRs also manifested clear phenotypic and metabolic consequences. In particular, the genetically balanced, stable reciprocal translocations lines showed immediate enhanced reproductive fitness relative to wild type. Our results underscore the profound impact of CRs on gene expression in nascent allopolyploids with wide-ranging phenotypic and metabolic consequences, suggesting CRs are an important source of genetic variation that can be exploited for crop breeding.

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.

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.

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.

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.