Comparative transcriptomics among floral organs of the basal eudicot Eschscholzia californica as reference for floral evolutionary developmental studies

California poppy (Eschscholzia californica), the Papaveraceae golden girl model organism for evodevo and specialized metabolism

California poppy or golden poppy (Eschscholzia californica) is the iconic state flower of California, with native ranges from Northern California to Southwestern Mexico. It grows well as an ornamental plant in Mediterranean climates, but it might be invasive in many parts of the world. California poppy was also highly prized by Native Americans for its medicinal value, mainly due to its various specialized metabolites, especially benzylisoquinoline alkaloids (BIAs). As a member of the Ranunculales, the sister lineage of core eudicots it occupies an interesting phylogenetic position. California poppy has a short-lived life cycle but can be maintained as a perennial. It has a comparatively simple floral and vegetative morphology. Several genetic resources, including options for genetic manipulation and a draft genome sequence have been established already with many more to come. Efficient cell and tissue culture protocols are established to study secondary metabolite biosynthesis and its regulation. Here, we review the use of California poppy as a model organism for plant genetics, with particular emphasis on the evolution of development and BIA biosynthesis. In the future, California poppy may serve as a model organism to combine two formerly separated lines of research: the regulation of morphogenesis and the regulation of secondary metabolism. This can provide insights into how these two integral aspects of plant biology interact with each other.

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.

Transcriptome profiling of flower buds of male-sterile lines provides new insights into male sterility mechanism in alfalfa

BACKGROUND

The use of heterosis to produce hybrid seeds is a challenge to breeding for improved crop yield. In previous studies, we isolated a male sterile alfalfa hybrid and successfully obtained a genetically stable alfalfa male sterile line through backcrossing, henceforth named MS-4. In this study, we used RNA-seq technology to analyze the transcriptome profiles of the male sterile line (MS-4) and the male fertile line (MF) of alfalfa to elucidate the mechanism of male sterility.

RESULTS

We screened a total of 11,812 differentially expressed genes (DEGs) from both MS-4 and MF lines at three different stages of anther development. Gene Ontology (GO), and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses revealed that these DEGs are mainly involved in processes such as energy metabolism, lipid and amino acid metabolism, carbohydrate metabolism, in addition to cell synthesis and aging. The results from protein-protein interaction (PPI) network analysis showed that the ribosomal protein (MS.Gene25178) was the core gene in the network. We also found that transcriptional regulation was an influential factor in the development of anthers.

CONCLUSIONS

Our findings provide new insights into understanding of the fertility changes in the male sterile (MS-4) of alfalfa.

Comparative transcriptome between male fertile and male sterile alfalfa (Medicago varia)

Male sterility is an important factor in improving crop quality and yield through heterosis breeding. In this study, we analyzed the transcriptomes of male fertile (MF) and male sterile (MS) alfalfa flower buds using the Illumina HiSeq™ 4000 platform. A total of 54.05 million clean reads were generated and assembled into 65,777 unigenes with an average length of 874 bp. The differentially expressed genes (DEGs) between the MF and MS flowers at three stages of pollen development were identified, and there were 3832, 5678 and 5925 DEGs respectively in stages 1, 2 and 3. GO and KEGG functional enrichment analysis revealed 12, 12, 6 and 12 key branch-point genes involved in circadian rhythm, transcription factors, pollen development and flavonoid biosynthesis. Our findings provide novel insights into the mechanism of male sterility in alfalfa.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s12298-021-01026-x.

Insights into triterpene synthesis and unsaturated fatty-acid accumulation provided by chromosomal-level genome analysis of Akebia trifoliata subsp. australis

