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

Genome-Wide Identification of Trehalose-6-phosphate Synthase (TPS) Gene Family Reveals the Potential Role in Carbohydrate Metabolism in Peach

Trehalose-6-phosphate synthase (TPS) is essential for plant growth and development, linking trehalose-6-phosphate (T6P) to carbon metabolism. However, little is known about the TPS gene family in peaches and their potential roles in regulating carbohydrates in peach fruit. In this study, nine TPS genes were identified in the peach genome and named according to the homologous genes in Arabidopsis. Phylogenetic analysis showed that three subfamilies were identified, including TPSI, TPSII-1, and TPSII-2, which were also consistent with gene structure analysis. Considerable cis-elements were enriched in the promoters, including plant hormone-related elements. Tissue-specific analysis showed that these TPS genes were mainly expressed in leaves, stems, and fruit, showing different expression patterns for each gene. In addition, during fruit development, the content of trehalose-6-phosphate (T6P) was positively correlated with the expression of PpTPS7a and negatively with sucrose non-fermenting-1-related kinase 1 (SnRK1) activity. Transient overexpression and silencing of PpTPS7a in peach fruit validated its function in regulating T6P content and SnRK1 activity.

Functional Divergence Analysis of AGL6 Genes in Prunus mume

The AGAMOUS-LIKE6 (AGL6) lineage is an important clade of MADS-box transcription factors that play essential roles in floral organ development. The genome of Prunus mume contains two homoeologous AGL6 genes that are replicated as gene fragments. In this study, two AGL6 homologs, PmAGL6-1 and PmAGL6-2, were cloned from P. mume and then functionally characterized. Sequence alignment and phylogenetic analyses grouped both genes into the AGL6 lineage. The expression patterns and protein-protein interaction patterns showed significant differences between the two genes. However, the ectopic expression of the two genes in Arabidopsis thaliana resulted in similar phenotypes, including the promotion of flowering, alteration of floral organ structure, participation in the formation of the floral meristem and promotion of pod bending. Therefore, gene duplication has led to some functional divergence of PmAGL6-1 and PmAGL6-2 but their functions are similar. We thus speculated that AGL6 genes play a crucial role in flower development in P. mume.

Molecular evolution and structural analyses of proteins involved in metabolic pathways of volatile organic compounds in Petunia hybrida (Solanaceae)

The association between plants and their pollinators is essential for increasing the diversity in angiosperms. Morphological and physiological traits, mainly floral scent, can influence the pollination dynamics and select pollinators for each plant species. In this work, we studied two proteins involved in producing volatile organic compounds in plants, conyferyl alcohol acyltransferase (CFAT) and benzoyl-CoA:benzyl alcohol/phenyl ethanol benzoyl transferase (BPBT) genes. We aimed to understand these proteins with respect to evolutionary and structural aspects and functions in Solanaceae using phylogenetic methods and comparative molecular modeling. We used Bayesian inference to describe the proteins' evolutionary history using Petunia x hybrida as a query to search for homologs in the Solanaceae family. Theoretical 3D models were obtained for both proteins using Panicum virgatum as a template. The phylogenetic tree included several different enzymes with diverse biological roles in Solanaceae, displaying the transferase domain. We identified only one sequence of CFAT in the databases, which belongs to Petunia x hybrida, and found several BPBT sequences from the genera Nicotiana, Solanum, and Capsicum. The 3D structures of CFAT and BPBT have two different domains, and we have identified the amino acid residues essential for the enzymatic activity and interaction with substrates.

Diverse function of the PISTILLATA, APETALA 3, and AGAMOUS-like MADS-box genes involved in the floral development in Alpinia hainanensis (Zingiberaceae)

Zingiberaceae is the vital clue and key node in the decreased process of fertile stamens in Zingiberales, helping to understand the evolution of the ginger families. This study focuses on Alpinia hainanensis to investigate the function of B- and C-class MADS-box genes in floral development. The introns size of two B-class genes AhPI and AhAP3, and one C-class gene AhAG are quite variable. By contrast, the positions of the corresponding introns are conserved, resulting in a similar exon size in homologs. The typical region 70 bp-CCAATCA element was not found in the second intron of AhAG compared to AG homologs. The subcellular localization showed that AhAP3 was in both intranuclear and extranuclear. The heterodimer was formed between APETALA3 and PISTILLATA but not between the B- and C-class proteins using Y2H and BiFC. The 35S::AhAG heterologous transformed Arabidopsis had curly and smaller rosette leaves with early flowering. Floral organs had no homeotic conversion, albeit sepals and petals reduced in size. Siliques development was affected and displayed wrinkled and shorter. By contrast, 35S::AhAP3 and 35S::AhPI did not show any modified phenotype in transgenic Arabidopsis thaliana. We first proposed the model for Alpinia flower development. MADS-box transcription factor binding at particular genomic locations and interaction with partners may be crucial for the development of the floral organ.

