Genome-wide identification of MADS-box gene family in sacred lotus (Nelumbo nucifera) identifies a SEPALLATA homolog gene involved in floral development

Genome-wide identification and expression pattern analysis of MIKC-Type MADS-box genes in Chionanthus retusus, an androdioecy plant

BACKGROUND

The MADS-box gene family is widely distributed in the plant kingdom, and its members typically encoding transcription factors to regulate various aspects of plant growth and development. In particular, the MIKC-type MADS-box genes play a crucial role in the determination of floral organ development and identity recognition. As a type of androdioecy plant, Chionanthus retusus have unique gender differentiation. Manifested as male individuals with only male flowers and female individuals with only bisexual flowers. However, due to the lack of reference genome information, the characteristics of MIKC-type MADS-box genes in C. retusus and its role in gender differentiation of C. retusus remain largely unknown. Therefore, it is necessary to identify and characterize the MADS-box gene family within the genome of the C. retusus.

RESULTS

In this study, we performed a genome-wide identification and analysis of MIKC-type MADS-box genes in C. retusus (2n = 2x = 46), utilizing the latest reference genome, and studied its expression pattern in individuals of different genders. As a result, we identified a total of 61 MIKC-type MADS-box genes in C. retusus. 61 MIKC-type MADS-box genes can be divided into 12 subfamilies and distributed on 18 chromosomes. Genome collinearity analysis revealed their conservation in evolution, while gene structure, domains and motif analysis indicated their conservation in structure. Finally, based on their expression patterns in floral organs of different sexes, we have identified that CrMADS45 and CrMADS60 may potentially be involved in the gender differentiation of C. retusus.

CONCLUSIONS

Our studies have provided a general understanding of the conservation and characteristics of the MIKC-type MADS-box genes family in C. retusus. And it has been demonstrated that members of the AG subfamily, CrMADS45 and CrMADS60, may play important roles in the gender differentiation of C. retusus. This provides a reference for future breeding efforts to improve flower types in C. retusus and further investigate the role of MIKC-type MADS-box genes in gender differentiation.

Genome-Wide Identification and Co-Expression Networks of WOX Gene Family in Nelumbo nucifera

WUSCHEL-related homeobox (WOX) genes are a class of plant-specific transcription factors, regulating the development of multiple tissues. However, the genomic characterizations and expression patterns of WOX genes have not been analyzed in lotus. In this study, 15 NnWOX genes were identified based on the well-annotated reference genome of lotus. According to the phylogenetic analysis, the NnWOX genes were clustered into three clades, i.e., ancient clade, intermediate clade, and WUS clade. Except for the conserved homeobox motif, we further found specific motifs of NnWOX genes in different clades and divergence gene structures, suggesting their distinct functions. In addition, two NnWOX genes in the ancient clade have conserved expression patterns and other NnWOX genes exhibit different expression patterns in lotus tissues, suggesting a low level of functional redundancy in lotus WOX genes. Furthermore, we constructed the gene co-expression networks for each NnWOX gene. Based on weighted gene co-expression network analysis (WGCNA), ten NnWOX genes and their co-expressed genes were assigned to the modules that were significantly related to the cotyledon and seed coat. We further performed RT-qPCR experiments, validating the expression levels of ten NnWOX genes in the co-expression networks. Our study reveals comprehensive genomic features of NnWOX genes in lotus, providing a solid basis for further function studies.

OfBFT genes play an essential role in the proliferate flower formation of Osmanthus fragrans

Floral organ development is one of the most vital events in flowering plants and is closely related to ornamental properties. The proliferate flower (a new branch or flower occurring in the centre of a flower) in plants is an interesting type, while the specific molecular mechanism remains largely unknown. Osmanthus fragrans 'Tianxiang Taige' has two different flower morphologies: normal flower and proliferate flower. Phenotypic observation suggested that a normal flower was composed of calyx, petal, stamen and pistil (reduced to leaf-like carpel). While in proliferate flower, the leaf-like carpel continued to grow and was replaced by a new branch. Paraffin section indicated that the re-growth of leaf carpels might be the main reason for proliferate flower formation. Transcriptome sequencing of normal and proliferate flower was performed, and the expression levels of related genes were analysed. Among the differentially expressed genes, OfBFT-a and OfBFT-b had differential expression during the proliferate flower formation process. The expression patterns revealed that both OfBFT-a and OfBFT-b were highly accumulated in carpels, and were significantly downregulated during the proliferate flower development process. Subcellular localization indicated that OfBFT-a and OfBFT-b proteins were located in the nucleus. Functional studies in 'Tianxiang Taige' and Arabidopsis showed that OfBFT-a and OfBFT-b had important roles in floral organ development, especially the proliferate flower formation process by downregulating the accumulation of AG and SEP3 homologous genes. These results may shed new light on the study of proliferate flower formation and flower morphology breeding in flowering plants.

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 and Expression Analysis of the MADS Gene Family in Tulips (Tulipa gesneriana)

