Evolutionary trends in the floral transcriptome: insights from one of the basalmost angiosperms, the water lily Nuphar advena (Nymphaeaceae)

Model Species to Investigate the Origin of Flowers

The angiosperms, or flowering plants, arose at least 135 million years ago (Ma) and rapidly diversified to form over 300,000 species alive today. This group appears, however, to have separated from its closest living relatives, the extant gymnosperms, much earlier: over 300 Ma. Representatives of basally-diverging angiosperm lineages are of key importance to studies aimed at reconstructing the most recent common ancestor of living angiosperms, including its morphological, anatomical, eco-physiological and molecular aspects. Furthermore, evo-devo comparisons of angiosperms with living gymnosperms may help to determine how the many novel aspects of angiosperms, including those of the flower, first came about. This chapter reviews literature on the origin of angiosperms and focusses on basally-diverging angiosperms and gymnosperms that show advantages as potential experimental models, reviewing information and protocols for the use of these species in an evo-devo context. The final section suggests a means by which data from living and fossil groups could be integrated to better elucidate evolutionary events that took place on the long stem-lineage that apparently preceded the radiation of living angiosperms.

NBLAST: a graphical user interface-based two-way BLAST software with a dot plot viewer

BLAST, a basic bioinformatics tool for searching local sequence similarity, has been one of the most widely used bioinformatics programs since its introduction in 1990. Users generally use the web-based NCBI-BLAST program for BLAST analysis. However, users with large sequence data are often faced with a problem of upload size limitation while using the web-based BLAST program. This proves inconvenient as scientists often want to run BLAST on their own data, such as transcriptome or whole genome sequences. To overcome this issue, we developed NBLAST, a graphical user interface-based BLAST program that employs a two-way system, allowing the use of input sequences either as "query" or "target" in the BLAST analysis. NBLAST is also equipped with a dot plot viewer, thus allowing researchers to create custom database for BLAST and run a dot plot similarity analysis within a single program. It is available to access to the NBLAST with http://nbitglobal.com/nblast.

Disentangling historical signal and pollinator selection on the micromorphology of flowers: an example from the floral epidermis of the Nymphaeaceae

The family Nymphaeaceae includes most of the diversity among the ANA-grade angiosperms. Among the species of this family, floral structures and pollination strategies vary. The genus Victoria, as well as subgenera Lotos and Hydrocallis in Nymphaea, present night-blooming, scented flowers pollinated by scarab beetles. Such similar pollination strategies have led to macromorphological similarities among the flowers of these species, which could be interpreted as homologies or convergences based on different phylogenetic hypotheses about the relationships of these groups. We employed scanning electron microscopy of floral epidermis for seven species of the Nymphaeaceae with contrasting pollination biology to identify the main characters of the floral organs and the potential homologous nature of the structures involved in pollinator attraction. Moreover, we used transmission electron microscopy to observe ultrastructure of papillate-conical epidermis in the stamen of Victoria cruziana. We then tested the phylogenetic or ecological distribution of these traits using both consensus network approaches and ancestral state reconstruction on fixed phylogenies. Our results show that the night-blooming flowers present different specialisations in their epidermis, with V. cruziana presenting the most elaborate floral anatomy. We also identify for the first time the presence of conical-papillate cells in the order Nymphaeales. The epidermal characters tend to reflect phylogenetic relationships more than convergence due to pollinator selection. These results point to an independent and parallel evolution of scarab pollination in Nymphaeaceae and demonstrate the promise of floral anatomy as a phylogenetic marker. Moreover, they indicate a degree of sophistication in the anatomical basis of cantharophilous flowers in the Nymphaeales that diverges from the most simplistic views of floral evolution in the angiosperms.

A Global View of Transcriptome Dynamics During Male Floral Bud Development in Populus tomentosa