Genome-wide identification of MIKC-type genes related to stamen and gynoecium development in Liriodendron

The organogenesis and development of reproductive organs, i.e., stamen and gynoecium, are important floral characteristics that are closely related to pollinators and reproductive fitness. As a genus from Magnoliaceae, Liriodendron has only two relict species: L. chinense and L. tulipifera. Despite the similar flower shapes of these species, their natural seed-setting rates differ significantly, implying interspecies difference in floral organogenesis and development. MADS-box genes, which participate in floral organogenesis and development, remain unexplored in Liriodendron. Here, to explore the interspecies difference in floral organogenesis and development and identify MADS-box genes in Liriodendron, we examined the stamen and gynoecium primordia of the two Liriodendron species by scanning electron microscopy combined with paraffin sectioning, and then collected two types of primordia for RNA-seq. A total of 12 libraries were constructed and 42,268 genes were identified, including 35,269 reference genes and 6,999 new genes. Monoterpenoid biosynthesis was enriched in L. tulipifera. Genome-wide analysis of 32 MADS-box genes was conducted, including phylogenetic trees, exon/intron structures, and conserved motif distributions. Twenty-six genes were anchored on 17 scaffolds, and six new genes had no location information. The expression profiles of MIKC-type genes via RT-qPCR acrossing six stamen and gynoecium developmental stages indicates that the PI-like, AG/STK-like, SEP-like, and SVP-like genes may contribute to the species-specific differentiation of the organogenesis and development of reproductive organs in Liriodendron. Our findings laid the groundwork for the future exploration of the mechanism underlying on the interspecific differences in reproductive organ development and fitness in Liriodendron.

Comprehensive multiomics analysis reveals key roles of NACs in plant growth and development and its environmental adaption mechanism by regulating metabolite pathways

Abnormal environmental conditions induce polyploidization and exacerbate vulnerability to agricultural production. Polyploidization is a pivotal event for plant adaption to stress and the expansion of transcription factors. NACs play key roles in plant stress resistance and growth and development, but the adaptive mechanism of NACs during plant polyploidization remain to be explored. Here, we identified and analyzed NACs from 15 species and found that the expansion of NACs was contributed by polyploidization. The regulatory networks were systematically analyzed based on polyomics. NACs might influence plant phenotypes and were correlated with amino acids acting as nitrogen source, indicating that NACs play a vital role in plant development. More importantly, in quinoa and Arabidopsis thaliana, NACs enabled plants to resist stress by regulating flavonoid pathways, and the universality was further confirmed by the Arabidopsis population. Our study provides a cornerstone for future research into improvement of important agronomic traits by transcription factors in a changing global environment.

Can the anatomy of abnormal flowers elucidate relationships of the androecial members in the ginger (Zingiberaceae)?

Background

Interpretation of the floral structure of Zingiberaceae has long concentrated on the relationships of the androecial members. It suggested that labellum is composed of two structures rather than three or five, and glands are interpreted either as gynoecial part or as androecial members.

Methods

Serial sections were used to observe the vasculature of normal and two-staminate flowers in Alpinia intermedia ‘shengzhen’. Floral diagrams were drawn to interpret the morphological category of the floral organs and the relationships of the androecial members. Androecial vascular bundles were associated with carpellary dorsal bundles (CDBs) and parietal bundles (PBs) in a Zingiberales phylogeny setting using ancestral state reconstruction.

Results

Anatomical observations demonstrate that the fertile stamen(s) incorporate parietal bundles both in normal and two-staminate flowers. The three appendages represent the three members of the outer whorl of the androecium, while the labellum represents the inner whorl of the androecium in the two-staminate flower. Reconstruction of the origin of the vascular system in the androecium suggests that the outer whorl of androecium receives its vascular supply from the CDBs, and the inner whorl of androecium receives from the PBs in both the basal banana group and the more derived ginger clade.

Conclusions

The present study adds to a growing body of literature suggesting that the anatomy of abnormal flowers may not provide enough evidence for elucidating the relationships of the androecial members, and help us to better understand how the vascular system is constructed during the androecial petaloidy evolution.

Conservation and divergence of ancestral AGAMOUS/SEEDSTICK subfamily genes from the basal angiosperm Magnolia wufengensis

