Genome-wide identification, phylogenetic and expression pattern analysis of MADS-box family genes in foxtail millet (Setaria italica)

Tribe Paniceae Cereals with Different Ploidy Levels: Setaria italica, Panicum miliaceum, and Echinochloa esculenta

Plants have repeatedly undergone whole-genome duplication during their evolutionary history. Even in modern plants, there is diversity in ploidy within and between species, providing a snapshot of the evolutionary turnover of ploidy. Here, I will review the diversity of ploidy and the evolution of the genome constitution, focusing on the millet species Setaria italica, Panicum miliaceum, and Echinochloa esculenta. These are all historically important cereal crops that have been domesticated in East Asia. They all display a basic chromosome set of nine, but they are diploid, tetraploid, and hexaploid, respectively. The timing of ploidy is different among the millet species, as is the extent of gene family expansion and gene loss. There also exists complex subgenomic evolution in the wild species within each genus. These three millet species and their related wild species are suitable models for elucidating the molecular evolution and diversity of genome duplication by comparative genomic analysis.

Realizing visionary goals for the International Year of Millet (IYoM): accelerating interventions through advances in molecular breeding and multiomics resources

MAIN CONCLUSION

Leveraging advanced breeding and multi-omics resources is vital to position millet as an essential "nutricereal resource," aligning with IYoM goals, alleviating strain on global cereal production, boosting resilience to climate change, and advancing sustainable crop improvement and biodiversity. The global challenges of food security, nutrition, climate change, and agrarian sustainability demand the adoption of climate-resilient, nutrient-rich crops to support a growing population amidst shifting environmental conditions. Millets, also referred to as "Shree Anna," emerge as a promising solution to address these issues by bolstering food production, improving nutrient security, and fostering biodiversity conservation. Their resilience to harsh environments, nutritional density, cultural significance, and potential to enhance dietary quality index made them valuable assets in global agriculture. Recognizing their pivotal role, the United Nations designated 2023 as the "International Year of Millets (IYoM 2023)," emphasizing their contribution to climate-resilient agriculture and nutritional enhancement. Scientific progress has invigorated efforts to enhance millet production through genetic and genomic interventions, yielding a wealth of advanced molecular breeding technologies and multi-omics resources. These advancements offer opportunities to tackle prevailing challenges in millet, such as anti-nutritional factors, sensory acceptability issues, toxin contamination, and ancillary crop improvements. This review provides a comprehensive overview of molecular breeding and multi-omics resources for nine major millet species, focusing on their potential impact within the framework of IYoM. These resources include whole and pan-genome, elucidating adaptive responses to abiotic stressors, organelle-based studies revealing evolutionary resilience, markers linked to desirable traits for efficient breeding, QTL analysis facilitating trait selection, functional gene discovery for biotechnological interventions, regulatory ncRNAs for trait modulation, web-based platforms for stakeholder communication, tissue culture techniques for genetic modification, and integrated omics approaches enabled by precise application of CRISPR/Cas9 technology. Aligning these resources with the seven thematic areas outlined by IYoM catalyzes transformative changes in millet production and utilization, thereby contributing to global food security, sustainable agriculture, and enhanced nutritional consequences.

Omics approaches to understand the MADS-box gene family in common bean (Phaseolus vulgaris L.) against drought stress

MADS-box genes are known to play important roles in diverse aspects of growth/devolopment and stress response in several plant species. However, no study has yet examined about MADS-box genes in P. vulgaris. In this study, a total of 79 PvMADS genes were identified and classified as type I and type II according to the phylogenetic analysis. While both type I and type II PvMADS classes were found to contain the MADS domain, the K domain was found to be present only in type II PvMADS proteins, in agreement with the literature. All chromosomes of the common bean were discovered to contain PvMADS genes and 17 paralogous gene pairs were identified. Only two of them were tandemly duplicated gene pairs (PvMADS-19/PvMADS-23 and PvMADS-20/PvMADS-24), and the remaining 15 paralogous gene pairs were segmentally duplicated genes. These duplications were found to play an important role in the expansion of type II PvMADS genes. Moreover, the RNAseq and RT-qPCR analyses showed the importance of PvMADS genes in response to drought stress in P. vulgaris.

