The Transcriptional Landscape of Polyploid Wheats and Their Diploid Ancestors during Embryogenesis and Grain Development

Spatiotemporal transcriptomics reveals key gene regulation for grain yield and quality in wheat

BACKGROUND

Cereal grain size and quality are critical agronomic traits in crop production. Wheat grain development is governed by intricate regulatory networks that require precise spatiotemporal coordination of gene expression to establish functional compartments in different cell types.

RESULTS

Here, we perform a spatial transcriptomics study covering the early stages of wheat grain development, from 4 to 12 days after pollination. We classify the grain into 10 distinct cell types and identify 192 marker genes associated with them. WGCNA analysis reveals that highly expressed genes in different cell types exhibit distinct enrichment patterns, significantly influencing grain development and filling. Through co-expression and motif analyses, we identify a specific group of genes that may regulate wheat grain development, including TaABI3-B1, a transcription factor specifically expressed in the embryo and surrounding endosperm, which negatively affects embryo and grain size.

CONCLUSIONS

This study presents a comprehensive spatiotemporal transcriptional dataset for understanding wheat grain development. Additionally, it identifies key genetic resources with potential applications for improving wheat yield.

A telomere-to-telomere genome assembly coupled with multi-omic data provides insights into the evolution of hexaploid bread wheat

The complete assembly of vast and complex plant genomes, like the hexaploid wheat genome, remains challenging. Here we present CS-IAAS, a comprehensive telomere-to-telomere (T2T) gap-free Triticum aestivum L. genome, encompassing 14.51 billion base pairs and featuring all 21 centromeres and 42 telomeres. Annotation revealed 90.8 Mb additional centromeric satellite arrays and 5,611 rDNA units. Genome-wide rearrangements, centromeric elements, transposable element expansion and segmental duplications were deciphered during tetraploidization and hexaploidization, providing a comprehensive understanding of wheat subgenome evolution. Among them, transposable element insertions during hexaploidization greatly influenced gene expression balances, thus increasing the genome plasticity of transcriptional levels. Additionally, we generated 163,329 full-length cDNA sequences and proteomic data that helped annotate 141,035 high-confidence protein-coding genes. The complete T2T reference genome (CS-IAAS), along with its transcriptome and proteome, represents a significant step in our understanding of wheat genome complexity and provides insights for future wheat research and breeding.

Temporal Gene Expression Profiles From Pollination to Seed Maturity in Sorghum Provide Core Candidates for Engineering Seed Traits

Sorghum (Sorghum bicolor (L.) Moench) is a highly nutritional multipurpose millet crop. However, the genetic and molecular regulatory mechanisms governing sorghum grain development and the associated agronomic traits remain unexplored. In this study, we performed a comprehensive transcriptomic analysis of pistils collected 1-2 days before pollination, and developing seeds collected -2, 10, 20 and 30 days after pollination of S. bicolor variety M35-1. Out of 31 337 genes expressed in these stages, 12 804 were differentially expressed in the consecutive stages of seed development. These exhibited 10 dominant expression patterns correlated with the distinct pathways and gene functions. Functional analysis, based on the pathway mapping, transcription factor enrichment and orthology, delineated the key patterns associated with pollination, fertilization, early seed development, grain filling and seed maturation. Furthermore, colocalization with previously reported quantitative trait loci (QTLs) for grain weight/size revealed 48 differentially expressed genes mapping to these QTL regions. Comprehensive literature mining integrated with QTL mapping and expression data shortlisted 25, 17 and 8 core candidates for engineering grain size, starch and protein content, respectively.

20 years of the Bio-Analytic Resource for Plant Biology

