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.

A scalable phenotyping approach for female floral organ development and senescence in the absence of pollination in wheat.

In the absence of pollination, female reproductive organs senesce, leading to an irrevocable loss in the reproductive potential of the flower, which directly affects seed set. In self-pollinating crops like wheat (Triticum aestivum), the post-anthesis viability of unpollinated carpels has been overlooked, despite its importance for hybrid seed production systems. To advance our knowledge of carpel development in the absence of pollination, we created a high-throughput phenotyping approach to quantify stigma and ovary morphology. We demonstrate the suitability of the approach, which uses light-microscopy imaging and machine learning, for the analysis of floral organ traits in field-grown plants using fresh and fixed samples. We show that the unpollinated carpel undergoes a well-defined initial growth phase, followed by a peak phase in which stigma area reaches its maximum and the radial expansion of the ovary slows, and a final deterioration phase. These developmental dynamics were consistent across years and could be used to classify male-sterile cultivars. This phenotyping approach provides a new tool for examining carpel development, which we hope will advance research into female fertility of wheat.

Combining Traditional Mutagenesis with New High-Throughput Sequencing and Genome Editing to Reveal Hidden Variation in Polyploid Wheat.

Induced mutations have been used to generate novel variation for breeding purposes since the early 1900s. However, the combination of this old technology with the new capabilities of high-throughput sequencing has resulted in powerful reverse genetic approaches in polyploid crops. Sequencing genomes or exomes of large mutant populations can generate extensive databases of mutations for most genes. These mutant collections, together with genome editing, are being used in polyploid species to combine mutations in all copies of a gene (homoeologs), and to expose phenotypic variation that was previously hidden by functional redundancy among homoeologs. This redundancy is more extensive in recently formed polyploids such as wheat, which can now benefit from the deployment of useful recessive mutations previously identified in its diploid relatives. Sequenced mutant populations and genome editing have changed the paradigm of what is possible in functional genetic analysis of wheat.

Pathogen-induced biosynthetic pathways encode defense-related molecules in bread wheat.

Wheat is a widely grown food crop that suffers major yield losses due to attack by pests and pathogens. A better understanding of biotic stress responses in wheat is thus of major importance. The recently assembled bread wheat genome coupled with extensive transcriptomic resources provides unprecedented new opportunities to investigate responses to pathogen challenge. Here, we analyze gene coexpression networks to identify modules showing consistent induction in response to pathogen exposure. Within the top pathogen-induced modules, we identify multiple clusters of physically adjacent genes that correspond to six pathogen-induced biosynthetic pathways that share a common regulatory network. Functional analysis reveals that these pathways, all of which are encoded by biosynthetic gene clusters, produce various different classes of compounds­namely, flavonoids, diterpenes, and triterpenes, including the defense-related compound ellarinacin. Through comparative genomics, we also identify associations with the known rice phytoalexins momilactones, as well as with a defense-related gene cluster in the grass model plant Brachypodium distachyon. Our results significantly advance the understanding of chemical defenses in wheat and open up avenues for enhancing disease resistance in this agriculturally important crop. They also exemplify the power of transcriptional networks to discover the biosynthesis of chemical defenses in plants with large, complex genomes.

Root angle is controlled by EGT1 in cereal crops employing an antigravitropic mechanism.

Root angle in crops represents a key trait for efficient capture of soil resources. Root angle is determined by competing gravitropic versus antigravitropic offset (AGO) mechanisms. Here we report a root angle regulatory gene termed ENHANCED GRAVITROPISM1 (EGT1) that encodes a putative AGO component, whose loss-of-function enhances root gravitropism. Mutations in barley and wheat EGT1 genes confer a striking root phenotype, where every root class adopts a steeper growth angle. EGT1 encodes an F-box and Tubby domain-containing protein that is highly conserved across plant species. Haplotype analysis found that natural allelic variation at the barley EGT1 locus impacts root angle. Gravitropic assays indicated that Hvegt1 roots bend more rapidly than wild-type. Transcript profiling revealed Hvegt1 roots deregulate reactive oxygen species (ROS) homeostasis and cell wall-loosening enzymes and cofactors. ROS imaging shows that Hvegt1 root basal meristem and elongation zone tissues have reduced levels. Atomic force microscopy measurements detected elongating Hvegt1 root cortical cell walls are significantly less stiff than wild-type. In situ analysis identified HvEGT1 is expressed in elongating cortical and stele tissues, which are distinct from known root gravitropic perception and response tissues in the columella and epidermis, respectively. We propose that EGT1 controls root angle by regulating cell wall stiffness in elongating root cortical tissue, counteracting the gravitropic machinery's known ability to bend the root via its outermost tissues. We conclude that root angle is controlled by EGT1 in cereal crops employing an antigravitropic mechanism.

