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

Cloning and functional analysis of ZmMADS42 gene in maize

Maize (Zea mays L.) is the most important cereal crop in the world. Flowering period and photoperiod play important roles in the reproductive development of maize. This study, investigated ZmMADS42, a gene that is highly expressed in the shoot apical meristem. Agrobacterium infection was used to successfully obtain overexpressed ZmMADS42 plants. Fluorescence quantitative PCR revealed that the expression of the ZmMADS42 gene in the shoot apical meristem of transgenic plants was 2.8 times higher than that of the wild-type(WT). In addition, the expression of the ZmMADS42 gene in the endosperm was 2.4 times higher than that in the wild-type. The seed width of the T2 generation increased by 5.35%, whereas the seed length decreased by 7.78% compared with that of the wild-type. Dissection of the shoot tips of transgenic and wild-type plants from the 7-leaf stage to the 9-leaf stage revealed that the transgenic plants entered the differentiation stage earlier and exhibited more tassel meristems during their vegetative growth period. The mature transgenic plants were approximately 20 cm shorter in height and had a lower panicle position than the wild-type plants. Comparing the flowering period, the tasseling, powdering, and silking stages of the transgenic plants occurred 10 days earlier than those of the wild-type plants. The results showed that the ZmMADS42 gene played a significant role in regulating the flowering period and plant height of maize.

Characterization and validation of TaAGL66, a gene related to fertility conversion of wheat in the presence of Aegilops kotschyi cytoplasm

MAIN CONCLUSION

TaAGL66, a MADS-box transcription factor highly expressed in fertile anthers of KTM3315A, regulates anther and/or pollen development, as well as male fertility in wheat with Aegilops kotschyi cytoplasm. Male sterility, as a string of sophisticated biological processes in higher plants, is commonly regulated by transcription factors (TFs). Among them, MADS-box TFs are mainly participated in the processes of floral organ formation and pollen development, which are tightly related to male sterility, but they have been little studied in the reproductive development in wheat. In our study, TaAGL66, a gene that was specifically expressed in spikes and highly expressed in fertile anthers, was identified by RNA sequencing and the expression profiles data of these genes, and qRT-PCR analyses, which was localized to the nucleus. Silencing of TaAGL66 under fertility condition in KTM3315A, a thermo-sensitive male sterile line with Ae. kotschyi cytoplasm, displayed severe fertility reduction, abnormal anther dehiscence, defective pollen development, decreased viability, and low seed-setting. It can be concluded that TaAGL66 plays an important role in wheat pollen development in the presence of Ae. kotschyi cytoplasm, providing new insights into the utilization of male sterility.

A SEPALLATA MADS-Box Transcription Factor, SlMBP21, Functions as a Negative Regulator of Flower Number and Fruit Yields in Tomato

MADS-box transcription factors act as the crucial regulators in plant organ differentiation. Crop yields are highly influenced by the flower number and fruit growth. However, flower identification is a very complex biological process, which involves many cascade regulations. The molecular mechanisms underlying the genetic regulation of flower identification in cultivated plants, such as tomato, are intricate and require further exploration. In this study, we investigated the vital function of a SEPALLATA (SEP) MADS-box gene, SlMBP21, in tomato sympodial inflorescence meristem (SIM) development for the conversion from SIMs to floral meristems (FMs). SlMBP21 transcripts were primarily accumulated in young inflorescence meristem, flowers, sepals, and abscission zones. The Ailsa Craig (AC++) tomato plants with suppressed SlMBP21 mRNA levels using RNAi exhibited a large increase in flower number and fruit yields in addition to enlarged sepals and inhibited abscission zone development. Scanning electron microscopy (SEM) revealed that the maturation of inflorescence meristems (IMs) was repressed in SlMBP21-RNAi lines. RNA-seq and qRT-PCR analyses showed that numerous genes related to the flower development, plant hormone signal transduction, cell cycle, and cell proliferation et al. were dramatically changed in SlMBP21-RNAi lines. Yeast two-hybrid assay exhibited that SlMBP21 can respectively interact with SlCMB1, SFT, JOINTLESS, and MC, which play key roles in inflorescence meristems or FM development. In summary, our data demonstrate that SlMBP21 functions as a key regulator in SIM development and the conversion from SIMs to FMs, through interacting with other regulatory proteins to control the expression of related genes.

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

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

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

