A Survey of MIKC Type MADS-Box Genes in Non-seed Plants: Algae, Bryophytes, Lycophytes and Ferns

OsMADS6-OsMADS32 and REP1 control palea cellular heterogeneity and morphogenesis in rice

Precise regulation of cell proliferation and differentiation is vital for organ morphology. Rice palea, serving as sepal, comprises two distinct regions: the marginal region (MRP) and body of palea (BOP), housing heterogeneous cell populations, which makes it an ideal system for studying organ morphogenesis. We report that the transcription factor (TF) REP1 promotes epidermal cell proliferation and differentiation in the BOP, resulting in hard silicified protrusion cells, by regulating the cyclin-dependent kinase gene, OsCDKB1;1. Conversely, TFs OsMADS6 and OsMADS32 are expressed exclusively in the MRP, where they limit cell division rates by inhibiting OsCDKB2;1 expression and promote endoreduplication, yielding elongated epidermal cells. Furthermore, reciprocal inhibition between the OsMADS6-OsMADS32 complex and REP1 fine-tunes the balance between cell division and differentiation during palea morphogenesis. We further show the functional conservation of these organ identity genes in heterogeneous cell growth in Arabidopsis, emphasizing a critical framework for controlling cellular heterogeneity in organ morphogenesis.

The Shoot Apical Meristem: An Evolutionary Molding of Higher Plants

The shoot apical meristem (SAM) gives rise to the aerial structure of plants by producing lateral organs and other meristems. The SAM is responsible for plant developmental patterns, thus determining plant morphology and, consequently, many agronomic traits such as the number and size of fruits and flowers and kernel yield. Our current understanding of SAM morphology and regulation is based on studies conducted mainly on some angiosperms, including economically important crops such as maize (Zea mays) and rice (Oryza sativa), and the model species Arabidopsis (Arabidopsis thaliana). However, studies in other plant species from the gymnosperms are scant, making difficult comparative analyses that help us understand SAM regulation in diverse plant species. This limitation prevents deciphering the mechanisms by which evolution gave rise to the multiple plant structures within the plant kingdom and determines the conserved mechanisms involved in SAM maintenance and operation. This review aims to integrate and analyze the current knowledge of SAM evolution by combining the morphological and molecular information recently reported from the plant kingdom.

Nuclear translocation of OsMADS25 facilitated by OsNAR2.1 in reponse to nitrate signals promotes rice root growth by targeting OsMADS27 and OsARF7

Nitrate is an important nitrogen source and signaling molecule that regulates plant growth and development. Although several components of the nitrate signaling pathway have been identified, the detailed mechanisms are still unclear. Our previous results showed that OsMADS25 can regulate root development in response to nitrate signals, but the mechanism is still unknown. Here, we try to answer two key questions: how does OsMADS25 move from the cytoplasm to the nucleus, and what are the direct target genes activated by OsMADS25 to regulate root growth after it moves to the nucleus in response to nitrate? Our results demonstrated that OsMADS25 moves from the cytoplasm to the nucleus in the presence of nitrate in an OsNAR2.1-dependent manner. Chromatin immunoprecipitation sequencing, chromatin immunoprecipitation qPCR, yeast one-hybrid, and luciferase experiments showed that OsMADS25 directly activates the expression of OsMADS27 and OsARF7, which are reported to be associated with root growth. Finally, OsMADS25-RNAi lines, the Osnar2.1 mutant, and OsMADS25-RNAi Osnar2.1 lines exhibited significantly reduced root growth compared with the wild type in response to nitrate supply, and expression of OsMADS27 and OsARF7 was significantly suppressed in these lines. Collectively, these results reveal a new mechanism by which OsMADS25 interacts with OsNAR2.1. This interaction is required for nuclear accumulation of OsMADS25, which promotes OsMADS27 and OsARF7 expression and root growth in a nitrate-dependent manner.

Structural evidence for MADS-box type I family expansion seen in new assemblies of Arabidopsis arenosa and A. lyrata

