SiMADS34, an E-class MADS-box transcription factor, regulates inflorescence architecture and grain yield in Setaria italica

Efficient identification of QTL for agronomic traits in foxtail millet (Setaria italica) using RTM- and MLM-GWAS

KEY MESSAGE

Eleven QTLs for agronomic traits were identified by RTM- and MLM-GWAS, putative candidate genes were predicted and two markers for grain weight were developed and validated. Foxtail millet (Setaria italica), the second most cultivated millet crop after pearl millet, is an important grain crop in arid regions. Seven agronomic traits of 408 diverse foxtail millet accessions from 15 provinces in China were evaluated in three environments. They were clustered into two divergent groups based on genotypic data using ADMIXTURE, which was highly consistent with their geographical distribution. Two models for genome-wide association studies (GWAS), namely restricted two-stage multi-locus multi-allele (RTM)-GWAS and mixed linear model (MLM)-GWAS, were used to dissect the genetic architecture of the agronomic traits based on 13,723 SNPs. Eleven quantitative trait loci (QTLs) for seven traits were identified using two models (RTM- and MLM-GWAS). Among them, five were considered stable QTLs that were identified in at least two environments using MLM-GWAS. One putative candidate gene (SETIT_006045mg, Chr4: 744,701-746,852) that can enhance grain weight per panicle was identified based on homologous gene comparison and gene expression analysis and was validated by haplotype analysis of 330 accessions with high-depth (10×) resequencing data (unpublished). In addition, homologous gene comparison and haplotype analysis identified one putative foxtail millet ortholog (SETIT_032906mg, Chr2: 5,020,600-5,029,771) with rice affecting the target traits. Two markers (cGWP6045 and kTGW2906) were developed and validated and can be used for marker-assisted selection of foxtail millet with high grain weight. The results provide a fundamental resource for foxtail millet genetic research and breeding and demonstrate the power of integrating RTM- and MLM-GWAS approaches as a complementary strategy for investigating complex traits in foxtail millet.

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.

DbMADS regulates carotenoid metabolism by repressing two carotenogenic genes in the green alga Dunaliella sp. FACHB-847

MADS transcription factors are involved in the regulation of fruit development and carotenoid metabolism in plants. However, whether and how carotenoid accumulation is regulated by algal MADS are largely unknown. In this study, we first used functional complementation to confirm the functional activity of phytoene synthase from the lutein-rich Dunaliella sp. FACHB-847 (DbPSY), the key rate-limiting enzyme in the carotenoid biosynthesis. Promoters of DbPSY and DbLcyB (lycopene β-cyclase) possessed multiple cis-acting elements such as light-, UV-B-, dehydration-, anaerobic-, and salt-responsive elements, W-box, and C-A-rich-G-box (MADS-box). Meanwhile, we isolated one nucleus-localized MADS transcription factor (DbMADS), belonging to type I MADS gene. Three carotenogenic genes, DbPSY, DbLcyB, and DbBCH (β-carotene hydroxylase) genes were upregulated at later stages, which was well correlated with the carotenoid accumulation. In contrast, DbMADS gene was highly expressed at lag phase with low carotenoid accumulation. Yeast one-hybrid assay and dual-luciferase reporter assay demonstrated that DbMADS could directly bind to the promoters of two carotenogenic genes, DbPSY and DbLcyB, and repress their transcriptions. This study suggested that DbMADS may act as a negative regulator of carotenoid biosynthesis by repressing DbPSY and DbLcyB at the lag phase, which provide new insights into the regulatory mechanisms of carotenoid metabolism in Dunaliella.

Molecular evolution, diversification, and expression assessment of MADS gene family in Setaria italica, Setaria viridis, and Panicum virgatum

KEY MESSAGE

