Transcriptome Association Identifies Regulators of Wheat Spike Architecture

Ploidy variation induces butterfly effect on chromatin topology in wheat

Polyploidy is a prominent driver of plant diversification, accompanied with dramatic chromosomal rearrangement and epigenetic changes that affect gene expression. How chromatin interactions within and between subgenomes adapt to ploidy transition remains poorly understood. We generate open chromatin interaction maps for natural hexaploid wheat (AABBDD), extracted tetraploid wheat (AABB), diploid wheat progenitor Aegilops tauschii (DD) and resynthesized hexaploid wheat (RHW, AABBDD). Thousands of intra- and interchromosomal loops are de novo established or disappeared in AB subgenomes after separation of D subgenome, in which 37-95% of novel loops are lost again in RHW after merger of D genome. Interestingly, more than half of novel loops are formed by cascade reactions that are triggered by disruption of chromatin interaction between AB and D subgenomes. The interaction repressed genes in RHW relative to DD are expression suppressed, resulting in more balanced expression of the three homoeologs in RHW. The interaction levels of cascade anchors are decreased step-by-step. Leading single nucleotide polymorphisms of yield- and plant architecture-related quantitative trait locus are significantly enriched in cascade anchors. The expression of 116 genes interacted with these anchors are significantly correlated with the corresponding traits. Our findings reveal trans-regulation of intrachromosomal loops by interchromosomal interactions during genome merger and separation in polyploid species.

TabHLH27 orchestrates root growth and drought tolerance to enhance water use efficiency in wheat

Cultivating high-yield wheat under limited water resources is crucial for sustainable agriculture in semiarid regions. Amid water scarcity, plants activate drought response signaling, yet the delicate balance between drought tolerance and development remains unclear. Through genome-wide association studies and transcriptome profiling, we identified a wheat atypical basic helix-loop-helix (bHLH) transcription factor (TF), TabHLH27-A1, as a promising quantitative trait locus candidate for both relative root dry weight and spikelet number per spike in wheat. TabHLH27-A1/B1/D1 knock-out reduced wheat drought tolerance, yield, and water use efficiency (WUE). TabHLH27-A1 exhibited rapid induction with polyethylene glycol (PEG) treatment, gradually declining over days. It activated stress response genes such as TaCBL8-B1 and TaCPI2-A1 while inhibiting root growth genes like TaSH15-B1 and TaWRKY70-B1 under short-term PEG stimulus. The distinct transcriptional regulation of TabHLH27-A1 involved diverse interacting factors such as TaABI3-D1 and TabZIP62-D1. Natural variations of TabHLH27-A1 influence its transcriptional responses to drought stress, with TabHLH27-A1Hap-II associated with stronger drought tolerance, larger root system, more spikelets, and higher WUE in wheat. Significantly, the excellent TabHLH27-A1Hap-II was selected during the breeding process in China, and introgression of TabHLH27-A1Hap-II allele improved drought tolerance and grain yield, especially under water-limited conditions. Our study highlights TabHLH27-A1's role in balancing root growth and drought tolerance, providing a genetic manipulation locus for enhancing WUE in wheat.

Recent advances in exploring transcriptional regulatory landscape of crops

Crop breeding entails developing and selecting plant varieties with improved agronomic traits. Modern molecular techniques, such as genome editing, enable more efficient manipulation of plant phenotype by altering the expression of particular regulatory or functional genes. Hence, it is essential to thoroughly comprehend the transcriptional regulatory mechanisms that underpin these traits. In the multi-omics era, a large amount of omics data has been generated for diverse crop species, including genomics, epigenomics, transcriptomics, proteomics, and single-cell omics. The abundant data resources and the emergence of advanced computational tools offer unprecedented opportunities for obtaining a holistic view and profound understanding of the regulatory processes linked to desirable traits. This review focuses on integrated network approaches that utilize multi-omics data to investigate gene expression regulation. Various types of regulatory networks and their inference methods are discussed, focusing on recent advancements in crop plants. The integration of multi-omics data has been proven to be crucial for the construction of high-confidence regulatory networks. With the refinement of these methodologies, they will significantly enhance crop breeding efforts and contribute to global food security.

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.

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.

HOMEOBOX2, the paralog of SIX-ROWED SPIKE1/HOMEOBOX1, is dispensable for barley spikelet development