Evolutionary analysis of MADS-box genes in buckwheat species and functional study of FdMADS28 in flavonoid metabolism

The MADS-box gene family is a transcription factor family that is widely expressed in plants. It controls secondary metabolic processes in plants and encourages the development of tissues like roots and flowers. However, the phylogenetic analysis and evolutionary model of MADS-box genes in Fagopyrum species has not been reported yet. This study identified the MADS-box genes of three buckwheat species at the whole genome level, and conducted systematic evolution and physicochemical analysis. The results showed that these genes can be divided into four subfamilies, with fragment duplication being the main way for the gene family expansion. During the domestication process from golden buckwheat to tartary buckwheat and the common buckwheat, the Ka/Ks ratio indicated that most members of the family experienced strong purification selection pressure, and with individual gene pairs experiencing positive selection. In addition, we combined the expression profile data of the MADS genes, mGWAS data, and WGCNA data to mine genes FdMADS28/48/50 that may be related to flavonoid metabolism. The results also showed that overexpression of FdMADS28 could increase rutin content by decreasing Kaempferol pathway content in hairy roots, and increase the resistance and growth of hairy roots to PEG and NaCl. This study systematically analyzed the evolutionary relationship of MADS-box genes in the buckwheat species, and elaborated on the expression patterns of MADS genes in different tissues under biotic and abiotic stresses, laying an important theoretical foundation for further elucidating their role in flavonoid metabolism.

Genome-Wide Characterization and Phylogenetic and Stress Response Expression Analysis of the MADS-Box Gene Family in Litchi (Litchi chinensis Sonn.)

The MADS-box protein is an important transcription factor in plants and plays an important role in regulating the plant abiotic stress response. In this study, a total of 94 MADS-box genes were predicted in the litchi genome, and these genes were widely distributed on all the chromosomes. The LcMADS-box gene family was divided into six subgroups (Mα, Mβ, Mγ, Mδ, MIKC, and UN) based on their phylogenetical relationships with Arabidopsis, and the closely linked subgroups exhibited more similarity in terms of motif distribution and intron/exon numbers. Transcriptome analysis indicated that LcMADS-box gene expression varied in different tissues, which can be divided into universal expression and specific expression. Furthermore, we further validated that LcMADS-box genes can exhibit different responses to various stresses using quantitative real-time PCR (qRT-PCR). Moreover, physicochemical properties, subcellular localization, collinearity, and cis-acting elements were also analyzed. The findings of this study provide valuable insights into the MADS-box gene family in litchi, specifically in relation to stress response. The identification of hormone-related and stress-responsive cis-acting elements in the MADS-box gene promoters suggests their involvement in stress signaling pathways. This study contributes to the understanding of stress tolerance mechanisms in litchi and highlights potential regulatory mechanisms underlying stress responses.

Genome-wide identification of the MADS-box gene family in Avena sativa and its role in photoperiod-insensitive oat

Background

Traditional spring-summer sown oat is a typical long-day crop that cannot head under short-day conditions. The creation of photoperiod-insensitive oats overcomes this limitation. MADS-box genes are a class of transcription factors involved in plant flowering signal transduction regulation. Previous transcriptome studies have shown that MADS-box genes may be related to the oat photoperiod.

Methods

Putative MADS-box genes were identified in the whole genome of oat. Bioinformatics methods were used to analyze their classification, conserved motifs, gene structure, evolution, chromosome localization, collinearity and cis-elements. Ten representative genes were further screened via qRT‒PCR analysis under short days.

Results

