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.

An autoactive NB-LRR gene causes Rht13 dwarfism in wheat.

Semidwarfing genes have greatly increased wheat yields globally, yet the widely used gibberellin (GA)-insensitive genes Rht-B1b and Rht-D1b have disadvantages for seedling emergence. Use of the GA-sensitive semidwarfing gene Rht13 avoids this pleiotropic effect. Here, we show that Rht13 encodes a nucleotide-binding site/leucine-rich repeat (NB-LRR) gene. A point mutation in the semidwarf Rht-B13b allele autoactivates the NB-LRR gene and causes a height reduction comparable with Rht-B1b and Rht-D1b in diverse genetic backgrounds. The autoactive Rht-B13b allele leads to transcriptional up-regulation of pathogenesis-related genes including class III peroxidases associated with cell wall remodeling. Rht13 represents a new class of reduced height (Rht) gene, unlike other Rht genes, which encode components of the GA signaling or metabolic pathways. This discovery opens avenues to use autoactive NB-LRR genes as semidwarfing genes in a range of crop species, and to apply Rht13 in wheat breeding programs using a perfect genetic marker.

Applying the latest advances in genomics and phenomics for trait discovery in polyploid wheat.

Improving traits in wheat has historically been challenging due to its large and polyploid genome, limited genetic diversity and in-field phenotyping constraints. However, within recent years many of these barriers have been lowered. The availability of a chromosome-level assembly of the wheat genome now facilitates a step-change in wheat genetics and provides a common platform for resources, including variation data, gene expression data and genetic markers. The development of sequenced mutant populations and gene-editing techniques now enables the rapid assessment of gene function in wheat directly. The ability to alter gene function in a targeted manner will unmask the effects of homoeolog redundancy and allow the hidden potential of this polyploid genome to be discovered. New techniques to identify and exploit the genetic diversity within wheat wild relatives now enable wheat breeders to take advantage of these additional sources of variation to address challenges facing food production. Finally, advances in phenomics have unlocked rapid screening of populations for many traits of interest both in greenhouses and in the field. Looking forwards, integrating diverse data types, including genomic, epigenetic and phenomics data, will take advantage of big data approaches including machine learning to understand trait biology in wheat in unprecedented detail.

An improved assembly and annotation of the allohexaploid wheat genome identifies complete families of agronomic genes and provides genomic evidence for chromosomal translocations.

Advances in genome sequencing and assembly technologies are generating many high-quality genome sequences, but assemblies of large, repeat-rich polyploid genomes, such as that of bread wheat, remain fragmented and incomplete. We have generated a new wheat whole-genome shotgun sequence assembly using a combination of optimized data types and an assembly algorithm designed to deal with large and complex genomes. The new assembly represents >78% of the genome with a scaffold N50 of 88.8 kb that has a high fidelity to the input data. Our new annotation combines strand-specific Illumina RNA-seq and Pacific Biosciences (PacBio) full-length cDNAs to identify 104,091 high-confidence protein-coding genes and 10,156 noncoding RNA genes. We confirmed three known and identified one novel genome rearrangements. Our approach enables the rapid and scalable assembly of wheat genomes, the identification of structural variants, and the definition of complete gene models, all powerful resources for trait analysis and breeding of this key global crop.

Blurring the boundaries between cereal crops and model plants.

The cereal crops rice (Oryza sativa), maize (Zea mays ssp. mays) and wheat (Triticum aestivum) provide half of the food eaten by humankind. However, understanding their biology has proved challenging due to their large size, long lifecycle and large genomes. The model plant Arabidopsis thaliana avoids these practical problems and has provided fundamental understanding of plant biology, however not all of this knowledge is directly transferrable to cereals. Recent developments in gene editing, speed breeding and genome assembly techniques mean that the challenges associated with working with the major cereal crops can be overcome. Resources such as mutant collections and genome sequences are now available for these crops, making them attractive experimental systems with which to make discoveries that are directly applicable to increasing crop production.

Identification of Transcription Factors Regulating Senescence in Wheat through Gene Regulatory Network Modelling.

