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.

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.

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.

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.

A heat-shock inducible system for flexible gene expression in cereals

BACKGROUND

Functional characterisation of genes using transgenic methods is increasingly common in cereal crops. Yet standard methods of gene over-expression can lead to undesirable developmental phenotypes, or even embryo lethality, due to ectopic gene expression. Inducible expression systems allow the study of such genes by preventing their expression until treatment with the specific inducer. When combined with the Cre-Lox recombination system, inducible promoters can be used to initiate constitutive expression of a gene of interest. Yet while these systems are well established in dicot model plants, like Arabidopsis thaliana, they have not yet been implemented in grasses.

RESULTS

Here we present an irreversible heat-shock inducible system developed using Golden Gate-compatible components which utilises Cre recombinase to drive constitutive gene expression in barley and wheat. We show that a heat shock treatment of 38 °C is sufficient to activate the construct and drive expression of the gene of interest. Modulating the duration of heat shock controls the density of induced cells. Short durations of heat shock cause activation of the construct in isolated single cells, while longer durations lead to global construct activation. The system can be successfully activated in multiple tissues and at multiple developmental stages and shows no activation at standard growth temperatures (~ 20 °C).

CONCLUSIONS

This system provides an adaptable framework for use in gene functional characterisation in cereal crops. The developed vectors can be easily adapted for specific genes of interest within the Golden Gate cloning system. By using an environmental signal to induce activation of the construct, the system avoids pitfalls associated with consistent and complete application of chemical inducers. As with any inducible system, care must be taken to ensure that the expected construct activation has indeed taken place.

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.

LYS3 encodes a prolamin-box-binding transcription factor that controls embryo growth in barley and wheat

Mutations at the LYS3 locus in barley have multiple effects on grain development, including an increase in embryo size and a decrease in endosperm starch content. The gene underlying LYS3 was identified by genetic mapping and mutations in this gene were identified in all four barley lys3 alleles. LYS3 encodes a transcription factor called Prolamin Binding Factor (PBF). Its role in controlling embryo size was confirmed using wheat TILLING mutants. To understand how PBF controls embryo development, we studied its spatial and temporal patterns of expression in developing grains. The PBF gene is expressed in both the endosperm and the embryos, but the timing of expression in these organs differs. PBF expression in wild-type embryos precedes the onset of embryo enlargement in lys3 mutants, suggesting that PBF suppresses embryo growth. We predicted the down-stream target genes of PBF in wheat and found them to be involved in a wide range of biological processes, including organ development and starch metabolism. Our work suggests that PBF may influence embryo size and endosperm starch synthesis via separate gene control networks.

•

LYS3 encodes a transcription factor called Prolamin Binding Factor (PBF) that is expressed in grains only.

•

Wheat and barley LYS3/PBF mutants have enlarged embryos suggesting that this gene suppresses embryo growth.

•

The predicted targets of wheat PBF are involved in a wide range of biological processes, including organ development.

A roadmap for gene functional characterisation in crops with large genomes: Lessons from polyploid wheat

Understanding the function of genes within staple crops will accelerate crop improvement by allowing targeted breeding approaches. Despite their importance, a lack of genomic information and resources has hindered the functional characterisation of genes in major crops. The recent release of high-quality reference sequences for these crops underpins a suite of genetic and genomic resources that support basic research and breeding. For wheat, these include gene model annotations, expression atlases and gene networks that provide information about putative function. Sequenced mutant populations, improved transformation protocols and structured natural populations provide rapid methods to study gene function directly. We highlight a case study exemplifying how to integrate these resources. This review provides a helpful guide for plant scientists, especially those expanding into crop research, to capitalise on the discoveries made in Arabidopsis and other plants. This will accelerate the improvement of crops of vital importance for food and nutrition security.

A Co-Expression Network in Hexaploid Wheat Reveals Mostly Balanced Expression and Lack of Significant Gene Loss of Homeologous Meiotic Genes Upon Polyploidization

