Isolation and phylogenetic footprinting analysis of the 5'-regulatory region of the floral homeotic gene OrcPI from Orchis italica (Orchidaceae)

Homeologs of Brassica SOC1, a central regulator of flowering time, are differentially regulated due to partitioning of evolutionarily conserved transcription factor binding sites in promoters

Evolution of Brassica genome post-polyploidization reveals asymmetrical genome fractionation and copy number variation. Herein, we describe the impact of promoter divergence among SUPPRESSOR OF OVEREXPRESSION OF CONSTANS1 (SOC1) homeologs on expression and function in Brassica spp. SOC1, a regulated floral pathway integrator, is conserved as 3 redundant homeologs in diploid Brassicas. Even with high sequence identity within coding regions (92.8-100%), the spatio-temporal expression patterns of 9 SOC1 homologs in B. juncea and B. nigra indicates regulatory divergence. While LF and MF2 SOC1 homeologs are upregulated during floral transition, MF1 is barely expressed. Also, MF2 homeolog levels do not decline post-flowering, unlike LF. To investigate the underlying source of divergence, we analyzed the sequence and phylogeny of all reported (22) and isolated (21) upstream regions of Brassica SOC1. Full length upstream regions (4712-19189 bp) reveal 5 ubiquitously conserved ancestral Blocks, harboring binding sites of 18 TFs (TFBSs) characterized in Arabidopsis thaliana. The orthologs of these TFBSs are differentially conserved among Brassica SOC1 homeologs, imparting expression divergence. No crucial TFBSs are exclusively lost from LF_SOC1 promoter, while MF1_SOC1 has lost NF-Y binding site crucial for SOC1 activation by CONSTANS. MF2_SOC1 homeologs have lost important TFBSs (SEP3, AP1 and SMZ), responsible for SOC1 repression post-flowering. BjuAALF_SOC1 promoter (proximal 2 kb) shows ubiquitous reporter expression in B. juncea cv. Varuna transgenics, while BjuAAMF1_SOC1 promoter shows absence of reporter expression, validating the impact of TFBS divergence. Conservation of the original primary protein sequence is discovered in B. rapa homeologs (46) of 18 TFs. Co-regulation pattern of these TFs appeared similar for B. rapa LF and MF2 SOC1 homeologs; MF1 shows significant variation. Strong regulatory association is recorded for AP1, AP2, SEP3, FLC and CONSTANS/NF-Y, highlighting their importance in homeolog-specific SOC1 regulation. Correlation of B. juncea AP1, AP2 and FLC expression with SOC1 homeologs also complies with the TFBS differences. We thus conclude that redundant SOC1 loci contribute differentially to cumulative expression of SOC1 due to divergent selection of ancestral TFBSs.

The MADS-box genes expressed in the inflorescence of Orchis italica (Orchidaceae)

The Orchidaceae family, which is one of the most species-rich flowering plant families, includes species with highly diversified and specialized flower shapes. The aim of this study was to analyze the MADS-box genes expressed in the inflorescence of Orchis italica, a wild Mediterranean orchid species. MADS-box proteins are transcription factors involved in various plant biological processes, including flower development. In the floral tissues of O. italica, 29 MADS-box genes are expressed that are classified as both class I and II. Class I MADS-box genes include one Mβ-type gene, thereby confirming the presence of this type of MADS-box genes in orchids. The class II MIKC* gene is highly expressed in the column, which is consistent with the conserved function of the MIKC* genes in gametophyte development. In addition, homologs of the SOC, SVP, ANR1, AGL12 and OsMADS32 genes are expressed. Compared with previous knowledge on class II MIKC^C genes of O. italica involved in the ABCDE model of flower development, the number of class B and D genes has been confirmed. In addition, 4 class A (AP1/FUL) transcripts, 2 class E (SEP) transcripts, 2 new class C (AG) transcripts and 1 new AGL6 transcript have been identified. Within the AP1/FUL genes, the sequence divergence, relaxation of purifying selection and expression profiles suggest a possible functional diversification within these orchid genes. The detection of only two SEP transcripts in O. italica, in contrast with the 4 genes found in other orchids, suggests that only two SEP genes could be present in the subfamily Orchidoideae. The expression pattern of the MIKC^C genes of O. italica indicates that low levels at the boundary of the domain of a given MADS-box gene can overlap with the expression of genes belonging to a different functional A-E class in the adjacent domain, thereby following a “fading borders” model.