This paper sheds light on the evolution and expression patterns of MADS genes in Setaria and Panicum virgatum. SiMADS51 and SiMADS64 maybe involved in the ABA-dependent pathway of drought response. The MADS gene family is a key regulatory factor family that controls growth, reproduction, and response to abiotic stress in plants. However, the molecular evolution of this family is rarely reported. Here, a total of 265 MADS genes were identified in Setaria italica (foxtail millet), Setaria viridis (green millet), and Panicum virgatum (switchgrass) and analyzed by bioinformatics, including physicochemical characteristics, subcellular localization, chromosomal position and duplicate, motif distribution, genetic structure, genetic evolvement, and expression patterns. Phylogenetic analysis was used to categorize these genes into M and MIKC types. The distribution of motifs and gene structure were similar for the corresponding types. According to a collinearity study, the MADS genes have been mostly conserved during evolution. The principal cause of their expansion is segmental duplication. However, the MADS gene family tends to shrink in foxtail millet, green millet, and switchgrass. The MADS genes were subjected to purifying selection, but several positive selection sites were also identified in three species. And most of the promoters of MADS genes contain cis-elements related to stress and hormonal response. RNA-seq and quantitative Real-time PCR (qRT-PCR) analysis also were examined. SiMADS genes expression levels are considerably changed in reaction to various treatments, following qRT-PCR analysis. This sheds fresh light on the evolution and expansion of the MADS family in foxtail millet, green millet, and switchgrass, and lays the foundation for further research on their functions.

Genome-wide identification of the MIKCc-type MADS-box gene family in Chrysanthemum lavandulifolium reveals their roles in the capitulum development

Chrysanthemum ×morifolium is well known throughout the world for its diverse and exquisite flower types. However, due to the complicated genetic background of C. ×morifolium, it is difficult to understand the molecular mechanism of its flower development. And it limits the molecular breeding of improving chrysanthemum flower types. C. ×morifolium has the typical radial capitulum, and many researches showed that the members of the MIKCc-type MADS box gene family play a key role in the formation and development of the capitulum. However, it has been difficult to isolate the important MIKCc and investigate their roles in this process due to the lack of genomic information in chrysanthemum. Here, we identified MIKCc-type MADS box genes at whole genome-wide level in C. lavandulifolium, a diploid species closely related to C. ×morifolium, and investigated their roles in capitulum development by gene expression pattern analysis and protein interaction analysis. A total of 40 ClMIKCc were identified and were phylogenetically grouped into 12 clades. Members of all clades showed different enriched expression patterns during capitulum formation. We speculate that the E-class genes in C. lavandulifolium underwent subfunctionalization because they have a significantly expanded, more diverse expression patterns, and specifically tissue expression than AtSEPs. Meanwhile, we detected the C-class expressed in disc floret corolla, which could be the clue to explore the morphological differences between disc and ray floret corolla. In addition, the potential roles of some MIKCcs in complex inflorescence formation were explored by comparing the number and phylogenetic relationship of MIKCc subfamily members in Asteraceae with different capitulum types. Members of the FLC branch in Asteraceae were found to be possibly related to the differentiation and development of the ray floret.

Identification and haplotype analysis of SiCHLI: a gene for yellow-green seedling as morphological marker to accelerate foxtail millet (Setaria italica) hybrid breeding

We cloned and developed functional markers for the SiCHLI gene, which is responsible for the yellow-green color of leaves in foxtail millet, a frequently used marker trait in the hybrid breeding of foxtail millet by using bulked segregant analysis sequencing and haplotype analysis on the F₂ and core-collected nature populations. The color of leaves has been widely used as a marker for the hybrid breeding of foxtail millet; however, few related gene have been cloned to date. Here, we used two F₂ populations generated from crosses between the highly male-sterile material 125A with yellow-green leaves, and CG58 and S410, which have green leaves, to identify the genes underlying the yellow-green color of the leaves of foxtail millet. The leaves of 125A seedlings were yellow-green, but they became green at the heading stage. The content of chlorophyll a and chlorophyll b was lower, the number of thylakoid lamellae and grana was reduced, and the chloroplasts was more rounded in 125A than in S410 at the yellow-green leaf stage; however, no differences were observed between 125A and S410 in these traits and photosynthetic at the heading stage. Bulked segregant analysis and map-based cloning revealed that the SiCHLI gene is responsible for the leaf colors of 125A. A nonsynonymous mutation (C/T) in exon 3 causes yellow-green leaves in 125A at the seedling stage. Haplotype analysis of the SiCHLI gene in 596 core collected accessions revealed a new haplotype associated with high photosynthetic metabolic potential at the heading and mature stages, which could be used to enhance sterile lines with yellow-green leaves. We developed a functional marker that will facilitate the identification of foxtail millet accessions with the different types of yellow-green leaves. Generally, our study provides new genetic resources to guide the future marker-assisted or target-base editing in foxtail millet hybrid breeding.