Arabidopsis thaliana diverged from A. arenosa and A. lyrata at least 6 million years ago. The three species differ by genome-wide polymorphisms and morphological traits. The species are to a high degree reproductively isolated, but hybridization barriers are incomplete. A special type of hybridization barrier is based on the triploid endosperm of the seed, where embryo lethality is caused by endosperm failure to support the developing embryo. The MADS-box type I family of transcription factors is specifically expressed in the endosperm and has been proposed to play a role in endosperm-based hybridization barriers. The gene family is well known for its high evolutionary duplication rate, as well as being regulated by genomic imprinting. Here we address MADS-box type I gene family evolution and the role of type I genes in the context of hybridization. Using two de-novo assembled and annotated chromosome-level genomes of A. arenosa and A. lyrata ssp. petraea we analyzed the MADS-box type I gene family in Arabidopsis to predict orthologs, copy number, and structural genomic variation related to the type I loci. Our findings were compared to gene expression profiles sampled before and after the transition to endosperm cellularization in order to investigate the involvement of MADS-box type I loci in endosperm-based hybridization barriers. We observed substantial differences in type-I expression in the endosperm of A. arenosa and A. lyrata ssp. petraea, suggesting a genetic cause for the endosperm-based hybridization barrier between A. arenosa and A. lyrata ssp. petraea.

Updated Phylogeny and Protein Structure Predictions Revise the Hypothesis on the Origin of MADS-box Transcription Factors in Land Plants

MADS-box transcription factors (TFs), among the first TFs extensively studied, exhibit a wide distribution across eukaryotes and play diverse functional roles. Varying by domain architecture, MADS-box TFs in land plants are categorized into Type I (M-type) and Type II (MIKC-type). Type I and II genes have been considered orthologous to the SRF and MEF2 genes in animals, respectively, presumably originating from a duplication before the divergence of eukaryotes. Here, we exploited the increasing availability of eukaryotic MADS-box sequences and reassessed their evolution. While supporting the ancient duplication giving rise to SRF- and MEF2-types, we found that Type I and II genes originated from the MEF2-type genes through another duplication in the most recent common ancestor (MRCA) of land plants. Protein structures predicted by AlphaFold2 and OmegaFold support our phylogenetic analyses, with plant Type I and II TFs resembling the MEF2-type structure, rather than SRFs. We hypothesize that the ancestral SRF-type TFs were lost in the MRCA of Archaeplastida (the kingdom Plantae sensu lato). The retained MEF2-type TFs acquired a Keratin-like domain and became MIKC-type before the divergence of Streptophyta. Subsequently in the MRCA of land plants, M-type TFs evolved from a duplicated MIKC-type precursor through loss of the Keratin-like domain, leading to the Type I clade. Both Type I and II TFs expanded and functionally differentiated in concert with the increasing complexity of land plant body architecture. The recruitment of these originally stress-responsive TFs into developmental programs, including those underlying reproduction, may have facilitated the adaptation to the terrestrial environment.

Cracking the Floral Quartet Code: How Do Multimers of MIKCC-Type MADS-Domain Transcription Factors Recognize Their Target Genes?

MADS-domain transcription factors (MTFs) are involved in the control of many important processes in eukaryotes. They are defined by the presence of a unique and highly conserved DNA-binding domain, the MADS domain. MTFs bind to double-stranded DNA as dimers and recognize specific sequences termed CArG boxes (such as 5'-CC(A/T)₆GG-3') and similar sequences that occur hundreds of thousands of times in a typical flowering plant genome. The number of MTF-encoding genes increased by around two orders of magnitude during land plant evolution, resulting in roughly 100 genes in flowering plant genomes. This raises the question as to how dozens of different but highly similar MTFs accurately recognize the cis-regulatory elements of diverse target genes when the core binding sequence (CArG box) occurs at such a high frequency. Besides the usual processes, such as the base and shape readout of individual DNA sequences by dimers of MTFs, an important sublineage of MTFs in plants, termed MIKC^C-type MTFs (M^C-MTFs), has evolved an additional mechanism to increase the accurate recognition of target genes: the formation of heterotetramers of closely related proteins that bind to two CArG boxes on the same DNA strand involving DNA looping. M^C-MTFs control important developmental processes in flowering plants, ranging from root and shoot to flower, fruit and seed development. The way in which M^C-MTFs bind to DNA and select their target genes is hence not only of high biological interest, but also of great agronomic and economic importance. In this article, we review the interplay of the different mechanisms of target gene recognition, from the ordinary (base readout) via the extravagant (shape readout) to the idiosyncratic (recognition of the distance and orientation of two CArG boxes by heterotetramers of M^C-MTFs). A special focus of our review is on the structural prerequisites of M^C-MTFs that enable the specific recognition of target genes.