Senescence is a tightly regulated developmental program coordinated by transcription factors. Identifying these transcription factors in crops will provide opportunities to tailor the senescence process to different environmental conditions and regulate the balance between yield and grain nutrient content. Here, we use ten time points of gene expression data along with gene network modeling to identify transcription factors regulating senescence in polyploid wheat (Triticum aestivum). We observe two main phases of transcriptional changes during senescence: early down-regulation of housekeeping functions and metabolic processes followed by up-regulation of transport and hormone-related genes. These two phases are largely conserved with Arabidopsis (Arabidopsis thaliana), although the individual genes underlying these changes are often not orthologous. We have identified transcription factor families associated with these early and later waves of differential expression. Using gene regulatory network modeling, we identified candidate transcription factors that may control senescence. Using independent, publicly available datasets, we found that the most highly ranked candidate genes in the network were enriched for senescence-related functions compared with all genes in the network. We validated the function of one of these candidate transcription factors in senescence using wheat chemically induced mutants. This study lays the groundwork to understand the transcription factors that regulate senescence in polyploid wheat and exemplifies the integration of time-series data with publicly available expression atlases and networks to identify candidate regulatory genes.

Overgrowth mutants determine the causal role of gibberellin GA2oxidaseA13 in Rht12 dwarfism of wheat.

The induced dwarf mutant Rht12 was previously shown to have agronomic potential to replace the conventional DELLA mutants Rht-B1b/Rht-D1b in wheat. The Rht12 dwarfing gene is not associated with reduced coleoptile length (unlike the DELLA mutants) and it is dominant, characteristics which are shared with the previously characterized dwarfing genes Rht18 and Rht14. Using the Rht18/Rht14 model, a gibberellin (GA) 2-oxidase gene was identified in the Rht12 region on chromosome 5A. A screen for suppressor mutants in the Rht12 background identified tall overgrowth individuals that were shown to contain loss-of-function mutations in GA2oxidaseA13, demonstrating the role of this gene in the Rht12 dwarf phenotype. It was concluded that Rht12, Rht18, and Rht14 share the same height-reducing mechanism through the increased expression of GA 2-oxidase genes. Some of the overgrowth mutants generated in this study were semi-dwarf and taller than the original Rht12 dwarf, providing breeders with new sources of agronomically useful dwarfism.

Applying genomic resources to accelerate wheat biofortification.

Wheat has low levels of the micronutrients iron and zinc in the grain, which contributes to 2 billion people suffering from micronutrient deficiency globally. While wheat flour is commonly fortified during processing, an attractive and more sustainable solution is biofortification, which could improve micronutrient content in the human diet, without the sustainability issues and costs associated with conventional fortification. Although many studies have used quantitative trait loci mapping and genome-wide association to identify genetic loci to improve micronutrient contents, recent developments in genomics offer an opportunity to accelerate marker discovery and use gene-focussed approaches to engineer improved micronutrient content in wheat. The recent publication of a high-quality wheat genome sequence, alongside gene expression atlases, variation datasets and sequenced mutant populations, provides a foundation to identify genetic loci and genes controlling micronutrient content in wheat. We discuss how novel genomic resources can identify candidate genes for biofortification, integrating knowledge from other cereal crops, and how these genes can be tested using gene editing, transgenic and TILLING approaches. Finally, we highlight key challenges remaining to develop wheat cultivars with high levels of iron and zinc.

Uncovering hidden variation in polyploid wheat.

Comprehensive reverse genetic resources, which have been key to understanding gene function in diploid model organisms, are missing in many polyploid crops. Young polyploid species such as wheat, which was domesticated less than 10,000 y ago, have high levels of sequence identity among subgenomes that mask the effects of recessive alleles. Such redundancy reduces the probability of selection of favorable mutations during natural or human selection, but also allows wheat to tolerate high densities of induced mutations. Here we exploited this property to sequence and catalog more than 10 million mutations in the protein-coding regions of 2,735 mutant lines of tetraploid and hexaploid wheat. We detected, on average, 2,705 and 5,351 mutations per tetraploid and hexaploid line, respectively, which resulted in 35-40 mutations per kb in each population. With these mutation densities, we identified an average of 23-24 missense and truncation alleles per gene, with at least one truncation or deleterious missense mutation in more than 90% of the captured wheat genes per population. This public collection of mutant seed stocks and sequence data enables rapid identification of mutations in the different copies of the wheat genes, which can be combined to uncover previously hidden variation. Polyploidy is a central phenomenon in plant evolution, and many crop species have undergone recent genome duplication events. Therefore, the general strategy and methods developed herein can benefit other polyploid crops.