The Integration of Genome-Wide Association Study and Homology Analysis to Explore the Genomic Regions and Candidate Genes for Panicle-Related Traits in Foxtail Millet

Panicle traits are important factors affecting yield, and their improvement has long been a critical goal in foxtail millet breeding. In order to understand the genetic basis of panicle formation, a large-scale genome-wide association study (GWAS) was performed in this study for six panicle-related traits based on 706,646 high-polymorphism SNP loci in 407 accessions. As a result, 87 quantitative trait loci (QTL) regions with a physical distance of less than 100 kb were detected to be associated with these traits in three environments. Among them, 27 core regions were stably detected in at least two environments. Based on rice-foxtail millet homologous comparison, expression, and haplotype analysis, 27 high-confidence candidate genes in the QTL regions, such as Si3g11200 (OsDER1), Si1g27910 (OsMADS6), Si7g27560 (GS5), etc., affected panicle-related traits by involving multiple plant growth regulator pathways, a photoperiod response, as well as panicle and grain development. Most of these genes showed multiple effects on different panicle-related traits, such as Si3g11200 affecting all six traits. In summary, this study clarified a strategy based on the integration of GWAS, a homologous comparison, and haplotype analysis to discover the genomic regions and candidate genes for important traits in foxtail millet. The detected QTL regions and candidate genes could be further used for gene clone and marker-assisted selection in foxtail millet breeding.

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.

Conuping BSA-Seq and RNA-Seq Reveal the Molecular Pathway and Genes Associated with the Plant Height of Foxtail Millet (Setaria italica)

Foxtail millet (Setaria italica) plays an important role in C4 crop research and agricultural development in arid areas due to its short growth period, drought tolerance, and barren tolerance. Exploration of the dwarfing mechanism and the dwarf genes of foxtail millet can provide a reference for dwarf breeding and dwarf research of other C4 crops. In this study, genetic analysis was performed using phenotypic data, candidate genes were screened by bulk segregant analysis sequencing (BSA-Seq); differentially expressed genes and metabolic pathways in different strains of high samples were analyzed by RNA sequencing (RNA-Seq). The association analysis of BSA-Seq and RNA-Seq further narrowed the candidate range. As a result, a total of three quantitative trait loci (QTLs) and nine candidate genes related to plant height were obtained on chromosomes I and IX. Based on the functional prediction of the candidate genes, we propose a hypothetical mechanism for the formation of millet dwarfing, in which, metabolism and MAPK signaling play important roles in the formation of foxtail millet plant height.

Identification and expression analysis of the MADS-box genes of Kentucky bluegrass during inflorescence development