The HD-ZIP class I transcription factor Homeobox 1 (HvHOX1), also known as Vulgare Row-type Spike 1 (VRS1) or Six-rowed Spike 1, regulates lateral spikelet fertility in barley (Hordeum vulgare L.). It was shown that HvHOX1 has a high expression only in lateral spikelets, while its paralog HvHOX2 was found to be expressed in different plant organs. Yet, the mechanistic functions of HvHOX1 and HvHOX2 during spikelet development are still fragmentary. Here, we show that compared with HvHOX1, HvHOX2 is more highly conserved across different barley genotypes and Hordeum species, hinting at a possibly vital but still unclarified biological role. Using bimolecular fluorescence complementation, DNA-binding, and transactivation assays, we validate that HvHOX1 and HvHOX2 are bona fide transcriptional activators that may potentially heterodimerize. Accordingly, both genes exhibit similar spatiotemporal expression patterns during spike development and growth, albeit their mRNA levels differ quantitatively. We show that HvHOX1 delays the lateral spikelet meristem differentiation and affects fertility by aborting the reproductive organs. Interestingly, the ancestral relationship of the two genes inferred from their co-expressed gene networks suggested that HvHOX1 and HvHOX2 might play a similar role during barley spikelet development. However, CRISPR-derived mutants of HvHOX1 and HvHOX2 demonstrated the suppressive role of HvHOX1 on lateral spikelets, while the loss of HvHOX2 does not influence spikelet development. Collectively, our study shows that through the suppression of reproductive organs, lateral spikelet fertility is regulated by HvHOX1, whereas HvHOX2 is dispensable for spikelet development in barley.

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.

TaLBD41 interacts with TaNAC2 to regulate nitrogen uptake and metabolism in response to nitrate availability

Nitrate is the main source of nitrogen (N) available to plants and also is a signal that triggers complex regulation of transcriptional networks to modulate a wide variety of physiological and developmental responses in plants. How plants adapt to soil nitrate fluctuations is a complex process involving a fine-tuned response to nitrate provision and N starvation, the molecular mechanisms of which remain largely uncharted. Here, we report that the wheat transcription factor TaLBD41 interacts with the nitrate-inducible transcription factor TaNAC2 and is repressed by nitrate provision. Electrophoretic mobility shift assay and dual-luciferase system show that the TaLBD41-NAC2 interaction confers homeostatic coordination of nitrate uptake, reduction, and assimilation by competitively binding to TaNRT2.1, TaNR1.2, and TaNADH-GOGAT. Knockdown of TaLBD41 expression enhances N uptake and assimilation, increases spike number, grain yield, and nitrogen harvest index under different N supply conditions. We also identified an elite haplotype of TaLBD41-2B associated with increased spike number and grain yield. Our study uncovers a novel mechanism underlying the interaction between two transcription factors in mediating wheat adaptation to nitrate availability by antagonistically regulating nitrate uptake and assimilation, providing a potential target for designing varieties with efficient N use in wheat (Triticum aestivum).

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.

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.

Deciphering the evolution and complexity of wheat germplasm from a genomic perspective

Bread wheat provides an essential fraction of the daily calorific intake for humanity. Due to its huge and complex genome, progress in studying on the wheat genome is substantially trailed behind those of the other two major crops, rice and maize, for at least a decade. With rapid advances in genome assembling and reduced cost of high-throughput sequencing, emerging de novo genome assemblies of wheat and whole-genome sequencing data are leading to a paradigm shift in wheat research. Here, we review recent progress in dissecting the complex genome and germplasm evolution of wheat since the release of the first high-quality wheat genome. New insights have been gained in the evolution of wheat germplasm during domestication and modern breeding progress, genomic variations at multiple scales contributing to the diversity of wheat germplasm, and complex transcriptional and epigenetic regulations of functional genes in polyploid wheat. Genomics databases and bioinformatics tools meeting the urgent needs of wheat genomics research are also summarized. The ever-increasing omics data, along with advanced tools and well-structured databases, are expected to accelerate deciphering the germplasm and gene resources in wheat for future breeding advances.

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.

Dynamic gene regulatory networks improving spike fertility through regulation of floret primordia fate in wheat

The developmental process of spike is critical for spike fertility through affecting floret primordia fate in wheat; however, the genetic regulation of this dynamic and complex developmental process remains unclear. Here, we conducted a high temporal-resolution analysis of spike transcriptomes and monitored the number and morphology of floret primordia within spike. The development of all floret primordia in a spike was clearly separated into three distinct phases: differentiation, pre-dimorphism and dimorphism. Notably, we identified that floret primordia with meiosis ability at the pre-dimorphism phase usually develop into fertile floret primordia in the next dimorphism phase. Compared to control, increasing plant space treatment achieved the maximum increasement range (i.e., 50%) in number of fertile florets by accelerating spike development. The process of spike fertility improvement was directed by a continuous and dynamic regulatory network involved in transcription factor and genes interaction. This was based on the coordination of genes related to heat shock protein and jasmonic acid biosynthesis during differentiation phase, and genes related to lignin, anthocyanin and chlorophyll biosynthesis during dimorphism phase. The multi-dimensional association with high temporal-resolution approach reported here allows rapid identification of genetic resource for future breeding studies to realise the maximum spike fertility potential in more cereal crops.

Cell signaling in the shoot apical meristem