Life cycle and functional genomics of the unicellular red alga Galdieria for elucidating algal and plant evolution and industrial use

Sexual reproduction has not been observed in unicellular red algae and Glaucophyceae, early branching groups in Archaeplastida, in which red algae and Viridiplantae independently evolved multicellular sexual life cycles. The finding of sexual reproduction in the unicellular red alga Galdieria provides information on the missing link of life cycle evolution in Archaeplastida. In addition, the metabolic plasticity, the polyextremophilic features, a relatively small genome, transcriptome data for the diploid and haploid, and the genetic modification tools developed here provide a useful platform for understanding the evolution of Archaeplastida, photosynthesis, metabolism, and environmental adaptation. For biotechnological use of the information and tools of Galdieria, the newly found cell wall–less haploid makes cell disruption less energy/cost intensive than the cell-walled diploid.

Sexual reproduction is widespread in eukaryotes; however, only asexual reproduction has been observed in unicellular red algae, including Galdieria, which branched early in Archaeplastida. Galdieria possesses a small genome; it is polyextremophile, grows either photoautotrophically, mixotrophically, or heterotrophically, and is being developed as an industrial source of vitamins and pigments because of its high biomass productivity. Here, we show that Galdieria exhibits a sexual life cycle, alternating between cell-walled diploid and cell wall–less haploid, and that both phases can proliferate asexually. The haploid can move over surfaces and undergo self-diploidization or generate heterozygous diploids through mating. Further, we prepared the whole genome and a comparative transcriptome dataset between the diploid and haploid and developed genetic tools for the stable gene expression, gene disruption, and selectable marker recycling system using the cell wall–less haploid. The BELL/KNOX and MADS-box transcription factors, which function in haploid-to-diploid transition and development in plants, are specifically expressed in the haploid and diploid, respectively, and are involved in the haploid-to-diploid transition in Galdieria, providing information on the missing link of the sexual life cycle evolution in Archaeplastida. Four actin genes are differently involved in motility of the haploid and cytokinesis in the diploid, both of which are myosin independent and likely reflect ancestral roles of actin. We have also generated photosynthesis-deficient mutants, such as blue-colored cells, which were depleted in chlorophyll and carotenoids, for industrial pigment production. These features of Galdieria facilitate the understanding of the evolution of algae and plants and the industrial use of microalgae.

Damaged Dickinsonia specimens provide clues to Ediacaran vendobiont biology

Recently reported specimens of the enigmatic Ediacaran fossil Dickinsonia from Russia show damage and repair that provides evidence of how they grew, and of their biological affinities. Marginal and terminal areas of wilting deformation are necrotic zones separating regenerated growth, sometimes on two divergent axes, rather than a single axis. Necrotic zones of damage to Dickinsonia are not a thick scar or callus, like a wound or amputation. Nor are they smooth transitions to a regenerated tail or arm. The wilted necrotic zone is most like damage by freezing, salt, or sunburn of leaves and lichens, compatible with evidence of terrestrial habitat from associated frigid and gypsic paleosols. Dickinsonia did not regrow by postembryonic addition of modules from a subterminal or patterned growth zone as in earthworms, myriapods, trilobites, crustaceans, and lizards. Rather Dickinsonia postembryonic regrowth from sublethal damage was from microscopic apical and lateral meristems, as in plants and lichens. Considered as fungal, Dickinsonia, and perhaps others of Class Vendobionta, were more likely Glomeromycota or Mucoromycotina, rather than Ascomycota or Basidiomycota.