MADS-box genes play vital roles in multiple biological processes of plants growth and development, especially inflorescence development. In the present study, a comprehensive investigation into the identification and classification of MADS-box genes in Kentucky bluegrass (Poa pratensis) has not been reported. Here, based on the transcriptome of inflorescence, we identified 44 PpMADS-box genes, and gave an overview of the physicochemical properties, phylogeny, protein structures, and potential functions of the proteins encoded by these genes through various bioinformatics software for the first time. Analysis of physicochemical properties revealed that most PpMADS-box were alkaline proteins and possessed similar conserved motifs. Additionally, it was demonstrated that 33 PpMADS-box proteins without signal peptide, leading peptide, transmembrane structure and located in the nucleus were not transported or secreted, so directly played transcriptional regulatory roles in the nucleus. Then, peptide sequences BLAST search and analysis of phylogenetic relationships with MADS-box proteins of P. pratensis, Arabidopsis thaliana, and Oryza sativa were performed. It was found that 44 PpMADS-box proteins were separated into 33 MIKC-type (3 BS, 1 AGL17, 8 AP3/P2, 3 AP1, 5 SEP, 6 SOC and 7 AG genes, respectvely) and 11 type I-type, which include 7 Mγ and 4 Mα. Furthermore, the relative expression levels of the selected 12 genes (MADS3, 15, 16, 17, 18, 20, 24, 27, 30, 36, 38 and 40) at the booting stage, pre-anthesis, anthesis, post-anthesis, and seed filling stage of inflorescences, as well as leaves and roots of the corresponding stages of inflorescences were analyzed, showing that most PpMADS-box genes were highly expressed mainly in young leaves and later inflorescences, and had complex patters in roots. Morever, except for PpMADS30 being highly expressed in the leaves, others were significantly highly expressed in inflorescence and/ or roots, demonstrating PpMADS-box genes also regulate leaves and roots development in plant. This study provides valuable insights into the MADS-box family genes in Kentucky bluegrass and its potential functional characteristics, expression pattern, and evolution in floral organogenesis and even reproduction development. @media print { .ms-editor-squiggler { display:none !important; } } .ms-editor-squiggler { all: initial; display: block !important; height: 0px !important; width: 0px !important; }.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s12298-022-01216-1.

Machine-Learning-Based Genome-Wide Association Studies for Uncovering QTL Underlying Soybean Yield and Its Components

A genome-wide association study (GWAS) is currently one of the most recommended approaches for discovering marker-trait associations (MTAs) for complex traits in plant species. Insufficient statistical power is a limiting factor, especially in narrow genetic basis species, that conventional GWAS methods are suffering from. Using sophisticated mathematical methods such as machine learning (ML) algorithms may address this issue and advance the implication of this valuable genetic method in applied plant-breeding programs. In this study, we evaluated the potential use of two ML algorithms, support-vector machine (SVR) and random forest (RF), in a GWAS and compared them with two conventional methods of mixed linear models (MLM) and fixed and random model circulating probability unification (FarmCPU), for identifying MTAs for soybean-yield components. In this study, important soybean-yield component traits, including the number of reproductive nodes (RNP), non-reproductive nodes (NRNP), total nodes (NP), and total pods (PP) per plant along with yield and maturity, were assessed using a panel of 227 soybean genotypes evaluated at two locations over two years (four environments). Using the SVR-mediated GWAS method, we were able to discover MTAs colocalized with previously reported quantitative trait loci (QTL) with potential causal effects on the target traits, supported by the functional annotation of candidate gene analyses. This study demonstrated the potential benefit of using sophisticated mathematical approaches, such as SVR, in a GWAS to complement conventional GWAS methods for identifying MTAs that can improve the efficiency of genomic-based soybean-breeding programs.

Genome-wide identification of quantitative trait loci for morpho-agronomic and yield-related traits in foxtail millet (Setaria italica) across multi-environments

Foxtail millet (Setaria italica) is an ideal model of genetic system for functional genomics of the Panicoideae crop. Identification of QTL responsible for morpho-agronomic and yield-related traits facilitates dissection of genetic control and breeding in cereal crops. Here, based on a Yugu1 × Longgu7 RIL population and genome-wide resequencing data, an updated linkage map harboring 2297 bin and 74 SSR markers was constructed, spanning 1315.1 cM with an average distance of 0.56 cM between adjacent markers. A total of 221 QTL for 17 morpho-agronomic and yield-related traits explaining 5.5 ~ 36% of phenotypic variation were identified across multi-environments. Of these, 109 QTL were detected in two to nine environments, including the most stable qLMS6.1 harboring a promising candidate gene Seita.6G250500, of which 70 were repeatedly identified in different trials in the same geographic location, suggesting that foxtail millet has more identical genetic modules under the similar ecological environment. One hundred-thirty QTL with overlapping intervals formed 22 QTL clusters. Furthermore, six superior recombinant inbred lines, RIL35, RIL48, RIL77, RIL80, RIL115 and RIL125 with transgressive inheritance and enrichment of favorable alleles in plant height, tiller, panicle morphology and yield related-traits were screened by hierarchical cluster. These identified QTL, QTL clusters and superior lines lay ground for further gene-trait association studies and breeding practice in foxtail millet.

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.