Distinct from animals, plants maintain organogenesis from specialized tissues termed meristems throughout life. In the shoot apex, the shoot apical meristem (SAM) produces all aerial organs, such as leaves, from its periphery. For this, the SAM needs to precisely balance stem cell renewal and differentiation, which is achieved through dynamic zonation of the SAM, and cell signaling within functional domains is key for SAM functions. The WUSCHEL-CLAVATA feedback loop plays a key role in SAM homeostasis, and recent studies have uncovered new components, expanding our understanding of the spatial expression and signaling mechanism. Advances in polar auxin transport and signaling have contributed to knowledge of the multifaceted roles of auxin in the SAM and organogenesis. Finally, single-cell techniques have expanded our understanding of the cellular functions within the shoot apex at single-cell resolution. In this review, we summarize the most up-to-date understanding of cell signaling in the SAM and focus on the multiple levels of regulation of SAM formation and maintenance.

An intron-located single nucleotide variation of TaGS5-3D is related to wheat grain size through accumulating intron retention transcripts

KEY MESSAGE

Thirty-three stable QTL for 13 yield-related traits across ten environments were identified in the PD34/MY47 RIL population, and a candidate gene TaGS5-3D in Qmt.nwafu.3D was preliminarily identified to affect grain-related traits through accumulation of specific transcripts. Dissecting the genetic basis of yield-related traits is pivotal for improvement of wheat yield potential. In this study, a recombinant inbred line (RIL) population genotyped by SNP markers was used to detect quantitative trait loci (QTL) related to yield-related traits in ten environments. A total of 102 QTL were detected, including 33 environmentally stable QTL and 69 putative QTL. Among them, Qmt.nwafu.3D was identified as a pleiotropic QTL in the physical interval of 149.77-154.11 Mb containing a potential candidate gene TaGS5-3D. An SNP (T > C) was detected in its ninth intron, and TaGS5-3D-C was validated as a superior allele associated with larger grains using a CAPS marker. Interestingly, we found that TaGS5-3D-C was closely related to significantly up-regulated expression of intron-retained transcript (TaGS5-3D-PD34.1), while TaGS5-3D-T was related to dominant expression of normal splicing transcript (TaGS5-3D-MY47.1). Our results indicated that alternative splicing associated with the SNP T/C could be involved in the regulation of grain-related traits, laying a foundation for the functional analysis of TaGS5-3D and its greater potential application in high-yield wheat breeding.

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.

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.

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

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

Genome-wide identification and expression analysis disclose the pivotal PHOSPHATIDYLETHANOLAMINE BINDING PROTEIN members that may be utilized for yield improvement of Chenopodium quinoa

Quinoa (Chenopodium quinoa) is a prospective orphan crop that needs yield improvement. Previous studies indicate PHOSPHATIDYLETHANOLAMINE BINDING PROTEIN (PEBP) family genes are highly associated with the key agronomic traits of crops. Characterizing the pivotal PEBP genes will speed up the domestication and yield improvement of quinoa. Previous investigations on PEBP genes of Chenopodium species indicated that, the PEBP genes, despite in the same subclade, may have experienced functional diversification. Especially, the allotetraploidy (AABB) and numerous segmental duplications and chromosomal rearrangements in quinoa make it more difficult to understand the functions of PEBP genes. More recently, 6 quinoa FT subfamily genes were predicted to be related to flowering of quinoa. However, investigation on the whole PEBP family members is still lacking. In this study, we obtained 23 PEBP genes, including 5 MFT, 11 FTL and 7 TFL genes. We found 7 orthologous gene pairs, from sub-genome A and sub-genome B, respectively, showing collinearities with sugar beet. Evolution analysis on PEBP genes of two quinoa sub-genomes, sugar beet and relatives of diploid ancestors indicated that, the reasons for gene duplication events varied and 4 tandem duplications are the major reason for PEBP family expansion. Tissue-specific expression analysis suggested that expression patterns are mostly differing between orthologous gene pairs. Analysis on gene expressions at 6 stages suggested the possible positive roles of CqFTL1/CqFTL2, CqFTL5, CqFTL8, CqFTL6/CqFTL9 and CqTFL6/CqTFL7, and negative roles of CqTFL1/CqTFL2/CqTFL3, CqTFL4/CqTFL5 in inflorescence branching. Expression analysis in ABA-treated seed, in combination with the cis-acting element distribution analysis, indicated that CqMFT2, CqMFT3 and CqMFT4 may regulate seed germination via ABA signaling. Observations on responses to night break and photoperiod changes highlighted the roles of CqFTL5 and CqFTL8 under short day, and CqFTL6 under long day for quinoa flowering. Further, co-expression network analysis indicated that 64 transcription factors may act upstream of CqFTL5 and CqFTL8 to regulate flowering. Together, this study will help us identify the pivotal PEBP genes that may be utilized for quinoa breeding in future.

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

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

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

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