KEY MESSAGE

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

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.

A co-located QTL for seven spike architecture-related traits shows promising breeding use potential in common wheat (Triticum aestivum L.)

KEY MESSAGE

A co-located novel QTL for TFS, FPs, FMs, FFS, FFPs, KWS, and KWPs with potential of improving wheat yield was identified and validated. Spike-related traits, including fertile florets per spike (FFS), kernel weight per spike (KWS), total florets per spike (TFS), florets per spikelet (FPs), florets in the middle spikelet (FMs), fertile florets per spikelet (FFPs), and kernel weight per spikelet (KWPs), are key traits in improving wheat yield. In the present study, quantitative trait loci (QTL) for these traits evaluated under various environments were detected in a recombinant inbred line population (msf/Chuannong 16) mainly genotyped using the 16 K SNP array. Ultimately, we identified 60 QTL, but only QFFS.sau-MC-1A for FFS was a major and stably expressed QTL. It was located on chromosome arm 1AS, where loci for TFS, FPs, FMs, FFS, FFPs, KWS, and KWPs were also simultaneously co-mapped. The effect of QFFS.sau-MC-1A was further validated in three independent segregating populations using a Kompetitive Allele-Specific PCR marker. For the co-located QTL, QFFS.sau-MC-1A, the presence of a positive allele from msf was associate with increases for all traits: + 12.29% TFS, + 10.15% FPs, + 13.97% FMs, + 17.12% FFS, + 14.75% FFPs, + 22.17% KWS, and + 19.42% KWPs. Furthermore, pleiotropy analysis showed that the positive allele at QFFS.sau-MC-1A simultaneously increased the spike length, spikelet number per spike, and thousand-kernel weight. QFFS.sau-MC-1A represents a novel QTL for marker-assisted selection with the potential for improving wheat yield. Four genes, TraesCS1A03G0012700, TraesCS1A03G0015700, TraesCS1A03G0016000, and TraesCS1A03G0016300, which may affect spike development, were predicted in the physical interval harboring QFFS.sau-MC-1A. Our results will help in further fine mapping QFFS.sau-MC-1A and be useful for improving wheat yield.

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

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

MADS1-regulated lemma and awn development benefits barley yield

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

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.

Grain yield improvement by genome editing of TaARF12 that decoupled peduncle and rachis development trajectories via differential regulation of gibberellin signalling in wheat

Plant breeding is constrained by trade-offs among different agronomic traits by the pleiotropic nature of many genes. Genes that contribute to two or more favourable traits with no penalty on yield are rarely reported, especially in wheat. Here, we describe the editing of a wheat auxin response factor TaARF12 by using CRISPR/Cas9 that rendered shorter plant height with larger spikes. Changes in plant architecture enhanced grain number per spike up to 14.7% with significantly higher thousand-grain weight and up to 11.1% of yield increase under field trials. Weighted Gene Co-Expression Network Analysis (WGCNA) of spatial-temporal transcriptome profiles revealed two hub genes: RhtL1, a DELLA domain-free Rht-1 paralog, which was up-regulated in peduncle, and TaNGR5, an organ size regulator that was up-regulated in rachis, in taarf12 plants. The up-regulation of RhtL1 in peduncle suggested the repression of GA signalling, whereas up-regulation of TaNGR5 in spike may promote GA response, a working model supported by differential expression patterns of GA biogenesis genes in the two tissues. Thus, TaARF12 complemented plant height reduction with larger spikes that gave higher grain yield. Manipulation of TaARF12 may represent a new strategy in trait pyramiding for yield improvement in wheat.

