FRIZZY PANICLE drives supernumerary spikelets in bread wheat

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.

Photoperiod-1 regulates the wheat inflorescence transcriptome to influence spikelet architecture and flowering time

Photoperiod insensitivity has been selected by breeders to help adapt crops to diverse environments and farming practices. In wheat, insensitive alleles of Photoperiod-1 (Ppd-1) relieve the requirement of long daylengths to flower by promoting expression of floral promoting genes early in the season; however, these alleles also limit yield by reducing the number and fertility of grain-producing florets through processes that are poorly understood. Here, we performed transcriptome analysis of the developing inflorescence using near-isogenic lines that contain either photoperiod-insensitive or null alleles of Ppd-1, during stages when spikelet number is determined and floret development initiates. We report that Ppd-1 influences the stage-specific expression of genes with roles in auxin signaling, meristem identity, and protein turnover, and analysis of differentially expressed transcripts identified bZIP and ALOG transcription factors, namely PDB1 and ALOG1, which regulate flowering time and spikelet architecture. These findings enhance our understanding of genes that regulate inflorescence development and introduce new targets for improving yield potential.

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.

Identification and map-based cloning of an EMS-induced mutation in wheat gene TaSP1 related to spike architecture

KEY MESSAGE

A candidate gene TaSP1 related to spike shape was cloned, and the gene-specific marker was developed to efficiently track the superior haplotype in common wheat. Spike shape, an important factor that affects wheat grain yield, is mainly defined by spike length (SPL), spikelet number (SPN), and compactness. Zhoumai32 mutant 1160 (ZM1160), a mutant obtained from ethyl methane sulfonate (EMS) treatment of hexaploid wheat variety Zhoumai32, was used to identify and clone the candidate gene that conditioned the spike shape. Genetic analysis of an F2 population derived from a cross of ZM1160 and Bainong207 suggested that the compact spike shape in ZM1160 was controlled by a single recessive gene, and therefore, the mutated gene was designated as Tasp1. With polymorphic markers identified through bulked segregant analysis (BSA), the gene was mapped to a 2.65-cM interval flanked by markers YZU0852 and MIS46239 on chromosome 7D, corresponding to a 0.42-Mb physical interval of Chinese spring (CS) reference sequences (RefSeq v1.0). To fine map TaSP1, 15 and seven recombinants were, respectively, screened from 1599 and 1903 F3 plants derived from the heterozygous F2 plants. Finally, TaSP1 was delimited to a 21.9 Kb (4,870,562 to 4,892,493 bp) Xmis48123-Xmis48104 interval. Only one high-confidence gene TraesCS7D02G010200 was annotated in this region, which encodes an unknown protein with a putative vWA domain. Quantitative reverse transcription PCR (qRT-PCR) analysis showed that TraesCS7D02G010200 was mainly expressed in the spike. Haplotype analysis of 655 wheat cultivars using the candidate gene-specific marker Xg010200p2 identified a superior haplotype TaSP1b with longer spike and more spikelet number. TaSP1 is beneficial to the improvement in wheat spike shape.

Genetic insights into superior grain number traits: a QTL analysis of wheat-Agropyron cristatum derivative pubing3228

BACKGROUND

Agropyron cristatum (L.) is a valuable genetic resource for expanding the genetic diversity of common wheat. Pubing3228, a novel wheat-A. cristatum hybrid germplasm, exhibits several desirable agricultural traits, including high grain number per spike (GNS). Understanding the genetic architecture of GNS in Pubing3228 is crucial for enhancing wheat yield. This study aims to analyze the specific genetic regions and alleles associated with high GNS in Pubing3228.

METHODS

The study employed a recombination inbred line (RIL) population derived from a cross between Pubing3228 and Jing4839 to investigate the genetic regions and alleles linked to high GNS. Quantitative Trait Loci (QTL) analysis and candidate gene investigation were utilized to explore these traits.

RESULTS

A total of 40 QTLs associated with GNS were identified across 16 chromosomes, accounting for 4.25-17.17% of the total phenotypic variation. Five QTLs (QGns.wa-1D, QGns.wa-5 A, QGns.wa-7Da.1, QGns.wa-7Da.2 and QGns.wa-7Da.3) accounter for over 10% of the phenotypic variation in at least two environments. Furthermore, 94.67% of the GNS QTL with positive effects originated from Pubing3228. Candidate gene analysis of stable QTLs identified 11 candidate genes for GNS, including a senescence-associated protein gene (TraesCS7D01G148000) linked to the most significant SNP (AX-108,748,734) on chromosome 7D, potentially involved in reallocating nutrients from senescing tissues to developing seeds.

CONCLUSION

This study provides new insights into the genetic mechanisms underlying high GNS in Pubing3228, offering valuable resources for marker-assisted selection in wheat breeding to enhance yield.

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.

Grain yield trade-offs in spike-branching wheat can be mitigated by elite alleles affecting sink capacity and post-anthesis source activity

Introducing variations in inflorescence architecture, such as the 'Miracle-Wheat' (Triticum turgidum convar. compositum (L.f.) Filat.) with a branching spike, has relevance for enhancing wheat grain yield. However, in the spike-branching genotypes, the increase in spikelet number is generally not translated into grain yield advantage because of reduced grains per spikelet and grain weight. Here, we investigated if such trade-offs might be a function of source-sink strength by using 385 recombinant inbred lines developed by intercrossing the spike-branching landrace TRI 984 and CIRNO C2008, an elite durum (T. durum L.) cultivar; they were genotyped using the 25K array. Various plant and spike architectural traits, including flag leaf, peduncle, and spike senescence rate, were phenotyped under field conditions for 2 consecutive years. On chromosome 5AL, we found a new modifier QTL for spike branching, branched headt3 (bht-A3), which was epistatic to the previously known bht-A1 locus. Besides, bht-A3 was associated with more grains per spikelet and a delay in flag leaf senescence rate. Importantly, favourable alleles, viz. bht-A3 and grain protein content (gpc-B1) that delayed senescence, are required to improve grain number and grain weight in the spike-branching genotypes. In summary, achieving a balanced source-sink relationship might minimize grain yield trade-offs in Miracle-Wheat.