Crossover-active regions of the wheat genome are distinguished by DMC1, the chromosome axis, H3K27me3, and signatures of adaptation.

The hexaploid bread wheat genome comprises over 16 gigabases of sequence across 21 chromosomes. Meiotic crossovers are highly polarized along the chromosomes, with elevation in the gene-dense distal regions and suppression in the Gypsy retrotransposon-dense centromere-proximal regions. We profiled the genomic landscapes of the meiotic recombinase DMC1 and the chromosome axis protein ASY1 in wheat and investigated their relationships with crossovers, chromatin state, and genetic diversity. DMC1 and ASY1 chromatin immunoprecipitation followed by sequencing (ChIP-seq) revealed strong co-enrichment in the distal, crossover-active regions of the wheat chromosomes. Distal ChIP-seq enrichment is consistent with spatiotemporally biased cytological immunolocalization of DMC1 and ASY1 close to the telomeres during meiotic prophase I. DMC1 and ASY1 ChIP-seq peaks show significant overlap with genes and transposable elements in the Mariner and Mutator superfamilies. However, DMC1 and ASY1 ChIP-seq peaks were detected along the length of each chromosome, including in low-crossover regions. At the fine scale, crossover elevation at DMC1 and ASY1 peaks and genes correlates with enrichment of the Polycomb histone modification H3K27me3. This indicates a role for facultative heterochromatin, coincident with high DMC1 and ASY1, in promoting crossovers in wheat and is reflected in distalized H3K27me3 enrichment observed via ChIP-seq and immunocytology. Genes with elevated crossover rates and high DMC1 and ASY1 ChIP-seq signals are overrepresented for defense-response and immunity annotations, have higher sequence polymorphism, and exhibit signatures of selection. Our findings are consistent with meiotic recombination promoting genetic diversity, shaping host-pathogen co-evolution, and accelerating adaptation by increasing the efficiency of selection.

FIGL1 prevents aberrant chromosome associations and fragmentation and limits crossovers in polyploid wheat meiosis.

Meiotic crossovers (COs) generate genetic diversity and are crucial for viable gamete production. Plant COs are typically limited to 1-3 per chromosome pair, constraining the development of improved varieties, which in wheat is exacerbated by an extreme distal localisation bias. Advances in wheat genomics and related technologies provide new opportunities to investigate, and possibly modify, recombination in this important crop species. Here, we investigate the disruption of FIGL1 in tetraploid and hexaploid wheat as a potential strategy for modifying CO frequency/position. We analysed figl1 mutants and virus-induced gene silencing lines cytogenetically. Genetic mapping was performed in the hexaploid. FIGL1 prevents abnormal meiotic chromosome associations/fragmentation in both ploidies. It suppresses class II COs in the tetraploid such that CO/chiasma frequency increased 2.1-fold in a figl1 msh5 quadruple mutant compared with a msh5 double mutant. It does not appear to affect class I COs based on HEI10 foci counts in a hexaploid figl1 triple mutant. Genetic mapping in the triple mutant suggested no significant overall increase in total recombination across examined intervals but revealed large increases in specific individual intervals. Notably, the tetraploid figl1 double mutant was sterile but the hexaploid triple mutant was moderately fertile, indicating potential utility for wheat breeding.

A two-gene strategy increases iron and zinc concentrations in wheat flour, improving mineral bioaccessibility.