A major QTL simultaneously increases the number of spikelets per spike and thousand-kernel weight in a wheat line

KEY MESSAGE

A novel and stably expressed QTL QSNS.sicau-SSY-7A for spikelet number per spike in wheat without negative effects on thousand-kernel weight was identified and validated in different genetic backgrounds. Spikelet number per spike (SNS) is an important determinant of yield in wheat. In the present study, we combined bulked segregant analysis (BSA) and the wheat 660 K single-nucleotide polymorphism (SNP) array to rapidly identify genomic regions associated with SNS from a recombinant inbred line (RIL) population derived from a cross between the wheat lines S849-8 and SY95-71. A genetic map was constructed using Kompetitive Allele Specific PCR markers in the SNP-enriched region on the long arm of chromosome 7A. A major and stably expressed QTL, QSNS.sicau-SSY-7A, was detected in multiple environments. It was located in a 1.6 cM interval on chromosome arm 7AL flanked by the markers AX-109983514 and AX-109820548. This QTL explained 6.86-15.72% of the phenotypic variance, with LOD values ranging from 3.66 to 8.66. Several genes associated with plant growth and development were identified in the interval where QSNS.sicau-SSY-7A was located on the 'Chinese Spring' wheat and wild emmer reference genomes. Furthermore, the effects of QSNS.sicau-SSY-7A and WHEAT ORTHOLOG OFAPO1(WAPO1) on SNS were analyzed. Interestingly, QSNS.sicau-SSY-7A significantly increased SNS without negative effects on thousand-kernel weight, anthesis date and plant height, demonstrating its great potential for breeding aimed at improving grain yield. Taken together, these results indicate that QSNS.sicau-SSY-7A is a promising locus for yield improvement, and its linkage markers are helpful for fine mapping and molecular breeding.

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.

A high-resolution genotype-phenotype map identifies the TaSPL17 controlling grain number and size in wheat

BACKGROUND

Large-scale genotype-phenotype association studies of crop germplasm are important for identifying alleles associated with favorable traits. The limited number of single-nucleotide polymorphisms (SNPs) in most wheat genome-wide association studies (GWASs) restricts their power to detect marker-trait associations. Additionally, only a few genes regulating grain number per spikelet have been reported due to sensitivity of this trait to variable environments.

RESULTS

We perform a large-scale GWAS using approximately 40 million filtered SNPs for 27 spike morphology traits. We detect 132,086 significant marker-trait associations and the associated SNP markers are located within 590 associated peaks. We detect additional and stronger peaks by dividing spike morphology into sub-traits relative to GWAS results of spike morphology traits. We propose that the genetic dissection of spike morphology is a powerful strategy to detect signals for grain yield traits in wheat. The GWAS results reveal that TaSPL17 positively controls grain size and number by regulating spikelet and floret meristem development, which in turn leads to enhanced grain yield per plant. The haplotypes at TaSPL17 indicate geographical differentiation, domestication effects, and breeding selection.

CONCLUSION

Our study provides valuable resources for genetic improvement of spike morphology and a fast-forward genetic solution for candidate gene detection and cloning in wheat.

Genetic basis of geographical differentiation and breeding selection for wheat plant architecture traits

BACKGROUND

Plant architecture associated with increased grain yield and adaptation to the local environments is selected during wheat (Triticum aestivum) breeding. The internode length of individual stems and tiller length of individual plants are important for the determination of plant architecture. However, few studies have explored the genetic basis of these traits.

RESULTS