AGAMOUS/SEEDSTICK (AG/STK) subfamily genes play crucial roles in the reproductive development of plants. However, most of our current knowledge of AG/STK subfamily genes is restricted to core eudicots and grasses, and the knowledge of ancestral exon-intron structures, expression patterns, protein-protein interaction patterns and functions of AG/STK subfamily genes remains unclear. To determine these, we isolated AG/STK subfamily genes (MawuAG1, MawuAG2 and MawuSTK) from a woody basal angiosperm Magnolia wufengensis (Magnoliaceae). MawuSTK arose from the gene duplication event occurring before the diversification of extant angiosperms, and MawuAG1 and MawuAG2 may result from a gene duplication event occurring before the divergence of Magnoliaceae and Lauraceae. Gene duplication led to apparent diversification in their expression and interaction patterns. It revealed that expression in both stamens and carpels likely represents the ancestral expression profiles of AG lineage genes, and expression of STK-like genes in stamens may have been lost soon after the appearance of the STK lineage. Moreover, AG/STK subfamily proteins may have immediately established interactions with the SEPALLATA (SEP) subfamily proteins following the emergence of the SEP subfamily; however, their interactions with the APETALA1/FRUITFULL subfamily proteins or themselves differ from those found in monocots and basal and core eudicots. MawuAG1 plays highly conserved roles in the determinacy of stamen, carpel and ovule identity, while gene duplication contributed to the functional diversification of MawuAG2 and MawuSTK. In addition, we investigated the evolutionary history of exon-intron structural changes of the AG/STK subfamily, and a novel splice-acceptor mode (GUU-AU) and the convergent evolution of N-terminal extension in the euAG and PLE subclades were revealed for the first time. These results further advance our understanding of ancestral AG/STK subfamily genes in terms of phylogeny, exon-intron structures, expression and interaction patterns, and functions, and provide strong evidence for the significance of gene duplication in the expansion and evolution of the AG/STK subfamily.

Gene duplication led to divergence of expression patterns, protein-protein interaction patterns and floral development functions of AGL6-like genes in the basal angiosperm Magnolia wufengensis (Magnoliaceae)

The MADS-box family genes play critical roles in the regulation of growth and development of flowering plants. AGAMOUS-LIKE 6 (AGL6)-like genes are one of the most enigmatic subfamilies of the MADS-box family because of highly variable expression patterns and ambiguous functions, which have long puzzled researchers. A lot of AGL6 homologs have been identified from gymnosperms and angiosperms. However, only a few have been characterized, especially for basal angiosperm taxa. Magnolia wufengensis is a woody basal angiosperm from the family Magnoliaceae. In the current study, the phylogenesis, expression and protein-protein interaction (PPI) patterns, and functions of two AGL6 homologs from M. wufengensis, MawuAGL6-1 and MawuAGL6-2, were analyzed. Phylogenetic analysis indicated that the two AGL6 duplicates may have arisen by gene duplication before the divergence of Magnoliaceae and Lauraceae, with the diversification of their expression and PPI patterns after gene duplication. Functional analysis revealed that, in addition to common functions in accelerating flowering, MawuAGL6-1 might be responsible for flower meristem determinacy, while MawuAGL6-2 is preferentially recruited to regulate tepal morphogenesis. These findings further advance our understanding of the evolution of phylogenesis, expression, interaction and functions of AGL6 lineage genes from basal angiosperms, as well as the entire AGL6 lineage genes, and the significance of AGL6 lineage genes in the evolution and biological diversity.

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.

Floral MADS-box protein interactions in the early diverging angiosperm Aristolochia fimbriata Cham. (Aristolochiaceae: Piperales)

Floral identity MADS-box A, B, C, D, E, and AGL6 class genes are predominantly single copy in Magnoliids, and predate the whole genome duplication (WGD) events in monocots and eudicots. By comparison with the model species Arabidopsis thaliana, the expression patterns of B-, C-, and D-class genes in stamen, carpel, and ovules are conserved in Aristolochia fimbriata, whereas A-, E-class, and AGL6 genes have different expression patterns. Nevertheless, the interactions of these proteins that act through multimeric complexes remain poorly known in early divergent angiosperms. This study evaluates protein interactions among all floral MADS-box A. fimbriata proteins using the Yeast Two Hybrid System (Y2H). We found no homodimers and less heterodimers formed by AfimFUL when compared to AfimAGL6, which allowed us to suggest AGL6 homodimers in combination with AfimSEP2 as the most likely tetramer in sepal identity. We found AfimAP3-AfimPI obligate heterodimers and AfimAG-AfimSEP2 protein interactions intact suggesting conserved stamen and carpel tetrameric complexes in A. fimbriata. We observed a broader interaction partner set for AfimSEP2 than for its paralog AfimSEP1. We show conserved and exclusive MADS-box protein interactions in A. fimbriata in comparison with other eudicot and monocot model species in order to establish plesiomorphic MADS-box protein floral networks in angiosperms.