To investigate the cold response mechanism and low temperature regulation of flowering in tulips, this study identified 32 MADS-box transcription factor family members in tulips based on full-length transcriptome sequencing, named TgMADS1-TgMADS32. Phylogenetic analysis revealed that these genes can be divided into two classes: type I and type II. Structural analysis showed that TgMADS genes from different subfamilies have a similar distribution of conserved motifs. Quantitative real-time PCR results demonstrated that some TgMADS genes (e.g., TgMADS3, TgMADS15, TgMADS16, and TgMADS19) were significantly upregulated in buds and stems under cold conditions, implying their potential involvement in the cold response of tulips. In summary, this study systematically identified MADS family members in tulips and elucidated their evolutionary relationships, gene structures, and cold-responsive expression patterns, laying the foundation for further elucidating the roles of these transcription factors in flowering and the cold adaptability of tulips.

Genome-Wide Identified MADS-Box Genes in Prunus campanulata 'Plena' and Theirs Roles in Double-Flower Development

The MADS-box gene family plays key roles in flower induction, floral initiation, and floral morphogenesis in flowering plants. To understand their functions in the double-flower formation of Prunus campanulata 'Plena' (hereafter referred to as PCP), which is an excellent flowering cherry cultivar, we performed genome-wide identification of the MADS-box gene family. In this study, 71 MADS-box genes were identified and grouped into the Mα, Mβ, Mγ and MIKC subfamilies according to their structures and phylogenetic relationships. All 71 MADS-box genes were located on eight chromosomes of PCP. Analysis of the cis-acting elements in the promoter region of MADS-box genes indicated that they were associated mainly with auxin, abscisic acid, gibberellin, MeJA (methyl jasmonate), and salicylic acid responsiveness, which may be involved in floral development and differentiation. By observing the floral organ phenotype, we found that the double-flower phenotype of PCP originated from petaloid stamens. The analysis of MIKC-type MADS-box genes in PCP vegetative and floral organs by qRT-PCR revealed six upregulated genes involved in petal development and three downregulated genes participating in stamen identity. Comparative analysis of petaloid stamens and normal stamens also indicated that the expression level of the AG gene (PcMADS40) was significantly reduced. Thus, we speculated that these upregulated and downregulated genes, especially PcMADS40, may lead to petaloid stamen formation and thus double flowers. This study lays a theoretical foundation for MADS-box gene identification and classification and studying the molecular mechanism underlying double flowers in other ornamental plants.

Genome-Wide Identification of the MADS-Box Gene Family during Male and Female Flower Development in Chayote (Sechium edule)

The MADS-box gene plays an important role in plant growth and development. As an important vegetable of Cucurbitaceae, chayote has great edible and medicinal value. So far, there is little molecular research on chayote, and there are no reports on the MADS-box transcription factor of chayote. In this study, the MADS-box gene family of chayote was analyzed for the first time, and a total of 70 MADS-box genes were identified, including 14 type I and 56 type II MICK MADS genes. They were randomly distributed on 13 chromosomes except for chromosome 11. The light response element, hormone response element and abiotic stress response element were found in the promoter region of 70 MADS genes, indicating that the MADS gene can regulate the growth and development of chayote, resist abiotic stress, and participate in hormone response; GO and KEGG enrichment analysis also found that SeMADS genes were mainly enriched in biological regulation and signal regulation, which further proved the important role of MADS-box gene in plant growth and development. The results of collinearity showed that segmental duplication was the main driving force of MADS gene expansion in chayote. RNA-seq showed that the expression levels of SeMADS06, SeMADS13, SeMADS26, SeMADS28, SeMADS36 and SeMADS37 gradually increased with the growth of chayote, indicating that these genes may be related to the development of root tubers of 'Tuershao'. The gene expression patterns showed that 12 SeMADS genes were specifically expressed in the male flower in 'Tuershao' and chayote. In addition, SeMADS03 and SeMADS52 may be involved in regulating the maturation of male flowers of 'Tuershao' and chayote. SeMADS21 may be the crucial gene in the development stage of the female flower of 'Tuershao'. This study laid a theoretical foundation for the further study of the function of the MADS gene in chayote in the future.

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.