AM22, a novel synthetic microRNA, inhibits the proliferation of colorectal cancer cells by targeting core binding factor subunit β (CBFB)

Our previous studies have revealed the important roles of the nonseed regions of microRNAs (miRNAs) in gene regulation, which provided novel insight into the development of miRNA analogs for cancer therapy. Here, we altered each nucleotide in the nonseed region of miR-34a and obtained novel synthetic miRNA analogs. Among them, AM22, with a base alteration from G to C at the 17^th nucleotide of miR-34a, showed extensive antiproliferative activity against several colorectal tumor cell lines and achieved effective inhibition of core binding factor subunit β (CBFB) expression. Subsequent investigations demonstrated that AM22 directly targeted CBFB by binding to its 3'-untranslated region (3'-UTR). Inhibition of CBFB showed obvious antiproliferative activity on HCT-116 and SW620 cells. Furthermore, the antiproliferative effects of AM22 on these cells were also measured in xenograft mouse models. In conclusion, this study identified AM22 as a potential antitumor miRNA by targeting CBFB and provided a new design approach for miRNA-based cancer treatment by changing the nonseed region of miRNA.

Ectopic expression of Torenia fournieri TCP8 and TCP13 alters the leaf and petal phenotypes in Arabidopsis thaliana

Teosinte branched1/cycloidea/proliferating cell factor (TCP) transcription factors (TFs) are essential for regulating plant developmental processes, which is still largely unknown in Torenia fournieri (T. fournieri), a widely used horticultural flower. In this study, we used a de novo transcriptome assembly method to predict the TCP transcription factors in T. fournieri. In total, 15 out of 21 predicted T. fournieri TCPs (TfTCPs) were isolated and verified with Sanger sequencing. Phylogenetic analysis showed that these 15 TfTCPs could be classified into two major classes. Most of these TfTCPs were expressed in floral buds, flowers, or leaves, suggesting an important role in developmental regulation in these tissues. Moreover, TfTCP8 and TfTCP13, the homologues of the Arabidopsis thaliana TCP5-like transcription factor, were able to bind to the conserved Class II TCP binding motifs and are localized to the nucleus, indicating that TfTCP8 and TfTCP13 act as transcriptional regulators. In agreement with the overexpression phenotype of AtTCP5, ectopic expression of TfTCP8 and TfTCP13 resulted in narrow leaves and the small petal phenotype in Arabidopsis, suggesting that these two TfTCPs potentially regulate leaf or flower shape in T. fournieri.

CsubMADS1, a lag phase transcription factor, controls development of polar eukaryotic microalga Coccomyxa subellipsoidea C-169

MADS-box transcription factors (TFs) have not been functionally delineated in microalgae. In this study, the role of CsubMADS1 from microalga Coccomyxa subellipsoidea C-169 has been explored. Unlike Type II MADS-box proteins of seed plants with MADS, Intervening, K-box, and C domains, CsubMADS1 only has MADS and Intervening domains. It forms a group with MADS TFs from algae in the phylogenetic tree within the Type II MIKC^C clade. CsubMADS1 is expressed strongly in the lag phase of growth. The CsubMADS1 monomer does not have a specific localization in the nucleus, and it forms homodimers to localize exclusively in the nucleus. The monomer has two nuclear localization signals (NLSs): an N-terminal NLS and an internal NLS. The internal NLS is functional, and the homodimer requires two NLSs for specific nuclear localization. Overexpression (OX) of CsubMADS1 slows down the growth of the culture and leads to the creation of giant polyploid multinucleate cells, resembling autospore mother cells. This implies that the release of autospores from autospore mother cells may be delayed. Thus, in wild-type (WT) cells, CsubMADS1 may play a crucial role in slowing down growth during the lag phase. Due to starvation in 2-month-old colonies on solid media, the WT colonies produce mucilage, whereas OX colonies produce significantly less mucilage. Thus, CsubMADS1 also negatively regulates stress-induced mucilage production and probably plays a role in stress tolerance during the lag phase. Taken together, our results reveal that CsubMADS1 is a key TF involved in the development and stress tolerance of this polar microalga.