Comparative analysis of whole flower transcriptomes in the Zingiberales

The advancement of next generation sequencing technologies (NGS) has revolutionized our ability to generate large quantities of data at a genomic scale. Despite great challenges, these new sequencing technologies have empowered scientists to explore various relevant biological questions on non-model organisms, even in the absence of a complete sequenced reference genome. Here, we analyzed whole flower transcriptome libraries from exemplar species across the monocot order Zingiberales, using a comparative approach in order to gain insight into the evolution of the molecular mechanisms underlying flower development in the group. We identified 4,153 coding genes shared by all floral transcriptomes analyzed, and 1,748 genes that are only retrieved in the Zingiberales. We also identified 666 genes that are unique to the ginger lineage, and 2,001 that are only found in the banana group, while in the outgroup species Dichorisandra thyrsiflora J.C. Mikan (Commelinaceae) we retrieved 2,686 unique genes. It is possible that some of these genes underlie lineage-specific molecular mechanisms of floral diversification. We further discuss the nature of these lineage-specific datasets, emphasizing conserved and unique molecular processes with special emphasis in the Zingiberales. We also briefly discuss the strengths and shortcomings of de novo assembly for the study of developmental processes across divergent taxa from a particular order. Although this comparison is based exclusively on coding genes, with particular emphasis in transcription factors, we believe that the careful study of other regulatory mechanisms, such as non-coding RNAs, might reveal new levels of complexity, which were not explored in this work.

Transcriptome profile analysis reveals the regulation mechanism of floral sex differentiation in Jatropha curcas L

The seeds of Jatropha curcas contain a high percentage of biodiesel. However, low seed yield which was limited by its poor female flowers was a bottleneck for its utilization. Here, we compared the transcriptomic profiles of five different samples during floral sex differentiation stages using Illumina Hiseq 4000. Our results showed that hundreds of differentially expressed genes (DEGs) were detected in floral sex initiation period, but thousands of DEGs were involved in the stamens and ovules development process. Moreover, the DEGs were mainly shown up-regulation in male floral initiation, but mainly down-regulation in female floral initiation. Male floral initiation was associated with the flavonoid biosynthesis pathway while female floral initiation was related to the phytohormone signal transduction pathway. Cytokinin (CTK) signaling triggered the initiation of female floral primordium, thereafter other phytohormones co-promoted the female floral development. In addition, the floral organ identity genes played important roles in floral sex differentiation process and displayed a general conservation of the ABCDE model in J. curcas. To the best of our knowledge, this data is the first comprehensive analysis of the underlying regulatory mechanism and the related genes during floral sex differentiation in J. curcas, which help in engineering high-yielding varieties of J. curcas.

Change of Fate and Staminodial Laminarity as Potential Agents of Floral Diversification in the Zingiberales

The evolution of floral morphology in the monocot order Zingiberales shows a trend in which androecial whorl organs are progressively modified into variously conspicuous "petaloid" structures with differing degrees of fertility. Petaloidy of androecial members results from extensive laminarization of an otherwise radially symmetric structure. The genetic basis of the laminarization of androecial members has been addressed through recent candidate gene studies focused on understanding the spatiotemporal expression patterns of genes known to be necessary to floral organ formation. Here, we explore the correlation between gene duplication events and floral and inflorescence morphological diversification across the Zingiberales by inferring ancestral character states and gene copy number using the most widely accepted phylogenetic hypotheses. Our results suggest that the duplication and differential loss of GLOBOSA (GLO) copies is correlated with a change in the degree of the laminarization of androecial members. We also find an association with increased diversification in most families. We hypothesize that retention of paralogs in flower development genes could have led to a developmental shift affecting androecial organs with potential adaptive consequences, thus favoring diversification in some lineages but not others.

Imbalanced positive selection maintains the functional divergence of duplicated DIHYDROKAEMPFEROL 4-REDUCTASE genes

Gene duplication could be beneficial by functional division but might increase the risk of genetic load. The dynamics of duplicated paralogs number could involve recombination, positive selection, and functional divergence. Duplication of DIHYDROFLAVONOL 4-REDUCTASE (DFR) has been reported in several organisms and may have been retained by escape from adaptive conflict (EAC). In this study, we screened the angiosperm DFR gene focusing on a diversified genus Scutellaria to investigate how these duplicated genes are retained. We deduced that gene duplication involved multiple independent events in angiosperms, but the duplication of DFR was before the divergence of Scutellaria. Asymmetric positive selective pressures resulted in different evolutionary rates between the duplicates. Different numbers of regulatory elements, differential codon usages, radical amino acid changes, and differential gene expressions provide evidences of functional divergence between the two DFR duplicates in Scutellaria, implying adaptive subfunctionalization between duplicates. The discovery of pseudogenes accompanying a reduced replacement rate in one DFR paralogous gene suggested possibly leading to “loss of function” due to dosage imbalance after the transient adaptive subfunctionalization in the early stage of duplication. Notwithstanding, episodic gene duplication and functional divergence may be relevant to the diversification of ecological function of DFR gene in Scutellaria.

