Wheat VRN1, FUL2 and FUL3 play critical and redundant roles in spikelet development and spike determinacy

Phenotypic plasticity of flowering time and plant height related traits in wheat

BACKGROUND

Climate changes pose challenges to crop production. However, the causes of phenotypic differences across environments remain unclear.

RESULTS

Here, heading date (HD), flowering date (FD), and plant height (PH) were measured along with four environmental factors (day length (DL), growing degree days (GDD), precipitation (PRCP), and photothermal ratio (PTR)) to investigate the genetic basis of phenotypic plasticity of these traits in 616 wheat accessions using genome-wide association studies. Regarding quantitative trait locus-by-environment interactions (QEIs), five known and three candidate genes for HD, six known and seven candidate genes for FD, and four known and eighteen candidate genes for PH were identified. For the genes associated with phenotypic plasticity, 10 genes exhibited responsiveness to alterations in diverse environmental conditions according to transcriptome data; haplotype effects of 33 genes were identified as significantly correlated with the changes in environmental factors; six candidate genes were identified as hub genes in the gene network, possibly influencing other genes and causing the phenotypic plasticity. And over-dominant effects can explain over 50% the genetic variance of phenotypic plasticity. More importantly, one FD/HD candidate gene (TraesCS4A01G180700) and two PH candidate genes (TraesCS5B01G054800 and TraesCS2A01G539400) partly explain the phenotypic plasticity for the FD/HD and PH traits, respectively. In addition, the potential utilization of these genes in wheat breeding was discussed.

CONCLUSIONS

This study elucidated the genetic basis of phenotypic differences caused by environments and provided a foundation for addressing the impact of climate change on crop production.

Transcription factors BnaC09.FUL and BnaC06.WIP2 antagonistically regulate flowering time under long-day conditions in Brassica napus

Appropriate flowering time in rapeseed (Brassica napus L.) is vital for preventing losses from weather, diseases, and pests. However, the molecular basis of its regulation remains largely unknown. Here, a genome-wide association study identifies BnaC09.FUL, a MADS-box transcription factor, as a promising candidate gene regulating flowering time in B. napus. BnaC09.FUL expression increases sharply in B. napus shoot apices near bolting. BnaC09.FUL overexpression results in early flowering, while BnaFUL mutation causes delayed flowering in B. napus. A zinc finger transcription factor, BnaC06.WIP2, is identified as an interaction partner of BnaC09.FUL, and BnaC06.WIP2 overexpression delays flowering in B. napus, with RNA sequencing revealing its influence on the expression of many flowering-associated genes. We further demonstrate that BnaC06.WIP2 directly represses the expression of BnaA05.SOC1, BnaC03.SOC1, BnaC04.SOC1, BnaC06.FT, BnaA06.LFY, BnaC07.FUL, BnaA08.CAL, and BnaC03.CAL and indirectly inhibits the expression of other flowering time-related genes. Genetic and molecular investigations highlight the antagonistic relationship between BnaC09.FUL and BnaC06.WIP2 in regulating the flowering time in B. napus through direct regulation of the expression of BnaC03.SOC1, BnaA08.CAL, and BnaC03.CAL. Overall, our findings provide a mechanism by which the BnaC09.FUL-BnaC06.WIP2 transcriptional regulatory module controls the flowering time in B. napus.

QTL analysis of spike traits and KASP marker validation for basal sterile spikelet number in wheat

Spike traits play a pivotal role in determining wheat yield. In this study, two recombinant inbred line (RIL) populations were employed to identify quantitative trait loci (QTLs) for spike length (SL), spikelet number per spike (SNS), and basal sterile spikelet number (BSSN). A total of 30 QTLs were identified, including 7 major-effect QTLs that exhibited stability across multiple environments. Among these, two novel QTLs were discovered: QBssn.A-3 A.3 and QBssn.AV-5B.2. Furthermore, it was found that the combination of these two QTLs with other major effect QTLs had no impact on thousand-kernel weight (TKW). For QBssn.A-3 A.3, a Kompetitive Allele-Specific PCR (KASP) marker was developed and validated based on its physical location. Additionally, candidate genes within the intervals of QBssn.A-3 A.3 and QBssn.AV-5B.2 were predicted, based on these genes primarily associated with sterile spikelet number per spike. This study provides a foundation for future map-baesd cloning of genes related to spike traits, and for the breeding of high-yield wheat varieties.

Historic rewiring of grass flowering time pathways and implications for crop improvement under climate change

Grasses are fundamental to human survival, providing a large percentage of our calories, fuel, and fodder for livestock, and an enormous global carbon sink. A particularly important part of the grass plant is the grain-producing inflorescence that develops in response to both internal and external signals that converge at the shoot tip to influence meristem behavior. Abiotic signals that trigger reproductive development vary across the grass family, mostly due to the unique ecological and phylogenetic histories of each clade. The time it takes a grass to flower has implications for its ability to escape harsh environments, while also indirectly affecting abiotic stress tolerance, inflorescence architecture, and grain yield. Here, we synthesize recent insights into the evolution of grass flowering time in response to past climate change, particularly focusing on genetic convergence in underlying traits. We then discuss how and why the rewiring of a shared ancestral flowering pathway affects grass yields, and outline ways in which researchers are using this and other information to breed higher yielding, climate-proof cereal crops.

Decoupling the pleiotropic effects of VRT-A2 during reproductive development enhances wheat grain length and weight

