Genome interplay in the grain transcriptome of hexaploid bread wheat

A gain-of-function mutation at the C-terminus of FT-D1 promotes heading by interacting with 14-3-3A and FDL6 in wheat

Vernalization and photoperiod pathways converging at FT1 control the transition to flowering in wheat. Here, we identified a gain-of-function mutation in FT-D1 that results in earlier heading date (HD), and shorter plant height and spike length in the gamma ray-induced eh1 wheat mutant. Knockout of the wild-type and overexpression of the mutated FT-D1 indicate that both alleles are functional to affect HD and plant height. Protein interaction assays demonstrated that the frameshift mutation in FT-D1eh1 exon 3 led to gain-of-function interactions with 14-3-3A and FDL6, thereby enabling the formation of florigen activation complex (FAC) and consequently activating a flowering-related transcriptomic programme. This mutation did not affect FT-D1eh1 interactions with TaNaKR5 or TaFTIP7, both of which could modulate HD, potentially via mediating FT-D1 translocation to the shoot apical meristem. Furthermore, the 'Segment B' external loop is essential for FT-D1 interaction with FDL6, while residue Y85 is required for interactions with TaNaKR5 and TaFTIP7. Finally, the flowering regulatory hub gene, ELF5, was identified as the FT-D1 regulatory target. This study illustrates FT-D1 function in determining wheat HD with a suite of interaction partners and provides genetic resources for tuning HD in elite wheat lines.

A Quantitative Computational Framework for Allopolyploid Single-Cell Data Integration and Core Gene Ranking in Development

Polyploidization drives regulatory and phenotypic innovation. How the merger of different genomes contributes to polyploid development is a fundamental issue in evolutionary developmental biology and breeding research. Clarifying this issue is challenging because of genome complexity and the difficulty in tracking stochastic subgenome divergence during development. Recent single-cell sequencing techniques enabled probing subgenome-divergent regulation in the context of cellular differentiation. However, analyzing single-cell data suffers from high error rates due to high dimensionality, noise, and sparsity, and the errors stack up in polyploid analysis due to the increased dimensionality of comparisons between subgenomes of each cell, hindering deeper mechanistic understandings. In this study, we develop a quantitative computational framework, called "pseudo-genome divergence quantification" (pgDQ), for quantifying and tracking subgenome divergence directly at the cellular level. Further comparing with cellular differentiation trajectories derived from single-cell RNA sequencing data allows for an examination of the relationship between subgenome divergence and the progression of development. pgDQ produces robust results and is insensitive to data dropout and noise, avoiding high error rates due to multiple comparisons of genes, cells, and subgenomes. A statistical diagnostic approach is proposed to identify genes that are central to subgenome divergence during development, which facilitates the integration of different data modalities, enabling the identification of factors and pathways that mediate subgenome-divergent activity during development. Case studies have demonstrated that applying pgDQ to single-cell and bulk tissue transcriptomic data promotes a systematic and deeper understanding of how dynamic subgenome divergence contributes to developmental trajectories in polyploid evolution.

Understanding grain development in the Poaceae family by comparing conserved and distinctive pathways through omics studies in wheat and maize

The Poaceae family, commonly known as the grass family, encompasses a diverse group of crops that play an essential role in providing food, fodder, biofuels, environmental conservation, and cultural value for both human and environmental well-being. Crops in Poaceae family are deeply intertwined with human societies, economies, and ecosystems, making it one of the most significant plant families in the world. As the major reservoirs of essential nutrients, seed grain of these crops has garnered substantial attention from researchers. Understanding the molecular and genetic processes that controls seed formation, development and maturation can provide insights for improving crop yield, nutritional quality, and stress tolerance. The diversity in photosynthetic pathways between C3 and C4 plants introduces intriguing variations in their physiological and biochemical processes, potentially affecting seed development. In this review, we explore recent studies performed with omics technologies, such as genomics, transcriptomics, proteomics and metabolomics that shed light on the mechanisms underlying seed development in wheat and maize, as representatives of C3 and C4 plants respectively, providing insights into their unique adaptations and strategies for reproductive success.

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.

RNA-Seq transcriptome profiling of immature grain wheat is a technique for understanding comparative modeling of baking quality

Improving the baking quality is a primary challenge in the wheat flour production value chain, as baking quality represents a crucial factor in determining its overall value. In the present study, we conducted a comparative RNA-Seq analysis on the high baking quality mutant "O-64.1.10" genotype and its low baking quality wild type "Omid" cultivar to recognize potential genes associated with bread quality. The cDNA libraries were constructed from immature grains that were 15 days post-anthesis, with an average of 16.24 and 18.97 million paired-end short-read sequences in the mutant and wild-type, respectively. A total number of 733 transcripts with differential expression were identified, 585 genes up-regulated and 188 genes down-regulated in the "O-64.1.10" genotype compared to the "Omid". In addition, the families of HSF, bZIP, C2C2-Dof, B3-ARF, BES1, C3H, GRF, HB-HD-ZIP, PLATZ, MADS-MIKC, GARP-G2-like, NAC, OFP and TUB were appeared as the key transcription factors with specific expression in the "O-64.1.10" genotype. At the same time, pathways related to baking quality were identified through Kyoto Encyclopedia of Genes and Genomes. Collectively, we found that the endoplasmic network, metabolic pathways, secondary metabolite biosynthesis, hormone signaling pathway, B group vitamins, protein pathways, pathways associated with carbohydrate and fat metabolism, as well as the biosynthesis and metabolism of various amino acids, have a great deal of potential to play a significant role in the baking quality. Ultimately, the RNA-seq results were confirmed using quantitative Reverse Transcription PCR for some hub genes such as alpha-gliadin, low molecular weight glutenin subunit and terpene synthase (gibberellin) and as a resource for future study, 127 EST-SSR primers were generated using RNA-seq data.

An Accurate Representation of the Number of bZIP Transcription Factors in the Triticum aestivum (Wheat) Genome and the Regulation of Functional Genes during Salt Stress

Climate change is dramatically increasing the overall area of saline soils around the world, which is increasing by approximately two million hectares each year. Soil salinity decreases crop yields and, thereby, makes farming less profitable, potentially causing increased poverty and hunger in many areas. A solution to this problem is increasing the salt tolerance of crop plants. Transcription factors (TFs) within crop plants represent a key to understanding salt tolerance, as these proteins play important roles in the regulation of functional genes linked to salt stress. The basic leucine zipper (bZIP) TF has a well-documented role in the regulation of salt tolerance. To better understand how bZIP TFs are linked to salt tolerance, we performed a genome-wide analysis in wheat using the Chinese spring wheat genome, which has been assembled by the International Wheat Genome Sequencing Consortium. We identified 89 additional bZIP gene sequences, which brings the total of bZIP gene sequences in wheat to 237. The majority of these 237 sequences included a single bZIP protein domain; however, different combinations of five other domains also exist. The bZIP proteins are divided into ten subfamily groups. Using an in silico analysis, we identified five bZIP genes (ABF2, ABF4, ABI5, EMBP1, and VIP1) that were involved in regulating salt stress. By scrutinizing the binding properties to the 2000 bp upstream region, we identified putative functional genes under the regulation of these TFs. Expression analyses of plant tissue that had been treated with or without 100 mM NaCl revealed variable patterns between the TFs and functional genes. For example, an increased expression of ABF4 was correlated with an increased expression of the corresponding functional genes in both root and shoot tissues, whereas VIP1 downregulation in root tissues strongly decreased the expression of two functional genes. Identifying strategies to sustain the expression of the functional genes described in this study could enhance wheat's salt tolerance.

Innovative computational tools provide new insights into the polyploid wheat genome

Bread wheat (Triticum aestivum) is an important crop and serves as a significant source of protein and calories for humans, worldwide. Nevertheless, its large and allopolyploid genome poses constraints on genetic improvement. The complex reticulate evolutionary history and the intricacy of genomic resources make the deciphering of the functional genome considerably more challenging. Recently, we have developed a comprehensive list of versatile computational tools with the integration of statistical models for dissecting the polyploid wheat genome. Here, we summarize the methodological innovations and applications of these tools and databases. A series of step-by-step examples illustrates how these tools can be utilized for dissecting wheat germplasm resources and unveiling functional genes associated with important agronomic traits. Furthermore, we outline future perspectives on new advanced tools and databases, taking into consideration the unique features of bread wheat, to accelerate genomic-assisted wheat breeding.