SEP-like genes of Gossypium hirsutum promote flowering via targeting different loci in a concentration-dependent manner

SEP genes are famous for their function in the morphological novelty of bisexual flowers. Although the diverse functions of SEP genes were reported, only the regulatory mechanisms underlying floral organ development have been addressed. In this study, we identified SEP-like genes in Gossypium and found that SEP3 genes were duplicated in diploid cotton varieties. GhSEP4.1 and GhSEP4.2 were abundantly transcribed in the shoot apical meristem (SAM), but only GhSEP4.2 was expressed in the leaf vasculature. The expression pattern of GhSEPs in floral organs was conserved with that of homologs in Arabidopsis, except for GhSEP2 that was preponderantly expressed in ovules and fibers. The overexpression and silencing of each single GhSEP gene suggested their distinct role in promoting flowering via direct binding to GhAP1 and GhLFY genomic regions. The curly leaf and floral defects in overexpression lines with a higher expression of GhSEP genes revealed the concentration-dependent target gene regulation of GhSEP proteins. Moreover, GhSEP proteins were able to dimerize and interact with flowering time regulators. Together, our results suggest the dominant role of GhSEP4.2 in leaves to promote flowering via GhAP1-A04, and differently accumulated GhSEP proteins in the SAM alternately participate in forming the dynamic tetramer complexes to target at the different loci of GhAP1 and GhLFY to maintain reproductive growth. The regulatory roles of cotton SEP genes reveal their conserved and diversified functions.

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.

Identification of the NAC Transcription Factors and Their Function in ABA and Salinity Response in Nelumbo nucifera

Nelumbo nucifera Gaertn. is an important perennial aquatic herb that has high ornamental, edible, medicinal, and economic value, being widely distributed and used in China. The NAC superfamily (NAM, ATAF1/2, CUC2) plays critical roles in plant growth, development, and response to abiotic and biotic stresses. Though there have been a few reports about NAC genes in lotus, systematic analysis is still relatively lacking. The present study aimed to characterize all the NAC genes in the lotus and obtain better insights on the NnNACs in response to salt stress by depending on ABA signaling. Here, 97 NAC genes were identified by searching the whole lotus genome based on the raw HMM models of the conserved NAM domain and NAC domain. They were characterized by bioinformatics analysis and divided into 18 subgroups based on the phylogenetic tree. Cis-element analysis demonstrated that NAC genes are responsive to biotic and abiotic stresses, light, low temperature, and plant hormones. Meanwhile, NAC genes had tissue expression specificity. qRT-PCR analysis indicated that NAC genes could be upregulated or downregulated by NaCl treatment, ABA, and fluoridone. In addition, NAC016, NAC025, and NAC070, whose encoding genes were significantly induced by NaCl and ABA, were located in the nucleus. Further analysis showed the three NAC proteins had transcriptional activation capabilities. The co-expression network analysis reflected that NAC proteins may form complexes with other proteins to play a role together. Our study provides a theoretical basis for further research to be conducted on the regulatory mechanisms of salinity resistance in the lotus.

Studies on Lotus Genomics and the Contribution to Its Breeding

Lotus (Nelumbo nucifera), under the Nelumbonaceae family, is one of the relict plants possessing important scientific research and economic values. Because of this, much attention has been paid to this species on both its biology and breeding among the scientific community. In the last decade, the genome of lotus has been sequenced, and several high-quality genome assemblies are available, which have significantly facilitated functional genomics studies in lotus. Meanwhile, re-sequencing of the natural and genetic populations along with different levels of omics studies have not only helped to classify the germplasm resources but also to identify the domestication of selected regions and genes controlling different horticultural traits. This review summarizes the latest progress of all these studies on lotus and discusses their potential application in lotus breeding.

Comprehensive Analysis for GRF Transcription Factors in Sacred Lotus (Nelumbo nucifera)

Sacred lotus (Nelumbo nucifera) is an aquatic perennial plant with essential food, ornamental, and pharmacological value. Growth-regulating factor (GRF) is a transcription factor (TF) family that plays an important role in regulating the growth and development of plants. In this study, a comprehensive analysis of the GRF family in N. nucifera was performed, and its role in N. nucifera development was studied. A total of eight GRF genes were identified in the N. nucifera genome. Phylogenetic analysis divided the 38 GRF genes into six clades, while the NuGRFs only contained five clades. The analyses of gene structures, motifs, and cis-acting regulatory elements of the GRF gene family were performed. In addition, the chromosome location and collinearity were analyzed. The expression pattern based on transcriptomic data and real-time reverse transcription-quantitative PCR (qRT-PCR) revealed that the GRF genes were expressed in multiple organs and were abundant in actively growing tissues, and the expression levels decreased as the age of N. nucifera increased. Then, 3D structures of the NuGRF proteins were predicted by homology modeling. Finally, the subcellular localization of GRF1 was ascertained in the tobacco leaf through a vector. Therefore, this study provides a comprehensive overview of the GRF TF family in N. nucifera.