Comparative analysis of Phytophthora genomes reveals oomycete pathogenesis in crops

The oomycete genus Phytophthora includes devastating plant pathogens that are found in almost all ecosystems. We sequenced the genomes of two quarantined Phytophthora species–P. fragariae and P. rubi. Comparing these Phytophthora species and related genera allowed reconstruction of the phylogenetic relationships within the genus Phytophthora and revealed Phytophthora genomic features associated with infection and pathogenicity. We found that several hundred Phytophthora genes are putatively inherited from red algae, but Phytophthora does not have vestigial plastids originating from phototrophs. The horizontally-transferred Phytophthora genes are abundant transposons that “transmit” exogenous gene to Phytophthora species thus bring about the gene recombination possibility. Several expansion events of Phytophthora gene families associated with cell wall biogenesis can be used as mutational targets to elucidate gene function in pathogenic interactions with host plants. This work enhanced the understanding of Phytophthora evolution and will also be helpful for the design of phytopathological control strategies.

All type II classic MADS-box genes in the lycophyte Selaginella moellendorffii are broadly yet discretely expressed in vegetative and reproductive tissues

The MADS-box genes constitute a large transcription factor family that appear to have evolved by duplication and diversification of function. Two types of MADS-box genes are distinguished throughout eukaryotes, types I and II. Type II classic MADS-box genes, also known as MIKC-type, are key developmental regulators in flowering plants and are particularly well-studied for their role in floral organ specification. However, very little is known about the role that these genes might play outside of the flowering plants. We investigated the evolution of type II classic MADS-box genes across land plants by performing a maximum likelihood analysis with a particular focus on lycophytes. Here, we present the expression patterns of all three type II classic MADS-box homologs throughout plant development in the lycophyte Selaginella moellendorffii: SmMADS1, SmMADS3, and SmMADS6. We used scanning electron microscopy and histological analyses to define stages of sporangia development in S. moellendorffii. We performed phylogenetic analyses of this gene lineage across land plants and found that lycophyte sequences appeared before the multiple duplication events that gave rise to the major MADS-box gene lineages in seed plants. Our expression analyses by in situ hybridization show that all type II classic MADS-box genes in S. moellendorffii have broad but distinct patterns of expression in vegetative and reproductive tissues, where SmMADS1 and SmMADS6 only differ during late sporangia development. The broad expression during S. moellendorffii development suggests that MADS-box genes have undergone neofunctionalization and subfunctionalization after duplication events in seed plants.

Bibenzyls and bisbybenzyls of bryophytic origin as promising source of novel therapeutics: pharmacology, synthesis and structure-activity

BACKGROUND

The amphibian, non-vascular, gametophyte-dominant, bio-indicator class, bryophytes; with their wide ranges of habitat have attained importance due to their promising medicinal attributions and therapeutic role; mostly aided by presence of aromatic bibenzyl and bisbybenzyl class of compounds. Bibenzyls are steroidal ethane derivatives, resembling the structural moiety of bioactive dihydro-stilbenoids or iso-quinoline alkaloids. These stress triggered secondary metabolites are the by-products of the flavonoid biosynthetic pathway. Different classes of bryophytes (Bryophyta, Marchantiophyta and Anthocerotophyta) possess different subtypes of bibenzyls and dimeric bisbibenzyls. Among the liverwort, hornwort and mosses, former one is mostly enriched with bibenzyl type constituents as per the extensive study conducted for phytochemical deposit. Considering macrocyclic and acyclic group of bibenzyls and bisbybenzyls, generally marchantin type compounds are reported vividly for significant biological activity that includes neuro-nephro-cardio-protection besides anti-allergic, anti-microbial, anti-apoptotic and cytotoxic activities studied on in-vitro and in-vivo models or on cell lines.

RESULT

The critical analysis of reported chemical and pharmaceutical attributions of bibenzyls and bis-bibenzyls yielded detailed report on this compound class along with their application, mode of action, natural source, techniques of synthesis, extraction procedure, isolation and characterization. Further, the structure activity relationship studies and bioactivity of bibenzyls derived from non-bryophytic origin were also summarized.