Polyploidization has played an important role in plant evolution. However, upon polyploidization, the process of meiosis must adapt to ensure the proper segregation of increased numbers of chromosomes to produce balanced gametes. It has been suggested that meiotic gene (MG) duplicates return to a single copy following whole genome duplication to stabilize the polyploid genome. Therefore, upon the polyploidization of wheat, a hexaploid species with three related (homeologous) genomes, the stabilization process may have involved rapid changes in content and expression of MGs on homeologous chromosomes (homeologs). To examine this hypothesis, sets of candidate MGs were identified in wheat using co-expression network analysis and orthology informed approaches. In total, 130 RNA-Seq samples from a range of tissues including wheat meiotic anthers were used to define co-expressed modules of genes. Three modules were significantly correlated with meiotic tissue samples but not with other tissue types. These modules were enriched for GO terms related to cell cycle, DNA replication, and chromatin modification and contained orthologs of known MGs. Overall, 74.4% of genes within these meiosis-related modules had three homeologous copies which was similar to other tissue-related modules. Amongst wheat MGs identified by orthology, rather than co-expression, the majority (93.7%) were either retained in hexaploid wheat at the same number of copies (78.4%) or increased in copy number (15.3%) compared to ancestral wheat species. Furthermore, genes within meiosis-related modules showed more balanced expression levels between homeologs than genes in non-meiosis-related modules. Taken together, our results do not support extensive gene loss nor changes in homeolog expression of MGs upon wheat polyploidization. The construction of the MG co-expression network allowed identification of hub genes and provided key targets for future studies.

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.

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.

Identification of a Dominant Chlorosis Phenotype Through a Forward Screen of the Triticum turgidum cv. Kronos TILLING Population

Durum wheat (Triticum turgidum) derives from a hybridization event approximately 400,000 years ago which led to the creation of an allotetraploid genome. The evolutionary recent origin of durum wheat means that its genome has not yet been fully diploidised. As a result, many of the genes present in the durum genome act in a redundant fashion, where loss-of-function mutations must be present in both gene copies to observe a phenotypic effect. Here, we use a novel set of induced variation within the cv. Kronos TILLING population to identify a locus controlling a dominant, environmentally dependent chlorosis phenotype. We carried out a forward screen of the sequenced cv. Kronos TILLING lines for senescence phenotypes and identified a line with a dominant early senescence and chlorosis phenotype. Mutant plants contained less chlorophyll throughout their development and displayed premature flag leaf senescence. A segregating population was classified into discrete phenotypic groups and subjected to bulked-segregant analysis using exome capture followed by next-generation sequencing. This allowed the identification of a single region on chromosome 3A, Yellow Early Senescence 1 (YES-1), which was associated with the mutant phenotype. While this phenotype was consistent across 4 years of field trials in the United Kingdom, the mutant phenotype was not observed when grown in Davis, CA (United States). To obtain further SNPs for fine-mapping, we isolated chromosome 3A using flow sorting and sequenced the entire chromosome. By mapping these reads against both the cv. Chinese Spring reference sequence and the cv. Kronos assembly, we could identify high-quality, novel EMS-induced SNPs in non-coding regions within YES-1 that were previously missed in the exome capture data. This allowed us to fine-map YES-1 to 4.3 Mb, containing 59 genes. Our study shows that populations containing induced variation can be sources of novel dominant variation in polyploid crop species, highlighting their importance in future genetic screens. We also demonstrate the value of using cultivar-specific genome assemblies alongside the gold-standard reference genomes particularly when working with non-coding regions of the genome. Further fine-mapping of the YES-1 locus will be pursued to identify the causal SNP underpinning this dominant, environmentally dependent phenotype.

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.

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.

Genome-Wide Transcription During Early Wheat Meiosis Is Independent of Synapsis, Ploidy Level, and the Ph1 Locus