Akebia trifoliata subsp. australis is a well-known medicinal and potential woody oil plant in China. The limited genetic information available for A. trifoliata subsp. australis has hindered its exploitation. Here, a high-quality chromosome-level genome sequence of A. trifoliata subsp. australis is reported. The de novo genome assembly of 682.14 Mb was generated with a scaffold N50 of 43.11 Mb. The genome includes 25,598 protein-coding genes, and 71.18% (485.55 Mb) of the assembled sequences were identified as repetitive sequences. An ongoing massive burst of long terminal repeat (LTR) insertions, which occurred ~1.0 million years ago, has contributed a large proportion of LTRs in the genome of A. trifoliata subsp. australis. Phylogenetic analysis shows that A. trifoliata subsp. australis is closely related to Aquilegia coerulea and forms a clade with Papaver somniferum and Nelumbo nucifera, which supports the well-established hypothesis of a close relationship between basal eudicot species. The expansion of UDP-glucoronosyl and UDP-glucosyl transferase gene families and β-amyrin synthase-like genes and the exclusive contraction of terpene synthase gene families may be responsible for the abundant oleanane-type triterpenoids in A. trifoliata subsp. australis. Furthermore, the acyl-ACP desaturase gene family, including 12 stearoyl-acyl-carrier protein desaturase (SAD) genes, has expanded exclusively. A combined transcriptome and fatty-acid analysis of seeds at five developmental stages revealed that homologs of SADs, acyl-lipid desaturase omega fatty acid desaturases (FADs), and oleosins were highly expressed, consistent with the rapid increase in the content of fatty acids, especially unsaturated fatty acids. The genomic sequences of A. trifoliata subsp. australis will be a valuable resource for comparative genomic analyses and molecular breeding.

The complete chloroplast genome of Aquilegia barnebyi, a basal eudicot species

Aquilegia barnebyi, belonging to the genus Aquilegia (Ranunculaceae), is a member of basal eudicot species. In this study, we obtained the complete chloroplast (cp) genome of A. barnebyi. The genome size is 161,954 bp with a GC content of 38.98%. A total of 113 unique genes including 79 protein-coding genes, 30 tRNA genes, four rRNA genes were annotated. The large single-copy region and small single-copy region contains 91,250 bp and 17,359 bp, respectively. The inverted repeat regions are 26,671 bp in length. The phyologenetic analysis indicated that A. barnebyi had a close relationship with A. coerulea. And four species in genus Aquilegia formed a monophyletic group with high support value. The availability of A. barnebyi cp genomic resources will greatly helpful for taxonomy, phylogeny and conservation genetic studies of basal eudicot specie.

Genome-wide analysis of the MADS-box gene family in polyploid cotton (Gossypium hirsutum) and in its diploid parental species (Gossypium arboreum and Gossypium raimondii)

The MADS-box gene family encodes transcription factors that share a highly conserved domain known to bind to DNA. Members of this family control various processes of development in plants, from root formation to fruit ripening. In this work, a survey of diploid (Gossypium raimondii and Gossypium arboreum) and tetraploid (Gossypium hirsutum) cotton genomes found a total of 147, 133 and 207 MADS-box genes, respectively, distributed in the MIKC, Mα, Mβ, Mγ, and Mδ subclades. A comparative phylogenetic analysis among cotton species, Arabidopsis, poplar and grapevine MADS-box homologous genes allowed us to evaluate the evolution of each MADS-box lineage in cotton plants and identify sequences within well-established subfamilies. Chromosomal localization and phylogenetic analysis revealed that G. raimondii and G. arboreum showed a conserved evolution of the MIKC subclade and a distinct pattern of duplication events in the Mα, Mγ and Mδ subclades. Additionally, G. hirsutum showed a combination of its parental subgenomes followed by a distinct evolutionary history including gene gain and loss in each subclade. qPCR analysis revealed the expression patterns of putative homologs in the AP1, AP3, AGL6, SEP4, AGL15, AG, AGL17, TM8, SVP, SOC and TT16 subfamilies of G. hirsutum. The identification of putative cotton orthologs is discussed in the light of evolution and gene expression data from other plants. This analysis of the MADS-box genes in Gossypium species opens an avenue to understanding the origin and evolution of each gene subfamily within diploid and polyploid species and paves the way for functional studies in cotton species.

Global gene expression defines faded whorl specification of double flower domestication in Camellia