Temporal-Spatial Transcriptome Analyses Provide Insights into the Development of Petaloid Androecium in Canna indica

Canna indica (Zingiberales) is one of the most important ornamental species characterized with beautiful petaloid staminodes, which are considered to evolve from stamens. However, the genetic basis for the development of petaloid staminodes remains unclear largely because the genomic sequences are not available. By using RNA-Seq, we sequenced the transcripts in the flower of C. indica, and quantified the temporal gene expressions in flower primordium and differentiated flower, as well as the spatial gene expressions in petal and petaloid staminode. In total, 118,869 unigenes were assembled, among which 67,299 unigenes were annotated. Quantification analysis identified the differentially expressed genes in the temporal and spatial two comparisons, based on which, Gene Ontology enrichment analysis highlighted the representative terms in each sample, such as specification of organ number in flower primordium, growth in differentiated flower, secondary cell wall biogenesis in petal and cell division in petaloid staminode. Among the 51 analyzed MADS-box unigenes, 37 were up-regulated in differentiated flower compared with those in flower primordium. A-class unigenes were expressed higher in petal than in petaloid staminode, and C-class unigenes were expressed oppositely, whereas B-class unigenes demonstrated close expression levels in these two organs, indicating that petaloid staminode retains stamen identity to some degree. In situ hybridization provided more detailed expression patterns of these unigenes, and revealed the extended expression of B-class to the carpel at later stages when the style turned flat. These results constitute a preliminary basis for the study of flower development in C. indica and can be applied in further study of the evolution of Zingiberales.

The developmental and genetic bases of apetaly in Bocconia frutescens (Chelidonieae: Papaveraceae)

BACKGROUND

Bocconia and Macleaya are the only genera of the poppy family (Papaveraceae) lacking petals; however, the developmental and genetic processes underlying such evolutionary shift have not yet been studied.

RESULTS

We studied floral development in two species of petal-less poppies Bocconia frutescens and Macleaya cordata as well as in the closely related petal-bearing Stylophorum diphyllum. We generated a floral transcriptome of B. frutescens to identify MADS-box ABCE floral organ identity genes expressed during early floral development. We performed phylogenetic analyses of these genes across Ranunculales as well as RT-PCR and qRT-PCR to assess loci-specific expression patterns. We found that petal-to-stamen homeosis in petal-less poppies occurs through distinct developmental pathways. Transcriptomic analyses of B. frutescens floral buds showed that homologs of all MADS-box genes are expressed except for the APETALA3-3 ortholog. Species-specific duplications of other ABCE genes in B. frutescens have resulted in functional copies with expanded expression patterns than those predicted by the model.

CONCLUSIONS

Petal loss in B. frutescens is likely associated with the lack of expression of AP3-3 and an expanded expression of AGAMOUS. The genetic basis of petal identity is conserved in Ranunculaceae and Papaveraceae although they have different number of AP3 paralogs and exhibit dissimilar floral groundplans.

The developmental and genetic bases of apetaly in Bocconia frutescens (Chelidonieae: Papaveraceae)

Background

Bocconia and Macleaya are the only genera of the poppy family (Papaveraceae) lacking petals; however, the developmental and genetic processes underlying such evolutionary shift have not yet been studied.

Results

We studied floral development in two species of petal-less poppies Bocconiafrutescens and Macleayacordata as well as in the closely related petal-bearing Stylophorum diphyllum. We generated a floral transcriptome of B. frutescens to identify MADS-box ABCE floral organ identity genes expressed during early floral development. We performed phylogenetic analyses of these genes across Ranunculales as well as RT-PCR and qRT-PCR to assess loci-specific expression patterns. We found that petal-to-stamen homeosis in petal-less poppies occurs through distinct developmental pathways. Transcriptomic analyses of B. frutescens floral buds showed that homologs of all MADS-box genes are expressed except for the APETALA3-3 ortholog. Species-specific duplications of other ABCE genes in B. frutescens have resulted in functional copies with expanded expression patterns than those predicted by the model.

Conclusions

Petal loss in B. frutescens is likely associated with the lack of expression of AP3-3 and an expanded expression of AGAMOUS. The genetic basis of petal identity is conserved in Ranunculaceae and Papaveraceae although they have different number of AP3 paralogs and exhibit dissimilar floral groundplans.