Dissection and validation of a promising QTL controlling spikelet number on 5B in bread wheat

KEY MESSAGE

Five environmentally stable QTLs for spikelet number per spike and days to heading were identified using a high-genetic map containing 95,444 SNPs, among which QSns.ucas-5B was validated using residual heterozygous line at multiple environments. Spikelet number per spike (SNS) and days to heading (DTH) play pivotal roles in the improvement of wheat yield. In this study, a high-density genetic map for a recombinant inbred lines (RILs) population derived from Zhengnong 17 (ZN17) and Yangbaimai (YBM) was constructed using 95,444 single-nucleotide polymorphism (SNP) markers from the Wheat660K SNP array. Our study identified a total of five environmentally stable QTLs for SNS and DTH, one of which was named QSns.ucas-5B, with a physical interval of approximately 545.4-552.1 Mb on the 5BL chromosome arm. Importantly, the elite haplotype within QSns.ucas-5B showed a consistent and positive effect on SNS, grain number and weight per spike, without extending the days to heading. These findings provide a foundation for future efforts to map and clone the gene(s) responsible for QSns.ucas-5B and further indicate the potential application of the developed and validated InDel marker of QSns.ucas-5B for molecular breeding purposes, aimed at improving wheat grain yield.

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.

Wheat bZIPC1 interacts with FT2 and contributes to the regulation of spikelet number per spike

KEY MESSAGE

The wheat transcription factor bZIPC1 interacts with FT2 and affects spikelet and grain number per spike. We identified a natural allele with positive effects on these two economically important traits. Loss-of-function mutations and natural variation in the gene FLOWERING LOCUS T2 (FT2) in wheat have previously been shown to affect spikelet number per spike (SNS). However, while other FT-like wheat proteins interact with bZIP-containing transcription factors from the A-group, FT2 does not interact with any of them. In this study, we used a yeast-two-hybrid screen with FT2 as bait and identified a grass-specific bZIP-containing transcription factor from the C-group, designated here as bZIPC1. Within the C-group, we identified four clades including wheat proteins that show Y2H interactions with different sets of FT-like and CEN-like encoded proteins. bZIPC1 and FT2 expression partially overlap in the developing spike, including the inflorescence meristem. Combined loss-of-function mutations in bZIPC-A1 and bZIPC-B1 (bzipc1) in tetraploid wheat resulted in a drastic reduction in SNS with a limited effect on heading date. Analysis of natural variation in the bZIPC-B1 (TraesCS5B02G444100) region revealed three major haplotypes (H1-H3), with the H1 haplotype showing significantly higher SNS, grain number per spike and grain weight per spike than both the H2 and H3 haplotypes. The favorable effect of the H1 haplotype was also supported by its increased frequency from the ancestral cultivated tetraploids to the modern tetraploid and hexaploid wheat varieties. We developed markers for the two non-synonymous SNPs that differentiate the bZIPC-B1b allele in the H1 haplotype from the ancestral bZIPC-B1a allele present in all other haplotypes. These diagnostic markers are useful tools to accelerate the deployment of the favorable bZIPC-B1b allele in pasta and bread wheat breeding programs.

Deciphering genetic basis of developmental and agronomic traits by integrating high-throughput optical phenotyping and genome-wide association studies in wheat

Dissecting the genetic basis of complex traits such as dynamic growth and yield potential is a major challenge in crops. Monitoring the growth throughout growing season in a large wheat population to uncover the temporal genetic controls for plant growth and yield-related traits has so far not been explored. In this study, a diverse wheat panel composed of 288 lines was monitored by a non-invasive and high-throughput phenotyping platform to collect growth traits from seedling to grain filling stage and their relationship with yield-related traits was further explored. Whole genome re-sequencing of the panel provided 12.64 million markers for a high-resolution genome-wide association analysis using 190 image-based traits and 17 agronomic traits. A total of 8327 marker-trait associations were detected and clustered into 1605 quantitative trait loci (QTLs) including a number of known genes or QTLs. We identified 277 pleiotropic QTLs controlling multiple traits at different growth stages which revealed temporal dynamics of QTLs action on plant development and yield production in wheat. A candidate gene related to plant growth that was detected by image traits was further validated. Particularly, our study demonstrated that the yield-related traits are largely predictable using models developed based on i-traits and provide possibility for high-throughput early selection, thus to accelerate breeding process. Our study explored the genetic architecture of growth and yield-related traits by combining high-throughput phenotyping and genotyping, which further unravelled the complex and stage-specific contributions of genetic loci to optimize growth and yield in wheat.

Delayed development of basal spikelets in wheat explains their increased floret abortion and rudimentary nature

Large differences exist in the number of grains per spikelet across an individual wheat (Triticum aestivum L.) spike. The central spikelets produce the highest number of grains, while apical and basal spikelets are less productive, and the most basal spikelets are commonly only developed in rudimentary form. Basal spikelets are delayed in initiation, yet they continue to develop and produce florets. The precise timing or the cause of their abortion remains largely unknown. Here, we investigated the underlying causes of basal spikelet abortion using shading applications in the field. We found that basal spikelet abortion is likely to be the consequence of complete floret abortion, as both occur concurrently and have the same response to shading treatments. We detected no differences in assimilate availability across the spike. Instead, we show that the reduced developmental age of basal florets pre-anthesis is strongly associated with their increased abortion. Using the developmental age pre-abortion, we were able to predict final grain set per spikelet across the spike, alongside the characteristic gradient in the number of grains from basal to central spikelets. Future efforts to improve spikelet homogeneity across the spike could thus focus on improving basal spikelet establishment and increasing floret development rates pre-abortion.

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.

Identification and validation of two major QTLs for spikelet number per spike in wheat (Triticum aestivum L.)

