A Genome-wide View of Transcriptome Dynamics During Early Spike Development in Bread Wheat

Identification of molecular markers and candidate regions associated with grain number per spike in Pubing3228 using SLAF-BSA

Grain number per spike, a pivotal agronomic trait dictating wheat yield, lacks a comprehensive understanding of its underlying mechanism in Pubing3228, despite the identification of certain pertinent genes. Thus, our investigation sought to ascertain molecular markers and candidate regions associated with grain number per spike through a high-density genetic mapping approach that amalgamates site-specific amplified fragment sequencing (SLAF-seq) and bulked segregation analysis (BSA). To facilitate this, we conducted a comparative analysis of two wheat germplasms, Pubing3228 and Jing4839, known to exhibit marked discrepancies in spike shape. By leveraging this methodology, we successfully procured 2,810,474 SLAF tags, subsequently resulting in the identification of 187,489 single nucleotide polymorphisms (SNPs) between the parental strains. We subsequently employed the SNP-index association algorithm alongside the extended distribution (ED) association algorithm to detect regions associated with the trait. The former algorithm identified 24 trait-associated regions, whereas the latter yielded 70. Remarkably, the intersection of these two algorithms led to the identification of 25 trait-associated regions. Amongst these regions, we identified 399 annotated genes, including three genes harboring non-synonymous mutant SNP loci. Notably, the APETALA2 (AP2) transcription factor families, which exhibited a strong correlation with spike type, were also annotated. Given these findings, it is plausible to hypothesize that these genes play a critical role in determining spike shape. In summation, our study contributes significant insights into the genetic foundation of grain number per spike. The molecular markers and candidate regions we have identified can be readily employed for marker-assisted breeding endeavors, ultimately leading to the development of novel wheat cultivars possessing enhanced yield potential. Furthermore, conducting further functional analyses on the identified genes will undoubtedly facilitate a comprehensive elucidation of the underlying mechanisms governing spike development in wheat.

Integrated transcriptomic and metabolomic analyses elucidate the mechanism of flavonoid biosynthesis in the regulation of mulberry seed germination under salt stress

Seed propagation is the main method of mulberry expansion in China, an important economic forest species. However, seed germination is the most sensitive stage to various abiotic stresses, especially salinity stress. To reveal the molecular regulatory mechanism of mulberry seed germination under salt stress, flavonoid metabolomics and transcriptomics analyses were performed on mulberry seeds germinated under 50 and 100 mmol/L NaCl stress. Analysis of the flavonoid metabolome revealed that a total of 145 differential flavonoid metabolites (DFMs) were classified into 9 groups, 40 flavonols, 32 flavones, 16 chalcones and 14 flavanones. Among them, 61.4% (89) of the DFMs accumulated continuously with increasing salt concentration, reaching the highest level at a 100 mmol/L salt concentration; these DFMs included quercetin-3-O-glucoside (isoquercitrin), kaempferol (3,5,7,4'-tetrahydroxyflavone), quercetin-7-O-glucoside, taxifolin (dihydroquercetin) and apigenin (4',5,7-trihydroxyflavone), indicating that these flavonoids may be key metabolites involved in the response to salt stress. Transcriptional analysis identified a total of 3055 differentially expressed genes (DEGs), most of which were enriched in flavonoid biosynthesis (ko00941), phenylpropanoid biosynthesis (ko00940) and biosynthesis of secondary metabolites (ko01110). Combined analysis of flavonoid metabolomic and transcriptomic data indicated that phenylalanine ammonia-lyase (PAL), 4-coumarate-CoA ligase (4CL), chalcone synthase (CHS), flavonol synthase (FLS), bifunctional dihydroflavonol 4-reductase/flavanone 4-reductase (DFR) and anthocyanidin reductase (ANR) were the key genes involved in flavonoid accumulation during mulberry seed germination under 50 and 100 mmol/L NaCl stress. In addition, three transcription factors, MYB, bHLH and NAC, were involved in the regulation of flavonoid accumulation under salt stress. The results of quantitative real-time PCR (qRT‒PCR) validation showed that the expression levels of 11 DEGs, including 7 genes involved in flavonoid biosynthesis, under different salt concentrations were consistent with the transcriptomic data, and parallel reaction monitoring (PRM) results showed that the expression levels of 6 key enzymes (proteins) involved in flavonoid synthesis were consistent with the accumulation of flavonoids. This study provides a new perspective for investigating the regulatory role of flavonoid biosynthesis in the regulation of mulberry seed germination under salt stress at different concentrations.

Auxins and grass shoot architecture: how the most important hormone makes the most important plants

Cereals are a group of grasses cultivated by humans for their grain. It is from these cereal grains that the majority of all calories consumed by humans are derived. The production of these grains is the result of the development of a series of hierarchical reproductive structures that form the distinct shoot architecture of the grasses. Being spatiotemporally complex, the coordination of grass shoot development is tightly controlled by a network of genes and signals, including the key phytohormone auxin. Hormonal manipulation has therefore been identified as a promising potential approach to increasing cereal crop yields and therefore ultimately global food security. Recent work translating the substantial body of auxin research from model plants into cereal crop species is revealing the contribution of auxin biosynthesis, transport, and signalling to the development of grass shoot architecture. This review discusses this still-maturing knowledge base and examines the possibility that changes in auxin biology could have been a causative agent in the evolution of differences in shoot architecture between key grass species, or could underpin the future selective breeding of cereal crops.

Deciphering spike architecture formation towards yield improvement in wheat