VEGETATIVE TO REPRODUCTIVE TRANSITION 2 (VRT-A2) is a subspecies-forming gene that confers the long-glume and large-grain traits of tetraploid Polish wheat (Triticum polonicum; AABB) and hexaploid Xinjiang rice wheat (T. petropavlovskyi; AABBDD). Transcriptional activation of VRT-A2 due to a natural sequence variation in its Intron-1 region significantly enhances grain weight but also causes some basal spikelets to fail to completely develop, thus decreasing grain number per spike and yield. This yield penalty has presented a challenge for the use of VRT-A2 in breeding high-yield wheat. Here, we report the characterization of 2 regulatory modules that fine-tune VRT-A2 expression in bread wheat (T. aestivum): (i) the APETALA2/Ethylene Responsive Factor (AP2/ERF)-type transcription factor MULTI-FLORET SPIKELET1 (TaMFS1) represses VRT-A2 expression by recruiting a transcriptional corepressor and a histone deacetylase and (ii) the STRUCTURE-SPECIFIC RECOGNITION PROTEIN 1 (TaSSRP1) facilitates VRT-A2 activation by assembling Mediator and further RNA polymerase II. Deleting TaMFS1 triggered moderate upregulation of VRT-A2 results in significantly increased grain weight without the yield penalty. Our study thus provides a feasible strategy for overcoming the tradeoffs of pleotropic genes by editing their upstream transcriptional regulators.

From Images to Loci: Applying 3D Deep Learning to Enable Multivariate and Multitemporal Digital Phenotyping and Mapping the Genetics Underlying Nitrogen Use Efficiency in Wheat

The selection and promotion of high-yielding and nitrogen-efficient wheat varieties can reduce nitrogen fertilizer application while ensuring wheat yield and quality and contribute to the sustainable development of agriculture; thus, the mining and localization of nitrogen use efficiency (NUE) genes is particularly important, but the localization of NUE genes requires a large amount of phenotypic data support. In view of this, we propose the use of low-altitude aerial photography to acquire field images at a large scale, generate 3-dimensional (3D) point clouds and multispectral images of wheat plots, propose a wheat 3D plot segmentation dataset, quantify the plot canopy height via combination with PointNet++, and generate 4 nitrogen utilization-related vegetation indices via index calculations. Six height-related and 24 vegetation-index-related dynamic digital phenotypes were extracted from the digital phenotypes collected at different time points and fitted to generate dynamic curves. We applied height-derived dynamic numerical phenotypes to genome-wide association studies of 160 wheat cultivars (660,000 single-nucleotide polymorphisms) and found that we were able to locate reliable loci associated with height and NUE, some of which were consistent with published studies. Finally, dynamic phenotypes derived from plant indices can also be applied to genome-wide association studies and ultimately locate NUE- and growth-related loci. In conclusion, we believe that our work demonstrates valuable advances in 3D digital dynamic phenotyping for locating genes for NUE in wheat and provides breeders with accurate phenotypic data for the selection and breeding of nitrogen-efficient wheat varieties.

Regulation of tillering and panicle branching in rice and wheat

Branching is a critical aspect of plant architecture that significantly impacts the yield and adaptability of staple cereal crops like rice and wheat. Cereal crops develop tillers during the vegetative stage and panicle or spike branches during the reproductive stage, respectively, both of which are significantly impacted by hormones and genetic factors. Tillering and panicle branching are closely interconnected and exhibit high environmental plasticity. Here, we summarize the recent progress in genetic, hormonal, and environmental factors regulation in the branching of rice and wheat. This review not only provides a comprehensive overview of the current knowledge on branching mechanisms in rice and wheat, but also explores the prospects for future research aimed at optimizing crop architecture for enhanced productivity.

Tae-miR396b regulates TaGRFs in spikes of three wheat spike mutants

Tillering and spike differentiation are key agronomic traits for wheat (Triticum aestivum L.) production. Numerous studies have shown that miR396 and growth-regulating factor genes (GRFs) are involved in growth and development of different plant organs. Previously, we have reported that wheat miR396b (tae-miR396b) and their targets TaGRFs (T. aestivum GRFs) play important roles in regulating wheat tillering. This study was to investigate the regulatory roles of tae-miR396b and TaGRFs played during wheat spike development. Wheat cultivar Guomai 301 (wild type, WT) and its three sipke mutants dwarf round spike mutant (drs), apical spikelet sterility mutant (ass) and prematurely terminated spike differentiation mutant (ptsd1) were studied. Three homeologous genes of tae-miR396b on the long arms of chromosomes 6A, 6B, and 6D were identified, and they encoded the same mature miRNA. Complementary sequences of mature tae-miR396b were identified in 23 TaGRFs, indicating they were the target genes of tae-miR396b. Tae-miR396b had different regulatory effects on TaGRFs between Guomai 301 and its mutants. TaGRF2-7A was confirmed to be the target gene of tae-miR396b by molecular interaction assay. The expression levels of tae-miR396b and TaGRFs were different between WT and mutants drs, ass and ptsd1 at the floret primordium visible (S1), the two awns/spikelet reaching apical meristem of the spikelet (S2), and the green anther stage (S3). The expression level of tae-miR396b in WT was significantly higher than that in mutants drs and ass. The most TaGRFs were negatively regulated by tae-miR396b. The abnormal expressions of TaGRF1 (6A, 6D), TaGRF2 (7A, 7B, 7D), TaGRF4 (6A, 6B), TaGRF5 (4A, 7A, 7D), and TaGRF10 (6A, 6B, 6D) were important causes for abnormal spike development in the three mutants. This study laid foundation for further elucidating functions of tae-miR396b and TaGRFs underlying wheat spike development. Regulating tae-miR396b and TaGRFs will be a new approach for wheat high yield breeding.

Identification and validation of two quantitative trait loci showing pleiotropic effect on multiple grain-related traits in bread wheat (Triticum aestivum L.)

KEY MESSAGE

