Global transcriptome analysis and identification of a CONSTANS-like gene family in the orchid Erycina pusilla

Evolution of major flowering pathway integrators in Orchidaceae

The Orchidaceae is a mega-diverse plant family with ca. 29,000 species with a large variety of life forms that can colonize transitory habitats. Despite this diversity, little is known about their flowering integrators in response to specific environmental factors. During the reproductive transition in flowering plants a vegetative apical meristem (SAM) transforms into an inflorescence meristem (IM) that forms bracts and flowers. In model grasses, like rice, a flowering genetic regulatory network (FGRN) controlling reproductive transitions has been identified, but little is known in the Orchidaceae. In order to analyze the players of the FRGN in orchids, we performed comprehensive phylogenetic analyses of CONSTANS-like/CONSTANS-like 4 (COL/COL4), FLOWERING LOCUS D (FD), FLOWERING LOCUS C/FRUITFULL (FLC/FUL) and SUPRESSOR OF OVEREXPRESSION OF CONSTANS 1 (SOC1) gene lineages. In addition to PEBP and AGL24/SVP genes previously analyzed, here we identify an increase of orchid homologs belonging to COL4, and FUL gene lineages in comparison with other monocots, including grasses, due to orchid-specific gene lineage duplications. Contrariwise, local duplications in Orchidaceae are less frequent in the COL, FD and SOC1 gene lineages, which points to a retention of key functions under strong purifying selection in essential signaling factors. We also identified changes in the protein sequences after such duplications, variation in the evolutionary rates of resulting paralogous clades and targeted expression of isolated homologs in different orchids. Interestingly, vernalization-response genes like VERNALIZATION1 (VRN1) and FLOWERING LOCUS C (FLC) are completely lacking in orchids, or alternatively are reduced in number, as is the case of VERNALIZATION2/GHD7 (VRN2). Our findings point to non-canonical factors sensing temperature changes in orchids during reproductive transition. Expression data of key factors gathered from Elleanthus auratiacus, a terrestrial orchid in high Andean mountains allow us to characterize which copies are actually active during flowering. Altogether, our data lays down a comprehensive framework to assess gene function of a restricted number of homologs identified more likely playing key roles during the flowering transition, and the changes of the FGRN in neotropical orchids in comparison with temperate grasses.

Combined Transcriptomic and Metabolomic Approach Revealed a Relationship between Light Control, Photoprotective Pigments, and Lipid Biosynthesis in Olives

Olive possesses excellent nutritional and economic values for its main healthy products. Among them, a high content of antioxidant compounds, balanced during the ripening process, are produced under genetic and environmental control, resulting in high variability among cultivars. The genes involved in these complex pathways are mainly known, but despite many studies which indicated the key role of light quality and quantity for the synthesis of many metabolites in plants, limited information on these topics is available in olive. We carried out a targeted gene expression profiling in three olive cultivars, Cellina di Nardò, Ruveia, and Salella, which were selected for their contrasting oleic acid and phenolic content. The -omics combined approach revealed a direct correlation between a higher expression of the main flavonoid genes and the high content of these metabolites in 'Cellina di Nardò'. Furthermore, it confirmed the key role of FAD2-2 in the linoleic acid biosynthesis. More interestingly, in all the comparisons, a co-regulation of genes involved in photoperception and circadian clock machinery suggests a key role of light in orchestrating the regulation of these pathways in olive. Therefore, the identified genes in our analyses might represent a useful tool to support olive breeding, although further investigations are needed.

Transcriptome mining of hormonal and floral integrators in the leafless flowers of three cymbidium orchids