The total number of spikelets (TSPN) and the number of fertile spikelets (FSPN) affect the final number of grains per spikelet in wheat. This study constructed a high-density genetic map using 55K single nucleotide polymorphism (SNP) arrays from a population of 152 recombinant inbred lines (RIL) from crossing the wheat accessions 10-A and B39. Twenty-four quantitative trait loci (QTLs) for TSPN and 18 QTLs for FSPN were localized based on the phenotype in 10 environments in 2019-2021. Two major QTLs, QTSPN/QFSPN.sicau-2D.4 (34.43-47.43 Mb) and QTSPN/QFSPN.sicau-2D.5(32.97-34.43 Mb), explained 13.97%-45.90% of phenotypic variation. Linked kompetitive allele-specific PCR (KASP) markers further validated these two QTLs and revealed that QTSPN.sicau-2D.4 had less effect on TSPN than QTSPN.sicau-2D.5 in 10-A×BE89 (134 RILs) and 10-A×Chuannong 16 (192 RILs) populations, and one population of Sichuan wheat (233 accessions). The alleles combination haplotype 3 with the allele from 10-A of QTSPN/QFSPN.sicau-2D.5 and the allele from B39 of QTSPN.sicau-2D.4 resulted in the highest number of spikelets. In contrast, the allele from B39 for both loci resulted in the lowest number of spikelets. Using bulk-segregant analysis-exon capture sequencing, six SNP hot spots that included 31 candidate genes were identified in the two QTLs. We identified Ppd-D1a from B39 and Ppd-D1d from 10-A and further analyzed Ppd-D1 variation in wheat. These results identified loci and molecular markers with potential utility for wheat breeding and laid a foundation for further fine mapping and cloning of the two loci.

Current Progress and Future Prospect of Wheat Genetics Research towards an Enhanced Nitrogen Use Efficiency

To improve the yield and quality of wheat is of great importance for food security worldwide. One of the most effective and significant approaches to achieve this goal is to enhance the nitrogen use efficiency (NUE) in wheat. In this review, a comprehensive understanding of the factors involved in the process of the wheat nitrogen uptake, assimilation and remobilization of nitrogen in wheat were introduced. An appropriate definition of NUE is vital prior to its precise evaluation for the following gene identification and breeding process. Apart from grain yield (GY) and grain protein content (GPC), the commonly recognized major indicators of NUE, grain protein deviation (GPD) could also be considered as a potential trait for NUE evaluation. As a complex quantitative trait, NUE is affected by transporter proteins, kinases, transcription factors (TFs) and micro RNAs (miRNAs), which participate in the nitrogen uptake process, as well as key enzymes, circadian regulators, cross-talks between carbon metabolism, which are associated with nitrogen assimilation and remobilization. A series of quantitative genetic loci (QTLs) and linking markers were compiled in the hope to help discover more efficient and useful genetic resources for breeding program. For future NUE improvement, an exploration for other criteria during selection process that incorporates morphological, physiological and biochemical traits is needed. Applying new technologies from phenomics will allow high-throughput NUE phenotyping and accelerate the breeding process. A combination of multi-omics techniques and the previously verified QTLs and molecular markers will facilitate the NUE QTL-mapping and novel gene identification.

Identification and validation of quantitative trait loci for fertile spikelet number per spike and grain number per fertile spikelet in bread wheat (Triticum aestivum L.)

A major and stable QTL for fertile spikelet number per spike and grain number per fertile spikelet identified in a 4.96-Mb interval on chromosome 2A was validated in different genetic backgrounds. Fertile spikelet number per spike (FSN) and grain number per fertile spikelet (GNFS) contribute greatly to wheat yield improvement. To detect quantitative trait loci (QTL) associated with FSN and GNFS, we used a recombinant inbred line population crossed by Zhongkemai 13F10 and Chuanmai 42 in eight environments. Two Genomic regions associated with FSN were detected on chromosomes 2A and 6A using bulked segregant exome sequencing analysis. After the genetic linkage maps were constructed, four QTL QFsn.cib-2A, QFsn.cib-6A, QGnfs.cib-2A and QGnfs.cib-6A were identified in three or more environments. Among them, two major QTL QFsn.cib-2A (LOD = 4.67-9.34, PVE = 6.66-13.05%) and QGnfs.cib-2A (LOD = 5.27-11.68, PVE = 7.95-16.71%) were detected in seven and six environments, respectively. They were co-located in the same region, namely QFsn/Gnfs.cib-2A. The developed linked Kompetitive Allele Specific PCR (KASP) markers further validated this QTL in a different genetic background. QFsn/Gnfs.cib-2A showed pleiotropic effects on grain number per spike (GNS) and spike compactness (SC), and had no effect on grain weight. Since QFsn/Gnfs.cib-2A might be a new locus, it and the developed KASP markers can be used in wheat breeding. According to haplotype analysis, QFsn/Gnfs.cib-2A was identified as a target of artificial selection during wheat improvement. Based on haplotype analysis, sequence differences, spatiotemporal expression patterns, and gene annotation, the potential candidate genes for QFsn/Gnfs.cib-2A were predicted. These results provide valuable information for fine mapping and cloning gene(s) underlying QFsn/Gnfs.cib-2A.

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.

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.

Translational research from models to crops: comparative genomics for plant breeding

The concept of translational research, which originated in the medical field in the 1980s, consists in improving the efficient transfer of research results obtained in a species (which can be considered as a model or pivot) to all the species for which these results are of interest for its improvement in Agriculture. In this context, comparative genomics is an important tool for translational research, effectively identifying genes controlling common functions between species. Editing and phenotyping tools must thus allow the functional validation of the gene conserved within the species for which the knowledge has been extrapolated, that is to say transferred, and the identification of the best alleles and associated genotypes for exploitation in current breeding programs.

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.

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 control of branching patterns in grass inflorescences

Inflorescence branching in the grasses controls the number of florets and hence the number of seeds. Recent data on the underlying genetics come primarily from rice and maize, although new data are accumulating in other systems as well. This review focuses on a window in developmental time from the production of primary branches by the inflorescence meristem through to the production of glumes, which indicate the transition to producing a spikelet. Several major developmental regulatory modules appear to be conserved among most or all grasses. Placement and development of primary branches are controlled by conserved auxin regulatory genes. Subtending bracts are repressed by a network including TASSELSHEATH4, and axillary branch meristems are regulated largely by signaling centers that are adjacent to but not within the meristems themselves. Gradients of SQUAMOSA-PROMOTER BINDING-like and APETALA2-like proteins and their microRNA regulators extend along the inflorescence axis and the branches, governing the transition from production of branches to production of spikelets. The relative speed of this transition determines the extent of secondary and higher order branching. This inflorescence regulatory network is modified within individual species, particularly as regards formation of secondary branches. Differences between species are caused both by modifications of gene expression and regulators and by presence or absence of critical genes. The unified networks described here may provide tools for investigating orphan crops and grasses other than the well-studied maize and rice.

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.

