MADS-Box gene diversity in seed plants 300 million years ago

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.

Systematic analysis of MADS-box gene family in the U's triangle species and targeted mutagenesis of BnaAG homologs to explore its role in floral organ identity in Brassica napus

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.

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

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

Genome-Wide Analysis of the MADS-Box Gene Family in Holoparasitic Plants (Balanophora subcupularis and Balanophora fungosa var. globosa)

MADS-box is an important transcription factor family that is involved in the regulation of various stages of plant growth and development, especially flowering regulation and flower development. Being a holoparasitic plant, the body structure of Balanophoraceae has changed dramatically over time, and its vegetative and reproductive organs have been extensively modified, with rudimentary flower organs. Meanwhile, extraordinary gene losses have been identified in holoparasitic plants compared with autotrophs. Our study reveals that the MADS-box gene family contracted sharply in Balanophora subcupularis and Balanophora fungosa var. globosa, and some subfamilies were lost, exhibiting reduced redundancy in both. The genes that functioned in the transition from the vegetative to floral production stages suffered a significant loss, but the ABCE model genes remained intact. We further investigated genes related to flowering regulation in B. subcupularis and B. fungosa var. globosa, vernalization and autonomous ways of regulating flowering time remained comparatively integrated, while genes in photoperiod and circadian clock pathways were almost lost. Convergent gene loss in flowering regulation occurred in Balanophora and another holoparasitic plant Sapria himalayana (Rafflesiaceae). The genome-wide analysis of the MADS-box gene family in Balanophora species provides valuable information for understanding the classification, gene loss pattern, and flowering regulation mechanism of MADS-box gene family in parasitic plants.

Functional identification of the different regions in B-class floral homeotic MADS-box proteins IiAP3 and IiPI from Isatis indigotica

APETALA3 (AP3) and PISTILLATA (PI) are B-class MADS-box floral homeotic genes of Arabidopsis and are involved in specifying the identity of petals and stamens. In the present work, IiAP3 and IiPI, the respective orthologous genes of AP3 and PI, were cloned from Isatis indigotica. By expressing in ap3-6 and pi-1 homozygous mutant and in wild-type Arabidopsis under the control of AP3 promoter or CaMV 35S promoter, we demonstrated that IiAP3 and IiPI were functionally equivalent to AP3 and PI of Arabidopsis. Referring to previous reports and the research results in the present work, expression patterns of AP3 and PI homologs are not the same in different angiosperms possessing diverse floral structures. It suggests that the alterations in expression may contribute to the changing morphology of flowers. To further determine the relationship between IiAP3 and IiPI, the coding sequences of the different structural regions in these two proteins were swapped with each other, and the data collected from transgenic Arabidopsis plants of the chimeric constructs suggested that MADS domain was irreplaceable for the function of IiAP3, K domain of IiAP3 was involved in specifying the identity of stamens, K domain of IiPI was mainly related to the formation of petals, and C-terminal region of IiPI was involved in characterization of stamens. In addition, a complete KC region of these two proteins was more effective in phenotypic complementation of the mutants.

Transcriptional Regulation of Pine Male and Female Cone Initiation and Development: Key Players Identified Through Comparative Transcriptomics

With long reproductive timescales, large complex genomes, and a lack of reliable reference genomes, understanding gene function in conifers is extremely challenging. Consequently, our understanding of which genetic factors influence the development of reproductive structures (cones) in monoecious conifers remains limited. Genes with inferred roles in conifer reproduction have mostly been identified through homology and phylogenetic reconstruction with their angiosperm counterparts. We used RNA-sequencing to generate transcriptomes of the early morphological stages of cone development in the conifer species Pinus densiflora and used these to gain a deeper insight into the transcriptional changes during male and female cone development. Paired-end Illumina sequencing was used to generate transcriptomes from non-reproductive tissue and male and female cones at four time points with a total of 382.82 Gbp of data generated. After assembly and stringent filtering, a total of 37,164 transcripts were retrieved, of which a third were functionally annotated using the Mercator plant pipeline. Differentially expressed gene (DEG) analysis resulted in the identification of 172,092 DEGs in the nine tissue types. This, alongside GO gene enrichment analyses, pinpointed transcripts putatively involved in conifer reproductive structure development, including co-orthologs of several angiosperm flowering genes and several that have not been previously reported in conifers. This study provides a comprehensive transcriptome resource for male and early female cone development in the gymnosperm species Pinus densiflora. Characterisation of this resource has allowed the identification of potential key players and thus provides valuable insights into the molecular regulation of reproductive structure development in monoecious conifers.

Circlize package in R and Analytic Hierarchy Process (AHP): Contribution values of ABCDE and AGL6 genes in the context of floral organ development

The morphological diversity of floral organs can largely be attributed to functional divergence in the MADS-box gene family. Nonetheless, research based on the ABCDE model has yet to conclusively determine whether the AGAMOUS-LIKE 6 (AGL6) subgroup has a direct influence on floral organ development. In the current study, the ABCDE model was used to quantify the contributions of ABCDE and AGL6 genes in the emergence of floral organs. We determined that the flower formation contribution values of the ABCDE and AGL6 genes were as follows: A gene, 0.192; B gene, 0.231; CD gene, 0.192; E gene, 0.385; and AGL6, 0.077. As AGL6 does not directly influence floral structure formation, the contribution value of AGL6 to flower formation was low. Furthermore, the gradient values of the floral organs were as follows: sepals, 0.572; petals, 1.606; stamens, 2.409; and carpels, 2.288. We also performed detailed analysis of the ABCDE and AGL6 genes using the Circlize package in R. Our results suggest that these genes likely emerged in one of two orders: 1) B genes→CD genes→AGL6→E genes→A genes; or 2) B genes→CD genes→AGL6/E genes→A genes. We use the analytic hierarchy process (AHP) to prove the contribution values and gradient values of floral organs. This is the first study to understand the contribution values of ABCDE and AGL6 genes using the AHP and the Circlize package in R.

A Novel Embryo Phenotype Associated With Interspecific Hybrid Weakness in Rice Is Controlled by the MADS-Domain Transcription Factor OsMADS8

A novel hybrid weakness gene, DTE9, associated with a dark tip embryo (DTE) trait, was observed in CR6078, an introgression line derived from a cross between the Oryza sativa spp. japonica "Hwayeong" (HY) and the wild relative Oryza rufipogon. CR6078 seeds exhibit protruding embryos and flowers have altered inner floral organs. DTE9 was also associated with several hybrid weakness symptoms including decreased grain weight. Map-based cloning and transgenic approaches revealed that DTE9 is an allele of OsMADS8, a MADS-domain transcription factor. Genetic analysis indicated that two recessive complementary genes were responsible for the expression of the DTE trait. No sequence differences were observed between the two parental lines in the OsMADS8 coding region; however, numerous single nucleotide polymorphisms were detected in the promoter and intronic regions. We generated overexpression (OX) and RNA interference (RNAi) transgenic lines of OsMADS8 in HY and CR6078, respectively. The OsMADS8-OX lines showed the dark tip embryo phenotype, whereas OsMADS8-RNAi recovered the normal embryo phenotype. Changes in gene expression, including of ABCDE floral homeotic genes, were observed in the OsMADS8-OX and OsMADS8-RNAi lines. Overexpression of OsMADS8 led to decreased expression of OsEMF2b and ABA signaling-related genes including OsVP1/ABI3. HY seeds showed higher ABA content than CR6078 seeds, consistent with OsMADS8/DTE9 regulating the expression of genes related ABA catabolism in CR6078. Our results suggest that OsMADS8 is critical for floral organ determination and seed germination and that these effects are the result of regulation of the expression of OsEMF2b and its role in ABA signaling and catabolism.