Double flowers in cultivated camellias are divergent in floral patterns which present a rich resource for demonstrating molecular modifications influenced by the human demands. Despite the key principle of ABCE model in whorl specification, the underlying mechanism of fine-tuning double flower formation remains largely unclear. Here a comprehensive comparative transcriptomics interrogation of gene expression among floral organs of wild type and “formal double” and “anemone double” is presented. Through a combination of transcriptome, small RNA and “degradome” sequencing, we studied the regulatory gene expression network underlying the double flower formation. We obtained the differentially expressed genes between whorls in wild and cultivated Camellia. We showed that the formation of double flowers tends to demolish gene expression canalization of key functions; the faded whorl specification mechanism was fundamental under the diverse patterns of double flowers. Furthermore, we identified conserved miRNA-targets regulations in the control of double flowers, and we found that miR172-AP2, miR156-SPLs were critical regulatory nodes contributing to the diversity of double flower forms. This work highlights the hierarchical patterning of global gene expression in floral development, and supports the roles of “faded ABC model” mechanism and miRNA-targets regulations underlying the double flower domestication.

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.

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

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

Perigone Lobe Transcriptome Analysis Provides Insights into Rafflesia cantleyi Flower Development

Rafflesia is a biologically enigmatic species that is very rare in occurrence and possesses an extraordinary morphology. This parasitic plant produces a gigantic flower up to one metre in diameter with no leaves, stem or roots. However, little is known about the floral biology of this species especially at the molecular level. In an effort to address this issue, we have generated and characterised the transcriptome of the Rafflesia cantleyi flower, and performed a comparison with the transcriptome of its floral bud to predict genes that are expressed and regulated during flower development. Approximately 40 million sequencing reads were generated and assembled de novo into 18,053 transcripts with an average length of 641 bp. Of these, more than 79% of the transcripts had significant matches to annotated sequences in the public protein database. A total of 11,756 and 7,891 transcripts were assigned to Gene Ontology categories and clusters of orthologous groups respectively. In addition, 6,019 transcripts could be mapped to 129 pathways in Kyoto Encyclopaedia of Genes and Genomes Pathway database. Digital abundance analysis identified 52 transcripts with very high expression in the flower transcriptome of R. cantleyi. Subsequently, analysis of differential expression between developing flower and the floral bud revealed a set of 105 transcripts with potential role in flower development. Our work presents a deep transcriptome resource analysis for the developing flower of R. cantleyi. Genes potentially involved in the growth and development of the R. cantleyi flower were identified and provide insights into biological processes that occur during flower development.

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.

Transcriptomic Analysis Using Olive Varieties and Breeding Progenies Identifies Candidate Genes Involved in Plant Architecture

Plant architecture is a critical trait in fruit crops that can significantly influence yield, pruning, planting density and harvesting. Little is known about how plant architecture is genetically determined in olive, were most of the existing varieties are traditional with an architecture poorly suited for modern growing and harvesting systems. In the present study, we have carried out microarray analysis of meristematic tissue to compare expression profiles of olive varieties displaying differences in architecture, as well as seedlings from their cross pooled on the basis of their sharing architecture-related phenotypes. The microarray used, previously developed by our group has already been applied to identify candidates genes involved in regulating juvenile to adult transition in the shoot apex of seedlings. Varieties with distinct architecture phenotypes and individuals from segregating progenies displaying opposite architecture features were used to link phenotype to expression. Here, we identify 2252 differentially expressed genes (DEGs) associated to differences in plant architecture. Microarray results were validated by quantitative RT-PCR carried out on genes with functional annotation likely related to plant architecture. Twelve of these genes were further analyzed in individual seedlings of the corresponding pool. We also examined Arabidopsis mutants in putative orthologs of these targeted candidate genes, finding altered architecture for most of them. This supports a functional conservation between species and potential biological relevance of the candidate genes identified. This study is the first to identify genes associated to plant architecture in olive, and the results obtained could be of great help in future programs aimed at selecting phenotypes adapted to modern cultivation practices in this species.

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

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

The 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.

A de novo floral transcriptome reveals clues into Phalaenopsis orchid flower development