Tracing the Evolution of the SEPALLATA Subfamily across Angiosperms Associated with Neo- and Sub-Functionalization for Reproductive and Agronomically Relevant Traits

SEPALLATA transcription factors (SEP TFs) have been extensively studied in angiosperms as pivotal components of virtually all the MADS-box tetrameric complex master regulators of floral organ identities. However, there are published reports that suggest that some SEP members also regulate earlier reproductive events, such as inflorescence meristem determinacy and inflorescence architecture, with potential for application in breeding programs in crops. The SEP subfamily underwent a quite complex pattern of duplications during the radiation of the angiosperms. Taking advantage of the many whole genomic sequences now available, we present a revised and expanded SEP phylogeny and link it to the known functions of previously characterized genes. This snapshot supports the evidence that the major SEP3 clade is highly specialized for the specification of the three innermost floral whorls, while its sister LOFSEP clade is functionally more versatile and has been recruited for diverse roles, such as the regulation of extra-floral bract formation and inflorescence determinacy and shape. This larger pool of angiosperm SEP genes confirms previous evidence that their evolution was driven by whole-genome duplications rather than small-scale duplication events. Our work may help to identify those SEP lineages that are the best candidates for the improvement of inflorescence traits, even in far distantly related crops.

Transcriptome analysis and mining of genes related to shade tolerance in foxtail millet (Setaria italica (L.) P. Beauv.)

A stereo interplanting system with foxtail millet beneath chestnut trees is an effective planting method to raise the utilization of land in chestnut orchards, increase yields and improve quality of chestnut nuts. Consequently, exploration of genes involved in shade tolerance response in foxtail millet and breeding shade-tolerant varieties have become urgent issues. In this study, RNA-seq of leaf samples from two shade-tolerant varieties and three shade-intolerant varieties of foxtail millet at the booting stage was performed. Comparisons between the varieties revealed that 70 genes were commonly differentially expressed. Moreover, the ratio of net photosynthetic rate under shaded environment to that under light environment could be used as an indicator of shade tolerance. Subsequently, weighted gene co-expression network analysis was employed to construct a co-expression network and modules were correlated with this ratio. A total of 375 genes were identified as potentially relevant to shade tolerance, among which nine genes were also present in the 70 differentially expressed genes, which implied that they were good candidates for genes involved in shade tolerance. Our results provide valuable resources for elucidation of the molecular mechanisms underlying shade tolerance and will contribute to breeding of shade-tolerant foxtail millet that are adapted to the shaded environment under chestnut trees.

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.

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.

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.

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.

Combining QTL mapping and gene co-expression network analysis for prediction of candidate genes and molecular network related to yield in wheat

BACKGROUND

Wheat (Triticum aestivum L.) is an important cereal crop. Increasing grain yield for wheat is always a priority. Due to the complex genome of hexaploid wheat with 21 chromosomes, it is difficult to identify underlying genes by traditional genetic approach. The combination of genetics and omics analysis has displayed the powerful capability to identify candidate genes for major quantitative trait loci (QTLs), but such studies have rarely been carried out in wheat. In this study, candidate genes related to yield were predicted by a combined use of linkage mapping and weighted gene co-expression network analysis (WGCNA) in a recombinant inbred line population.

RESULTS

QTL mapping was performed for plant height (PH), spike length (SL) and seed traits. A total of 68 QTLs were identified for them, among which, 12 QTLs were stably identified across different environments. Using RNA sequencing, we scanned the 99,168 genes expression patterns of the whole spike for the recombinant inbred line population. By the combined use of QTL mapping and WGCNA, 29, 47, 20, 26, 54, 46 and 22 candidate genes were predicted for PH, SL, kernel length (KL), kernel width, thousand kernel weight, seed dormancy, and seed vigor, respectively. Candidate genes for different traits had distinct preferences. The known PH regulation genes Rht-B and Rht-D, and the known seed dormancy regulation genes TaMFT can be selected as candidate gene. Moreover, further experiment revealed that there was a SL regulatory QTL located in an interval of about 7 Mbp on chromosome 7A, named TaSL1, which also involved in the regulation of KL.

CONCLUSIONS

A combination of QTL mapping and WGCNA was applied to predicted wheat candidate genes for PH, SL and seed traits. This strategy will facilitate the identification of candidate genes for related QTLs in wheat. In addition, the QTL TaSL1 that had multi-effect regulation of KL and SL was identified, which can be used for wheat improvement. These results provided valuable molecular marker and gene information for fine mapping and cloning of the yield-related trait loci in the future.

A wiring diagram to integrate physiological traits of wheat yield potential