Dietary deficiencies of iron and zinc cause human malnutrition that can be mitigated by biofortified staple crops. Conventional breeding approaches to increase grain mineral concentrations in wheat (Triticum aestivum L.) have had only limited success, and our understanding of the genetic and physiological barriers to altering this trait is incomplete. Here we demonstrate that a transgenic approach combining endosperm-specific expression of the wheat VACUOLAR IRON TRANSPORTER gene TaVIT2-D with constitutive expression of the rice (Oryza sativa) NICOTIANAMINE SYNTHASE gene OsNAS2 significantly increases the total concentration of zinc and relocates iron to white-flour fractions. In two distinct bread wheat cultivars, we show that the so called VIT-NAS construct led to a two-fold increase in zinc in wholemeal flour, to ∼50 µg g-1. Total iron was not significantly increased, but redistribution within the grain resulted in a three-fold increase in iron in highly pure, roller-milled white flour, to ∼25 µg g-1. Interestingly, expression of OsNAS2 partially restored iron translocation to the aleurone, which is iron depleted in grain overexpressing TaVIT2 alone. A greater than three-fold increase in the level of the natural plant metal chelator nicotianamine in the grain of VIT-NAS lines corresponded with improved iron and zinc bioaccessibility in white flour. The growth of VIT-NAS plants in the greenhouse was indistinguishable from untransformed controls. Our results provide insights into mineral translocation and distribution in wheat grain and demonstrate that the individual and combined effects of the two transgenes can enhance the nutritional quality of wheat beyond what is possible by conventional breeding.

Ectopic expression of Triticum polonicum VRT-A2 underlies elongated glumes and grains in hexaploid wheat in a dosage-dependent manner.

Flower development is an important determinant of grain yield in crops. In wheat (Triticum spp.), natural variation for the size of spikelet and floral organs is particularly evident in Triticum turgidum ssp. polonicum (also termed Triticum polonicum), a tetraploid subspecies of wheat with long glumes, lemmas, and grains. Using map-based cloning, we identified VEGETATIVE TO REPRODUCTIVE TRANSITION 2 (VRT2), which encodes a MADS-box transcription factor belonging to the SHORT VEGETATIVE PHASE family, as the gene underlying the T. polonicum long-glume (P1) locus. The causal P1 mutation is a sequence rearrangement in intron-1 that results in ectopic expression of the T. polonicum VRT-A2 allele. Based on allelic variation studies, we propose that the intron-1 mutation in VRT-A2 is the unique T. polonicum subspecies-defining polymorphism, which was later introduced into hexaploid wheat via natural hybridizations. Near-isogenic lines differing for the P1 locus revealed a gradient effect of P1 across spikelets and within florets. Transgenic lines of hexaploid wheat carrying the T. polonicum VRT-A2 allele show that expression levels of VRT-A2 are highly correlated with spike, glume, grain, and floral organ length. These results highlight how changes in expression profiles, through variation in cis-regulation, can affect agronomic traits in a dosage-dependent manner in polyploid crops.

ENHANCED GRAVITROPISM 2 encodes a STERILE ALPHA MOTIF-containing protein that controls root growth angle in barley and wheat.

The root growth angle defines how roots grow toward the gravity vector and is among the most important determinants of root system architecture. It controls water uptake capacity, nutrient use efficiency, stress resilience, and, as a consequence, yield of crop plants. We demonstrated that the egt2 (enhanced gravitropism 2) mutant of barley exhibits steeper root growth of seminal and lateral roots and an auxin-independent higher responsiveness to gravity compared to wild-type plants. We cloned the EGT2 gene by a combination of bulked-segregant analysis and whole genome sequencing. Subsequent validation experiments by an independent CRISPR/Cas9 mutant allele demonstrated that egt2 encodes a STERILE ALPHA MOTIF domain-containing protein. In situ hybridization experiments illustrated that EGT2 is expressed from the root cap to the elongation zone. We demonstrated the evolutionary conserved role of EGT2 in root growth angle control between barley and wheat by knocking out the EGT2 orthologs in the A and B genomes of tetraploid durum wheat. By combining laser capture microdissection with RNA sequencing, we observed that seven expansin genes were transcriptionally down-regulated in the elongation zone. This is consistent with a role of EGT2 in this region of the root where the effect of gravity sensing is executed by differential cell elongation. Our findings suggest that EGT2 is an evolutionary conserved regulator of root growth angle in barley and wheat that could be a valuable target for root-based crop improvement strategies in cereals.

High expression of the MADS-box gene VRT2 increases the number of rudimentary basal spikelets in wheat.