MicroRNA-resistant alleles of HOMEOBOX DOMAIN-2 modify inflorescence branching and increase grain protein content of wheat

Plant and inflorescence architecture determine the yield potential of crops. Breeders have harnessed natural diversity for inflorescence architecture to improve yields, and induced genetic variation could provide further gains. Wheat is a vital source of protein and calories; however, little is known about the genes that regulate the development of its inflorescence. Here, we report the identification of semidominant alleles for a class III homeodomain-leucine zipper transcription factor, HOMEOBOX DOMAIN-2 (HB-2), on wheat A and D subgenomes, which generate more flower-bearing spikelets and enhance grain protein content. These alleles increase HB-2 expression by disrupting a microRNA 165/166 complementary site with conserved roles in plants; higher HB-2 expression is associated with modified leaf and vascular development and increased amino acid supply to the inflorescence during grain development. These findings enhance our understanding of genes that control wheat inflorescence development and introduce an approach to improve the nutritional quality of grain.

High-resolution detection of quantitative trait loci for seven important yield-related traits in wheat (Triticum aestivum L.) using a high-density SLAF-seq genetic map

BACKGROUND

Yield-related traits including thousand grain weight (TGW), grain number per spike (GNS), grain width (GW), grain length (GL), plant height (PH), spike length (SL), and spikelet number per spike (SNS) are greatly associated with grain yield of wheat (Triticum aestivum L.). To detect quantitative trait loci (QTL) associated with them, 193 recombinant inbred lines derived from two elite winter wheat varieties Chuanmai42 and Chuanmai39 were employed to perform QTL mapping in six/eight environments.

RESULTS

A total of 30 QTLs on chromosomes 1A, 1B, 1D, 2A, 2B, 2D, 3A, 4A, 5A, 5B, 6A, 6D, 7A, 7B and 7D were identified. Among them, six major QTLs QTgw.cib-6A.1, QTgw.cib-6A.2, QGw.cib-6A, QGl.cib-3A, QGl.cib-6A, and QSl.cib-2D explaining 5.96-23.75% of the phenotypic variance were detected in multi-environments and showed strong and stable effects on corresponding traits. Three QTL clusters on chromosomes 2D and 6A containing 10 QTLs were also detected, which showed significant pleiotropic effects on multiple traits. Additionally, three Kompetitive Allele Specific PCR (KASP) markers linked with five of these major QTLs were developed. Candidate genes of QTgw.cib-6A.1/QGl.cib-6A and QGl.cib-3A were analyzed based on the spatiotemporal expression patterns, gene annotation, and orthologous search.

CONCLUSIONS

Six major QTLs for TGW, GL, GW and SL were detected. Three KASP markers linked with five of these major QTLs were developed. These QTLs and KASP markers will be useful for elucidating the genetic architecture of grain yield and developing new wheat varieties with high and stable yield in wheat.

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.

Elbaiomy R, Chen L, Wang Z, Jiao J, Zhu J, Zhou W, Chen B, Soaud SA, Abbas M, Lin N, El-Sappah AH. Evolution analysis of FRIZZY PANICLE (FZP) orthologs explored the mutations in DNA coding sequences in the grass family (Poaceae)

FRIZZY PANICLE (FZP), an essential gene that controls spikelet differentiation and development in the grass family (Poaceae), prevents the formation of axillary bud meristems and is closely associated with crop yields. It is unclear whether the FZP gene or its orthologs were selected during the evolutionary process of grass species, which possess diverse spike morphologies. In the present study, we adopted bioinformatics methods for the evolutionary analysis of FZP orthologs in species of the grass family. Thirty-five orthologs with protein sequences identical to that of the FZP gene were identified from 29 grass species. Analysis of conserved domains revealed that the AP2/ERF domains were highly conserved with almost no amino acid mutations. However, species of the tribe Triticeae, genus Oryza, and C4 plants exhibited more significant amino acid mutations in the acidic C-terminus region. Results of the phylogenetic analysis showed that the 29 grass species could be classified into three groups, namely, Triticeae, Oryza, and C4 plants. Within the Triticeae group, the FZP genes originating from the same genome were classified into the same sub-group. When selection pressure analysis was performed, significant positive selection sites were detected in species of the Triticeae and Oryza groups. Our results show that the FZP gene was selected during the grass family's evolutionary process, and functional divergence may have already occurred among the various species. Therefore, researchers investigating the FZP gene's functions should take note of the possible presence of various roles in other grass species.

Meta-QTLs, ortho-meta-QTLs and candidate genes for grain yield and associated traits in wheat (Triticum aestivum L.)

In wheat, 2852 major QTLs of 8998 QTLs available for yield and related traits were used for meta-analysis; 141 meta-QTLs were identified, which included 13 breeder's MQTLs and 24 ortho-MQTLs; 1202 candidate genes and 50 homologues of genes for yield from other cereals were also identified. Meta-QTL analysis was conducted using 2852 of the 8998 known QTLs, retrieved from 230 reports published during 1999-2020 (including 19 studies on tetraploid wheat) for grain yield (GY) and the following ten component traits: (i) grain weight (GWei), (ii) grain morphology-related traits (GMRTs), (iii) grain number (GN), (iv) spikes-related traits (SRTs), (v) plant height (PH), (vi) tiller number (TN), (vii) harvest index (HI), (viii) biomass yield (BY), (ix) days to heading/flowering and maturity (DTH/F/M), and (x) grain filling duration (GFD). The study resulted in the identification of 141 meta-QTLs (MQTLs), with an average confidence interval (CI) of 1.4 cM as against a CI of > 12.1 cM (8.8 fold reduction) in the QTLs that were used. The corresponding physical length of CI ranged from 0.01 Mb to 661.9 Mb (mean, 31.5 Mb). Seventy-seven (77) of these 141 MQTLs overlapped marker-trait associations (MTAs) reported in genome-wide association studies. Also, 63 MQTLs (each based on at least 10 QTLs) were considered stable and robust, with 13 MQTLs described as breeder's MQTLs (selected based on small CI, large LOD, and high level of phenotypic variation explained). Thirty-five yield-related genes from rice, barley, and maize were also utilized to identify 50 wheat homologues in MQTLs. Further, the use of synteny and collinearity allowed the identification of 24 ortho-MQTLs which were common among the wheat, barley, rice, and maize. The results of the present study should prove useful for wheat breeding and future basic research in cereals including wheat, barley, rice, and maize.