Histone acetylation accompanied with promoter sequences displaying differential expression profiles of B-class MADS-box genes for phalaenopsis floral morphogenesis

Five B-class MADS-box genes, including four APETALA3 (AP3)-like PeMADS2∼5 and one PISTILLATA (PI)-like PeMADS6, specify the spectacular flower morphology in orchids. The PI-like PeMADS6 ubiquitously expresses in all floral organs. The four AP3-like genes, resulted from two duplication events, express ubiquitously at floral primordia and early floral organ stages, but show distinct expression profiles at late floral organ primordia and floral bud stages. Here, we isolated the upstream sequences of PeMADS2∼6 and studied the regulatory mechanism for their distinct gene expression. Phylogenetic footprinting analysis of the 1.3-kb upstream sequences of AP3-like PeMADS2∼5 showed that their promoter regions have sufficiently diverged and contributed to their subfunctionalization. The amplified promoter sequences of PeMADS2∼6 could drive beta-glucuronidase (GUS) gene expression in all floral organs, similar to their expression at the floral primordia stage. The promoter sequence of PeMADS4, exclusively expressed in lip and column, showed a 1.6∼3-fold higher expression in lip/column than in sepal/petal. Furthermore, we noted a 4.9-fold increase in histone acetylation (H3K9K14ac) in the translation start region of PeMADS4 in lip as compared in petal. All these results suggest that the regulation via the upstream sequences and increased H3K9K14ac level may act synergistically to display distinct expression profiles of the AP3-like genes at late floral organ primordia stage for Phalaenopsis floral morphogenesis.

De novo transcriptome assembly from inflorescence of Orchis italica: analysis of coding and non-coding transcripts

The floral transcriptome of Orchis italica, a wild orchid species, was obtained using Illumina RNA-seq technology and specific de novo assembly and analysis tools. More than 100 million raw reads were processed resulting in 132,565 assembled transcripts and 86,079 unigenes with an average length of 606 bp and N50 of 956 bp. Functional annotation assigned 38,984 of the unigenes to records present in the NCBI non-redundant protein database, 32,161 of them to Gene Ontology terms, 15,775 of them to Eukaryotic Orthologous Groups (KOG) and 7,143 of them to Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways. The in silico expression analysis based on the Fragments Per Kilobase of transcript per Million mapped reads (FPKM) was confirmed by real-time RT-PCR experiments on 10 selected unigenes, which showed high and statistically significant positive correlation with the RNA-seq based expression data. The prediction of putative long non-coding RNAs was assessed using two different software packages, CPC and Portrait, resulting in 7,779 unannotated unigenes that matched the threshold values for both of the analyses. Among the predicted long non-coding RNAs, one is the homologue of TAS3, a long non-coding RNA precursor of trans-acting small interfering RNAs (ta-siRNAs). The differential expression pattern observed for the selected putative long non-coding RNAs suggests their possible functional role in different floral tissues.

The MADS and the Beauty: Genes Involved in the Development of Orchid Flowers

Since the time of Darwin, biologists have studied the origin and evolution of the Orchidaceae, one of the largest families of flowering plants. In the last two decades, the extreme diversity and specialization of floral morphology and the uncoupled rate of morphological and molecular evolution that have been observed in some orchid species have spurred interest in the study of the genes involved in flower development in this plant family. As part of the complex network of regulatory genes driving the formation of flower organs, the MADS-box represents the most studied gene family, both from functional and evolutionary perspectives. Despite the absence of a published genome for orchids, comparative genetic analyses are clarifying the functional role and the evolutionary pattern of the MADS-box genes in orchids. Various evolutionary forces act on the MADS-box genes in orchids, such as diffuse purifying selection and the relaxation of selective constraints, which sometimes reveals a heterogeneous selective pattern of the coding and non-coding regions. The emerging theory regarding the evolution of floral diversity in orchids proposes that the diversification of the orchid perianth was a consequence of duplication events and changes in the regulatory regions of the MADS-box genes, followed by sub- and neo-functionalization. This specific developmental-genetic code is termed the "orchid code."