As crop yields are pushed closer to biophysical limits, achieving yield gains becomes increasingly challenging and will require more insight into deterministic pathways to yields. Here, we propose a wiring diagram as a platform to illustrate the interrelationships of the physiological traits that impact wheat yield potential and to serve as a decision support tool for crop scientists. The wiring diagram is based on the premise that crop yield is a function of photosynthesis (source), the investment of assimilates into reproductive organs (sinks) and the underlying processes that enable expression of both. By illustrating these linkages as coded wires, the wiring diagram can show connections among traits that may not have been apparent, and can inform new research hypotheses and guide crosses designed to accumulate beneficial traits and alleles in breeding. The wiring diagram can also serve to create an ever-richer common point of reference for refining crop models in the future.

Genomic variants affecting homoeologous gene expression dosage contribute to agronomic trait variation in allopolyploid wheat

Allopolyploidy greatly expands the range of possible regulatory interactions among functionally redundant homoeologous genes. However, connection between the emerging regulatory complexity and expression and phenotypic diversity in polyploid crops remains elusive. Here, we use diverse wheat accessions to map expression quantitative trait loci (eQTL) and evaluate their effects on the population-scale variation in homoeolog expression dosage. The relative contribution of cis- and trans-eQTL to homoeolog expression variation is strongly affected by both selection and demographic events. Though trans-acting effects play major role in expression regulation, the expression dosage of homoeologs is largely influenced by cis-acting variants, which appear to be subjected to selection. The frequency and expression of homoeologous gene alleles showing strong expression dosage bias are predictive of variation in yield-related traits, and have likely been impacted by breeding for increased productivity. Our study highlights the importance of genomic variants affecting homoeolog expression dosage in shaping agronomic phenotypes and points at their potential utility for improving yield in polyploid crops.

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

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

A transcriptome-based association study of growth, wood quality, and oleoresin traits in a slash pine breeding population

Slash pine (Pinus elliottii Engelm.) is an important timber and resin species in the United States, China, Brazil and other countries. Understanding the genetic basis of these traits will accelerate its breeding progress. We carried out a genome-wide association study (GWAS), transcriptome-wide association study (TWAS) and weighted gene co-expression network analysis (WGCNA) for growth, wood quality, and oleoresin traits using 240 unrelated individuals from a Chinese slash pine breeding population. We developed high quality 53,229 single nucleotide polymorphisms (SNPs). Our analysis reveals three main results: (1) the Chinese breeding population can be divided into three genetic groups with a mean inbreeding coefficient of 0.137; (2) 32 SNPs significantly were associated with growth and oleoresin traits, accounting for the phenotypic variance ranging from 12.3% to 21.8% and from 10.6% to 16.7%, respectively; and (3) six genes encoding PeTLP, PeAP2/ERF, PePUP9, PeSLP, PeHSP, and PeOCT1 proteins were identified and validated by quantitative real time polymerase chain reaction for their association with growth and oleoresin traits. These results could be useful for tree breeding and functional studies in advanced slash pine breeding program.

Slash pine is an important source of original timber and resin production on commercial forest plantations. It is necessary to implement precise breeding strategies to improve timber quality and resin yield. However, little is known about the species’ molecular genetic basis. Using a transcriptome dataset with sequencing from 240 individuals in the slash pine population, we combined multiple approaches (based on gene variation, expression variation and co-expression network) to dissect the genetic structure for slash pine major breeding traits. We found that the research population could be divided into three genetic groups with a mean heterozygosity of 0.2246. We also found that six genes with important functions in slash pine resin synthesis and timber formation through association studies. Four new SNPs associatation with the average ring width were also discovered. Our results provide new insights into the molecular genetic basis of important traits in slash pine and provide a comprehensive method for association analyses of conifer tree species with large genome.

Rice SEPALLATA genes OsMADS5 and OsMADS34 cooperate to limit inflorescence branching by repressing the TERMINAL FLOWER1-like gene RCN4

The spatiotemporal control of meristem identity is critical for determining inflorescence architecture, and thus yield, of cereal plants. However, the precise mechanisms underlying inflorescence and spikelet meristem determinacy in cereals are still largely unclear. We have generated loss-of-function and overexpression mutants of the paralogous OsMADS5 and OsMADS34 genes in rice (Oryza sativa), and analysed their panicle phenotypes. Using chromatin immunoprecipitation, electrophoretic mobility-shift and dual-luciferase assays, we have also identified RICE CENTRORADIALIS 4 (RCN4), a TFL1-like gene, as a direct downstream target of both OsMADS proteins, and have analysed RCN4 mutants. The osmads5 osmads34 mutant lines had significantly enhanced panicle branching with increased secondary, and even tertiary and quaternary, branches, compared to wild-type (WT) and osmads34 plants. The osmads34 mutant phenotype could largely be rescued by also knocking out RCN4. Moreover, transgenic panicles overexpressing RCN4 had significantly increased branching, and initiated development of c. 7× more spikelets than WT. Our results reveal a role for OsMADS5 in panicle development, and show that OsMADS5 and OsMADS34 play similar functions in limiting branching and promoting the transition to spikelet meristem identity, in part by repressing RCN4 expression. These findings provide new insights to better understand the molecular regulation of rice inflorescence architecture.

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

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