The role of pollination in controlling Ginkgo biloba ovule development

Generally, in gymnosperms, pollination and fertilization events are temporally separated and the developmental processes leading the switch from ovule integument into seed coat are still unknown. The single ovule integument of Ginkgo biloba acquires the typical characteristics of the seed coat long before the fertilization event. In this study, we investigated whether pollination triggers the transformation of the ovule integument into the seed coat. Transcriptomics and metabolomics analyses performed on ovules just prior and after pollination lead to the identification of changes occurring in Ginkgo ovules during this specific time. A morphological atlas describing the developmental stages of ovule development is presented. The metabolic pathways involved in the lignin biosynthesis and in the production of fatty acids are activated upon pollination, suggesting that the ovule integument starts its differentiation into a seed coat before the fertilization. Omics analyses allowed an accurate description of the main changes that occur in Ginkgo ovules during the pollination time frame, suggesting the crucial role of the pollen arrival on the progression of ovule development.

Genome-Wide Identification MIKC-Type MADS-Box Gene Family and Their Roles during Development of Floral Buds in Wheel Wingnut (Cyclocarya paliurus)

MADS-box transcription factors (TFs) have fundamental roles in regulating floral organ formation and flowering time in flowering plants. In order to understand the function of MIKC-type MADS-box family genes in Cyclocarya paliurus (Batal.) Iljinskaja, we first implemented a genome-wide analysis of MIKC-type MADS-box genes in C. paliurus. Here, the phylogenetic relationships, chromosome location, conserved motif, gene structure, promoter region, and gene expression profile were analyzed. The results showed that 45 MIKC-type MADS-box were divided into 14 subfamilies: BS (3), AGL12 (1), AP3-PI (3), MIKC* (3), AGL15 (3), SVP (5), AGL17 (2), AG (3), TM8 (1), AGL6 (2), SEP (5), AP1-FUL (6), SOC1 (7), and FLC (1). The 43 MIKC-type MADS-box genes were distributed unevenly in 14 chromosomes, but two members were mapped on unanchored scaffolds. Gene structures were varied in the same gene family or subfamily, but conserved motifs shared similar distributions and sequences. The element analysis in promoters’ regions revealed that MIKC-type MADS-box family genes were associated with light, phytohormone, and temperature responsiveness, which may play important roles in floral development and differentiation. The expression profile showed that most MIKC-type MADS-box genes were differentially expressed in six tissues (specifically expressed in floral buds), and the expression patterns were also visibly varied in the same subfamily. CpaF1st24796 and CpaF1st23405, belonging to AP3-PI and SEP subfamilies, exhibited the high expression levels in PA-M and PG-F, respectively, indicating their functions in presenting heterodichogamy. We further verified the MIKC-type MADS-box gene expression levels on the basis of transcriptome and qRT-PCR analysis. This study would provide a theoretical basis for classification, cloning, and regulation of flowering mechanism of MIKC-type MADS-box genes in C. paliurus.

Revealing the developmental dynamics in male strobilus transcriptome of Gnetum luofuense using nanopore sequencing technology

Gnetum is a pantropical distributed gymnosperm genus. As being dioecious, Gnetum species apply female and male strobili to attract and provide nutrition to insect pollinators. Due to its unique gross morphology, a Gnetum male strobilus receives much attention in previous taxonomic and evolutionary studies. However, underlying molecular mechanisms that control male strobilus development and pollination adaptation have not been well studied. In the present study, nine full-length transcriptomes were sequenced from three developmental stages of the G. luofuense male strobili using Oxford Nanopore Technologies. In addition, weighted gene co-expression network analysis (WGCNA), and RT-qPCR analysis were performed. Our results show that a total of 3138 transcription factors and 466 long non-coding RNAs (lncRNAs) were identified, and differentially expressed lncRNAs and TFs reveal a dynamic pattern during the male strobilus development. Our results show that MADS-box and Aux/IAA TFs were differentially expressed at the three developmental stages, suggesting their important roles in the regulation of male strobilus development of G. luofuense. Results of WGCNA analysis and annotation of differentially expressed transcripts corroborate that the male strobilus development of G. luofuense is closely linked to plant hormone changes, photosynthesis, pollination drop secretion and reproductive organ defense. Our results provide a valuable resource for understanding the molecular mechanisms that drive organ evolution and pollination biology in Gnetum.

Genome-wide analysis of the MADS-Box gene family in Chrysanthemum

MADS-box family transcription factors play key roles in various developmental processes in plants. Here, we identified 108 MADS-box genes in the genome of chrysanthemum (Chrysanthemum nankingense). We classified these genes based on their phylogenetic relationships with MADS-box genes in Arabidopsis thaliana and lettuce (Lactuca sativa). Type I genes were subdivided into classes Mα (19 genes), Mβ (12 genes), and Mγ (10 genes), and type II genes were subdivided into classes MIKC^C (64 genes) and MIKC* (3 genes). The MIKC^C class genes were further divided into 16 subclasses that included genes described in the ABCDE flower development model. Each group of MADS-box genes showed a specific pattern of conserved protein motifs and exon-intron structure. We analyzed the expression levels of each MADS-box gene in root, stem, leaf, flower bud, disc floret, and ray floret tissues. Subfamilies AGL18, FLC, and SVP contained more members in chrysanthemum. The asterid-specific TM8 subfamily and eleven Asteraceae Specific-MADS CnMADS genes were present in chrysanthemum. Chrysanthemum is the lacking members of the AGL15 subfamily. Among the genes responsible for the ABCDE model, B-class genes were expanded in chrysanthemum with three AP3 and four PI genes. One AP3 homolog functions in marginal ray floret development, whereas the two other AP3 homologs function in the development of the central disc floret. Two of the four PI genes are expressed in chrysanthemum, specifically in both types of florets. The results of this study lay the foundation for further studies of the roles of MADS-box genes in flower development in chrysanthemum and of the evolution of MADS-box genes in plants.

The genomic architecture of the sex-determining region and sex-related metabolic variation in Ginkgobiloba

Sex differences and evolutionary differences are critical biological issues. Ginkgo is an ancient lineage of dioecious gymnosperms with special value for studying the mechanism of sex determination in plants. However, the major genetic basic underlying sex chromosomes remains to be uncovered. In this study, we identify the sex-determining region of Ginkgo and locate it to the area from megabases 48 to 75 on chromosome 2. We find that the male sex-determining region of Ginkgo contains more than 200 genes, including four MADS-box genes, demonstrating that the Ginkgo sex determination system is of the XY type. We also find that genetic sex differences result in specialized flavonoid metabolism and regulation in each sex. These findings establish a foundation for revealing the molecular mechanism of sexual dimorphism and promoting the development of the Ginkgo industry.

Evolution of SHORT VEGETATIVE PHASE (SVP) genes in Rosaceae: Implications of lineage-specific gene duplication events and function diversifications with respect to their roles in processes other than bud dormancy