Evolution of inflorescence branch modifications in cereal crops

Grasses are ubiquitous in our daily lives, with gramineous cereal crops such as maize, rice, and wheat constituting a large proportion of our daily staple food intake. Evolutionary forces, especially over the past ∼20 million years, have shaped grass adaptability, inflorescence architecture, and reproductive success. Here, we provide basic information on grass evolution and inflorescence structures mainly related to two inflorescence types: branched panicle- and spike-type inflorescences, the latter of which has highly modified branching. We summarize and compare known genetic pathways underlying each infloresecence type and discuss how the maize RAMOSA, rice ABERRANT PANICLE ORGANIZATION, and Triticeae COMPOSITUM pathways are regulated. Our analyses might lay the foundation for understanding species-specific gene regulatory networks that could result in improved sink capacities.

The barley mutant multiflorus2.b reveals quantitative genetic variation for new spikelet architecture

Spikelet indeterminacy and supernumerary spikelet phenotypes in barley multiflorus2.b mutant show polygenic inheritance. Genetic analysis of multiflorus2.b revealed major QTLs for spikelet determinacy and supernumerary spikelet phenotypes on 2H and 6H chromosomes. Understanding the genetic basis of yield forming factors in small grain cereals is of extreme importance, especially in the wake of stagnation of further yield gains in these crops. One such yield forming factor in these cereals is the number of grain-bearing florets produced per spikelet. Wild-type barley (Hordeum vulgare L.) spikelets are determinate structures, and the spikelet axis (rachilla) degenerates after producing single floret. In contrast, the rachilla of wheat (Triticum ssp.) spikelets, which are indeterminate, elongates to produce up to 12 florets. In our study, we characterized the barley spikelet determinacy mutant multiflorus2.b (mul2.b) that produced up to three fertile florets on elongated rachillae of lateral spikelets. Apart from the lateral spikelet indeterminacy (LS-IN), we also characterized the supernumerary spikelet phenotype in the central spikelets (CS-SS) of mul2.b. Through our phenotypic and genetic analyses, we identified two major QTLs on chromosomes 2H and 6H, and two minor QTLs on 3H for the LS-IN phenotype. For, the CS-SS phenotype, we identified one major QTL on 6H, and a minor QTL on 5H chromosomes. Notably, the 6H QTLs for CS-SS and LS-IN phenotypes co-located with each other, potentially indicating that a single genetic factor might regulate both phenotypes. Thus, our in-depth phenotyping combined with genetic analyses revealed the quantitative nature of the LS-IN and CS-SS phenotypes in mul2.b, paving the way for cloning the genes underlying these QTLs in the future.

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.

Brachypodium distachyon UNICULME4 and LAXATUM-A are redundantly required for development

In cultivated grasses, tillering, leaf, and inflorescence architecture, as well as abscission ability, are major agronomical traits. In barley (Hordeum vulgare), maize (Zea mays), rice (Oryza sativa), and brachypodium (Brachypodium distachyon), NOOT-BOP-COCH-LIKE (NBCL) genes are essential regulators of vegetative and reproductive development. Grass species usually possess 2-4 NBCL copies and until now a single study in O. sativa showed that the disruption of all NBCL genes strongly altered O. sativa leaf development. To improve our understanding of the role of NBCL genes in grasses, we extended the study of the two NBCL paralogs BdUNICULME4 (CUL4) and BdLAXATUM-A (LAXA) in the nondomesticated grass B. distachyon. For this, we applied reversed genetics and generated original B. distachyon single and double nbcl mutants by clustered regularly interspaced short palindromic repeats - CRISPR associated protein 9 (CRISPR-Cas9) approaches and genetic crossing between nbcl targeting induced local lesions in genomes (TILLING) mutants. Through the study of original single laxa CRISPR-Cas9 null alleles, we validated functions previously proposed for LAXA in tillering, leaf patterning, inflorescence, and flower development and also unveiled roles for these genes in seed yield. Furthermore, the characterization of cul4laxa double mutants revealed essential functions for nbcl genes in B. distachyon development, especially in the regulation of tillering, stem cell elongation and secondary cell wall composition as well as for the transition toward the reproductive phase. Our results also highlight recurrent antagonist interactions between NBCLs occurring in multiple aspects of B. distachyon development.

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.

Statistical estimates of multiple transcription factors binding in the model plant genomes based on ChIP-seq data

The development of high-throughput genomic sequencing coupled with chromatin immunoprecipitation technologies allows studying the binding sites of the protein transcription factors (TF) in the genome scale. The growth of data volume on the experimentally determined binding sites raises qualitatively new problems for the analysis of gene expression regulation, prediction of transcription factors target genes, and regulatory gene networks reconstruction. Genome regulation remains an insufficiently studied though plants have complex molecular regulatory mechanisms of gene expression and response to environmental stresses. It is important to develop new software tools for the analysis of the TF binding sites location and their clustering in the plant genomes, visualization, and the following statistical estimates. This study presents application of the analysis of multiple TF binding profiles in three evolutionarily distant model plant organisms. The construction and analysis of non-random ChIP-seq binding clusters of the different TFs in mammalian embryonic stem cells were discussed earlier using similar bioinformatics approaches. Such clusters of TF binding sites may indicate the gene regulatory regions, enhancers and gene transcription regulatory hubs. It can be used for analysis of the gene promoters as well as a background for transcription networks reconstruction. We discuss the statistical estimates of the TF binding sites clusters in the model plant genomes. The distributions of the number of different TFs per binding cluster follow same power law distribution for all the genomes studied. The binding clusters in Arabidopsis thaliana genome were discussed here in detail.