Wheat is the most widely grown crop globally, providing 20% of the daily consumed calories and protein content around the world. With the growing global population and frequent occurrence of extreme weather caused by climate change, ensuring adequate wheat production is essential for food security. The architecture of the inflorescence plays a crucial role in determining the grain number and size, which is a key trait for improving yield. Recent advances in wheat genomics and gene cloning techniques have improved our understanding of wheat spike development and its applications in breeding practices. Here, we summarize the genetic regulation network governing wheat spike formation, the strategies used for identifying and studying the key factors affecting spike architecture, and the progress made in breeding applications. Additionally, we highlight future directions that will aid in the regulatory mechanistic study of wheat spike determination and targeted breeding for grain yield improvement.

Transcriptome Dynamics during Spike Differentiation of Wheat Reveal Amazing Changes in Cell Wall Metabolic Regulators

Coordinated cell proliferation and differentiation result in the complex structure of the inflorescence in wheat. It exhibits unique differentiation patterns and structural changes at different stages, which have attracted the attention of botanists studying the dynamic regulation of its genes. Our research aims to understand the molecular mechanisms underlying the regulation of spike development genes at different growth stages. We conducted RNA-Seq and qRT-PCR evaluations on spikes at three stages. Our findings revealed that genes associated with the cell wall and carbohydrate metabolism showed high expression levels between any two stages throughout the entire process, suggesting their regulatory role in early spike development. Furthermore, through transgenic experiments, we elucidated the role of the cell wall regulator gene in spike development regulation. These research results contribute to identifying essential genes associated with the morphology and development of wheat spike tissue.

Boosting wheat functional genomics via an indexed EMS mutant library of KN9204

A better understanding of wheat functional genomics can improve targeted breeding for better agronomic traits and environmental adaptation. However, the lack of gene-indexed mutants and the low transformation efficiency of wheat limit in-depth gene functional studies and genetic manipulation for breeding. In this study, we created a library for KN9204, a popular wheat variety in northern China, with a reference genome, transcriptome, and epigenome of different tissues, using ethyl methyl sulfonate (EMS) mutagenesis. This library contains a vast developmental diversity of critical tissues and transition stages. Exome capture sequencing of 2090 mutant lines using KN9204 genome-designed probes revealed that 98.79% of coding genes had mutations, and each line had an average of 1383 EMS-type SNPs. We identified new allelic variations for crucial agronomic trait-related genes such as Rht-D1, Q, TaTB1, and WFZP. We tested 100 lines with severe mutations in 80 NAC transcription factors (TFs) under drought and salinity stress and identified 13 lines with altered sensitivity. Further analysis of three lines using transcriptome and chromatin accessibility data revealed hundreds of direct NAC targets with altered transcription patterns under salt or drought stress, including SNAC1, DREB2B, CML16, and ZFP182, factors known to respond to abiotic stress. Thus, we have generated and indexed a KN9204 EMS mutant library that can facilitate functional genomics research and offer resources for genetic manipulation of wheat.

A framework for improving wheat spike development and yield based on the master regulatory TOR and SnRK gene systems

The low rates of yield gain in wheat breeding programs create an ominous situation for the world. Amongst the reasons for this low rate are issues manifested in spike development that result in too few spikelets, fertile florets, and therefore grains being produced. Phases in spike development are particularly sensitive to stresses of various kinds and origins, and these are partly responsible for the deficiencies in grain production and slow rates of gain in yield. The diversity of developmental processes, stresses, and the large numbers of genes involved make it particularly difficult to prioritize approaches in breeding programs without an overarching, mechanistic framework. Such a framework, introduced here, is provided around the master regulator target of rapamycin and sucrose non-fermenting-1 (SNF1)-related protein kinase complexes and their control by trehalose-6-phosphate and other molecules. Being master regulators of the balance between growth and growth inhibition under stress, these provide genetic targets for creating breakthroughs in yield enhancement. Examples of potential targets and experimental approaches are described.

A 'wiring diagram' for sink strength traits impacting wheat yield potential

Identifying traits for improving sink strength is a bottleneck to increasing wheat yield. The interacting processes determining sink strength and yield potential are reviewed and visualized in a set of 'wiring diagrams', covering critical phases of development (and summarizing known underlying genetics). Using this framework, we reviewed and assembled the main traits determining sink strength and identified research gaps and potential hypotheses to be tested for achieving gains in sink strength. In pre-anthesis, grain number could be increased through: (i) enhanced spike growth associated with optimized floret development and/or a reduction in specific stem-internode lengths and (ii) improved fruiting efficiency through an accelerated rate of floret development, improved partitioning between spikes, or optimized spike cytokinin levels. In post-anthesis, grain, sink strength could be augmented through manipulation of grain size potential via ovary size and/or endosperm cell division and expansion. Prospects for improving spike vascular architecture to support all rapidly growing florets, enabling the improved flow of assimilate, are also discussed. Finally, we considered the prospects for enhancing grain weight realization in relation to genetic variation in stay-green traits as well as stem carbohydrate remobilization. The wiring diagrams provide a potential workspace for breeders and crop scientists to achieve yield gains in wheat and other field crops.

Integrated metabolomics and transcriptomic analysis of the flavonoid regulatory networks in Sorghum bicolor seeds

Background

The objective of this study was to reveal the flavonoid biosynthesis pathway in white (Z6), red (Z27) and black (HC4) seeds of the sweet sorghum (Sorghum bicolor) using metabolomics and transcriptomics, to identify different flavonoid metabolites, and to analyze the differentially expressed genes involved in flavonoid biosynthesis.

Results