Phalaenopsis has a zygomorphic floral structure, including three outer tepals, two lateral inner tepals and a highly modified inner median tepal called labellum or lip; however, the regulation of its organ development remains unelucidated. We generated RNA-seq reads with the Illumina platform for floral organs of the Phalaenopsis wild-type and peloric mutant with a lip-like petal. A total of 43,552 contigs were obtained after de novo assembly. We used differentially expressed gene profiling to compare the transcriptional changes in floral organs for both the wild-type and peloric mutant. Pair-wise comparison of sepals, petals and labellum between peloric mutant and its wild-type revealed 1,838, 758 and 1,147 contigs, respectively, with significant differential expression. PhAGL6a (CUFF.17763), PhAGL6b (CUFF.17763.1), PhMADS1 (CUFF.36625.1), PhMADS4 (CUFF.25909) and PhMADS5 (CUFF.39479.1) were significantly upregulated in the lip-like petal of the peloric mutant. We used real-time PCR analysis of lip-like petals, lip-like sepals and the big lip of peloric mutants to confirm the five genes' expression patterns. PhAGL6a, PhAGL6b and PhMADS4 were strongly expressed in the labellum and significantly upregulated in lip-like petals and lip-like sepals of peloric-mutant flowers. In addition, PhAGL6b was significantly downregulated in the labellum of the big lip mutant, with no change in expression of PhAGL6a. We provide a comprehensive transcript profile and functional analysis of Phalaenopsis floral organs. PhAGL6a PhAGL6b, and PhMADS4 might play crucial roles in the development of the labellum in Phalaenopsis. Our study provides new insights into how the orchid labellum differs and why the petal or sepal converts to a labellum in Phalaenopsis floral mutants.

Comparative transcriptome analysis of the petal degeneration mutant pdm in Chinese cabbage (Brassica campestris ssp. pekinensis) using RNA-Seq

Flowering, which plays a crucial role in the growth and development of flowering plants, is a crucial point from vegetative growth to reproductive growth. The goal of this study was to examine the differences between the transcriptomes of the Chinese cabbage mutant pdm and the corresponding wild-type line 'FT'. We performed transcriptome analysis on mRNA isolated from flower buds of pdm and 'FT' using Illumina RNA sequencing (RNA-Seq) data. A total of 117 differentially expressed genes (DEGs) were detected. Among the DEGs, we identified a number of genes involved in floral development and flowering, including an F-box protein gene, EARLY FLOWERING 4 (ELF4), and transcription factors BIGPETAL (BPE) and MYB21 (v-myb avian myeloblastosis viral oncogene homolog); differential expression of these genes could potentially explain the difference in the flowers between pdm and 'FT'. In addition, the expression patterns of 20 DEGs, including 12 floral development and flowering-related genes and eight randomly selected genes, were validated by qRT-PCR, and the results were highly concordant with the RNA-Seq results. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes pathway enrichment analyses were performed to better understand the functions of these DEGs. We also identified a large number of single nucleotide polymorphism and insertion/deletion markers, which will be a rich resource for future marker development and breeding research in Chinese cabbage. Also, our analysis revealed numerous novel transcripts and alternative splicing events. The transcriptome analysis provides valuable information for furthering our understanding of the molecular mechanisms that regulate the flowering process, and establishes a solid foundation for future genetic and functional genomic studies in Chinese cabbage.

De novo transcriptome analysis of petal senescence in Gardenia jasminoides Ellis

BACKGROUND

The petal senescence of ethylene insensitive species has not been investigated thoroughly while little is known about the temporal and tissue specific expression patterns of transcription factors (TFs) in this developmental process. Even less is known on flower senescence of the ornamental pot plant Gardenia jasminoides, a non climacteric flower with significant commercial value.

RESULTS