Flowering is the most studied ornamental trait in orchids where long vegetative phase may span up to three years. Cymbidium orchids produce beautiful flowers with astonishing shapes and pleasant scent. However, an unusually long vegetative phase is a major drawback to their ornamental value. We observed that under certain culture conditions, three cymbidium species (Cymbidium ensifolium, C. goeringii and C. sinense) skipped vegetative growth phase and directly flowered within six months, that could be a breakthrough for future orchids with limited vegetative growth. Hormonal and floral regulators could be the key factors arresting vegetative phase. Therefore, transcriptomic analyses were performed for leafless flowers and normal vegetative leaves to ascertain differentially expressed genes (DEGs) related to hormones (auxin, cytokinin, gibberellin, abscisic acid and ethylene), floral integrators and MADS-box genes. A significant difference of cytokinin and floral regulators was observed among three species as compared to other hormones. The MADS-box genes were significantly expressed in the leafless flowers of C. sinense as compared to other species. Among the key floral regulators, CONSTANS and AGAMOUS-like genes showed the most differential expression in the leafless flowers as compared to leaves where the expression was negligible. However, CONSTANS also showed downregulation. Auxin efflux carriers were mainly downregulated in the leafless flowers of C. ensifolium and C. sinense, while they were upregulated in C. goeringii. Moreover, gibberellin and cytokinin genes were also downregulated in C. ensifolium and C. sinense flowers, while they were upregulated in C. goeringii, suggesting that species may vary in their responses. The data mining thus, outsources the valuable information to direct future research on orchids at industrial levels.

Genetic insights into the regulatory pathways for continuous flowering in a unique orchid Arundina graminifolia

BACKGROUND

Manipulation of flowering time and frequency of blooming is key to enhancing the ornamental value of orchids. Arundina graminifolia is a unique orchid that flowers year round, although the molecular basis of this flowering pattern remains poorly understood.

RESULTS

We compared the A. graminifolia transcriptome across tissue types and floral developmental stages to elucidate important genetic regulators of flowering and hormones. Clustering analyses identified modules specific to floral transition and floral morphogenesis, providing a set of candidate regulators for the floral initiation and timing. Among candidate floral homeotic genes, the expression of two FT genes was positively correlated with flower development. Assessment of the endogenous hormone levels and qRT-PCR analysis of 32 pathway-responsive genes supported a role for the regulatory networks in floral bud control in A. graminifolia. Moreover, WGCNA showed that flowering control can be delineated by modules of coexpressed genes; especially, MEgreen presented group of genes specific to flowering.

CONCLUSIONS

Candidate gene selection coupled with hormonal regulators brings a robust source to understand the intricate molecular regulation of flowering in precious orchids.

Morphological and Molecular Characterization of Orchid Fruit Development

Efficient seed dispersal in flowering plants is enabled by the development of fruits, which can be either dehiscent or indehiscent. Dehiscent fruits open at maturity to shatter the seeds, while indehiscent fruits do not open and the seeds are dispersed in various ways. The diversity in fruit morphology and seed shattering mechanisms is enormous within the flowering plants. How these different fruit types develop and which molecular networks are driving fruit diversification is still largely unknown, despite progress in eudicot model species. The orchid family, known for its astonishing floral diversity, displays a huge variation in fruit dehiscence types, which have been poorly investigated. We undertook a combined approach to understand fruit morphology and dehiscence in different orchid species to get more insight into the molecular network that underlies orchid fruit development. We describe fruit development in detail for the epiphytic orchid species Erycina pusilla and compare it to two terrestrial orchid species: Cynorkis fastigiata and Epipactis helleborine. Our anatomical analysis provides further evidence for the split carpel model, which explains the presence of three fertile and three sterile valves in most orchid species. Interesting differences were observed in the lignification patterns of the dehiscence zones. While C. fastigiata and E. helleborine develop a lignified layer at the valve boundaries, E. pusilla fruits did not lignify at these boundaries, but formed a cuticle-like layer instead. We characterized orthologs of fruit-associated MADS-domain transcription factors and of the Arabidopsis dehiscence-related genes INDEHISCENT (IND)/HECATE 3 (HEC3), REPLUMLESS (RPL) and SPATULA (SPT)/ALCATRAZ (ALC) in E. pusilla, and found that the key players of the eudicot fruit regulatory network appear well-conserved in monocots. Protein-protein interaction studies revealed that MADS-domain complexes comprised of FRUITFULL (FUL), SEPALLATA (SEP) and AGAMOUS (AG) /SHATTERPROOF (SHP) orthologs can also be formed in E. pusilla, and that the expression of HEC3, RPL, and SPT can be associated with dehiscence zone development similar to Arabidopsis. Our expression analysis also indicates differences, however, which may underlie fruit divergence.