We analyzed the metabolomics and transcriptomics data of sweet sorghum seeds. Six hundred and fifty-one metabolites including 171 flavonoids were identified in three samples. Integrated analysis of transcriptomics and metabolomics showed that 8 chalcone synthase genes (gene19114, gene19115, gene19116, gene19117, gene19118, gene19120, gene19122 and gene19123) involved in flavonoid biosynthesis, were identified and play central role in change of color. Six flavanone including homoeriodictyol, naringin, prunin, naringenin, hesperetin and pinocembrin were main reason for the color difference.

Conclusions

Our results provide valuable information on the flavonoid metabolites and the candidate genes involved in the flavonoid biosynthesis pathway in sweet sorghum seeds.

The Central Circadian Clock Protein TaCCA1 Regulates Seedling Growth and Spike Development in Wheat (Triticum aestivum L.)

The biological functions of the circadian clock on growth and development have been well elucidated in model plants, while its regulatory roles in crop species, especially the roles on yield-related traits, are poorly understood. In this study, we characterized the core clock gene CIRCADIAN CLOCK-ASSOCIATED 1 (CCA1) homoeologs in wheat and studied their biological functions in seedling growth and spike development. TaCCA1 homoeologs exhibit typical diurnal expression patterns, which are positively regulated by rhythmic histone modifications including histone H3 lysine 4 trimethylation (H3K4me3), histone H3 lysine 9 acetylation (H3K9Ac), and histone H3 lysine 36 trimethylation (H3K36me3). TaCCA1s are preferentially located in the nucleus and tend to form both homo- and heterodimers. TaCCA1 overexpression (TaCCA1-OE) transgenic wheat plants show disrupted circadian rhythmicity coupling with reduced chlorophyll and starch content, as well as biomass at seedling stage, also decreased spike length, grain number per spike, and grain size at the ripening stage. Further studies using DNA affinity purification followed by deep sequencing [DNA affinity purification and sequencing (DAP-seq)] indicated that TaCCA1 preferentially binds to sequences similarly to "evening elements" (EE) motif in the wheat genome, particularly genes associated with photosynthesis, carbon utilization, and auxin homeostasis, and decreased transcriptional levels of these target genes are observed in TaCCA1-OE transgenic wheat plants. Collectively, our study provides novel insights into a circadian-mediated mechanism of gene regulation to coordinate photosynthetic and metabolic activities in wheat, which is important for optimal plant growth and crop yield formation.

Population Genomics and Haplotype Analysis in Bread Wheat Identify a Gene Regulating Glume Pubescence

Glume hairiness or pubescence is an important morphological trait with high heritability to distinguish/characterize wheat and is related to the resistance to biotic and abiotic stresses. Hg1 (formerly named Hg) on chromosome arm 1AS controlled glume hairiness in wheat. Its genetic analysis and mapping have been widely studied, yet more useful and accurate information for fine mapping of Hg1 and identification of its candidate gene is lacking. The cloning of this gene has not yet been reported for the large complex wheat genome. Here, we performed a GWAS between SNP markers and glume pubescence (Gp) in a wheat population with 352 lines and further demonstrated the gene expression and haplotype analysis approach for isolating the Hg1 gene. One gene, TraesCSU02G143200 (TaELD1-1A), encoding glycosyltransferase-like ELD1/KOBITO 1, was identified as the most promising candidate gene of Hg1. The gene annotation, expression pattern, function SNP variation, haplotype analysis, and co-expression analysis in floral organ (spike) development indicated that it is likely to be involved in the regulation of glume pubescence. Our study demonstrates the importance of high-quality reference genomes and annotation information, as well as bioinformatics analysis, for gene cloning in wheat.

Capturing Wheat Phenotypes at the Genome Level

Recent technological advances in next-generation sequencing (NGS) technologies have dramatically reduced the cost of DNA sequencing, allowing species with large and complex genomes to be sequenced. Although bread wheat (Triticum aestivum L.) is one of the world's most important food crops, efficient exploitation of molecular marker-assisted breeding approaches has lagged behind that achieved in other crop species, due to its large polyploid genome. However, an international public-private effort spanning 9 years reported over 65% draft genome of bread wheat in 2014, and finally, after more than a decade culminated in the release of a gold-standard, fully annotated reference wheat-genome assembly in 2018. Shortly thereafter, in 2020, the genome of assemblies of additional 15 global wheat accessions was released. As a result, wheat has now entered into the pan-genomic era, where basic resources can be efficiently exploited. Wheat genotyping with a few hundred markers has been replaced by genotyping arrays, capable of characterizing hundreds of wheat lines, using thousands of markers, providing fast, relatively inexpensive, and reliable data for exploitation in wheat breeding. These advances have opened up new opportunities for marker-assisted selection (MAS) and genomic selection (GS) in wheat. Herein, we review the advances and perspectives in wheat genetics and genomics, with a focus on key traits, including grain yield, yield-related traits, end-use quality, and resistance to biotic and abiotic stresses. We also focus on reported candidate genes cloned and linked to traits of interest. Furthermore, we report on the improvement in the aforementioned quantitative traits, through the use of (i) clustered regularly interspaced short-palindromic repeats/CRISPR-associated protein 9 (CRISPR/Cas9)-mediated gene-editing and (ii) positional cloning methods, and of genomic selection. Finally, we examine the utilization of genomics for the next-generation wheat breeding, providing a practical example of using in silico bioinformatics tools that are based on the wheat reference-genome sequence.

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

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

Genetic Dissection of Three Major Quantitative Trait Loci for Spike Compactness and Length in Bread Wheat (Triticum aestivum L.)