To obtain a comprehensive overview of the dynamic transcriptome during male floral bud development in Populus tomentosa, high-throughput RNA-seq was conducted during eight flowering-related stages. Among the 109,212 de novo assembled unigenes, 6,959 were differentially expressed during the eight stages. The overrepresented classed of genes identified by Gene Ontology (GO) enrichment included 'response to environmental stimuli' and 'plant-type spore development'. One-third of the differentially expressed genes were transcription factors (TFs). Several genes and gene families were analyzed in depth, including MADS-box TFs, Squamosa promoter binding protein-like family, receptor-like kinases, FLOWERING LOCUS T/TERMINAL-FLOWER-LIKE 1 family, key genes involved in anther and tapetum development, as well as LEAFY, WUSCHEL and CONSTANS. The results provided new insights into the roles of these and other well known gene families during the annual flowering cycle. To explore the mechanisms regulating poplar flowering, a weighted gene co-expression network was constructed using 98 floral-related genes involved in flower meristem identity and flower development. Many modules of co-expressed genes and hub genes were identified, such as APETALA1 and HUA1. This work provides many new insights on the annual flowering cycle in a perennial plant, and a major new resource for plant biology and biotechnology.

Evolution of floral diversity: genomics, genes and gamma

A salient feature of flowering plant diversification is the emergence of a novel suite of floral features coinciding with the origin of the most species-rich lineage, Pentapetalae. Advances in phylogenetics, developmental genetics and genomics, including new analyses presented here, are helping to reconstruct the specific evolutionary steps involved in the evolution of this clade. The enormous floral diversity among Pentapetalae appears to be built on a highly conserved ground plan of five-parted (pentamerous) flowers with whorled phyllotaxis. By contrast, lability in the number and arrangement of component parts of the flower characterize the early-diverging eudicot lineages subtending Pentapetalae. The diversification of Pentapetalae also coincides closely with ancient hexaploidy, referred to as the gamma whole-genome triplication, for which the phylogenetic timing, mechanistic details and molecular evolutionary consequences are as yet not fully resolved. Transcription factors regulating floral development often persist in duplicate or triplicate in gamma-derived genomes, and both individual genes and whole transcriptional programmes exhibit a shift from broadly overlapping to tightly defined expression domains in Pentapetalae flowers. Investigations of these changes associated with the origin of Pentapetalae can lead to a more comprehensive understanding of what is arguably one of the most important evolutionary diversification events within terrestrial plants.This article is part of the themed issue 'Evo-devo in the genomics era, and the origins of morphological diversity'.

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.

Perianth organs in Nymphaeaceae: comparative study on epidermal and structural characters

The perianth organs of six species of Nymphaeaceae, representing Euryale, Nymphaea and Victoria, were studied on the basis of macroscopical, micromorphological, and anatomical characters. The aims were to determine whether perianth is differentiated among tepal whorls considering the presence of sepaloid and petaloid characters, and to evaluate the occurrence of both features in individual tepals. Selected perianth series were examined macroscopically, with light microscopy, and scanning electron microscopy. Osmophores were detected using neutral red and Sudan. In all tepals examined, stomata and hydropotes were present on the abaxial and adaxial surfaces. These are anomocytic or stephanocytic; hydropotes of irregular type are also present. The outer series of tepals display morpho-anatomical characters in most part related with photosynthetic and protective functions. Osmophore activity is very scarce and petaloid epidermal morphology is present only in N. lotus, thus allowing interpretation of this whorl as primarily sepaloid. The second series exhibits both petal-like and sepal-like characters; in N. amazonum and N. gardneriana sepaloid and petaloid group of cells are present on the abaxial surface of individual tepals. Therefore, this whorl is transitional between the outer and the innermost ones. Both the morpho-anatomy and presence of osmophore activity indicate that the innermost series is entirely petaloid. Inner tepals of E. ferox, N. alba, and V. cruziana share the presence of epidermal cells with predominantly smooth cuticle, whereas those of N. amazonum, N. gardneriana, and N. lotus share a cuticular ornamentation consisting of numerous papillae on each cell. Morphological characters of the perianth epidermis are in some respects congruent with the molecular phylogeny of Nymphaeaceae. Our results support the co-expression of sepaloidy and petaloidy within individual tepals and the mosaic model of perianth evolution proposed for the angiosperms.

Comparative transcriptome analyses of flower development in four species of Achimenes (Gesneriaceae)

BACKGROUND

Flowers have an amazingly diverse display of colors and shapes, and these characteristics often vary significantly among closely related species. The evolution of diverse floral form can be thought of as an adaptive response to pollination and reproduction, but it can also be seen through the lens of morphological and developmental constraints. To explore these interactions, we use RNA-seq across species and development to investigate gene expression and sequence evolution as they relate to the evolution of the diverse flowers in a group of Neotropical plants native to Mexico-magic flowers (Achimenes, Gesneriaceae).

