Evolution of unisexual flowers involves extreme changes in floral development. Spinach is one of the species to discern the formation and evolution of dioecy. MADS-box gene family is involved in regulation of floral organ identity and development and in many other plant developmental processes. However, there is no systematic analysis of MADS-box family genes in spinach. A comprehensive genome-wide analysis and transcriptome profiling of MADS-box genes were undertaken to understand their involvement in unisexual flower development at different stages in spinach. In total, 54 MADS-box genes found to be unevenly located across 6 chromosomes and can be divided into type I and type II genes. Twenty type I MADS-box genes are subdivided into Mα, Mβ and Mγ subgroups. While thirty-four type II SoMADSs consist of 3 MIKC*, and 31 MIKCC -type genes including sixteen floral homeotic MADS-box genes that are orthologous to the proposed Arabidopsis ABCDE model of floral organ identity determination, were identified in spinach. Gene structure, motif distribution, physiochemical properties, gene duplication and collinearity analyses for these genes are performed in detail. Promoters of both types of SoMADS genes contain mainly MeJA and ABA response elements. Expression profiling indicated that MIKCc genes exhibited more dynamic and intricate expression patterns compared to M-type genes and the majority of type-II genes AP1, SVP, and SOC1 sub-groups showed female flower-biased expression profiles, suggesting their role in carpel development, while PI showed male-biased expression throughout flower developmental stages, suggesting their role in stamen development. These results provide genomic resources and insights into spinach dioecious flower development and expedite spinach improvement.

The development of floral organs, crucial for the establishment of floral symmetry and morphology in higher plants, is regulated by MADS-box genes. In sunflower, the capitulum is comprised of ray and disc florets with various floral organs. In the sunflower long petal mutant (lpm), the abnormal disc (ray-like) floret possesses prolongated petals and degenerated stamens, resulting in a transformation from zygomorphic to actinomorphic symmetry. In this study, we investigated the effect of MADS-box genes on floral organs, particularly on petals, using WT and lpm plants as materials. Based on our RNA-seq data, 29 MADS-box candidate genes were identified, and their roles on floral organ development, especially in petals, were explored, by analyzing the expression levels in various tissues in WT and lpm plants through RNA-sequencing and qPCR. The results suggested that HaMADS3, HaMADS7, and HaMADS8 could regulate petal development in sunflower. High levels of HaMADS3 that relieved the inhibition of cell proliferation, together with low levels of HaMADS7 and HaMADS8, promoted petal prolongation and maintained the morphology of ray florets. In contrast, low levels of HaMADS3 and high levels of HaMADS7 and HaMADS8 repressed petal extension and maintained the morphology of disc florets. Their coordination may contribute to the differentiation of disc and ray florets in sunflower and maintain the balance between attracting pollinators and producing offspring. Meanwhile, Pearson correlation analysis between petal length and expression levels of MADS-box genes further indicated their involvement in petal prolongation. Additionally, the analysis of cis-acting elements indicated that these three MADS-box genes may regulate petal development and floral symmetry establishment by regulating the expression activity of HaCYC2c. Our findings can provide some new understanding of the molecular regulatory network of petal development and floral morphology formation, as well as the differentiation of disc and ray florets in sunflower.

Prunus conradinae, a valuable flowering cherry belonging to the Rosaceae family subgenus Cerasus and endemic to China, has high economic and ornamental value. However, a high-quality P. conradinae genome is unavailable, which hinders our understanding of its genetic relationships and phylogenesis, and ultimately, the possibility of mining of key genes for important traits. Herein, we have successfully assembled a chromosome-scale P. conradinae genome, identifying 31,134 protein-coding genes, with 98.22% of them functionally annotated. Furthermore, we determined that repetitive sequences constitute 46.23% of the genome. Structural variation detection revealed some syntenic regions, inversions, translocations, and duplications, highlighting the genetic diversity and complexity of Cerasus. Phylogenetic analysis demonstrated that P. conradinae is most closely related to P. campanulata, from which it diverged ~ 19.1 million years ago (Mya). P. avium diverged earlier than P. cerasus and P. conradinae. Similar to the other Prunus species, P. conradinae underwent a common whole-genome duplication event at ~ 138.60 Mya. Furthermore, 79 MADS-box members were identified in P. conradinae, accompanied by the expansion of the SHORT VEGETATIVE PHASE subfamily. Our findings shed light on the complex genetic relationships, and genome evolution of P. conradinae and will facilitate research on the molecular breeding and functions of key genes related to important horticultural and economic characteristics of subgenus Cerasus.