The wheat AGL6-like MADS-box gene is a master regulator for floral organ identity and a target for spikelet meristem development manipulation

The AGAMOUS-LIKE6 (AGL6)-like genes are ancient MADS-box genes and are functionally studied in a few model plants. The knowledge of these genes in wheat remains limited. Here, by studying a 'double homoeolog mutant' of the AGL6 gene in tetraploid wheat, we showed that AGL6 was required for the development of all four whorls of floral organs with dosage-dependent effect on floret fertility. Yeast two-hybrid analyses detected interactions of AGL6 with all classes of MADS-box proteins in the ABCDE model for floral organ development. AGL6 was found to interact with several additional proteins, including the G protein β and γ (DEP1) subunits. Analysis of the DEP1-B mutant showed a significant reduction in spikelet number per spike in tetraploid wheat, while overexpression of AGL6 in common wheat increased the spikelet number per spike and hence the grain number per spike. RNA-seq analysis identified the regulation of several meristem activity genes by AGL6, such as FUL2 and TaMADS55. Our work therefore extensively updated the wheat ABCDE model and proposed an alternative approach to improve wheat grain yield by manipulating the AGL6 gene.

Interactions between SQUAMOSA and SHORT VEGETATIVE PHASE MADS-box proteins regulate meristem transitions during wheat spike development

Inflorescence architecture is an important determinant of crop productivity. The number of spikelets produced by the wheat inflorescence meristem (IM) before its transition to a terminal spikelet (TS) influences the maximum number of grains per spike. Wheat MADS-box genes VERNALIZATION 1 (VRN1) and FRUITFULL 2 (FUL2) (in the SQUAMOSA-clade) are essential to promote the transition from IM to TS and for spikelet development. Here we show that SQUAMOSA genes contribute to spikelet identity by repressing MADS-box genes VEGETATIVE TO REPRODUCTIVE TRANSITION 2 (VRT2), SHORT VEGETATIVE PHASE 1 (SVP1), and SVP3 in the SVP clade. Constitutive expression of VRT2 resulted in leafy glumes and lemmas, reversion of spikelets to spikes, and downregulation of MADS-box genes involved in floret development, whereas the vrt2 mutant reduced vegetative characteristics in spikelets of squamosa mutants. Interestingly, the vrt2 svp1 mutant showed similar phenotypes to squamosa mutants regarding heading time, plant height, and spikelets per spike, but it exhibited unusual axillary inflorescences in the elongating stem. We propose that SQUAMOSA-SVP interactions are important to promote heading, formation of the TS, and stem elongation during the early reproductive phase, and that downregulation of SVP genes is then necessary for normal spikelet and floral development. Manipulating SVP and SQUAMOSA genes can contribute to engineering spike architectures with improved productivity.

Reassessing the evolution of the 1-deoxy-D-xylulose 5-phosphate synthase family suggests a possible novel function for the DXS class 3 proteins

The methylerythritol 4-phosphate (MEP) pathway is of paramount importance for generating plastidial isoprenoids. The first enzyme of the MEP pathway, 1-deoxy-D-xylulose-5-phosphate synthase (DXS), catalyzes a flux-controlling step. In plants the DXS gene family is composed of three distinct classes with non-redundant functions. Although the DXS1 and DXS2 subfamilies have been well characterized, the DXS3 subfamily has been considerably understudied. Here, we carried out in silico and functional analyses to better understand the DXS3 class. Our phylogenetic analysis showed high variation in copy number among the different DXS classes, with the apparent absence of DXS1 class in some species. We found that DXS3 subfamily emerged later than DXS1 and DXS2 and it is under less intense purifying selection. Furthermore, in the DXS3 subfamily critical amino acids positions in the thiamine pyrophosphate binding pocket are not conserved. We demonstrated that the DXS3 proteins from Arabidopsis, Maize, and Rice lack functional DXS activity. Moreover, the Arabidopsis DXS3 protein displayed distinctive sub-organellar chloroplast localization not observed in any DXS1 or DXS2 proteins. Co-expression analysis of the DXS3 from Arabidopsis showed that, unlike DXS1 and DXS2 proteins, it co-expresses with genes related to post-embryonic development and reproduction and not with primary metabolism and isoprenoid synthesis.

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

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

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

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

Interactive Effects of Lanthanum and Calcium on Cadmium Accumulation in Wheat with Special Reference to TaNramp5 Expression Regulated by Calmodulin

