These two authors contributed equally to this work. SUMMARYThe establishment of meristematic domains with different transcriptional activity is essential for many developmental processes. The asymmetry of the Antirrhinum majus flower is established by transcription factors with an asymmetric pattern of activity. To understand how this asymmetrical pattern is established, we studied the molecular mechanism through which the dorsal MYB protein RADIALIS (RAD) restricts the activity of the MYB transcription factor DIVARICATA (DIV) to the ventral region of the flower meristem. We show that RAD competes with DIV for binding with other MYB-like proteins, termed DRIF1 and DRIF2 (DIVand-RAD-interacting-factors). DRIF1 and DIV interact to form a protein complex that binds to the DIV-DNA consensus region, suggesting that the DRIFs act as co-regulators of DIV transcriptional activity. In the presence of RAD, the interaction between DRIF1 and DIV bound to DNA is disrupted. Moreover, the DRIFs are sequestered in the cytoplasm by RAD, thus, preventing or reducing the formation of DRIF-DIV heterodimers in the nuclei. Our results suggest that in the dorsal region of the Antirrhinum flower meristem the dorsal protein RAD antagonises the activity of the ventral identity protein DIV in a subcellular competition for a DRIF protein promoting the establishment of the asymmetric pattern of gene activity in the Antirrhinum flower.

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.

BackgroundCork oak (Quercus suber) is one of the rare trees with the ability to produce cork, a material widely used to make wine bottle stoppers, flooring and insulation materials, among many other uses. The molecular mechanisms of cork formation are still poorly understood, in great part due to the difficulty in studying a species with a long life-cycle and for which there is scarce molecular/genomic information. Cork oak forests are of great ecological importance and represent a major economic and social resource in Southern Europe and Northern Africa. However, global warming is threatening the cork oak forests by imposing thermal, hydric and many types of novel biotic stresses. Despite the economic and social value of the Q. suber species, few genomic resources have been developed, useful for biotechnological applications and improved forest management.ResultsWe generated in excess of 7 million sequence reads, by pyrosequencing 21 normalized cDNA libraries derived from multiple Q. suber tissues and organs, developmental stages and physiological conditions. We deployed a stringent sequence processing and assembly pipeline that resulted in the identification of ~159,000 unigenes. These were annotated according to their similarity to known plant genes, to known Interpro domains, GO classes and E.C. numbers. The phylogenetic extent of this ESTs set was investigated, and we found that cork oak revealed a significant new gene space that is not covered by other model species or EST sequencing projects. The raw data, as well as the full annotated assembly, are now available to the community in a dedicated web portal at http://www.corkoakdb.org.ConclusionsThis genomic resource represents the first trancriptome study in a cork producing species. It can be explored to develop new tools and approaches to understand stress responses and developmental processes in forest trees, as well as the molecular cascades underlying cork differentiation and disease response.

Four putative AGP genes were identified that are preferentially expressed in the male flower compared with the female flower. The putative Arabidopsis thaliana orthologues of these genes are associated with preferential expression in pollen, suggesting that the AGPs probably play a significant role in cork oak reproduction.

Monoecious species provide an excellent system to study the specific determinants that underlie male and female flower development. Quercus suber is a monoecious species with unisexual flowers at inception. Despite the overall importance of this and other tree species with a similar reproductive habit, little is known regarding the mechanisms involved in the development of their male and female flowers. Here, we have characterised members of the ABCDE MADS-box gene family of Q. suber. The temporal expression of these genes was found to be sex-biased. The B-class genes, in particular, are predominantly, or exclusively (in the case of QsPISTILLATA), expressed in the male flowers. Functional analysis in Arabidopsis suggests that the B-class genes have their function conserved. The identification of sex-biased gene expression plus the identification of unusual protein-protein interactions suggest that the floral organ identity of Q. suber may be under control of specific changes in the dynamics of the ABCDE model. This study constitutes a major step towards the characterisation of the mechanisms involved in reproductive organ identity in a monoecious tree with a potential contribution towards the knowledge of conserved developmental mechanisms in other species with a similar sex habit.