QKl/Tgw/Gns.yaas-2D associates with KL, TGW, and GNS, and QKl/Tgw.yaas-5A associates with KL and TGW. Significantly pleiotropic and additive effects of these two QTL were validated. The YM5 allele both at QKl/Tgw/Gns.yaas-2D and QKl/Tgw.yaas-5A was proved to be the best allelic combination for improving yield potential. Kernel length (KL), kernel width (KW), thousand grain weight (TGW), and grain number per spike (GNS) play important roles in the yield improvement of wheat. In this study, one recombinant inbred line (RIL) derived from a cross between Yangmai 5 (YM5) and Yanzhan 1 (YZ1) was used to identify quantitative trait loci (QTL) associated with KL, KW, TGW, and GNS across three years. Two pleiotropic QTL namely QKl/Tgw/Gns.yaas-2D and QKl/Tgw.yaas-5A were located in two genomic regions on chromosomes 2D and 5A, respectively. Breeder-friendly Kompetitive Allele-Specific PCR (KASP) markers for QKl/Tgw/Gns.yaas-2D and QKl/Tgw.yaas-5A were developed and validated in a set of 246 wheat cultivars/lines. Analysis of allelic combinations indicated that the YM5 allele both at QKl/Tgw/Gns.yaas-2D and QKl/Tgw.yaas-5A is probably the best one to promote TGW, GNS, and grain weight per spike. Based on the analysis of gene annotation, sequence variations, expression patterns, and GO enrichment, twenty-five and twenty-four candidate genes of QKl/Tgw/Gns.yaas-2D and QKl/Tgw.yaas-5A, respectively, were identified. These results provide the basis of fine-mapping the target QTL and marker-assisted selection in wheat yield-breeding programs.

Natural alleles of LEAFY and WAPO1 interact to regulate spikelet number per spike in wheat

KEY MESSAGE

Specific combinations of LFY and WAPO1 natural alleles maximize spikelet number per spike in wheat. Spikelet number per spike (SNS) is an important yield component in wheat that determines the maximum number of grains that can be formed in a wheat spike. In wheat, loss-of-function mutations in LEAFY (LFY) or its interacting protein WHEAT ORTHOLOG OF APO1 (WAPO1) significantly reduce SNS by reducing the rate of formation of spikelet meristems. In previous studies, we identified a natural amino acid change in WAPO1 (C47F) that significantly increases SNS in hexaploid wheat. In this study, we searched for natural variants in LFY that were associated with differences in SNS and detected significant effects in the LFY-B region in a nested association mapping population. We generated a large mapping population and confirmed that the LFY-B polymorphism R80S is linked with the differences in SNS, suggesting that LFY-B is the likely causal gene. A haplotype analysis revealed two amino acid changes P34L and R80S, which were both enriched during wheat domestication and breeding suggesting positive selection. We also explored the interactions between the LFY and WAPO1 natural variants for SNS using biparental populations and identified significant interaction, in which the positive effect of the 80S and 34L alleles from LFY-B was only detected in the WAPO-A1 47F background but not in the 47C background. Based on these results, we propose that the allele combination WAPO-A1-47F/LFY-B 34L 80S can be used in wheat breeding programs to maximize SNS and increase grain yield potential in wheat.

Distinct roles of H3K27me3 and H3K36me3 in vernalization response, maintenance, and resetting in winter wheat

Winter plants rely on vernalization, a crucial process for adapting to cold conditions and ensuring successful reproduction. However, understanding the role of histone modifications in guiding the vernalization process in winter wheat remains limited. In this study, we investigated the transcriptome and chromatin dynamics in the shoot apex throughout the life cycle of winter wheat in the field. Two core histone modifications, H3K27me3 and H3K36me3, exhibited opposite patterns on the key vernalization gene VERNALIZATION1 (VRN1), correlating with its induction during cold exposure. Moreover, the H3K36me3 level remained high at VRN1 after cold exposure, which may maintain its active state. Mutations in FERTILIZATION-INDEPENDENT ENDOSPERM (TaFIE) and SET DOMAIN GROUP 8/EARLY FLOWERING IN SHORT DAYS (TaSDG8/TaEFS), components of the writer complex for H3K27me3 and H3K36me3, respectively, affected flowering time. Intriguingly, VRN1 lost its high expression after the cold exposure memory in the absence of H3K36me3. During embryo development, VRN1 was silenced with the removal of active histone modifications in both winter and spring wheat, with selective restoration of H3K27me3 in winter wheat. The mutant of Tafie-cr-87, a component of H3K27me3 "writer" complex, did not influence the silence of VRN1 during embryo development, but rather attenuated the cold exposure requirement of winter wheat. Integrating gene expression with H3K27me3 and H3K36me3 patterns identified potential regulators of flowering. This study unveils distinct roles of H3K27me3 and H3K36me3 in controlling vernalization response, maintenance, and resetting in winter wheat.

Genome-wide identification of the E-class gene family in wheat: evolution, expression, and interaction

Introduction

Wheat (Triticum aestivum L.) is among themost important crop worldwide. Given a growing population and changing climate, enhancing wheat yield is of great importance. Yield is closely associated with flower and spike development, and E-class genes play important roles in the flower and kernel development of plants. Currently, the absence of systematic analysis on the E gene family hinders our comprehension of their roles in plant growth and development.

Methods

Identify E-class genes based on homologous sequence searches. Analyze the identified E-class genes through a series of gene family analyses. Determine the expression levels of wheat E-class genes by searching public databases. Validate the functions of these genes by transforming them into Arabidopsis. Finally, determine the interactions between the genes through yeast two-hybrid experiments.

Results