RNA-Seq-Based WGCNA and Association Analysis Reveal the Key Regulatory Module and Genes Responding to Salt Stress in Wheat Roots

Soil salinization is the main abiotic stressor faced by crops. An improved understanding of the transcriptional response to salt stress in roots, the organ directly exposed to a high salinity environment, can inform breeding strategies to enhance tolerance and increase crop yield. Here, RNA-sequencing was performed on the roots of salt-tolerant wheat breeding line CH7034 at 0, 1, 6, 24, and 48 h after NaCl treatment. Based on transcriptome data, a weighted gene co-expression network analysis (WGCNA) was constructed, and five gene co-expression modules were obtained, of which the blue module was correlated with the time course of salt stress at 1 and 48 h. Two GO terms containing 249 differentially expressed genes (DEGs) related to osmotic stress response and salt-stress response were enriched in the blue module. These DEGs were subsequently used for association analysis with a set of wheat germplasm resources, and the results showed that four genes, namely a Walls Are Thin 1-related gene (TaWAT), an aquaporin gene (TaAQP), a glutathione S-transfer gene (TaGST), and a zinc finger gene (TaZFP), were associated with the root salt-tolerance phenotype. Using the four candidate genes as hub genes, a co-expression network was constructed with another 20 DEGs with edge weights greater than 0.6. The network showed that TaWAT and TaAQP were mainly co-expressed with fifteen interacting DEGs 1 h after salt treatment, while TaGST and TaZFP were mainly co-expressed with five interacting DEGs 48 h after salt treatment. This study provides key modules and candidate genes for understanding the salt-stress response mechanism in wheat roots.

Identification of candidate regulators of the response to early heat stress in climate-adapted wheat landraces via transcriptomic and co-expression network analyses

Introduction

Climate change is likely to lead to not only increased global temperatures but also a more variable climate where unseasonal periods of heat stress are more prevalent. This has been evidenced by the observation of spring-time temperatures approaching 40°C in some of the main spring-wheat producing countries, such as the USA, in recent years. With an optimum growth temperature of around 20°C, wheat is particularly prone to damage by heat stress. A warming climate with increasingly common fluctuations in temperature therefore threatens wheat crops and subsequently the lives and livelihoods of billions of people who depend on the crop for food. To futureproof wheat against a variable climate, a better understanding of the response to early heat stress is required.

Methods

Here, we utilised DESeq2 to identify 7,827 genes which were differentially expressed in wheat landraces after early heat stress exposure. Candidate hub genes, which may regulate the transcriptional response to early heat stress, were identified via weighted gene co-expression network analysis (WGCNA), and validated by qRT-PCR.

Results

Two of the most promising candidate hub genes (TraesCS3B02G409300 and TraesCS1B02G384900) may downregulate the expression of genes involved in the drought, salinity, and cold responses-genes which are unlikely to be required under heat stress-as well as photosynthesis genes and stress hormone signalling repressors, respectively. We also suggest a role for a poorly characterised sHSP hub gene (TraesCS4D02G212300), as an activator of the heat stress response, potentially inducing the expression of a vast suite of heat shock proteins and transcription factors known to play key roles in the heat stress response.

Discussion

The present work represents an exploratory examination of the heat-induced transcriptional change in wheat landrace seedlings and identifies several candidate hub genes which may act as regulators of this response and, thus, may be targets for breeders in the production of thermotolerant wheat varieties.

Type I MADS-box transcription factor TaMADS-GS regulates grain size by stabilizing cytokinin signalling during endosperm cellularization in wheat

Grain size is one of the important traits in wheat breeding programs aimed at improving yield, and cytokinins, mainly involved in cell division, have a positive impact on grain size. Here, we identified a novel wheat gene TaMADS-GS encoding type I MADS-box transcription factor, which regulates the cytokinins signalling pathway during early stages of grain development to modulate grain size and weight in wheat. TaMADS-GS is exclusively expressed in grains at early stage of seed development and its knockout leads to delayed endosperm cellularization, smaller grain size and lower grain weight. TaMADS-GS protein interacts with the Polycomb Repressive Complex 2 (PRC2) and leads to repression of genes encoding cytokinin oxidase/dehydrogenases (CKXs) stimulating cytokinins inactivation by mediating accumulation of the histone H3 trimethylation at lysine 27 (H3K27me3). Through the screening of a large wheat germplasm collection, an elite allele of the TaMADS-GS exhibits higher ability to repress expression of genes inactivating cytokinins and a positive correlation with grain size and weight, thus representing a novel marker for breeding programs in wheat. Overall, these findings support the relevance of TaMADS-GS as a key regulator of wheat grain size and weight.

A papain-like cysteine protease-released small signal peptide confers wheat resistance to wheat yellow mosaic virus

Wheat yellow mosaic virus (WYMV), a soil-borne pathogen, poses a serious threat to global wheat production. Here, we identify a WYMV resistance gene, TaRD21A, that belongs to the papain-like cysteine protease family. Through genetic manipulation of TaRD21A expression, we establish its positive role in the regulation of wheat to WYMV resistance. Furthermore, our investigation shows that the TaRD21A-mediated plant antiviral response relies on the release of a small peptide catalyzed by TaRD21A protease activity. To counteract wheat resistance, WYMV-encoded nuclear inclusion protease-a (NIa) suppress TaRD21A activity to promote virus infection. In resistant cultivars, a natural variant of TaRD21A features a glycine-to-threonine substitution and this substitution enables the phosphorylation of threonine, thereby weakening the interaction between NIa and TaRD21A, reinforcing wheat resistance against WYMV. Our study not only unveils a WYMV resistance gene but also offers insights into the intricate mechanisms underpinning resistance against WYMV.

LHP1-mediated epigenetic buffering of subgenome diversity and defense responses confers genome plasticity and adaptability in allopolyploid wheat

Polyploidization is a major driver of genome diversification and environmental adaptation. However, the merger of different genomes may result in genomic conflicts, raising a major question regarding how genetic diversity is interpreted and regulated to enable environmental plasticity. By analyzing the genome-wide binding of 191 trans-factors in allopolyploid wheat, we identified like heterochromatin protein 1 (LHP1) as a master regulator of subgenome-diversified genes. Transcriptomic and epigenomic analyses of LHP1 mutants reveal its role in buffering the expression of subgenome-diversified defense genes by controlling H3K27me3 homeostasis. Stripe rust infection releases latent subgenomic variations by eliminating H3K27me3-related repression. The simultaneous inactivation of LHP1 homoeologs by CRISPR-Cas9 confers robust stripe rust resistance in wheat seedlings. The conditional repression of subgenome-diversified defenses ensures developmental plasticity to external changes, while also promoting neutral-to-non-neutral selection transitions and adaptive evolution. These findings establish an LHP1-mediated buffering system at the intersection of genotypes, environments, and phenotypes in polyploid wheat. Manipulating the epigenetic buffering capacity offers a tool to harness cryptic subgenomic variations for crop improvement.

Genome-Wide Association Study of Asian and European Common Wheat Accessions for Yield-Related Traits and Stripe Rust Resistance

Identifying novel loci of yield-related traits and resistance to stripe rust (caused by Puccinia striiformis f. sp. tritici) in wheat will help in breeding wheat that can meet projected demands in diverse environmental and agricultural practices. We performed a genome-wide association study with 24,767 single nucleotide polymorphisms (SNPs) in 180 wheat accessions that originated in 16 Asian or European countries between latitudes 30°N and 45°N. We detected seven accessions with desirable yield-related traits and 42 accessions that showed stable, high degrees of stripe rust resistance in multienvironment field assessments. A marker-trait association analysis of yield-related traits detected 18 quantitative trait loci (QTLs) in at least two test environments and two QTLs related to stripe rust resistance in at least three test environments. Five of these QTLs were identified as potentially novel QTLs by comparing their physical locations with those of known QTLs in the Chinese Spring (CS) reference genome RefSeq v1.1 published by the International Wheat Genome Sequencing Consortium; two were for spike length, one was for grain number per spike, one was for spike number, and one was for stripe rust resistance at the adult plant stage. We also identified 14 candidate genes associated with the five novel QTLs. These QTLs and candidate genes will provide breeders with new germplasm and can be used to conduct marker-assisted selection in breeding wheat with improved yield and stripe rust resistance.