RESULTS

The assembled transcriptomes contain between 29,000 and 42,000 genes expressed during development. We combine sequence orthology and coexpression clustering with analyses of protein evolution to identify candidate genes for roles in floral form evolution. Over 25% of transcripts captured were distinctive to Achimenes and overrepresented by genes involved in transcription factor activity. Using a model-based clustering approach we find dynamic, temporal patterns of gene expression among species. Selection tests provide evidence of positive selection in several genes with roles in pigment production, flowering time, and morphology. Combining these approaches to explore genes related to flower color and flower shape, we find distinct patterns that correspond to transitions of floral form among Achimenes species.

CONCLUSIONS

The floral transcriptomes developed from four species of Achimenes provide insight into the mechanisms involved in the evolution of diverse floral form among closely related species with different pollinators. We identified several candidate genes that will serve as an important and useful resource for future research. High conservation of sequence structure, patterns of gene coexpression, and detection of positive selection acting on few genes suggests that large phenotypic differences in floral form may be caused by genetic differences in a small set of genes. Our characterized floral transcriptomes provided here should facilitate further analyses into the genomics of flower development and the mechanisms underlying the evolution of diverse flowers in Achimenes and other Neotropical Gesneriaceae.

Evolving Ideas on the Origin and Evolution of Flowers: New Perspectives in the Genomic Era

The origin of the flower was a key innovation in the history of complex organisms, dramatically altering Earth's biota. Advances in phylogenetics, developmental genetics, and genomics during the past 25 years have substantially advanced our understanding of the evolution of flowers, yet crucial aspects of floral evolution remain, such as the series of genetic and morphological changes that gave rise to the first flowers; the factors enabling the origin of the pentamerous eudicot flower, which characterizes ∼70% of all extant angiosperm species; and the role of gene and genome duplications in facilitating floral innovations. A key early concept was the ABC model of floral organ specification, developed by Elliott Meyerowitz and Enrico Coen and based on two model systems,Arabidopsis thalianaandAntirrhinum majus Yet it is now clear that these model systems are highly derived species, whose molecular genetic-developmental organization must be very different from that of ancestral, as well as early, angiosperms. In this article, we will discuss how new research approaches are illuminating the early events in floral evolution and the prospects for further progress. In particular, advancing the next generation of research in floral evolution will require the development of one or more functional model systems from among the basal angiosperms and basal eudicots. More broadly, we urge the development of "model clades" for genomic and evolutionary-developmental analyses, instead of the primary use of single "model organisms." We predict that new evolutionary models will soon emerge as genetic/genomic models, providing unprecedented new insights into floral evolution.

Water lilies as emerging models for Darwin's abominable mystery

Water lilies are not only highly favored aquatic ornamental plants with cultural and economic importance but they also occupy a critical evolutionary space that is crucial for understanding the origin and early evolutionary trajectory of flowering plants. The birth and rapid radiation of flowering plants has interested many scientists and was considered 'an abominable mystery' by Charles Darwin. In searching for the angiosperm evolutionary origin and its underlying mechanisms, the genome of Amborella has shed some light on the molecular features of one of the basal angiosperm lineages; however, little is known regarding the genetics and genomics of another basal angiosperm lineage, namely, the water lily. In this study, we reviewed current molecular research and note that water lily research has entered the genomic era. We propose that the genome of the water lily is critical for studying the contentious relationship of basal angiosperms and Darwin's 'abominable mystery'. Four pantropical water lilies, especially the recently sequenced Nymphaea colorata, have characteristics such as small size, rapid growth rate and numerous seeds and can act as the best model for understanding the origin of angiosperms. The water lily genome is also valuable for revealing the genetics of ornamental traits and will largely accelerate the molecular breeding of water lilies.

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.

The floral transcriptome of Eucalyptus grandis