MADS-box genes that are homologous to Arabidopsis SHORT VEGETATIVE PHASE (SVP) have been shown to play key roles in the regulation of bud dormancy in perennial species, particularly in the deciduous fruit trees of Rosaceae. However, their evolutionary profiles in Rosaceae have not yet been analyzed systematically. Here, The SVP genes were found to be significantly expanded in Rosaceae when compared with annual species from Brassicaceae. Phylogenetic analysis showed that Rosaceae SVP genes could be classified into five clades, namely, SVP1, SVP2-R1, SVP2-R2, SVP2-R3 and SVP3. The SVP1 clade genes were retained in most of the species, whereas the SVP2-R2 and SVP2-R3 clades were found to be Maleae- and Amygdaleae-specific (Both of the lineages belong to Amygdaloideae), respectively, and SVP2-R1 was Rosoideae-specific in Rosaceae. Furthermore, 10 lineage-specific gene duplication (GD) events (GD1-10) were proposed for the expansion of SVP genes, suggesting that the expansion and divergence of Rosaceae SVP genes were mainly derived by lineage-specific manner during evolution. Moreover, tandem and segmental duplications were the major reasons for the expansion of SVP genes, and interestingly, tandem duplications, a well-known evolutionary feature of SVP genes, were found to be mainly Amygdaloideae-specific. Sequence alignment, selection pressure, and cis-acting element analysis suggested large functional innovations and diversification of SVP genes in different lineages of Rosaceae. Finally, the different growth cycle of Rosa multiflora and their novel expression patterns of RmSVP genes provided new insights into the functional diversification of SVP genes in terms of their roles in processes other than bud dormancy.

The Floral Repressor GmFLC-like Is Involved in Regulating Flowering Time Mediated by Low Temperature in Soybean

Soybean is an important crop that is grown worldwide. Flowering time is a critical agricultural trait determining successful reproduction and yields. For plants, light and temperature are important environmental factors that regulate flowering time. Soybean is a typical short-day (SD) plant, and many studies have elucidated the fine-scale mechanisms of how soybean responds to photoperiod. Low temperature can delay the flowering time of soybean, but little is known about the detailed mechanism of how temperature affects soybean flowering. In this study, we isolated GmFLC-like from soybean, which belongs to the FLOWERING LOCUS C clade of the MADS-box family and is intensely expressed in soybean leaves. Heterologous expression of GmFLC-like results in a delayed-flowering phenotype in Arabidopsis. Additional experiments revealed that GmFLC-like is involved in long-term low temperature-triggered late flowering by inhibiting FT gene expression. In addition, yeast one-hybrid, dual-luciferase reporter assay, and electrophoretic mobility shift assay revealed that the GmFLC-like protein could directly repress the expression of FT2a by physically interacting with its promoter region. Taken together, our results revealed that GmFLC-like functions as a floral repressor involved in flowering time during treatments with various low temperature durations. As the only the FLC gene in soybean, GmFLC-like was meaningfully retained in the soybean genome over the course of evolution, and this gene may play an important role in delaying flowering time and providing protective mechanisms against sporadic and extremely low temperatures.

Pollination Drop Proteome and Reproductive Organ Transcriptome Comparison in Gnetum Reveals Entomophilous Adaptation

Gnetum possesses morphologically bisexual but functionally unisexual reproductive structures that exude sugary pollination drops to attract insects. Previous studies have revealed that the arborescent species (G. gnemon L.) and the lianoid species (G. luofuense C.Y.Cheng) possess different pollination syndromes. This study compared the proteome in the pollination drops of these two species using label-free quantitative techniques. The transcriptomes of fertile reproductive units (FRUs) and sterile reproductive units (SRUs) for each species were furthermore compared using Illumina Hiseq sequencing, and integrated proteomic and transcriptomic analyses were subsequently performed. Our results show that the differentially expressed proteins between FRUs and SRUs were involved in carbohydrate metabolism, the biosynthesis of amino acids and ovule defense. In addition, the differentially expressed genes between the FRUs and SRUs (e.g., MADS-box genes) were engaged in reproductive development and the formation of pollination drops. The integrated protein-transcript analyses revealed that FRUs and their exudates were relatively conservative while the SRUs and their exudates were more diverse, probably functioning as pollinator attractants. The evolution of reproductive organs appears to be synchronized with changes in the pollination drop proteome of Gnetum, suggesting that insect-pollinated adaptations are not restricted to angiosperms but also occur in gymnosperms.

Genome-Wide Identification, Gene Structure and Expression Analysis of the MADS-Box Gene Family Indicate Their Function in the Development of Tobacco (Nicotiana tabacum L.)

MADS-box genes play a pivotal role in various processes, including floral and seed development, controlling flowering time, regulation of fruits ripening, and respond to abiotic and biotic stressors in planta. Tobacco (Nicotiana tabacum) has been widely used as a model plant for analyzing the gene function, however, there has been less information on the regulation of flowering, and the associated genes. In the present study, a total of 168 NtMADS-box genes were identified from tobacco, and their phylogenetic relationship, chromosome locations, and gene structures were further analyzed. NtMADS-box genes can be clustered into four sub-families of Mα, Mγ, MIKC*, and MIKCC. A total of 111 NtMADS-box genes were distributed on 20 chromosomes, and 57 NtMADS-box genes were located on the unanchored scaffolds due to the complex and incomplete assembly of the tobacco genome. Expression profiles of NtMADS-box genes by microarray from 23 different tissues indicated that members in different NtMADS-box gene subfamilies might play specific roles in the growth and flower development, and the transcript levels of 24 NtMADS-box genes were confirmed by quantitative real-time PCR. Importantly, overexpressed NtSOC1/NtMADS133 could promote early flowering and dwarfism in transgenic tobacco plants. Therefore, our findings provide insights on the characterization of NtMADS-box genes to further study their functions in plant development.

Genome-wide Analysis of the MADS-Box Gene Family in Watermelon

MADS-box genes comprise a family of transcription factors that function in the growth and development of plants. To obtain insights into their evolution in watermelon (Citrullus lanatus), we carried out a genome-wide analysis and identified 39 MADS-box genes. These genes were classified into MIKC^c (25), MIKC^*(3), Mα (5), Mβ (3), and Mγ (3) clades according to their phylogenetic relationship with Arabidopsis thaliana and Cucumis sativus; moreover, these 25 genes in the MIKC clade could be classified into 13 subfamilies, and the Flowering Locus C (FLC) subfamily is absent in watermelon. Analysis of the conserved gene motifs showed similar motifs among clades. Continuing chromosomal localizations analysis indicated that MADS-box genes were distributed across 11 chromosomes in watermelon, and these genes were conditioned to be differentially expressed during plant growth and development. This research provides information that will aid further investigations into the evolution of the MADS-box gene family in plants.

Origination and selection of ABCDE and AGL6 subfamily MADS-box genes in gymnosperms and angiosperms

BACKGROUND