Co-expression network analysis of diverse wheat landraces reveals markers of early thermotolerance and a candidate master regulator of thermotolerance genes

Triticum aestivum L. (bread wheat) is a crop relied upon by billions of people around the world, as a major source of both income and calories. Rising global temperatures, however, pose a genuine threat to the livelihood of these people, as wheat growth and yields are extremely vulnerable to damage by heat stress. Here we present the YoGI wheat landrace panel, comprising 342 accessions that show remarkable phenotypic and genetic diversity thanks to their adaptation to different climates. We quantified the abundance of 110 790 transcripts from the panel and used these data to conduct weighted co-expression network analysis and to identify hub genes in modules associated with abiotic stress tolerance. We found that the expression of three hub genes, all heat-shock proteins (HSPs), were significantly correlated with early thermotolerance in a validation panel of landraces. These hub genes belong to the same module, with one (TraesCS4D01G207500.1) being a candidate master-regulator potentially controlling the expression of the other two hub genes, as well as a suite of other HSPs and heat-stress transcription factors (HSFs). In this work, therefore, we identify three validated hub genes, the expression of which can serve as markers of thermotolerance during early development, and suggest that TraesCS4D01G207500.1 is a potential master regulator of HSP and HSF expression - presenting the YoGI landrace panel as an invaluable tool for breeders wishing to determine and introduce novel alleles into modern varieties, for the production of climate-resilient crops.

Design, execution, and interpretation of plant RNA-seq analyses

Genomics has transformed our understanding of the genetic architecture of traits and the genetic variation present in plants. Here, we present a review of how RNA-seq can be performed to tackle research challenges addressed by plant sciences. We discuss the importance of experimental design in RNA-seq, including considerations for sampling and replication, to avoid pitfalls and wasted resources. Approaches for processing RNA-seq data include quality control and counting features, and we describe common approaches and variations. Though differential gene expression analysis is the most common analysis of RNA-seq data, we review multiple methods for assessing gene expression, including detecting allele-specific gene expression and building co-expression networks. With the production of more RNA-seq data, strategies for integrating these data into genetic mapping pipelines is of increased interest. Finally, special considerations for RNA-seq analysis and interpretation in plants are needed, due to the high genome complexity common across plants. By incorporating informed decisions throughout an RNA-seq experiment, we can increase the knowledge gained.

Transcriptomic and co-expression network analyses on diverse wheat landraces identifies candidate master regulators of the response to early drought

Introduction

Over four billion people around the world rely on bread wheat (Triticum aestivum L.) as a major constituent of their diet. The changing climate, however, threatens the food security of these people, with periods of intense drought stress already causing widespread wheat yield losses. Much of the research into the wheat drought response has centred on the response to drought events later in development, during anthesis or grain filling. But as the timing of periods of drought stress become increasingly unpredictable, a more complete understanding of the response to drought during early development is also needed.

Methods

Here, we utilized the YoGI landrace panel to identify 10,199 genes which were differentially expressed under early drought stress, before weighted gene co-expression network analysis (WGCNA) was used to construct a co-expression network and identify hub genes in modules particularly associated with the early drought response.

Results

Of these hub genes, two stood out as novel candidate master regulators of the early drought response - one as an activator (TaDHN4-D1; TraesCS5D02G379200) and the other as a repressor (uncharacterised gene; TraesCS3D02G361500).

Discussion

As well as appearing to coordinate the transcriptional early drought response, we propose that these hub genes may be able to regulate the physiological early drought response due to potential control over the expression of members of gene families well-known for their involvement in the drought response in many plant species, namely dehydrins and aquaporins, as well as other genes seemingly involved in key processes such as, stomatal opening, stomatal closing, stomatal morphogenesis and stress hormone signalling.

TaTPP-7A positively feedback regulates grain filling and wheat grain yield through T6P-SnRK1 signalling pathway and sugar-ABA interaction

Grain size and filling are two key determinants of grain thousand-kernel weight (TKW) and crop yield, therefore they have undergone strong selection since cereal was domesticated. Genetic dissection of the two traits will improve yield potential in crops. A quantitative trait locus significantly associated with wheat grain TKW was detected on chromosome 7AS flanked by a simple sequence repeat marker of Wmc17 in Chinese wheat 262 mini-core collection by genome-wide association study. Combined with the bulked segregant RNA-sequencing (BSR-seq) analysis of an F₂ genetic segregation population with extremely different TKW traits, a candidate trehalose-6-phosphate phosphatase gene located at 135.0 Mb (CS V1.0), designated as TaTPP-7A, was identified. This gene was specifically expressed in developing grains and strongly influenced grain filling and size. Overexpression (OE) of TaTPP-7A in wheat enhanced grain TKW and wheat yield greatly. Detailed analysis revealed that OE of TaTPP-7A significantly increased the expression levels of starch synthesis- and senescence-related genes involved in abscisic acid (ABA) and ethylene pathways. Moreover, most of the sucrose metabolism and starch regulation-related genes were potentially regulated by SnRK1. In addition, TaTPP-7A is a crucial domestication- and breeding-targeted gene and it feedback regulates sucrose lysis, flux, and utilization in the grain endosperm mainly through the T6P-SnRK1 pathway and sugar-ABA interaction. Thus, we confirmed the T6P signalling pathway as the central regulatory system for sucrose allocation and source-sink interactions in wheat grains and propose that the trehalose pathway components have great potential to increase yields in cereal crops.

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.

Genome-Wide Analysis of MYB Transcription Factors in the Wheat Genome and Their Roles in Salt Stress Response

Large and rapidly increasing areas of salt-affected soils are posing major challenges for the agricultural sector. Most fields used for the important food crop Triticum aestivum (wheat) are expected to be salt-affected within 50 years. To counter the associated problems, it is essential to understand the molecular mechanisms involved in salt stress responses and tolerance, thereby enabling their exploitation in the development of salt-tolerant varieties. The myeloblastosis (MYB) family of transcription factors are key regulators of responses to both biotic and abiotic stress, including salt stress. Thus, we used the Chinese spring wheat genome assembled by the International Wheat Genome Sequencing Consortium to identify putative MYB proteins (719 in total). Protein families (PFAM) analysis of the MYB sequences identified 28 combinations of 16 domains in the encoded proteins. The most common consisted of MYB_DNA-binding and MYB-DNA-bind_6 domains, and five highly conserved tryptophans were located in the aligned MYB protein sequence. Interestingly, we found and characterized a novel 5R-MYB group in the wheat genome. In silico studies showed that MYB transcription factors MYB3, MYB4, MYB13 and MYB59 are involved in salt stress responses. qPCR analysis confirmed upregulation of the expression of all these MYBs in both roots and shoots of the wheat variety BARI Gom-25 (except MYB4, which was downregulated in roots) under salt stress. Moreover, we identified nine target genes involved in salt stress that are regulated by the four MYB proteins, most of which have cellular locations and are involved in catalytic and binding activities associated with various cellular and metabolic processes.

Non-Rabl chromosome organization in endoreduplicated nuclei of barley embryo and endosperm tissues

Rabl organization is a type of interphase chromosome arrangement with centromeres and telomeres clustering at opposite nuclear poles. Here, we analyzed nuclear morphology and chromosome organization in cycling and endoreduplicated nuclei isolated from embryo and endosperm tissues of developing barley seeds. We show that endoreduplicated nuclei have an irregular shape, less sister chromatid cohesion at 5S rDNA loci, and a reduced amount of centromeric histone CENH3. While the chromosomes of the embryo and endosperm nuclei are initially organized in Rabl configuration, the centromeres and telomeres are intermingled within the nuclear space in the endoreduplicated nuclei with an increasing endoreduplication level. Such a loss of chromosome organization suggests that Rabl configuration is introduced and further reinforced by mitotic divisions in barley cell nuclei in a tissue- and seed age-dependent manner.