Photoperiod response-related gene SiCOL1 contributes to flowering in sesame

BACKGROUND

Sesame is a major oilseed crop which is widely cultivated all around the world. Flowering, the timing of transition from vegetative to reproductive growth, is one of the most important events in the life cycle of sesame. Sesame is a typical short-day (SD) plant and its flowering is largely affected by photoperiod. However, the flowering mechanism in sesame at the molecular level is still not very clear. Previous studies showed that the CONSTANS (CO) gene is the crucial photoperiod response gene which plays a center role in duration of the plant vegetative growth.

RESULTS

In this study, the CO-like (COL) genes were identified and characterized in the sesame genome. Two homologs of the CO gene in the SiCOLs, SiCOL1 and SiCOL2, were recognized and comprehensively analyzed. However, sequence analysis showed that SiCOL2 lacked one of the B-box motifs. In addition, the flowering time of the transgenic Arabidopsis lines with overexpressed SiCOL2 were longer than that of SiCOL1, indicating that SiCOL1 was more likely to be the potential functional homologue of CO in sesame. Expression analysis revealed that SiCOL1 had high expressed levels before flowering in leaves and exhibited diurnal rhythmic expression in both SD and long-day (LD) conditions. In total, 16 haplotypes of SiCOL1 were discovered in the sesame collections from Asia. However, the mutated haplotypes did not express under both SD and LD conditions and was regarded as a nonfunctional allele. Notably, the sesame landraces from high-latitude regions harboring nonfunctional alleles of SiCOL1 flowered much earlier than landraces from low-latitude regions under LD condition, and adapted to the northernmost regions of sesame cultivation. The result indicated that sesame landraces from high-latitude regions might have undergone artificial selection to adapt to the LD environment.

CONCLUSIONS

Our results suggested that SiCOL1 might contribute to regulation of flowering in sesame and natural variations in SiCOL1 were probably related to the expansion of sesame cultivation to high-latitude regions. The results could be used in sesame breeding and in broadening adaptation of sesame varieties to new regions.

Comparative transcriptomics provides insight into the molecular basis of species diversification of section Trigonopedia (Cypripedium) on the Qinghai-Tibetan Plateau

Deceptive pollination is key to the species richness of Orchidaceae. However, the genetic basis of species diversification is still under study. Section Trigonopedia is a monophyletic clade of genus Cypripedium distributed in the southwest of China. The species of this section are pollinated by different flies. Pollinator differentiation makes section Trigonopedia an ideal group for studying the genetic basis underlying species diversification. Here, we sequenced the transcriptomes of eight species of the genus Cypripedium, including six co-flowering species of section Trigonopedia and two species outside this section as an outgroup. We reconstructed the phylogeny of the section with the combined 1572 single-copy genes extracted from the eight species and produced a highly resolved tree of the section. Furthermore, we combined substitution rate estimation and differential expression analysis to identify candidate genes, including genes related to floral scent synthesis and environmental adaptation, involved in species differentiation. Field investigations showed that these species have adapted to different habitats. We propose that the species diversification in this section is initiated by floral scent differentiation, followed by habitat differentiation, finally leading to speciation. This study sheds novel light on the diversification of closely related orchid species in the Qinghai-Tibetan region.

Post genomics era for orchid research