The morphological diversity of flower organs is closely related to functional divergence within the MADS-box gene family. Bryophytes and seedless vascular plants have MADS-box genes but do not have ABCDE or AGAMOUS-LIKE6 (AGL6) genes. ABCDE and AGL6 genes belong to the subgroup of MADS-box genes. Previous works suggest that the B gene was the first ABCDE and AGL6 genes to emerge in plant but there are no mentions about the probable origin time of ACDE and AGL6 genes. Here, we collected ABCDE and AGL6 gene 381 protein sequences and 361 coding sequences from gymnosperms and angiosperms and reconstructed a complete Bayesian phylogeny of these genes. In this study, we want to clarify the probable origin time of ABCDE and AGL6 genes is a great help for understanding the role of the formation of the flower, which can decipher the forming order of MADS-box genes in the future.

RESULTS

These genes appeared to have been under purifying selection and their evolutionary rates are not significantly different from each other. Using the Bayesian evolutionary analysis by sampling trees (BEAST) tool, we estimated that: the mutation rate of the ABCDE and AGL6 genes was 2.617 × 10-3 substitutions/site/million years, and that B genes originated 339 million years ago (MYA), CD genes originated 322 MYA, and A genes shared the most recent common ancestor with E/AGL6 296 MYA, respectively.

CONCLUSIONS

The phylogeny of ABCDE and AGL6 genes subfamilies differed. The APETALA1 (AP1 or A gene) subfamily clustered into one group. The APETALA3/PISTILLATA (AP3/PI or B genes) subfamily clustered into two groups: the AP3 and PI clades. The AGAMOUS/SHATTERPROOF/SEEDSTICK (AG/SHP/STK or CD genes) subfamily clustered into a single group. The SEPALLATA (SEP or E gene) subfamily in angiosperms clustered into two groups: the SEP1/2/4 and SEP3 clades. The AGL6 subfamily clustered into a single group. Moreover, ABCDE and AGL6 genes appeared in the following order: AP3/PI → AG/SHP/STK → AGL6/SEP/AP1. In this study, we collected candidate sequences from gymnosperms and angiosperms. This study highlights important events in the evolutionary history of the ABCDE and AGL6 gene families and clarifies their evolutionary path.

Evolution and functional divergence of MADS-box genes in Pyrus

MADS-box transcription factors widely regulate all aspects of plant growth including development and reproduction. Although the MADS-box gene family genes have been extensively characterized in many plants, they have not been studied in closely related species. In this study, 73 and 74 MADS-box genes were identified in European pear (Pyrus communis) and Chinese pear (Pyrus bretschneideri), respectively. Based on the phylogenetic relationship, these genes could be clustered into five groups (Mα, Mβ, Mr, MIKC^C, MIKC*) and the MIKC^C group was further categorized into 10 subfamilies. The distribution of MADS-box genes on each chromosome was significantly nonrandom. Thirty-seven orthologs, twenty-five PcpMADS (P. communis MADS-box) paralogs and nineteen PbrMADS (P. bretschneideri MADS-box) paralogs were predicted. Among these paralogous genes, two pairs arose from tandem duplications (TD), nineteen from segmental duplication (SD) events and twenty-three from whole genome duplication (WGD) events, indicating SD/WGD events led to the expansion of MADS-box gene family. The MADS-box genes expression profiles in pear fruits indicated functional divergence and neo-functionalization or sub-functionalization of some orthologous genes originated from a common ancestor. This study provided a useful reference for further analysis the mechanisms of species differentiation and biodiversity formation among closely related species.

MIKCC-type MADS-box genes in Rosa chinensis: the remarkable expansion of ABCDE model genes and their roles in floral organogenesis

MIKC^C-type MADS-box (MIKC^C) genes encode transcription factors that have crucial roles in controlling floral organogenesis and flowering time in plants. Although this gene family has been well characterized in many plant species, its evolutionary and comprehensive functional analysis in rose is lacking. In this study, 58 non-redundant MIKC^C uni-transcripts were extensively identified from rose transcriptomes. Phylogenetic analysis placed these genes into 12 clades with their Arabidopsis and strawberry counterparts, and revealed that ABCDE model (including AP1/FUL, AP3/PI, AG, and SEP clades), and SOC1 and AGL6 clade genes have remarkably expanded in Rosa chinensis, whereas genes from the FLC and AGL17 clades were undetectable. Sequence alignments suggest that the AP3/PI clade may contribute to more specific functions in rose due to a high variation of amino acid residues within its MADS-box domains. A comparative analysis of gene expression in specific floral organ differentiation stages and floral organs between R. chinensis cv. Old Blush and the closely related mutant genotype R. chinensis cv. Viridiflora (floral organs mutated into leaf-like structures) further revealed the roles of ABCDE model genes during floral organogenesis in rose. Analysis of co-expression networks provided an overview of the regulatory mechanisms of rose MIKC^C genes and shed light on both the prominent roles of AP3/PI clade genes in floral organogenesis and the roles of RcAGL19, RcAGL24, and RcSOC1 in regulating floral transition in rose. Our analyses provide an overall insight of MIKC^C genes in rose and their potential roles in floral organogenesis.

Evolutionary Analysis of MIKCc-Type MADS-Box Genes in Gymnosperms and Angiosperms

MIKC^c-type MADS-box genes encode transcription factors that control floral organ morphogenesis and flowering time in flowering plants. Here, in order to determine when the subfamilies of MIKC^c originated and their early evolutionary trajectory, we sampled and analyzed the genomes and large-scale transcriptomes representing all the orders of gymnosperms and basal angiosperms. Through phylogenetic inference, the MIKC^c-type MADS-box genes were subdivided into 14 monophyletic clades. Among them, the gymnosperm orthologs of AGL6, SEP, AP1, GMADS, SOC1, AGL32, AP3/PI, SVP, AGL15, ANR1, and AG were identified. We identified and characterized the origin of a novel subfamily GMADS within gymnosperms but lost orthologs in monocots and Brassicaceae. ABCE model prototype genes were relatively conserved in terms of gene number in gymnosperms, but expanded in angiosperms, whereas SVP, SOC1, and GMADS had dramatic expansions in gymnosperms but conserved in angiosperms. Our results provided the most detailed evolutionary history of all MIKC^c gene clades in gymnosperms and angiosperms. We proposed that although the near complete set of MIKC^c genes had evolved in gymnosperms, the duplication and expressional transition of ABCE model MIKC^c genes in the ancestor of angiosperms triggered the first flower.

Genome-Wide Identification of the MIKC-Type MADS-Box Gene Family in Gossypium hirsutum L. Unravels Their Roles in Flowering

Cotton is one of the major world oil crops. Cottonseed oil meets the increasing demand of fried food, ruminant feed, and renewable bio-fuels. MADS intervening keratin-like and C-terminal (MIKC)-type MADS-box genes encode transcription factors that have crucial roles in various plant developmental processes. Nevertheless, this gene family has not been characterized, nor its functions investigated, in cotton. Here, we performed a comprehensive analysis of MIKC-type MADS genes in the tetraploid Gossypium hirsutum L., which is the most widely cultivated cotton species. In total, 110 GhMIKC genes were identified and phylogenetically classified into 13 subfamilies. The Flowering locus C (FLC) subfamily was absent in the Gossypium hirsutum L. genome but is found in Arabidopsis and Vitis vinifera L. Among the genes, 108 were distributed across the 13 A and 12 of the D genome's chromosomes, while two were located in scaffolds. GhMIKCs within subfamilies displayed similar exon/intron characteristics and conserved motif compositions. According to RNA-sequencing, most MIKC genes exhibited high flowering-associated expression profiles. A quantitative real-time PCR analysis revealed that some crucial MIKC genes determined the identities of the five flower organs. Furthermore, the overexpression of GhAGL17.9 in Arabidopsis caused an early flowering phenotype. Meanwhile, the expression levels of the flowering-related genes CONSTANS (CO), LEAFY (LFY) and SUPPRESSOR OF OVEREXPRESSION OF CONSTANS1 (SOC1) were significantly increased in these lines. These results provide useful information for future studies of GhMIKCs' regulation of cotton flowering.