Fifteen E-class genes (TaEs) were identified in common wheat. Nine E-class genes were detected in five ancestral/closely related species, including one in Aegilops tauschii (AtE), one in T. Urartu (TuEs), two in T. turgidum (TtEs), two in T. dicoccoides (TdEs), and three in T. spelta (TsEs). The 24 E-class genes were classified into three subgroups using a phylogenetic approach. All genes were highly expressed in spikes, and most were only highly expressed at the floret meristem stage. The effects of TaSEP5-A on flowering and growth cycles were confirmed in homologous mutants and transgenic Arabidopsis thaliana. The E-class genes were able to regulate the growth cycle of Arabidopsis. Finally, we confirmed the interactions between TaSEP5-A and other wheat E-class genes based on yeast two-hybrid assays.

Discussion

Our findings provide information regarding the E-class genes in wheat and will potentially promote the application of these genes in wheat improvement.

Independent genetic factors control floret number and spikelet number in Triticum turgidum ssp

Wheat grain yield is a complex trait resulting from a trade-off among many distinct components. During wheat evolution, domestication events and then modern breeding have strongly increased the yield potential of wheat plants, by enhancing spike fertility. To address the genetic bases of spike fertility in terms of spikelet number per spike and floret number per spikelet, a population of 110 recombinant inbred lines (RILS) obtained crossing a Triticum turgidum ssp. durum cultivar (Latino) and a T. dicoccum accession (MG5323) was exploited. Being a modern durum and a semi-domesticated genotype, respectively, the two parents differ for spike architecture and fertility, and thus the corresponding RIL population is the ideal genetic material to dissect genetic bases of yield components. The RIL population was phenotyped in four environments. Using a high-density SNP genetic map and taking advantage of several genome sequencing available for Triticeae, a total of 94 QTLs were identified for the eight traits considered; these QTLs were further reduced to 17 groups, based on their genetic and physical co-location. QTLs controlling floret number per spikelet and spikelet number per spike mapped in non-overlapping chromosomal regions, suggesting that independent genetic factors determine these fertility-related traits. The physical intervals of QTL groups were considered for possible co-location with known genes functionally involved in spike fertility traits and with yield-related QTLs previously mapped in tetraploid wheat. The most interesting result concerns a QTL group on chromosome 5B, associated with spikelet number per spike, since it could host genes still uncharacterized for their association to spike fertility. Finally, we identified two different regions where the trade-off between fertility related traits and kernel weight is overcome. Further analyses of these regions could pave the way for a future identification of new genetic loci contributing to fertility traits essential for yield improvement in durum wheat.

LEAFY and WAPO1 jointly regulate spikelet number per spike and floret development in wheat

In wheat, the transition of the inflorescence meristem to a terminal spikelet (IM→TS) determines the spikelet number per spike (SNS), an important yield component. In this study, we demonstrate that the plant-specific transcription factor LEAFY (LFY) physically and genetically interacts with WHEAT ORTHOLOG OF APO1 (WAPO1) to regulate SNS and floret development. Loss-of-function mutations in either or both genes result in significant and similar reductions in SNS, as a result of a reduction in the rate of spikelet meristem formation per day. SNS is also modulated by significant genetic interactions between LFY and the SQUAMOSA MADS-box genes VRN1 and FUL2, which promote the IM→TS transition. Single-molecule fluorescence in situ hybridization revealed a downregulation of LFY and upregulation of the SQUAMOSA MADS-box genes in the distal part of the developing spike during the IM→TS transition, supporting their opposite roles in the regulation of SNS in wheat. Concurrently, the overlap of LFY and WAPO1 transcription domains in the developing spikelets contributes to normal floret development. Understanding the genetic network regulating SNS is a necessary first step to engineer this important agronomic trait.

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.

Evolution of major flowering pathway integrators in Orchidaceae

The Orchidaceae is a mega-diverse plant family with ca. 29,000 species with a large variety of life forms that can colonize transitory habitats. Despite this diversity, little is known about their flowering integrators in response to specific environmental factors. During the reproductive transition in flowering plants a vegetative apical meristem (SAM) transforms into an inflorescence meristem (IM) that forms bracts and flowers. In model grasses, like rice, a flowering genetic regulatory network (FGRN) controlling reproductive transitions has been identified, but little is known in the Orchidaceae. In order to analyze the players of the FRGN in orchids, we performed comprehensive phylogenetic analyses of CONSTANS-like/CONSTANS-like 4 (COL/COL4), FLOWERING LOCUS D (FD), FLOWERING LOCUS C/FRUITFULL (FLC/FUL) and SUPRESSOR OF OVEREXPRESSION OF CONSTANS 1 (SOC1) gene lineages. In addition to PEBP and AGL24/SVP genes previously analyzed, here we identify an increase of orchid homologs belonging to COL4, and FUL gene lineages in comparison with other monocots, including grasses, due to orchid-specific gene lineage duplications. Contrariwise, local duplications in Orchidaceae are less frequent in the COL, FD and SOC1 gene lineages, which points to a retention of key functions under strong purifying selection in essential signaling factors. We also identified changes in the protein sequences after such duplications, variation in the evolutionary rates of resulting paralogous clades and targeted expression of isolated homologs in different orchids. Interestingly, vernalization-response genes like VERNALIZATION1 (VRN1) and FLOWERING LOCUS C (FLC) are completely lacking in orchids, or alternatively are reduced in number, as is the case of VERNALIZATION2/GHD7 (VRN2). Our findings point to non-canonical factors sensing temperature changes in orchids during reproductive transition. Expression data of key factors gathered from Elleanthus auratiacus, a terrestrial orchid in high Andean mountains allow us to characterize which copies are actually active during flowering. Altogether, our data lays down a comprehensive framework to assess gene function of a restricted number of homologs identified more likely playing key roles during the flowering transition, and the changes of the FGRN in neotropical orchids in comparison with temperate grasses.

Finding the right balance: The enduring role of florigens during cereal inflorescence development and their influence on fertility