Among 300,000 species in angiosperms, Orchidaceae containing 30,000 species is one of the largest families. Almost every habitats on earth have orchid plants successfully colonized, and it indicates that orchids are among the plants with significant ecological and evolutionary importance. So far, four orchid genomes have been sequenced, including Phalaenopsis equestris, Dendrobium catenatum, Dendrobium officinale, and Apostaceae shengen. Here, we review the current progress and the direction of orchid research in the post genomics era. These include the orchid genome evolution, genome mapping (genome-wide association analysis, genetic map, physical map), comparative genomics (especially receptor-like kinase and terpene synthase), secondary metabolomics, and genome editing.

Characterising genes associated with flowering time in carrot (Daucus carota L.) using transcriptome analysis

Carrot is generally regarded as a biennial plant with an obligatory vernalization requirement. Early spring cultivation makes plants vulnerable to premature bolting, which results in a loss of commercial value. However, our knowledge of flowering time genes and flowering mechanisms in carrot remain limited. Bolting behavior of D. carota ssp. carota 'Songzi', a wild species sensitive to flower induction by vernalization and photoperiod, and orange cultivar 'Amsterdam forcing', and their offspring were investigated in different growing conditions. We performed RNA-seq to identify the flowering time genes, and digital gene expression (DGE) analysis to examine their expression levels. The circadian patterns of related genes were identified by qPCR. The results showed bolting behavior of carrot was influenced by low temperature, illumination intensity and photoperiod. A total of 45 flowering time-related unigenes were identified, which were classified into five categories including photoperiod, vernalization, autonomous and gibberellin pathway, and floral integrators. Homologs of LATE ELONGATED HYPOCOTYL (LHY) and CONSTANS-LIKE 2 (COL2) were more highly expressed under short day condition than under long day condition. Homologs of COL2, CONSTANS-LIKE 5 (COL5), SUPPRESSION OF OVEREXPRESSION OF CONSTANS 1 (SOC1), FLOWERING LOCUS C (FLC) and GIBBERELLIC ACID INSENSITIVE (GAI) were differentially expressed between 'Songzi' and 'Amsterdam forcing'. The homolog of COL2 (Dct43207) was repressed by light, but that of COL5 (Dct20940) was induced. A preliminary model of genetic network controlling flowering time was constructed by associating the results of DGE analysis with correlation coefficients between genes. This study provides useful information for further investigating the genetic mechanism of flowering in carrot.

GBS-based single dosage markers for linkage and QTL mapping allow gene mining for yield-related traits in sugarcane

BACKGROUND

Sugarcane (Saccharum spp.) is predominantly an autopolyploid plant with a variable ploidy level, frequent aneuploidy and a large genome that hampers investigation of its organization. Genetic architecture studies are important for identifying genomic regions associated with traits of interest. However, due to the genetic complexity of sugarcane, the practical applications of genomic tools have been notably delayed in this crop, in contrast to other crops that have already advanced to marker-assisted selection (MAS) and genomic selection. High-throughput next-generation sequencing (NGS) technologies have opened new opportunities for discovering molecular markers, especially single nucleotide polymorphisms (SNPs) and insertion-deletion (indels), at the genome-wide level. The objectives of this study were to (i) establish a pipeline for identifying variants from genotyping-by-sequencing (GBS) data in sugarcane, (ii) construct an integrated genetic map with GBS-based markers plus target region amplification polymorphisms and microsatellites, (iii) detect QTLs related to yield component traits, and (iv) perform annotation of the sequences that originated the associated markers with mapped QTLs to search putative candidate genes.

RESULTS