Global transcriptome analysis of allopolyploidization reveals large-scale repression of the D-subgenome in synthetic hexaploid wheat

Synthetic hexaploid wheat (SHW) lines are created as pre-breeding germplasm to diversify the D subgenome of hexaploid wheat and capitalize upon the untapped genetic diversity of the Aegilops tauschii gene pool. However, the phenotypes observed in the Ae. tauschii parents are not always recovered in the SHW lines, possibly due to inter-subgenome interactions. To elucidate this post-polyploidization genome reprogramming phenomenon, we performed RNA-seq of four SHW lines and their corresponding tetraploid and diploid parents, across ten tissues and three biological replicates. Homoeologue expression bias (HEB) analysis using more than 18,000 triads suggests massive suppression of homoeoalleles of the D subgenome in SHWs. Comparative transcriptome analysis of the whole-genome gene set further corroborated this finding. Alternative splicing analysis of the high-confidence genes indicates an additional layer of complexity where all five splice events are identified, and retained intron is predominant. Homoeologue expression upon resynthesis of hexaploid wheat has implications to the usage and handling of this germplasm in breeding as it relates to capturing the effects of epistatic interaction across subgenomes upon polyploidization. Special considerations must be given to this germplasm in pre-breeding activities to consider the extent of the inter-subgenome interactions on gene expression and their impact on traits for crop improvement.

Dissection of Developmental Programs and Regulatory Modules Directing Endosperm Transfer Cell and Aleurone Identity in the Syncytial Endosperm of Barley

Endosperm development in barley starts with the formation of a multinucleate syncytium, followed by cellularization in the ventral part of the syncytium generating endosperm transfer cells (ETCs) as first differentiating subdomain, whereas aleurone (AL) cells will originate from the periphery of the enclosing syncytium. Positional signaling in the syncytial stage determines cell identity in the cereal endosperm. Here, we performed a morphological analysis and employed laser capture microdissection (LCM)-based RNA-seq of the ETC region and the peripheral syncytium at the onset of cellularization to dissect developmental and regulatory programs directing cell specification in the early endosperm. Transcriptome data revealed domain-specific characteristics and identified two-component signaling (TCS) and hormone activities (auxin, ABA, ethylene) with associated transcription factors (TFs) as the main regulatory links for ETC specification. On the contrary, differential hormone signaling (canonical auxin, gibberellins, cytokinin) and interacting TFs control the duration of the syncytial phase and timing of cellularization of AL initials. Domain-specific expression of candidate genes was validated by in situ hybridization and putative protein-protein interactions were confirmed by split-YFP assays. This is the first transcriptome analysis dissecting syncytial subdomains of cereal seeds and provides an essential framework for initial endosperm differentiation in barley, which is likely also valuable for comparative studies with other cereal crops.

Form follows function in Triticeae inflorescences

Grass inflorescences produce grains, which are directly connected to our food. In grass crops, yields are mainly affected by grain number and weight; thus, understanding inflorescence shape is crucially important for cereal crop breeding. In the last two decades, several key genes controlling inflorescence shape have been elucidated, thanks to the availability of rich genetic resources and powerful genomics tools. In this review, we focus on the inflorescence architecture of Triticeae species, including the major cereal crops wheat and barley. We summarize recent advances in our understanding of the genetic basis of spike branching, and spikelet and floret development in the Triticeae. Considering our changing climate and its impacts on cereal crop yields, we also discuss the future orientation of research.

Insights into the Bioinformatics and Transcriptional Analysis of the Elongator Complexes (ELPs) Gene Family of Wheat: TaELPs Contribute to Wheat Abiotic Stress Tolerance and Leaf Senescence

Elongator complexes (ELPs) are the protein complexes that promote transcription through histone acetylation in eukaryotic cells and interact with elongating RNA polymerase II (RNAPII). ELPs' role in plant growth and development, signal transduction, and response to biotic and abiotic stresses have been confirmed in model plants. However, the functions of the wheat ELP genes are not well documented. The present study identified 18 members of the ELPs from the wheat genome with a homology search. Further, bioinformatics and transcription patterns in response to different stress conditions were analyzed to dissect their potential regulatory mechanisms in wheat. Gene duplication analysis showed that 18 pairs of ELP paralogous genes were derived from segmental duplication, which was divided into six clades by protein phylogenetic and cluster analysis. The orthologous analysis of wheat TaELP genes showed that TaELP genes may have evolved from orthologous genes of other plant species or closely related plants. Moreover, a variety of cis-acting regulatory elements (CAREs) related to growth and development, hormone response, and biotic and abiotic stresses were identified in the TaELPs' promoter regions. The qRT-PCR analysis showed that the transcription of TaELPs was induced under hormone, salt, and drought stress and during leaf senescence. The TaELP2 gene was silenced with BSMV-VIGS, and TaELP2 was preliminarily verified to be involved in the regulation of wheat leaf senescence. Overall, TaELP genes might be regulated by hormone signaling pathways and response to abiotic stress and leaf senescence, which could be investigated further as potential candidate genes for wheat abiotic stress tolerance and yield improvement.

QTL cluster analysis and marker development for kernel traits based on DArT markers in spring bread wheat (Triticum aestivum L.)

Genetic dissection of yield component traits including kernel characteristics is essential for the continuous improvement in wheat yield. In the present study, one recombinant inbred line (RIL) F₆ population derived from a cross between Avocet and Chilero was used to evaluate the phenotypes of kernel traits of thousand-kernel weight (TKW), kernel length (KL), and kernel width (KW) in four environments at three experimental stations during the 2018-2020 wheat growing seasons. The high-density genetic linkage map was constructed with the diversity arrays technology (DArT) markers and the inclusive composite interval mapping (ICIM) method to identify the quantitative trait loci (QTLs) for TKW, KL, and KW. A total of 48 QTLs for three traits were identified in the RIL population on the 21 chromosomes besides 2A, 4D, and 5B, accounting for 3.00%-33.85% of the phenotypic variances. Based on the physical positions of each QTL, nine stable QTL clusters were identified in the RILs, and among these QTL clusters, TaTKW-1A was tightly linked to the DArT marker interval 3950546-1213099, explaining 10.31%-33.85% of the phenotypic variances. A total of 347 high-confidence genes were identified in a 34.74-Mb physical interval. TraesCS1A02G045300 and TraesCS1A02G058400 were among the putative candidate genes associated with kernel traits, and they were expressed during grain development. Moreover, we also developed high-throughput kompetitive allele-specific PCR (KASP) markers of TaTKW-1A, validated in a natural population of 114 wheat varieties. The study provides a basis for cloning the functional genes underlying the QTL for kernel traits and a practical and accurate marker for molecular breeding.

Molecular prospective on the wheat grain development

Wheat grain development is an important biological process to determine grain yield and quality, which is controlled by the interplay of genetic, epigenetic, and environmental factors. Wheat grain development has been extensively characterized at the phenotypic and genetic levels. The advent of innovative molecular technologies allows us to characterize genes, proteins, and regulatory factors involved in wheat grain development, which have enhanced our understanding of the wheat seed development process. However, wheat is an allohexaploid with a large genome size, the molecular mechanisms underlying the wheat grain development have not been well understood as those in diploids. Understanding grain development, and how it is regulated, is of fundamental importance for improving grain yield and quality through conventional breeding or genetic engineering. Herein, we review the current discoveries on the molecular mechanisms underlying wheat grain development. Notably, only a handful of genes that control wheat grain development have, thus far, been well characterized, their interplay underlying the grain development remains elusive. The synergistic network-integrated genomics and epigenetics underlying wheat grain development and how the subgenome divergence dynamically and precisely regulates wheat grain development are unknown.

Low-affinity SPL binding sites contribute to subgenome expression divergence in allohexaploid wheat