We initiated a de novo transcriptome study to investigate the petal senescence in four developmental stages of cut gardenia flowers considering that the visible symptoms of senescence appear within 4 days of flower opening. De novo assembly of transcriptome sequencing resulted in 102,263 contigs with mean length of 360 nucleotides that generated 57,503 unigenes. These were further clustered into 20,970 clusters and 36,533 singletons. The comparison of the consecutive developmental stages resulted in 180 common, differentially expressed unigenes. A large number of Simple Sequence Repeats were also identified comprising a large number of dinucleotides and trinucleotides. The prevailing families of differentially expressed TFs comprise the AP2/EREBP, WRKY and the bHLH. There are 81 differentially expressed TFs when the symptoms of flower senescence become visible with the most prevailing being the WRKY family with 19 unigenes. No other WRKY TFs had been identified up to now in petal senescence of ethylene insensitive species. A large number of differentially expressed genes were identified at the initiation of visible symptoms of senescence compared to the open flower stage indicating a significant shift in the expression profiles which might be coordinated by up-regulated and/or down-regulated TFs. The expression of 16 genes that belong to the TF families of WRKY, bHLH and the ethylene sensing pathway was validated using qRT--PCR.

CONCLUSION

This de novo transcriptome analysis resulted in the identification of TFs with specific temporal expression patterns such as two WRKYs and one bHLH, which might play the role of senescence progression regulators. Further research is required to investigate their role in gardenia flowers in order to develop tools to delay petal senescence.

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.

Dissecting organ-specific transcriptomes through RNA-sequencing

BACKGROUND

Organ-specific gene expression contains rich information about in vivo biological processes. This kind of information, previously gathered through microarray profiling, has been proven fruitful to the understanding of specific mutants, regulatory events, signaling, and development. With the advent of the next-generation sequencing (NGS) of RNAs, more quantitative and detailed information of gene expressions than previously available can now be collected for each organ or organ developmental stages. The combination of an object-oriented experimental design and an efficient treatment of the high volume information generated through a NGS platform may offer a powerful tool for inferring previously intractable developmental processes.

RESULTS

We collected transcriptomic data over a Solexa/Illumina platform on samples of Ipomoea leaf, sepal, and petals (at three developmental stages), and presented a method for analyzing transcriptomic variations within and between organs. We demonstrated that in vivo signals of transcriptomes can be retrieved de novo through the NGS techniques, proper data handling, bioinformatic tools, and the current understanding of molecular networks. We found that numbers of transcribed genes from both nuclear and chloroplast genomes decreased by the same order of leaf → sepal →petal. Petal resembled leaf in cell division patterns and abundance level of commonly expressed organelle genes. Its chloroplast transcripts constituted a subset of those in leaf. Moreover, reconstructions of multiple metabolic networks for each organ enabled inferences of substance flow, providing transcript evidence for the path of sucrose in leaf to anthocyanin synthesis in petal.

CONCLUSION

Our results attest that developmental transcriptomes are highly informative for exploring connections between morphological traits and the associated molecular networks. Significant hypotheses have been developed, including that the petal is a derived organ of leaf and that its color can be modified by fluctuations of substance flow within the associated metabolic networks among organs.

Separate-channel analysis of two-channel microarrays: recovering inter-spot information

BACKGROUND

Two-channel (or two-color) microarrays are cost-effective platforms for comparative analysis of gene expression. They are traditionally analysed in terms of the log-ratios (M-values) of the two channel intensities at each spot, but this analysis does not use all the information available in the separate channel observations. Mixed models have been proposed to analyse intensities from the two channels as separate observations, but such models can be complex to use and the gain in efficiency over the log-ratio analysis is difficult to quantify. Mixed models yield test statistics for the null distributions can be specified only approximately, and some approaches do not borrow strength between genes.

RESULTS

This article reformulates the mixed model to clarify the relationship with the traditional log-ratio analysis, to facilitate information borrowing between genes, and to obtain an exact distributional theory for the resulting test statistics. The mixed model is transformed to operate on the M-values and A-values (average log-expression for each spot) instead of on the log-expression values. The log-ratio analysis is shown to ignore information contained in the A-values. The relative efficiency of the log-ratio analysis is shown to depend on the size of the intraspot correlation. A new separate channel analysis method is proposed that assumes a constant intra-spot correlation coefficient across all genes. This approach permits the mixed model to be transformed into an ordinary linear model, allowing the data analysis to use a well-understood empirical Bayes analysis pipeline for linear modeling of microarray data. This yields statistically powerful test statistics that have an exact distributional theory. The log-ratio, mixed model and common correlation methods are compared using three case studies. The results show that separate channel analyses that borrow strength between genes are more powerful than log-ratio analyses. The common correlation analysis is the most powerful of all.