CONCLUSION

This review encompasses prospective biological application of botanical reservoir of this primarily ignored, primeval land plant group where recent technical advances has paved the way for qualitative and quantitative isolation and estimation of novel compounds as well as marker components to study their impact on environment, as bio-control agents and as key leads in future drug designing. Graphical abstract.

Radial or Bilateral? The Molecular Basis of Floral Symmetry

In the plant kingdom, the flower is one of the most relevant evolutionary novelties. Floral symmetry has evolved multiple times from the ancestral condition of radial to bilateral symmetry. During evolution, several transcription factors have been recruited by the different developmental pathways in relation to the increase of plant complexity. The MYB proteins are among the most ancient plant transcription factor families and are implicated in different metabolic and developmental processes. In the model plant Antirrhinum majus, three MYB transcription factors (DIVARICATA, DRIF, and RADIALIS) have a pivotal function in the establishment of floral dorsoventral asymmetry. Here, we present an updated report of the role of the DIV, DRIF, and RAD transcription factors in both eudicots and monocots, pointing out their functional changes during plant evolution. In addition, we discuss the molecular models of the establishment of flower symmetry in different flowering plants.

Genomic and Transcriptional Analysis of Banana Ovate Family Proteins Reveals Their Relationship with Fruit Development and Ripening

Ovate Family Proteins (OFPs) belong to a plant-specific transcription factor family. They have been found to have significant roles in growth and development in Arabidopsis and tomato; however, little is known regarding their role in banana. Thus, a genome-wide study of OFP genes in banana was conducted for the first time in the present study. The results demonstrated that 49 OFP family members are unequally distributed across 11 chromosomes. Phylogenetic analysis grouped these genes into two subfamilies and eight subgroups, which was confirmed by the conserved motif and gene structure analysis. Furthermore, MaOFPs genes duplicates were found to have originated from whole-genome duplication (WGD). The expression patterns of the genes in the various tissues and at different fruit development and ripening stages in the BaXi Jiao (BX) and Feng Jiao (FJ), banana cultivars were elucidated using transcriptome analysis. Using co-expression network analysis, MaOFP1 was found to interact not only with MaMADS36 but also with hormone response proteins. These findings improve our understanding of the functions of MaOFPs genes in the control of plant hormone signal transduction pathways during banana growth and ripening, which should inform the genetic improvement of important agricultural characters.

Structural Basis for Plant MADS Transcription Factor Oligomerization

MADS transcription factors (TFs) are DNA binding proteins found in almost all eukaryotes that play essential roles in diverse biological processes. While present in animals and fungi as a small TF family, the family has dramatically expanded in plants over the course of evolution, with the model flowering plant, Arabidopsis thaliana, possessing over 100 type I and type II MADS TFs. All MADS TFs contain a core and highly conserved DNA binding domain called the MADS or M domain. Plant MADS TFs have diversified this domain with plant-specific auxiliary domains. Plant type I MADS TFs have a highly diverse and largely unstructured Carboxy-terminal (C domain), whereas type II MADS have added oligomerization domains, called Intervening (I domain) and Keratin-like (K domain), in addition to the C domain. In this mini review, we describe the overall structure of the type II "MIKC" type MADS TFs in plants, with a focus on the K domain, a critical oligomerization module. We summarize the determining factors for oligomerization and provide mechanistic insights on how secondary structural elements are required for oligomerization capability and specificity. Using MADS TFs that are involved in flower organ specification as an example, we provide case studies and homology modeling of MADS TFs complex formation. Finally, we highlight outstanding questions in the field.

Ectopic expression of the AaFUL1 gene identified in Anthurium andraeanum affected floral organ development and seed fertility in tobacco