Expression divergence caused by genetic variation and crosstalks among subgenomes of the allohexaploid bread wheat (Triticum aestivum. L., BBAADD) is hypothesized to increase its adaptability and/or plasticity. However, the molecular basis of expression divergence remains unclear. Squamosa promoter-binding protein-like (SPL) transcription factors are critical for a wide array of biological processes. In this study, we constructed expression regulatory networks by combining DAP-seq for 40 SPLs, ATAC-seq, and RNA-seq. Our findings indicate that a group of low-affinity SPL binding regions (SBRs) were targeted by diverse SPLs and caused different sequence preferences around the core GTAC motif. The SBRs including the low-affinity ones are evolutionarily conserved, enriched GWAS signals related to important agricultural traits. However, those SBRs are highly diversified among the cis-regulatory regions (CREs) of syntenic genes, with less than 8% SBRs coexisting in triad genes, suggesting that CRE variations are critical for subgenome differentiations. Knocking out of TaSPL7A/B/D and TaSPL15A/B/D subfamily further proved that both high- and low-affinity SBRs played critical roles in the differential expression of genes regulating tiller number and spike sizes. Our results have provided baseline data for downstream networks of SPLs and wheat improvements and revealed that CRE variations are critical sources for subgenome divergence in the allohexaploid wheat.

Transposable elements orchestrate subgenome-convergent and -divergent transcription in common wheat

The success of common wheat as a global staple crop was largely attributed to its genomic diversity and redundancy due to the merge of different genomes, giving rise to the major question how subgenome-divergent and -convergent transcription is mediated and harmonized in a single cell. Here, we create a catalog of genome-wide transcription factor-binding sites (TFBSs) to assemble a common wheat regulatory network on an unprecedented scale. A significant proportion of subgenome-divergent TFBSs are derived from differential expansions of particular transposable elements (TEs) in diploid progenitors, which contribute to subgenome-divergent transcription. Whereas subgenome-convergent transcription is associated with balanced TF binding at loci derived from TE expansions before diploid divergence. These TFBSs have retained in parallel during evolution of each diploid, despite extensive unbalanced turnover of the flanking TEs. Thus, the differential evolutionary selection of paleo- and neo-TEs contribute to subgenome-convergent and -divergent regulation in common wheat, highlighting the influence of TE repertory plasticity on transcriptional plasticity in polyploid.

Identification of Glu-D1 Alleles and Novel Marker-Trait Associations for Flour Quality and Grain Yield Traits under Heat-Stress Environments in Wheat Lines Derived from Diverse Accessions of Aegilops tauschii

Heat stress during grain filling is considered one of the major abiotic factors influencing wheat grain yield and quality in arid and semi-arid regions. We studied the effect of heat stress on flour quality and grain yield at moderate and continuous heat stress under natural field conditions using 147 lines of wheat multiple synthetic derivatives (MSD) containing Aegilops tauschii introgressions. The study aimed to identify the marker-trait associations (MTAs) for the quality traits and grain yield under heat-stress conditions and identify stress-resilient germplasm-combining traits for good flour quality and grain yield. The MSD lines showed considerable genetic variation for quality traits and grain yield under heat-stress conditions; some lines performed better than the recurrent parent, Norin 61. We identified two MSD lines that consistently maintained relative performance (RP) values above 100% for grain yield and dough strength. We found the presence of three high-molecular-weight glutenin subunits (HMW-GSs) at the Glu-D1 locus derived from Ae. tauschii, which were associated with stable dough strength across the four environments used in this study. These HMW-GSs could be potentially useful in applications for future improvements of end-use quality traits targeting wheat under severe heat stress. A total of 19,155 high-quality SNP markers were used for the genome-wide association analysis and 251 MTAs were identified, most of them on the D genome, confirming the power of the MSD panel as a platform for mining and exploring the genes of Ae. tauschii. We identified the MTAs for dough strength under heat stress, which simultaneously control grain yield and relative performance for dough strength under heat-stress/optimum conditions. This study proved that Ae. tauschii is an inexhaustible resource for genetic mining, and the identified lines and pleiotropic MTAs reported in this study are considered a good resource for the development of resilient wheat cultivars that combine both good flour quality and grain yield under stress conditions using marker-assisted selection.

Brachypodium: 20 years as a grass biology model system; the way forward?

It has been 20 years since Brachypodium distachyon was suggested as a model grass species, but ongoing research now encompasses the entire genus. Extensive Brachypodium genome sequencing programmes have provided resources to explore the determinants and drivers of population diversity. This has been accompanied by cytomolecular studies to make Brachypodium a platform to investigate speciation, polyploidisation, perenniality, and various aspects of chromosome and interphase nucleus organisation. The value of Brachypodium as a functional genomic platform has been underscored by the identification of key genes for development, biotic and abiotic stress, and cell wall structure and function. While Brachypodium is relevant to the biofuel industry, its impact goes far beyond that as an intriguing model to study climate change and combinatorial stress.

Consensus genomic regions associated with grain protein content in hexaploid and tetraploid wheat

A meta-analysis of QTLs associated with grain protein content (GPC) was conducted in hexaploid and tetraploid wheat to identify robust and stable meta-QTLs (MQTLs). For this purpose, as many as 459 GPC-related QTLs retrieved from 48 linkage-based QTL mapping studies were projected onto the newly developed wheat consensus map. The analysis resulted in the prediction of 57 MQTLs and 7 QTL hotspots located on all wheat chromosomes (except chromosomes 1D and 4D) and the average confidence interval reduced 2.71-fold in the MQTLs and QTL hotspots compared to the initial QTLs. The physical regions occupied by the MQTLs ranged from 140 bp to 224.02 Mb with an average of 15.2 Mb, whereas the physical regions occupied by QTL hotspots ranged from 1.81 Mb to 36.03 Mb with a mean of 8.82 Mb. Nineteen MQTLs and two QTL hotspots were also found to be co-localized with 45 significant SNPs identified in 16 previously published genome-wide association studies in wheat. Candidate gene (CG) investigation within some selected MQTLs led to the identification of 705 gene models which also included 96 high-confidence CGs showing significant expressions in different grain-related tissues and having probable roles in GPC regulation. These significantly expressed CGs mainly involved the genes/gene families encoding for the following proteins: aminotransferases, early nodulin 93, glutamine synthetases, invertase/pectin methylesterase inhibitors, protein BIG GRAIN 1-like, cytochrome P450, glycosyl transferases, hexokinases, small GTPases, UDP-glucuronosyl/UDP-glucosyltransferases, and EamA, SANT/Myb, GNAT, thioredoxin, phytocyanin, and homeobox domains containing proteins. Further, eight genes including GPC-B1, Glu-B1-1b, Glu-1By9, TaBiP1, GSr, TaNAC019-A, TaNAC019-D, and bZIP-TF SPA already known to be associated with GPC were also detected within some of the MQTL regions confirming the efficacy of MQTLs predicted during the current study.

Transcriptional regulation of photomorphogenesis in seedlings of Brassica napus under different light qualities

This study displayed the transcriptional regulation network of key regulators and downstream pathway in seedling morphogenesis of Brassica napus under different light quality. Plants undergo photomorphogenesis upon the presence of light, mediated by different light (e.g., blue, red, and far-red) signaling pathways. Although the light signaling pathway has been well documented in Arabidopsis, the underlying mechanisms were studied to a less extent in other plant species including Brassica napus. In this study, we investigated the effect of different light qualities (white, blue, red, and far-red light) on the hypocotyl elongation in B. napus, and performed the transcriptomic analysis of seedlings in response to different light qualities. The results showed that hypocotyl elongation was slightly inhibited by red light, while it was strongly inhibited by blue/far-red light. Transcriptome analysis identified 9748 differentially expressed genes (DEGs) among treatments. Gene ontology (GO) enrichment analysis of DEGs showed that light-responsive and photosynthesis-related genes were highly expressed in response to blue/far-red light rather than in red light. Furthermore, the key genes in light signaling (i.e., PHYB, HY5, HYH, HFR1, and PIF3) exhibited distinct expression patterns between blue/far-red and red light treatments. In addition, subgenome dominant expression of homoeologous genes were observed for some genes, such as PHYA, PHYB, HFR1, and BBXs. The current study displayed a comprehensive dissection of light-mediated transcriptional regulation network, including light signaling, phytohormone, and cell elongation/modification, which improved the understanding on the underlying mechanism of light-regulated hypocotyl growth in B. napus.

New insight into the sucrose biosynthesis pathway from genome-wide identification, gene expression analysis, and subcellular localization in hexaploid wheat (Triticum aestivum L.)