Polyploidization is a fundamental process in plant evolution. One of the biggest challenges faced by a new polyploid is meiosis, particularly discriminating between multiple related chromosomes so that only homologous chromosomes synapse and recombine to ensure regular chromosome segregation and balanced gametes. Despite its large genome size, high DNA repetitive content and similarity between homoeologous chromosomes, hexaploid wheat completes meiosis in a shorter period than diploid species with a much smaller genome. Therefore, during wheat meiosis, mechanisms additional to the classical model based on DNA sequence homology, must facilitate more efficient homologous recognition. One such mechanism could involve exploitation of differences in chromosome structure between homologs and homoeologs at the onset of meiosis. In turn, these chromatin changes, can be expected to be linked to transcriptional gene activity. In this study, we present an extensive analysis of a large RNA-seq data derived from six different genotypes: wheat, wheat-rye hybrids and newly synthesized octoploid triticale, both in the presence and absence of the Ph1 locus. Plant material was collected at early prophase, at the transition leptotene-zygotene, when the telomere bouquet is forming and synapsis between homologs is beginning. The six genotypes exhibit different levels of synapsis and chromatin structure at this stage; therefore, recombination and consequently segregation, are also different. Unexpectedly, our study reveals that neither synapsis, whole genome duplication nor the absence of the Ph1 locus are associated with major changes in gene expression levels during early meiotic prophase. Overall wheat transcription at this meiotic stage is therefore highly resilient to such alterations, even in the presence of major chromatin structural changes. Further studies in wheat and other polyploid species will be required to reveal whether these observations are specific to wheat meiosis.

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.

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.

The transcriptional landscape of polyploid wheat

The coordinated expression of highly related homoeologous genes in polyploid species underlies the phenotypes of many of the world's major crops. Here we combine extensive gene expression datasets to produce a comprehensive, genome-wide analysis of homoeolog expression patterns in hexaploid bread wheat. Bias in homoeolog expression varies between tissues, with ~30% of wheat homoeologs showing nonbalanced expression. We found expression asymmetries along wheat chromosomes, with homoeologs showing the largest inter-tissue, inter-cultivar, and coding sequence variation, most often located in high-recombination distal ends of chromosomes. These transcriptionally dynamic genes potentially represent the first steps toward neo- or subfunctionalization of wheat homoeologs. Coexpression networks reveal extensive coordination of homoeologs throughout development and, alongside a detailed expression atlas, provide a framework to target candidate genes underpinning agronomic traits in wheat.

Impact of transposable elements on genome structure and evolution in bread wheat

BACKGROUND

Transposable elements (TEs) are major components of large plant genomes and main drivers of genome evolution. The most recent assembly of hexaploid bread wheat recovered the highly repetitive TE space in an almost complete chromosomal context and enabled a detailed view into the dynamics of TEs in the A, B, and D subgenomes.

RESULTS

The overall TE content is very similar between the A, B, and D subgenomes, although we find no evidence for bursts of TE amplification after the polyploidization events. Despite the near-complete turnover of TEs since the subgenome lineages diverged from a common ancestor, 76% of TE families are still present in similar proportions in each subgenome. Moreover, spacing between syntenic genes is also conserved, even though syntenic TEs have been replaced by new insertions over time, suggesting that distances between genes, but not sequences, are under evolutionary constraints. The TE composition of the immediate gene vicinity differs from the core intergenic regions. We find the same TE families to be enriched or depleted near genes in all three subgenomes. Evaluations at the subfamily level of timed long terminal repeat-retrotransposon insertions highlight the independent evolution of the diploid A, B, and D lineages before polyploidization and cases of concerted proliferation in the AB tetraploid.

CONCLUSIONS

Even though the intergenic space is changed by the TE turnover, an unexpected preservation is observed between the A, B, and D subgenomes for features like TE family proportions, gene spacing, and TE enrichment near genes.

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