Electronic supplementary material

The online version of this article (doi:10.1186/s13227-016-0054-6) contains supplementary material, which is available to authorized users.

Evolutionary dynamics and biogeography of Musaceae reveal a correlation between the diversification of the banana family and the geological and climatic history of Southeast Asia

Tropical Southeast Asia, which harbors most of the Musaceae biodiversity, is one of the most species-rich regions in the world. Its high degree of endemism is shaped by the region's tectonic and climatic history, with large differences between northern Indo-Burma and the Malayan Archipelago. Here, we aim to find a link between the diversification and biogeography of Musaceae and geological history of the Southeast Asian subcontinent. The Musaceae family (including five Ensete, 45 Musa and one Musella species) was dated using a large phylogenetic framework encompassing 163 species from all Zingiberales families. Evolutionary patterns within Musaceae were inferred using ancestral area reconstruction and diversification rate analyses. All three Musaceae genera - Ensete, Musa and Musella - originated in northern Indo-Burma during the early Eocene. Musa species dispersed from 'northwest to southeast' into Southeast Asia with only few back-dispersals towards northern Indo-Burma. Musaceae colonization events of the Malayan Archipelago subcontinent are clearly linked to the geological and climatic history of the region. Musa species were only able to colonize the region east of Wallace's line after the availability of emergent land from the late Miocene onwards.

Molecular evolution analysis of WUSCHEL-related homeobox transcription factor family reveals functional divergence among clades in the homeobox region

Gene families have been shown to play important roles in plant evolution and are associated with diversification and speciation. Genes of WUSCHEL-related homeobox family of transcription factors have important functions in plant development and are correlated with the appearance of evolutionary novelties. There are several published studies related to this family, but little is known about the relationships among the main clades in the phylogeny and the molecular evolution of the family. In this study, we obtained a well-resolved Bayesian phylogenetic tree establishing the relationships among the main clades and determining the position of Selaginella moellendorffii WOX genes. Moreover, a correlation was identified between the number of genes in the genomes and the events of whole-genome duplications. The intron-exon structure is more consistent across the modern clade, which appeared more recently in the WOX evolutionary history, and coincides with the development of higher complexity in plant species. No positive selection was detected among sites through the branches in the tree. However, with regard to the main clades, functional divergence among certain amino acids in the homeodomain region was found. Relaxed purifying selection could be the main driving force in the evolution of these genes and in agreement with some genes have been demonstrated to be functionally redundant.

Revisiting the Zingiberales: using multiplexed exon capture to resolve ancient and recent phylogenetic splits in a charismatic plant lineage

The Zingiberales are an iconic order of monocotyledonous plants comprising eight families with distinctive and diverse floral morphologies and representing an important ecological element of tropical and subtropical forests. While the eight families are demonstrated to be monophyletic, phylogenetic relationships among these families remain unresolved. Neither combined morphological and molecular studies nor recent attempts to resolve family relationships using sequence data from whole plastomes has resulted in a well-supported, family-level phylogenetic hypothesis of relationships. Here we approach this challenge by leveraging the complete genome of one member of the order, Musa acuminata, together with transcriptome information from each of the other seven families to design a set of nuclear loci that can be enriched from highly divergent taxa with a single array-based capture of indexed genomic DNA. A total of 494 exons from 418 nuclear genes were captured for 53 ingroup taxa. The entire plastid genome was also captured for the same 53 taxa. Of the total genes captured, 308 nuclear and 68 plastid genes were used for phylogenetic estimation. The concatenated plastid and nuclear dataset supports the position of Musaceae as sister to the remaining seven families. Moreover, the combined dataset recovers known intra- and inter-family phylogenetic relationships with generally high bootstrap support. This is a flexible and cost effective method that gives the broader plant biology community a tool for generating phylogenomic scale sequence data in non-model systems at varying evolutionary depths.

Prevalent Exon-Intron Structural Changes in the APETALA1/FRUITFULL, SEPALLATA, AGAMOUS-LIKE6, and FLOWERING LOCUS C MADS-Box Gene Subfamilies Provide New Insights into Their Evolution