BACKGROUND: Lonicera japonica Thunb. is widely used in traditional Chinese medicine. Medicinal L. japonica mainly consists of dried flower buds and partially opened flowers, thus flowers are an important quality indicator. MADS-box genes encode transcription factors that regulate flower development. However, little is known about these genes in L. japonica.
RESULTS: In this study, 48 MADS-box genes were identified in L. japonica, including 20 Type-I genes (8 Mα, 2 Mβ, and 10 Mγ) and 28 Type-II genes (26 MIKCc and 2 MIKC*). The Type-I and Type-II genes differed significantly in gene structure, conserved domains, protein structure, chromosomal distribution, phylogenesis, and expression pattern. Type-I genes had a simpler gene structure, lacked the K domain, had low protein structure conservation, were tandemly distributed on the chromosomes, had more frequent lineage-specific duplications, and were expressed at low levels. In contrast, Type-II genes had a more complex gene structure; contained conserved M, I, K, and C domains; had highly conserved protein structure; and were expressed at high levels throughout the flowering period. Eleven floral homeotic MADS-box genes that are orthologous to the proposed Arabidopsis ABCDE model of floral organ identity determination, were identified in L. japonica. By integrating expression pattern and protein interaction data for these genes, we developed a possible model for floral organ identity determination.
CONCLUSIONS: This study genome-widely identified and characterized the MADS-box gene family in L. japonica. Eleven floral homeotic MADS-box genes were identified and a possible model for floral organ identity determination was also developed. This study contributes to our understanding of the MADS-box gene family and its possible involvement in floral organ development in L. japonica.

Normal floral organ development in rice is necessary for grain formation. Many MADS-box family genes that belong to ABCDE model have been widely implicated in rice flower development. The LAX1 allele encodes a plant-specific basic helix-loop-helix (bHLH) transcription factor, which is the main regulator of axillary meristem formation in rice. However, the molecular mechanisms of LAX1 allele together with MADS-box family genes underlying palea development have not been reported. We found a short palea mutant plant in a population of indica rice variety 9311 treated with cobalt 60. We report the map-based cloning and characterization of lax1-7, identified as a new mutant allele of the LAX1 locus, and the role of its wild-type allele LAX1 in rice palea development. Through complementary experiments, combined with genetic and molecular biological analyses, the function of the LAX1 allele was determined. We showed that LAX1 allele is expressed specifically in young spikelets and encodes a nucleus-localized protein. In vitro and in vivo experiments revealed that the LAX1 protein physically interacts with OsMADS1, OsMADS6 and OsMADS7. The LAX1 allele is pleiotropic for the maintenance of rice palea identity via cooperation with MADS-box genes and other traits, including axillary meristem initiation, days to heading, plant height, panicle length and spikelet fertility.

Eucommia ulmoides Oliver is a typical dioecious plant endemic to China that has great medicinal and economic value. Here, we report a high-quality chromosome-level female genome of E. ulmoides obtained by PacBio and Hi-C technologies. The size of the female genome assembly was 1.01 Gb with 17 pseudochromosomes and 31,665 protein coding genes. In addition, Hi-C technology was used to reassemble the male genome released in 2018. The reassembled male genome was 1.24 Gb with the superscaffold N50 (48.30 Mb), which was increased 25.69 times, and the number of predicted genes increased by 11,266. Genome evolution analysis indicated that E. ulmoides has undergone two whole-genome duplication events before the divergence of female and male, including core eudicot γ whole-genome triplication event (γ-WGT) and a recent whole genome duplication (WGD) at approximately 27.3 million years ago (Mya). Based on transcriptome analysis, EuAP3 and EuAG may be the key genes involved in regulating the sex differentiation of E. ulmoides. Pathway analysis showed that the high expression of ω-3 fatty acid desaturase coding gene EU0103017 was an important reason for the high α-linolenic acid content in E. ulmoides. The genome of female and male E. ulmoides presented here is a valuable resource for the molecular biological study of sex differentiation of E. ulmoides and also will provide assistance for the breeding of superior varieties.

