FRIZZY PANICLE defines a regulatory hub for simultaneously controlling spikelet formation and awn elongation in bread wheat

The OsNAC41-RoLe1-OsAGAP module promotes root development and drought resistance in upland rice

Drought is a major environmental stress limiting crop yields worldwide. Upland rice (Oryza sativa) has evolved complex genetic mechanisms for adaptative growth under drought stress. However, few genetic variants that mediate drought resistance in upland rice have been identified, and little is known about the evolution of this trait during rice domestication. In this study, using a genome-wide association study we identified ROOT LENGTH 1 (RoLe1) that controls rice root length and drought resistance. We found that a G-to-T polymorphism in the RoLe1 promoter causes increased binding of the transcription factor OsNAC41 and thereby enhanced expression of RoLe1. We further showed that RoLe1 interacts with OsAGAP, an ARF-GTPase activating protein involved in auxin-dependent root development, and interferes with its function to modulate root development. Interestingly, RoLe1 could enhance crop yield by increasing the seed-setting rate under moderate drought conditions. Genomic evolutionary analysis revealed that a newly arisen favorable allelic variant, proRoLe1-526T, originated from the midwest Asia and was retained in upland rice during domestication. Collectively, our study identifies an OsNAC41-RoLe1-OsAGAP module that promotes upland rice root development and drought resistance, providing promising genetic targets for molecular breeding of drought-resistant rice varieties.

The TaWAK2-TaNAL1-TaDST pathway regulates leaf width via cytokinin signaling in wheat

Leaves play a crucial role in photosynthesis and respiration, ultimately affecting the final grain yield of crops, including wheat (Triticum aestivum L.); however, the molecular mechanisms underlying wheat leaf development remain largely unknown. Here, we isolated a narrow-leaf gene, TaWAK2-A, through a map-based cloning strategy. TaWAK2-A encodes a wall-associated kinase (WAK), for which a single Ala-to-Val amino acid substitution reduces the protein stability, leading to a narrow-leaf phenotype in wheat. Further investigation suggests that TaWAK2 directly interacts with and phosphorylates TaNAL1, a trypsin-like serine/cysteine protease. The phosphorylated TaNAL1 is then involved in the degradation of the zinc finger transcription factor TaDST, which acts as a repressor of leaf expansion by activating the expression of the cytokinin oxidase gene TaCKX9 and triggering in vivo cytokinin degradation. Therefore, our findings elucidate a signaling cascade involving TaWAK2-TaNAL1-TaDST that sheds light on the regulation of wheat leaf development.

Genomic insights into the modifications of spike morphology traits during wheat breeding

Over the past century, environmental changes have significantly impacted wheat spike morphology, crucial for adaptation and grain yield. However, the changes in wheat spike modifications during this period remain largely unknown. This study examines 16 spike morphology traits in 830 accessions released from 1900 to 2020. It finds that spike weight, grain number per spike (GN), and thousand kernel weight have significantly increased, while spike length has no significant change. The increase in fertile spikelets is due to fewer degenerated spikelets, resulting in a higher GN. Genome-wide association studies identified 49,994 significant SNPs, grouped into 293 genomic regions. The accumulation of favorable alleles in these genomic regions indicates the genetic basis for modification in spike morphology traits. Genetic network analysis of these genomic regions reveals the genetic basis for phenotypic correlations among spike morphology traits. The haplotypes of the identified genomic regions display obvious geographical differentiation in global accessions and environmental adaptation over the past 120 years. In summary, we reveal the genetic basis of adaptive evolution and the interactions of spike morphology, offering valuable resources for the genetic improvement of spike morphology to enhance environmental adaptation.

Functional Characterization of Accessible Chromatin in Common Wheat

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

Mapping QTLs with additive and epistatic effects for awn length and their effects on kernel-related traits in common wheat

Introduction

Wheat awns are crucial determinants of wheat yield due to their capacity to photosynthesize and exchange gas. Understanding the genetic basis of awn length (AL) is essential for improving wheat yield in molecular breeding programs.

Methods

In this study, quantitative trait loci (QTLs) of AL were analyzed using recombinant inbred line (RIL) mapping population referred to as YY-RILs, which was derived from a cross between Yannong 15 (YN15) and Yannong 1212 (YN1212).

Results and discussion

Seven putative additive QTLs and 30 pairwise epistatic QTLs for AL were identified. Among them, five novel additive QTLs (except qAl-2A and qAl-5A.2) and 30 novel pairwise epistatic QTLs were identified. qAl-5A.1 was repeatedly identified in all five environment datasets, which was considered to be one novel stable QTL for AL with minor additive effects. eqAl-2B.2-2 significantly interacted with eight loci and could be of great importance in regulating awn development. The genes associated with the major stable QTL of qAl-5A.2 and the minor stable QTL of qAl-2A were B1 and WFZP-A, respectively. Awn lengths exhibited significant genetic correlations with kernel weight and kernels per spike, which could affect grain protein content to a lesser extent. This study enhances our understanding of the genetic basis of awn development and identifies novel genes as well as markers for future genetic improvement of wheat yield.

Identification, validation and candidate gene analysis of major QTL for Supernumerary spikelets in wheat

BACKGROUND

The number of spikelets per spike is a key trait that affects the yield of bread wheat (Triticum aestivum L.). Identification of the QTL for spikelets per spike and its genetic effects that could be used in molecular assistant breeding in the future.

RESULTS