We used four pseudo-references to align the GBS reads. Depending on the reference, from 3,433 to 15,906 high-quality markers were discovered, and half of them segregated as single-dose markers (SDMs) on average. In addition to 7,049 non-redundant SDMs from GBS, 629 gel-based markers were used in a subsequent linkage analysis. Of 7,678 SDMs, 993 were mapped. These markers were distributed throughout 223 linkage groups, which were clustered in 18 homo(eo)logous groups (HGs), with a cumulative map length of 3,682.04 cM and an average marker density of 3.70 cM. We performed QTL mapping of four traits and found seven QTLs. Our results suggest the presence of a stable QTL across locations. Furthermore, QTLs to soluble solid content (BRIX) and fiber content (FIB) traits had markers linked to putative candidate genes.

CONCLUSIONS

This study is the first to report the use of GBS for large-scale variant discovery and genotyping of a mapping population in sugarcane, providing several insights regarding the use of NGS data in a polyploid, non-model species. The use of GBS generated a large number of markers and still enabled ploidy and allelic dosage estimation. Moreover, we were able to identify seven QTLs, two of which had great potential for validation and future use for molecular breeding in sugarcane.

Transcriptome-wide analysis of the MADS-box gene family in the orchid Erycina pusilla

Orchids exhibit a range of unique flower shapes and are a valuable ornamental crop. MADS-box transcription factors are key regulatory components in flower initiation and development. Changing the flower shape and flowering time can increase the value of the orchid in the ornamental horticulture industry. In this study, 28 MADS-box genes were identified from the transcriptome database of the model orchid Erycina pusilla. The full-length genomic sequences of these MADS-box genes were obtained from BAC clones. Of these, 27 were MIKC-type EpMADS (two truncated forms) and one was a type I EpMADS. Eleven EpMADS genes contained introns longer than 10 kb. Phylogenetic analysis classified the 24 MIKC(c) genes into nine subfamilies. Three specific protein motifs, AG, FUL and SVP, were identified and used to classify three subfamilies. The expression profile of each EpMADS gene correlated with its putative function. The phylogenetic analysis was highly correlated with the protein domain identification and gene expression results. Spatial expression of EpMADS6, EpMADS12 and EpMADS15 was strongly detected in the inflorescence meristem, floral bud and seed via in situ hybridization. The subcellular localization of the 28 EpMADS proteins was also investigated. Although EpMADS27 lacks a complete MADS-box domain, EpMADS27-YFP was localized in the nucleus. This characterization of the orchid MADS-box family genes provides useful information for both orchid breeding and studies of flowering and evolution.

Digital Gene Expression Analysis Based on De Novo Transcriptome Assembly Reveals New Genes Associated with Floral Organ Differentiation of the Orchid Plant Cymbidium ensifolium

Cymbidium ensifolium belongs to the genus Cymbidium of the orchid family. Owing to its spectacular flower morphology, C. ensifolium has considerable ecological and cultural value. However, limited genetic data is available for this non-model plant, and the molecular mechanism underlying floral organ identity is still poorly understood. In this study, we characterize the floral transcriptome of C. ensifolium and present, for the first time, extensive sequence and transcript abundance data of individual floral organs. After sequencing, over 10 Gb clean sequence data were generated and assembled into 111,892 unigenes with an average length of 932.03 base pairs, including 1,227 clusters and 110,665 singletons. Assembled sequences were annotated with gene descriptions, gene ontology, clusters of orthologous group terms, the Kyoto Encyclopedia of Genes and Genomes, and the plant transcription factor database. From these annotations, 131 flowering-associated unigenes, 61 CONSTANS-LIKE (COL) unigenes and 90 floral homeotic genes were identified. In addition, four digital gene expression libraries were constructed for the sepal, petal, labellum and gynostemium, and 1,058 genes corresponding to individual floral organ development were identified. Among them, eight MADS-box genes were further investigated by full-length cDNA sequence analysis and expression validation, which revealed two APETALA1/AGL9-like MADS-box genes preferentially expressed in the sepal and petal, two AGAMOUS-like genes particularly restricted to the gynostemium, and four DEF-like genes distinctively expressed in different floral organs. The spatial expression of these genes varied distinctly in different floral mutant corresponding to different floral morphogenesis, which validated the specialized roles of them in floral patterning and further supported the effectiveness of our in silico analysis. This dataset generated in our study provides new insights into the molecular mechanisms underlying floral patterning of Cymbidium and supports a valuable resource for molecular breeding of the orchid plant.