An annotated reference sequence representing the hexaploid bread wheat genome in 21 pseudomolecules has been analyzed to identify the distribution and genomic context of coding and noncoding elements across the A, B, and D subgenomes. With an estimated coverage of 94% of the genome and containing 107,891 high-confidence gene models, this assembly enabled the discovery of tissue- and developmental stage-related coexpression networks by providing a transcriptome atlas representing major stages of wheat development. Dynamics of complex gene families involved in environmental adaptation and end-use quality were revealed at subgenome resolution and contextualized to known agronomic single-gene or quantitative trait loci. This community resource establishes the foundation for accelerating wheat research and application through improved understanding of wheat biology and genomics-assisted breeding.

Genome-Wide Sequence and Expression Analysis of the NAC Transcription Factor Family in Polyploid Wheat

Many important genes in agriculture correspond to transcription factors (TFs) that regulate a wide range of pathways from flowering to responses to disease and abiotic stresses. In this study, we identified 5776 TFs in hexaploid wheat (Triticum aestivum) and classified them into gene families. We further investigated the NAC family exploring the phylogeny, C-terminal domain (CTD) conservation, and expression profiles across 308 RNA-seq samples. Phylogenetic trees of NAC domains indicated that wheat NACs divided into eight groups similar to rice (Oryza sativa) and barley (Hordeum vulgare). CTD motifs were frequently conserved between wheat, rice, and barley within phylogenetic groups; however, this conservation was not maintained across phylogenetic groups. Three homeologous copies were present for 58% of NACs, whereas evidence of single homeolog gene loss was found for 33% of NACs. We explored gene expression patterns across a wide range of developmental stages, tissues, and abiotic stresses. We found that more phylogenetically related NACs shared more similar expression patterns compared to more distant NACs. However, within each phylogenetic group there were clades with diverse expression profiles. We carried out a coexpression analysis on all wheat genes and identified 37 modules of coexpressed genes of which 23 contained NACs. Using gene ontology (GO) term enrichment, we obtained putative functions for NACs within coexpressed modules including responses to heat and abiotic stress and responses to water: these NACs may represent targets for breeding or biotechnological applications. This study provides a framework and data for hypothesis generation for future studies on NAC TFs in wheat.

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.

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.

Systematic Investigation of FLOWERING LOCUS T-Like Poaceae Gene Families Identifies the Short-Day Expressed Flowering Pathway Gene, TaFT3 in Wheat (Triticum aestivum L.)

To date, a small number of major flowering time loci have been identified in the related Triticeae crops, bread wheat (Triticum aestivum), durum wheat (T. durum), and barley (Hordeum vulgare). Natural genetic variants at these loci result in major phenotypic changes which have adapted crops to the novel environments encountered during the spread of agriculture. The polyploid nature of bread and durum wheat means that major flowering time loci in which recessive alleles confer adaptive advantage in related diploid species have not been readily identified. One such example is the PPD-H2 flowering time locus encoded by FLOWERING LOCUS T 3 (HvFT3) in the diploid crop barley, for which recessive mutant alleles confer delayed flowering under short day (SD) photoperiods. In autumn-sown barley, such alleles aid the repression of flowering over the winter, which help prevent the development of cold-sensitive floral organs until the onset of inductive long day (LD) photoperiods the following spring. While the identification of orthologous loci in wheat could provide breeders with alternative mechanisms to fine tune flowering time, systematic identification of wheat orthologs of HvFT3 has not been reported. Here, we characterize the FT gene families in six Poaceae species, identifying novel members in all taxa investigated, as well as FT3 homoeologs from the A, B and D genomes of hexaploid (TaFT3) and tetraploid wheat. Sequence analysis shows TaFT3 homoeologs display high similarity to the HvFT3 coding region (95-96%) and predicted protein (96-97%), with conservation of intron/exon structure across the five cereal species investigated. Genetic mapping and comparative analyses in hexaploid and tetraploid wheat find TaFT3 homoeologs map to the long arms of the group 1 chromosomes, collinear to HvFT3 in barley and FT3 orthologs in rice, foxtail millet and brachypodium. Genome-specific expression analyses show FT3 homoeologs in tetraploid and hexaploid wheat are upregulated under SD photoperiods, but not under LDs, analogous to the expression of HvFT3. Collectively, these results indicate that functional wheat orthologs of HvFT3 have been identified. The molecular resources generated here provide the foundation for engineering a novel major flowering time locus in wheat using forward or reverse genetics approaches.

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 Grain Filling Is Limited by Grain Filling Capacity rather than the Duration of Flag Leaf Photosynthesis: A Case Study Using NAM RNAi Plants