The MADS Box Genes ABS, SHP1, and SHP2 Are Essential for the Coordination of Cell Divisions in Ovule and Seed Coat Development and for Endosperm Formation in Arabidopsis thaliana

Seed formation is a pivotal process in plant reproduction and dispersal. It begins with megagametophyte development in the ovule, followed by fertilization and subsequently coordinated development of embryo, endosperm, and maternal seed coat. Two closely related MADS-box genes, SHATTERPROOF 1 and 2 (SHP1 and SHP2) are involved in specifying ovule integument identity in Arabidopsis thaliana. The MADS box gene ARABIDOPSIS BSISTER (ABS or TT16) is required, together with SEEDSTICK (STK) for the formation of endothelium, part of the seed coat and innermost tissue layer formed by the maternal plant. Little is known about the genetic interaction of SHP1 and SHP2 with ABS and the coordination of endosperm and seed coat development. In this work, mutant and expression analysis shed light on this aspect of concerted development. Triple tt16 shp1 shp2 mutants produce malformed seedlings, seed coat formation defects, fewer seeds, and mucilage reduction. While shp1 shp2 mutants fail to coordinate the timely development of ovules, tt16 mutants show less peripheral endosperm after fertilization. Failure in coordinated division of the innermost integument layer in early ovule stages leads to inner seed coat defects in tt16 and tt16 shp1 shp2 triple mutant seeds. An antagonistic action of ABS and SHP1/SHP2 is observed in inner seed coat layer formation. Expression analysis also indicates that ABS represses SHP1, SHP2, and FRUITFUL expression. Our work shows that the evolutionary conserved Bsister genes are required not only for endothelium but also for endosperm development and genetically interact with SHP1 and SHP2 in a partially antagonistic manner.

A transcriptomics investigation into pine reproductive organ development

The development of reproductive structures in gymnosperms is still poorly studied because of a lack of genomic information and useful genetic tools. The hermaphroditic reproductive structure derived from unisexual gymnosperms is an even less studied aspect of seed plant evolution. To extend our understanding of the molecular mechanism of hermaphroditism and the determination of sexual identity of conifer reproductive structures in general, unisexual and bisexual cones from Pinus tabuliformis were profiled for gene expression using 60K microarrays. Expression patterns of genes during progression of sexual cone development were analysed using RNA-seq. The results showed that, overall, the transcriptomes of male structures in bisexual cones were more similar to those of female cones. However, the expression of several MADS-box genes in the bisexual cones was similar to that of male cones at the more juvenile developmental stage, while despite these expression shifts, male structures of bisexual cones and normal male cones were histologically indistinguishable and cone development was continuous. This study represents a starting point for in-depth analysis of the molecular regulation of cone development and also the origin of hermaphroditism in pine.

Transcriptome-wide analysis of the MADS-box gene family in the orchid Erycina pusilla

Orchids exhibit a range of unique flower shapes and are a valuable ornamental crop. MADS-box transcription factors are key regulatory components in flower initiation and development. Changing the flower shape and flowering time can increase the value of the orchid in the ornamental horticulture industry. In this study, 28 MADS-box genes were identified from the transcriptome database of the model orchid Erycina pusilla. The full-length genomic sequences of these MADS-box genes were obtained from BAC clones. Of these, 27 were MIKC-type EpMADS (two truncated forms) and one was a type I EpMADS. Eleven EpMADS genes contained introns longer than 10 kb. Phylogenetic analysis classified the 24 MIKC(c) genes into nine subfamilies. Three specific protein motifs, AG, FUL and SVP, were identified and used to classify three subfamilies. The expression profile of each EpMADS gene correlated with its putative function. The phylogenetic analysis was highly correlated with the protein domain identification and gene expression results. Spatial expression of EpMADS6, EpMADS12 and EpMADS15 was strongly detected in the inflorescence meristem, floral bud and seed via in situ hybridization. The subcellular localization of the 28 EpMADS proteins was also investigated. Although EpMADS27 lacks a complete MADS-box domain, EpMADS27-YFP was localized in the nucleus. This characterization of the orchid MADS-box family genes provides useful information for both orchid breeding and studies of flowering and evolution.

Strategies to improve low copy transgenic events in Agrobacterium-mediated transformation of maize

Transgenic plants containing low copy transgene insertion free of vector backbone are highly desired for many biotechnological applications. We have investigated two different strategies for increasing the percentage of low copy events in Agrobacterium-mediated transformation experiments in maize. One of the strategies is to use a binary vector with two separate T-DNAs, one T-DNA containing an intact E.coli manA gene encoding phosphomannose isomerase (PMI) as selectable marker gene cassette and another T-DNA containing an RNAi cassette of PMI sequences. By using this strategy, low copy transgenic events containing the transgenes were increased from 43 to 60 % in maize. An alternate strategy is using selectable marker gene cassettes containing regulatory or coding sequences derived from essential plant genes such as 5-enolpyruvylshikimate-3-phosphate synthase (EPSPS) or MADS box transcription factor. In this paper we demonstrate that higher percentage of low copy transgenic events can be obtained in Agrobacterium-mediated maize transformation experiments using both strategies. We propose that the above two strategies can be used independently or in combination to increase transgenic events that contain low copy transgene insertion in Agrobacterium-mediated transformation experiments.

The significance and scope of evolutionary developmental biology: a vision for the 21st century

Evolutionary developmental biology (evo-devo) has undergone dramatic transformations since its emergence as a distinct discipline. This paper aims to highlight the scope, power, and future promise of evo-devo to transform and unify diverse aspects of biology. We articulate key questions at the core of eleven biological disciplines-from Evolution, Development, Paleontology, and Neurobiology to Cellular and Molecular Biology, Quantitative Genetics, Human Diseases, Ecology, Agriculture and Science Education, and lastly, Evolutionary Developmental Biology itself-and discuss why evo-devo is uniquely situated to substantially improve our ability to find meaningful answers to these fundamental questions. We posit that the tools, concepts, and ways of thinking developed by evo-devo have profound potential to advance, integrate, and unify biological sciences as well as inform policy decisions and illuminate science education. We look to the next generation of evolutionary developmental biologists to help shape this process as we confront the scientific challenges of the 21st century.

The identification of novel PMADS3 interacting proteins indicates a role in post-transcriptional control