In this study, four recombinant inbred line (RIL) populations were generated and used, having YuPi branching wheat (YP), with Supernumerary Spikelets (SS) phenotype, as a common parent. QTL (QSS.sicau-2 A and QSS.sicau-2D) related to SS trait were mapped on chromosomes 2 A and 2D through bulked segregant exome sequencing (BSE-Seq). Fourteen molecular markers were further developed within the localization interval, and QSS.sicau-2 A was narrowed to 3.0 cM covering 7.6 Mb physical region of the reference genome, explaining 13.7 - 15.9% the phenotypic variance. Similarly, the QSS.sicau-2D was narrowed to 1.8 cM covering 2.4 Mb physical region of the reference genome, and it explained 27.4 - 32.9% the phenotypic variance. These two QTL were validated in three different genetic backgrounds using the linked markers. QSS.sicau-2 A was identified as WFZP-A, and QSS.sicau-2D was identified a novel locus, different to the previously identified WFZP-D. Based on the gene expression patterns, gene annotation and sequence analysis, TraesCS2D03G0260700 was predicted to be a potential candidate gene for QSS.sicau-2D.

CONCLUSION

Two significant QTL for SS, namely QSS.sicau-2 A and QSS.sicau-2D were identified in multiple environments were identified and their effect in diverse genetic populations was assessed. QSS.sicau-2D is a novel QTL associated with the SS trait, with TraesCS2D03G0260700 predicted as its candidate gene.

Identification and validation of two major QTL for grain number per spike on chromosomes 2B and 2D in bread wheat (Triticum aestivum L.)

KEY MESSAGE

Major QTL for grain number per spike were identified on chromosomes 2B and 2D. Haplotypes and candidate genes of QGns.cib-2B.1 were analyzed. Grain number per spike (GNS) is one of the main components of wheat yield. Genetic dissection of their regulatory factors is essential to improve the yield potential. In present study, a recombinant inbred line population comprising 180 lines developed from the cross between a high GNS line W7268 and a cultivar Chuanyu12 was employed to identify quantitative trait loci (QTL) associated with GNS across six environments. Two major QTL, QGns.cib-2B.1 and QGns.cib-2D.1, were detected in at least four environments with the phenotypic variations of 12.99-27.07% and 8.50-13.79%, respectively. And significant interactions were observed between the two major QTL. In addition, QGns.cib-2B.1 is a QTL cluster for GNS, grain number per spikelet and fertile tiller number, and they were validated in different genetic backgrounds using Kompetitive Allele Specific PCR (KASP) markers. QGns.cib-2B.1 showed pleotropic effects on other yield-related traits including plant height, spike length, and spikelet number per spike, but did not significantly affect thousand grain weight which suggested that it might be potentially applicable in breeding program. Comparison analysis suggested that QGns.cib-2B.1 might be a novel QTL. Furthermore, haplotype analysis of QGns.cib-2B.1 indicated that it is a hot spot of artificial selection during wheat improvement. Based on the expression patterns, gene annotation, orthologs analysis and sequence variations, the candidate genes of QGns.cib-2B.1 were predicted. Collectively, the major QTL and KASP markers reported here provided a wealth of information for the genetic basis of GNS and grain yield improvement.

Dissecting the molecular basis of spike traits by integrating gene regulatory networks and genetic variation in wheat

Spike architecture influences both grain weight and grain number per spike, which are the two major components of grain yield in bread wheat (Triticum aestivum L.). However, the complex wheat genome and the influence of various environmental factors pose challenges in mapping the causal genes that affect spike traits. Here, we systematically identified genes involved in spike trait formation by integrating information on genomic variation and gene regulatory networks controlling young spike development in wheat. We identified 170 loci that are responsible for variations in spike length, spikelet number per spike, and grain number per spike through genome-wide association study and meta-QTL analyses. We constructed gene regulatory networks for young inflorescences at the double ridge stage and the floret primordium stage, in which the spikelet meristem and the floret meristem are predominant, respectively, by integrating transcriptome, histone modification, chromatin accessibility, eQTL, and protein-protein interactome data. From these networks, we identified 169 hub genes located in 76 of the 170 QTL regions whose polymorphisms are significantly associated with variation in spike traits. The functions of TaZF-B1, VRT-B2, and TaSPL15-A/D in establishment of wheat spike architecture were verified. This study provides valuable molecular resources for understanding spike traits and demonstrates that combining genetic analysis and developmental regulatory networks is a robust approach for dissection of complex traits.

Systematic identification of wheat spike developmental regulators by integrated multi-omics, transcriptional network, GWAS, and genetic analyses