Here, we conduct a genome-wide association study (GWAS) to dissect the genetic basis of geographical differentiation of these traits in 306 worldwide wheat accessions including both landraces and traditional varieties. We determine the changes of haplotypes for the associated genomic regions in frequency in 831 wheat accessions that are either introduced from other countries or developed in China from last two decades. We identify 83 loci that are associated with one trait, while the remaining 247 loci are pleiotropic. We also find 163 associated loci are under strong selective sweep. GWAS results demonstrate independent regulation of internode length of individual stems and consistent regulation of tiller length of individual plants. This makes it possible to obtain ideal haplotype combinations of the length of four internodes. We also find that the geographical distribution of the haplotypes explains the observed differences in internode length among the worldwide wheat accessions.

CONCLUSION

This study provides insights into the genetic basis of plant architecture. It will facilitate gene functional analysis and molecular design of plant architecture for breeding.

The Genomic Variation and Differentially Expressed Genes on the 6P Chromosomes in Wheat-Agropyron cristatum Addition Lines 5113 and II-30-5 Confer Different Desirable Traits

Wild relatives of wheat are essential gene pools for broadening the genetic basis of wheat. Chromosome rearrangements and genomic variation in alien chromosomes are widespread. Knowledge of the genetic variation between alien homologous chromosomes is valuable for discovering and utilizing alien genes. In this study, we found that 5113 and II-30-5, two wheat-A. cristatum 6P addition lines, exhibited considerable differences in heading date, grain number per spike, and grain weight. Genome resequencing and transcriptome analysis revealed significant differences in the 6P chromosomes of the two addition lines, including 143,511 single-nucleotide polymorphisms, 62,103 insertion/deletion polymorphisms, and 757 differentially expressed genes. Intriguingly, genomic variations were mainly distributed in the middle of the chromosome arms and the proximal centromere region. GO and KEGG analyses of the variant genes and differentially expressed genes showed the enrichment of genes involved in the circadian rhythm, carbon metabolism, carbon fixation, and lipid metabolism, suggesting that the differential genes on the 6P chromosome are closely related to the phenotypic differences. For example, the photosynthesis-related genes PsbA, PsbT, and YCF48 were upregulated in II-30-5 compared with 5113. ACS and FabG are related to carbon fixation and fatty acid biosynthesis, respectively, and both carried modification variations and were upregulated in 5113 relative to II-30-5. Therefore, this study provides important guidance for cloning desirable genes from alien homologous chromosomes and for their effective utilization in wheat improvement.

Large-scale population structure and genetic architecture of agronomic traits of garlic

Garlic, an asexually propagated crop, is the second important bulb crop after the onion and is used as a vegetable and medicinal plant. Abundant and diverse garlic resources have been formed over thousands of years of cultivation. However, genome variation, population structure and genetic architecture of garlic agronomic traits were still not well elucidated. Here, 1 100 258 single nucleotide polymorphisms (SNPs) were identified using genotyping-by-sequencing in 606 garlic accessions collected from 43 countries. Population structure, principal component and phylogenetic analysis showed that these accessions were divided into five subpopulations. Twenty agronomic traits, including above-ground growth traits, bulb-related and bolt-related traits in two consecutive years were implemented in a genome-wide association study. In total, 542 SNPs were associated with these agronomic traits, among which 188 SNPs were repeatedly associated with more than two traits. One SNP (chr6: 1896135972) was repeatedly associated with ten traits. These associated SNPs were located within or near 858 genes, 56 of which were transcription factors. Interestingly, one non-synonymous SNP (Chr4: 166524085) in ribosomal protein S5 was repeatedly associated with above-ground growth and bulb-related traits. Additionally, gene ontology enrichment analysis of candidate genes for genomic selection regions between complete-bolting and non-bolting accessions showed that these genes were significantly enriched in 'vegetative to reproductive phase transition of meristem', 'shoot system development', 'reproductive process', etc. These results provide valuable information for the reliable and efficient selection of candidate genes to achieve garlic genetic improvement and superior varieties.

Global Transcriptomic Analyses Provide New Insight into the Molecular Mechanisms of Endocarp Formation and Development in Iron Walnut (Juglans sigillata Dode)