PMADS3, a known MADS-box transcriptional factor and a C-class gene for floral development, plays dual roles in controlling the identity of inner floral organs and the termination of flower meristems in petunia. In this study, it was confirmed by bimolecular fluorescence complementation (BiFC) assays that the PMADS3 protein can interact individually with E-class proteins FBP2, FBP5, FBP9 and PMADS12. A yeast two-hybrid cDNA library was screened using the entire PMADS3 as bait, and this identified further potential interaction candidates. Two novel genes, PheIF3f and PhAGO10, were isolated, and suggested to regulate mRNA and translational processes according to the analysis of protein functional domains and subcellular localization predictions. Notably, the PhAGO10 protein belongs to the Argonaute family, members of which are major players in small-RNA-guided gene silencing processes via mRNA cleavage or translational inhibition. The results of yeast two-hybrid and BiFC assays indicated that PheIF3f and PhAGO10 could interact with PMADS3. Our findings indicate that the C-class gene PMADS3 potentially participates in post-transcriptional control, as well as transcriptional regulation.

My favourite flowering image: a cob of pod corn

For good reasons scientists usually do not report the personal circumstances of their work when publishing their results. This means, however, that the scientific facts being reported may not accurately reflect the personal importance of the respective work for the individual scientists. Pictures of pod corn (or Tunicate maize) have been on my mind for much of my life, through good and through bad times. This is why...

MADS goes genomic in conifers: towards determining the ancestral set of MADS-box genes in seed plants

BACKGROUND AND AIMS

MADS-box genes comprise a gene family coding for transcription factors. This gene family expanded greatly during land plant evolution such that the number of MADS-box genes ranges from one or two in green algae to around 100 in angiosperms. Given the crucial functions of MADS-box genes for nearly all aspects of plant development, the expansion of this gene family probably contributed to the increasing complexity of plants. However, the expansion of MADS-box genes during one important step of land plant evolution, namely the origin of seed plants, remains poorly understood due to the previous lack of whole-genome data for gymnosperms.

METHODS

The newly available genome sequences of Picea abies, Picea glauca and Pinus taeda were used to identify the complete set of MADS-box genes in these conifers. In addition, MADS-box genes were identified in the growing number of transcriptomes available for gymnosperms. With these datasets, phylogenies were constructed to determine the ancestral set of MADS-box genes of seed plants and to infer the ancestral functions of these genes.

KEY RESULTS

Type I MADS-box genes are under-represented in gymnosperms and only a minimum of two Type I MADS-box genes have been present in the most recent common ancestor (MRCA) of seed plants. In contrast, a large number of Type II MADS-box genes were found in gymnosperms. The MRCA of extant seed plants probably possessed at least 11-14 Type II MADS-box genes. In gymnosperms two duplications of Type II MADS-box genes were found, such that the MRCA of extant gymnosperms had at least 14-16 Type II MADS-box genes.

CONCLUSIONS

The implied ancestral set of MADS-box genes for seed plants shows simplicity for Type I MADS-box genes and remarkable complexity for Type II MADS-box genes in terms of phylogeny and putative functions. The analysis of transcriptome data reveals that gymnosperm MADS-box genes are expressed in a great variety of tissues, indicating diverse roles of MADS-box genes for the development of gymnosperms. This study is the first that provides a comprehensive overview of MADS-box genes in conifers and thus will provide a framework for future work on MADS-box genes in seed plants.

Structural basis for the oligomerization of the MADS domain transcription factor SEPALLATA3 in Arabidopsis

In plants, MADS domain transcription factors act as central regulators of diverse developmental pathways. In Arabidopsis thaliana, one of the most central members of this family is SEPALLATA3 (SEP3), which is involved in many aspects of plant reproduction, including floral meristem and floral organ development. SEP3 has been shown to form homo and heterooligomeric complexes with other MADS domain transcription factors through its intervening (I) and keratin-like (K) domains. SEP3 function depends on its ability to form specific protein-protein complexes; however, the atomic level determinants of oligomerization are poorly understood. Here, we report the 2.5-Å crystal structure of a small portion of the intervening and the complete keratin-like domain of SEP3. The domains form two amphipathic alpha helices separated by a rigid kink, which prevents intramolecular association and presents separate dimerization and tetramerization interfaces comprising predominantly hydrophobic patches. Mutations to the tetramerization interface demonstrate the importance of highly conserved hydrophobic residues for tetramer stability. Atomic force microscopy was used to show SEP3-DNA interactions and the role of oligomerization in DNA binding and conformation. Based on these data, the oligomerization patterns of the larger family of MADS domain transcription factors can be predicted and manipulated based on the primary sequence.

Transcriptome-wide profiling and expression analysis of transcription factor families in a liverwort, Marchantia polymorpha

BACKGROUND

Transcription factors (TFs) are vital elements that regulate transcription and the spatio-temporal expression of genes, thereby ensuring the accurate development and functioning of an organism. The identification of TF-encoding genes in a liverwort, Marchantia polymorpha, offers insights into TF organization in the members of the most basal lineages of land plants (embryophytes). Therefore, a comparison of Marchantia TF genes with other land plants (monocots, dicots, bryophytes) and algae (chlorophytes, rhodophytes) provides the most comprehensive view of the rates of expansion or contraction of TF genes in plant evolution.

RESULTS

In this study, we report the identification of TF-encoding transcripts in M. polymorpha for the first time, as evidenced by deep RNA sequencing data. In total, 3,471 putative TF encoding transcripts, distributed in 80 families, were identified, representing 7.4% of the generated Marchantia gametophytic transcriptome dataset. Overall, TF basic functions and distribution across families appear to be conserved when compared to other plant species. However, it is of interest to observe the genesis of novel sequences in 24 TF families and the apparent termination of 2 TF families with the emergence of Marchantia. Out of 24 TF families, 6 are known to be associated with plant reproductive development processes. We also examined the expression pattern of these TF-encoding transcripts in six male and female developmental stages in vegetative and reproductive gametophytic tissues of Marchantia.

CONCLUSIONS

The analysis highlighted the importance of Marchantia, a model plant system, in an evolutionary context. The dataset generated here provides a scientific resource for TF gene discovery and other comparative evolutionary studies of land plants.

Molecular control of normal and acrocona mutant seed cone development in Norway spruce (Picea abies) and the evolution of conifer ovule-bearing organs

Reproductive organs in seed plants are morphologically divergent and their evolutionary history is often unclear. The mechanisms controlling their development have been extensively studied in angiosperms but are poorly understood in conifers and other gymnosperms. Here, we address the molecular control of seed cone development in Norway spruce, Picea abies. We present expression analyses of five novel MADS-box genes in comparison with previously identified MADS and LEAFY genes at distinct developmental stages. In addition, we have characterized the homeotic transformation from vegetative shoot to female cone and associated changes in regulatory gene expression patterns occurring in the acrocona mutant. The analyses identified genes active at the onset of ovuliferous and ovule development and identified expression patterns marking distinct domains of the ovuliferous scale. The reproductive transformation in acrocona involves the activation of all tested genes normally active in early cone development, except for an AGAMOUS-LIKE6/SEPALLATA (AGL6/SEP) homologue. This absence may be functionally associated with the nondeterminate development of the acrocona ovule-bearing scales. Our morphological and gene expression analyses give support to the hypothesis that the modern cone is a complex structure, and the ovuliferous scale the result of reductions and compactions of an ovule-bearing axillary short shoot in cones of Paleozoic conifers.

Crystallization studies of the keratin-like domain from Arabidopsis thaliana SEPALLATA 3