Spikelets are the fundamental building blocks of Poaceae inflorescences, and their development and branching patterns determine the various inflorescence architectures and grain yield of grasses. In wheat (Triticum aestivum), the central spikelets produce the most and largest grains, while spikelet size gradually decreases acropetally and basipetally, giving rise to the characteristic lanceolate shape of wheat spikes. The acropetal gradient corresponds with the developmental age of spikelets; however, the basal spikelets are developed first, and the cause of their small size and rudimentary development is unclear. Here, we adapted G&T-seq, a low-input transcriptomics approach, to characterize gene expression profiles within spatial sections of individual spikes before and after the establishment of the lanceolate shape. We observed larger differences in gene expression profiles between the apical, central, and basal sections of a single spike than between any section belonging to consecutive developmental time points. We found that SHORT VEGETATIVE PHASE MADS-box transcription factors, including VEGETATIVE TO REPRODUCTIVE TRANSITION 2 (VRT-A2), are expressed highest in the basal section of the wheat spike and display the opposite expression gradient to flowering E-class SEPALLATA1 genes. Based on multi-year field trials and transgenic lines, we show that higher expression of VRT-A2 in the basal sections of the spike is associated with increased numbers of rudimentary basal spikelets. Our results, supported by computational modeling, suggest that the delayed transition of basal spikelets from vegetative to floral developmental programs results in the lanceolate shape of wheat spikes. This study highlights the value of spatially resolved transcriptomics to gain insights into developmental genetics pathways of grass inflorescences.

Delayed development of basal spikelets in wheat explains their increased floret abortion and rudimentary nature.

Large differences exist in the number of grains per spikelet across an individual wheat (Triticum aestivum L.) spike. The central spikelets produce the highest number of grains, while apical and basal spikelets are less productive, and the most basal spikelets are commonly only developed in rudimentary form. Basal spikelets are delayed in initiation, yet they continue to develop and produce florets. The precise timing or the cause of their abortion remains largely unknown. Here, we investigated the underlying causes of basal spikelet abortion using shading applications in the field. We found that basal spikelet abortion is likely to be the consequence of complete floret abortion, as both occur concurrently and have the same response to shading treatments. We detected no differences in assimilate availability across the spike. Instead, we show that the reduced developmental age of basal florets pre-anthesis is strongly associated with their increased abortion. Using the developmental age pre-abortion, we were able to predict final grain set per spikelet across the spike, alongside the characteristic gradient in the number of grains from basal to central spikelets. Future efforts to improve spikelet homogeneity across the spike could thus focus on improving basal spikelet establishment and increasing floret development rates pre-abortion.

Identification of a novel SNP in the miR172 binding site of Q homoeolog AP2L-D5 is associated with spike compactness and agronomic traits in wheat (Triticum aestivum L.).

KEY MESSAGE: This study found that the compact spike locus of ANK-15 is on chromosome 5D instead of 2B. We have identified a new allele of AP2L-D5 as the candidate causal polymorphism. Spike architecture is a key determinant of wheat yield, a crop which supports much of the human diet but whose yield gains are stagnating. Spike architecture mutants offer opportunities to identify genetic factors contributing to inflorescence development. Here, we investigate the locus underlying the compact spike phenotype of mutant line ANK-15 by conducting mRNA-sequencing and genetic mapping using ANK-15 and its non-compact spike near-isogenic line Novosibirskaya 67 (N67). Previous literature has placed the compact spike locus of ANK-15 to chromosome 2B. However, based on the single nucleotide polymorphisms (SNPs) identified using mRNA-seq data, we were unable to detect polymorphisms between N67 and ANK-15 in the putative chromosome 2B region. We performed differential expression analysis of developing rachis and found that AP2L-D5, the D homoeolog of the domestication Q gene, is upregulated in ANK-15 in comparison to N67. ANK-15 carries a SNP in the microRNA172 binding site of AP2L-D5, which is predicted to lead to higher expression of AP2L-D5 due to decreased miRNA172-mediated degradation. Furthermore, we performed genetic mapping using an ANK-15 × N67 F2 population and found a single quantitative trait locus on chromosome 5D coinciding with the position of AP2L-D5. This result suggests that AP2L-D5 is likely the underlying causal gene for the compact spike phenotype in ANK-15. We performed a field trial to investigate the effect of the AP2L-D5 allele on agronomic traits and found that the AP2L-D5 allele from ANK-15 is associated with a significant reduction in height, increased thousand grain weight (TGW), and increased grain width.

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.