Flowering is a vital process in a plant's lifecycle and variation for flowering-time has helped cereals adapt to diverse environments. Much cereal research has focused on understanding how flowering signals, or florigens, regulate the floral transition and timing of ear emergence. However, flowering genes also perform an enduring role during inflorescence development, with genotypes that elicit a weaker flowering signal producing more elaborately branched inflorescences with extra floret-bearing spikelets. While this outcome indicates that variable expression of flowering genes could boost yield potential, further analysis has shown that dampened florigen levels can compromise fertility, negating the benefit of extra grain-producing sites. Here, we discuss ways that florigens contribute to early and late inflorescence development, including their influence on branch/spikelet architecture and fertility. We propose that a deeper understanding of the role for florigens during inflorescence development could be used to balance the effects of florigens throughout flowering to improve productivity.

Genome-wide association study and genomic selection of spike-related traits in bread wheat

KEY MESSAGE

Identification of 337 stable MTAs for wheat spike-related traits improved model accuracy, and favorable alleles of MTA259 and MTA64 increased grain weight and yield per plant. Wheat (Triticum aestivum L.) is one of the three primary global, staple crops. Improving spike-related traits in wheat is crucial for optimizing spike and plant morphology, ultimately leading to increased grain yield. Here, we performed a genome-wide association study using a dataset of 24,889 high-quality unique single-nucleotide polymorphisms (SNPs) and phenotypic data from 314 wheat accessions across eight diverse environments. In total, 337 stable and significant marker-trait associations (MTAs) related to spike-related traits were identified. MTA259 and MTA64 were consistently detected in seven and six environments, respectively. The presence of favorable alleles associated with MTA259 and MTA64 significantly reduced wheat spike exsertion length and spike length, while enhancing thousand kernel weight and yield per plant. Combined gene expression and network analyses identified TraesCS6D03G0692300 and TraesCS6D03G0692700 as candidate genes for MTA259 and TraesCS2D03G0111700 and TraesCS2D03G0112500 for MTA64. The identified MTAs significantly improved the prediction accuracy of each model compared with using all the SNPs, and the random forest model was optimal for genome selection. Additionally, the eight stable and major MTAs, including MTA259, MTA64, MTA66, MTA94, MTA110, MTA165, MTA180, and MTA164, were converted into cost-effective and efficient detection markers. This study provided valuable genetic resources and reliable molecular markers for wheat breeding programs.

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.

The Function of Florigen in the Vegetative-to-Reproductive Phase Transition in and around the Shoot Apical Meristem

Plants undergo a series of developmental phases throughout their life-cycle, each characterized by specific processes. Three critical features distinguish these phases: the arrangement of primordia (phyllotaxis), the timing of their differentiation (plastochron) and the characteristics of the lateral organs and axillary meristems. Identifying the unique molecular features of each phase, determining the molecular triggers that cause transitions and understanding the molecular mechanisms underlying these transitions are keys to gleaning a complete understanding of plant development. During the vegetative phase, the shoot apical meristem (SAM) facilitates continuous leaf and stem formation, with leaf development as the hallmark. The transition to the reproductive phase induces significant changes in these processes, driven mainly by the protein FT (FLOWERING LOCUS T) in Arabidopsis and proteins encoded by FT orthologs, which are specified as 'florigen'. These proteins are synthesized in leaves and transported to the SAM, and act as the primary flowering signal, although its impact varies among species. Within the SAM, florigen integrates with other signals, culminating in developmental changes. This review explores the central question of how florigen induces developmental phase transition in the SAM. Future research may combine phase transition studies, potentially revealing the florigen-induced developmental phase transition in the SAM.

Analysis of pea mutants reveals the conserved role of FRUITFULL controlling the end of flowering and its potential to boost yield

Monocarpic plants have a single reproductive phase in their life. Therefore, flower and fruit production are restricted to the length of this period. This reproductive strategy involves the regulation of flowering cessation by a coordinated arrest of the growth of the inflorescence meristems, optimizing resource allocation to ensure seed filling. Flowering cessation appears to be a regulated phenomenon in all monocarpic plants. Early studies in several species identified seed production as a major factor triggering inflorescence proliferative arrest. Recently, genetic factors controlling inflorescence arrest, in parallel to the putative signals elicited by seed production, have started to be uncovered in Arabidopsis, with the MADS-box gene FRUITFULL (FUL) playing a central role in the process. However, whether the genetic network regulating arrest is also at play in other species is completely unknown. Here, we show that this role of FUL is not restricted to Arabidopsis but is conserved in another monocarpic species with a different inflorescence structure, field pea, strongly suggesting that the network controlling the end of flowering is common to other plants. Moreover, field trials with lines carrying mutations in pea FUL genes show that they could be used to boost crop yield.

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.

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.

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

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

INDETERMINATE1-mediated expression of FT family genes is required for proper timing of flowering in Brachypodium distachyon

The transition to flowering is a major developmental switch in plants. In many temperate grasses, perception of indicators of seasonal change, such as changing day-length and temperature, leads to expression of FLOWERING LOCUS T1 (FT1) and FT-Like (FTL) genes that are essential for promoting the transition to flowering. However, little is known about the upstream regulators of FT1 and FTL genes in temperate grasses. Here, we characterize the monocot-specific gene INDETERMINATE1 (BdID1) in Brachypodium distachyon and demonstrate that BdID1 is a regulator of FT family genes. Mutations in ID1 impact the ability of the short-day (SD) vernalization, cold vernalization, and long-day (LD) photoperiod pathways to induce certain FTL genes. BdID1 is required for upregulation of FTL9 (FT-LIKE9) expression by the SD vernalization pathway, and overexpression of FTL9 in an id1 background can partially restore the delayed flowering phenotype of id1. We show that BdID1 binds in vitro to the promoter region of FTL genes suggesting that ID1 directly activates FTL expression. Transcriptome analysis shows that BdID1 is required for FT1, FT2, FTL12, and FTL13 expression under inductive LD photoperiods, indicating that BdID1 is a regulator of the FT gene family. Moreover, overexpression of FT1 in the id1 background results in rapid flowering similar to overexpressing FT1 in the wild type, demonstrating that BdID1 is upstream of FT family genes. Interestingly, ID1 negatively regulates a previously uncharacterized FTL gene, FTL4, and we show that FTL4 is a repressor of flowering. Thus, BdID1 is critical for proper timing of flowering in temperate grasses.