Analyses of MADS-box Genes Suggest HvMADS56 to Regulate Lateral Spikelet Development in Barley

MADS-box transcription factors are crucial regulators of inflorescence and flower development in plants. Therefore, the recent interest in this family has received much attention in plant breeding programs due to their impact on plant development and inflorescence architecture. The aim of this study was to investigate the role of HvMADS-box genes in lateral spikelet development in barley (Hordeum vulgare L.). A set of 30 spike-contrasting barley lines were phenotypically and genotypically investigated under controlled conditions. We detected clear variations in the spike and spikelet development during the developmental stages among the tested lines. The lateral florets in the deficiens and semi-deficiens lines were more reduced than in two-rowed cultivars except cv. Kristina. Interestingly, cv. Kristina, int-h.43 and int-i.39 exhibited the same behavior as def.5, def.6, semi-def.1, semi-def.8 regarding development and showed reduced lateral florets size. In HOR1555, HOR7191 and HOR7041, the lateral florets continued their development, eventually setting seeds. In contrast, lateral florets in two-rowed barley stopped differentiating after the awn primordia stage giving rise to lateral floret sterility. At harvest, the lines tested showed large variation for all central and lateral spikelet-related traits. Phylogenetic analysis showed that more than half of the 108 MADS-box genes identified are highly conserved and are expressed in different barley tissues. Re-sequence analysis of a subset of these genes showed clear polymorphism in either SNPs or in/del. Variation in HvMADS56 correlated with altered lateral spikelet morphology. This suggests that HvMADS56 plays an important role in lateral spikelet development in barley.

What are the regulatory targets for intervention in assimilate partitioning to improve crop yield and resilience?

Sucrose utilisation for the synthesis of cellular components involved in growth and development and the accumulation of biomass determines diversity in the plant kingdom; sucrose utilisation and partitioning also underpin crop yields. As a complex process the use of sucrose for the partitioning of plant products for yield is decided by the interaction of several regulatory hubs and the integration of metabolism and development. Understanding the regulation of assimilate partitioning has been a grand challenge in plant and crop science. There are emerging examples of genes and processes that appear important for assimilate partitioning that underpin yield in crops and which are amenable to intervention. Enzymes of carbon metabolism were some of the first targets in attempts to modify assimilate partitioning at the beginning (source) and end (sink) of the whole plant assimilate partitioning process. Metabolic enzymes are subject to regulatory and homeostatic mechanisms, a key factor to consider in modifying assimilate partitioning. Trehalose 6-phosphate, as a sucrose signal, may represent a special case in its ability to regulate and coordinate source and sink processes. This review summarises recent progress in understanding the underlying regulators of assimilate partitioning and the current and potentially most promising routes to crop yield enhancement with a main focus on cereals. A framework for how source-sink may regulate whole plant assimilate partitioning involving a few key elements and the central importance of reproductive development is presented.

Multi-locus genome-wide association mapping for spike-related traits in bread wheat (Triticum aestivum L.)

BACKGROUND

Bread wheat (Triticum aestivum L.) is one of the most important cereal food crops for the global population. Spike-layer uniformity (the consistency of the spike distribution in the vertical space)-related traits (SLURTs) are quantitative and have been shown to directly affect yield potential by modifying the plant architecture. Therefore, these parameters are important breeding targets for wheat improvement. The present study is the first genome-wide association study (GWAS) targeting SLURTs in wheat. In this study, a set of 225 diverse spring wheat accessions were used for multi-locus GWAS to evaluate SLURTs, including the number of spikes per plant (NSPP), spike length (SL), number of spikelets per spike (NSPS), grain weight per spike (GWPS), lowest tiller height (LTH), spike-layer thickness (SLT), spike-layer number (SLN) and spike-layer uniformity (SLU).

RESULTS

In total, 136 significant marker trait associations (MTAs) were identified when the analysis was both performed individually and combined for two environments. Twenty-nine MTAs were detected in environment one, 48 MTAs were discovered in environment two and 59 MTAs were detected using combined data from the two environments. Altogether, 15 significant MTAs were found for five traits in one of the two environments, and four significant MTAs were detected for the two traits, LTH and SLU, in both environments i.e. E1, E2 and also in combined data from the two environments. In total, 279 candidate genes (CGs) were identified, including Chaperone DnaJ, ABC transporter-like, AP2/ERF, SWEET sugar transporter, as well as genes that have previously been associated with wheat spike development, seed development and grain yield.

CONCLUSIONS

The MTAs detected through multi-locus GWAS will be useful for improving SLURTs and thus yield in wheat production through marker-assisted and genomic selection.

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.

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

Grain yield in bread wheat (Triticum aestivum L.) is largely determined by inflorescence architecture. Zang734 is an endemic Tibetan wheat variety that exhibits a rare triple spikelet (TRS) phenotype with significantly increased spikelet/floret number per spike. However, the molecular basis underlying this specific spike morphology is completely unknown. Through map-based cloning, the causal genes for TRS trait in Zang734 were isolated. Furthermore, using CRISPR/Cas9-based gene mutation, transcriptome sequencing and protein-protein interaction, the downstream signalling networks related to spikelet formation and awn elongation were defined. Results showed that the null mutation in WFZP-A together with deletion of WFZP-D led to the TRS trait in Zang734. More interestingly, WFZP plays a dual role in simultaneously repressing spikelet formation gene TaBA1 and activating awn development genes, basically through the recruitments of chromatin remodelling elements and the Mediator complex. Our findings provide insights into the molecular bases by which WFZP suppresses spikelet formation but promotes awn elongation and, more importantly, define WFZP-D as a favourable gene for high-yield crop breeding.

Genome-wide identification of loci modifying spike-branching in tetraploid wheat