Lanthanum (La), calcium (Ca), and cadmium (Cd) have similar physical and chemical properties because of their similar ionic radius. Although the interactions between La, Ca, and Cd have been frequently reported in plants, few studies have investigated the interactive effects of La and Ca on the growth and Cd accumulation in plants. Therefore, we investigated the interactive effects of La and Ca on the growth and Cd accumulation in wheat under Cd exposure by a hydroponic experiment. The results indicated that wheat growth was significantly affected by La-Cd and La-Ca interactions. The accumulation of Cd in wheat was significantly affected by La-Ca and La-Cd interactions and La-Ca-Cd interplay. Correlation analysis indicated that Ca deficiency stimulated La to promote wheat growth and mitigate Cd toxicity. Simultaneously, a low Ca supply stimulated La to decrease Cd accumulation in wheat and induce TaNramp5 expression. In addition, Cd accumulation in wheat was significantly affected by the W7-La interaction and W7-La-Ca interplay. All of the results suggested that La, Ca, and Cd probably share the same binding sites in calmodulin (TaCaM) and La could affect Cd accumulation in wheat by interacting with TaCaM and then downregulating the expression of TaNramp5.

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.

Using natural variation to understand plant responses to iron availability

Iron bioavailability varies dramatically between soil types across the globe. This has given rise to high levels of natural variation in plant iron responses, allowing members of even a single species to thrive across a wide range of soil types. In recent years we have seen the use of genome-wide association analysis to identify natural variants underlying plant responses to changes in iron availability in both Arabidopsis and important crop species. These studies have provided insights into which genes have been important in shaping local adaptation to iron availability in different plant species and have allowed the discovery of novel regulators and mechanisms, not previously identified using mutagenesis approaches. Furthermore, these studies have allowed the identification of markers that can be used to accelerate breeding of future elite varieties with increased resilience to iron stress and improved nutritional quality. The studies highlighted here show that, in addition to studying plant responses to iron alone, it is important to consider these responses within the context of plant nutrition more broadly and to also consider iron regulation in relation to additional traits of agronomic importance such as yield and disease resistance.

Thermal Stresses in Maize: Effects and Management Strategies

Climate change can decrease the global maize productivity and grain quality. Maize crop requires an optimal temperature for better harvest productivity. A suboptimal temperature at any critical stage for a prolonged duration can negatively affect the growth and yield formation processes. This review discusses the negative impact of temperature extremes (high and low temperatures) on the morpho-physiological, biochemical, and nutritional traits of the maize crop. High temperature stress limits pollen viability and silks receptivity, leading to a significant reduction in seed setting and grain yield. Likewise, severe alterations in growth rate, photosynthesis, dry matter accumulation, cellular membranes, and antioxidant enzyme activities under low temperature collectively limit maize productivity. We also discussed various strategies with practical examples to cope with temperature stresses, including cultural practices, exogenous protectants, breeding climate-smart crops, and molecular genomics approaches. We reviewed that identified quantitative trait loci (QTLs) and genes controlling high- and low temperature stress tolerance in maize could be introgressed into otherwise elite cultivars to develop stress-tolerant cultivars. Genome editing has become a key tool for developing climate-resilient crops. Moreover, challenges to maize crop improvement such as lack of adequate resources for breeding in poor countries, poor communication among the scientists of developing and developed countries, problems in germplasm exchange, and high cost of advanced high-throughput phenotyping systems are discussed. In the end, future perspectives for maize improvement are discussed, which briefly include new breeding technologies such as transgene-free clustered regularly interspaced short palindromic repeat (CRISPR)/CRISPR-associated (Cas)-mediated genome editing for thermo-stress tolerance in maize.

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.

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

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

Coevolution in Hybrid Genomes: Nuclear-Encoded Rubisco Small Subunits and Their Plastid-Targeting Translocons Accompanying Sequential Allopolyploidy Events in Triticum

The Triticum/Aegilops complex includes hybrid species resulting from homoploid hybrid speciation and allopolyploid speciation. Sequential allotetra- and allohexaploidy events presumably result in two challenges for the hybrids, which involve 1) cytonuclear stoichiometric disruptions caused by combining two diverged nuclear genomes with the maternal inheritance of the cytoplasmic organellar donor; and 2) incompatibility of chimeric protein complexes with diverged subunits from nuclear and cytoplasmic genomes. Here, we describe coevolution of nuclear rbcS genes encoding the small subunits of Rubisco (ribulose 1,5-bisphosphate carboxylase/oxygenase) and nuclear genes encoding plastid translocons, which mediate recognition and translocation of nuclear-encoded proteins into plastids, in allopolyploid wheat species. We demonstrate that intergenomic paternal-to-maternal gene conversion specifically occurred in the genic region of the homoeologous rbcS3 gene from the D-genome progenitor of wheat (abbreviated as rbcS3D) such that it encodes a maternal-like or B-subgenome-like SSU3D transit peptide in allohexaploid wheat but not in allotetraploid wheat. Divergent and limited interaction between SSU3D and the D-subgenomic TOC90D translocon subunit is implicated to underpin SSU3D targeting into the chloroplast of hexaploid wheat. This implicates early selection favoring individuals harboring optimal maternal-like organellar SSU3D targeting in hexaploid wheat. These data represent a novel dimension of cytonuclear evolution mediated by organellar targeting and transportation of nuclear proteins.