Iron walnut (Juglans sigillata Dode) is a native species in southwestern China that exhibits variation in both fruit morphology and shell thickness. However, the underlying molecular processes controlling hardened endocarp development in walnut has not yet been reported. Here, we generated transcriptional profiles of iron walnut endocarp at three developmental stages using “Dapao”, the most common commercial variety. Using pairwise comparisons between these three stages, a total of 8555 non-redundant differentially expressed genes (DEGs) were identified, and more than one-half of the total DEGs exhibited significant differential expression in stage I as compared with stage II or stage III, suggesting that the first stage may ultimately determine the final characteristics of the mature walnut shell. Furthermore, in the clustering analysis of the above DEGs, 3682, 2349, and 2388 genes exhibited the highest expression in stages I, II, and III, respectively. GO enrichment analysis demonstrated that the major transcriptional variation among the three developmental stages was caused by differences in cell growth, plant hormones, metabolic process, and phenylpropanoid metabolism. Namely, using the tissue-specific expression analysis and a gene co-expression network, we identified MADS-box transcription factor JsiFBP2 and bHLH transcription factor JsibHLH94 as candidate regulators of endocarp formation in the early stage, and JsiNAC56 and JsiMYB78 might play key roles in regulating the lignification process of endocarp in the late stage. This study provides useful information for further research to dissect the molecular mechanisms governing the shell formation and development of iron walnut.

The Torreya grandis genome illuminates the origin and evolution of gymnosperm-specific sciadonic acid biosynthesis

Torreya plants produce dry fruits with assorted functions. Here, we report the 19-Gb chromosome-level genome assembly of T. grandis. The genome is shaped by ancient whole-genome duplications and recurrent LTR retrotransposon bursts. Comparative genomic analyses reveal key genes involved in reproductive organ development, cell wall biosynthesis and seed storage. Two genes encoding a C₁₈ Δ⁹-elongase and a C₂₀ Δ⁵-desaturase are identified to be responsible for sciadonic acid biosynthesis and both are present in diverse plant lineages except angiosperms. We demonstrate that the histidine-rich boxes of the Δ⁵-desaturase are crucial for its catalytic activity. Methylome analysis reveals that methylation valleys of the T. grandis seed genome harbor genes associated with important seed activities, including cell wall and lipid biosynthesis. Moreover, seed development is accompanied by DNA methylation changes that possibly fuel energy production. This study provides important genomic resources and elucidates the evolutionary mechanism of sciadonic acid biosynthesis in land plants.

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.

Genome-Wide Identification of MADS-Box Family Genes in Safflower (Carthamus tinctorius L.) and Functional Analysis of CtMADS24 during Flowering

Safflower is an important economic crop with a plethora of industrial and medicinal applications around the world. The bioactive components of safflower petals are known to have pharmacological activity that promotes blood circulation and reduces blood stasis. However, fine-tuning the genetic mechanism of flower development in safflower is still required. In this study, we report the genome-wide identification of MADS-box transcription factors in safflower and the functional characterization of a putative CtMADS24 during vegetative and reproductive growth. In total, 77 members of MADS-box-encoding genes were identified from the safflower genome. The phylogenetic analysis divided CtMADS genes into two types and 15 subfamilies. Similarly, bioinformatic analysis, such as of conserved protein motifs, gene structures, and cis-regulatory elements, also revealed structural conservation of MADS-box genes in safflower. Furthermore, the differential expression pattern of CtMADS genes by RNA-seq data indicated that type II genes might play important regulatory roles in floral development. Similarly, the qRT-PCR analysis also revealed the transcript abundance of 12 CtMADS genes exhibiting tissue-specific expression in different flower organs. The nucleus-localized CtMADS24 of the AP1 subfamily was validated by transient transformation in tobacco using GFP translational fusion. Moreover, CtMADS24-overexpressed transgenic Arabidopsis exhibited early flowering and an abnormal phenotype, suggesting that CtMADS24 mediated the expression of genes involved in floral organ development. Taken together, these findings provide valuable information on the regulatory role of CtMADS24 during flower development in safflower and for the selection of important genes for future molecular breeding programs.