It has been proposed that delayed leaf senescence can extend grain filling duration and thus increase yields in cereal crops. We found that wheat (Triticum aestivum) NAM RNAi plants with delayed senescence carried out 40% more flag leaf photosynthesis after anthesis than control plants, but had the same rate and duration of starch accumulation during grain filling and the same final grain weight. The additional photosynthate available in NAM RNAi plants was in part stored as fructans in the stems, whereas stem fructans were remobilised during grain filling in control plants. In both genotypes, activity of starch synthase was limiting for starch synthesis in the later stages of grain filling. We suggest that in order to realise the potential yield gains offered by delayed leaf senescence, this trait should be combined with increased grain filling capacity.

Biofortification of wheat grain with iron and zinc: integrating novel genomic resources and knowledge from model crops

Wheat, like many other staple cereals, contains low levels of the essential micronutrients iron and zinc. Up to two billion people worldwide suffer from iron and zinc deficiencies, particularly in regions with predominantly cereal-based diets. Although wheat flour is commonly fortified during processing, an attractive and more sustainable solution is biofortification, which requires developing new varieties of wheat with inherently higher iron and zinc content in their grains. Until now most studies aimed at increasing iron and zinc content in wheat grains have focused on discovering natural variation in progenitor or related species. However, recent developments in genomics and transformation have led to a step change in targeted research on wheat at a molecular level. We discuss promising approaches to improve iron and zinc content in wheat using knowledge gained in model grasses. We explore how the latest resources developed in wheat, including sequenced genomes and mutant populations, can be exploited for biofortification. We also highlight the key research and practical challenges that remain in improving iron and zinc content in wheat.

Arabidopsis plants perform arithmetic division to prevent starvation at night

Photosynthetic starch reserves that accumulate in Arabidopsis leaves during the day decrease approximately linearly with time at night to support metabolism and growth. We find that the rate of decrease is adjusted to accommodate variation in the time of onset of darkness and starch content, such that reserves last almost precisely until dawn. Generation of these dynamics therefore requires an arithmetic division computation between the starch content and expected time to dawn. We introduce two novel chemical kinetic models capable of implementing analog arithmetic division. Predictions from the models are successfully tested in plants perturbed by a night-time light period or by mutations in starch degradation pathways. Our experiments indicate which components of the starch degradation apparatus may be important for appropriate arithmetic division. Our results are potentially relevant for any biological system dependent on a food reserve for survival over a predictable time period. DOI:http://dx.doi.org/10.7554/eLife.00669.001.

Arabidopsis plants perform arithmetic division to prevent starvation at night

Photosynthetic starch reserves that accumulate in Arabidopsis leaves during the day decrease approximately linearly with time at night to support metabolism and growth. We find that the rate of decrease is adjusted to accommodate variation in the time of onset of darkness and starch content, such that reserves last almost precisely until dawn. Generation of these dynamics therefore requires an arithmetic division computation between the starch content and expected time to dawn. We introduce two novel chemical kinetic models capable of implementing analog arithmetic division. Predictions from the models are successfully tested in plants perturbed by a night-time light period or by mutations in starch degradation pathways. Our experiments indicate which components of the starch degradation apparatus may be important for appropriate arithmetic division. Our results are potentially relevant for any biological system dependent on a food reserve for survival over a predictable time period. DOI:http://dx.doi.org/10.7554/eLife.00669.001.

Multiple Arabidopsis genes primed for recruitment into C₄ 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.