Spike compactness (SC) and length (SL) are the components of spike morphology and are strongly related to grain yield in wheat (Triticum aestivum L.). To investigate quantitative trait loci (QTL) associated with SC and SL, a recombinant inbred lines (RIL) population derived from the cross of Bailangmai (BLM, a Tibet landrace) and Chuanyu 20 (CY20, an improved variety) was employed in six environments. Three genomic regions responsible for SC and SL traits were identified on chromosomes 2A and 2D using bulked segregant exome sequencing (BSE-Seq). By constructing genetic maps, six major QTL were repeatedly detected in more than four environments and the best linear unbiased estimation (BLUE) datasets, explaining 7.00-28.56% of the phenotypic variation and the logarithm of the odd (LOD) score varying from 2.50 to 13.22. They were co-located on three loci, designed as QSc/Sl.cib-2AS, QSc/Sl.cib-2AL, and QSc/Sl.cib-2D, respectively. Based on the flanking markers, their interactions and effects on the corresponding trait and other agronomic traits were also analyzed. Comparison analysis showed that QSc/Sl.cib-2AS and QSc/Sl.cib-2AL were possibly two novel loci for SC and SL. QSc/Sl.cib-2AS and QSc/Sl.cib-2D showed pleiotropic effects on plant height and grain morphology, while QSc/Sl.cib-2AL showed effects on spikelet number per spike (SNS) and grain width (GW). Based on the gene annotation, orthologous search, and spatiotemporal expression patterns of genes, TraesCS2A03G0410600 and TraesCS2A03G0422300 for QSc/Sl.cib-2AS, and TraesCS2D03G1129300 and TraesCS2D03G1131500 for QSc/Sl.cib-2D were considered as potential candidate genes, respectively. These results will be useful for fine mapping and developing new varieties with high yield in the future.

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.

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.

Wheat Escapes Low Light Stress by Altering Pollination Types

Although low light stress seriously affects florets fertility and grain number during the reproductive period, crops can be fertilized by heterologous pollen to alleviate the reduction of grain number. However, wheat is strongly autogamous, how to change to outcross after low light remains unclear. To understand the mechanisms of this change process, an approach combined morphological, physiological, and transcriptomic analyses was performed under low light stress imposed at the young microspore stage the booting stage from tetrad to uni-nucleate microspores stage. The results showed that low light stress caused pollen abortion, and the unfertilized ovary is fertilized by heterologous pollen after floret opening. Compared to control, the opening angle of lemma and glume were increased by 11.6-48.6 and 48.4-78.5%, respectively. The outcross of stressed wheat compensated for the 2.1-18.0% of grain number loss. During this process, phytohormones played an important role. Jasmonic acid (JA) and methyl jasmonate (MeJA) levels in spikelets were increased. Meanwhile, lignin and cellulose content decreased, and genes associated with cell wall related GO terms were enriched. Among the differentially expressed genes (DEGs), were identified 88-710 transcription factors genes, of which some homologs in Arabidopsis are proposed to function in lignin and cellulose, influencing the glume and lemma opening. Our finding can provide new insight into a survival mechanism to set seeds through pollination way alteration in the absence of self-fertilization after the stress of adversity.

Identification and validation of two major QTLs for spike compactness and length in bread wheat (Triticum aestivum L.) showing pleiotropic effects on yield-related traits

Two major and stable QTLs for spike compactness and length were detected and validated in multiple genetic backgrounds and environments, and their pleiotropic effects on yield-related traits were analyzed. Spike compactness (SC) and length (SL) are greatly associated with wheat (Triticum aestivum L.) grain yield. To detect quantitative trait loci (QTL) associated with SC and SL, two biparental populations derived from crosses of Chuanmai42/Kechengmai1 and Chuanmai42/Chuannong16 were employed to perform QTL mapping in five environments. A total of 34 QTLs were identified, in which six major QTLs were repeatedly detected in more than four environments and the best linear unbiased prediction datasets, explaining 7.13-33.6% of phenotypic variation. These major QTLs were co-located in two genomic regions on chromosome 5A and 6A, namely QSc/Sl.cib-5A and QSc/Sl.cib-6A, respectively. By developing kompetitive allele-specific PCR (KASP) markers that linked to them, the two loci were validated in different genetic backgrounds, and their interactions were also analyzed. Comparison analysis showed that QSc/Sl.cib-5A was not Vrn-A1 and Q, and QSc/Sl.cib-6A was likely a new locus for SC and SL. Both QSc/Sl.cib-5A and QSc/Sl.cib-6A had pleiotropic effects on other yield-related traits including plant height, thousand grain weight and grain length. Therefore, the two loci combined with the developed KASP markers might be potentially applicable in wheat breeding. Furthermore, based on the spatiotemporal expression patterns, gene annotation, orthologous search and sequence differences, TraesCS5A01G301400 and TraesCS6A01G090300 were considered as potential candidates for QSc/Sl.cib-5A and QSc/Sl.cib-6A, respectively. These results provided valuable information for fine mapping and cloning of the two loci in the future.

Reproductive Stage Drought Tolerance in Wheat: Importance of Stomatal Conductance and Plant Growth Regulators

Drought stress requires plants to adjust their water balance to maintain tissue water levels. Isohydric plants (‘water-savers’) typically achieve this through stomatal closure, while anisohydric plants (‘water-wasters’) use osmotic adjustment and maintain stomatal conductance. Isohydry or anisohydry allows plant species to adapt to different environments. In this paper we show that both mechanisms occur in bread wheat (Triticum aestivum L.). Wheat lines with reproductive drought-tolerance delay stomatal closure and are temporarily anisohydric, before closing stomata and become isohydric at higher threshold levels of drought stress. Drought-sensitive wheat is isohydric from the start of the drought treatment. The capacity of the drought-tolerant line to maintain stomatal conductance correlates with repression of ABA synthesis in spikes and flag leaves. Gene expression profiling revealed major differences in the drought response in spikes and flag leaves of both wheat lines. While the isohydric drought-sensitive line enters a passive growth mode (arrest of photosynthesis, protein translation), the tolerant line mounts a stronger stress defence response (ROS protection, LEA proteins, cuticle synthesis). The drought response of the tolerant line is characterised by a strong response in the spike, displaying enrichment of genes involved in auxin, cytokinin and ethylene metabolism/signalling. While isohydry may offer advantages for longer term drought stress, anisohydry may be more beneficial when drought stress occurs during the critical stages of wheat spike development, ultimately improving grain yield.