Organ-enriched gene expression during floral morphogenesis in wild barley

Floral morphology varies considerably between dicots and monocots. The ABCDE model explaining how floral organ development is controlled was formulated using core eudicots and applied to grass crops. Barley (Hordeum. vulgare) has unique floral morphogenesis. Wild barley (H. vulgare ssp. spontaneum), which is the immediate ancestor of cultivated barley (H. vulgare ssp. vulgare), contains a rich reservoir of genetic diversity. However, the wild barley genes involved in floral organ development are still relatively uncharacterized. In this study, we generated an organ-specific transcriptome atlas for wild barley floral organs. Genome-wide transcription profiles indicated that 22 838 protein-coding genes were expressed in at least one organ. These genes were grouped into seven clusters according to the similarities in their expression patterns. Moreover, 5619 genes exhibited organ-enriched expression, 677 of which were members of 47 transcription factor families. Gene ontology analyses suggested that the functions of the genes with organ-enriched expression influence the biological processes in floral organs. The co-expression regulatory network showed that the expression of 690 genes targeted by MADS-box proteins was highly positively correlated with the expression of ABCDE model genes during floral morphogenesis. Furthermore, the expression of 138 genes was specific to the wild barley OUH602 genome and not the Morex genome; most of these genes were highly expressed in the glume, awn, lemma, and palea. This study revealed the global gene expression patterns underlying floral morphogenesis in wild barley. On the basis of the study findings, a molecular mechanism controlling floral morphology in barley was proposed.

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.

FUL homologous gene CmFL1 is involved in regulating flowering time and floret numbers in Chrysanthemum morifolium

Flowering time and floret numbers are important ornamental characteristics of chrysanthemums that control their adaptability and inflorescence morphology, respectively. The FRUITFULL (FUL) gene plays a key role in inducing flowering and inflorescence meristem development. In this study, we isolated a homolog of the MADS-box gene FUL, CmFUL-Like 1 (CmFL1), from chrysanthemum inflorescence buds. Quantitative RT-PCR and in situ analyses showed that CmFL1 was strongly expressed in young inflorescence buds. Overexpression of CmFL1 caused early flowering while co-suppression expression of CmFL1 increased the number of florets. Furthermore, the floral promoting factors CmSOC1, CmFDL1, and CmLFY were up-regulated in the shoot tips of transgenic plants. In addition, RNA-seq analysis of the transgenic plants suggested that certain differentially expressed genes (DEGs)-such as MADS-box, homeobox family, and ethylene pathway genes-may be involved in the inflorescence meristem development. GO pathway enrichment analysis showed that the differentially transcribed genes enriched the representation of the carbohydrate metabolic pathway, which is critical for flower development. Overall, our findings revealed the conserved function of CmFL1 in controlling flowering time along with a novel function in regulating the number of florets.

Multi-locus genome-wide association studies reveal the genetic architecture of Fusarium head blight resistance in durum wheat

Durum wheat is more susceptible to Fusarium head blight (FHB) than other types or classes of wheat. The disease is one of the most devastating in wheat; it reduces yield and end-use quality and contaminates the grain with fungal mycotoxins such as deoxynivalenol (DON). A panel of 265 Canadian and European durum wheat cultivars, as well as breeding and experimental lines, were tested in artificially inoculated field environments (2019-2022, inclusive) and two greenhouse trials (2019 and 2020). The trials were assessed for FHB severity and incidence, visual rating index, Fusarium-damaged kernels, DON accumulation, anthesis or heading date, maturity date, and plant height. In addition, yellow pigment and protein content were analyzed for the 2020 field season. To capture loci underlying FHB resistance and related traits, GWAS was performed using single-locus and several multi-locus models, employing 13,504 SNPs. Thirty-one QTL significantly associated with one or more FHB-related traits were identified, of which nine were consistent across environments and associated with multiple FHB-related traits. Although many of the QTL were identified in regions previously reported to affect FHB, the QTL QFhb-3B.2, associated with FHB severity, incidence, and DON accumulation, appears to be novel. We developed KASP markers for six FHB-associated QTL that were consistently detected across multiple environments and validated them on the Global Durum Panel (GDP). Analysis of allelic diversity and the frequencies of these revealed that the lines in the GDP harbor between zero and six resistance alleles. This study provides a comprehensive assessment of the genetic basis of FHB resistance and DON accumulation in durum wheat. Accessions with multiple favorable alleles were identified and will be useful genetic resources to improve FHB resistance in durum breeding programs through marker-assisted recurrent selection and gene stacking.

Sequencing 4.3 million mutations in wheat promoters to understand and modify gene expression