Anthurium andraeanum is a high-grade potted flower that enjoys global popularity. Its floral organs have been substantially modified, and its ornamental value is based on its petaloid bracts. MADS-box gene products are important transcription factors that control plant development. In particular, the APETALA1 (AP1)/FRUITFULL (FUL) family of MADS-box genes plays a key role in flowering transitions and out-whorl floral organ identity specification. In this report, one FUL-like gene was cloned from Anthurium andraeanum and named AaFUL1 after bioinformatics identification. Subsequent subcellular localization experiments confirmed that the AaFUL1 protein was located in the nucleus, and data obtained from an expression analysis indicated that the relative expression level of AaFUL1 was the highest in bracts and inflorescences, while its expression was relatively low in stems and roots. Next, an AaFUL1 overexpression vector was constructed and ectopically expressed in tobacco. The transformants did not show any early flowering phenotype, but the average internode length of the inflorescence branch was significantly higher than that observed in the control, and its petal color had substantially faded. The morphology of the petal and pistil was clearly changed, the fruit was deformed, and the seed was largely aborted. These data indicate that even though the sequence of AaFUL1 is relatively conserved, its function differs from that of other orthologs, and the FUL subfamily of MADS-box transcription factors may have taken on new functions during the evolution processes. The results of this experiment enrich our knowledge of FUL transcription factors in monocotyledon 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.

Dating Whole Genome Duplication in Ceratopteris thalictroides and Potential Adaptive Values of Retained Gene Duplicates

Whole-genome duplications (WGDs) are widespread in plants and frequently coincide with global climatic change events, such as the Cretaceous–Tertiary (KT) extinction event approximately 65 million years ago (mya). Ferns have larger genomes and higher chromosome numbers than seed plants, which likely resulted from multiple rounds of polyploidy. Here, we use diploid and triploid material from a model fern species, Ceratopteris thalictroides, for the detection of WGDs. High-quality RNA-seq data was used to infer the number of synonymous substitutions per synonymous site (Ks) between paralogs; Ks age distribution and absolute dating approach were used to determine the age of WGD events. Evidence of an ancient WGD event with a Ks peak value of approximately 1.2 was obtained for both samples; however, the Ks frequency distributions varied significantly. Importantly, we dated the WGD event at 51–53 mya, which coincides with the Paleocene-Eocene Thermal Maximum (PETM), when the Earth became warmer and wetter than any other period during the Cenozoic. Duplicate genes were preferentially retained for specific functions, such as environment response, further support that the duplicates may have promoted quick adaption to environmental changes and potentially resulted in evolutionary success, especially for pantropical species, such as C. thalictroides, which exhibits higher temperature tolerance.

Asymmetric birth and death of type I and type II MADS-box gene subfamilies in the rubber tree facilitating laticifer development

The rubber tree (Hevea brasiliensis Muell. Arg.) is a rubber producing crop and contains specialized laticifers. MADS-box genes are a family of transcription factor genes that regulate plant development, especially floral organ and gametophyte development. 97 MADS-box genes were identified in the rubber tree through transcriptomes and genome mining. 93.8% of the genes were mapped onto the genome scaffolds in correspondence to the coverage (93.8%) of current version of sequenced genome. Phylogenetic analysis indicates that type II MADS-box genes have been more actively duplicated than their orthologous genes in Arabidopsis and rice, so that most (70, 72.2%) of the MADS-box genes in the rubber tree belong to type II subfamily. This is a high percentage compared to those in Arabidopsis (43.7%) and rice (56.8%). Moreover, 69 out of 70 type II genes in the rubber tree are transcribed, and they are mostly predominantly expressed in flowers, but some genes are predominantly expressed in laticifers, suggesting their roles in both flower and laticifer development. The number of type I genes in the rubber tree is only 27 (27.8%), a much smaller number compared to their orthologous genes in Arabidopsis (56.3%) and rice (43.2%). At the same time, most of the type I genes (55.6%, 15) in the rubber tree are silent and are probably pseudogenes. The high birth rate and low death rate of type II genes and low birth rate and high death rate of type I genes may corresponds to special developmental requirements in the rubber tree, e.g. the development of laticifer system for biosynthesis of cis-polyisoprene, the rubber. Moreover, atypical MIKC* factors (e.g. HbMADS1 in S-clade, and HbMADS20 in P-clade) are identified. These genes are diverged to typical MIKC* genes in sequences and facilitate functions required in laticifer development and rubber biosynthesis, which is not necessary in Arabidopsis and rice.