The establishment of new interactions between transcriptional regulators increases the regulatory diversity that drives phenotypic novelty. To understand how such interactions evolve, we have studied a regulatory module (DDR) composed by three MYB-like proteins: DIVARICATA (DIV), RADIALIS (RAD), and DIV-and-RAD-Interacting Factor (DRIF). The DIV and DRIF proteins form a transcriptional complex that is disrupted in the presence of RAD, a small interfering peptide, due to the formation of RAD–DRIF dimers. This dynamic interaction result in a molecular switch mechanism responsible for the control of distinct developmental processes in plants. Here, we have determined how the DDR regulatory module was established by analyzing the origin and evolution of the DIV, DRIF, and RAD protein families and the evolutionary history of their interactions. We show that duplications of a pre-existing MYB domain originated the DIV and DRIF protein families in the ancestral lineage of green algae, and, later, the RAD family in seed plants. Intraspecies interactions between the MYB domains of DIV and DRIF proteins are detected in green algae, whereas the earliest evidence of an interaction between DRIF and RAD proteins occurs in the gymnosperms, coincident with the establishment of the RAD family. Therefore, the DDR module evolved in a stepwise progression with the DIV–DRIF transcription complex evolving prior to the antagonistic RAD–DRIF interaction that established the molecular switch mechanism. Our results suggest that the successive rearrangement and divergence of a single protein domain can be an effective evolutionary mechanism driving new protein interactions and the establishment of novel regulatory modules.

The MYB transcription factors DIVARICATA (DIV), DIV-and-RAD-Interacting-Factor (DRIF), and the small interfering peptide RADIALIS (RAD) can interact, forming a regulatory module that controls different plant developmental processes. In the snapdragon Antirrhinum majus, this module, together with the TCP transcription factor CYCLOIDEA (CYC), is responsible for the establishment of floral dorsoventral asymmetry. The spatial gene expression pattern of the OitDIV, OitDRIF, and OitRAD homologs of Orchis italica, an orchid with zygomorphic flowers, has suggested a possible conserved role of these genes in bilateral symmetry of the orchid flower. Here, we have identified four DRIF genes of orchids and have reconstructed their genomic organization and evolution. In addition, we found snapdragon transcriptional cis-regulatory elements of DIV and RAD loci generally conserved within the corresponding orchid orthologues. We have tested the biochemical interactions among OitDIV, OitDRIF1, and OitRAD of O. italica, showing that OitDRIF1 can interact both with OitDIV and OitRAD, whereas OitDIV and OitRAD do not directly interact, as in A. majus. The analysis of the quantitative expression profile of these MYB genes revealed that in zygomorphic orchid flowers, the DIV, DRIF1, and RAD transcripts are present at higher levels in the lip than in lateral inner tepals, whereas in peloric orchid flowers they show similar expression levels. These results indicate that MYB transcription factors could have a role in shaping zygomorphy of the orchid flower, potentially enriching the underlying orchid developmental code.

The understanding of the molecular mechanisms responsible for the making of a unisexual flower has been a long-standing quest in plant biology. Plants with male and female flowers can be divided mainly into two categories: dioecious and monoecious, and both sexual systems co-exist in nature in ca of 10% of the angiosperms. The establishment of male and female traits has been extensively described in a hermaphroditic flower and requires the interplay of networks, directly and indirectly related to the floral organ identity genes including hormonal regulators, transcription factors, microRNAs, and chromatin-modifying proteins. Recent transcriptomic studies have been uncovering the molecular processes underlying the establishment of unisexual flowers and there are many parallelisms between monoecious, dioecious, and hermaphroditic individuals. Here, we review the paper entitled “Comparative transcriptomic analysis of male and female flowers of monoecious Quercus suber” published in 2014 in the Frontiers of Plant Science (volume 5 |Article 599) and discussed it in the context of recent studies with other dioecious and monoecious plants that utilized high-throughput platforms to obtain transcriptomic profiles of male and female unisexual flowers. In some unisexual flowers, the developmental programs that control organ initiation fail and male or female organs do not form, whereas in other species, organ initiation and development occur but they abort or arrest during different species-specific stages of differentiation. Therefore, a direct comparison of the pathways responsible for the establishment of unisexual flowers in different species are likely to reveal conserved modules of gene regulatory hubs involved in stamen or carpel development, as well as differences that reflect the different stages of development in which male and/or female organ arrest or loss-of-function occurs.