CONCLUSIONS

The common correlation method proposed in this article for separate-channel analysis of two-channel microarray data is no more difficult to apply in practice than the traditional log-ratio analysis. It provides an intuitive and powerful means to conduct analyses and make comparisons that might otherwise not be possible.

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.

Virus-induced gene silencing of the alkaloid-producing basal eudicot model plant Eschscholzia californica (California Poppy)

Eschscholzia californica (California poppy), a member of the basal eudicot family of the Papaveraceae, is an important species to study alkaloid biosynthesis and the effect of alkaloids on plant metabolism. More recently, it has also been developed as a model system to study the evolution of plant morphogenesis. While progress has been made towards establishing methods for generating genetically modified cell culture lines, transcriptome data and gene expression analysis, the stable transformation and subsequent regeneration of transgenic plants has proven extremely time consuming and difficult. Here, we describe in detail a method to transiently down regulate expression of a target gene by virus-induced gene silencing (VIGS) and the subsequent analysis of the VIGS treated plants. VIGS in E. californica allows for the study of gene function within 2 to 3 weeks after inoculation, and the method proves very efficient, enabling the rapid analysis of gene functions.

Developmental and evolutionary diversity of plant MADS-domain factors: insights from recent studies

Members of the MADS-box transcription factor family play essential roles in almost every developmental process in plants. Many MADS-box genes have conserved functions across the flowering plants, but some have acquired novel functions in specific species during evolution. The analyses of MADS-domain protein interactions and target genes have provided new insights into their molecular functions. Here, we review recent findings on MADS-box gene functions in Arabidopsis and discuss the evolutionary history and functional diversification of this gene family in plants. We also discuss possible mechanisms of action of MADS-domain proteins based on their interactions with chromatin-associated factors and other transcriptional regulators.

De novo sequencing and characterization of the floral transcriptome of Dendrocalamus latiflorus (Poaceae: Bambusoideae)

Background

Transcriptome sequencing can be used to determine gene sequences and transcript abundance in non-model species, and the advent of next-generation sequencing (NGS) technologies has greatly decreased the cost and time required for this process. Transcriptome data are especially desirable in bamboo species, as certain members constitute an economically and culturally important group of mostly semelparous plants with remarkable flowering features, yet little bamboo genomic research has been performed. Here we present, for the first time, extensive sequence and transcript abundance data for the floral transcriptome of a key bamboo species, Dendrocalamus latiflorus, obtained using the Illumina GAII sequencing platform. Our further goal was to identify patterns of gene expression during bamboo flower development.

Results

Approximately 96 million sequencing reads were generated and assembled de novo, yielding 146,395 high quality unigenes with an average length of 461 bp. Of these, 80,418 were identified as putative homologs of annotated sequences in the public protein databases, of which 290 were associated with the floral transition and 47 were related to flower development. Digital abundance analysis identified 26,529 transcripts differentially enriched between two developmental stages, young flower buds and older developing flowers. Unigenes found at each stage were categorized according to their putative functional categories. These sequence and putative function data comprise a resource for future investigation of the floral transition and flower development in bamboo species.

Conclusions

Our results present the first broad survey of a bamboo floral transcriptome. Although it will be necessary to validate the functions carried out by these genes, these results represent a starting point for future functional research on D. latiflorus and related species.

Plant 'evo-devo' goes genomic: from candidate genes to regulatory networks

Plant development gives rise to a staggering complexity of morphological structures with different shapes, colors, and functions. Understanding the evolution of control mechanisms that underlie developmental processes provides insights into causes of morphological diversity and, therefore, is of great interest to biologists. New genomic resources and techniques enable biologists to assess for the first time the evolution of developmental regulatory networks at a global scale. Here, we address the question of how comparative regulatory genomics can be used to reveal the evolutionary dynamics of control networks linked to morphological evolution in plants.