Wheat is an important contributor to global food security, and further improvements are required to feed a growing human population. Functional genetics and genomics tools can help us to understand the function of different genes and to engineer beneficial changes. In this study, we used a promoter capture assay to sequence 2-kb regions upstream of all high-confidence annotated genes from 1,513 mutagenized plants from the tetraploid wheat variety Kronos. We identified 4.3 million induced mutations with an accuracy of 99.8%, resulting in a mutation density of 41.9 mutations per kb. We also remapped Kronos exome capture reads to Chinese Spring RefSeq v1.1, identified 4.7 million mutations, and predicted their effects on annotated genes. Using these predictions, we identified 59% more nonsynonymous substitutions and 49% more truncation mutations than in the original study. To show the biological value of the promoter dataset, we selected two mutations within the promoter of the VRN-A1 vernalization gene. Both mutations, located within transcription factor binding sites, significantly altered VRN-A1 expression, and one reduced the number of spikelets per spike. These publicly available sequenced mutant datasets provide rapid and inexpensive access to induced variation in the promoters and coding regions of most wheat genes. These mutations can be used to understand and modulate gene expression and phenotypes for both basic and commercial applications, where limited governmental regulations can facilitate deployment. These mutant collections, together with gene editing, provide valuable tools to accelerate functional genetic studies in this economically important crop.

TaSPL17s act redundantly with TaSPL14s to control spike development and their elite haplotypes may improve wheat grain yield

Wheat is a staple crop for the world's population, and there is constant pressure to improve grain yield, which is largely determined by plant architecture. SQUAMOSA promotor-binding protein-like (SPL) genes have been widely studied in rice, including their effects on plant architecture, grain development, and grain yield. However, the function of SPL homologous genes in wheat has not been well investigated. In this study, TaSPL14s and TaSPL17s, wheat's closest orthologous of OsSPL14, were functionally investigated using gene-editing assays, revealing that these genes redundantly influence plant height, tiller number, spike length, and thousand-grain weight (TGW). Bract outgrowth was frequently observed in the hexa-mutant, occasionally in the quintuple mutant but never in the wild type. Transcriptome analysis revealed that the expression of many spike development-associated genes was altered in taspl14taspl17 hexa-mutants compared to that in the wild type. In addition, we analyzed the sequence polymorphisms of TaSPL14s and TaSPL17s among wheat germplasm and found superior haplotypes of TaSPL17-A and TaSPL17-D with significantly higher TGW, which had been positively selected during wheat breeding. Accordingly, dCAPS and KASP markers were developed for TaSPL17-A and TaSPL17-D, respectively, providing a novel insight for molecular marker-assisted breeding in wheat. Overall, our results highlight the role of TaSPLs in regulating plant architecture and their potential application for wheat grain yield improvement through molecular breeding.

Molecular Cytological Analysis and Specific Marker Development in Wheat-Psathyrostachys huashanica Keng 3Ns Additional Line with Elongated Glume

Psathyrostachys huashanica Keng (2n = 2x = 14, NsNs) is an excellent gene resource for wheat breeding, which is characterized by early maturity, low plant height, and disease resistance. The wheat-P. huashanica derivatives were created by the elite genes of P. huashanica and permeate into common wheat through hybridization. Among them, a long-glume material 20JH1155 was identified, with larger grains and longer spike than its parents. In the present study, the methods of cytological observation, GISH, and sequential FISH analysis showed that 20JH1155 contained 21 pairs of wheat chromosomes and a pair of P. huashanica. There were some differences in 5A and 7B chromosomes between 20JH1155 and parental wheat 7182. Molecular marker, FISH, and sequence cloning indicated 20JH1155 alien chromosomes were 3Ns of P. huashanica. In addition, differentially expressed genes during immature spikelet development of 20JH1155 and 7182 and predicted transcription factors were obtained by transcriptome sequencing. Moreover, a total of 7 makers derived from Ph#3Ns were developed from transcriptome data. Taken together, the wheat-P. huashanica derived line 20JH1155 provides a new horizon on distant hybridization of wheat and accelerates the utilization of genes of P. huashanica.

Genome-edited TaTFL1-5 mutation decreases tiller and spikelet numbers in common wheat

Tillering is a critical agronomic trait of wheat (Triticum aestivum L.) that determines the shoot architecture and affects grain yield. TERMINAL FLOWER 1 (TFL1), encoding a phosphatidylethanolamine-binding protein, is implicated in the transition to flowering and shoot architecture in plant development. However, the roles of TFL1 homologs is little known in wheat development. CRISPR/Cas9-mediated targeted mutagenesis was used in this study to generate a set of wheat (Fielder) mutants with single, double or triple-null tatfl1-5 alleles. The wheat tatfl1-5 mutations decreased the tiller number per plant in the vegetative growth stage and the effective tiller number per plant and spikelet number per spike at maturity in the field. RNA-seq analysis showed that the expression of the auxin signaling-related and cytokinin signaling-related genes was significantly changed in the axillary buds of tatfl1-5 mutant seedlings. The results suggested that wheat TaTFL1-5s were implicated in tiller regulation by auxin and cytokinin signaling.

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.

Independent recruitment of FRUITFULL-like transcription factors in the convergent origins of vernalization-responsive grass flowering

Flowering in response to low temperatures (vernalization) has evolved multiple times independently across angiosperms as an adaptation to match reproductive development with the short growing season of temperate habitats. Despite the context of a generally conserved flowering time network, evidence suggests that the genes underlying vernalization responsiveness are distinct across major plant clades. Whether different or similar mechanisms underlie vernalization-induced flowering at narrower (e.g., family-level) phylogenetic scales is not well understood. To test the hypothesis that vernalization responsiveness has evolved convergently in temperate species of the grass family (Poaceae), we carried out flowering time experiments with and without vernalization in several representative species from different subfamilies. We then determined the likelihood that vernalization responsiveness evolved through parallel mechanisms by quantifying the response of Pooideae vernalization pathway FRUITFULL (FUL)-like genes to extended periods of cold. Our results demonstrate that vernalization-induced flowering has evolved multiple times independently in at least five grass subfamilies, and that different combinations of FUL-like genes have been recruited to this pathway on several occasions.

Tuning of meristem maturation rate increases yield in multiple Triticum aestivum cultivars