In total, sixteen AsMADS genes were identified and grouped into nine subfamilies. The domains, conserved motifs and gene structures of all AsMADS genes were conserved. All members contained light-responsive elements. Using the photoperiod-insensitive oat MENGSIYAN4HAO (MSY4) and spring-summer sown oat HongQi2hao (HQ2) as materials, qRT‒PCR analysis was used to analyze the AsMADS gene at different panicle differentiation stages under short-day conditions. Compared with HQ2, AsMADS3, AsMADS8, AsMADS11, AsMADS13, and AsMADS16 were upregulated from the initial stage to the branch differentiation stage in MSY4, while AsMADS12 was downregulated. qRT‒PCR analysis was also performed on the whole panicle differentiation stages in MSY4 under short-day conditions, the result showed that the expression levels of AsMADS9 and AsMADS11 gradually decreased. Based on the subfamily to which these genes belong, the above results indicated that AsMADS genes, especially SVP, SQUA and Mα subfamily members, regulated panicle development in MSY4 by responding to short-days. This work provides a foundation for revealing the function of the AsMADS gene family in the oat photoperiod pathway.

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 analysis of the MADS-box gene family involved in salt and waterlogging tolerance in barley (Hordeum vulgare L.)

MADS-box transcription factors are crucial members of regulatory networks underlying multiple developmental pathways and abiotic stress regulatory networks in plants. Studies on stress resistance-related functions of MADS-box genes are very limited in barley. To gain insight into this gene family and elucidate their roles in salt and waterlogging stress resistance, we performed genome-wide identification, characterization and expression analysis of MADS-box genes in barley. A whole-genome survey of barley revealed 83 MADS-box genes, which were categorized into type I (Mα, Mβ and Mγ) and type II (AP1, SEP1, AGL12, STK, AGL16, SVP and MIKC*) lineages based on phylogeny, protein motif structure. Twenty conserved motifs were determined and each HvMADS contained one to six motifs. We also found tandem repeat duplication was the driven force for HvMADS gene family expansion. Additionally, the co-expression regulatory network of 10 and 14 HvMADS genes was predicted in response to salt and waterlogging stress, and we proposed HvMADS11,13 and 35 as candidate genes for further exploration of the functions in abiotic stress. The extensive annotations and transcriptome profiling reported in this study ultimately provides the basis for MADS functional characterization in genetic engineering of barley and other gramineous crops.

Genome-wide identification and expression analysis of MADS-box transcription factors reveal their involvement in sex determination of hardy rubber tree (Eucommia ulmoides oliv.)

Eucommia ulmoides is a famous rubber-producing and medicinal tree species that produces unisexual flowers on separate individuals from the earliest stage of stamen/pistil primordium formation. To explore the genetic regulation pathway of sex in E. ulmoides, comprehensive genome-wide analyses and tissue-/sex-specific transcriptome comparisons of MADS-box transcription factors were performed for the first time in this work. Quantitative real-time PCR technique was employed to further validate the expression of genes that are assigned to floral organ ABCDE model. A total of 66 non-redundant E. ulmoides MADS-box (EuMADS) genes were identified, they were classified into Type I (M-type, 17 genes) and Type II (MIKC, 49 genes). Complex protein-motif composition, exon-intron structure and phytohormone-response cis-elements were detected in MIKC-EuMADS genes. Furthermore, 24 differentially-expressed EuMADS genes (DEGs) between male and female flowers, and two DEGs between male and female leaves were revealed. Amongst the 14 floral organ ABCDE model-related genes, there were 6 (A/B/C/E-class) and 5 (A/D/E-class) genes displayed male- and female-biased expression respectively. In particular, one B-class gene EuMADS39 and one A-class gene EuMADS65 were almost exclusively expressed in male trees, no matter in flower or leaf tissues. Collectively, these results suggested a critical role of MADS-box transcription factors in sex determination of E. ulmoides, which is conducive to decoding the molecular regulation mechanism of sex in E. ulmoides.

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.

Genome-Wide Identification, Evolution, and Expression Characterization of the Pepper (Capsicum spp.) MADS-box Gene Family