Genomic innovation for crop improvement.

Crop production needs to increase to secure future food supplies, while reducing its impact on ecosystems. Detailed characterization of plant genomes and genetic diversity is crucial for meeting these challenges. Advances in genome sequencing and assembly are being used to access the large and complex genomes of crops and their wild relatives. These have helped to identify a wide spectrum of genetic variation and permitted the association of genetic diversity with diverse agronomic phenotypes. In combination with improved and automated phenotyping assays and functional genomic studies, genomics is providing new foundations for crop-breeding systems.

Genome sequence and genetic diversity of European ash trees.

Ash trees (genus Fraxinus, family Oleaceae) are widespread throughout the Northern Hemisphere, but are being devastated in Europe by the fungus Hymenoscyphus fraxineus, causing ash dieback, and in North America by the herbivorous beetle Agrilus planipennis. Here we sequence the genome of a low-heterozygosity Fraxinus excelsior tree from Gloucestershire, UK, annotating 38,852 protein-coding genes of which 25% appear ash specific when compared with the genomes of ten other plant species. Analyses of paralogous genes suggest a whole-genome duplication shared with olive (Olea europaea, Oleaceae). We also re-sequence 37 F. excelsior trees from Europe, finding evidence for apparent long-term decline in effective population size. Using our reference sequence, we re-analyse association transcriptomic data, yielding improved markers for reduced susceptibility to ash dieback. Surveys of these markers in British populations suggest that reduced susceptibility to ash dieback may be more widespread in Great Britain than in Denmark. We also present evidence that susceptibility of trees to H. fraxineus is associated with their iridoid glycoside levels. This rapid, integrated, multidisciplinary research response to an emerging health threat in a non-model organism opens the way for mitigation of the epidemic.

Copy number variation of TdDof controls solid-stemmed architecture in wheat.

Stem solidness is an important agronomic trait of durum (Triticum turgidum L. var. durum) and bread (Triticum aestivum L.) wheat that provides resistance to the wheat stem sawfly. This dominant trait is conferred by the SSt1 locus on chromosome 3B. However, the molecular identity and mechanisms underpinning stem solidness have not been identified. Here, we demonstrate that copy number variation of TdDof, a gene encoding a putative DNA binding with one finger protein, controls the stem solidness trait in wheat. Using map-based cloning, we localized TdDof to within a physical interval of 2.1 Mb inside the SSt1 locus. Molecular analysis revealed that hollow-stemmed wheat cultivars such as Kronos carry a single copy of TdDof, whereas solid-stemmed cultivars such as CDC Fortitude carry multiple identical copies of the gene. Deletion of all TdDof copies from CDC Fortitude resulted in the loss of stem solidness, whereas the transgenic overexpression of TdDof restored stem solidness in the TdDof deletion mutant pithless1 and conferred stem solidness in Kronos. In solid-stemmed cultivars, increased TdDof expression was correlated with the down-regulation of genes whose orthologs have been implicated in programmed cell death (PCD) in other species. Anatomical and histochemical analyses revealed that hollow-stemmed lines had stronger PCD-associated signals in the pith cells compared to solid-stemmed lines, which suggests copy number-dependent expression of TdDof could be directly or indirectly involved in the negative regulation of PCD. These findings provide opportunities to manipulate stem development in wheat and other monocots for agricultural or industrial purposes.

The Triticum ispahanicum elongated glume locus P2 maps to chromosome 6A and is associated with the ectopic expression of SVP-A1.