AP1/FUL, SEP, AGL6, and FLC subfamily genes play important roles in flower development. The phylogenetic relationships among them, however, have been controversial, which impedes our understanding of the origin and functional divergence of these genes. One possible reason for the controversy may be the problems caused by changes in the exon-intron structure of genes, which, according to recent studies, may generate non-homologous sites and hamper the homology-based sequence alignment. In this study, we first performed exon-by-exon alignments of these and three outgroup subfamilies (SOC1, AG, and STK). Phylogenetic trees reconstructed based on these matrices show improved resolution and better congruence with species phylogeny. In the context of these phylogenies, we traced evolutionary changes of exon-intron structures in each subfamily. We found that structural changes have occurred frequently following gene duplication and speciation events. Notably, exons 7 and 8 (if present) suffered more structural changes than others. With the knowledge of exon-intron structural changes, we generated more reasonable alignments containing all the focal subfamilies. The resulting trees showed that the SEP subfamily is sister to the monophyletic group formed by AP1/FUL and FLC subfamily genes and that the AGL6 subfamily forms a sister group to the three abovementioned subfamilies. Based on this topology, we inferred the evolutionary history of exon-intron structural changes among different subfamilies. Particularly, we found that the eighth exon originated before the divergence of AP1/FUL, FLC, SEP, and AGL6 subfamilies and degenerated in the ancestral FLC-like gene. These results provide new insights into the origin and evolution of the AP1/FUL, FLC, SEP, and AGL6 subfamilies.

Analysis of the floral MADS-box genes from monocotyledonous Trilliaceae species indicates the involvement of SEPALLATA3-like genes in sepal-petal differentiation

The evolution of greenish sepals from petaloid outer tepals has occurred repeatedly in various lineages of non-grass monocots. Studies in distinct monocot species showed that the evolution of sepals could be explained by the ABC model; for example, the defect of B-class function in the outermost whorl was linked to the evolution of sepals. Here, floral MADS-box genes from three sepal-bearing monocotyledonous Trilliaceae species, Trillium camschatcense, Paris verticillata, and Kinugasa japonica were examined. Unexpectedly, expression of not only A- but also B-class genes was detected in the sepals of all three species. Although the E-class gene is generally expressed across all floral whorls, no expression was detected in sepals in the three species examined here. Overexpression of the E-class SEPALLATA3-like gene from T. camschatcense (TcamSEP) in Arabidopsis thaliana produced phenotypes identical to those reported for orthologs in other monocots. Additionally, yeast hybrid experiments indicated that TcamSEP could form a higher-order complex with an endogenous heterodimer of B-class APETALA3/DEFICIENS-like (TcamDEF) and PISTILLATA/GLOBOSA-like (TcamGLO) proteins. These results suggest a conserved role for Trilliaceae SEPALLATA3-like genes in functionalization of the B-class genes, and that a lack of SEPALLATA3-like gene expression in the outermost whorl may be related to the formation of greenish sepals.

Interactions of OsMADS1 with Floral Homeotic Genes in Rice Flower Development

During reproductive development, rice plants develop unique flower organs which determine the final grain yield. OsMADS1, one of SEPALLATA-like MADS-box genes, has been unraveled to play critical roles in rice floral organ identity specification and floral meristem determinacy. However, the molecular mechanisms underlying interactions of OsMADS1 with other floral homeotic genes in regulating flower development remains largely elusive. In this work, we studied the genetic interactions of OsMADS1 with B-, C-, and D-class genes along with physical interactions among their proteins. We show that the physical and genetic interactions between OsMADS1 and OsMADS3 are essential for floral meristem activity maintenance and organ identity specification; while OsMADS1 physically and genetically interacts with OsMADS58 in regulating floral meristem determinacy and suppressing spikelet meristem reversion. We provided important genetic evidence to support the neofunctionalization of two rice C-class genes (OsMADS3 and OsMADS58) during flower development. Gene expression profiling and quantitative RT-PCR analyses further revealed that OsMADS1 affects the expression of many genes involved in floral identity and hormone signaling, and chromatin immunoprecipitation (ChIP)-PCR assay further demonstrated that OsMADS17 is a direct target gene of OsMADS1. Taken together, these results reveal that OsMADS1 has diversified regulatory functions in specifying rice floral organ and meristem identity, probably through its genetic and physical interactions with different floral homeotic regulators.

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.

Adaptation in flower form: a comparative evodevo approach

Evolutionary developmental biology (evodevo) attempts to explain how the process of organismal development evolves, utilizing a comparative approach to investigate changes in developmental pathways and processes that occur during the evolution of a given lineage. Evolutionary genetics uses a population approach to understand how organismal changes in form or function are linked to underlying genetics, focusing on changes in gene and genotype frequencies within populations and the fixation of genotypic variation into traits that define species or evoke speciation events. Microevolutionary processes, including mutation, genetic drift, natural selection and gene flow, can provide the foundation for macroevolutionary patterns observed as morphological evolution and adaptation. The temporal element linking microevolutionary processes to macroevolutionary patterns is development: an organism's genotype is converted to phenotype by ontogenetic processes. Because selection acts upon the phenotype, the connection between evolutionary genetics and developmental evolution becomes essential to understanding adaptive evolution in organismal form and function. Here, we discuss how developmental genetic studies focused on key developmental processes could be linked within a comparative framework to study the developmental genetics of adaptive evolution, providing examples from research on two key processes of plant evodevo - floral symmetry and organ fusion - and their role in the adaptation of floral form.