As a step toward functional annotation of genes required for floral initiation and development within the Eucalyptus genome, we used short read sequencing to analyze transcriptomes of floral buds from early and late developmental stages, and compared these with transcriptomes of diverse vegetative tissues, including leaves, roots, and stems. A subset of 4807 genes (13% of protein-coding genes) were differentially expressed between floral buds of either stage and vegetative tissues. A similar proportion of genes were differentially expressed among all tissues. A total of 479 genes were differentially expressed between early and late stages of floral development. Gene function enrichment identified 158 gene ontology classes that were overrepresented in floral tissues, including 'pollen development' and 'aromatic compound biosynthetic process'. At least 40 floral-dominant genes lacked functional annotations and thus may be novel floral transcripts. We analyzed several genes and gene families in depth, including 49 putative biomarkers of floral development, the MADS-box transcription factors, 'S-domain'-receptor-like kinases, and selected gene family members with phosphatidylethanolamine-binding protein domains. Expanded MADS-box gene subfamilies in Eucalyptus grandis included SUPPRESSOR OF OVEREXPRESSION OF CO 1 (SOC1), SEPALLATA (SEP) and SHORT VEGETATIVE PHASE (SVP) Arabidopsis thaliana homologs. These data provide a rich resource for functional and evolutionary analysis of genes controlling eucalypt floral development, and new tools for breeding and biotechnology.

Flower diversity and angiosperm diversification

The flower itself, which comprises most of the evolutionary innovations of flowering plants, bears special significance for understanding the origin and diversification of angiosperms. The sudden origin of angiosperms in the fossil record poses unanswered questions on both the origins of flowering plants and their rapid spread and diversification. Central to these questions is the role that the flower, and floral diversity, played. Recent clarifications of angiosperm phylogeny provide the foundation for investigating evolutionary transitions in floral features and the underlying genetic mechanisms of stasis and change. The general features of floral diversity can best be addressed by considering key patterns of variation: an undifferentiated versus a differentiated perianth; elaboration of perianth organs in size and color; merosity of the flower; and phyllotaxy of floral organs. Various models of gene expression now explain the regulation of floral organization and floral organ identity; the best understood are the ABC(E) model and its modifications, but other gene systems are important in specific clades and require further study. Furthermore, the propensity for gene and genome duplications in angiosperms provides abundant raw material for novel floral features--emphasizing the importance of understanding the conservation and diversification of gene lineages and functions in studies of macroevolution.

Comparative transcriptomic analysis of male and female flowers of monoecious Quercus suber

Monoecious species provide a comprehensive system to study the developmental programs underlying the establishment of female and male organs in unisexual flowers. However, molecular resources for most monoecious non-model species are limited, hampering our ability to study the molecular mechanisms involved in flower development of these species. The objective of this study was to identify differentially expressed genes during the development of male and female flowers of the monoecious species Quercus suber, an economically important Mediterranean tree. Total RNA was extracted from different developmental stages of Q. suber flowers. Non-normalized cDNA libraries of male and female flowers were generated using 454 pyrosequencing technology producing a total of 962,172 high-quality reads with an average length of 264 nucleotides. The assembly of the reads resulted in 14,488 contigs for female libraries and 10,438 contigs for male libraries. Comparative analysis of the transcriptomes revealed genes differentially expressed in early and late stages of development of female and male flowers, some of which have been shown to be involved in pollen development, in ovule formation and in flower development of other species with a monoecious, dioecious, or hermaphroditic sexual system. Moreover, we found differentially expressed genes that have not yet been characterized and others that have not been previously shown to be implicated in flower development. This transcriptomic analysis constitutes a major step toward the characterization of the molecular mechanisms involved in flower development in a monoecious tree with a potential contribution toward the knowledge of conserved developmental mechanisms in other species.

Comparative transcriptional profiling provides insights into the evolution and development of the zygomorphic flower of Vicia sativa (Papilionoideae)

Background

Vicia sativa (the common vetch) possesses a predominant zygomorphic flower and belongs to the subfamily Papilionoideae, which is related to Arabidopsis thaliana in the eurosid II clade of the core eudicots. Each vetch flower consists of 21 concentrically arranged organs: the outermost five sepals, then five petals and ten stamens, and a single carpel in the center.

Methodology/Principal Findings

We explored the floral transcriptome to examine a genome-scale genetic model of the zygomorphic flower of vetch. mRNA was obtained from an equal mixture of six floral organs, leaves and roots. De novo assembly of the vetch transcriptome using Illumina paired-end technology produced 71,553 unigenes with an average length of 511 bp. We then compared the expression changes in the 71,553 unigenes in the eight independent organs through RNA-Seq Quantification analysis. We predominantly analyzed gene expression patterns specific to each floral organ and combinations of floral organs that corresponded to the traditional ABC model domains. Comparative analyses were performed in the floral transcriptomes of vetch and Arabidopsis, and genomes of vetch and Medicago truncatula.