MADS-box is a class of transcriptional regulators that are ubiquitous in plants and plays important roles in the process of plant growth and development. Identification and analysis of blueberry MADS-box genes can lay a foundation for their function investigations. In the present study, 249 putative MADS-box genes were identified in the blueberry genome. Those MADS-box genes were distributed on 47 out of 48 chromosomes. The phylogenetic and evolutionary analyses showed that blueberry MADS-box genes were divided into 131 type I members and 118 type II members. The type I genes contained an average of 1.89 exons and the type II genes contained an average of 7.83 exons. Motif analysis identified 15 conserved motifs, of which 4 were related to the MADS domain and 3 were related to the K-box domain. A variety of cis-acting elements were found in the promoter region of the blueberry MADS-box gene, indicating that the MADS-box gene responded to various hormones and environmental alterations. A total of 243 collinear gene pairs were identified, most of which had a Ka/Ks value of less than 1. Nine genes belonging to SEP, AP3/PI, and AGL6 subfamilies were screened based on transcriptomic data. The expression patterns of those nine genes were also verified using quantitative PCR, suggesting that VcMADS6, VcMADS35, VcMADS44, VcMADS58, VcMADS125, VcMADS188, and VcMADS212 had potential functions in blueberry fruit ripening. The results of this study provide references for an in-depth understanding of the biological function of the blueberry MADS-box genes and the mechanism of blueberry fruit ripening.

MADS-box transcription factors play an important role in regulating floral organ development and participate in environmental responses. To date, the MADS-box gene family has been widely identified in Brassica rapa (B. rapa), Brassica oleracea (B. oleracea), and Brassica napus (B. napus); however, there are no analogous reports in Brassica nigra (B. nigra), Brassica juncea (B. juncea), and Brassica carinata (B. carinata). In this study, a whole-genome survey of the MADS-box gene family was performed for the first time in the triangle of U species, and a total of 1430 MADS-box genes were identified. Based on the phylogenetic relationship and classification of MADS-box genes in Arabidopsis thaliana (A. thaliana), 1430 MADS-box genes were categorized as M-type subfamily (627 genes), further divided into Mα, Mβ, Mγ, and Mδ subclades, and MIKC-type subfamily (803 genes), further classified into 35 subclades. Gene structure and conserved protein motifs of MIKC-type MADS-box exhibit diversity and specificity among different subclades. Comparative analysis of gene duplication events and syngenic gene pairs among different species indicated that polyploidy is beneficial for MIKC-type gene expansion. Analysis of transcriptome data within diverse tissues and stresses in B. napus showed tissue-specific expression of MIKC-type genes and a broad response to various abiotic stresses, particularly dehydration stress. In addition, four representative floral organ mutants (wtl, feml, aglf-2, and aglf-1) in the T0 generation were generated by editing four AGAMOUS (BnaAG) homoeologs in B. napus that enriched the floral organ variant phenotype. In brief, this study provides useful information for investigating the function of MADS-box genes and contributes to revealing the regulatory mechanisms of floral organ development in the genetic improvement of new varieties.

Japanese cultivated gentians are perennial plants that flower in early summer to late autumn in Japan, depending on the cultivar. Several flowering-related genes, including GtFT1 and GtTFL1, are known to be involved in regulating flowering time, but many such genes remain unidentified. In this study, we obtained transcriptome profiling data using the Gentiana triflora cultivar \'Maciry\', which typically flowers in late July. We conducted deep RNA sequencing analysis using gentian plants grown under natural field conditions for three months before flowering. To investigate diurnal changes, the plants were sampled at 4 h intervals over 24 h. Using these transcriptome data, we determined the expression profiles of leaves based on homology searches against the Flowering-Interactive Database of Arabidopsis. In particular, we focused on transcription factor genes, belonging to the BBX and MADS-box families, and analyzed their developmental and diurnal variation. The expression levels of representative BBX genes were also analyzed under long- and short-day conditions using in-vitro-grown seedlings, and the expression patterns of some BBX genes differed. Clustering analysis revealed that the transcription factor genes were coexpressed with GtFT1. Overall, these expression profiles will facilitate further analysis of the molecular mechanisms underlying the control of flowering time in gentians.

Sex determination in dioecious plants has been broadly and progressively studied with the blooming of genome sequencing and editing techniques. This provides us with a great opportunity to explore the evolution and genetic mechanisms underlining the sex-determining system in dioecious plants. In this study, comprehensively reviewing advances in sex-chromosomes, sex-determining genes, and floral MADS-box genes in dioecious plants, we proposed a convergent model that governs plant dioecy across divergent species using a cascade regulation pathway connecting sex-determining genes and MADS-box genes e.g., B-class genes. We believe that this convergent mechanism of sex determination in dioecious plants will shed light on our understanding of gene regulation and evolution of plant dioecy. Perspectives concerning the evolutionary pathway of plant dioecy are also suggested.