Genome-wide profiling of long noncoding RNAs involved in wheat spike development

BACKGROUND

Long noncoding RNAs (lncRNAs) have been shown to play important roles in the regulation of plant growth and development. Recent transcriptomic analyses have revealed the gene expression profiling in wheat spike development, however, the possible regulatory roles of lncRNAs in wheat spike morphogenesis remain largely unclear.

RESULTS

Here, we analyzed the genome-wide profiling of lncRNAs during wheat spike development at six stages, and identified a total of 8,889 expressed lncRNAs, among which 2,753 were differentially expressed lncRNAs (DE lncRNAs) at various developmental stages. Three hundred fifteen differentially expressed cis- and trans-regulatory lncRNA-mRNA pairs comprised of 205 lncRNAs and 279 genes were predicted, which were found to be mainly involved in the stress responses, transcriptional and enzymatic regulations. Moreover, the 145 DE lncRNAs were predicted as putative precursors or target mimics of miRNAs. Finally, we identified the important lncRNAs that participate in spike development by potentially targeting stress response genes, TF genes or miRNAs.

CONCLUSIONS

This study outlines an overall view of lncRNAs and their possible regulatory networks during wheat spike development, which also provides an alternative resource for genetic manipulation of wheat spike architecture and thus yield.

Positional-based cloning 'fail-safe' approach is overpowered by wheat chromosome structural variation

Positional-based cloning is a foundational method for understanding the genes and gene networks that control valuable agronomic traits such as grain yield components. In this study, we sought to positionally clone the causal genetic variant of a 1000-grain weight (TGW) quantitative trait loci (QTL) on wheat (Triticum aestivum L.) chromosome arm 5AL. We developed heterogenous inbred families (HIFs) (>5,000 plants) for enhanced genotypic resolution and fine-mapped the QTL to a 10-Mbp region. The transcriptome of developing grains from positive and negative control HIF haplotypes revealed presence-absence chromosome arm 5AS structural variation and unexpectedly no differential expression of genes within the chromosome arm 5AL candidate region. Evaluation of genomic, transcriptomic, and phenotypic data, and predicted function of genes, identified that the 5AL QTL was the result of strong linkage disequilibrium (LD) with chromosome arm 5AS presence or absence (HIF r² = 0.91). Structural variation is common in wheat, and our results highlight that the redundant polyploid genome's masking of such variation is a significant barrier to positional cloning. We propose recommendations for more efficient and robust detection of structural variation, including transitioning from a single nucleotide polymorphism (SNP) to a haplotype-based approach to identify positional cloning targets. We also present nine candidate genes for grain yield components based on chromosome arm 5AS presence or absence, which may unveil hidden variation of homoeolog dosage-dependent genes across the group five chromosome short arms. Taken together, our discovery demonstrates the phenotypic resiliency of polyploid genomic structural variation and highlights a considerable challenge to routine positional cloning in wheat.

Wheat FRIZZY PANICLE activates VERNALIZATION1-A and HOMEOBOX4-A to regulate spike development in wheat

Kernel number per spike determined by the spike or inflorescence development is one important agricultural trait for wheat yield that is critical for global food security. While a few important genes for wheat spike development were identified, the genetic regulatory mechanism underlying supernumerary spikelets (SSs) is still unclear. Here, we cloned the wheat FRIZZY PANICLE (WFZP) gene from one local wheat cultivar. WFZP is specifically expressed at the sites where the spikelet meristem and floral meristem are initiated, which differs from the expression patterns of its homologs FZP/BD1 in rice and maize, indicative of its functional divergence during species differentiation. Moreover, WFZP directly activates VERNALIZATION1 (VRN1) and wheat HOMEOBOX4 (TaHOX4) to regulate the initiation and development of spikelet. The haplotypes analysis showed that the favourable alleles of WFZP associated with spikelet number per spike (SNS) were preferentially selected during breeding. Our findings provide insights into the molecular and genetic mechanisms underlying wheat spike development and characterize the WFZP as elite resource for wheat molecular breeding with enhanced crop yield.

Identification of eight QTL controlling multiple yield components in a German multi-parental wheat population, including Rht24, WAPO-A1, WAPO-B1 and genetic loci on chromosomes 5A and 6A

KEY MESSAGE