MADS domain transcription factors play roles throughout the whole lifecycle of plants from seeding to flowering and fruit-bearing. However, systematic research into MADS-box genes of the economically important vegetable crop pepper (Capsicum spp.) is still lacking. We identified 174, 207, and 72 MADS-box genes from the genomes of C. annuum, C. baccatum, and C. chinense, respectively. These 453 MADS-box genes were divided into type I (Mα, Mβ, Mγ) and type II (MIKC* and MIKC^C) based on their phylogenetic relationships. Collinearity analysis identified 144 paralogous genes and 195 orthologous genes in the three Capsicum species, and 70, 114, and 10 MADS-box genes specific to C. annuum, C. baccatum, and C. chinense, respectively. Comparative genomic analysis highlighted functional differentiation among homologous MADS-box genes during pepper evolution. Tissue expression analysis revealed three main expression patterns: highly expressed in roots, stems, leaves, and flowers (CaMADS93/CbMADS35/CcMADS58); only expressed in roots; and specifically expressed in flowers (CaMADS26/CbMADS31/CcMADS11). Protein interaction network analysis showed that type II CaMADS mainly interacted with proteins related to flowering pathway and flower organ development. This study provides the basis for an in-depth study of the evolutionary features and biological functions of pepper MADS-box genes.

Genome-wide identification and expression analysis of the MADS-box gene family during female and male flower development in Juglans mandshurica

The MADS-box gene family plays a crucial role in multiple developmental processes of plants, especially in floral organ specification and the regulation of fruit development and ripening. Juglans mandshurica is a precious fruit material whose quality and yield are determined by floral organ development. The molecular mechanism of J. mandshurica female and male flower development depending on MADS-box genes remains unclear. In our study, 67 JmMADS genes were identified and unevenly distributed on 15 of 16 J. mandshurica chromosomes. These genes were divided into two types [type I (Mα, Mγ, Mδ) and type II (MIKC)]. The gene structure and motif analyses showed that most genes belonging to the same type had similar gene structures and conserved motifs. The analysis of syntenic relationships showed that MADS-box genes in J. mandshurica, J. sigillata, and J. regia exhibited the highest homology and great collinearity. Analysis of cis-acting elements showed that JmMADS gene promoter regions contained light, stress and hormone response cis-acting elements. The gene expression patterns demonstrated that 30 and 26 JmMADS genes were specifically expressed in the female and male flowers, respectively. In addition, 12 selected genes common to J. mandshurica female and male flowers were significantly upregulated at the mature stage and were used to validate the reliability of the transcriptome data using quantitative real-time PCR. This comprehensive and systematic analysis of J. mandshurica MADS-box genes lays a foundation for future studies on MADS-box gene family functions.

Identification of MADS-Box Transcription Factors in Iris laevigata and Functional Assessment of IlSEP3 and IlSVP during Flowering

Iris laevigata is ideal for gardening and landscaping in northeast China because of its beautiful flowers and strong cold resistance. However, the short length of flowering time (2 days for individual flowers) greatly limits its applications. Molecular breeding and engineering hold high potential for producing I. laevigata of desirable flowering properties. A prerequisite is to identify and characterize key flowering control genes, the identity of which remains largely unknown in I. laevigata due to the lack of genome information. To fill this knowledge gap, we used sequencing data of the I. laevigata transcriptome to identify MADS-box gene-encoding transcription factors that have been shown to play key roles in developmental processes, including flowering. Our data revealed 41 putative MADS-box genes, which consisted of 8 type I (5 Mα and 3 Mβ, respectively) and 33 type II members (2 MIKC* and 31 MIKCC, respectively). We then selected IlSEP3 and IlSVP for functional studies and found that both are localized to the nucleus and that they interact physically in vitro. Ectopic expression of IlSEP3 in Arabidopsis resulted in early flowering (32 days) compared to that of control plants (36 days), which could be mediated by modulating the expression of FT, SOC1, AP1, SVP, SPL3, VRN1, and GA20OX. By contrast, plants overexpressing IlSVP were phenotypically similar to that of wild type. Our functional validation of IlSEP3 was consistent with the notion that SEP3 promotes flowering in multiple plant species and indicated that IlSEP3 regulates flowering in I. laevigata. Taken together, this work provided a systematic identification of MADS-box genes in I. laevigata and demonstrated that the flowering time of I. laevigata can be genetically controlled by altering the expression of key MADS-box genes.