Sucrose, the main synthesized product and transported form of photoassimilates, moves from leaves to support plant growth and storage substance biosynthesis occurring in the heterotrophic sink organs. Enhancing sucrose biosynthesis efficiency is a top priority for crop yield breeding programs. However, the molecular mechanism of sucrose biosynthesis is still elusive, especially in wheat. We performed transcriptome sequencing, subcellular localization, and bioinformatics analysis to identify key isoforms and metabolic branches associated with sucrose biosynthesis in wheat. Our analysis identified 45 genes from 13 families that exhibited high expression in young leaves with an evident diurnal change. The carbon flux from photoassimilates to sucrose was divided into two pathways. In the cytoplasm, assimilates initiating at phosphotrioses (TPs) exported by TaTPT1 from chloroplasts flowed along the TP-Sucrose branch formed by TaALD6, TaFBP5, TacPGI, TacPGM, TaUGP1, TaSPS5, and TaSPP1. Intermediates either from the Calvin cycle or TP-Sucrose branch were converted into ADPGlc to synthesize the simple starch, which was transiently degraded by a series of enzymes, including TaBAM4 and TaSEX4 in the chloroplast. Similar to the starch-biosynthesis branch in endosperms, the TP-Sucrose branch is the most prominent in leaves because each reaction can be catalyzed by at least one highly expressed isoform with expected cytosolic localization. The key isoforms and major branches highlighted in the wheat sucrose biosynthesis pathway expand our molecular understanding of crop sucrose biosynthesis and provide clues to increase wheat yield by enhancing the sucrose synthesis efficiency of leaves.

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.

Asymmetric subgenomic chromatin architecture impacts on gene expression in resynthesized and natural allopolyploid Brassica napus

Although asymmetric subgenomic epigenetic modification and gene expression have been revealed in the successful establishment of allopolyploids, the changes in chromatin accessibility and their relationship with epigenetic modifications and gene expression are poorly understood. Here, we synthetically analyzed chromatin accessibility, four epigenetic modifications and gene expression in natural allopolyploid Brassica napus, resynthesized allopolyploid B. napus, and diploid progenitors (B. rapa and B. oleracea). “Chromatin accessibility shock” occurred in both allopolyploidization and natural evolutionary processes, and genic accessible chromatin regions (ACRs) increased after allopolyploidization. ACRs associated with H3K27me3 modifications were more accessible than those with H3K27ac or H3K4me3. Although overall chromatin accessibility may be defined by H3K27me3, the enrichment of H3K4me3 and H3K27ac and depletion of DNA methylation around transcriptional start sites up-regulated gene expression. Moreover, we found that subgenome C_n exhibited higher chromatin accessibility than A_n, which depended on the higher chromatin accessibility of C_n-unique genes but not homologous genes.

Changes in chromatin accessibility occuring during the process of allopolyploidization of Brassica napus are analysed using ATAC and ChIPseq, with differences in asymmetric chromatin accessibility between subgenomes of B. napus investigated.

The molecular basis of cereal grain proteostasis

Storage proteins deposited in the endosperm of cereal grains are both a nitrogen reserve for seed germination and seedling growth and a primary protein source for human nutrition. Detailed surveys of the patterns of storage protein accumulation in cereal grains during grain development have been undertaken, but an in-depth understanding of the molecular mechanisms that regulate these patterns is still lacking. Accumulation of storage proteins in cereal grains involves a series of subcellular compartments, a set of energy-dependent events that compete with other cellular processes, and a balance of protein synthesis and protein degradation rates at different times during the developmental process. In this review, we focus on the importance of rates in cereal grain storage protein accumulation during grain development and outline the potential implications and applications of this information to accelerate modern agriculture breeding programmes and optimize energy use efficiency in proteostasis.

Unravelling consensus genomic regions associated with quality traits in wheat using meta-analysis of quantitative trait loci

Meta-analysis in wheat for three major quality traits identified 110 meta-QTL (MQTL) with reduced confidence interval (CI). Five GWAS validated MQTL (viz., 1A.1, 1B.2, 3B.4, 5B.2, and 6B.2), each involving more than 20 initial QTL and reduced CI (95%) (< 2 cM), were selected for quality breeding programmes. Functional characterization including candidate gene mining and expression analysis discovered 44 high confidence candidate genes associated with quality traits. A meta-analysis of quantitative trait loci (QTL) associated with dough rheology properties, nutritional traits, and processing quality traits was conducted in wheat. For this purpose, as many as 2458 QTL were collected from 50 interval mapping studies published during 2013-2020. Of the total QTL, 1126 QTL were projected onto the consensus map saturated with 249,603 markers which led to the identification of 110 meta-QTL (MQTL). These MQTL exhibited an 18.84-fold reduction in the average CI compared to the average CI of the initial QTL (ranging from 14.87 to 95.55 cM with an average of 40.35 cM). Of the 110, 108 MQTL were physically anchored to the wheat reference genome, including 51 MQTL verified with marker-trait associations (MTAs) reported from earlier genome-wide association studies. Candidate gene (CG) mining allowed the identification of 2533 unique gene models from the MQTL regions. In-silico expression analysis discovered 439 differentially expressed gene models with > 2 transcripts per million expressions in grains and related tissues, which also included 44 high-confidence CGs involved in the various cellular and biochemical processes related to quality traits. Nine functionally characterized wheat genes associated with grain protein content, high-molecular-weight glutenin, and starch synthase enzymes were also found to be co-localized with some of the MQTL. Synteny analysis between wheat and rice MQTL regions identified 23 wheat MQTL syntenic to 16 rice MQTL associated with quality traits. Furthermore, 64 wheat orthologues of 30 known rice genes were detected in 44 MQTL regions. Markers flanking the MQTL identified in the present study can be used for marker-assisted breeding and as fixed effects in the genomic selection models for improving the prediction accuracy during quality breeding. Wheat orthologues of rice genes and other CGs available from MQTLs can be promising targets for further functional validation and to better understand the molecular mechanism underlying the quality traits in wheat.

Targeting Induced Local Lesions in the Wheat DEMETER and DRE2 Genes, Responsible for Transcriptional Derepression of Wheat Gluten Proteins in the Developing Endosperm

Wheat is a major source of energy and nutrition worldwide, but it is also a primary cause of frequent diet-induced health issues, specifically celiac disease, for which the only effective therapy so far is strict dietary abstinence from gluten-containing grains. Wheat gluten proteins are grouped into two major categories: high-molecular-weight glutenin subunits (HMWgs), vital for mixing and baking properties, and gliadins plus low-molecular-weight glutenin subunits (LMWgs) that contain the overwhelming majority of celiac-causing epitopes. We put forth a hypothesis that eliminating gliadins and LMWgs while retaining HMWgs might allow the development of reduced-immunogenicity wheat genotypes relevant to most gluten-sensitive individuals. This hypothesis stems from the knowledge that the molecular structures and regulatory mechanisms of the genes encoding the two groups of gluten proteins are quite different, and blocking one group's transcription, without affecting the other's, is possible. The genes for gliadins and LMWgs have to be de-methylated by 5-methylcytosine DNA glycosylase/lyase (DEMETER) and an iron-sulfur (Fe-S) cluster biogenesis enzyme (DRE2) early during endosperm development to permit their transcription. In this study, a TILLING (Targeting Induced Local Lesions IN Genomes) approach was undertaken to identify mutations in the homoeologous DEMETER (DME) and DRE2 genes in common and durum wheat. Lines with mutations in these genes were obtained that displayed reduced content of immunogenic gluten proteins while retaining essential baking properties. Although our data at first glance suggest new possibilities for treating celiac disease and are therefore of medical and agronomical interest, it also shows that inducing mutations in the DME and DRE2 genes analyzed here affected pollen viability and germination. Hence there is a need to develop other approaches in the future to overcome this undesired effect.

Meta-analysis reveals consensus genomic regions associated with multiple disease resistance in wheat (Triticum aestivum L.)