Conclusions/Significance

Our comparative analysis of vetch and Arabidopsis showed that the vetch flowers conform to a strict ABC model. We analyzed the evolution and expression of the TCP gene family in vetch at a whole-genome level, and several unigenes specific to three different vetch petals, which might offer some clues toward elucidating the molecular mechanisms underlying floral zygomorphy. Our results provide the first insights into the genome-scale molecular regulatory network that controls the evolution and development of the zygomorphic flower in Papilionoideae.

Inference of genome duplications from age distributions revisited

Whole-genome duplications (WGDs), thought to facilitate evolutionary innovations and adaptations, have been uncovered in many phylogenetic lineages. WGDs are frequently inferred from duplicate age distributions, where they manifest themselves as peaks against a small-scale duplication background. However, the interpretation of duplicate age distributions is complicated by the use of K(S), the number of synonymous substitutions per synonymous site, as a proxy for the age of paralogs. Two particular concerns are the stochastic nature of synonymous substitutions leading to increasing uncertainty in K(S) with increasing age since duplication and K(S) saturation caused by the inability of evolutionary models to fully correct for the occurrence of multiple substitutions at the same site. K(S) stochasticity is expected to erode the signal of older WGDs, whereas K(S) saturation may lead to artificial peaks in the distribution. Here, we investigate the consequences of these effects on K(S)-based age distributions and WGD inference by simulating the evolution of duplicated sequences according to predefined real age distributions and re-estimating the corresponding K(S) distributions. We show that, although K(S) estimates can be used for WGD inference far beyond the commonly accepted K(S) threshold of 1, K(S) saturation effects can cause artificial peaks at higher ages. Moreover, K(S) stochasticity and saturation may lead to confounded peaks encompassing multiple WGD events and/or saturation artifacts. We argue that K(S) effects need to be properly accounted for when inferring WGDs from age distributions and that the failure to do so could lead to false inferences.

[Expression analysis of the Medicago truncatula floral specific expression genes]

The expression of genes specific to floral organ is important for the floral organ formation and development in Medicago truncatula. Screening of the genes specifically expressed in M. truncatula flowers and comparing the expression patterns of their orthologous homologous genes among different model plants can provide novel insights into the functions of these genes in controlling the floral organ development in M. truncatula. According to the expression profile data of PISTILLATA (PI), we screened 97 genes specifically expressed in M. truncatula floral organs (ratio≥10 and Z≥7.9). Their homolog genes were also identified in Arabidopsis thaliana, soybean (Glycine max L.), Lotus japonicus, and rice (Oryza sativa L.). The results of comparing the gene expression levels, the gene expression patterns, and the gene functions among these species indicated that the expression variation of the orthologous homolog genes was small in the kindred species and was great in distant species. Furthermore, we compared the cis-acting regulatory elements of the genes, which had large expression variation among different plants. These results suggest that the great discrepancy of the orthologous homolog gene expression caused by the different character of cis-element in the promoter region.

Conservation and canalization of gene expression during angiosperm diversification accompany the origin and evolution of the flower

The origin and rapid diversification of the angiosperms (Darwin's "Abominable Mystery") has engaged generations of researchers. Here, we examine the floral genetic programs of phylogenetically pivotal angiosperms (water lily, avocado, California poppy, and Arabidopsis) and a nonflowering seed plant (a cycad) to obtain insight into the origin and subsequent evolution of the flower. Transcriptional cascades with broadly overlapping spatial domains, resembling the hypothesized ancestral gymnosperm program, are deployed across morphologically intergrading organs in water lily and avocado flowers. In contrast, spatially discrete transcriptional programs in distinct floral organs characterize the more recently derived angiosperm lineages represented by California poppy and Arabidopsis. Deep evolutionary conservation in the genetic programs of putatively homologous floral organs traces to those operating in gymnosperm reproductive cones. Female gymnosperm cones and angiosperm carpels share conserved genetic features, which may be associated with the ovule developmental program common to both organs. However, male gymnosperm cones share genetic features with both perianth (sterile attractive and protective) organs and stamens, supporting the evolutionary origin of the floral perianth from the male genetic program of seed plants.