Characterization of Phytohormones and Transcriptomic Profiling of the Female and Male Inflorescence Development in Manchurian Walnut (Juglans mandshurica Maxim.)

Flowers are imperative reproductive organs and play a key role in the propagation of offspring, along with the generation of several metabolic products in flowering plants. In Juglans mandshurica, the number and development of flowers directly affect the fruit yield and subsequently its commercial value. However, owing to the lack of genetic information, there are few studies on the reproductive biology of Juglans mandshurica, and the molecular regulatory mechanisms underlying the development of female and male inflorescence remain unclear. In this study, phytohormones and transcriptomic sequencing analyses at the three stages of female and male inflorescence growth were performed to understand the regulatory functions underlying flower development. Gibberellin is the most dominant phytohormone that regulates flower development. In total, 14,579 and 7188 differentially expressed genes were identified after analyzing the development of male and female flowers, respectively, wherein, 3241 were commonly expressed. Enrichment analysis for significantly enriched pathways suggested the roles of MAPK signaling, phytohormone signal transduction, and sugar metabolism. Genes involved in floral organ transition and flowering were obtained and analyzed; these mainly belonged to the M-type MADS-box gene family. Three flowering-related genes (SOC1/AGL20, ANT, and SVP) strongly interacted with transcription factors in the co-expression network. Two key CO genes (CO3 and CO1) were identified in the photoperiod pathway. We also identified two GA20xs genes, one SVP gene, and five AGL genes (AGL8, AGL9, AGL15, AGL19, and AGL42) that contributed to flower development. The findings are expected to provide a genetic basis for the studies on the regulatory networks and reproductive biology in inflorescence development for J. mandshurica.

An ortholog of the MADS-box gene SEPALLATA3 regulates stamen development in the woody plant Jatropha curcas

Overexpression of JcSEP3 causes defective stamen development in Jatropha curcas, in which brassinosteroid and gibberellin signaling pathways may be involved. SEPALLATAs (SEPs), the class E genes of the ABCE model, are required for floral organ determination. In this study, we investigated the role of the JcSEP3 gene in floral organ development in the woody plant Jatropha curcas. Transgenic Jatropha plants overexpressing JcSEP3 displayed abnormal phenotypes such as deficient anthers and pollen, as well as free stamen filaments, whereas JcSEP3-RNA interference (RNAi) transgenic plants had no obvious phenotypic changes, suggesting that JcSEP3 is redundant with other JcSEP genes in Jatropha. Moreover, we compared the transcriptomes of wild-type plants, JcSEP3-overexpressing, and JcSEP3-RNAi transgenic plants. In the JcSEP3-overexpressing transgenic plants, we discovered 25 upregulated genes involved in anther and pollen development, as well as 12 induced genes in brassinosteroid (BR) and gibberellin (GA) signaling pathways. These results suggest that JcSEP3 directly or indirectly regulates stamen development, concomitant with the regulation of BR and GA signaling pathways. Our findings help to understand the roles of SEP genes in stamen development in perennial woody plants.

The MADS-box gene PpPI is a key regulator of the double-flower trait in peach

Double flower is an invaluable trait in ornamental peach, but the mechanism underlying its development remains largely unknown. Here, we report the roles of ABCE model genes in double flower development in peach. A total of nine ABCE regulatory genes, including eight MADS-box genes and one AP2/EREBP gene, were identified in the peach genome. Subcellular localization assay showed that all the ABCE proteins were localized in the nucleus. Four genes, PpAP1, PpAP3, PpSEP3, and PpPI, showed a difference in expression levels between single and double flowers. Ectopic overexpression of PpPI increased petal number in Arabidopsis, while transgenic lines overexpressing PpAP3 or PpSEP3 were morphologically similar to wild-type. Ectopic overexpression of PpAP1 resulted in a significant decrease in the number of basal leaves and caused early flowering. These results suggest that PpPI is likely crucial for double flower development in peach. In addition, double flowers have petaloid sepals and stamens, and single flower could occasionally change to be double flower by converting stamens to petals in peach, suggesting that the double-flower trait is likely to have evolved from an ancestral single-flower structure. Our results provide new insights into mechanisms underlying the double-flower trait in peach.