In wheat, meta-QTLs (MQTLs) and candidate genes (CGs) were identified for multiple disease resistance (MDR). For this purpose, information was collected from 58 studies for mapping QTLs for resistance to one or more of the five diseases. As many as 493 QTLs were available from these studies, which were distributed in five diseases as follows: septoria tritici blotch (STB) 126 QTLs; septoria nodorum blotch (SNB), 103 QTLs; fusarium head blight (FHB), 184 QTLs; karnal bunt (KB), 66 QTLs; and loose smut (LS), 14 QTLs. Of these 493 QTLs, only 291 QTLs could be projected onto a consensus genetic map, giving 63 MQTLs. The CI of the MQTLs ranged from 0.04 to 15.31 cM with an average of 3.09 cM per MQTL. This is a ~ 4.39 fold reduction from the CI of QTLs, which ranged from 0 to 197.6 cM, with a mean of 13.57 cM. Of 63 MQTLs, 60 were anchored to the reference physical map of wheat (the physical interval of these MQTLs ranged from 0.30 to 726.01 Mb with an average of 74.09 Mb). Thirty-eight (38) of these MQTLs were verified using marker-trait associations (MTAs) derived from genome-wide association studies. As many as 874 CGs were also identified which were further investigated for differential expression using data from five transcriptome studies, resulting in 194 differentially expressed candidate genes (DECGs). Among the DECGs, 85 genes had functions previously reported to be associated with disease resistance. These results should prove useful for fine mapping and cloning of MDR genes and marker-assisted breeding.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-022-01282-z.

Conservation and Divergence of SQUAMOSA-PROMOTER BINDING PROTEIN-LIKE (SPL) Gene Family between Wheat and Rice

The SQUAMOSA-PROMOTER BINDING PROTEIN-LIKE (SPL) gene family affects plant architecture, panicle structure, and grain development, representing key genes for crop improvements. The objective of the present study is to utilize the well characterized SPLs' functions in rice to facilitate the functional genomics of TaSPL genes. To achieve these goals, we combined several approaches, including genome-wide analysis of TaSPLs, comparative genomic analysis, expression profiling, and functional study of TaSPL3 in rice. We established the orthologous relationships of 56 TaSPL genes with the corresponding OsSPLs, laying a foundation for the comparison of known SPL functions between wheat and rice. Some TaSPLs exhibited different spatial-temporal expression patterns when compared to their rice orthologs, thus implicating functional divergence. TaSPL2/6/8/10 were identified to respond to different abiotic stresses through the combination of RNA-seq and qPCR expression analysis. Additionally, ectopic expression of TaSPL3 in rice promotes heading dates, affects leaf and stem development, and leads to smaller panicles and decreased yields per panicle. In conclusion, our work provides useful information toward cataloging of the functions of TaSPLs, emphasized the conservation and divergence between TaSPLs and OsSPLs, and identified the important SPL genes for wheat improvement.

Genome-Wide Analysis of Serine Hydroxymethyltransferase Genes in Triticeae Species Reveals That TaSHMT3A-1 Regulates Fusarium Head Blight Resistance in Wheat

Serine hydroxymethyltransferase (SHMT) plays a pivotal role in cellular one-carbon, photorespiration pathways and it influences the resistance to biotic and abiotic stresses. However, the function of SHMT proteins in wheat remains largely unexplored. In the present study, SHMT genes in five Triticeae species, Oryza sativa, and four dicotyledon species were identified based on whole genome information. The origin history of the target gene was traced by micro-collinearity analysis. Gene expression patterns of TaSHMTs in different tissues, various biotic stresses, exogenous hormones, and two biotic stresses were determined by Quantitative real-time reverse transcription polymerase chain reaction (qRT-PCR). The function of the selected TaSHMT3A-1 was studied by barley stripe mosaic virus-induced gene silencing in common wheat Bainong207. A total of 64 SHMT members were identified and further classified into two main classes based on the structure of SHMT proteins. The gene structure and motif composition analyses revealed that SHMTs kept relatively conserved within the same subclasses. Interestingly, there was a gene, TdSHMT7B-1, on chromosome 7B of Triticum dicoccoides, but there was no SHMT gene on chromosome 7 of other analyzed Triticeae species; TdSHMT7B-1 had fewer exons and conserved motifs than the genes in the same subclass, suggesting that the gene of TdSHMT7B-1 has a notable evolutionary progress. The micro-collinearity relationship showed that no homologs of TaSHMT3A-1 and its two neighboring genes were found in the collinearity region of Triticum urartu, and there were 27 genes inserted into the collinearity region of T. urartu. Furthermore, qRT-PCR results showed that TaSHMT3A-1 was responsive to abiotic stresses (NaCl and cold), abscisic acid, methyl jasmonate, and hydrogen peroxide. Significantly, upon Fusarium graminearum infection, the expression of TaSHMT3A-1 was highly upregulated in resistant cultivar Sumai3. More importantly, silencing of TaSHMT3A-1 compromises Fusarium head blight resistance in common wheat Bainong207. Our new findings suggest that the TaSHMT3A-1 gene in wheat plays an important role in resistance to Fusarium head blight. This provides a valuable reference for further study on the function of this gene family.

Transcriptome Analysis of Developing Wheat Grains at Rapid Expanding Phase Reveals Dynamic Gene Expression Patterns

Simple Summary

Understanding the regulatory mechanism underlying grain development is essential for wheat improvement. The early grain expanding phase boasts critical biological events like embryogenesis and initiation of grain filling. RNA sequencing analysis of this developmental stage revealed dynamic expressions of genes related to cell division, starch biosynthesis, and hormone biosynthesis. An unbalanced expression among triads may play critical roles as shown by multiple enriched metabolic pathways. Our work demonstrated complex regulation mechanisms in early grain development and provided useful information for future wheat improvement.

Abstract

Grain development, as a vital process in the crop’s life cycle, is crucial for determining crop quality and yield. The wheat grain expanding phase is the early process involving the rapid morphological changes and initiation of grain filling. However, little is known about the molecular basis of grain development at this stage. Here, we provide a time-series transcriptome profile of developing wheat grain at 0, 2, 4, 6, 8, and 10 days after pollination of the wheat landrace Chinese Spring. A total of 26,892 differentially expressed genes, including 1468 transcription factors, were found between adjacent time points. Co-expression cluster analysis and Gene Ontology enrichment revealed dynamic expressions of cell division and starch biosynthesis related structural genes and transcription factors. Moreover, diverse, differential and drastically varied expression trends of the key genes related to hormone metabolism were identified. Furthermore, ~30% of triads showed unbalanced expression patterns enriching for genes in multiple pivotal metabolic pathways. Hormone metabolism related genes, such as YUC10 (YUCCA flavin-containing monooxygenase 10), AOS2 (allene oxide synthase 2), CYP90D2 (cytochrome P450 90D2), and CKX1 (cytokinin dehydrogenase 1), were dominantly contributed by A or D homoeologs of the triads. Our study provided a systematic picture of transcriptional regulation of wheat grains at the early grain expanding phase which should deepen our understanding of wheat grain development and help in wheat yield improvement.

Gene expression profiling reveals transcription factor networks and subgenome bias during Brassica napus seed development

We profiled the global gene expression landscape across the reproductive lifecycle of Brassica napus. Comparative analysis of this nascent amphidiploid revealed the contribution of each subgenome to plant reproduction. Whole-genome transcription factor networks identified BZIP11 as a transcriptional regulator of early B. napus seed development. Knockdown of BZIP11 using RNA interference resulted in a similar reduction in gene activity of predicted gene targets, and a reproductive-lethal phenotype. Global mRNA profiling revealed lower accumulation of Cⁿ subgenome transcripts relative to the Aⁿ subgenome. Subgenome-specific transcription factor networks identified distinct transcription factor families enriched in each of the Aⁿ and Cⁿ subgenomes early in seed development. Analysis of laser-microdissected seed subregions further reveal subgenome expression dynamics in the embryo, endosperm and seed coat of early stage seeds. Transcription factors predicted to be regulators encoded by the Aⁿ subgenome are expressed primarily in the seed coat, whereas regulators encoded by the Cⁿ subgenome were expressed primarily in the embryo. Data suggest subgenome bias are characteristic features of the B. napus seed throughout development, and that such bias might not be universal across the embryo, endosperm and seed coat of the developing seed. Transcriptional networks spanning both the Aⁿ and Cⁿ genomes of the whole B. napus seed can identify valuable targets for seed development research and that -omics level approaches to studying gene regulation in B. napus can benefit from both broad and high-resolution analyses.