Quantitative trait locus (QTL) mapping of 15 yield component traits in a German multi-founder population identified eight QTL each controlling ≥2 phenotypes, including the genetic loci Rht24, WAPO-A1 and WAPO-B1. Grain yield in wheat (Triticum aestivum L.) is a polygenic trait representing the culmination of many developmental processes and their interactions with the environment. Toward maintaining genetic gains in yield potential, 'reductionist approaches' are commonly undertaken by which the genetic control of yield components, that collectively determine yield, are established. Here we use an eight-founder German multi-parental wheat population to investigate the genetic control and phenotypic trade-offs between 15 yield components. Increased grains per ear was significantly positively correlated with the number of fertile spikelets per ear and negatively correlated with the number of infertile spikelets. However, as increased grain number and fertile spikelet number per ear were significantly negatively correlated with thousand grain weight, sink strength limitations were evident. Genetic mapping identified 34 replicated quantitative trait loci (QTL) at two or more test environments, of which 24 resolved into eight loci each controlling two or more traits-termed here 'multi-trait QTL' (MT-QTL). These included MT-QTL associated with previously cloned genes controlling semi-dwarf plant stature, and with the genetic locus Reduced height 24 (Rht24) that further modulates plant height. Additionally, MT-QTL controlling spikelet number traits were located to chromosome 7A encompassing the gene WHEAT ORTHOLOG OF APO1 (WAPO-A1), and to its homoeologous location on chromosome 7B containing WAPO-B1. The genetic loci identified in this study, particularly those that potentially control multiple yield components, provide future opportunities for the targeted investigation of their underlying genes, gene networks and phenotypic trade-offs, in order to underpin further genetic gains in yield.

Comparative Analysis Based on Transcriptomics and Metabolomics Data Reveal Differences between Emmer and Durum Wheat in Response to Nitrogen Starvation

Mounting evidence indicates the key role of nitrogen (N) on diverse processes in plant, including development and defense. Using a combined transcriptomics and metabolomics approach, we studied the response of seedlings to N starvation of two different tetraploid wheat genotypes from the two main domesticated subspecies: emmer and durum wheat. We found that durum wheat exhibits broader and stronger response in comparison to emmer as seen from the expression pattern of both genes and metabolites and gene enrichment analysis. They showed major differences in the responses to N starvation for transcription factor families, emmer showed differential reduction in the levels of primary metabolites while durum wheat exhibited increased levels of most of them to N starvation. The correlation-based networks, including the differentially expressed genes and metabolites, revealed tighter regulation of metabolism in durum wheat in comparison to emmer. We also found that glutamate and γ-aminobutyric acid (GABA) had highest values of centrality in the metabolic correlation network, suggesting their critical role in the genotype-specific response to N starvation of emmer and durum wheat, respectively. Moreover, this finding indicates that there might be contrasting strategies associated to GABA and glutamate signaling modulating shoot vs. root growth in the two different wheat subspecies.

Identification and analysis of a differentially expressed wheat RING-type E3 ligase in spike primordia development during post-vernalization

We identified a RING-type E3 ligase (TaBAH1) protein in winter wheat that targets TaSAHH1 for degradation and might be involved in primordia development by regulating targeted protein degradation. Grain yield per spike in wheat (Triticum aestivum), is mainly determined prior to flowering during mature primordia development; however, the genes involved in primordia development have yet to be characterized. In this study, we demonstrated that, after vernalization for 50 days at 4 °C, there was a rapid acceleration in primordia development to the mature stages in the winter wheat cultivars Keumgang and Yeongkwang compared with the Chinese Spring cultivar. Although Yeongkwang flowers later than Keumgang under normal condition, it has the same heading time and reaches the WS9 stage of floral development after vernalization for 50 days. Using RNA sequencing, we identified candidate genes associated with primordia development in cvs. Keumgang and Yeongkwang, that are differentially expressed during wheat reproductive stages. Among these, the RING-type E3 ligase TaBAH1 (TraesCS5B01G373000) was transcriptionally upregulated between the double-ridge (WS2.5) stage and later stages of floret primordia development (WS10) after vernalization. Transient expression analysis indicated that TaBAH1 was localized to the plasma membrane and nucleus and was characterized by self-ubiquitination activity. Furthermore, we found that TaBAH1 interacts with TaSAHH1 to mediate its polyubiquitination and degradation through a 26S proteasomal pathway. Collectively, the findings of this study indicate that TaBAH1 might play a prominent role in post-vernalization floret primordia development.

Wheat omics: Classical breeding to new breeding technologies

Wheat is an important cereal crop, and its significance is more due to compete for dietary products in the world. Many constraints facing by the wheat crop due to environmental hazardous, biotic, abiotic stress and heavy matters factors, as a result, decrease the yield. Understanding the molecular mechanism related to these factors is significant to figure out genes regulate under specific conditions. Classical breeding using hybridization has been used to increase the yield but not prospered at the desired level. With the development of newly emerging technologies in biological sciences i.e., marker assisted breeding (MAB), QTLs mapping, mutation breeding, proteomics, metabolomics, next-generation sequencing (NGS), RNA_sequencing, transcriptomics, differential expression genes (DEGs), computational resources and genome editing techniques i.e. (CRISPR cas9; Cas13) advances in the field of omics. Application of new breeding technologies develops huge data; considerable development is needed in bioinformatics science to interpret the data. However, combined omics application to address physiological questions linked with genetics is still a challenge. Moreover, viroid discovery opens the new direction for research, economics, and target specification. Comparative genomics important to figure gene of interest processes are further discussed about considering the identification of genes, genomic loci, and biochemical pathways linked with stress resilience in wheat. Furthermore, this review extensively discussed the omics approaches and their effective use. Integrated plant omics technologies have been used viroid genomes associated with CRISPR and CRISPR-associated Cas13a proteins system used for engineering of viroid interference along with high-performance multidimensional phenotyping as a significant limiting factor for increasing stress resistance in wheat.

The soybean (Glycine max L.) cytokinin oxidase/dehydrogenase multigene family; Identification of natural variations for altered cytokinin content and seed yield