Genome-Wide Identification of MIKCc-Type MADS-Box Family Gene and Floral Organ Transcriptome Characterization in Ma Bamboo (Dendrocalamus latiflorus Munro)

Most bamboos die after flowering, and the molecular mechanisms responsible for flowering is poorly understood. The MIKCc-type MADS-box family gene is involved in the flowering process. To explore the mechanism of the MIKCc-type MADS-box gene and phytohormone regulation in the flowering of Dendrocalamus latiflorus Munro (D. latiflorus), characterized by extremely rapid growth and widely cultivated woody bamboo, we initially did a genome-wide analysis of the MIKCc-type MADS-box gene in D. latiflorus. In the meantime, transcriptome analysis was performed using the floral organs. A total of 170 MIKCc-Type MADS-Box genes were identified and divided into 15 categories. The cis-acting element analysis in promoters regions revealed that MIKC-type MADS-box family genes were associated with hormones, including auxin, abscisic acid (ABA), gibberellin (GA) and jasmonic acid (JA), which was found at 79, 476, 96, 486 sites and cover 61, 103, 73, 128 genes. Genome synteny analysis showed subgenome AA and BB were better than CC and obtained 49, 40, 39 synteny genes compared with Oryza sativa (O. sativa). In transcriptome analysis of floral organs, the enriched pathway from DEGs included circadian, vernalization and gibberellin pathways associated with the flowering process. We found that the jasmonic acid synthesis gene is highly expressed in the pistil, which may be the cause of Ma bamboo pollen abortion. The expression profile showed that most MIKC-type MADS-box genes exhibited high expression in flower organs. The consequences of this study will provide insight into the irregular flowering and low pollen counts of Ma bamboo.

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.

Genome-wide analysis of MADS-box gene family in kiwifruit (Actinidia chinensis var. chinensis) and their potential role in floral sex differentiation

Kiwifruit (Actinidia chinensis Planch.) is a functionally dioecious plant, which displays diverse morphology in male and female flowers. MADS-box is an ancient and huge gene family that plays a key role in plant floral organ differentiation. In this study, we have identified 89 MADS-box genes from A. chinensis Red 5 genome. These genes are distributed on 26 chromosomes and are classified into type I (21 genes) and type II (68 genes). Overall, type II AcMADS-box genes have more complex structures than type I with more exons, protein domains, and motifs, indicating that type II genes may have more diverse functions. Gene duplication analysis showed that most collinearity occurred in type II AcMADS-box genes, which was consistent with a large number of type II genes. Analysis of cis-acting elements in promoters showed that AcMADS-box genes are mainly associated with light and phytohormone responsiveness. The expression profile of AcMADS-box genes in different tissues showed that most genes were highly expressed in flowers. Further, the qRT-PCR analysis of the floral organ ABCDE model-related genes in male and female flowers revealed that AcMADS4, AcMADS56, and AcMADS70 were significantly expressed in female flowers. It indicated that those genes may play an important role in the sex differentiation of kiwifruit. This work provided a comprehensive analysis of the AcMADS-box genes and may help facilitate our understanding of the sex differentiation regulatory mechanism in kiwifruit.

Transcriptional signatures of wheat inflorescence development

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

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.

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

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

Proteomic and Phosphoproteomic Analyses Reveal a Complex Network Regulating Pollen Abortion and Potential Candidate Proteins in TCMS Wheat