Protein turnover in the developing Triticum aestivum grain

Protein abundance in cereal grains is determined by the relative rates of protein synthesis and protein degradation during grain development but quantitation of these rates is lacking. Through combining in vivo stable isotope labelling and in-depth quantitative proteomics, we have measured the turnover of 1400 different types of proteins during wheat grain development. We demonstrate that there is a spatiotemporal pattern to protein turnover rates which explain part of the variation in protein abundances that is not attributable to differences in wheat gene expression. We show that c. 20% of total grain adenosine triphosphate (ATP) production is used for grain proteome biogenesis and maintenance, and nearly half of this budget is invested exclusively in storage protein synthesis. We calculate that 25% of newly synthesized storage proteins are turned over during grain development rather than stored. This approach to measure protein turnover rates at proteome scale reveals how different functional categories of grain proteins accumulate, calculates the costs of protein turnover during wheat grain development and identifies the most and the least stable proteins in the developing wheat grain.

Spatial and Temporal Transcriptomic Heredity and Asymmetry in an Artificially Constructed Allotetraploid Wheat (AADD)

Polyploidy, or whole-genome duplication (WGD), often induces dramatic changes in gene expression due to "transcriptome shock. " However, questions remain about how allopolyploidy (the merging of multiple nuclear genomes in the same nucleus) affects gene expression within and across multiple tissues and developmental stages during the initial foundation of allopolyploid plants. Here, we systematically investigated the immediate effect of allopolyploidy on gene expression variation in an artificial allopolyploidy system consisting of a constructed allotetraploid wheat (AADD genome, accession AT2) and its diploid progenitors Triticum urartu and Aegilops tauschii. We performed comprehensive RNA sequencing of 81 samples from different genotypes, tissues, and developmental stages. First, we found that intrinsic interspecific differences between the diploid parents played a major role in establishing the expression architecture of the allopolyploid. Nonetheless, allopolyploidy per se also induced dramatic and asymmetric patterns of differential gene expression between the subgenomes, and genes from the D subgenome exhibited a more drastic response. Second, analysis of homoeolog expression bias (HEB) revealed that the D subgenome exhibited significant expression bias and that de novo-generated HEB was attributed mainly to asymmetrical differential gene expression. Homoeolog-specific expression (HSE) analyses showed that the cis-only regulatory pattern was predominant in AT2, reflecting significant divergence between the parents. Co-expression network analysis revealed that homoeolog expression connectivity (HEC) was significantly correlated with sequence divergence in cis elements between subgenomes. Interestingly, allopolyploidy-induced reconstruction of network modules was also associated with different HSE patterns. Finally, a transcriptome atlas of spike development demonstrated that the phenotypic similarity of AT2 to T. urartu may be attributed to the combination of relatively stable expression of A-subgenome genes and drastic downregulation of their D-subgenome homoeologs. These findings provide a broad, multidimensional characterization of allopolyploidy-induced transcriptomic responses and suggest that allopolyploidy can have immediate and complex regulatory effects on the expression of nuclear genes.

Genome-Wide Identification and Expression Profile of OSCA Gene Family Members in Triticum aestivum L

Hyperosmolality and various other stimuli can trigger an increase in cytoplasmic-free calcium concentration ([Ca²⁺]_cyt). Members of the Arabidopsis thaliana (L.) reduced hyperosmolality-gated calcium-permeable channels (OSCA) gene family are reported to be involved in sensing extracellular changes to trigger hyperosmolality-induced [Ca²⁺]_cyt increases and controlling stomatal closure during immune signaling. Wheat (Triticum aestivum L.) is a very important food crop, but there are few studies of its OSCA gene family members. In this study, 42 OSCA members were identified in the wheat genome, and phylogenetic analysis can divide them into four clades. The members of each clade have similar gene structures, conserved motifs, and domains. TaOSCA genes were predicted to be regulated by cis-acting elements such as STRE, MBS, DRE1, ABRE, etc. Quantitative PCR results showed that they have different expression patterns in different tissues. The expression profiles of 15 selected TaOSCAs were examined after PEG (polyethylene glycol), NaCl, and ABA (abscisic acid) treatment. All 15 TaOSCA members responded to PEG treatment, while TaOSCA12/-39 responded simultaneously to PEG and ABA. This study informs research into the biological function and evolution of TaOSCA and lays the foundation for the breeding and genetic improvement of wheat.

[Ribosomal profiling as a tool for studying translation in plants: main results, problems and future prospects]

The expression of eukaryotic genes can be regulated at several stages, including the translation of mRNA. It is known that the structure of mRNA can affect both the efficiency of interaction with the translation apparatus in general and the choice of translation initiation sites. To study the translated fraction of the transcriptome, experimental methods of analysis were developed, the most informative of which is ribosomal profiling (RP, Ribo-seq). Originally developed for use in yeast systems, this method has been adapted for research in translation mechanisms in many plant species. This technology includes the isolation of the polysomal fraction and high-performance sequencing of a pool of mRNA fragments associated with ribosomes. Comparing the results of transcript coverage with reads obtained using the ribosome profiling with the transcriptional efficiency of genes allows the translation efficiency to be evaluated for each transcript. The exact positions of ribosomes determined on mRNA sequences allow determining the translation of open reading frames and switching between the translation of several reading frames - a phenomenon in which two or more overlapping frames are read from one mRNA and different proteins are synthesized. The advantage of this method is that it provides quantitative estimates of ribosome coverage of mRNA and can detect relatively rare translation events. Using this technology, it was possible to identify and classify plant genes by the type of regulation of their expression at the transcription, translation, or both levels. Features of the mRNA structure that affect translation levels have been revealed: the formation of G2 quadruplexes and the presence of specific motifs in the 5'-UTR region, GC content, the presence of alternative translation starts, and the influence of uORFs on the translation of downstream mORFs. In this review, we briefly reviewed the RP methodology and the prospects for its application to study the structural and functional organization and regulation of plant gene expression.

Asymmetric Divergence in Transmitted SNPs of DNA Replication/Transcription and Their Impact on Gene Expression in Polyploid Brassica napus

The marked increase in plant genomic data has provided valuable resources for investigating the dynamic evolution of duplicate genes in polyploidy. Brassica napus is an ideal model species for investigating polyploid genome evolution. The present study comprehensively analyzed DNA and RNA variation of two representative B. napus inbredlines, Zhongshuang11 and Zhongyou821, and we investigated gene expression levels of A_n and C_n subgenomes in multiple tissues of the two lines. The distribution of transmitted single nucleotide polymorphisms (SNPs) was significantly different in two subgenomes of B. napus. Gene expression levels were significantly negatively correlated with number of variations in replication and transcription of the corresponding genes, but were positively correlated with the ratios of transmitted SNPs from DNA to RNA. We found a higher density of SNP variation in A_n than that in C_n during DNA replication and more SNPs were transmitted to RNA during transcription, which may contribute to A_n expression dominance. These activities resulted in asymmetrical gene expression in polyploid B. napus. The SNPs transmitted from DNA to RNA could be an important complement feature in comparative genomics, and they may play important roles in asymmetrical genome evolution in polyploidy.

Exploring the diversity of promoter and 5'UTR sequences in ancestral, historic and modern wheat

A data set of promoter and 5'UTR sequences of homoeo-alleles of 459 wheat genes that contribute to agriculturally important traits in 95 ancestral and commercial wheat cultivars is presented here. The high-stringency myBaits technology used made individual capture of homoeo-allele promoters possible, which is reported here for the first time. Promoters of most genes are remarkably conserved across the 83 hexaploid cultivars used with <7 haplotypes per promoter and 21% being identical to the reference Chinese Spring. InDels and many high-confidence SNPs are located within predicted plant transcription factor binding sites, potentially changing gene expression. Most haplotypes found in the Watkins landraces and a few haplotypes found in Triticum monococcum, germplasms hitherto not thought to have been used in modern wheat breeding, are already found in many commercial hexaploid wheats. The full data set which is useful for genomic and gene function studies and wheat breeding is available at https://rrescloud.rothamsted.ac.uk/index.php/s/DMCFDu5iAGTl50u/authenticate.