In higher plants, the MADS-box genes encode a large family of transcription factors (TFs) involved in key developmental processes, most notably plant reproduction, flowering and floral organ development. SEPALLATA 3 (SEP3) is a member of the MADS TF family and plays a role in the development of the floral organs through the formation of multiprotein complexes with other MADS-family TFs. SEP3 is divided into four domains: the M (MADS) domain, involved in DNA binding and dimerization, the I (intervening) domain, a short domain involved in dimerization, the K (keratin-like) domain important for multimeric MADS complex formation and the C (C-terminal) domain, a largely unstructured region putatively important for higher-order complex formation. The entire K domain along with a portion of the I and C domains of SEP3 was crystallized using high-throughput robotic screening followed by optimization. The crystals belonged to space group P2(1)2(1)2, with unit-cell parameters a = 123.44, b = 143.07, c = 49.83 Å, and a complete data set was collected to 2.53 Å resolution.

Molecular evolution and patterns of duplication in the SEP/AGL6-like lineage of the Zingiberales: a proposed mechanism for floral diversification

The diversity of floral forms in the plant order Zingiberales has evolved through alterations in floral organ morphology. One striking alteration is the shift from fertile, filamentous stamens to sterile, laminar (petaloid) organs in the stamen whorls, attributed to specific pollination syndromes. Here, we examine the role of the SEPALLATA (SEP) genes, known to be important in regulatory networks underlying floral development and organ identity, in the evolution of development of the diverse floral organs phenotypes in the Zingiberales. Phylogenetic analyses show that the SEP-like genes have undergone several duplication events giving rise to multiple copies. Selection tests on the SEP-like genes indicate that the two copies of SEP3 have mostly evolved under balancing selection, probably due to strong functional restrictions as a result of their critical role in floral organ specification. In contrast, the two LOFSEP copies have undergone differential positive selection, indicating neofunctionalization. Reverse transcriptase-polymerase chain reaction, gene expression from RNA-seq data, and in situ hybridization analyses show that the recovered genes have differential expression patterns across the various whorls and organ types found in the Zingiberales. Our data also suggest that AGL6, sister to the SEP-like genes, may play an important role in stamen morphology in the Zingiberales. Thus, the SEP-like genes are likely to be involved in some of the unique morphogenetic patterns of floral organ development found among this diverse order of tropical monocots. This work contributes to a growing body of knowledge focused on understanding the role of gene duplications and the evolution of entire gene networks in the evolution of flower development.

Genome-wide survey and expression analysis of the MADS-box gene family in soybean

MADS-box genes encode important transcription factors in plants that are involved in many processes during plant growth and development. An investigation of the soybean genome revealed 106 putative MADS-box genes. These genes were classified into two classes, type I and type II, based on phylogenetic analysis. The soybean type II group has 72 members, which is higher than that of Arabidopsis, indicating that soybean type II genes have undergone a higher rate of duplication and/or a lower rate of gene loss after duplication. Soybean MADS-box genes are present on all chromosomes. Like Arabidopsis and rice MADS-box genes, soybean MADS-box genes expanded through tandem gene duplication and segmental duplication events. There are many duplicate genes distributed across the soybean genome, with two genomic regions, i.e., MADS-box gene hotspots, where MADS-box genes with high degrees of similarity are clustered. Analysis of high-throughput sequencing data from soybean at different developmental stages and in different tissues revealed that MADS-box genes are expressed in embryos of various stages and in floral buds. This expression pattern suggests that soybean MADS-box genes play an important role in soybean growth and floral development.

An EST dataset for Metasequoia glyptostroboides buds: the first EST resource for molecular genomics studies in Metasequoia

The "living fossil" Metasequoia glyptostroboides Hu et Cheng, commonly known as dawn redwood or Chinese redwood, is the only living species in the genus and is valued for its essential oil and crude extracts that have great potential for anti-fungal activity. Despite its paleontological significance and economical value as a rare relict species, genomic resources of Metasequoia are very limited. In order to gain insight into the molecular mechanisms behind the formation of reproductive buds and the transition from vegetative phase to reproductive phase in Metasequoia, we performed sequencing of expressed sequence tags from Metasequoia vegetative buds and female buds. By using the 454 pyrosequencing technology, a total of 1,571,764 high-quality reads were generated, among which 733,128 were from vegetative buds and 775,636 were from female buds. These EST reads were clustered and assembled into 114,124 putative unique transcripts (PUTs) with an average length of 536 bp. The 97,565 PUTs that were at least 100 bp in length were functionally annotated by a similarity search against public databases and assigned with Gene Ontology (GO) terms. A total of 59 known floral gene families and 190 isotigs involved in hormone regulation were captured in the dataset. Furthermore, a set of PUTs differentially expressed in vegetative and reproductive buds, as well as SSR motifs and high confidence SNPs, were identified. This is the first large-scale expressed sequence tags ever generated in Metasequoia and the first evidence for floral genes in this critically endangered deciduous conifer species.

Live and let die - the B(sister) MADS-box gene OsMADS29 controls the degeneration of cells in maternal tissues during seed development of rice (Oryza sativa)

B_sister genes have been identified as the closest relatives of class B floral homeotic genes. Previous studies have shown that B_sister genes from eudicots are involved in cell differentiation during ovule and seed development. However, the complete function of B_sister genes in eudicots is masked by redundancy with other genes and little is known about the function of B_sister genes in monocots, and about the evolution of B_sister gene functions. Here we characterize OsMADS29, one of three MADS-box B_sister genes in rice. Our analyses show that OsMADS29 is expressed in female reproductive organs including the ovule, ovule vasculature, and the whole seed except for the outer layer cells of the pericarp. Knock-down of OsMADS29 by double-stranded RNA-mediated interference (RNAi) results in shriveled and/or aborted seeds. Histological analyses of the abnormal seeds at 7 days after pollination (DAP) indicate that the symplastic continuity, including the ovular vascular trace and the nucellar projection, which is the nutrient source for the filial tissue at early development stages, is affected. Moreover, degeneration of all the maternal tissues in the transgenic seeds, including the pericarp, ovular vascular trace, integuments, nucellar epidermis and nucellar projection, is blocked as compared to control plants. Our results suggest that OsMADS29 has important functions in seed development of rice by regulating cell degeneration of maternal tissues. Our findings provide important insights into the ancestral function of B_sister genes.

Selaginella Genome Analysis - Entering the "Homoplasy Heaven" of the MADS World

In flowering plants, arguably the most significant transcription factors regulating development are MADS-domain proteins, encoded by Type I and Type II MADS-box genes. Type II genes are divided into the MIKC(C) and MIKC* groups. In angiosperms, these types and groups play distinct roles in the development of female gametophytes, embryos, and seeds (Type I); vegetative and floral tissues in sporophytes (MIKC(C)); and male gametophytes (MIKC*), but their functions in other plants are largely unknown. The complete set of MADS-box genes has been described for several angiosperms and a moss, Physcomitrella patens. Our examination of the complete genome sequence of a lycophyte, Selaginella moellendorffii, revealed 19 putative MADS-box genes (13 Type I, 3 MIKC(C), and 3 MIKC*). Our results suggest that the most recent common ancestor of vascular plants possessed at least two Type I and two Type II genes. None of the S. moellendorffii MIKC(C) genes were identified as orthologs of any floral organ identity genes. This strongly corroborates the view that the clades of floral organ identity genes originated in a common ancestor of seed plants after the lineage that led to lycophytes had branched off, and that expansion of MIKC(C) genes in the lineage leading to seed plants facilitated the evolution of their unique reproductive organs. The number of MIKC* genes and the ratio of MIKC* to MIKC(C) genes is lower in S. moellendorffii and angiosperms than in P. patens, correlated with reduction of the gametophyte in vascular plants. Our data indicate that Type I genes duplicated and diversified independently within lycophytes and seed plants. Our observations on MADS-box gene evolution echo morphological evolution since the two lineages of vascular plants appear to have arrived independently at similar body plans. Our annotation of MADS-box genes in S. moellendorffii provides the basis for functional studies to reveal the roles of this crucial gene family in basal vascular plants.