Transcript and metabolite profiling in cell cultures of 18 plant species that produce benzylisoquinoline alkaloids

Benzylisoquinoline alkaloids (BIAs) are a large and diverse group of ~2500 specialized metabolites found predominantly in plants of the order Ranunculales. Research focused on BIA metabolism in a restricted number of plant species has identified many enzymes and cognate genes involved in the biosynthesis of compounds such as morphine, sanguinarine and berberine. However, the formation of most BIAs remains uncharacterized at the molecular biochemical level. Herein a compendium of sequence- and metabolite-profiling resources from 18 species of BIA-accumulating cell cultures was established, representing four related plant families. Our integrated approach consisted of the construction of EST libraries each containing approximately 3500 unigenes per species for a total of 58,787 unigenes. The EST libraries were manually triaged using known BIA-biosynthetic genes as queries to identify putative homologs with similar or potentially different functions. Sequence resources were analyzed in the context of the targeted metabolite profiles obtained for each cell culture using electrospray-ionization and collision-induced dissociation mass spectrometry. Fragmentation analysis was used for the identification or structural characterization coupled with the relative quantification of 72 BIAs, which establishes a key resource for future work on alkaloid biosynthesis. The metabolite profile obtained for each species provides a rational basis for the prediction of enzyme function in BIA metabolism. The metabolic frameworks assembled through the integration of transcript and metabolite profiles allow a comparison of BIA metabolism across several plant species and families. Taken together, these data represent an important tool for the discovery of BIA biosynthetic genes.

Virus-induced gene silencing (VIGS) in Cysticapnos vesicaria, a zygomorphic-flowered Papaveraceae (Ranunculales, basal eudicots)

BACKGROUND AND AIMS

Studies of evolutionary diversification in the basal eudicot family Papaveraceae, such as the transition from actinomorphy to zygomorphy, are hampered by the lack of comparative functional studies. So far, gene silencing methods are only available in the actinomorphic species Eschscholzia californica and Papaver somniferum. This study addresses the amenability of Cysticapnos vesicaria, a derived fumitory with zygomorphic flowers, to virus-induced gene silencing (VIGS), and describes vegetative and reproductive traits in this species.

METHODS

VIGS-mediated downregulation of the C. vesicaria PHYTOENE DESATURASE gene (CvPDS) and of the FLORICAULA gene CvFLO was carried out using Agrobacterium tumefaciens transfer of Tobacco rattle virus (TRV)-based vectors. Wild-type and vector-treated plants were characterized using reverse transcription-PCR (RT-PCR), in situ hybridization, and macroscopic and scanning electron microscopic imaging.

KEY RESULTS

Cysticapnos vesicaria germinates rapidly, can be grown at high density, has a short life cycle and is self-compatible. Inoculation of C. vesicaria with a CvPDS-VIGS vector resulted in strong photobleaching of green parts and reduction of endogenous CvPDS transcript levels. Gene silencing persisted during inflorescence development until fruit set. Inoculation of plants with CvFLO-VIGS affected floral phyllotaxis, symmetry and floral organ identities.

CONCLUSIONS

The high penetrance, severity and stability of pTRV-mediated silencing, including the induction of meristem-related phenotypes, make C. vesicaria a very promising new focus species for evolutionary-developmental (evo-devo) studies in the Papaveraceae. This now enables comparative studies of flower symmetry, inflorescence determinacy and other traits that diversified in the Papaveraceae.

AMS-dependent and independent regulation of anther transcriptome and comparison with those affected by other Arabidopsis anther genes

BACKGROUND

In flowering plants, the development of male reproductive organs is controlled precisely to achieve successful fertilization and reproduction. Despite the increasing knowledge of genes that contribute to anther development, the regulatory mechanisms controlling this process are still unclear.

RESULTS