The spike architecture of wheat plays a crucial role in determining grain number, making it a key trait for optimization in wheat breeding programs. In this study, we used a multi-omic approach to analyze the transcriptome and epigenome profiles of the young spike at eight developmental stages, revealing coordinated changes in chromatin accessibility and H3K27me3 abundance during the flowering transition. We constructed a core transcriptional regulatory network (TRN) that drives wheat spike formation and experimentally validated a multi-layer regulatory module involving TaSPL15, TaAGLG1, and TaFUL2. By integrating the TRN with genome-wide association studies, we identified 227 transcription factors, including 42 with known functions and 185 with unknown functions. Further investigation of 61 novel transcription factors using multiple homozygous mutant lines revealed 36 transcription factors that regulate spike architecture or flowering time, such as TaMYC2-A1, TaMYB30-A1, and TaWRKY37-A1. Of particular interest, TaMYB30-A1, downstream of and repressed by WFZP, was found to regulate fertile spikelet number. Notably, the excellent haplotype of TaMYB30-A1, which contains a C allele at the WFZP binding site, was enriched during wheat breeding improvement in China, leading to improved agronomic traits. Finally, we constructed a free and open access Wheat Spike Multi-Omic Database (http://39.98.48.156:8800/#/). Our study identifies novel and high-confidence regulators and offers an effective strategy for dissecting the genetic basis of wheat spike development, with practical value for wheat breeding.

Genetic dissection of major QTL for grain number per spike on chromosomes 5A and 6A in bread wheat (Triticum aestivum L.)

Grain number per spike (GNS) is a crucial component of grain yield and plays a significant role in improving wheat yield. To identify quantitative trait loci (QTL) associated with GNS, a recombinant inbred line (RIL) population derived from the cross of Zhongkemai 13F10 and Chuanmai 42 was employed to conduct QTL mapping across eight environments. Based on the bulked segregant exome sequencing (BSE-Seq), genomic regions associated with GNS were detected on chromosomes 5A and 6A. According to the constructed genetic maps, two major QTL QGns.cib-5A (LOD = 4.35-8.16, PVE = 8.46-14.43%) and QGns.cib-6A (LOD = 3.82-30.80, PVE = 5.44-12.38%) were detected in five and four environments, respectively. QGns.cib-6A is a QTL cluster for other seven yield-related traits. QGns.cib-5A and QGns.cib-6A were further validated using linked Kompetitive Allele Specific PCR (KASP) markers in different genetic backgrounds. QGns.cib-5A exhibited pleiotropic effects on productive tiller number (PTN), spike length (SL), fertile spikelet number per spike (FSN), and ratio of grain length to grain width (GL/GW) but did not significantly affect thousand grain weight (TGW). Haplotype analysis revealed that QGns.cib-5A and QGns.cib-6A were the targets of artificial selection during wheat improvement. Candidate genes for QGns.cib-5A and QGns.cib-6A were predicted by analyzing gene annotation, spatiotemporal expression patterns, and orthologous and sequence differences. These findings will be valuable for fine mapping and map-based cloning of genes underlying QGns.cib-5A and QGns.cib-6A.

MADS1-regulated lemma and awn development benefits barley yield

Floral organ shape and size in cereal crops can affect grain size and yield, so genes that regulate their development are promising breeding targets. The lemma, which protects inner floral organs, can physically constrain grain growth; while the awn, a needle-like extension of the lemma, creates photosynthate to developing grain. Although several genes and modules controlling grain size and awn/lemma growth in rice have been characterized, these processes, and the relationships between them, are not well understood for barley and wheat. Here, we demonstrate that the barley E-class gene HvMADS1 positively regulates awn length and lemma width, affecting grain size and weight. Cytological data indicates that HvMADS1 promotes awn and lemma growth by promoting cell proliferation, while multi-omics data reveals that HvMADS1 target genes are associated with cell cycle, phytohormone signaling, and developmental processes. We define two potential targets of HvMADS1 regulation, HvSHI and HvDL, whose knockout mutants mimic awn and/or lemma phenotypes of mads1 mutants. Additionally, we demonstrate that HvMADS1 interacts with APETALA2 (A-class) to synergistically activate downstream genes in awn/lemma development in barley. Notably, we find that MADS1 function remains conserved in wheat, promoting cell proliferation to increase awn length. These findings extend our understanding of MADS1 function in floral organ development and provide insights for Triticeae crop improvement strategies.

Boosting Triticeae crop grain yield by manipulating molecular modules to regulate inflorescence architecture: insights and knowledge from other cereal crops

One of the challenges for global food security is to reliably and sustainably improve the grain yield of cereal crops. One solution is to modify the architecture of the grain-bearing inflorescence to optimize for grain number and size. Cereal inflorescences are complex structures, with determinacy, branching patterns, and spikelet/floret growth patterns that vary by species. Recent decades have witnessed rapid advancements in our understanding of the genetic regulation of inflorescence architecture in rice, maize, wheat, and barley. Here, we summarize current knowledge on key genetic factors underlying the different inflorescence morphologies of these crops and model plants (Arabidopsis and tomato), focusing particularly on the regulation of inflorescence meristem determinacy and spikelet meristem identity and determinacy. We also discuss strategies to identify and utilize these superior alleles to optimize inflorescence architecture and, ultimately, improve crop grain yield.

Whole-genome sequencing of diverse wheat accessions uncovers genetic changes during modern breeding in China and the United States

Breeding has dramatically changed the plant architecture of wheat (Triticum aestivum), resulting in the development of high-yielding varieties adapted to modern farming systems. However, how wheat breeding shaped the genomic architecture of this crop remains poorly understood. Here, we performed a comprehensive comparative analysis of a whole-genome resequencing panel of 355 common wheat accessions (representing diverse landraces and modern cultivars from China and the United States) at the phenotypic and genomic levels. The genetic diversity of modern wheat cultivars was clearly reduced compared to landraces. Consistent with these genetic changes, most phenotypes of cultivars from China and the United States were significantly altered. Of the 21 agronomic traits investigated, 8 showed convergent changes between the 2 countries. Moreover, of the 207 loci associated with these 21 traits, more than half overlapped with genomic regions that showed evidence of selection. The distribution of selected loci between the Chinese and American cultivars suggests that breeding for increased productivity in these 2 regions was accomplished by pyramiding both shared and region-specific variants. This work provides a framework to understand the genetic architecture of the adaptation of wheat to diverse agricultural production environments, as well as guidelines for optimizing breeding strategies to design better wheat varieties.

CRISPR-mediated acceleration of wheat improvement: advances and perspectives

Common wheat (Triticum aestivum) is one of the most widely cultivated and consumed crops globally. In the face of limited arable land and climate changes, it is a great challenge to maintain current and increase future wheat production. Enhancing agronomic traits in wheat by introducing mutations across all three homoeologous copies of each gene has proven to be a difficult task due to its large genome with high repetition. However, clustered regularly interspaced short palindromic repeat (CRISPR)/CRISPR-associated nuclease (Cas) genome editing technologies offer a powerful means of precisely manipulating the genomes of crop species, thereby opening up new possibilities for biotechnology and breeding. In this review, we first focus on the development and optimization of the current CRISPR-based genome editing tools in wheat, emphasizing recent breakthroughs in precise and multiplex genome editing. We then describe the general procedure of wheat genome editing and highlight different methods to deliver the genome editing reagents into wheat cells. Furthermore, we summarize the recent applications and advancements of CRISPR/Cas technologies for wheat improvement. Lastly, we discuss the remaining challenges specific to wheat genome editing and its future prospects.

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.

Characterization and fine mapping analysis of a major stable QTL qKnps-4A for kernel number per spike in wheat

KEY MESSAGE

A major stable QTL for kernel number per spike was narrowed down to a 2.19-Mb region containing two potential candidate genes, and its effects on yield-related traits were characterized. Kernel number per spike (KNPS) in wheat is a key yield component. Dissection and characterization of major stable quantitative trait loci (QTLs) for KNPS would be of considerable value for the genetic improvement of yield potential using molecular breeding technology. We had previously reported a major stable QTL controlling KNPS, qKnps-4A. In the current study, primary fine-mapping analysis, based on the primary mapping population, located qKnps-4A to an interval of approximately 6.8-Mb from 649.0 to 655.8 Mb on chromosome 4A refering to 'Kenong 9204' genome. Further fine-mapping analysis based on a secondary mapping population narrowed qKnps-4A to an approximately 2.19-Mb interval from 653.72 to 655.91 Mb. Transcriptome sequencing, gene function annotation analysis and homologous gene related reports showed that TraesKN4A01HG38570 and TraesKN4A01HG38590 were most likely to be candidate genes of qKnps-4A. Phenotypic analysis based on paired near-isogenic lines in the target region showed that qKnps-4A increased KNPS mainly by increasing the number of central florets per spike. We also evaluated the effects of qKnps-4A on other yield-related traits. Moreover, we dissected the QTL cluster of qKnps-4A and qTkw-4A and proved that the phenotypic effects were probably due to close linkage of two or more genes rather than pleiotropic effects of a single gene. This study provides molecular marker resource for wheat molecular breeding designed to improve yield potential, and lay the foundation for gene functional analysis of qKnps-4A.

A promising QTL QSns.sau-MC-3D.1 likely superior to WAPO1 for the number of spikelets per spike of wheat shows no adverse effects on yield-related traits

KEY MESSAGE

A likely new locus QSns.sau-MC-3D.1 associated with SNS showing no negative effect on yield-related traits compared to WAPO1 was identified and validated in various genetic populations under multiple environments. The number of spikelets per spike (SNS) is one of the crucial factors determining wheat yield. Thus, improving our understanding of the genes that regulate SNS could help develop wheat varieties with higher yield. In this study, a recombinant inbred line (RIL) population (MC) containing 198 lines derived from a cross between msf and Chuannong 16 (CN16) was used to construct a genetic linkage map using the GenoBaits Wheat 16 K Panel. The genetic map contained 5,991 polymorphic SNP markers spanning 2,813.25 cM. A total of twelve QTL for SNS were detected, and two of them, i.e., QSns.sau-MC-3D.1 and QSns.sau-MC-7A, were stably expressed. QSns.sau-MC-3D.1 had high LOD values ranging from 4.99 to 11.06 and explained 9.71-16.75% of the phenotypic variation. Comparison of QSns.sau-MC-3D.1 with previously reported SNS QTL suggested that it is likely a novel one, and two kompetitive allele-specific PCR (KASP) markers were further developed. The positive effect of QSns.sau-MC-3D.1 was also validated in three biparental populations and a diverse panel containing 388 Chinese wheat accessions. Genetic analysis indicated that WHEAT ORTHOLOG OFAPO1 (WAPO1) was a candidate gene for QSns.sau-MC-7A. Pyramiding of QSns.sau-MC-3D.1 and WAP01 had a great additive effect increasing SNS by 7.10%. Correlation analysis suggested that QSns.sau-MC-3D.1 was likely independent of effective tiller number, plant height, spike length, anthesis date, and thousand kernel weight. However, the H2 haplotype of WAPO1 may affect effective tiller number and plant height. These results indicated that utilization of QSns.sau-MC-3D.1 should be given priority for wheat breeding. Geographical distribution analysis showed that the positive allele of QSns.nsau-MC-3D.1 was dominant in most wheat-producing regions of China, and it has been positively selected among modern cultivars released in China since the 1940s. Gene prediction, qRT-PCR analysis, and sequence alignment suggested that TraesCS3D03G0216800 may be the candidate gene of QSns.nsau-MC-3D.1. Taken together, these results enrich our understanding of the genetic basis of wheat SNS and will be useful for fine mapping and cloning of the gene underlying QSns.sau-MC-3D.1.

The TaTCP4/10-B1 cascade regulates awn elongation in wheat (Triticum aestivum L.)

Awns are important morphological markers for wheat and exert a strong physiological effect on wheat yield. The awn elongation suppressor B1 has recently been cloned through association and linkage analysis in wheat. However, the mechanism of awn inhibition centered around B1 remains to be clarified. Here, we identified an allelic variant in the coding region of B1 through analysis of re-sequencing data; this variant causes an amino acid substitution and premature termination, resulting in a long-awn phenotype. Transcriptome analysis indicated that B1 inhibited awn elongation by impeding cytokinin- and auxin-promoted cell division. Moreover, B1 directly repressed the expression of TaRAE2 and TaLks2, whose orthologs have been reported to promote awn development in rice or barley. More importantly, we found that TaTCP4 and TaTCP10 synergistically inhibited the expression of B1, and a G-to-A mutation in the B1 promoter attenuated its inhibition by TaTCP4/10. Taken together, our results reveal novel mechanisms of awn development and provide genetic resources for trait improvement in wheat.

Genetic architecture and molecular regulation of sorghum domestication

Over time, wild crops have been domesticated by humans, and the knowledge gained from parallel selection and convergent domestication-related studies in cereals has contributed to current techniques used in molecular plant breeding. Sorghum (Sorghum bicolor (L.) Moench) is the world's fifth-most popular cereal crop and was one of the first crops cultivated by ancient farmers. In recent years, genetic and genomic studies have provided a better understanding of sorghum domestication and improvements. Here, we discuss the origin, diversification, and domestication processes of sorghum based on archeological discoveries and genomic analyses. This review also comprehensively summarized the genetic basis of key genes related to sorghum domestication and outlined their molecular mechanisms. It highlights that the absence of a domestication bottleneck in sorghum is the result of both evolution and human selection. Additionally, understanding beneficial alleles and their molecular interactions will allow us to quickly design new varieties by further de novo domestication.

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.

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

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

Conservation and divergence: Regulatory networks underlying reproductive branching in rice and maize

BACKGROUND

Cereal crops are a major source of raw food and nutrition for humans worldwide. Inflorescence of cereal crops is their reproductive organ, which also contributes to crop productivity. The branching pattern in flowering plant species not only determines inflorescence architecture but also determines the grain yield. There are good reviews describing the grass inflorescence architecture contributing to the final grain yield. However, very few discuss the aspects of inflorescence branching.

AIM OF REVIEW

This review aimed at systematically and comprehensively summarizing the latest progress in the field of conservation and divergence of genetic regulatory network that controls inflorescence branching in maize and rice, provide strategies to efficiently utilize the achievements in reproductive branching for crop yield improvement, and suggest a potential regulatory network underlying the inflorescence branching and vegetative branching system.

KEY SCIENTIFIC CONCEPTS OF REVIEW

Inflorescence branching is the consequence of a series of developmental events including the initiation, outgrowth, determinacy, and identity of reproductive axillary meristems, and it is controlled by a complex functional hierarchy of genetic networks. Initially, we compared the inflorescence architecture of maize and rice; then, we reviewed the genetic regulatory pathways controlling the inflorescence meristem size, bud initiation, and outgrowth, and the key transition steps that shape the inflorescence branching in maize and rice; additionally, we summarized strategies to effectively apply the recent advances in inflorescence branching for crop yield improvement. Finally, we discussed how the newly discovered hormones coordinate the regulation of inflorescence branching and yield traits. Furthermore, we discussed the possible reason behind distinct regulatory pathways for vegetative and inflorescence branching.

Transcriptional signatures of wheat inflorescence development

In order to maintain global food security, it will be necessary to increase yields of the cereal crops that provide most of the calories and protein for the world's population, which includes common wheat (Triticum aestivum L.). An important wheat yield component is the number of grain-holding spikelets which form on the spike during inflorescence development. Characterizing the gene regulatory networks controlling the timing and rate of inflorescence development will facilitate the selection of natural and induced gene variants that contribute to increased spikelet number and yield. In the current study, co-expression and gene regulatory networks were assembled from a temporal wheat spike transcriptome dataset, revealing the dynamic expression profiles associated with the progression from vegetative meristem to terminal spikelet formation. Consensus co-expression networks revealed enrichment of several transcription factor families at specific developmental stages including the sequential activation of different classes of MIKC-MADS box genes. This gene regulatory network highlighted interactions among a small number of regulatory hub genes active during terminal spikelet formation. Finally, the CLAVATA and WUSCHEL gene families were investigated, revealing potential roles for TtCLE13, TtWOX2, and TtWOX7 in wheat meristem development. The hypotheses generated from these datasets and networks further our understanding of wheat inflorescence development.

Revisiting the origin and identity specification of the spikelet: A structural innovation in grasses (Poaceae)

Spikelets are highly specialized and short-lived branches and function as a constitutional unit of the complex grass inflorescences. A series of genetic, genomic, and developmental studies across different clades of the family have called for and permitted a synthesis on the regulation and evolution of spikelets, and hence inflorescence diversity. Here, we have revisited the identity specification of a spikelet, focusing on the diagnostic features of a spikelet from morphological, developmental, and molecular perspectives. Particularly, recent studies on a collection of barley (Hordeum vulgare L.), wheat (Triticum spp.), and rice (Oryza sativa L.) mutants have highlighted a set of transcription factors that are important in the control of spikelet identity and the patterning of floral parts of a spikelet. In addition, we have endeavored to clarify some puzzling issues on the (in)determinacy and modifications of spikelets over the course of evolution. Meanwhile, genomes of two sister taxa of the remaining grass species have again demonstrated the importance of genome duplication and subsequent gene losses on the evolution of spikelets. Accordingly, we argue that changes in the orthologs of spikelet-related genes could be critical for the development and evolution of the spikelet, an evolutionary innovation in the grass family. Likewise, the conceptual discussions on the regulation of a fundamental unit of compound inflorescences could be translated into other organismal groups where compound structures are similarly formed, permitting a comparative perspective on the control of biological complexity.

Molecular, cellular, and developmental foundations of grass diversity

Humans have cultivated grasses for food, feed, beverages, and construction materials for millennia. Grasses also dominate the landscape in vast parts of the world, where they have adapted morphologically and physiologically, diversifying to form ~12,000 species. Sequences of hundreds of grass genomes show that they are essentially collinear; nonetheless, not all species have the same complement of genes. Here, we focus on the molecular, cellular, and developmental bases of grain yield and dispersal-traits that are essential for domestication. Distinct genes, networks, and pathways were selected in different crop species, reflecting underlying genomic diversity. With increasing genomic resources becoming available in nondomesticated species, we anticipate advances in coming years that illuminate the ecological and economic success of the grasses.

The miR319/TaGAMYB3 module regulates plant architecture and improves grain yield in common wheat (Triticum aestivum)

Plant architecture is a key determinant of crop productivity and adaptation. The highly conserved microRNA319 (miR319) family functions in various biological processes, but little is known about how miR319 regulates plant architecture in wheat (Triticum aestivum). Here, we determined that the miR319/TaGAMYB3 module controls plant architecture and grain yield in common wheat. Repressing tae-miR319 using short tandem target mimics resulted in favorable plant architecture traits, including increased plant height, reduced tiller number, enlarged spikes and flag leaves, and thicker culms, as well as enhanced grain yield in field plot tests. Overexpressing tae-miR319 had the opposite effects on plant architecture and grain yield. Although both TaPCF8 and TaGAMYB3 were identified as miR319 target genes, genetic complementation assays demonstrated that only miR319-resistant TaGAMYB3 (rTaGAMYB3) abolished tae-miR319-mediated growth inhibition of flag leaves and spikes. TaGAMYB3 functions as a transcriptional activator of downstream genes, including TaPSKR1, TaXTH23, TaMADS5 and TaMADS51, by binding to their promoters. Furthermore, TaGAMYB3 physically interacts with TaBA1, an important regulator of spike development, to additively activate the transcription of downstream genes such as TaMADS5. Our findings provide insight into how the miR319/TaGAMYB3 module regulates plant architecture and improves grain yield in common wheat.

Improving bread wheat yield through modulating an unselected AP2/ERF gene

Crop breeding heavily relies on natural genetic variation. However, additional new variations are desired to meet the increasing human demand. Inflorescence architecture determines grain number per spike, a major determinant of bread wheat (Triticum aestivum L.) yield. Here, using Brachypodium distachyon as a wheat proxy, we identified DUO-B1, encoding an APETALA2/ethylene response factor (AP2/ERF) transcription factor, regulating spike inflorescence architecture in bread wheat. Mutations of DUO-B1 lead to mild supernumerary spikelets, increased grain number per spike and, importantly, increased yield under field conditions without affecting other major agronomic traits. DUO-B1 suppresses cell division and promotes the expression of BH^t/WFZP, whose mutations could lead to branched 'miracle-wheat'. Pan-genome analysis indicated that DUO-B1 has not been utilized in breeding, and holds promise to increase wheat yield further.

Genetic and transcriptomic dissection of an artificially induced paired spikelets mutant of wheat (Triticum aestivum L.)

Morphological, genetic and transcriptomic characterizations of an EMS-induced wheat paired spikelets (PS) mutant were performed. A novel qualitative locus WPS1 on chromosome 1D was identified. Grain yield of wheat is significantly associated with inflorescence or spike architecture. However, few genes related to wheat spike development have been identified and their underlying mechanisms are largely unknown. In this study, we characterized an ethyl methanesulfonate (EMS)-induced wheat mutant, wheat paired spikelets 1 (wps1). Unlike a single spikelet that usually develops at each node of rachis, a secondary spikelet appeared below the primary spikelet at most of the rachis nodes of wps1. The microscope observation showed that the secondary spikelet initiated later than the primary spikelet. Genetic analysis suggested that the PS of wps1 is controlled by a single dominant nuclear gene, designated WHEAT PAIRED SPIKELETS 1 (WPS1). Further RNA-seq based bulked segregant analysis and molecular marker mapping localized WPS1 in an interval of 208.18-220.92 Mb on the chromosome arm 1DL, which is different to known genes related to spike development in wheat. By using wheat omics data, TraesCS1D02G155200 encoding a HD-ZIP III transcription factor was considered as a strong candidate gene for WPS1. Transcriptomic analysis indicated that PS formation in wps1 is associated with auxin-related pathways and may be regulated by networks involving TB1, Ppd1, FT1, VRN1, etc. This study laid the solid foundation for further validation of the causal gene of WPS1 and explored its regulatory mechanism in PS formation and inflorescence development, which may benefit to kernel yield improvement of wheat based on optimization or design of spike architecture in the future.

The control of compound inflorescences: insights from grasses and legumes

A major challenge in biology is to understand how organisms have increased developmental complexity during evolution. Inflorescences, with remarkable variation in branching systems, are a fitting model to understand architectural complexity. Inflorescences bear flowers that may become fruits and/or seeds, impacting crop productivity and species fitness. Great advances have been achieved in understanding the regulation of complex inflorescences, particularly in economically and ecologically important grasses and legumes. Surprisingly, a synthesis is still lacking regarding the common or distinct principles underlying the regulation of inflorescence complexity. Here, we synthesize the similarities and differences in the regulation of compound inflorescences in grasses and legumes, and propose that the emergence of novel higher-order repetitive modules is key to the evolution of inflorescence complexity.

Translational repression of FZP mediated by CU-rich element/OsPTB interactions modulates panicle development in rice

Panicle development is an important determinant of the grain number in rice. A thorough characterization of the molecular mechanism underlying panicle development will lead to improved breeding of high-yielding rice varieties. Frizzy Panicle (FZP), a critical gene for panicle development, is regulated by OsBZR1 and OsARFs at the transcriptional stage. However, the translational modulation of FZP has not been reported. We reveal that the CU-rich elements (CUREs) in the 3' UTR of the FZP mRNA are crucial for efficient FZP translation. The knockout of CUREs in the FZP 3' UTR or the over-expression of the FZP 3' UTR fragment containing CUREs resulted in an increase in FZP mRNA translation efficiency. Moreover, the number of secondary branches (NSB) and the grain number per panicle (GNP) decreased in the transformed rice plants. The CUREs in the 3' UTR of FZP mRNA were verified as the targets of the polypyrimidine tract-binding proteins OsPTB1 and OsPTB2 in rice. Both OsPTB1 and OsPTB2 were highly expressed in young panicles. The knockout of OsPTB1/2 resulted in an increase in the FZP translational efficiency and a decrease in the NSB and GNP. Furthermore, the over-expression of OsPTB1/2 decreased the translation of the reporter gene fused to FZP 3' UTR in vivo and in vitro. These results suggest that OsPTB1/2 can mediate FZP translational repression by interacting with CUREs in the 3' UTR of FZP mRNA, leading to changes in the NSB and GNP. Accordingly, in addition to transcriptional regulation, FZP expression is also fine-tuned at the translational stage during rice panicle development.

Genetic Incorporation of the Favorable Alleles for Three Genes Associated With Spikelet Development in Wheat

The number of spikelets per spike is an important trait that directly affects grain yield in wheat. Three quantitative trait loci (QTLs) associated with spikelet nodes per spike (SNS) were mapped in a population of recombinant inbred lines generated from a cross between two advanced breeding lines of winter wheat based on the phenotypic variation evaluated over six locations/years. Two of the three QTLs are QSns.sxau-2A at the WHEATFRIZZY PANICLE (WFZP) loci and QSns.sxau-7A at the WHEAT ORTHOLOG OF APO1 (WAPO1) loci. The WFZP-A1b allele with a 14-bp deletion at QSns.sxau-2A was associated with increased spikelets per spike. WAPO-A1e, as a novel allele at WAPO1, were regulated at the transcript level that was associated with the SNS trait. The third SNS QTL, QSns.sxau-7D on chromosome 7D, was not associated with homoeologous WAPO-D1 or any other genes known to regulate SNS. The favorable alleles for each of WZFP-A1, WAPO-A1, and QSns.sxau-7D are identified and incorporated to increase up to 3.4 spikelets per spike in the RIL lines. Molecular markers for the alleles were developed. This study has advanced our understanding of the genetic basis of natural variation in spikelet development in wheat.

A single nucleotide deletion in the third exon of FT-D1 increases the spikelet number and delays heading date in wheat (Triticum aestivum L.)

The spikelet number and heading date are two crucial and correlated traits for yield in wheat. Here, a quantitative trait locus (QTL) analysis was conducted in F₈ recombinant inbred lines (RILs) derived from crossing two common wheats with different spikelet numbers. A total of 15 stable QTL influencing total spikelet number (TSN) and heading date (HD) were detected. Notably, FT-D1, a well-known flowering time gene in wheat, was located within the finely mapped interval of a major QTL on 7DS (QTsn/Hd.cau-7D). A causal indel of one G in the third exon of FT-D1 was significantly associated with total spikelet number and heading date. Consistently, CRISPR/Cas9 mutant lines with homozygous mutations in FT-D1 displayed an increase in total spikelet number and heading date when compared with wild type. Moreover, one simple and robust marker developed according to the polymorphic site of FT-D1 revealed that this one G indel had been preferentially selected to adapt to different environments. Collectively, these data provide further insights into the genetic basis of spikelet number and heading date, and the diagnostic marker of FT-D1 will be useful for marker-assisted pyramiding in wheat breeding.

A Multi-Omics Approach for Rapid Identification of Large Genomic Lesions at the Wheat Dense Spike (wds) Locus

Optimal spike architecture provides a favorable structure for grain development and yield improvement. However, the number of genes cloned to underlie wheat spike architecture is extremely limited. Here, we obtained a wheat dense spike mutant (wds) induced by 60Co treatment of a common wheat landrace Huangfangzhu that exhibited significantly reduced spike and grain lengths. The shortened spike length was caused by longitudinal reduction in number and length of rachis cells. We adopted a multi-omics approach to identify the genomic locus underlying the wds mutant. We performed Exome Capture Sequencing (ECS) and identified two large deletion segments, named 6BL.1 at 334.8∼424.3 Mb and 6BL.2, 579.4∼717.8 Mb in the wds mutant. RNA-seq analysis confirmed that genes located in these regions lost their RNA expression. We then found that the 6BL.2 locus was overlapping with a known spike length QTL, qSL6B.2. Totally, 499 genes were located within the deleted region and two of them were found to be positively correlated with long spike accessions but not the ones with short spike. One of them, TraesCS6B01G334600, a well-matched homolog of the rice OsBUL1 gene that works in the Brassinosteroids (BR) pathway, was identified to be involved in cell size and number regulation. Further transcriptome analysis of young spikes showed that hormone-related genes were enriched among differentially expressed genes, supporting TraesCS6B01G334600 as a candidate gene. Our work provides a strategy to rapid locate genetic loci with large genomic lesions in wheat and useful resources for future wheat study.

Genome-Wide Association Study of Grain Number in Common Wheat From Shanxi Under Different Water Regimes

Water availability is a crucial environmental factor on grain number in wheat, which is one of the important yield-related traits. In this study, a diverse panel of 282 wheat accessions were phenotyped for grain number per spike (GNS), spikelet number (SN), basal sterile spikelet number (BSSN), and apical sterile spikelet number (ASSN) under different water regimes across two growing seasons. Correlation analysis showed that GNS is significantly correlated with both SN and BSSN under two water regimes. A total of 9,793 single nucleotide polymorphism (SNP) markers from the 15 K wheat array were employed for genome-wide association study (GWAS). A total of 77 significant marker-trait associations (MTAs) for investigated traits as well as 8 MTAs for drought tolerance coefficient (DTC) were identified using the mixed linear model. Favored alleles for breeding were inferred according to their estimated effects on GNS, based on the mean difference of varieties. Frequency changes in favored alleles associated with GNS in modern varieties indicate there is still considerable genetic potential for their use as markers for genome selection of GNS in wheat breeding.

Analysis of Genetic Regions Related to Field Grain Number per Spike From Chinese Wheat Founder Parent Linfen 5064

Wheat founder parents have been important in the development of new wheat cultivars. Understanding the effects of specific genome regions on yield-related traits in founder variety derivatives can enable more efficient use of these genetic resources through molecular breeding. In this study, the genetic regions related to field grain number per spike (GNS) from the founder parent Linfen 5064 were analyzed using a doubled haploid (DH) population developed from a cross between Linfen 5064 and Nongda 3338. Quantitative trait loci (QTL) for five spike-related traits over nine experimental locations/years were identified, namely, total spikelet number per spike (TSS), base sterile spikelet number per spike (BSSS), top sterile spikelet number per spike (TSSS), fertile spikelet number per spike (FSS), and GNS. A total of 13 stable QTL explaining 3.91-19.51% of the phenotypic variation were found. The effect of six of these QTL, Qtss.saw-2B.1, Qtss.saw-2B.2, Qtss.saw-3B, Qfss.saw-2B.2, Qbsss.saw-5A.1, and Qgns.saw-1A, were verified by another DH population (Linfen 5064/Jinmai 47), which showed extreme significance (P < 0.05) in more than three environments. No homologs of reported grain number-related from grass species were found in the physical regions of Qtss.saw-2B.1 and Qtss.saw-3B, that indicating both of them are novel QTL, or possess novel-related genes. The positive alleles of Qtss.saw-2B.2 from Linfen 5064 have the larger effect on TSS (3.30%, 0.62) and have 66.89% in Chinese cultivars under long-term artificial selection. This study revealed three key regions for GNS in Linfen 5064 and provides insights into molecular marker-assisted breeding.

Genetic Features of Triticale–Wheat Hybrids with Vaviloid-Type Spike Branching

Vaviloid spike branching, also called sham ramification, is a typical trait of Triticum vavilovii Jakubz. and is characterized by a lengthening of the spikelet axis. In this article, we present the results of a study of three triticale–wheat hybrid lines with differences in terms of the manifestation of the vaviloid spike branching. Lines were obtained by crossing triticale with hexaploid wheat, T. aestivum var. velutinum. The parental triticale is a hybrid of synthetic wheat (T. durum × Ae. tauschii var. meyrei) with rye, S. cereale ssp. segetale. Line 857 has a karyotype corresponding to hexaploid wheat and has a spike morphology closest to normal, whereas Lines 808/1 and 844/4 are characterized by the greatest manifestation of vaviloid spike branching. In Lines 808/1 and 844/4, we found the substitution 2RL(2DL). The karyotypes of the latter lines differ in that a pair of telocentric chromosomes 2DS is detected in Line 808/1, and these telocentrics are fused into one unpaired chromosome in Line 844/4. Using molecular genetic analysis, we found a deletion of the wheat domestication gene Q located on 5AL in the three studied hybrid lines. The deletion is local since an analysis of the adjacent gene B1 showed the presence of this gene. We assume that the manifestation of vaviloid spike branching in two lines (808/1 and 844/4) is associated with a disturbance in the joint action of genes Q and AP2L2-2D, which is another important gene that determines spike morphology and is located on 2DL.

Characterization of the Durum Wheat-Aegilops tauschii 4D(4B) Disomic Substitution Line YL-443 With Superior Characteristics of High Yielding and Stripe Rust Resistance

Durum wheat is one of the important food and cash crops. The main goals in current breeding programs are improving its low yield potential, kernel characteristics, and lack of resistance or tolerance to some biotic and abiotic stresses. In this study, a nascent synthesized hexaploid wheat Lanmai/AT23 is used as the female parent in crosses with its AB genome donor Lanmai. A tetraploid line YL-443 with supernumerary spikelets and high resistance to stripe rust was selected out from the pentaploid F₇ progeny. Somatic analysis using multicolor fluorescence in situ hybridization (mc-FISH) revealed that this line is a disomic substitution line with the 4B chromosome pair of Lanmai replaced by the 4D chromosome pair of Aegilops tauschii AT23. Comparing with Lanmai, YL-443 shows an increase in the number of spikelets and florets per spike by 36.3 and 75.9%, respectively. The stripe rust resistance gene Yr28 carried on the 4D chromosome was fully expressed in the tetraploid background. The present 4D(4B) disomic substitution line YL-443 was distinguished from the previously reported 4D(4B) lines with the 4D chromosomes from Chinese Spring (CS). Our study demonstrated that YL-443 can be used as elite germplasm for durum wheat breeding targeting high yield potential and stripe rust resistance. The Yr28-specific PCR marker and the 4D chromosome-specific KASP markers together with its unique features of pubescent leaf sheath and auricles can be utilized for assisting selection in breeding.