MADS-Box Genes Are Key Components of Genetic Regulatory Networks Involved in Abiotic Stress and Plastic Developmental Responses in Plants

Plants, as sessile organisms, adapt to different stressful conditions, such as drought, salinity, extreme temperatures, and nutrient deficiency, via plastic developmental and growth responses. Depending on the intensity and the developmental phase in which it is imposed, a stress condition may lead to a broad range of responses at the morphological, physiological, biochemical, and molecular levels. Transcription factors are key components of regulatory networks that integrate environmental cues and concert responses at the cellular level, including those that imply a stressful condition. Despite the fact that several studies have started to identify various members of the MADS-box gene family as important molecular components involved in different types of stress responses, we still lack an integrated view of their role in these processes. In this review, we analyze the function and regulation of MADS-box gene family members in response to drought, salt, cold, heat, and oxidative stress conditions in different developmental processes of several plants. In addition, we suggest that MADS-box genes are key components of gene regulatory networks involved in plant responses to stress and plant developmental plasticity in response to seasonal changes in environmental conditions.

Evolution and expression analyses of the MADS-box gene family in Brassica napus

MADS-box transcription factors are important for plant growth and development, and hundreds of MADS-box genes have been functionally characterized in plants. However, less is known about the functions of these genes in the economically important allopolyploid oil crop, Brassica napus. We identified 307 potential MADS-box genes (BnMADSs) in the B. napus genome and categorized them into type I (Mα, Mβ, and Mγ) and type II (MADS DNA-binding domain, intervening domain, keratin-like domain, and C-terminal domain [MIKC]c and MIKC*) based on phylogeny, protein motif structure, and exon-intron organization. We identified one conserved intron pattern in the MADS-box domain and seven conserved intron patterns in the K-box domain of the MIKCc genes that were previously ignored and may be associated with function. Chromosome distribution and synteny analysis revealed that hybridization between Brassica rapa and Brassica oleracea, segmental duplication, and homologous exchange (HE) in B. napus were the main BnMADSs expansion mechanisms. Promoter cis-element analyses indicated that BnMADSs may respond to various stressors (drought, heat, hormones) and light. Expression analyses showed that homologous genes in a given subfamily or sister pair are highly conserved, indicating widespread functional conservation and redundancy. Analyses of BnMADSs provide a basis for understanding their functional roles in plant development.

A well-resolved fern nuclear phylogeny reveals the evolution history of numerous transcription factor families

Ferns account for 80% of nonflowering vascular plant species and are the sister lineage of seed plants. Recent molecular phylogenetics have greatly advanced understanding of fern tree of life, but relationships among some major lineages remain unclear. To better resolve the phylogenetic relationships of ferns, we generated transcriptomes from 125 ferns and two lycophytes, with three additional public datasets, to represent all 11 orders and 85% of families of ferns. Our nuclear phylogeny provides strong supports for the monophyly of all four subclasses and nearly all orders and families, and for relationships among these lineages. The only exception is Gleicheniales, which was highly supported as being paraphyletic with Dipteridaceae sister to a clade with Gleicheniaceae + Hymenophyllales. In addition, new and strongly supported phylogenetic relationships are found for suborders and families in Polypodiales. We provide the first dated fern phylogenomic tree using many nuclear genes from a large majority of families, with an estimate for separation of the ancestors of ferns and seed plants in early Devonian at ∼400 Mya and subsequent gradual divergences of fern orders from ∼380 to 200 Mya. Moreover, the newly obtained fern phylogeny provides a framework for gene family analyses, which indicate that the vast majority of transcription factor families found in seed plants were already present in the common ancestor of extant vascular plants. In addition, fern transcription factor genes show similar duplication patterns to those in seed plants, with some showing stable copy number and others displaying independent expansions in both ferns and seed plants. This study provides a robust phylogenetic and gene family evolution framework, as well as rich molecular resources for understanding the morphological and functional evolution in ferns.