Deficiencies in the formation and regulation of anther cuticle and tryphine contribute to male sterility in cotton PGMS line

BACKGROUND

Male sterility is a simple and efficient pollination control system that is widely exploited in hybrid breeding. In upland cotton, CCRI9106, a photosensitive genetic male sterile (PGMS) mutant isolated from CCRI040029, was reported of great advantages to cotton heterosis. However, little information concerning the male sterility of CCRI9106 is known. Here, comparative transcriptome analysis of CCRI9106 (the mutant, MT) and CCRI040029 (the wild type, WT) anthers in Anyang (long-day, male sterile condition to CCRI9106) was performed to reveal the potential male sterile mechanism of CCRI9106.

RESULTS

Light and electron microscopy revealed that the male sterility phenotype of MT was mainly attributed to irregularly exine, lacking tryphine and immature anther cuticle. Based on the cytological characteristics of MT anthers, anther RNA libraries (18 in total) of tetrad (TTP), late uninucleate (lUNP) and binucleate (BNP) stages in MT and WT were constructed for transcriptomic analysis, therefore revealing a total of 870,4 differentially expressed genes (DEGs). By performing gene expression pattern analysis and protein-protein interaction (PPI) networks construction, we found down-regulation of DEGs, which enriched by the lipid biosynthetic process and the synthesis pathways of several types of secondary metabolites such as terpenoids, flavonoids and steroids, may crucial to the male sterility phenotype of MT, and resulting in the defects of anther cuticle and tryphine, even the irregularly exine. Furthermore, several lipid-related genes together with ABA-related genes and MYB transcription factors were identified as hub genes via weighted gene co-expression network analysis (WGCNA). Additionally, the ABA content of MT anthers was reduced across all stages when compared with WT anthers. At last, genes related to the formation of anther cuticle and tryphine could activated in MT under short-day condition.

CONCLUSIONS

We propose that the down-regulation of genes related to the assembly of anther cuticle and tryphine may lead to the male sterile phenotype of MT, and MYB transcription factors together with ABA played key regulatory roles in these processes. The conversion of fertility in different photoperiods may closely relate to the functional expression of these genes. These findings contribute to elucidate the mechanism of male sterility in upland cotton.

Deficiencies in the formation and regulation of anther cuticle and tryphine contribute to male sterility in cotton PGMS line

Background

Male sterility is a simple and efficient pollination control system that is widely exploited in hybrid breeding. In upland cotton, CCRI9106, a photosensitive genetic male sterile (PGMS) mutant isolated from CCRI040029, was reported of great advantages to cotton heterosis. However, little information concerning the male sterility of CCRI9106 is known. Here, comparative transcriptome analysis of CCRI9106 (the mutant, MT) and CCRI040029 (the wild type, WT) anthers in Anyang (long-day, male sterile condition to CCRI9106) was performed to reveal the potential male sterile mechanism of CCRI9106.

Results

Light and electron microscopy revealed that the male sterility phenotype of MT was mainly attributed to irregularly exine, lacking tryphine and immature anther cuticle. Based on the cytological characteristics of MT anthers, anther RNA libraries (18 in total) of tetrad (TTP), late uninucleate (lUNP) and binucleate (BNP) stages in MT and WT were constructed for transcriptomic analysis, therefore revealing a total of 870,4 differentially expressed genes (DEGs). By performing gene expression pattern analysis and protein-protein interaction (PPI) networks construction, we found down-regulation of DEGs, which enriched by the lipid biosynthetic process and the synthesis pathways of several types of secondary metabolites such as terpenoids, flavonoids and steroids, may crucial to the male sterility phenotype of MT, and resulting in the defects of anther cuticle and tryphine, even the irregularly exine. Furthermore, several lipid-related genes together with ABA-related genes and MYB transcription factors were identified as hub genes via weighted gene co-expression network analysis (WGCNA). Additionally, the ABA content of MT anthers was reduced across all stages when compared with WT anthers. At last, genes related to the formation of anther cuticle and tryphine could activated in MT under short-day condition.

Conclusions

We propose that the down-regulation of genes related to the assembly of anther cuticle and tryphine may lead to the male sterile phenotype of MT, and MYB transcription factors together with ABA played key regulatory roles in these processes. The conversion of fertility in different photoperiods may closely relate to the functional expression of these genes. These findings contribute to elucidate the mechanism of male sterility in upland cotton.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-020-07250-1.