Breeding programs aim to improve crop yield and environmental stability for enhanced food security. The principal methodology in breeding for stable yield gain relies on the indirect selection of beneficial genetics by yield evaluation across diverse environmental conditions. This methodology requires substantial resources while delivering a slow pace of yield gain and environmental adaptation. Alternative methods are required to accelerate gain and adaptation, becoming even more imperative in a changing climate. New molecular tools and approaches can enable accelerated creation and deployment of multiple alleles of genes identified to control key traits. With the advent of tools that enable breeding by targeted allelic selection, identifying gene targets associated with an improved crop performance ideotype will become crucial. Previous studies have shown that altered photoperiod regimes increase yield in wheat (Triticum aestivum). In the current study, we have employed such treatments to study the resulting yield ideotype in five spring wheat cultivars. We found that the photoperiod treatment creates a yield ideotype arising from delayed spike establishment rates that are accompanied by increased early shoot expression of TARGET OF EAT1 (TaTOE1) genes. Genes identified in this way could be used for ideotype-based improve crop performance through targeted allele creation and selection in relevant environments.

Wheat genomic study for genetic improvement of traits in China

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

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.

Overexpression of the WAPO-A1 gene increases the number of spikelets per spike in bread wheat

Two homoeologous QTLs for number of spikelets per spike (SPS) were mapped on chromosomes 7AL and 7BL using two wheat MAGIC populations. Sets of lines contrasting for the QTL on 7AL were developed which allowed for the validation and fine mapping of the 7AL QTL and for the identification of a previously described candidate gene, WHEAT ORTHOLOG OF APO1 (WAPO1). Using transgenic overexpression in both a low and a high SPS line, we provide a functional validation for the role of this gene in determining SPS also in hexaploid wheat. We show that the expression levels of this gene positively correlate with SPS in multiple MAGIC founder lines under field conditions as well as in transgenic lines grown in the greenhouse. This work highlights the potential use of WAPO1 in hexaploid wheat for further yield increases. The impact of WAPO1 and SPS on yield depends on other genetic and environmental factors, hence, will require a finely balanced expression level to avoid the development of detrimental pleiotropic phenotypes.

Long-Amplicon Single-Molecule Sequencing Reveals Novel, Trait-Associated Variants of VERNALIZATION1 Homoeologs in Hexaploid Wheat

The gene VERNALIZATION1 (VRN1) is a key controller of vernalization requirement in wheat. The genome of hexaploid wheat (Triticum aestivum) harbors three homoeologous VRN1 loci on chromosomes 5A, 5B, and 5D. Structural sequence variants including small and large deletions and insertions and single nucleotide polymorphisms (SNPs) in the three homoeologous VRN1 genes not only play an important role in the control of vernalization requirement, but also have been reported to be associated with other yield related traits of wheat. Here we used single-molecule sequencing of barcoded long-amplicons to assay the full-length sequences (∼13 kbp plus 700 bp from the promoter sequence) of the three homoeologous VRN1 genes in a panel of 192 predominantly European winter wheat cultivars. Long read sequences revealed previously undetected duplications, insertions and single-nucleotide polymorphisms in the three homoeologous VRN1 genes. All the polymorphisms were confirmed by Sanger sequencing. Sequence analysis showed the predominance of the winter alleles vrn-A1, vrn-B1, and vrn-D1 across the investigated cultivars. Associations of SNPs and structural variations within the three VRN1 genes with 20 economically relevant traits including yield, nodal root-angle index and quality related traits were evaluated at the levels of alleles, haplotypes, and copy number variants. Cultivars carrying structural variants within VRN1 genes showed lower grain yield, protein yield and biomass compared to those with intact genes. Cultivars carrying a single vrn-A1 copy and a unique haplotype with a high number of SNPs were found to have elevated grain yield, kernels per spike and kernels per m2 along with lower grain sedimentation values. In addition, we detected a novel SNP polymorphism within the G-quadruplex region of the promoter of vrn-A1 that was associated with deeper roots in winter wheat. Our findings show that multiplex, single-molecule long-amplicon sequencing is a useful tool for detecting variants in target genes within large plant populations, and can be used to simultaneously assay sequence variants among target multiple gene homoeologs in polyploid crops. Numerous novel VRN1 haplotypes and alleles were identified that showed significantly associations to economically important traits. These polymorphisms were converted into PCR or KASP assays for use in marker-assisted breeding.

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.

The Triticum ispahanicum elongated glume locus P2 maps to chromosome 6A and is associated with the ectopic expression of SVP-A1

We propose the MADS-box transcription factor SVP-A1 as a promising candidate gene for the elongated glume locus P2, which maps to chromosome 6A instead of the previously proposed chromosome 7B. In rice and wheat, glume and floral organ length are positively correlated with grain size, making them an important target to increase grain size and potentially yield. The wheat subspecies Triticum ispahanicum is known to develop elongated glumes and floral organs as well as long grains. These multiple phenotypic effects are controlled by the P2 locus, which was previously mapped to wheat chromosome 7B. Using three mapping populations, we show that the long glume locus P2 does not map to chromosome 7B, but instead maps to a 1.68 Mbp interval on chromosome 6A. Within this interval, we identified SVP-A1, a MADS box transcription factor which is the direct ortholog of the maize gene underlying the 'pod corn' Tunicate locus and is a paralog to the T. polonicum elongated glume P1 gene. In T. ispahanicum, we identified a unique allele which has a 482-bp deletion in the SVP-A1 promoter and is associated with ectopic and higher expression of SVP-A1 in the elongated glumes and floral organs. We used near-isogenic lines (NILs) to show that P2 has a consistent positive effect on the length of glume, lemma, palea, spike and grain. Based on the mapping data, natural variation, biological function of SVP genes in cereals and expression analyses, we propose the MADS-box transcription factor SVP-A1 as a promising candidate for P2.

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

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

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.

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

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

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.