Identification of a CONSTANS homologous gene with distinct diurnal expression patterns in varied photoperiods in ramie (Boehmeria nivea L. Gaud)

Ramie is an important natural fiber and forage crop in China. Breeding of late- or non-flowering varieties, with higher vegetative yields, is an important goal in ramie improvement. However, the ramie genes involved in flowering regulation have not previously been identified. In model plants, such as rice and Arabidopsis, the CONSTANS (CO) and CONSTANS-like (COL) genes play key roles in flowering regulation. In the present study, six ramie COL genes (BnCOL1-BnCOL6) with a full-length open reading frame (ORF) were identified. Sequence alignment revealed that all six BnCOL proteins contained conserved CCT (CO, COL, TOC1) and B-box I domains, but that only four of these proteins contained the B-box II domain. Expression pattern analysis showed that BnCOL1-BnCOL6 were mainly expressed in the stem xylem, flowers, and leaves. Phylogenetic analysis classified the six newly identified BnCOL proteins, and also 16 COL proteins with known functions in other species, into three groups. The BnCOL2-encoded protein was assigned to the same group as the CO- and Hd1-encoded proteins, suggesting that this BnCOL2-encoded protein showed the highest level of homology with the CO/Hd1-encoded proteins. Photoperiodic experiments showed that BnCOL2 exhibited a diurnal expression pattern under long- and short-day conditions. Subcellular localization examination revealed that the BnCOL2 protein fused with YFP was localized in the nucleus. Because the homologous sequence and similar expression pattern between BnCOL2 and CO/Hd1, the BnCOL2 possibly has a role in flowering modulation, and can be used as a candidate gene for research in the flowering regulation of ramie.

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.

BeMADS1 is a key to delivery MADSs into nucleus in reproductive tissues-De novo characterization of Bambusa edulis transcriptome and study of MADS genes in bamboo floral development

BACKGROUND

The bamboo Bambusa edulis has a long juvenile phase in situ, but can be induced to flower during in vitro tissue culture, providing a readily available source of material for studies on reproductive biology and flowering. In this report, in vitro-derived reproductive and vegetative materials of B. edulis were harvested and used to generate transcriptome databases by use of two sequencing platforms: Illumina and 454. Combination of the two datasets resulted in high transcriptome quality and increased length of the sequence reads. In plants, many MADS genes control flower development, and the ABCDE model has been developed to explain how the genes function together to create the different whorls within a flower.

RESULTS

As a case study, published floral development-related OsMADS proteins from rice were used to search the B. edulis transcriptome datasets, identifying 16 B. edulis MADS (BeMADS). The BeMADS gene expression levels were determined qRT-PCR and in situ hybridization. Most BeMADS genes were highly expressed in flowers, with the exception of BeMADS34. The expression patterns of these genes were most similar to the rice homologs, except BeMADS18 and BeMADS34, and were highly similar to the floral development ABCDE model in rice. Transient expression of MADS-GFP proteins showed that only BeMADS1 entered leaf nucleus. BeMADS18, BeMADS4, and BeMADS1 were located in the lemma nucleus. When co-transformed with BeMADS1, BeMADS15, 16, 13, 21, 6, and 7 translocated to nucleus in lemmas, indicating that BeMADS1 is a key factor for subcellular localization of other BeMADS.

CONCLUSION

Our study provides abundant B. edulis transcriptome data and offers comprehensive sequence resources. The results, molecular materials and overall strategy reported here can be used for future gene identification and for further reproductive studies in the economically important crop of bamboo.

Global transcriptome analysis and identification of the flowering regulatory genes expressed in leaves of Lagerstroemia indica