Hotspots in the genomic architecture of field drought responses in wheat as breeding targets.

Wheat can adapt to most agricultural conditions across temperate regions. This success is the result of phenotypic plasticity conferred by a large and complex genome composed of three homoeologous genomes (A, B, and D). Although drought is a major cause of yield and quality loss in wheat, the adaptive mechanisms and gene networks underlying drought responses in the field remain largely unknown. Here, we addressed this by utilizing an interdisciplinary approach involving field water status phenotyping, sampling, and gene expression analyses. Overall, changes at the transcriptional level were reflected in plant spectral traits amenable to field-level physiological measurements, although changes in photosynthesis-related pathways were found likely to be under more complex post-transcriptional control. Examining homoeologous genes with a 1:1:1 relationship across the A, B, and D genomes (triads), we revealed a complex genomic architecture for drought responses under field conditions, involving gene homoeolog specialization, multiple gene clusters, gene families, miRNAs, and transcription factors coordinating these responses. Our results provide a new focus for genomics-assisted breeding of drought-tolerant wheat cultivars.

Conserved residues in the wheat (Triticum aestivum) NAM-A1 NAC domain are required for protein binding and when mutated lead to delayed peduncle and flag leaf senescence.

BACKGROUND: NAC transcription factors contain five highly conserved subdomains which are required for protein dimerisation and DNA binding. Few residues within these subdomains have been identified as essential for protein function, and fewer still have been shown to be of biological relevance in planta. Here we use a positive regulator of senescence in wheat, NAM-A1, to test the impact of missense mutations at specific, highly conserved residues of the NAC domain on protein function. RESULTS: We identified missense mutations in five highly conserved residues of the NAC domain of NAM-A1 in a tetraploid TILLING population. TILLING lines containing these mutations, alongside synonymous and non-conserved mutation controls, were grown under glasshouse conditions and scored for senescence. Four of the five mutations showed a significant and consistent delay in peduncle senescence but had no consistent effects on flag leaf senescence. All four mutant alleles with the delayed senescence phenotype also lost the ability to interact with the homoeolog NAM-B1 in a yeast two-hybrid assay. Two of these residues were previously shown to be involved in NAC domain function in Arabidopsis, suggesting conservation of residue function between species. Three of these four alleles led to an attenuated cell death response compared to wild-type NAM-A1 when transiently over-expressed in Nicotiana benthamiana. One of these mutations was further tested under field conditions, in which there was a significant and consistent delay in both peduncle and leaf senescence. CONCLUSIONS: We combined field and glasshouse studies of a series of mutant alleles with biochemical analyses to identify four residues of the NAC domain which are required for NAM-A1 function and protein interaction. We show that mutations in these residues lead to a gradient of phenotypes, raising the possibility of developing allelic series of mutations for traits of agronomic importance. We also show that mutations in NAM-A1 more severely impact peduncle senescence, compared to the more commonly studied flag leaf senescence, highlighting this as an area deserving of further study. The results from this integrated approach provide strong evidence that conserved residues within the functional domains of NAC transcription factors have biological significance in planta.

Final grain weight is not limited by the activity of key starch-synthesising enzymes during grain filling in wheat.

Since starch is by far the major component of the mature wheat grain, it has been assumed that variation in the capacity for starch synthesis during grain filling can influence final grain weight. We investigated this assumption by studying a total of 54 wheat genotypes including elite varieties and landraces that were grown in two successive years in fields in the east of England. The weight, water content, sugars, starch, and maximum catalytic activities of two enzymes of starch biosynthesis, ADP-glucose pyrophosphorylase and soluble starch synthase, were measured during grain filling. The relationships between these variables and the weights and starch contents of mature grains were analysed. Final grain weight showed few or no significant correlations with enzyme activities, sugar levels, or starch content during grain filling, or with starch content at maturity. We conclude that neither sugar availability nor enzymatic capacity for starch synthesis during grain filling significantly influenced final grain weight in our field conditions. We suggest that final grain weight may be largely determined by developmental processes prior to grain filling. Starch accumulation then fills the grain to a physical limit set by developmental processes. This conclusion is in accord with those from previous studies in which source or sink strength has been artificially manipulated.