Genetic control and phenotypic characterization of panicle architecture and grain yield-related traits in foxtail millet (Setaria italica)

Multi-environment QTL mapping identified 23 stable loci and 34 co-located QTL clusters for panicle architecture and grain yield-related traits, which provide a genetic basis for foxtail millet yield improvement. Panicle architecture and grain weight, both of which are influenced by genetic and environmental factors, have significant effects on grain yield potential. Here, we used a recombinant inbred line (RIL) population of 333 lines of foxtail millet, which were grown in 13 trials with varying environmental conditions, to identify quantitative trait loci (QTL) controlling nine agronomic traits related to panicle architecture and grain yield. We found that panicle weight, grain weight per panicle, panicle length, panicle diameter, and panicle exsertion length varied across different geographical locations. QTL mapping revealed 159 QTL for nine traits. Of the 159 QTL, 34 were identified in 2 to 12 environments, suggesting that the genetic control of panicle architecture in foxtail millet is sensitive to photoperiod and/or other environmental factors. Eighty-eight QTL controlling different traits formed 34 co-located QTL clusters, including the triple QTL cluster qPD9.2/qPL9.5/qPEL9.3, which was detected 23 times in 13 environments. Several candidate genes, including Seita.2G388700, Seita.3G136000, Seita.4G185300, Seita.5G241500, Seita.5G243100, Seita.9G281300, and Seita.9G342700, were identified in the genomic intervals of multi-environmental QTL or co-located QTL clusters. Using available phenotypic and genotype data, we conducted haplotype analysis for Seita.2G002300 and Seita.9G064000,which showed high correlations with panicle weight and panicle exsertion length, respectively. These results not only provided a basis for further fine mapping, functional studies and marker-assisted selection of traits related to panicle architecture in foxtail millet, but also provide information for comparative genomics analyses of cereal crops.

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

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

Genome-Wide DNA Polymorphism Analysis and Molecular Marker Development for the Setaria italica Variety "SSR41" and Positional Cloning of the Setaria White Leaf Sheath Gene SiWLS1

Genome-wide DNA polymorphism analysis and molecular marker development are important for forward genetics research and DNA marker-assisted breeding. As an ideal model system for Panicoideae grasses and an important minor crop in East Asia, foxtail millet (Setaria italica) has a high-quality reference genome as well as large mutant libraries based on the "Yugu1" variety. However, there is still a lack of genetic and mutation mapping tools available for forward genetics research on S. italica. Here, we screened another S. italica genotype, "SSR41", which is morphologically similar to, and readily cross-pollinates with, "Yugu1". High-throughput resequencing of "SSR41" identified 1,102,064 reliable single nucleotide polymorphisms (SNPs) and 196,782 insertions/deletions (InDels) between the two genotypes, indicating that these two genotypes have high genetic diversity. Of the 8,361 high-quality InDels longer than 20 bp that were developed as molecular markers, 180 were validated with 91.5% accuracy. We used "SSR41" and these developed molecular markers to map the white leaf sheath gene SiWLS1. Further analyses showed that SiWLS1 encodes a chloroplast-localized protein that is involved in the regulation of chloroplast development in bundle sheath cells in the leaf sheath in S. italica and is related to sensitivity to heavy metals. Our study provides the methodology and an important resource for forward genetics research on Setaria.