Isolation and characterization of a SEPALLATA-like gene, ZjMADS1, from marine angiosperm Zostera japonica

In flowering plants, floral homeotic MADS-box genes, which constitute a large multigene family, play important roles in the specification of floral organs as defined by the ABCDE model. In this study, a MADS-box gene, ZjMADS1, was isolated and characterized from the marine angiosperm Zostera japonica. The predicted length of the ZjMADS1 protein was 246 amino acids (AA), and the AA sequence was most similar to those of the SEPALLATA (SEP) subfamily, corresponding to E-function genes. Southern blot analysis suggested the presence of two SEP3-like genes in the Z. japonica genome. ZjMADS1 mRNA levels were extremely high in the spadices, regardless of the developmental stage, compared to other organs from the reproductive and vegetative shoots. These results suggest that the ZjMADS1 gene may be involved in spadix development in Z. japonica and act as an E-function gene in floral organ development in marine angiosperms.

Evolutionary and comparative analysis of MYB and bHLH plant transcription factors

The expansion of gene families encoding regulatory proteins is typically associated with the increase in complexity characteristic of multi-cellular organisms. The MYB and basic helix-loop-helix (bHLH) families provide excellent examples of how gene duplication and divergence within particular groups of transcription factors are associated with, if not driven by, the morphological and metabolic diversity that characterize the higher plants. These gene families expanded dramatically in higher plants; for example, there are approximately 339 and 162 MYB and bHLH genes, respectively, in Arabidopsis, and approximately 230 and 111, respectively, in rice. In contrast, the Chlamydomonas genome has only 38 MYB genes and eight bHLH genes. In this review, we compare the MYB and bHLH gene families from structural, evolutionary and functional perspectives. The knowledge acquired on the role of many of these factors in Arabidopsis provides an excellent reference to explore sequence-function relationships in crops and other plants. The physical interaction and regulatory synergy between particular sub-classes of MYB and bHLH factors is perhaps one of the best examples of combinatorial plant gene regulation. However, members of the MYB and bHLH families also interact with a number of other regulatory proteins, forming complexes that either activate or repress the expression of sets of target genes that are increasingly being identified through a diversity of high-throughput genomic approaches. The next few years are likely to witness an increasing understanding of the extent to which conserved transcription factors participate at similar positions in gene regulatory networks across plant species.

Molecular interactions of orthologues of floral homeotic proteins from the gymnosperm Gnetum gnemon provide a clue to the evolutionary origin of 'floral quartets'

Several lines of evidence suggest that the identity of floral organs in angiosperms is specified by multimeric transcription factor complexes composed of MADS-domain proteins. These bind to specific cis-regulatory elements ('CArG-boxes') of their target genes involving DNA-loop formation, thus constituting 'floral quartets'. Gymnosperms, angiosperms' closest relatives, contain orthologues of floral homeotic genes, but when and how the interactions constituting floral quartets were established during evolution has remained unknown. We have comprehensively studied the dimerization and DNA-binding of several classes of MADS-domain proteins from the gymnosperm Gnetum gnemon. Determination of protein-protein and protein-DNA interactions by yeast two-hybrid, in vitro pull-down and electrophoretic mobility shift assays revealed complex patterns of homo- and heterodimerization among orthologues of floral homeotic class B, class C and class E proteins and B(sister) proteins. Using DNase I footprint assays we demonstrate that both orthologues of class B with C proteins, and orthologues of class C proteins alone, but not orthologues of class B proteins alone can loop DNA in floral quartet-like complexes. This is in contrast to class B and class C proteins from angiosperms, which require other factors such as class E floral homeotic proteins to 'glue' them together in multimeric complexes. Our findings suggest that the evolutionary origin of floral quartet formation is based on the interaction of different DNA-bound homodimers, does not depend on class E proteins, and predates the origin of angiosperms.

The Arabidopsis B-sister MADS-box protein, GORDITA, represses fruit growth and contributes to integument development

The MADS-box family of transcription factors have diverse developmental roles in flower pattern formation, gametophyte cell division and fruit differentiation. The B-sister MADS-box proteins are most similar to the B-class floral homeotic proteins, and are expressed in female reproductive organs. The Arabidopsis B-sister MADS-box protein, TT16, is necessary for inner integument differentiation. We have functionally characterized the only other B-sister MADS-box gene in Arabidopsis, AGL63, renamed here as GORDITA (GOA). A loss-of-function mutation in goa or reduction of endogenous GOA expression results in larger fruits, illustrating its novel function in regulating fruit growth. Consistent with its function, GOA expression is detected in the walls of the valves and throughout the replum of the fruit. Our phenotypic and molecular analyses of 35S::GOA and goa plants show that GOA controls organ size via cell expansion. Further, functional studies of goa tt16 double mutants have shown their additive role in controlling seed coat development, and have revealed the importance of GOA expression in the outer integument. Together, our studies provide evidence of a new regulatory role for a B-sister MADS-box gene in the control of organ growth.

Functional conservation and diversification of class E floral homeotic genes in rice (Oryza sativa)

Mutant analyses in different eudicotyledonous flowering plants demonstrated that SEPALLATA-like MADS-box genes are required for the specification of sepals, petals, stamens and carpels, and for floral determinacy, thus defining class E floral organ identity genes. SEP-like genes encode MADS-domain transcription factors and constitute an angiosperm-specific gene clade whose members show remarkably different degrees of redundancy and sub-functionalization within eudicots. To better understand the evolutionary dynamics of SEP-like genes throughout the angiosperms we have knocked down SEP-like genes of rice (Oryza sativa), a distant relative of eudicots within the flowering plants. Plants affected in both OsMADS7 and OsMADS8 show severe phenotypes including late flowering, homeotic changes of lodicules, stamens and carpels into palea/lemma-like organs, and a loss of floral determinacy. Simultaneous knockdown of the four rice SEP-like genes OsMADS1, OsMADS5, OsMADS7 and OsMADS8, leads to homeotic transformation of all floral organs except the lemma into leaf-like organs. This mimics the phenotype observed with the sep1 sep2 sep3 sep4 quadruple mutant of Arabidopsis. Detailed analyses of the spatial and temporal mRNA expression and protein interaction patterns corresponding to the different rice SEP-like genes show strong similarities, but also gene-specific differences. These findings reveal conservation of SEP-like genes in specifying floral determinacy and organ identities since the separation of eudicots and monocots about 150 million years ago. However, they indicate also monocot-specific neo- and sub-functionalization events and hence underscore the evolutionary dynamics of SEP-like genes. Moreover, our findings corroborate the view that the lodicules of grasses are homologous to eudicot petals.