Flowering time is an important trait for ornamental plants, and flowering regulation has thus been both a focus of and challenge to researchers. Lagerstroemia indica is an important summer flowering tree in China and has been introduced abroad as a key parent of new cultivars; no previous reports have addressed the regulation of flowering time in this species. In this study, 28,567,778×2 reads were obtained from leaves of L. indica. A total of 37,325 unigenes were assembled with an average length of 849.56 bp, and 17,506 (46.90%) unigenes were significantly matched to known genes in the nr database of GenBank. The annotated sequences were clustered into putative functional categories using the Gene Ontology framework. Potential genes and their functions were predicted by the Cluster of Orthologous Groups analysis and Kyoto Encyclopedia of Genes and Genomes pathway mapping. A total of 115 unigenes related to flowering time control were discovered. Ten homologous genes of the CONSTANS-like (COL) gene family were identified based on transcript data. Phylogenetic analysis of the CONSTANS and COL genes from L. indica and other species grouped them into three clades. The transcriptome dataset and outcome of the analysis provide a valuable new resource for research on the functional genomics and molecular mechanisms of flowering control in L. indica.

The analysis of the inflorescence miRNome of the orchid Orchis italica reveals a DEF-like MADS-box gene as a new miRNA target

Plant microRNAs (miRNAs) are small, regulatory non-coding RNAs involved in a wide range of biological processes, from organ development to response to stimuli. In recent years, an increasing number of studies on model plant species have highlighted the evolutionary conservation of a high number of miRNA families and the existence of taxon-specific ones. However, few studies have examined miRNAs in non-model species such as orchids, which are characterized by highly diversified floral structures and pollination strategies. Therefore, we analysed a small RNA library of inflorescence tissue of the Mediterranean orchid Orchis italica to increase the knowledge on miRNAs in a non-model plant species. The high-throughput sequencing and analysis of a small RNA library of inflorescence of O. italica revealed 23 conserved and 161 putative novel miRNA families. Among the putative miRNA targets, experimental validation demonstrated that a DEF-like MADS-box transcript is cleaved by the homolog of miR5179 of O. italica. The presence of conserved miRNA families in the inflorescence of O. italica indicates that the basic developmental flower regulatory mechanisms mediated by miRNAs are maintained through evolution. Because, according to the "orchid code" theory, DEF-like genes exert a key function in the diversification of tepals and lip, the cleavage-mediated inhibitory activity of miR5179 on a OitaDEF-like transcript suggests that, in orchids, miRNAs play an important role in the diversification of the perianth organs.

Shedding some light over the floral metabolism by arum lily (Zantedeschia aethiopica) spathe de novo transcriptome assembly

Zantedeschia aethiopica is an evergreen perennial plant cultivated worldwide and commonly used for ornamental and medicinal purposes including the treatment of bacterial infections. However, the current understanding of molecular and physiological mechanisms in this plant is limited, in comparison to other non-model plants. In order to improve understanding of the biology of this botanical species, RNA-Seq technology was used for transcriptome assembly and characterization. Following Z. aethiopica spathe tissue RNA extraction, high-throughput RNA sequencing was performed with the aim of obtaining both abundant and rare transcript data. Functional profiling based on KEGG Orthology (KO) analysis highlighted contigs that were involved predominantly in genetic information (37%) and metabolism (34%) processes. Predicted proteins involved in the plant circadian system, hormone signal transduction, secondary metabolism and basal immunity are described here. In silico screening of the transcriptome data set for antimicrobial peptide (AMP) –encoding sequences was also carried out and three lipid transfer proteins (LTP) were identified as potential AMPs involved in plant defense. Spathe predicted protein maps were drawn, and suggested that major plant efforts are expended in guaranteeing the maintenance of cell homeostasis, characterized by high investment in carbohydrate, amino acid and energy metabolism as well as in genetic information.