expVIP: a Customizable RNA-seq Data Analysis and Visualization Platform.

The majority of transcriptome sequencing (RNA-seq) expression studies in plants remain underutilized and inaccessible due to the use of disparate transcriptome references and the lack of skills and resources to analyze and visualize these data. We have developed expVIP, an expression visualization and integration platform, which allows easy analysis of RNA-seq data combined with an intuitive and interactive interface. Users can analyze public and user-specified data sets with minimal bioinformatics knowledge using the expVIP virtual machine. This generates a custom Web browser to visualize, sort, and filter the RNA-seq data and provides outputs for differential gene expression analysis. We demonstrate expVIP's suitability for polyploid crops and evaluate its performance across a range of biologically relevant scenarios. To exemplify its use in crop research, we developed a flexible wheat (Triticum aestivum) expression browser (www.wheat-expression.com) that can be expanded with user-generated data in a local virtual machine environment. The open-access expVIP platform will facilitate the analysis of gene expression data from a wide variety of species by enabling the easy integration, visualization, and comparison of RNA-seq data across experiments.

Genomics as the key to unlocking the polyploid potential of wheat.

Polyploidy has played a central role in plant genome evolution and in the formation of new species such as tetraploid pasta wheat and hexaploid bread wheat. Until recently, the high sequence conservation between homoeologous genes, together with the large genome size of polyploid wheat, had hindered genomic analyses in this important crop species. In the past 5 yr, however, the advent of next-generation sequencing has radically changed the wheat genomics landscape. Here, we review a series of advances in genomic resources and tools for functional genomics that are shifting the paradigm of what is possible in wheat molecular genetics and breeding. We discuss how understanding the relationship between homoeologues can inform approaches to modulate the response of quantitative traits in polyploid wheat; we also argue that functional redundancy has 'locked up' a wide range of phenotypic variation in wheat. We explore how genomics provides key tools to inform targeted manipulation of multiple homoeologues, thereby allowing researchers and plant breeders to unlock the full polyploid potential of wheat.

Wheat NAC transcription factor NAC5-1 is a positive regulator of senescence.

Wheat (Triticum aestivum L.) is an important source of both calories and protein in global diets, but there is a trade-off between grain yield and protein content. The timing of leaf senescence could mediate this trade-off as it is associated with both declines in photosynthesis and nitrogen remobilization from leaves to grain. NAC transcription factors play key roles in regulating senescence timing. In rice, OsNAC5 expression is correlated with increased protein content and upregulated in senescing leaves, but the role of the wheat ortholog in senescence had not been characterized. We verified that NAC5-1 is the ortholog of OsNAC5 and that it is expressed in senescing flag leaves in wheat. To characterize NAC5-1, we combined missense mutations in NAC5-A1 and NAC5-B1 from a TILLING mutant population and overexpressed NAC5-A1 in wheat. Mutation in NAC5-1 was associated with delayed onset of flag leaf senescence, while overexpression of NAC5-A1 was associated with slightly earlier onset of leaf senescence. DAP-seq was performed to locate transcription factor binding sites of NAC5-1. Analysis of DAP-seq and comparison with other studies identified putative downstream target genes of NAC5-1 which could be associated with senescence. This work showed that NAC5-1 is a positive transcriptional regulator of leaf senescence in wheat. Further research is needed to test the effect of NAC5-1 on yield and protein content in field trials, to assess the potential to exploit this senescence regulator to develop high-yielding wheat while maintaining grain protein content.

Multiple Arabidopsis genes primed for recruitment into C4 photosynthesis.