Genetic modification of spike architecture is essential for improving wheat yield. Newly identified loci for the 'Miracle wheat' phenotype on chromosomes 1AS and 2BS have significant effects on spike traits. The wheat (Triticum ssp.) inflorescence, also known as a spike, forms an unbranched inflorescence in which the inflorescence meristem generates axillary spikelet meristems (SMs) destined to become sessile spikelets. Previously, we identified the putatively causative mutation in the branched head^t (bh^t) gene (TtBH-A1) of tetraploid wheat (T. turgidum convar. compositum (L.f.) Filat.) responsible for the loss of SM identity, converting the non-branching spike to a branched wheat spike. In the current study, we performed whole-genome quantitative trait loci (QTL) analysis using 146 recombinant inbred lines (RILs) derived from a cross between spike-branching wheat ('Miracle wheat') and an elite durum wheat cultivar showing broad phenotypic variation for spike architecture. Besides the previously found gene at the bh^t-A1 locus on the short arm of chromosome 2A, we also mapped two new modifier QTL for spike-branching on the short arm of chromosome 1A, termed bh^t-A2, and 2BS. Using biparental mapping population and GWAS in 302 diverse accessions, the 2BS locus was highly associated with coding sequence variation found at the homoeo-allele of TtBH-B1 (bh^t-B1). Thus, RILs that combined both bh^t-A1 and bh^t-B1 alleles showed additive genetic effects leading to increased penetrance and expressivity of the supernumerary spikelet and/or mini-spike formation.

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.

The Pharus latifolius genome bridges the gap of early grass evolution