In this study, we analyzed the transcriptome profiles of early anthers of sterile mutants aborted microspores (ams) and found that 1,368 genes were differentially expressed in ams compared to wild type anthers, affecting metabolism, transportation, ubiquitination and stress response. Moreover, the lack of significant enrichment of potential AMS binding sites (E-box) in the promoters of differentially expressed genes suggests both direct and indirect regulation for AMS-dependent regulation of anther transcriptome involving other transcription factors. Combining ams transcriptome profiles with those of two other sterile mutants, spl/nzz and ems1/exs, expression of 3,058 genes were altered in at least one mutant. Our investigation of expression patterns of major transcription factor families, such as bHLH, MYB and MADS, suggested that some closely related homologs of known anther developmental genes might also have similar functions. Additionally, comparison of expression levels of genes in different organs suggested that anther-preferential genes could play important roles in anther development.

CONCLUSION

Analysis of ams anther transcriptome and its comparison with those of spl/nzz and ems1/exs anthers uncovered overlapping and distinct sets of regulated genes, including those encoding transcription factors and other proteins. These results support an expanded regulatory network for early anther development, providing a series of hypotheses for future experimentation.

AMS-dependent and independent regulation of anther transcriptome and comparison with those affected by other Arabidopsis anther genes

BACKGROUND

In flowering plants, the development of male reproductive organs is controlled precisely to achieve successful fertilization and reproduction. Despite the increasing knowledge of genes that contribute to anther development, the regulatory mechanisms controlling this process are still unclear.

RESULTS

In this study, we analyzed the transcriptome profiles of early anthers of sterile mutants aborted microspores (ams) and found that 1,368 genes were differentially expressed in ams compared to wild type anthers, affecting metabolism, transportation, ubiquitination and stress response. Moreover, the lack of significant enrichment of potential AMS binding sites (E-box) in the promoters of differentially expressed genes suggests both direct and indirect regulation for AMS-dependent regulation of anther transcriptome involving other transcription factors. Combining ams transcriptome profiles with those of two other sterile mutants, spl/nzz and ems1/exs, expression of 3,058 genes were altered in at least one mutant. Our investigation of expression patterns of major transcription factor families, such as bHLH, MYB and MADS, suggested that some closely related homologs of known anther developmental genes might also have similar functions. Additionally, comparison of expression levels of genes in different organs suggested that anther-preferential genes could play important roles in anther development.

CONCLUSION

Analysis of ams anther transcriptome and its comparison with those of spl/nzz and ems1/exs anthers uncovered overlapping and distinct sets of regulated genes, including those encoding transcription factors and other proteins. These results support an expanded regulatory network for early anther development, providing a series of hypotheses for future experimentation.

Evolution of the YABBY gene family with emphasis on the basal eudicot Eschscholzia californica (Papaveraceae)

YABBY genes are seed plant-specific transcriptional regulators that are involved in diverse aspects of leaf, shoot and flower development. A series of duplications gave rise to five gene groups found throughout flowering plants. In Arabidopsis and other species, expression of two gene groups, CRABS CLAW and INNER NO OUTER, is restricted to floral organs. In contrast, members of the FILAMENTOUS FLOWER, YABBY2 and YABBY5 gene groups are also expressed in leaves and have been termed 'vegetative YABBYs'. How the five paralogue groups evolved and how their expression and function diversified have remained largely unresolved, precluding a reconstruction of the natural history of this gene family. Here, we report new genes from Eschscholzia californica (Ranunculales, Papaveraceae) that we use together with currently available database sequences in a comprehensive phylogenetic re-evaluation of the YABBY gene family. Multilayered Bayesian analysis covering seed plants allowed us to locate Eschscholzia YABBY sequences within the gene family phylogeny. We established that vegetative YABBYs do not form a monophyletic clade, and that CRABS CLAW and FILAMENTOUS FLOWER arose from a common ancestor gene. INNER NO OUTER genes are sister to that ancestral gene. We identified several conserved motifs outside of known amino acid domains that define all five angiosperm YABBY gene clades. Further, we inferred the evolution of gene expression and provide evidence for release of purifying constraint in certain branches of the gene family tree. Finally, we report expression patterns for five Eschscholzia YABBY genes consistent with functional conservation between early-diverged and core eudicots.

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.