KEY MESSAGE: We propose the MADS-box transcription factor SVP-A1 as a promising candidate gene for the elongated glume locus P2, which maps to chromosome 6A instead of the previously proposed chromosome 7B. In rice and wheat, glume and floral organ length are positively correlated with grain size, making them an important target to increase grain size and potentially yield. The wheat subspecies Triticum ispahanicum is known to develop elongated glumes and floral organs as well as long grains. These multiple phenotypic effects are controlled by the P2 locus, which was previously mapped to wheat chromosome 7B. Using three mapping populations, we show that the long glume locus P2 does not map to chromosome 7B, but instead maps to a 1.68 Mbp interval on chromosome 6A. Within this interval, we identified SVP-A1, a MADS box transcription factor which is the direct ortholog of the maize gene underlying the 'pod corn' Tunicate locus and is a paralog to the T. polonicum elongated glume P1 gene. In T. ispahanicum, we identified a unique allele which has a 482-bp deletion in the SVP-A1 promoter and is associated with ectopic and higher expression of SVP-A1 in the elongated glumes and floral organs. We used near-isogenic lines (NILs) to show that P2 has a consistent positive effect on the length of glume, lemma, palea, spike and grain. Based on the mapping data, natural variation, biological function of SVP genes in cereals and expression analyses, we propose the MADS-box transcription factor SVP-A1 as a promising candidate for P2.

Distribution of Puccinia striiformis f. sp. tritici Races and Virulence in Wheat Growing Regions of Kenya from 1970 to 2014.

Stripe rust, caused by the fungal pathogen Puccinia striiformis f. sp. tritici, is a major threat to wheat (Triticum spp.) production worldwide. The objective of this study was to determine the virulence of P. striiformis f. sp. tritici races prevalent in the main wheat growing regions of Kenya, which includes Mt. Kenya, Eastern Kenya, and the Rift Valley (Central, Southern, and Northern Rift). Fifty P. striiformis f. sp. tritici isolates collected from 1970 to 1992 and from 2009 to 2014 were virulence phenotyped with stripe rust differential sets, and 45 isolates were genotyped with sequence characterized amplified region (SCAR) markers to differentiate the isolates and identify aggressive strains PstS1 and PstS2. Virulence corresponding to stripe rust resistance genes Yr1, Yr2, Yr3, Yr6, Yr7, Yr8, Yr9, Yr17, Yr25, and Yr27 and the seedling resistance in genotype Avocet S were detected. Ten races were detected in the P. striiformis f. sp. tritici samples obtained from 1970 to 1992, and three additional races were detected from 2009 to 2014, with a single race being detected in both periods. The SCAR markers detected both Pst1 and Pst2 strains in the collection. Increasing P. striiformis f. sp. tritici virulence was found in the Kenyan P. striiformis f. sp. tritici population, and different P. striiformis f. sp. tritici race groups were found to dominate different wheat growing regions. Moreover, recent P. striiformis f. sp. tritici races in East Africa indicated possible migration of some race groups into Kenya from other regions. This study is important in elucidating P. striiformis f. sp. tritici evolution and virulence diversity and useful in breeding wheat cultivars with effective resistance to stripe rust.

The NLR-Annotator Tool Enables Annotation of the Intracellular Immune Receptor Repertoire.

Disease resistance genes encoding nucleotide-binding and leucine-rich repeat (NLR) intracellular immune receptor proteins detect pathogens by the presence of pathogen effectors. Plant genomes typically contain hundreds of NLR-encoding genes. The availability of the hexaploid wheat (Triticum aestivum) cultivar Chinese Spring reference genome allows a detailed study of its NLR complement. However, low NLR expression and high intrafamily sequence homology hinder their accurate annotation. Here, we developed NLR-Annotator, a software tool for in silico NLR identification independent of transcript support. Although developed for wheat, we demonstrate the universal applicability of NLR-Annotator across diverse plant taxa. We applied our tool to wheat and combined it with a transcript-validated subset of genes from the reference gene annotation to characterize the structure, phylogeny, and expression profile of the NLR gene family. We detected 3,400 full-length NLR loci, of which 1,560 were confirmed as expressed genes with intact open reading frames. NLRs with integrated domains mostly group in specific subclades. Members of another subclade predominantly locate in close physical proximity to NLRs carrying integrated domains, suggesting a paired helper function. Most NLRs (88%) display low basal expression (in the lower 10 percentile of transcripts). In young leaves subjected to biotic stress, we found up-regulation of 266 of the NLRs To illustrate the utility of our tool for the positional cloning of resistance genes, we estimated the number of NLR genes within the intervals of mapped rust resistance genes. Our study will support the identification of functional resistance genes in wheat to accelerate the breeding and engineering of disease-resistant varieties.