Cytokinins (CKs) play a fundamental role in regulating dynamics of organ source/sink relationships during plant development, including flowering and seed formation stages. As a result, CKs are key drivers of seed yield. The cytokinin oxidase/dehydrogenase (CKX) is one of the critical enzymes responsible for regulating plant CK levels by causing their irreversible degradation. Variation of CKX activity is significantly correlated with seed yield in many crop species while in soybean (Glycine max L.), the possible associations between CKX gene family members (GFMs) and yield parameters have not yet been assessed. In this study, 17 GmCKX GFMs were identified, and natural variations among GmCKX genes were probed among soybean cultivars with varying yield characteristics. The key CKX genes responsible for regulating CK content during seed filling stages of reproductive development were highlighted using comparative phylogenetics, gene expression analysis and CK metabolite profiling. Five of the seventeen identified GmCKX GFMs, showed natural variations in the form of single nucleotide polymorphisms (SNPs). The gene GmCKX7-1, with high expression during critical seed filling stages, was found to have a non-synonymous mutation (H105Q), on one of the active site residues, Histidine 105, previously reported to be essential for co-factor binding to maintain structural integrity of the enzyme. Soybean lines with this mutation had higher CK content and desired yield characteristics. The potential for marker-assisted selection based on the identified natural variation within GmCKX7-1, is discussed in the context of hormonal control that can result in higher soybean yield.

Stepwise increases in FT1 expression regulate seasonal progression of flowering in wheat (Triticum aestivum)

Flowering is regulated by genes that respond to changing daylengths and temperature, which have been well studied using controlled conditions; however, the molecular processes underpinning flowering in nature remain poorly understood. Here, we investigate the genetic pathways that coordinate flowering and inflorescence development of wheat (Triticum aestivum) as daylengths extend naturally in the field, using lines that contain variant alleles for the key photoperiod gene, Photoperiod-1 (Ppd-1). We found flowering involves a stepwise increase in the expression of FLOWERING LOCUS T1 (FT1), which initiates under day-neutral conditions of early spring. The incremental rise in FT1 expression is overridden in plants that contain a photoperiod-insensitive allele of Ppd-1, which hastens the completion of spikelet development and accelerates flowering time. The accelerated inflorescence development of photoperiod-insensitive lines is promoted by advanced seasonal expression of floral meristem identity genes. The completion of spikelet formation is promoted by FLOWERING LOCUS T2, which regulates spikelet number and is activated by Ppd-1. In wheat, flowering under natural photoperiods is regulated by stepwise increases in the expression of FT1, which responds dynamically to extending daylengths to promote early inflorescence development. This research provides a strong foundation to improve yield potential by fine-tuning the photoperiod-dependent control of inflorescence development.

High-Resolution Genome-wide Association Study Identifies Genomic Regions and Candidate Genes for Important Agronomic Traits in Wheat

Wheat (Triticum aestivum) is a major staple food crop worldwide. Genetic dissection of important agronomic traits is essential for continuous improvement of wheat yield to meet the demand of the world's growing population. We conducted a large-scale genome-wide association study (GWAS) using a panel of 768 wheat cultivars that were genotyped with 327 609 single-nucleotide polymorphisms generated by genotyping-by-sequencing and detected 395 quantitative trait loci (QTLs) for 12 traits under 7 environments. Among them, 273 QTLs were delimited to ≤1.0-Mb intervals and 7 of them are either known genes (Rht-D, Vrn-B1, and Vrn-D1) that have been cloned or known QTLs (TaGA2ox8, APO1, TaSus1-7B, and Rht12) that were previously mapped. Eight putative candidate genes were identified for three QTLs that enhance spike seed setting and grain size using gene expression data and were validated in three bi-parental populations. Protein sequence analysis identified 33 putative wheat orthologs that have high identity with rice genes in QTLs affecting similar traits. Large r² values for additive effects observed among the QTLs for most traits indicated that the phenotypes of these identified QTLs were highly predictable. Results from this study demonstrated that significantly increasing GWAS population size and marker density greatly improves detection and identification of candidate genes underlying a QTL, solidifying the foundation for large-scale QTL fine mapping, candidate gene validation, and developing functional markers for genomics-based breeding in wheat.

Tumorous Stem Development of Brassica Juncea: A Complex Regulatory Network of Stem Formation and Identification of Key Genes in Glucosinolate Biosynthesis

Stem mustard is a stem variety of mustard, an important Brassica vegetable. The formation and development of the tumorous stem, which is the key organ for the direct yield and quality, is a complex biological process involving morphogenesis, material accumulation and gene regulation. In this study, we demonstrated through anatomical studies that stem swelling is mainly dependent on the increase in the number of cells and the volume of parenchyma cells in the cortex and pith. To further understand transcript and metabolic changes during stem swelling, we obtained 27,901 differentially expressed genes, of which 671 were specifically detected using transcriptome sequencing technology in all four stages of stem swelling. Functional annotation identified enrichment for genes involved in photosynthesis, energy metabolism, cell growth, sulfur metabolism and glucosinolate biosynthesis. Glucosinolates are a group of nitrogen- and sulfur-containing secondary metabolites, which largely exist in the Cruciferous vegetables. HPLC analysis of the contents and components of glucosinolates in four different stem development stages revealed eight glucosinolates, namely, three aliphatic glucosinolates (sinigrin, glucoalyssin and gluconapin), four indole glucosinolates (4-hydroxyglucobrassicin, glucobrassicin, 4-methoxyglucobrassicin and neoglucobrassicin) and one aromatic glucosinolate (gluconasturtiin). All these types of glucosinolates showed a significant downward trend during the stem swelling period. The content of aliphatic glucosinolates was the highest, with sinigrin being the main component. In addition, qPCR was used to validate the expression of nine genes involved in glucosinolate biosynthesis. Most of these genes were down-regulated during stem swelling in qPCR, which is consistent with transcriptome data. These data provide a basic resource for further molecular and genetic research on Brassica juncea.

Heritable temporal gene expression patterns correlate with metabolomic seed content in developing hexaploid oat seed