Interpreting lemma and palea homologies: a point of view from rice floral mutants

In contrast to eudicot flowers which typically exhibit sepals and petals at their periphery, the flowers of grasses are distinguished by the presence of characteristic outer organs. In place of sepals, grasses have evolved the lemma and the palea, two bract-like structures that partially or fully enclose the inner reproductive organs. With little morphological similarities to sepals, whether the lemma and palea are part of the perianth or non-floral organs has been a longstanding debate. In recent years, comparative studies of floral mutants as well as the availability of whole genome sequences in many plant species have provided strong arguments in favor of the hypothesis of lemma and palea being modified sepals. In rice, a feature of the palea is the bending of its lateral region into a hook-shaped marginal structure. This allows the palea to lock into the facing lemma region, forming a close-fitting lemma-palea enclosure. In this article, we focus on the rice lemma and palea and review some of the key transcription factors involved in their development and functional specialization. Alternative interpretations of these organs are also addressed.

Functional conservation and divergence of four ginger AP1/AGL9 MADS-box genes revealed by analysis of their expression and protein-protein interaction, and ectopic expression of AhFUL gene in Arabidopsis

Alpinia genus are known generally as ginger-lilies for showy flowers in the ginger family, Zingiberaceae, and their floral morphology diverges from typical monocotyledon flowers. However, little is known about the functions of ginger MADS-box genes in floral identity. In this study, four AP1/AGL9 MADS-box genes were cloned from Alpinia hainanensis, and protein-protein interactions (PPIs) and roles of the four genes in floral homeotic conversion and in floral evolution are surveyed for the first time. AhFUL is clustered to the AP1 lineage, AhSEP4 and AhSEP3b to the SEP lineage, and AhAGL6-like to the AGL6 lineage. The four genes showed conserved and divergent expression patterns, and their encoded proteins were localized in the nucleus. Seven combinations of PPI (AhFUL-AhSEP4, AhFUL-AhAGL6-like, AhFUL-AhSEP3b, AhSEP4-AhAGL6-like, AhSEP4-AhSEP3b, AhAGL6-like-AhSEP3b, and AhSEP3b-AhSEP3b) were detected, and the PPI patterns in the AP1/AGL9 lineage revealed that five of the 10 possible combinations are conserved and three are variable, while conclusions cannot yet be made regarding the other two. Ectopic expression of AhFUL in Arabidopsis thaliana led to early flowering and floral organ homeotic conversion to sepal-like or leaf-like. Therefore, we conclude that the four A. hainanensis AP1/AGL9 genes show functional conservation and divergence in the floral identity from other MADS-box genes.

Gene duplication as a driver of plant morphogenetic evolution

Most duplicated genes (paralogs) are quickly erased during evolution, and only some are retained. Yet, gene and genome duplications are connected to the evolution of genetic and, in turn, morphological complexity. Plants are especially prone to experience polyploidizations and to enhance their gene repertoire after such events. Genes encoding proteins involved in transcriptional regulation are of especial interest since they are correlated with the occurrence of genome duplication events and with the rise of plant morphological complexity. Here, I review what we know about paralog retention as a driver for morphogenetic evolution of plants. The main focus is on the evolution of plant genes controlling development (morphogenetic transcription factors).

A method for extracting high-quality RNA from diverse plants for next-generation sequencing and gene expression analyses

PREMISE OF THE STUDY

To study gene expression in plants, high-quality RNA must be extracted in quantities sufficient for subsequent cDNA library construction. Field-based collections are often limited in quantity and quality of tissue and are typically preserved in RNAlater. Obtaining sufficient and high-quality yield from variously preserved samples is essential to studies of comparative biology. We present a protocol for the extraction of high-quality RNA from even the most recalcitrant plant tissues. •

METHODS AND RESULTS

Tissues from mosses, cycads, and angiosperm floral organs and leaves were preserved in RNAlater or frozen fresh at -80°C. Extractions were performed and quality was measured for yield and purity. •

CONCLUSIONS

This protocol results in the extraction of high-quality RNA from a variety of plant tissues representing vascular and nonvascular plants. RNA was used for cDNA synthesis to generate libraries for next-generation sequencing and for expression studies using quantitative PCR (qPCR) and semiquantitative reverse transcription PCR (RT-PCR).