C(4) photosynthesis occurs in the most productive crops and vegetation on the planet, and has become widespread because it allows increased rates of photosynthesis compared with the ancestral C(3) pathway. Leaves of C(4) plants typically possess complicated alterations to photosynthesis, such that its reactions are compartmented between mesophyll and bundle sheath cells. Despite its complexity, the C(4) pathway has arisen independently in 62 separate lineages of land plants, and so represents one of the most striking examples of convergent evolution known. We demonstrate that elements in untranslated regions (UTRs) of multiple genes important for C(4) photosynthesis contribute to the metabolic compartmentalization characteristic of a C(4) leaf. Either the 5' or the 3' UTR is sufficient for cell specificity, indicating that functional redundancy underlies this key aspect of C(4) gene expression. Furthermore, we show that orthologous PPDK and CA genes from the C(3) plant Arabidopsis thaliana are primed for recruitment into the C(4) pathway. Elements sufficient for M-cell specificity in C(4) leaves are also present in both the 5' and 3' UTRs of these C(3) A. thaliana genes. These data indicate functional latency within the UTRs of genes from C(3) species that have been recruited into the C(4) pathway. The repeated recruitment of pre-existing cis-elements in C(3) genes may have facilitated the evolution of C(4) photosynthesis. These data also highlight the importance of alterations in trans in producing a functional C(4) leaf, and so provide insight into both the evolution and molecular basis of this important type of photosynthesis.

Wheat NAM genes regulate the majority of early monocarpic senescence transcriptional changes including nitrogen remobilization genes.

Senescence enables the remobilization of nitrogen and micronutrients from vegetative tissues of wheat (Triticum aestivum L.) into the grain. Understanding the molecular players in this process will enable the breeding of wheat lines with tailored grain nutrient content. The NAC transcription factor NAM-B1 is associated with earlier senescence and higher levels of grain protein, iron, and zinc contents due to increased nutrient remobilization. To investigate how related NAM genes control nitrogen remobilization at the molecular level, we carried out a comparative transcriptomic study using flag leaves at 7 time points (3, 7, 10, 13, 15, 19, and 26 days after anthesis) in wild type and NAM RNA interference lines with reduced NAM gene expression. Approximately 2.5 times more genes were differentially expressed in wild type than NAM RNA interference plants during this early senescence time course (6,508 vs 2,605 genes). In both genotypes, differentially expressed genes were enriched for gene ontology terms related to photosynthesis, hormones, amino acid transport, and nitrogen metabolism. However, nitrogen metabolism genes including glutamine synthetase (GS1 and GS2), glutamate decarboxylase (GAD), glutamate dehydrogenase (GDH), and asparagine synthetase (ASN1) showed stronger or earlier differential expression in wild-type than in NAM RNA interference plants, consistent with higher nitrogen remobilization. The use of time course data identified the dynamics of NAM-regulated and NAM-independent gene expression changes during senescence and provides an entry point to functionally characterize the pathways regulating senescence and nutrient remobilization in wheat.

Transcription factor retention through multiple polyploidization steps in wheat.

Whole-genome duplication is widespread in plant evolutionary history and is followed by nonrandom gene loss to return to a diploid state. Across multiple angiosperm species, the retained genes tend to be dosage-sensitive regulatory genes such as transcription factors, yet data for younger polyploid species is sparse. Here, we analyzed the retention, expression, and genetic variation in transcription factors in the recent allohexaploid bread wheat (Triticum aestivum L.). By comparing diploid, tetraploid, and hexaploid wheat, we found that, following each of two hybridization and whole-genome duplication events, the proportion of transcription factors in the genome increased. Transcription factors were preferentially retained over other genes as homoeologous groups in tetraploid and hexaploid wheat. Across cultivars, transcription factor homoeologs contained fewer deleterious missense mutations than nontranscription factors, suggesting that transcription factors are maintained as three functional homoeologs in hexaploid wheat populations. Transcription factor homoeologs were more strongly coexpressed than nontranscription factors, indicating conservation of function between homoeologs. We found that the B3, MADS-M-type, and NAC transcription factor families were less likely to have three homoeologs present than other families, which was associated with low expression levels and high levels of tandem duplication. Together, our results show that transcription factors are preferentially retained in polyploid wheat genomes although there is variation between families. Knocking out one transcription factor homoeolog to alter gene dosage, using TILLING or CRISPR, could generate new phenotypes for wheat breeding.