Oat ranks sixth in world cereal production and has a higher content of health-promoting compounds compared with other cereals. However, there is neither a robust oat reference genome nor transcriptome. Using deeply sequenced full-length mRNA libraries of oat cultivar Ogle-C, a de novo high-quality and comprehensive oat seed transcriptome was assembled. With this reference transcriptome and QuantSeq 3' mRNA sequencing, gene expression was quantified during seed development from 22 diverse lines across six time points. Transcript expression showed higher correlations between adjacent time points. Based on differentially expressed genes, we identified 22 major temporal co-expression (TCoE) patterns of gene expression and revealed enriched gene ontology biological processes. Within each TCoE set, highly correlated transcripts, putatively commonly affected by genetic background, were clustered and termed genetic co-expression (GCoE) sets. Seventeen of the 22 TCoE sets had GCoE sets with median heritabilities higher than 0.50, and these heritability estimates were much higher than that estimated from permutation analysis, with no divergence observed in cluster sizes between permutation and non-permutation analyses. Linear regression between 634 metabolites from mature seeds and the PC1 score of each of the GCoE sets showed significantly lower p-values than permutation analysis. Temporal expression patterns of oat avenanthramides and lipid biosynthetic genes were concordant with previous studies of avenanthramide biosynthetic enzyme activity and lipid accumulation. This study expands our understanding of physiological processes that occur during oat seed maturation and provides plant breeders the means to change oat seed composition through targeted manipulation of key pathways.

Cytokinin dehydrogenase: a genetic target for yield improvement in wheat

The plant hormone group, the cytokinins, is implicated in both qualitative and quantitative components of yield. Cytokinins have opposing actions in shoot and root growth-actions shown to involve cytokinin dehydrogenase (CKX), the enzyme that inactivates cytokinin. We revise and provide unambiguous names for the CKX gene family members in wheat, based on the most recently released wheat genome database, IWGSC RefSeq v1.0 & v2.0. We review expression data of CKX gene family members in wheat, revealing tissue-specific gene family member expression as well as sub-genome-specific expression. Manipulation of CKX in cereals shows clear impacts on yield, root growth and orientation, and Zn nutrition, but this also emphasizes the necessity to unlink promotive effects on grain yield from negative effects of cytokinin on root growth and uptake of mineral nutrients, particularly Zn and Fe. Wheat is the most widely grown cereal crop globally, yet is under-research compared with rice and maize. We highlight gaps in our knowledge of the involvement of CKX for wheat. We also highlight the necessity for accurate analysis of endogenous cytokinins, acknowledging why this is challenging, and provide examples where inadequate analyses of endogenous cytokinins have led to unjustified conclusions. We acknowledge that the allohexaploid nature of bread wheat poses challenges in terms of uncovering useful mutations. However, we predict TILLING followed by whole-exome sequencing will uncover informative mutations and we indicate the potential for stacking mutations within the three genomes to modify yield components. We model a wheat ideotype based on CKX manipulation.

APETALA 2-like genes AP2L2 and Q specify lemma identity and axillary floral meristem development in wheat

The spikelet is the basic unit of the grass inflorescence. In tetraploid (Triticum turgidum) and hexaploid wheat (Triticum aestivum), the spikelet is a short indeterminate branch with two proximal sterile bracts (glumes) followed by a variable number of florets, each including a bract (lemma) with an axillary flower. Varying levels of miR172 and/or its target gene Q (AP2L5) result in gradual transitions of glumes to lemmas, and vice versa. Here, we show that AP2L5 and its related paralog AP2L2 play critical and redundant roles in the specification of axillary floral meristems and lemma identity. AP2L2, also targeted by miR172, displayed similar expression profiles to AP2L5 during spikelet development. Loss-of-function mutants in both homeologs of AP2L2 (henceforth ap2l2) developed normal spikelets, but ap2l2 ap2l5 double mutants generated spikelets with multiple empty bracts before transitioning to florets. The coordinated nature of these changes suggest an early role of these genes in floret development. Moreover, the flowers of ap2l2 ap2l5 mutants showed organ defects in paleas and lodicules, including the homeotic conversion of lodicules into carpels. Mutations in the miR172 target site of AP2L2 were associated with reduced plant height, more compact spikes, promotion of lemma-like characters in glumes and smaller lodicules. Taken together, our results show that the balance in the expression of miR172 and AP2-like genes is crucial for the correct development of spikelets and florets, and that this balance has been altered during the process of wheat and barley (Hordeum vulgare) domestication. The manipulation of this regulatory module provides an opportunity to modify spikelet architecture and improve grain yield.

A set of AP2-like genes is associated with inflorescence branching and architecture in domesticated rice

Rice yield is influenced by inflorescence size and architecture, and inflorescences from domesticated rice accessions produce more branches and grains. Neither the molecular control of branching nor the developmental differences between wild and domesticated rice accessions are fully understood. We surveyed phenotypes related to branching, size, and grain yield across 91 wild and domesticated African and Asian accessions. Characteristics related to axillary meristem identity were the main phenotypic differences between inflorescences from wild and domesticated accessions. We used whole transcriptome sequencing in developing inflorescences to measure gene expression before and after the transition from branching axillary meristems to determinate spikelet meristems. We identified a core set of genes associated with axillary meristem identity in Asian and African rice, and another set associated with phenotypic variability between wild and domesticated accessions. AP2/EREBP-like genes were enriched in both sets, suggesting that they are key factors in inflorescence branching and rice domestication. Our work has identified new candidates in the molecular control of inflorescence development and grain yield, and provides a detailed description of the effects of domestication on phenotype and gene expression.