Thermosensitive sterile lines are natural materials for exploring the effects of anther development on male fertility. To study the possible molecular mechanisms regulating protein activity during the induction of male sterility, proteomic and phosphoproteomic analyses with tandem mass tags (TMTs) were used to study the binucleate anther of the thermosensitive sterile wheat line YS3038. A total of 9072 proteins, including 5019 phosphoproteins, were identified. Enrichment analyses of differentially abundant proteins (DAPs) and phosphoproteins (DAPPs) in metabolic pathways showed that both were mainly related to energy metabolism. Soluble sugar and ATP content were significantly decreased, free fatty acid content was significantly increased, and ROS was abnormally accumulated in male sterile YS3038-A. In addition, 233 kinase–substrate pairs involved in potential phosphorylation control networks were predicted to regulate fertility. Candidate proteins were identified, and a quantitative real-time polymerase chain reaction (qRT-PCR) analysis was used to validate the TMT results. TaPDCD5 is likely to be involved in fertility conversion of YS3038 by barley stripe mosaic virus-induced gene silencing (BSMV-VIGS). Our data provide new insights into the mechanism of TCMS, which has value for identifying potential candidate proteins associated with the formation or abortion of pollen and promotion of wheat heterosis utilization.

InDels Identification and Association Analysis with Spike and Awn Length in Chinese Wheat Mini-Core Collection

Diversity surveys of germplasm are important for gaining insight into the genomic basis for crop improvement; especially InDels, which are poorly understood in hexaploid common wheat. Here, we describe a map of 89,923 InDels from exome sequencing of 262 accessions of a Chinese wheat mini-core collection. Population structure analysis, principal component analysis and selective sweep analysis between landraces and cultivars were performed. Further genome-wide association study (GWAS) identified five QTL (Quantitative Trait Loci) that were associated with spike length, two of them, on chromosomes 2B and 6A, were detected in 10 phenotypic data sets. Assisted with RNA-seq data, we identified 14 and 21 genes, respectively that expressed in spike and rachis within the two QTL regions that can be further investigated for candidate genes discovery. Moreover, InDels were found to be associated with awn length on chromosomes 5A, 6B and 4A, which overlapped with previously reported genetic loci B1 (Tipped 1), B2 (Tipped 2) and Hd (Hooded). One of the genes TaAGL6 that was previously shown to affect floral organ development was found at the B2 locus to affect awn length development. Our study shows that trait-associated InDels may contribute to wheat improvement and may be valuable molecular markers for future wheat breeding.

Comparative Analysis of the MADS-Box Genes Revealed Their Potential Functions for Flower and Fruit Development in Longan (Dimocarpus longan)

Longan (Dimocarpus longan Lour.) is an important economic crop widely planted in tropical and subtropical regions, and flower and fruit development play decisive effects on the longan yield and fruit quality formation. MCM1, AGAMOUS, DEFICIENS, Serum Response Factor (MADS)-box transcription factor family plays important roles for the flowering time, floral organ identity, and fruit development in plants. However, there is no systematic information of MADS-box family in longan. In this study, 114 MADS-box genes were identified from the longan genome, phylogenetic analysis divided them into type I (Mα, Mβ, Mγ) and type II (MIKC*, MIKC ^C ) groups, and MIKC ^C genes were further clustered into 12 subfamilies. Comparative genomic analysis of 12 representative plant species revealed the conservation of type II in Sapindaceae and analysis of cis-elements revealed that Dof transcription factors might directly regulate the MIKC ^C genes. An ABCDE model was proposed for longan based on the phylogenetic analysis and expression patterns of MADS-box genes. Transcriptome analysis revealed that MIKC ^C genes showed wide expression spectrums, particularly in reproductive organs. From 35 days after KClO₃ treatment, 11 MIKC genes were up-regulated, suggesting a crucial role in off-season flower induction, while DlFLC, DlSOC1, DlSVP, and DlSVP-LIKE may act as the inhibitors. The gene expression patterns of longan fruit development indicated that DlSTK, DlSEP1/2, and DlMADS53 could be involved in fruit growth and ripening. This paper carried out the whole genome identification and analysis of the longan MADS-box family for the first time, which provides new insights for further understanding its function in flowers and fruit.