The Bio-Analytic Resource for Plant Biology ('the BAR', at https://bar.utoronto.ca) is celebrating its 20th year in operation in 2025. The BAR encompasses and provides visualization tools for large 'omics data sets from plants. The BAR covers data from Arabidopsis, tomato, wheat, barley and 29 other plant species (with data for 2 others to be released soon). These data include nucleotide and protein sequence data, gene expression data, protein-protein and protein-DNA interactions, protein structures, subcellular localizations, and polymorphisms. The data are stored in more than 200 relational databases holding 186 GB of data and are presented to the researchers via web apps. These web apps provide data analysis and visualization tools. Some of the most popular tools are eFP ('electronic fluorescent pictograph') Browsers, ePlants and ThaleMine (an Arabidopsis-specific instance of InterMine). The BAR was designated a Global Core Biodata Resource in 2023. Like other GCBRs, the BAR has excellent operational stability, provides access without login requirement, and provides an API for researchers to be able to access BAR data programmatically. We present in this update a new overarching search tool called Gaia that permits easy access to all BAR data, powered by machine learning and artificial intelligence.

Deciphering the Transcriptional Regulatory Network Governing Starch and Storage Protein Biosynthesis in Wheat for Breeding Improvement

Starch and seed storage protein (SSP) composition profoundly impact wheat grain yield and quality. To unveil regulatory mechanisms governing their biosynthesis, transcriptome, and epigenome profiling is conducted across key endosperm developmental stages, revealing that chromatin accessibility, H3K27ac, and H3K27me3 collectively regulate SSP and starch genes with varying impact. Population transcriptome and phenotype analyses highlight accessible promoter regions' crucial role as a genetic variation resource, influencing grain yield and quality in a core collection of wheat accessions. Integration of time-serial RNA-seq and ATAC-seq enables the construction of a hierarchical transcriptional regulatory network governing starch and SSP biosynthesis, identifying 42 high-confidence novel candidates. These candidates exhibit overlap with genetic regions associated with grain size and quality traits, and their functional significance is validated through expression-phenotype association analysis among wheat accessions and loss-of-function mutants. Functional analysis of wheat abscisic acid insensitive 3-A1 (TaABI3-A1) with genome editing knock-out lines demonstrates its role in promoting SSP accumulation while repressing starch biosynthesis through transcriptional regulation. Excellent TaABI3-A1Hap1 with enhanced grain weight is selected during the breeding process in China, linked to altered expression levels. This study unveils key regulators, advancing understanding of SSP and starch biosynthesis regulation and contributing to breeding enhancement.

Functional Characterization of Accessible Chromatin in Common Wheat

Eukaryotic gene transcription is fine-tuned by precise spatiotemporal interactions between cis-regulatory elements (CREs) and trans-acting factors. However, how CREs individually or coordinated with epigenetic marks function in regulating homoeolog bias expression is still largely unknown in wheat. In this study, through comprehensively characterizing open chromatin coupled with DNA methylation in the seedling and spikelet of common wheat, we observed that differential chromatin openness occurred between the seedling and spikelet, which plays important roles in tissue development through regulating the expression of related genes or through the transcription factor (TF)-centered regulatory network. Moreover, we found that CHH methylation may act as a key determinant affecting the differential binding of TFs, thereby resulting in differential expression of target genes. In addition, we found that sequence variations in MNase hypersensitive sites (MHSs) result in the differential expression of key genes responsible for important agronomic traits. Thus, our study provides new insights into the roles of CREs in regulating tissue or homoeolog bias expression, and controlling important agronomic traits in common wheat. It also provides potential CREs for genetic and epigenetic manipulation toward improving desirable traits for wheat molecule breeding.

Voice from both sides: a molecular dialogue between transcriptional activators and repressors in seed-to-seedling transition and crop adaptation

High-quality seeds provide valuable nutrients to human society and ensure successful seedling establishment. During maturation, seeds accumulate storage compounds that are required to sustain seedling growth during germination. This review focuses on the epigenetic repression of the embryonic and seed maturation programs in seedlings. We begin with an extensive overview of mutants affecting these processes, illustrating the roles of core proteins and accessory components in the epigenetic machinery by comparing mutants at both phenotypic and molecular levels. We highlight how omics assays help uncover target-specific functional specialization and coordination among various epigenetic mechanisms. Furthermore, we provide an in-depth discussion on the Seed dormancy 4 (Sdr4) transcriptional corepressor family, comparing and contrasting their regulation of seed germination in the dicotyledonous species Arabidopsis and two monocotyledonous crops, rice and wheat. Finally, we compare the similarities in the activation and repression of the embryonic and seed maturation programs through a shared set of cis-regulatory elements and discuss the challenges in applying knowledge largely gained in model species to crops.

Recent advances in exploring transcriptional regulatory landscape of crops

Crop breeding entails developing and selecting plant varieties with improved agronomic traits. Modern molecular techniques, such as genome editing, enable more efficient manipulation of plant phenotype by altering the expression of particular regulatory or functional genes. Hence, it is essential to thoroughly comprehend the transcriptional regulatory mechanisms that underpin these traits. In the multi-omics era, a large amount of omics data has been generated for diverse crop species, including genomics, epigenomics, transcriptomics, proteomics, and single-cell omics. The abundant data resources and the emergence of advanced computational tools offer unprecedented opportunities for obtaining a holistic view and profound understanding of the regulatory processes linked to desirable traits. This review focuses on integrated network approaches that utilize multi-omics data to investigate gene expression regulation. Various types of regulatory networks and their inference methods are discussed, focusing on recent advancements in crop plants. The integration of multi-omics data has been proven to be crucial for the construction of high-confidence regulatory networks. With the refinement of these methodologies, they will significantly enhance crop breeding efforts and contribute to global food security.

Starch and storage protein dynamics in the developing and matured grains of durum wheat and diploid progenitor species

Durum wheat, less immunogenically intolerant than bread wheat, originates from diploid progenitors known for nutritional quality and stress tolerance. Present study involves the analysis of major grain parameters, viz. size, weight, sugar, starch, and protein content of Triticum durum (AABB genome) and its diploid progenitors, Triticum monococcum (AA genome) and Aegilops speltoides (BB genome). Samples were collected during 2-5 weeks after anthesis (WAA), and at maturity. The investigation revealed that T. durum displayed the maximum grain size and weight. Expression analysis of Grain Weight 2 (GW2) and Glutamine Synthase (GS2), negative and positive regulators of grain weight and size, respectively, revealed higher GW2 expression in Ae. speltoides and higher GS2 expression in T. durum. Further we explored total starch, sugar and protein content, observing higher levels of starch and sugar in durum wheat while AA genome species exhibited higher protein content dominated by the fractions of albumin/globulin. HPLC profiling revealed unique sub-fractions in all three genome species. Additionally, a comparative transcriptome analysis also corroborated with the starch and protein content in the grains. This study provides valuable insights into the genetic and biochemical distinctions among durum wheat and its diploid progenitors, offering a foundation for their nutritional composition.

Deciphering the differential expression patterns of yield-related negative regulators in hexaploid wheat cultivars and hybrids at different growth stages

BACKGROUND

Hexaploid bread wheat underwent a series of polyploidization events through interspecific hybridizations that conferred adaptive plasticity and resulted in duplication and neofunctionalization of major agronomic genes. The genetic architecture of polyploid wheat not only confers adaptive plasticity but also offers huge genetic diversity. However, the contribution of different gene copies (homeologs) encoded from different subgenomes (A, B, D) at different growth stages remained unexplored.

METHODS

In this study, hybrid of elite cultivars of wheat were developed via reciprocal crosses (cytoplasm swapping) and phenotypically evaluated. We assessed differential expression profiles of yield-related negative regulators in these cultivars and their F1 hybrids and identified various cis-regulatory signatures by employing bioinformatics tools. Furthermore, the preferential expression patterns of the syntenic triads encoded from A, B, and D subgenomes were assessed to decipher their functional redundancy at six different growth stages.

RESULTS

Hybrid progenies showed better heterosis such as up to 17% increase in the average number of grains and up to 50% increase in average thousand grains weight as compared to mid-parents. Based on the expression profiling, our results indicated significant dynamic transcriptional expression patterns, portraying the different homeolog-dominance at the same stage in the different cultivars and their hybrids. Albeit belonging to same syntenic triads, a dynamic trend was observed in the regulatory signatures of these genes that might be influencing their expression profiles.

CONCLUSION

These findings can substantially contribute and provide insights for the selective introduction of better cultivars into traditional and hybrid breeding programs which can be harnessed for the improvement of future wheat.

Reshaped DNA methylation cooperating with homoeolog-divergent expression promotes improved root traits in synthesized tetraploid wheat

Polyploidization is a major event driving plant evolution and domestication. However, how reshaped epigenetic modifications coordinate gene transcription to generate phenotypic variations during wheat polyploidization is currently elusive. Here, we profiled transcriptomes and DNA methylomes of two diploid wheat accessions (SlSl and AA) and their synthetic allotetraploid wheat line (SlSlAA), which displayed elongated root hair and improved root capability for nitrate uptake and assimilation after tetraploidization. Globally decreased DNA methylation levels with a reduced difference between subgenomes were observed in the roots of SlSlAA. DNA methylation changes in first exon showed strong connections with altered transcription during tetraploidization. Homoeolog-specific transcription was associated with biased DNA methylation as shaped by homoeologous sequence variation. The hypomethylated promoters showed significantly enriched binding sites for MYB, which may affect gene transcription in response to root hair growth. Two master regulators in root hair elongation pathway, AlCPC and TuRSL4, exhibited upregulated transcription levels accompanied by hypomethylation in promoter, which may contribute to the elongated root hair. The upregulated nitrate transporter genes, including NPFs and NRTs, also are significantly associated with hypomethylation, indicating an epigenetic-incorporated regulation manner in improving nitrogen use efficiency. Collectively, these results provided new insights into epigenetic changes in response to crop polyploidization and underscored the importance of epigenetic regulation in improving crop traits.

TabHLH95-TaNF-YB1 module promotes grain starch synthesis in bread wheat

Starch is the most abundant substance in wheat (Triticum aestivum L.) endosperm and provides the major carbohydrate energy for human daily life. Starch synthesis-related (SSR) genes are believed to be spatiotemporally specific, but their transcriptional regulation remains unclear in wheat. Here, we investigate the role of the basic helix-loop-helix (bHLH) transcription factor TabHLH95 in starch synthesis. TabHLH95 is preferentially expressed in the developing grains in wheat and encodes a nucleus localized protein without autoactivation activity. The Tabhlh95 knockout mutants display smaller grain size and less starch content than wild type, whereas overexpression of TabHLH95 enhances starch accumulation and significantly improves thousand grain weight. Transcriptome analysis reveals that the expression of multiple SSR genes is significantly reduced in the Tabhlh95 mutants. TabHLH95 binds to the promoters of ADP-glucose pyrophosphorylase large subunit 1 (AGPL1-1D/-1B), AGPL2-5D, and isoamylase (ISA1-7D) and enhances their transcription. Intriguingly, TabHLH95 interacts with the nuclear factor Y (NF-Y) family transcription factor TaNF-YB1, thereby synergistically regulating starch synthesis. These results suggest that the TabHLH95-TaNF-YB1 complex positively modulates starch synthesis and grain weight by regulating the expression of a subset of SSR genes, thus providing a good potential approach for genetic improvement of grain productivity in wheat.

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

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

Embryo Rescue in Plant Breeding

Embryo rescue (ER) techniques are among the oldest and most successful in vitro tissue culture protocols used with plant species. ER refers to a series of methods that promote the development of an immature or lethal embryo into a viable plant. Intraspecific, interspecific, or intergeneric crosses allow the introgression of important alleles of agricultural interest from wild species, such as resistance or tolerance to abiotic and biotic stresses or morphological traits in crops. However, pre-zygotic and post-zygotic reproductive barriers often present challenges in achieving successful hybridization. Pre-zygotic barriers manifest as incompatibility reactions that hinder pollen germination, pollen tube growth, or penetration into the ovule occurring in various tissues, such as the stigma, style, or ovary. To overcome these barriers, several strategies are employed, including cut-style or graft-on-style techniques, the utilization of mixed pollen from distinct species, placenta pollination, and in vitro ovule pollination. On the other hand, post-zygotic barriers act at different tissues and stages ranging from early embryo development to the subsequent growth and reproduction of the offspring. Many crosses among different genera result in embryo abortion due to the failure of endosperm development. In such cases, ER techniques are needed to rescue these hybrids. ER holds great promise for not only facilitating successful crosses but also for obtaining haploids, doubled haploids, and manipulating the ploidy levels for chromosome engineering by monosomic and disomic addition as well substitution lines. Furthermore, ER can be used to shorten the reproductive cycle and for the propagation of rare plants. Additionally, it has been repeatedly used to study the stages of embryonic development, especially in embryo-lethal mutants. The most widely used ER procedure is the culture of immature embryos taken and placed directly on culture media. In certain cases, the in vitro culture of ovule, ovaries or placentas enables the successful development of young embryos from the zygote stage to maturity.

The translational landscape of bread wheat during grain development

The dynamics of gene expression in crop grains has typically been investigated at the transcriptional level. However, this approach neglects translational regulation, a widespread mechanism that rapidly modulates gene expression to increase the plasticity of organisms. Here, we performed ribosome profiling and polysome profiling to obtain a comprehensive translatome data set of developing bread wheat (Triticum aestivum) grains. We further investigated the genome-wide translational dynamics during grain development, revealing that the translation of many functional genes is modulated in a stage-specific manner. The unbalanced translation between subgenomes is pervasive, which increases the expression flexibility of allohexaploid wheat. In addition, we uncovered widespread previously unannotated translation events, including upstream open reading frames (uORFs), downstream open reading frames (dORFs), and open reading frames (ORFs) in long noncoding RNAs, and characterized the temporal expression dynamics of small ORFs. We demonstrated that uORFs act as cis-regulatory elements that can repress or even enhance the translation of mRNAs. Gene translation may be combinatorially modulated by uORFs, dORFs, and microRNAs. In summary, our study presents a translatomic resource that provides a comprehensive and detailed overview of the translational regulation in developing bread wheat grains. This resource will facilitate future crop improvements for optimal yield and quality.

A multi-omic resource of wheat seed tissues for nutrient deposition and improvement for human health

As a globally important staple crop, wheat seeds provide us with nutrients and proteins. The trend of healthy dietary has become popular recently, emphasizing the consumption of whole-grain wheat products and the dietary benefits. However, the dynamic changes in nutritional profiles of different wheat seed regions (i.e., the embryo, endosperm and outer layers) during developmental stages and the molecular regulation have not been well studied. Here, we provide this multi-omic resource of wheat seeds and describe the generation, technical assessment and preliminary analyses. This resource includes a time-series RNA-seq dataset of the embryo, endosperm and outer layers of wheat seeds and their corresponding metabolomic dataset, covering the middle and late stages of seed development. Our RNA-seq experiments profile the expression of 63,708 genes, while the metabolomic data includes the abundance of 984 metabolites. We believe that this was the first reported transcriptome and metabolome dataset of wheat seeds that helps understand the molecular regulation of the deposition of beneficial nutrients and hence improvements for nutritional and processing quality traits.

Spatiotemporal transcriptomics and metabolic profiling provide insights into gene regulatory networks during lentil seed development

Lentil (Lens culinaris Medik.) is a nutritious legume with seeds rich in protein, minerals and an array of diverse specialized metabolites. The formation of a seed requires regulation and tight coordination of developmental programs to form the embryo, endosperm and seed coat compartments, which determines the structure and composition of mature seed and thus its end-use quality. Understanding the molecular and cellular events and metabolic processes of seed development is essential for improving lentil yield and seed nutritional value. However, such information remains largely unknown, especially at the seed compartment level. In this study, we generated high-resolution spatiotemporal gene expression profiles in lentil embryo, seed coat and whole seeds from fertilization through maturation. Apart from anatomic differences between the embryo and seed coat, comparative transcriptomics and weighted gene co-expression network analysis revealed embryo- and seed coat-specific genes and gene modules predominant in specific tissues and stages, which highlights distinct genetic programming. Furthermore, we investigated the dynamic profiles of flavonoid, isoflavone, phytic acid and saponin in seed compartments across seed development. Coupled with transcriptome data, we identified sets of candidate genes involved in the biosynthesis of these metabolites. The global view of the transcriptional and metabolic changes of lentil seed tissues throughout development provides a valuable resource for dissecting the genetic control of secondary metabolism and development of molecular tools for improving seed nutritional quality.

Chromatin accessibility landscapes revealed the subgenome-divergent regulation networks during wheat grain development

Development of wheat (Triticum aestivum L.) grain mainly depends on the processes of starch synthesis and storage protein accumulation, which are critical for grain yield and quality. However, the regulatory network underlying the transcriptional and physiological changes of grain development is still not clear. Here, we combined ATAC-seq and RNA-seq to discover the chromatin accessibility and gene expression dynamics during these processes. We found that the chromatin accessibility changes are tightly associated with differential transcriptomic expressions, and the proportion of distal ACRs was increased gradually during grain development. Specific transcription factor (TF) binding sites were enriched at different stages and were diversified among the 3 subgenomes. We further predicted the potential interactions between key TFs and genes related with starch and storage protein biosynthesis and found different copies of some key TFs played diversified roles. Overall, our findings have provided numerous resources and illustrated the regulatory network during wheat grain development, which would shed light on the improvement of wheat yields and qualities.

Supplementary Information

The online version contains supplementary material available at 10.1007/s42994-023-00095-8.

New Growth-Related Features of Wheat Grain Pericarp Revealed by Synchrotron-Based X-ray Micro-Tomography and 3D Reconstruction

Wheat (Triticum aestivum L.) is one of the most important crops as it provides 20% of calories and proteins to the human population. To overcome the increasing demand in wheat grain production, there is a need for a higher grain yield, and this can be achieved in particular through an increase in the grain weight. Moreover, grain shape is an important trait regarding the milling performance. Both the final grain weight and shape would benefit from a comprehensive knowledge of the morphological and anatomical determinism of wheat grain growth. Synchrotron-based phase-contrast X-ray microtomography (X-ray µCT) was used to study the 3D anatomy of the growing wheat grain during the first developmental stages. Coupled with 3D reconstruction, this method revealed changes in the grain shape and new cellular features. The study focused on a particular tissue, the pericarp, which has been hypothesized to be involved in the control of grain development. We showed considerable spatio-temporal diversity in cell shape and orientations, and in tissue porosity associated with stomata detection. These results highlight the growth-related features rarely studied in cereal grains, which may contribute significantly to the final grain weight and shape.

New cup out of old coffee: contribution of parental gene expression legacy to phenotypic novelty in coffee beans of the allopolyploid Coffea arabica L

BACKGROUND AND AIMS

Allopolyploidization is a widespread phenomenon known to generate novel phenotypes by merging evolutionarily distinct parental genomes and regulatory networks in a single nucleus. The objective of this study was to investigate the transcriptional regulation associated with phenotypic novelty in coffee beans of the allotetraploid Coffea arabica.

METHODS

A genome-wide comparative transcriptomic analysis was performed in C. arabica and its two diploid progenitors, C. canephora and C. eugenioides. Gene expression patterns and homeologue expression were studied on seeds at five different maturation stages. The involvement of homeologue expression bias (HEB) in specific traits was addressed both by functional enrichment analyses and by the study of gene expression in the caffeine and chlorogenic acid biosynthesis pathways.

KEY RESULTS

Expression-level dominance in C. arabica seed was observed for most of the genes differentially expressed between the species. Approximately a third of the genes analysed showed HEB. This proportion increased during seed maturation but the biases remained equally distributed between the sub-genomes. The relative expression levels of homeologues remained relatively constant during maturation and were correlated with those estimated in leaves of C. arabica and interspecific hybrids between C. canephora and C. eugenioides. Functional enrichment analyses performed on genes exhibiting HEB enabled the identification of processes potentially associated with physiological traits. The expression profiles of the genes involved in caffeine biosynthesis mirror the differences observed in the caffeine content of mature seeds of C. arabica and its parental species.

CONCLUSIONS

Neither of the two sub-genomes is globally preferentially expressed in C. arabica seeds, and homeologues appear to be co-regulated by shared trans-regulatory mechanisms. The observed HEBs are thought to be a legacy of gene expression differences inherited from diploid progenitor species. Pre-existing functional divergences between parental species appear to play an important role in controlling the phenotype of C. arabica seeds.

A wheat integrative regulatory network from large-scale complementary functional datasets enables trait-associated gene discovery for crop improvement

Gene regulation is central to all aspects of organism growth, and understanding it using large-scale functional datasets can provide a whole view of biological processes controlling complex phenotypic traits in crops. However, the connection between massive functional datasets and trait-associated gene discovery for crop improvement is still lacking. In this study, we constructed a wheat integrative gene regulatory network (wGRN) by combining an updated genome annotation and diverse complementary functional datasets, including gene expression, sequence motif, transcription factor (TF) binding, chromatin accessibility, and evolutionarily conserved regulation. wGRN contains 7.2 million genome-wide interactions covering 5947 TFs and 127 439 target genes, which were further verified using known regulatory relationships, condition-specific expression, gene functional information, and experiments. We used wGRN to assign genome-wide genes to 3891 specific biological pathways and accurately prioritize candidate genes associated with complex phenotypic traits in genome-wide association studies. In addition, wGRN was used to enhance the interpretation of a spike temporal transcriptome dataset to construct high-resolution networks. We further unveiled novel regulators that enhance the power of spike phenotypic trait prediction using machine learning and contribute to the spike phenotypic differences among modern wheat accessions. Finally, we developed an interactive webserver, wGRN (http://wheat.cau.edu.cn/wGRN), for the community to explore gene regulation and discover trait-associated genes. Collectively, this community resource establishes the foundation for using large-scale functional datasets to guide trait-associated gene discovery for crop improvement.

Beyond Transcript Concentrations: Quantifying Polyploid Expression Responses per Biomass, per Genome, and per Cell with RNA-Seq

RNA-seq has been used extensively to study expression responses to polyploidy. Most current methods for normalizing RNA-seq data yield estimates of transcript concentrations (transcripts per transcriptome). The implicit assumption of these normalization methods is that transcriptome size is equivalent between the samples being compared such that transcript concentrations are equivalent to transcripts per cell. In recent years, however, evidence has mounted that transcriptome size can vary dramatically in response to a range of factors including polyploidy and that such variation is ubiquitous. Where such variation exists, transcript concentration is often a poor or even misleading proxy for expression responses at other biologically relevant scales (e.g., expression per cell). Thus, it is important that transcriptomic studies of polyploids move beyond simply comparing transcript concentrations if we are to gain a complete understanding of how genome multiplication affects gene expression. I discuss this issue in more detail and summarize a suite of approaches that can leverage RNA-seq to quantify expression responses per genome, per cell, and per unit of biomass.

Toward kingdom-wide analyses of gene expression

Gene expression data for Archaeplastida are accumulating exponentially, with more than 300 000 RNA-sequencing (RNA-seq) experiments available for hundreds of species. The gene expression data stem from thousands of experiments that capture gene expression in various organs, tissues, cell types, (a)biotic perturbations, and genotypes. Advances in software tools make it possible to process all these data in a matter of weeks on modern office computers, giving us the possibility to study gene expression in a kingdom-wide manner for the first time. We discuss how the expression data can be accessed and processed and outline analyses that take advantage of cross-species analyses, allowing us to generate powerful and robust hypotheses about gene function and evolution.

Dynamic chromatin regulatory programs during embryogenesis of hexaploid wheat

BACKGROUND

Plant and animal embryogenesis have conserved and distinct features. Cell fate transitions occur during embryogenesis in both plants and animals. The epigenomic processes regulating plant embryogenesis remain largely elusive.

RESULTS

Here, we elucidate chromatin and transcriptomic dynamics during embryogenesis of the most cultivated crop, hexaploid wheat. Time-series analysis reveals stage-specific and proximal-distal distinct chromatin accessibility and dynamics concordant with transcriptome changes. Following fertilization, the remodeling kinetics of H3K4me3, H3K27ac, and H3K27me3 differ from that in mammals, highlighting considerable species-specific epigenomic dynamics during zygotic genome activation. Polycomb repressive complex 2 (PRC2)-mediated H3K27me3 deposition is important for embryo establishment. Later H3K27ac, H3K27me3, and chromatin accessibility undergo dramatic remodeling to establish a permissive chromatin environment facilitating the access of transcription factors to cis-elements for fate patterning. Embryonic maturation is characterized by increasing H3K27me3 and decreasing chromatin accessibility, which likely participates in restricting totipotency while preventing extensive organogenesis. Finally, epigenomic signatures are correlated with biased expression among homeolog triads and divergent expression after polyploidization, revealing an epigenomic contributor to subgenome diversification in an allohexaploid genome.

CONCLUSIONS

Collectively, we present an invaluable resource for comparative and mechanistic analysis of the epigenomic regulation of crop embryogenesis.

Differential Expression Feature Extraction (DEFE): A Case Study in Wheat FHB RNA-Seq Data Analysis

In differential gene expression data analysis, one objective is to identify groups of co-expressed genes from a large dataset in order to detect the association between such a group of genes and an experimental condition. This is often done through a clustering approach, such as k-means or bipartition hierarchical clustering, based on particular similarity measures in the grouping process. In such a dataset, the gene differential expression itself is an innate attribute that can be used in the feature extraction process. For example, in a dataset consisting of multiple treatments versus their controls, the expression of a gene in each treatment would have three possible behaviors, upregulated, downregulated, or unchanged. We present in this chapter, a differential expression feature extraction (DEFE) method by using a string consisting of three numerical values at each character to denote such behavior, i.e., 1 = up, 2 = down, and 0 = unchanged, which results in up to 3^B differential expression patterns across all B comparisons. This approach has been successfully applied in many research projects, and among these, we demonstrate the strength of DEFE in a case study on RNA-sequencing (RNA-seq) data analysis of wheat challenged with the phytopathogenic fungus, Fusarium graminearum. Combinations of multiple schemes of DEFE patterns revealed groups of genes putatively associated with resistance or susceptibility to FHB.

Asymmetric gene expression in grain development of reciprocal crosses between tetraploid and hexaploid wheats

Production of viable progeny from interploid crosses requires precise regulation of gene expression from maternal and paternal chromosomes, yet the transcripts contributed to hybrid seeds from polyploid parent species have rarely been explored. To investigate the genome-wide maternal and paternal contributions to polyploid grain development, we analyzed the transcriptomes of developing embryos, from zygote to maturity, alongside endosperm in two stages of development, using reciprocal crosses between tetraploid and hexaploid wheats. Reciprocal crosses between species with varied levels of ploidy displayed broad impacts on gene expression, including shifts in alternative splicing events in select crosses, as illustrated by active splicing events, enhanced protein synthesis and chromatin remodeling. Homoeologous gene expression was repressed on the univalent D genome in pentaploids, but this suppression was attenuated in crosses with a higher ploidy maternal parent. Imprinted genes were identified in endosperm and early embryo tissues, supporting predominant maternal effects on early embryogenesis. By systematically investigating the complex transcriptional networks in reciprocal-cross hybrids, this study presents a framework for understanding the genomic incompatibility and transcriptome shock that results from interspecific hybridization and uncovers the transcriptional impacts on hybrid seeds created from agriculturally-relevant polyploid species.

Transcriptome Analysis Reveals Potential Mechanism in Storage Protein Trafficking within Developing Grains of Common Wheat

Gluten proteins are the major storage protein fraction in the mature wheat grain. They are restricted to the starchy endosperm, which defines the viscoelastic properties of wheat dough. The synthesis of these storage proteins is controlled by the endoplasmic reticulum (ER) and is directed into the vacuole via the Golgi apparatus. In the present study, transcriptome analysis was used to explore the potential mechanism within critical stages of grain development of wheat cultivar "Shaannong 33" and its sister line used as the control (CK). Samples were collected at 10 DPA (days after anthesis), 14 DPA, 20 DPA, and 30 DPA for transcriptomic analysis. The comparative transcriptome analysis identified that a total of 18,875 genes were differentially expressed genes (DEGs) between grains of four groups "T10 vs. CK10, T14 vs. CK14, T20 vs. CK20, and T30 vs. CK30", including 2824 up-regulated and 5423 down-regulated genes in T30 vs. CK30. Further, the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment highlighted the maximum number of genes regulating protein processing in the endoplasmic reticulum (ER) during grain enlargement stages (10-20 DPA). In addition, KEGG database analysis reported 1362 and 788 DEGs involved in translation, ribosomal structure, biogenesis, flavonoid biosynthesis pathway and intracellular trafficking, secretion, and vesicular transport through protein processing within ER pathway (ko04141). Notably, consistent with the higher expression of intercellular storage protein trafficking genes at the initial 10 DPA, there was relatively low expression at later stages. Expression levels of nine randomly selected genes were verified by qRT-PCR, which were consistent with the transcriptome data. These data suggested that the initial stages of "cell division" played a significant role in protein quality control within the ER, thus maintaining the protein quality characteristics at grain maturity. Furthermore, our data suggested that the protein synthesis, folding, and trafficking pathways directed by a different number of genes during the grain enlargement stage contributed to the observed high-quality characteristics of gluten protein in Shaannong 33 (Triticum aestivum L.).

Chromosome painting reveals inter-chromosomal rearrangements and evolution of subgenome D of wheat

Aegilops species represent the most important gene pool for breeding bread wheat (Triticum aestivum). Thus, understanding the genome evolution, including chromosomal structural rearrangements and syntenic relationships among Aegilops species or between Aegilops and wheat, is important for both basic genome research and practical breeding applications. In the present study, we attempted to develop subgenome D-specific fluorescence in situ hybridization (FISH) probes by selecting D-specific oligonucleotides based on the reference genome of Chinese Spring. The oligo-based chromosome painting probes consisted of approximately 26 000 oligos per chromosome and their specificity was confirmed in both diploid and polyploid species containing the D subgenome. Two previously reported translocations involving two D chromosomes have been confirmed in wheat varieties and their derived lines. We demonstrate that the oligo painting probes can be used not only to identify the translocations involving D subgenome chromosomes, but also to determine the precise positions of chromosomal breakpoints. Chromosome painting of 56 accessions of Ae. tauschii from different origins led us to identify two novel translocations: a reciprocal 3D-7D translocation in two accessions and a complex 4D-5D-7D translocation in one accession. Painting probes were also used to analyze chromosomes from more diverse Aegilops species. These probes produced FISH signals in four different genomes. Chromosome rearrangements were identified in Aegilops umbellulata, Aegilops markgrafii, and Aegilops uniaristata, thus providing syntenic information that will be valuable for the application of these wild species in wheat breeding.

Wheat genomic study for genetic improvement of traits in China

Bread wheat (Triticum aestivum L.) is a major crop that feeds 40% of the world's population. Over the past several decades, advances in genomics have led to tremendous achievements in understanding the origin and domestication of wheat, and the genetic basis of agronomically important traits, which promote the breeding of elite varieties. In this review, we focus on progress that has been made in genomic research and genetic improvement of traits such as grain yield, end-use traits, flowering regulation, nutrient use efficiency, and biotic and abiotic stress responses, and various breeding strategies that contributed mainly by Chinese scientists. Functional genomic research in wheat is entering a new era with the availability of multiple reference wheat genome assemblies and the development of cutting-edge technologies such as precise genome editing tools, high-throughput phenotyping platforms, sequencing-based cloning strategies, high-efficiency genetic transformation systems, and speed-breeding facilities. These insights will further extend our understanding of the molecular mechanisms and regulatory networks underlying agronomic traits and facilitate the breeding process, ultimately contributing to more sustainable agriculture in China and throughout the world.

Spatiotemporal Transcriptomic Atlas of Developing Embryos and Vegetative Tissues in Flax

Flax (Linum usitatissimum L.) is an important multipurpose crop widely grown for oil and fiber. Despite recent advances in genomics, detailed gene activities during the important reproductive phase of its development are not well defined. In this study, we employed high-throughput RNA-sequencing methods to generate in-depth transcriptome profiles of flax tissues with emphasis on the reproductive phases of five key stages of embryogenesis (globular embryo, heart embryo, torpedo embryo, cotyledon embryo, and mature embryo), mature seed, and vegetative tissues viz. ovary, anther, and root. These datasets were used to establish the co-expression networks covering 36 gene modules based on the expression patterns for each gene through weighted gene co-expression network analysis (WGCNA). Functional interrogation with Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) of dominantly expressed genetic modules in tissues revealed pathways involved in the development of different tissues. Moreover, the essential genes in embryo development and synthesis of storage reserves were identified based on their dynamic expression patterns. Together, this comprehensive dataset for developing embryos, mature seeds and vegetative tissues provides new insights into molecular mechanisms of seed development with potential for flax crop improvement.

Spatio-temporal transcriptome profiling and subgenome analysis in Brassica napus

Brassica napus is an important oil crop and an allotetraploid species. However, the detailed analysis of gene function and homoeologous gene expression in all tissues at different developmental stages was not explored. In this study, we performed a global transcriptome analysis of 24 vegetative and reproductive tissues at six developmental stages (totally 111 tissues). These samples were clustered into eight groups. The gene functions of silique pericarp were similar to roots, stems and leaves. In particular, glucosinolate metabolic process was associated with root and silique pericarp. Genes involved in protein phosphorylation were often associated with stamen, anther and the early developmental stage of seeds. Transcription factor (TF) genes were more specific than structural genes. A total of 17 100 genes that were preferentially expressed in one tissue (tissue-preferred genes, TPGs), including 889 TFs (5.2%), were identified in the 24 tissues. Some TPGs were identified as hub genes in the co-expression network analysis, and some TPGs in different tissues were involved in different hormone pathways. About 67.0% of the homoeologs showed balanced expression, whereas biased expression of homoeologs was associated with structural divergence. In addition, the spatiotemporal expression of homoeologs was related to the presence of transposable elements (TEs) and regulatory elements (REs); more TEs and fewer REs in the promoters resulted in divergent expression in different tissues. This study provides a valuable transcriptional map for understanding the growth and development of B. napus, for identifying important genes for future crop improvement, and for exploring gene expression patterns in the B. napus.

Cloning and Functional Analysis of TaWRI1Ls, the Key Genes for Grain Fatty Acid Synthesis in Bread Wheat

WRINKLED1 (WRI1), an APETALA2 (AP2) transcription factor (TF), critically regulates the processes related to fatty acid synthesis, storage oil accumulation, and seed development in plants. However, the WRI1 genes remain unknown in allohexaploid bread wheat (Triticum aestivum L.). In this study, based on the sequence of Arabidopsis AtWRI1, two TaWRI1Ls genes of bread wheat, TaWRI1L1 and TaWRI1L2, were cloned. TaWRI1L2 was closely related to monocotyledons and clustered in one subgroup with AtWRI1, while TaWRI1L1 was clustered in another subgroup with AtWRI3 and AtWRI4. Both were expressed highly in the developmental grain, subcellular localized in the nucleus, and showed transcriptional activation activity. TaWRI1L2, rather than TaWRI1L1, promoted oil body accumulation and significantly increased triglyceride (TAG) content in tobacco leaves. Overexpression of TaWRI1L2 compensated for the functional loss of AtWRI1 in an Arabidopsis mutant and restored the wild-type phenotypes of seed shape, generation, and fatty acid synthesis and accumulation. Knockout of TaWRI1L2 reduced grain size, 1000 grain weight, and grain fatty acid synthesis in bread wheat. Conclusively, TaWRI1L2, rather than TaWRI1L1, was the key transcriptional factor in the regulation of grain fatty acid synthesis in bread wheat. This study lays a foundation for gene regulation and genetic manipulation of fatty acid synthesis in wheat genetic breeding programs.

Evolutionary divergence in embryo and seed coat development of U's Triangle Brassica species illustrated by a spatiotemporal transcriptome atlas

The economically valuable Brassica species include the six related members of U's Triangle. Despite the agronomic and economic importance of these Brassicas, the impacts of evolution and relatively recent domestication events on the genetic landscape of seed development have not been comprehensively examined in these species. Here we present a 3D transcriptome atlas for the six species of U's Triangle, producing a unique resource that captures gene expression data for the major subcompartments of the seed, from the unfertilized ovule to the mature embryo and seed coat. This comprehensive dataset for seed development in tetraploid and ancestral diploid Brassicas provides new insights into evolutionary divergence and expression bias at the gene and subgenome levels during the domestication of these valued crop species. Comparisons of gene expression associated with regulatory networks and metabolic pathways operating in the embryo and seed coat during seed development reveal differences in storage reserve accumulation and fatty acid metabolism among the six Brassica species. This study illustrates the genetic underpinnings of seed traits and the selective pressures placed on seed production, providing an immense resource for continued investigation of Brassica polyploid biology, genomics and evolution.

Genome-wide identification of seed storage protein gene regulators in wheat through coexpression analysis

Only a few transcriptional regulators of seed storage protein (SSP) genes have been identified in common wheat (Triticum aestivum L.). Coexpression analysis could be an efficient approach to characterize novel transcriptional regulators at the genome-scale considering the correlated expression between transcriptional regulators and target genes. As the A genome donor of common wheat, Triticum urartu is more suitable for coexpression analysis than common wheat considering the diploid genome and single gene copy. In this work, the transcriptome dynamics in endosperm of T. urartu throughout grain filling were revealed by RNA-Seq analysis. In the coexpression analysis, a total of 71 transcription factors (TFs) from 23 families were found to be coexpressed with SSP genes. Among these TFs, TuNAC77 enhanced the transcription of SSP genes by binding to cis-elements distributed in promoters. The homolog of TuNAC77 in common wheat, TaNAC77, shared an identical function, and the total SSPs were reduced by about 24% in common wheat when TaNAC77 was knocked down. This is the first genome-wide identification of transcriptional regulators of SSP genes in wheat, and the newly characterized transcriptional regulators will undoubtedly expand our knowledge of the transcriptional regulation of SSP synthesis.

GWAS identifies an ortholog of the rice D11 gene as a candidate gene for grain size in an international collection of hexaploid wheat

Grain size is a key agronomic trait that contributes to grain yield in hexaploid wheat. Grain length and width were evaluated in an international collection of 157 wheat accessions. These accessions were genetically characterized using a genotyping-by-sequencing (GBS) protocol that produced 73,784 single nucleotide polymorphism (SNP) markers. GBS-derived genotype calls obtained on Chinese Spring proved extremely accurate when compared to the reference (> 99.9%) and showed > 95% agreement with calls made at SNP loci shared with the 90 K SNP array on a subset of 71 Canadian wheat accessions for which both types of data were available. This indicates that GBS can yield a large amount of highly accurate SNP data in hexaploid wheat. The genetic diversity analysis performed using this set of SNP markers revealed the presence of six distinct groups within this collection. A GWAS was conducted to uncover genomic regions controlling variation for grain length and width. In total, seven SNPs were found to be associated with one or both traits, identifying three quantitative trait loci (QTLs) located on chromosomes 1D, 2D and 4A. In the vicinity of the peak SNP on chromosome 2D, we found a promising candidate gene (TraesCS2D01G331100), whose rice ortholog (D11) had previously been reported to be involved in the regulation of grain size. These markers will be useful in breeding for enhanced wheat productivity.

Conserved, divergent and heterochronic gene expression during Brachypodium and Arabidopsis embryo development

KEY MESSAGE

Developmental and transcriptomic analysis of Brachypodium embryogenesis and comparison with Arabidopsis identifies conserved and divergent phases of embryogenesis and reveals widespread heterochrony of developmental gene expression. Embryogenesis, transforming the zygote into the mature embryo, represents a fundamental process for all flowering plants. Current knowledge of cell specification and differentiation during plant embryogenesis is largely based on studies of the dicot model plant Arabidopsis thaliana. However, the major crops are monocots and the transcriptional programs associated with the differentiation processes during embryogenesis in this clade were largely unknown. Here, we combined analysis of cell division patterns with development of a temporal transcriptomic resource during embryogenesis of the monocot model plant Brachypodium distachyon. We found that early divisions of the Brachypodium embryo were highly regular, while later stages were marked by less stereotypic patterns. Comparative transcriptomic analysis between Brachypodium and Arabidopsis revealed that early and late embryogenesis shared a common transcriptional program, whereas mid-embryogenesis was divergent between species. Analysis of orthology groups revealed widespread heterochronic expression of potential developmental regulators between the species. Interestingly, Brachypodium genes tend to be expressed at earlier stages than Arabidopsis counterparts, which suggests that embryo patterning may occur early during Brachypodium embryogenesis. Detailed investigation of auxin-related genes shows that the capacity to synthesize, transport and respond to auxin is established early in the embryo. However, while early PIN1 polarity could be confirmed, it is unclear if an active response is mounted. This study presents a resource for studying Brachypodium and grass embryogenesis and shows that divergent angiosperms share a conserved genetic program that is marked by heterochronic gene expression.

Large-scale integration of meta-QTL and genome-wide association study discovers the genomic regions and candidate genes for yield and yield-related traits in bread wheat

Based on the large-scale integration of meta-QTL and Genome-Wide  Association Study, 76 high-confidence MQTL regions and 237 candidate genes that affected wheat yield and yield-related traits were discovered. Improving yield and yield-related traits are key goals in wheat breeding program. The integration of accumulated wheat genetic resources provides an opportunity to uncover important genomic regions and candidate genes that affect wheat yield. Here, a comprehensive meta-QTL analysis was conducted on 2230 QTL of yield-related traits obtained from 119 QTL studies. These QTL were refined into 145 meta-QTL (MQTL), and 89 MQTL were verified by GWAS with different natural populations. The average confidence interval (CI) of these MQTL was 2.92 times less than that of the initial QTL. Furthermore, 76 core MQTL regions with a physical distance less than 25 Mb were detected. Based on the homology analysis and expression patterns, 237 candidate genes in the MQTL involved in photoperiod response, grain development, multiple plant growth regulator pathways, carbon and nitrogen metabolism and spike and flower organ development were determined. A novel candidate gene TaKAO-4A was confirmed to be significantly associated with grain size, and a CAPS marker was developed based on its dominant haplotype. In summary, this study clarified a method based on the integration of meta-QTL, GWAS and homology comparison to reveal the genomic regions and candidate genes that affect important yield-related traits in wheat. This work will help to lay a foundation for the identification, transfer and aggregation of these important QTL or candidate genes in wheat high-yield breeding.

RNA-Seq-based DNA marker analysis of the genetics and molecular evolution of Triticeae species

The release of high-quality chromosome-level genome sequences of members of the Triticeae tribe has greatly facilitated genetic and genomic analyses of important crops such as wheat (Triticum aestivum) and barley (Hordeum vulgare). Due to the large diploid genome size of Triticeae plants (ca. 5 Gbp), transcript analysis is an important method for identifying genetic and genomic differences among Triticeae species. In this review, we summarize our results of RNA-Seq analyses of diploid wheat accessions belonging to the genera Aegilops and Triticum. We also describe studies of the molecular relationships among these accessions and provide insight into the evolution of common hexaploid wheat. DNA markers based on polymorphisms within species can be used to map loci of interest. Even though the genome sequence of diploid Aegilops tauschii, the D-genome donor of common wheat, has been released, the diploid barley genome continues to provide key information about the physical structures of diploid wheat genomes. We describe how a series of RNA-Seq analyses of wheat relatives has helped uncover the structural and evolutionary features of genomic and genetic systems in wild and cultivated Triticeae species.

Alternative splicing dynamics and evolutionary divergence during embryogenesis in wheat species

Among polyploid species with complex genomic architecture, variations in the regulation of alternative splicing (AS) provide opportunities for transcriptional and proteomic plasticity and the potential for generating trait diversities. However, the evolution of AS and its influence on grain development in diploid grass and valuable polyploid wheat crops are poorly understood. To address this knowledge gap, we developed a pipeline for the analysis of alternatively spliced transcript isoforms, which takes the high sequence similarity among polyploid wheat subgenomes into account. Through analysis of synteny and detection of collinearity of homoeologous subgenomes, conserved and specific AS events across five wheat and grass species were identified. A global analysis of the regulation of AS in diploid grass and polyploid wheat grains revealed diversity in AS events not only between the endosperm, pericarp and embryo overdevelopment, but also between subgenomes. Analysis of AS in homoeologous triads of polyploid wheats revealed evolutionary divergence between gene-level and transcript-level regulation of embryogenesis. Evolutionary age analysis indicated that the generation of novel transcript isoforms has occurred in young genes at a more rapid rate than in ancient genes. These findings, together with the development of comprehensive AS resources for wheat and grass species, advance understanding of the evolution of regulatory features of AS during embryogenesis and grain development in wheat.

Exploring the genic resources underlying metabolites through mGWAS and mQTL in wheat: From large-scale gene identification and pathway elucidation to crop improvement

Common wheat (Triticum aestivum L.) is a leading cereal crop, but has lagged behind with respect to the interpretation of the molecular mechanisms of phenotypes compared with other major cereal crops such as rice and maize. The recently available genome sequence of wheat affords the pre-requisite information for efficiently exploiting the potential molecular resources for decoding the genetic architecture of complex traits and identifying valuable breeding targets. Meanwhile, the successful application of metabolomics as an emergent large-scale profiling methodology in several species has demonstrated this approach to be accessible for reaching the above goals. One such productive avenue is combining metabolomics approaches with genetic designs. However, this trial is not as widespread as that for sequencing technologies, especially when the acquisition, understanding, and application of metabolic approaches in wheat populations remain more difficult and even arguably underutilized. In this review, we briefly introduce the techniques used in the acquisition of metabolomics data and their utility in large-scale identification of functional candidate genes. Considerable progress has been made in delivering improved varieties, suggesting that the inclusion of information concerning these metabolites and genes and metabolic pathways enables a more explicit understanding of phenotypic traits and, as such, this procedure could serve as an -omics-informed roadmap for executing similar improvement strategies in wheat and other species.

Reduced free asparagine in wheat grain resulting from a natural deletion of TaASN-B2: investigating and exploiting diversity in the asparagine synthetase gene family to improve wheat quality

BACKGROUND

Understanding the determinants of free asparagine concentration in wheat grain is necessary to reduce levels of the processing contaminant acrylamide in baked and toasted wheat products. Although crop management strategies can help reduce asparagine concentrations, breeders have limited options to select for genetic variation underlying this trait. Asparagine synthetase enzymes catalyse a critical step in asparagine biosynthesis in plants and, in wheat, are encoded by five homeologous gene triads that exhibit distinct expression profiles. Within this family, TaASN2 genes are highly expressed during grain development but TaASN-B2 is absent in some varieties.

RESULTS

Natural genetic diversity in the asparagine synthetase gene family was assessed in different wheat varieties revealing instances of presence/absence variation and other polymorphisms, including some predicted to affect the function of the encoded protein. The presence and absence of TaASN-B2 was determined across a range of UK and global common wheat varieties and related species, showing that the deletion encompassing this gene was already present in some wild emmer wheat genotypes. Expression profiling confirmed that TaASN2 transcripts were only detectable in the grain, while TaASN3.1 genes were highly expressed during the early stages of grain development. TaASN-A2 was the most highly expressed TaASN2 homeologue in most assayed wheat varieties. TaASN-B2 and TaASN-D2 were expressed at similar, lower levels in varieties possessing TaASN-B2. Expression of TaASN-A2 and TaASN-D2 did not increase to compensate for the absence of TaASN-B2, so total TaASN2 expression was lower in varieties lacking TaASN-B2. Consequently, free asparagine concentrations in field-produced grain were, on average, lower in varieties lacking TaASN-B2, although the effect was lost when free asparagine accumulated to very high concentrations as a result of sulphur deficiency.

CONCLUSIONS

Selecting wheat genotypes lacking the TaASN-B2 gene may be a simple and rapid way for breeders to reduce free asparagine concentrations in commercial wheat grain.

STARCH SYNTHASE 4 is required for normal starch granule initiation in amyloplasts of wheat endosperm

Starch granule initiation is poorly understood at the molecular level. The glucosyltransferase, STARCH SYNTHASE 4 (SS4), plays a central role in granule initiation in Arabidopsis leaves, but its function in cereal endosperms is unknown. We investigated the role of SS4 in wheat, which has a distinct spatiotemporal pattern of granule initiation during grain development. We generated TILLING mutants in tetraploid wheat (Triticum turgidum) that are defective in both SS4 homoeologs. The morphology of endosperm starch was examined in developing and mature grains. SS4 deficiency led to severe alterations in endosperm starch granule morphology. During early grain development, while the wild-type initiated single 'A-type' granules per amyloplast, most amyloplasts in the mutant formed compound granules due to multiple initiations. This phenotype was similar to mutants deficient in B-GRANULE CONTENT 1 (BGC1). SS4 deficiency also reduced starch content in leaves and pollen grains. We propose that SS4 and BGC1 are required for the proper control of granule initiation during early grain development that leads to a single A-type granule per amyloplast. The absence of either protein results in a variable number of initiations per amyloplast and compound granule formation.

Domestication of Crop Metabolomes: Desired and Unintended Consequences

The majority of the crops and vegetables of today were domesticated from their wild progenitors within the past 12 000 years. Considerable research effort has been expended on characterizing the genes undergoing positive and negative selection during the processes of crop domestication and improvement. Many studies have also documented how the contents of a handful of metabolites have been altered during human selection, but we are only beginning to unravel the true extent of the metabolic consequences of breeding. We highlight how crop metabolomes have been wittingly or unwittingly shaped by the processes of domestication, and highlight how we can identify new targets for metabolite engineering for the purpose of de novo domestication of crop wild relatives.

The Landscape of the Genomic Distribution and the Expression of the F-Box Genes Unveil Genome Plasticity in Hexaploid Wheat during Grain Development and in Response to Heat and Drought Stress

FBX proteins are subunits of the SCF complex (Skp1-cullin-FBX) belonging to the E3 ligase family, which is involved in the ubiquitin-proteasome 26S (UPS) pathway responsible for the post-translational protein turnover. By targeting, in a selective manner, key regulatory proteins for ubiquitination and 26S proteasome degradation, FBX proteins play a major role in plant responses to diverse developmental and stress conditions. Although studies on the genomic organization of the FBX gene family in various species have been reported, knowledge related to bread wheat (Triticum aestivum) is scarce and needs to be broadened. Using the latest assembly of the wheat genome, we identified 3670 TaFBX genes distributed non-homogeneously within the three subgenomes (A, B and D) and between the 21 chromosomes, establishing it as one of the richest gene families among plant species. Based on the presence of the five different chromosomal regions previously identified, the present study focused on the genomic distribution of the TaFBX family and the identification of differentially expressed genes during the embryogenesis stages and in response to heat and drought stress. Most of the time, when comparing the expected number of genes (taking into account the formal gene distribution on the entire wheat genome), the TaFBX family harbors a different pattern at the various stratum of observation (subgenome, chromosome, chromosomal regions). We report here that the local gene expansion of the TaFBX family must be the consequence of multiple and complex events, including tandem and small-scale duplications. Regarding the differentially expressed TaFBX genes, while the majority of the genes are localized in the distal chromosomal regions (R1 and R3), differentially expressed genes are more present in the interstitial regions (R2a and R2b) than expected, which could be an indication of the preservation of major genes in those specific chromosomal regions.

Developmental and genomic architecture of plant embryogenesis: from model plant to crops

Embryonic development represents an important reproductive phase of sexually reproducing plant species. The fusion of egg and sperm produces the plant zygote, a totipotent cell that, through cell division and cell identity specification in early embryogenesis, establishes the major cell lineages and tissues of the adult plant. The subsequent morphogenesis phase produces the full-sized embryo, while the late embryogenesis maturation process prepares the seed for dormancy and subsequent germination, ensuring continuation of the plant life cycle. In this review on embryogenesis, we compare the model eudicot Arabidopsis thaliana with monocot crops, focusing on genome activation, paternal and maternal regulation of early zygote development, and key organizers of patterning, such as auxin and WOX transcription factors. While the early stages of embryo development are apparently conserved among plant species, embryo maturation programs have diversified between eudicots and monocots. This diversification in crop species reflects the likely effects of domestication on seed quality traits that are determined during embryo maturation, and also assures seed germination in different environmental conditions. This review describes the most important features of embryonic development in plants, and the scope and applications of genomics in plant embryo studies.

Tissue-Specific Transcriptome Analysis Reveals Candidate Transcripts Associated with the Process of Programmed B Chromosome Elimination in Aegilops speltoides

Some eukaryotes exhibit dramatic genome size differences between cells of different organs, resulting from programmed elimination of chromosomes. Here, we present the first transcriptome analysis of programmed chromosome elimination using laser capture microdissection (LCM)-based isolation of the central meristematic region of Aegilops speltoides embryos where B chromosome (B) elimination occurs. The comparative RNA-seq analysis of meristematic cells of embryos with (Bplus) and without Bs (B0) allowed the identification of 14,578 transcript isoforms (35% out of 41,615 analyzed transcript isoforms) that are differentially expressed during the elimination of Bs. A total of 2908 annotated unigenes were found to be up-regulated in Bplus condition. These genes are either associated with the process of B chromosome elimination or with the presence of B chromosomes themselves. GO enrichment analysis categorized up-regulated transcript isoforms into 27 overrepresented terms related to the biological process, nine terms of the molecular function aspect and three terms of the cellular component category. A total of 2726 annotated unigenes were down-regulated in Bplus condition. Based on strict filtering criteria, 341 B-unique transcript isoforms could be identified in central meristematic cells, of which 70 were functionally annotated. Beside others, genes associated with chromosome segregation, kinetochore function and spindle checkpoint activity were retrieved as promising candidates involved in the process of B chromosome elimination.

The NAC transcription factor NAC019-A1 is a negative regulator of starch synthesis in wheat developing endosperm

Starch is a major component of wheat (Triticum aestivum L.) endosperm and is an important part of the human diet. The functions of many starch synthesis genes have been elucidated. However, little is known about their regulatory mechanisms in wheat. Here, we identified a novel NAC transcription factor, TaNAC019-A1 (TraesCS3A02G077900), that negatively regulates starch synthesis in wheat and rice (Oryza sativa L.) endosperms. TaNAC019-A1 was highly expressed in the endosperm of developing grains and encoded a nucleus-localized transcriptional repressor. Overexpression of TaNAC019-A1 in rice and wheat led to significantly reduced starch content, kernel weight, and kernel width. The TaNAC019-A1-overexpression wheat lines had smaller A-type starch granules and fewer B-type starch granules than wild-type. Moreover, TaNAC019-A1 could directly bind to the 'ACGCAG' motif in the promoter regions of ADP-glucose pyrophosphorylase small subunit 1 (TaAGPS1-A1, TraesCS7A02G287400) and TaAGPS1-B1 (TraesCS7B02G183300) and repress their expression, thereby inhibiting starch synthesis in wheat endosperm. One haplotype of TaNAC019-B1 (TaNAC019-B1-Hap2, TraesCS3B02G092800) was positively associated with thousand-kernel weight and underwent positive selection during the Chinese wheat breeding process. Our data demonstrate that TaNAC019-A1 is a negative regulator of starch synthesis in wheat endosperm and provide novel insight into wheat yield improvement.

Genome-wide survey of the amino acid transporter gene family in wheat (Triticum aestivum L.): Identification, expression analysis and response to abiotic stress

Amino acid transporters (AATs), which transport amino acids across cell membranes, play important roles in alleviating plant damage under stresses and in plant growth. To data, little is known about the AAT genes in wheat because of its complex genome. In this study, a total of 296 AAT genes were identified from the latest wheat genome sequence (IWGSC v1.1) and classified into 12 distinct subfamilies based upon their sequence composition and phylogenetic relationship. The expansion of the wheat AAT family was mainly the results of whole-genome duplication (WGD) and tandem events. The unequal expansion of different subfamilies brought new features to TaAATs. TaAATs were highly expressed and exhibited distinct expression patterns in different tissues. On the basis of homology and expression pattern analysis, we identified several wheat AAT family members that may affect grain quality. In addition, TaAAP3, TaATLa2 and TaATLb13 exhibited sustained expression in response to drought and high-temperature stress. These genes are involved in the response of wheat to abiotic stress by regulating the transport and distribution of amino acids. Overall, our results help to understand the complexity of TaAATs and provide a theoretical basis for further identification and utilization of AATs in wheat and other crop species.

Multi-Locus GWAS of Quality Traits in Bread Wheat: Mining More Candidate Genes and Possible Regulatory Network

In wheat breeding, improved quality traits, including grain quality and dough rheological properties, have long been a critical goal. To understand the genetic basis of key quality traits of wheat, two single-locus and five multi-locus GWAS models were performed for six grain quality traits and three dough rheological properties based on 19, 254 SNPs in 267 bread wheat accessions. As a result, 299 quantitative trait nucleotides (QTNs) within 105 regions were identified to be associated with these quality traits in four environments. Of which, 40 core QTN regions were stably detected in at least three environments, 19 of which were novel. Compared with the previous studies, these novel QTN regions explained smaller phenotypic variation, which verified the advantages of the multi-locus GWAS models in detecting important small effect QTNs associated with complex traits. After characterization of the function and expression in-depth, 67 core candidate genes involved in protein/sugar synthesis, histone modification and the regulation of transcription factor were observed to be associated with the formation of grain quality, which showed that multi-level regulations influenced wheat grain quality. Finally, a preliminary network of gene regulation that may affect wheat quality formation was inferred. This study verified the power and reliability of multi-locus GWAS methods in wheat quality trait research, and increased the understanding of wheat quality formation mechanisms. The detected QTN regions and candidate genes in this study could be further used for gene cloning and marker-assisted selection in high-quality breeding of bread wheat.