The grass family (Poaceae) includes all commercial cereal crops and is a major contributor to biomass in various terrestrial ecosystems. The ancestry of all grass genomes includes a shared whole-genome duplication (WGD), named rho (ρ) WGD, but the evolutionary significance of ρ-WGD remains elusive. We sequenced the genome of Pharus latifolius, a grass species (producing a true spikelet) in the subfamily Pharoideae, a sister lineage to the core Poaceae including the (Panicoideae, Arundinoideae, Chloridoideae, Micrairoideae, Aristidoideae, and Danthonioideae (PACMAD) and Bambusoideae, Oryzoideae, and Pooideae (BOP) clades. Our results indicate that the P. latifolius genome has evolved slowly relative to cereal grass genomes, as reflected by moderate rates of molecular evolution, limited chromosome rearrangements and a low rate of gene loss for duplicated genes. We show that the ρ-WGD event occurred approximately 98.2 million years ago (Ma) in a common ancestor of the Pharoideae and the PACMAD and BOP grasses. This was followed by contrasting patterns of diploidization in the Pharus and core Poaceae lineages. The presence of two FRIZZY PANICLE-like genes in P. latifolius, and duplicated MADS-box genes, support the hypothesis that the ρ-WGD may have played a role in the origin and functional diversification of the spikelet, an adaptation in grasses related directly to cereal yields. The P. latifolius genome sheds light on the origin and early evolution of grasses underpinning the biology and breeding of cereals.

A Heterozygous Genotype-Dependent Branched-Spike Wheat and the Potential Genetic Mechanism Revealed by Transcriptome Sequencing

Simple Summary

This paper reported a novel type of branched spike wheat from a natural mutation event. The branched spike was controlled by a heterozygous genotype. The genetic patterns showed that gametophytic male sterility probably contributed to the heterozygous genotype responsible for the branched spike trait. Expressional levels and Wheat FRIZZY PANICLE (WFZP) sequencing between the mutant with the branched spike and the wild-type with the normal spike showed that WFZP was not the causal gene for the branched spike. Data from transcriptome sequencing indicated that carbohydrate metabolism might be involved in the formation of the branched spike trait.

Abstract

Wheat (Triticum aestivum L.) spike architecture is an important trait associated with spike development and grain yield. Here, we report a naturally occurring wheat mutant with branched spikelets (BSL) from its wild-type YD-16, which has a normal spike trait and confers a moderate level of resistance to wheat Fusarium head blight (FHB). The lateral meristems positioned at the basal parts of the rachis node of the BSL mutant develop into ramified spikelets characterized as multiple spikelets. The BSL mutant shows three to four-day longer growth period but less 1000-grain weight than the wild type, and it becomes highly susceptible to FHB infection, indicating that the locus controlling the BSL trait may have undergone an intensively artificial and/or natural selection in modern breeding process. The self-pollinated descendants of the lines with the BSL trait consistently segregated with an equal ratio of branched and normal spikelets (NSL) wheat, and homozygotes with the BSL trait could not be achieved even after nine cycles of self-pollination. Distinct segregation patterns both from the self-pollinated progenies of the BSL plants and from the reciprocal crosses between the BSL plants with their sister NSL plants suggested that gametophytic male sterility was probably associated with the heterozygosity for the BSL trait. Transcriptome sequencing of the RNA bulks contrasting in the two types of spike trait at the heading stage indicated that the genes on chromosome 2D_S may be critical for the BSL trait formation since 329 out of 2540 differentially expressed genes (DEGs) were located on that chromosome, and most of them were down-regulated. Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis showed that carbohydrate metabolism may be involved in the BSL trait expression. This work provides valuable clues into understanding development and domestication of wheat spike as well as the association of the BSL trait with FHB susceptibility.

Molecular Insights into Inflorescence Meristem Specification for Yield Potential in Cereal Crops

Flowering plants develop new organs throughout their life cycle. The vegetative shoot apical meristem (SAM) generates leaf whorls, branches and stems, whereas the reproductive SAM, called the inflorescence meristem (IM), forms florets arranged on a stem or an axis. In cereal crops, the inflorescence producing grains from fertilized florets makes the major yield contribution, which is determined by the numbers and structures of branches, spikelets and florets within the inflorescence. The developmental progression largely depends on the activity of IM. The proper regulations of IM size, specification and termination are outcomes of complex interactions between promoting and restricting factors/signals. Here, we focus on recent advances in molecular mechanisms underlying potential pathways of IM identification, maintenance and differentiation in cereal crops, including rice (Oryza sativa), maize (Zea mays), wheat (Triticum aestivum), and barley (Hordeum vulgare), highlighting the researches that have facilitated grain yield by, for example, modifying the number of inflorescence branches. Combinatorial functions of key regulators and crosstalk in IM determinacy and specification are summarized. This review delivers the knowledge to crop breeding applications aiming to the improvements in yield performance and productivity.

Genome-wide association study reveals the genetic basis of yield- and quality-related traits in wheat

BACKGROUND

Identifying the loci and dissecting the genetic architecture underlying wheat yield- and quality-related traits are essential for wheat breeding. A genome-wide association study was conducted using a high-density 90 K SNP array to analyze the yield- and quality-related traits of 543 bread wheat varieties.

RESULTS

A total of 11,140 polymorphic SNPs were distributed on 21 chromosomes, including 270 significant SNPs associated with 25 yield- and quality-related traits. Additionally, 638 putative candidate genes were detected near the significant SNPs based on BLUP data, including three (TraesCS7A01G482000, TraesCS4B01G343700, and TraesCS6B01G295400) related to spikelet number per spike, diameter of the first internode, and grain volume. The three candidate genes were further analyzed using stage- and tissue- specific gene expression data derived from an RNA-seq analysis. These genes are promising candidates for enhancing yield- and quality-related traits in wheat.

CONCLUSIONS

The results of this study provide a new insight to understand the genetic basis of wheat yield and quality. Furthermore, the markers detected in this study may be applicable for marker-assisted selection in wheat breeding programs.

Genetic factors associated with favourable pollinator traits in the wheat cultivar Piko

Hybrid breeding in wheat has the potential to boost yields. An efficient hybrid seed production system requires elite pollinators; however, such germplasm is limited among modern cultivars. Piko, a winter wheat (Triticum aestivum L.) cultivar, has been identified as a superior pollinator and has been used in Europe. Piko has favourable pollinator traits for anther extrusion, anther length, pollen mass and hybrid seed set. However, the genetic factors responsible for Piko's favourable traits are largely unknown. Here, we report on the genetic analysis of a Piko-derived F2 mapping population. We confirmed that Piko's Rht-D1a allele for tall stature is associated with large anthers and high anther extrusion. However, Rht-D1 was not found to be associated with anther filament length, confirmed by near isogenic lines. Piko's photoperiod sensitive Ppd-B1b allele shows an association with increased spike length, more spikelets and spike architecture traits, while the insensitive Ppd-B1a allele is linked with high anther extrusion and larger anthers. We identified an anther extrusion quantitative trait locus (QTL) on chromosome 6A that showed significantly biased transmission of the favourable Piko allele amongst F2 progenies. The Piko allele is completely absent in the distal 6AS region and the central 6A region revealed a significantly lower ratio (<8%) of F2 with homozygous Piko alleles. Our study provided further evidence for the effects of Rht-D1 and Ppd-B1 loci on multiple pollinator traits and a novel anther extrusion QTL that exhibits segregation distortion.

Identification and Validation of a Novel Locus Controlling Spikelet Number in Bread Wheat (Triticum aestivum L.)

Spikelet number is an important target trait for wheat yield improvement. Thus, the identification and verification of novel quantitative trait locus (QTL)/genes controlling spikelet number are essential for dissecting the underlying molecular mechanisms and hence for improving grain yield. In the present study, we constructed a high-density genetic map for the Kechengmai1/Chuanmai42 doubled haploid (DH) population using 13,068 single-nucleotide polymorphism (SNP) markers from the Wheat 55K SNP array. A comparison between the genetic and physical maps indicated high consistence of the marker orders. Based on this genetic map, a total of 27 QTLs associated with total spikelet number per spike (TSN) and fertile spikelet number per spike (FSN) were detected on chromosomes 1B, 1D, 2B, 2D, 3D, 4A, 4D, 5A, 5B, 5D, 6A, 6B, and 7D in five environments. Among them, five QTLs on chromosome 2D, 3D, 5A, and 7D were detected in multiple environments and combined QTL analysis, explaining the phenotypic variance ranging from 3.64% to 23.28%. Particularly, QTsn/Fsn.cib-3D for TSN and FSN [phenotypic variation explained (PVE) = 5.97-23.28%, limit of detection (LOD) = 3.73-18.51] is probably a novel locus and located in a 4.5-cM interval on chromosome arm 3DL flanking by the markers AX-110914105 and AX-109429351. This QTL was further validated in other two populations with different genetic backgrounds using the closely linked Kompetitive Allele-Specific PCR (KASP) marker KASP_AX-110914105. The results indicated that QTsn/Fsn.cib-3D significantly increased the TSN (5.56-7.96%) and FSN (5.13-9.35%), which were significantly correlated with grain number per spike (GNS). We also preliminary analyzed the candidate genes within this locus by sequence similarity, spatial expression patterns, and collinearity analysis. These results provide solid foundation for future fine mapping and cloning of QTsn/Fsn.cib-3D. The developed and validated KASP markers could be utilized in molecular breeding aiming to increase the grain yield in wheat.

The 55K SNP-Based Exploration of QTLs for Spikelet Number Per Spike in a Tetraploid Wheat (Triticum turgidum L.) Population: Chinese Landrace "Ailanmai" × Wild Emmer

Spikelet number per spike (SNS) is the primary factor that determines wheat yield. Common wheat breeding reduces the genetic diversity among elite germplasm resources, leading to a detrimental effect on future wheat production. It is, therefore, necessary to explore new genetic resources for SNS to increase wheat yield. A tetraploid landrace "Ailanmai" × wild emmer wheat recombinant inbred line (RIL) population was used to construct a genetic map using a wheat 55K single- nucleotide polymorphism (SNP) array. The linkage map containing 1,150 bin markers with a total genetic distance of 2,411.8 cm was obtained. Based on the phenotypic data from the eight environments and best linear unbiased prediction (BLUP) values, five quantitative trait loci (QTLs) for SNS were identified, explaining 6.71-29.40% of the phenotypic variation. Two of them, QSns.sau-AM-2B.2 and QSns.sau-AM-3B.2, were detected as a major and novel QTL. Their effects were further validated in two additional F₂ populations using tightly linked kompetitive allele-specific PCR (KASP) markers. Potential candidate genes within the physical intervals of the corresponding QTLs were predicted to participate in inflorescence development and spikelet formation. Genetic associations between SNS and other agronomic traits were also detected and analyzed. This study demonstrates the feasibility of the wheat 55K SNP array developed for common wheat in the genetic mapping of tetraploid population and shows the potential application of wheat-related species in wheat improvement programs.